医学影像系统实验-CTsim图像重建实验.pptVIP

CT图像重建仿真实验;一、实验目的;三、实验原理;;;由于只是模拟X射线的投影过程，我们简化了问题。假设断面的结构如图1（Shepp-Logan）所示（各图元均为椭圆），各个椭圆表示了人体的不同的组织（内部是均匀的），分别有不同的线性衰减系数μ。那么，就可利用公式②来求某条X射线投影值。各个椭圆（组织）的线性衰减系数μ（Shepp-Logan图的各椭圆的位置、大小和线性衰减系数参见表2.1）是已知的，问题就是求X射线穿过椭圆时的行进距离。 设椭圆的长短轴为a，b；X射线与a的夹角为Φ；椭圆中心到X射线的距离为d。如下图所示。这样可由④⑤求得X射线穿过椭圆时的行进距离。再乘上各个椭圆的线性衰减系数μ后累加起来就可得到X射线的投影值。;为了对反投影重建后的图像进行评价，可对Shepp-Logan图进行光栅化。然后比较重建前后图像的不同。更多信息参见matlab 中phantom radon iradon ???程，自由软件CTSim。;四、实验方法与步骤;对体模进行平行线束（平移-旋转式）扫描投影。按process-projections,在projection parameter中选择geometry-parallel（平行束）, trace level中选择projections, 其他扫描参数使用默认参数。 软件进行体模扫描演示（如图2），扫描完成后，将各角度投影值按顺序排列在二维平面上（图3）。;选择reconstruct-filtered backprojection,在filtered backprojection parameter 中选择filter method 为FFT。 软件进行滤波反投影重建，选用的滤波方法为快速傅里叶变换，获得重建的sheep-logan体模图像（如图4）。;比较重建图像与原体模的差异。 在trace level中选择plot，重复1~3步，观察扫描过程中，每个投影值的变化。 在projection parameter中选择geometry-equiangule(扇形扫描)，重复1~3步，观察扫描过程和重建影像。其扫描过程如图5、投影图象如图6、重建图像如图7。比较投影图象、重建图像与平行线束扫描获得的图像的区别。;选用unit pulse体模，重复上述过程，观察重建图像中的伪影。 更换其它体模(自己建立模型)，进行上述方法，进行重建，观察重建方法与重建图象。 ;对象