医学影像系统实验PART II部分报告 U201012579 生医1002班刘志 - 图文

椭圆模型 Shepp-Logan头部

四、实验方法与步骤：

Ctsim软件是由Kevin M.Rosenberg等学者编写开发的CT重建仿真程序，其源代码开放。用户可选择Herman head、shepp-logan、unit pulse三种体模，进行滤波反投影重建和傅里叶变换重建。

1．选择要进行扫描的体模。打开Ctsim软件，按file-create phantom选择体模shepp-logan（如图1）。

图2．对体模进行平行线束（平移-旋转式）扫描投影。按process-projections,在projection parameter中选择geometry-parallel（平行束）, trace level中选择projections, 其他扫描参数使用默认参数。软件进行体模扫描演示（如图2），扫描完成后，将各角度投影值按顺序排列在二维平面上（图3）。

1

图3 图3．选择reconstruct-filtered backprojection,在filtered backprojection parameter 中选filter method 为FFT。软件进行滤波反投影重建，选用的滤波方法为快速傅里叶变换，获得重建的sheep-logan体模图像（如图4）。

2

4．比较重建图像与原体模的差异。

5．在trace level中选择plot，重复1~3步，观察扫描过程中，每个投影值的变化。

6．在projection parameter中选择geometry-equiangule(扇形扫描)，重复1~3步，观察扫描过程和重建影像。其扫描过程如图5、投影图象如图6、重建图像如图7。比较投影图象、重建图像与平行线束扫描获得的图像的区别。

图4

图5 图7 图6

7．选用unit pulse体模，重复上述过程，观察重建图像中的伪影。

8．更换其它体模，进行上述方法，进行重建，观察重建方法与重建图象。

五、实验结果显示及其分析：

(1) shepp-logan体模如下所示：

图一

对该图进行光栅化，得到的图形如下所示:

图二

(2)对体模进行平行线束（平移-旋转式）扫描投影过程中扫描演示图以及所得到的投影图像如下所示:

图三

图四

(3)用滤波反投影方法重建出的shepp-logan体模图像如下所示: