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 度角来形成投影结果，我们 可以看到最后的重建结果差别并不是很明显。