CT重建实验报告.doc

CT重建实验报告.doc:..平行光照射的CT重建艾文书220171807一、 实验目的认识并熟悉CT重建的原理和成像基础。通过编写CT图像重建程序,进一步了解CT重建过程,同时加强图像处理程序的编程训练。二、 实验软件:MATLAB三、 实验要求1•递交整个过程的源程序。比较原图和旋转360次(一度一次)2•釆用频率滤波反投影和时域卷积反投影得到重建结果。四、 ,物体对X线吸收起主要作用,在一个均匀的物体中,X线的衰减服从指数规律。在X线穿透人体器官或组织时,由于人体器官或组织是由多种物质成分和不同的密度构成的,所以各点对X线的吸收系数是不同的。将沿着X线束通过的物体分割成许多小单元体(体素),令每个体素的厚度相等(I)。设足够小,使得每个体素均匀,每个体素的吸收系数为常值,如果X线的入射强度10、透射强度I和物体体素的厚度I均为已知,沿着X线通过路径上的吸收系数之和P1+U2++un就可计算出来。吸收系数是一个物理量,它是CT影像中每个像素所对应的物质对X线线性平均衰减量大小的表示。再将图像面上各像素的CT值转换为灰度,就得到了图像面上的灰度分布,就是CT影像。CT重建过程可以采用肓接反投影和卷积反投影来实现。卷积反投影重建图像时,先把由检测器上获得的原始数据与一个滤波函数进行了卷积运算,得到各方向卷积的投影函数;然后再把它们从各方向进行反投影,即按其原路径平均分配到每一个矩阵元上,进行叠加后得到每一矩阵元的CT值;再经过适当处理后就可以得到被扫描物体的断层图像,卷积反投影可消除单纯的反投影产生的边缘失锐效应,补偿投影中的高频成分和降低投影中心密度,(重要步骤部分原理)Sinogram:所有9角度上的Projection:CT发出的X射线 =densityof25/Xw)□-densityof50eJ第一个Project!on五、MATLAB程序算法及结果塞辑器 发布 ;|+%,clear:% radon变换前的囹像 P=phantomCModifiedshepp-logan?,256):figure,imshow(P):title('原图’);%Radon变换:theta=linspace(0,179,360);R=radon(P,theta):%利用;Lradon变换用360个角度来重构囹像Il=iradon(R,theta):figure,imshow(Il):title(?360个角度重构图像');% Iradon %p=phantomCModifiedshepp-logan7,256):A=linspace(0,179,180);1R,xp>radon(p,A):% 方法1・频域廳波法 lien,ANJ=size(R):lenFilt=2nextpow2(2*size(R,1)):%滤波器的长度n=0:(lenFilt/2):w=2*pi*n/lenFilt;filt=Zw(end:-1:2)5:w(l:1:end-1)7];%