毕业论文-基于MATLAB的扇形束投影CT重建

1、表 44对比的PSNR值可知，当重建对象与重建参数相同时，等角型扇束重建出的断层图像质量明显优于等角型扇束，尤其是当旋转增角和探测器间距都很大的时候等距型就有更明显的优势。参考文献Radon J. MathJ. Physics. 1917, 21(6): 262-267.Bracewell R N. Image Reconstruction in Radio AstronomyJ. Implementation and Application. 1979, 13(7): 465-474.Cormak A M. Representation of a function by its line i

2、ntegrates with some radiological applicationsJ. Journal of Applied Physics. 1956, 34(9): 1635-1641.Herman G. T. Image reconstruction from projectionsM. New York: Academic Press, 1980.吕东辉. 变旋转中心扫描及BPF的局部特性研究D. 上海大学, 2006: 69-70. 张朝宗, 郭志平, 张鹏, 王贤刚等. 工业CT技术和原理M. 北京: 科学出版社, 2009.王黎明. 直接重建目标特征的CT算法研究D. 中

3、北大学, 2009: 54-55.崔屹. 数字图像处理技术与应用M. 北京:电子工业出版社, 1997.潘晋孝. 锥形束三维XCT重建算法研究D. 中北大学, 2007: 67-70.堵国樑. CT式行李安检系统反投影算法研究及其硬件设计D. 东南大学, 2010: 52-53.杨昆. 基于螺旋Micro-CT的成像质量优化研究D. 河北大学, 2013: 65-67.Henrik Turbell. Cone-beam reconstruction using filtered backprojectionD. Sweden: Department of Electrical Engineer

4、ing Link pings universitet, 2001: 23-34.庄天戈. CT原理与算法M. 上海: 上海交通大学出版社, 1992.庄天戈. 计算机在生物医学中的应用M. 北京: 科学出版社, 2000.致 谢至此，本科毕业设计的任务基本圆满完成。在进行毕业设计的这一个季度中，经历过许多意想不到的困难与挑战，但最终还是通过自己不懈的努力与同学老师的热心帮助顺利将其完成。回想这一段时间，感觉自己在科学研究收获很多。虽然不久后我将离开学校，但是我十感激关心我的老师与同学，感谢他们在此期间给予我的持久的帮助与鼓励。首先我感谢敖波老师对我的帮助，正是他的精心选题与对毕设严格要求才使我

5、顺利完成本次的毕业设计。同时，敖老师严谨治学的态度使我更加明白科研人员应有的责任与担当。其次，我要感谢检测楼的诸位老师和同学们，由于他们的热情帮助，使得我能一次又一次解决毕设期间遇到的问题。最后我要感谢110813班的同学们在我攥写毕业论文时给我的宝贵意见与建议。感谢昌航给我提供的进行毕业设计的机会与条件。也衷心感谢诸位参加评阅论文和答辩的专家、教授。附录 实现扇形束CT投影重建的MATLAB程序实现Shepp-logan脑图模型的扇形束CT重建clear; A = phantom(256); D=250;FS=1;OS=256;RI=1;SS=1;F=fanbeam(A,D,FanSenso

6、rGeometry,arc,. FanRotationIncrement,RI,. FanSensorSpacing,SS);AReconstructed=ifanbeam(F,D,FanSensorGeometry,arc,. FanRotationIncrement,RI,. FanSensorSpacing,SS,. FrequencyScaling,FS,. Filter,shepp-logan,. OutputSize,OS);m1=PSNR(AReconstructed,A);imwrite(AReconstructed,D:DocumentsMATLAB扇形束a0.1.JPG);

7、计算峰值信噪比（PSNR）function output=PSNR(input1,input2)m,n=size(input1);input1=double(input1);input2=double(input2);d=0;e=0;a=zeros(m,n);b=zeros(m,n);for i=1:m for j=1:n a(i,j)=2552; b(i,j)=input1(i,j)-input2(i,j); d=d+a(i,j); e=e+b(i,j)2; endendoutput=10*log10(d/e);利用MATLAB自建断层模型clear;n=256;A=zeros(256,256);A(20:236,20:236)=1;A(30:216,30:216)=0.2;A(150:190,50:100)=0.8;for x=1:n for y=1:n if (x-n/3).2/800+(y-n/3).2/2500=1 A(x,y)=1; end if (x-n/3).2/100+(y-n*2/3).2/100