CT重建实验报告.docxVIP

平行光照射的CT重建 艾文书 220171807 实验目的 认识并熟悉CT重建的原理和成像基础。通过编写CT图像重建程序，进一步了解CT重建过程，同时加强图像处理程序的编程训练。 实验软件：MATLAB 实验要求1.递交整个过程的源程序。比较原图和旋转360次（一度一次） 2.采用频率滤波反投影和时域卷积反投影得到重建结果。 四、算法原理及结果 1.CT重建原理 在CT成像中，物体对X线吸收起主要作用，在一个均匀的物体中，X线的衰减服从指数规律。在X线穿透人体器官或组织时，由于人体器官或组织是由多种物质成分和不同的密度构成的，所以各点对X线的吸收系数是不同的。将沿着X线束通过的物体分割成许多小单元体（体素），令每个体素的厚度相等(l)。设足够小，使得每个体素均匀，每个体素的吸收系数为常值，如果X线的入射强度I0、透射强度I和物体体素的厚度l均为已知，沿着X线通过路径上的吸收系数之和μ1＋μ2＋……＋μn就可计算出来。 吸收系数是一个物理量，它是CT影像中每个像素所对应的物质对X线线性平均衰减量大小的表示。再将图像面上各像素的CT值转换为灰度，就得到了图像面上的灰度分布，就是CT影像。 CT重建过程可以采用直接反投影和卷积反投影来实现。卷积反投影重建图像时，先把由检测器上获得的原始数据与一个滤波函数进行了卷积运算，得到各方向卷积的投影函数；然后再把它们从各方向进行反投影，即按其原路径平均分配到每一个矩阵元上，进行叠加后得到每一矩阵元的CT值；再经过适当处理后就可以得到被扫描物体的断层图像，卷积反投影可消除单纯的反投影产生的边缘失锐效应，补偿投影中的高频成分和降低投影中心密度，并保证重建图像边缘清晰和内部分布均匀。 2.图示CT重建原理（重要步骤部分原理） 五、MATLAB程序算法及结果 六、结果评估 本次实验，总体上实现了老师的要求，但是由于程序简单，对实现的效果并不是特别明显，在图形重建方面知识还是有待于加强与提高。对得到的结果，原图与360每旋转1度角来形成投影结果，我们可以看到最后的重建结果差别并不是很明显。