Author: Not specified Language: matlab
Description: Not specified Timestamp: 2017-10-11 04:56:04 +0000
View raw paste Reply
  1. % A1PartD
  2. %Comb phantom with resolution (-10:0.1:10)
  3. [xmat,ymat] = meshgrid(-10:0.1:10);
  4. h1 = (1/4).*exp(-(xmat.^2+ymat.^2))+abs((xmat.*sinc(xmat.^2)).*(sinc(ymat/16)/16));
  5. mesh(xmat, ymat, h1)
  6. title('Cyberdyne CT System PSF (h1) applied with Comb Phantom with Resolution of (-10:0.1:10)');
  7. h2 = (1/10).*exp(-((xmat-1).^2)/8)+(1/2).*exp(-((ymat-3).^2)/16);
  8. mesh(xmat,ymat,h2)
  9. title('Umbrella Corporation CT System PSF (h2) applied with Comb Phantom with Resolution of (-10:0.1:10)');
  10.  
  11. %Comb phantom with resolution (-20:0.1:20)
  12. [xmat,ymat] = meshgrid(-20:0.1:20);
  13. h1 = (1/4).*exp(-(xmat.^2+ymat.^2))+abs((xmat.*sinc(xmat.^2)).*(sinc(ymat/16)/16));
  14. mesh(xmat, ymat, h1)
  15. title('Cyberdyne CT System PSF (h1) applied with Comb Phantom with Resolution of (-20:0.1:20)');
  16. h2 = (1/10).*exp(-((xmat-1).^2)/8)+(1/2).*exp(-((ymat-3).^2)/16);
  17. mesh(xmat,ymat,h2)
  18. title('Umbrella Corporation CT System PSF (h2) applied with Comb Phantom with Resolution of (-20:0.1:20)');
  19.  
  20. %Comb phantom with resolution (-30:0.1:30)
  21. [xmat,ymat] = meshgrid(-30:0.1:30);
  22. h1 = (1/4).*exp(-(xmat.^2+ymat.^2))+abs((xmat.*sinc(xmat.^2)).*(sinc(ymat/16)/16));
  23. mesh(xmat, ymat, h1)
  24. title('Cyberdyne CT System PSF (h1) applied with Comb Phantom with Resolution of (-30:0.1:30)');
  25. h2 = (1/10).*exp(-((xmat-1).^2)/8)+(1/2).*exp(-((ymat-3).^2)/16);
  26. mesh(xmat,ymat,h2)
  27. title('Umbrella Corporation CT System PSF (h2) applied with Comb Phantom with Resolution of (-30:0.1:30)');
View raw paste Reply