投影数据采集与滤波反投影重建实验报告.docVIP

实验一 投影数据采集与滤波反投影重建 一、实验目的及要求： 用程序模拟X射线的投影，获得Shepp-Logan模型的投影数据。对获得的投影数据进行滤波反投影重建，获得Shepp-Logan模型的重建图像。 二、实验基本原理： X射线穿过人体时，人体的各种组织对X射线有不同程度的，即不同的组织有不同的线性衰减系数μ。假设强度为I0的X射线穿过均匀分布衰减系数为μ的物体，行进了x的距离，强度变为I，按Beer定理有 或 若物体时分段均匀的，系数分别是μ1、μ2、μ3、...,相应的长度为x1，x2，x3，...,则下式成立： 更一般的可用下面的积分式表示： 由于只是模拟X射线的投影过程，我们简化了问题。假设断面的结构如图1.1（Shepp-Logan）所示（各图元均为椭圆），各个椭圆表示了人体的不同的组织（内部是均匀的），分别有不同的线性衰减系数μ。那么，就可利用公式来求某条X射线投影值。各个椭圆（组织）的线性衰减系数μ（Shepp-Logan图的各椭圆的位置、大小和线性衰减系数参见表.1）是已知的，问题就是球X射线穿过椭圆时的行进距离。设椭圆的长短轴为a，b；X射线与b的夹角为Φ；椭圆中心到X射线的距离为d。如图1.2所示。这样可由求得X射线穿过椭圆时的行进距离。再乘上各个椭圆的线性衰减系数μ后累加起来就可得到X射线的投影值。 ⑤ ? 图1.1图1.2 表.1Shepp-L