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));
  6. mesh(xmat, ymat, h1)
  7. title('Cyberdyne CT System PSF (h1) applied with Comb Phantom with Resolution of (-10:0.1:10)');
  9. h2 = (1/10).*exp(-((xmat-1).^2)/8)+(1/2).*exp(-((ymat-3).^2)/16);
  10. mesh(xmat,ymat,h2)
  11. title('Umbrella Corporation CT System PSF (h2) applied with Comb Phantom with Resolution of (-10:0.1:10)');
  12.  
  13. %Comb phantom with resolution (-20:0.1:20)
  14. [xmat,ymat] = meshgrid(-20:0.1:20);
  15. h1 = (1/4).*exp(-(xmat.^2+ymat.^2))+abs((xmat.*sinc(xmat.^2)).*(sinc(ymat/16)/16));
  17. mesh(xmat, ymat, h1)
  18. title('Cyberdyne CT System PSF (h1) applied with Comb Phantom with Resolution of (-20:0.1:20)');
  20. h2 = (1/10).*exp(-((xmat-1).^2)/8)+(1/2).*exp(-((ymat-3).^2)/16);
  21. mesh(xmat,ymat,h2)
  22. title('Umbrella Corporation CT System PSF (h2) applied with Comb Phantom with Resolution of (-20:0.1:20)');
  23.  
  24. %Comb phantom with resolution (-30:0.1:30)
  25. [xmat,ymat] = meshgrid(-30:0.1:30);
  26. h1 = (1/4).*exp(-(xmat.^2+ymat.^2))+abs((xmat.*sinc(xmat.^2)).*(sinc(ymat/16)/16));
  28. mesh(xmat, ymat, h1)
  29. title('Cyberdyne CT System PSF (h1) applied with Comb Phantom with Resolution of (-30:0.1:30)');
  31. h2 = (1/10).*exp(-((xmat-1).^2)/8)+(1/2).*exp(-((ymat-3).^2)/16);
  32. mesh(xmat,ymat,h2)
  33. title('Umbrella Corporation CT System PSF (h2) applied with Comb Phantom with Resolution of (-30:0.1:30)');
View raw paste Reply