(完整)CT图像重建

1、（完整）（完整）CT图像重建- -back=ax-1；fork=1：N/2n=2*k-1；RL(forw)=-1/(n*pi*vstep)八2；RL(back)=RL(forw)；RL(forw+1)=0；RL(back-1)=0；forw=forw+2；back=back2；endfork=1:(180/angle+1)radon2(k，:)=conv(radon1(k，：)，RL)；endradonf=radon2(：，2*N：3*N)；%generateSLfunction%radon1=zeros(k,N),radon0；forv=1：(2*N+1)n=vN1；SL(v)=-2/(pi

2、八2*vstep八2*(4*nU1)；endfork=1:(180/angle+1)radon2(k，:)=conv(radon1(k，：),SL)；endradonf=radon2(:，2N：3*N)；figure(1)subplot(321)plot(1：(2N+1)，radon1(1，:)title(投影函数(已补零)subplot(323)plot(1:(2*N+1),SL)title(S-L卷积函数)subplot(325)plot(1：(4*N+1)，radon2(1,：)title(卷积结果)Xradon1=fft(radon1(1,：);subplot(322)plot(1:(

3、2N+1),abs(Xradon1)title(频谱)XRL=fft(SL)；subplot(324)plot(1:(2N+1),abs(XRL)Xradon2=fft(radon2(1,：);subplot(326)plot(1:(4*N+1),abs(Xradon2)%iradon%fork=1:(180/angle+1)theta=(k1)*angle*pi/180;C=N/2-(N-1)*(cos(theta)+sin(theta)/2;fori=1：NR=(i-1)*cos(theta)+C;n0=floor（R）；ifn00&n0（N+1）dot=Rn0;I（i,1，k）=（1do

4、t）*radonf（k,n0）+dot*radonf（k，n0+1）；elseI（i，1,k）=nan;endforj=2:NR=R+sin（theta）;n0=floor（R）；ifn00&n0（N+1）dot=R-n0；I（i,j，k）=（1-dot）radonf（k，n0）+dotradonf（k，n0+1）；elseI（i,j，k）=nan;endendendendIfinal=sum（I,3）;Gfinal=mat2gray（Ifinal）；Gfinal=imrotate（Gfinal，90）;figure（2）imshow（Gfinal）五、实验结果与分析反投影一般步骤为:程序流程

5、:卷积反投影法结果:重建图象原图象六、心得体会通过本次的实验，对CT图象重建的基本方法之一:卷积反投影，有了进一步的认识，在实验的过程中，采用的图象是经典的Shep-Logen头模型，得到的结果与原图象相比，有一定的差异，但影响不大.有待进一步的改进算法.七、参考文献1医用电子仪器及装备医学成像系统及放射治疗装置类，唐庆玉主编，清华大学电机系生物医学工程及仪器专业2医学成像系统，高上凯主编，清华大学出版社3G.T。赫尔曼著（严洪编译），由投影重建图象，科学出版社4董雏申、吴世法、王天童，卷积反投影法从x光图像重建三维轴对称图象，全国图基科学会议论文集，l989。5G,NHousfleldComputedMedicalImagingJournal0fComputedAssistedTomography，4（5）,65567419806赶荣椿等，数字图像处理导论，百北工业大学出版牡，西安，1995,p771797庄天戈，CT原理与算法），上海交通大学出版杜，1992,p31-338RadonJ，UberdieBestimmumgvonFunktionenmuchihreintegralwertelangsgewissserMannnm