CT重建-X射线断层成像仿真实验

今天来做一下医学成像系统大作业：CT重建-X射线断层成像仿真实验。

一. 对一幅仿体断层图像进行平行束扫描,获得并显示投影数据

1. 首先用phantom函数在matlab中建一个仿体断层图像，phantom函数具体用法如下：

（这个MATLAB函数生成一个头部仿体断层图像，可以用来测试radon和iradon或其他二维重建算法的数值精度）

其中， P是由一个大的椭圆(代表大脑)和几个小的椭圆(代表大脑中的特征)组成的灰度图像；Def指定要生成的头部仿体断层图像的类型，n指定仿体断层图像中的行数和列数。

P = phantom(E,n)生成一个用户定义的phantom，其中矩阵E的每一行指定图像中的一个椭圆。E有六列，每列包含椭圆的不同参数。

[P,E] =phantom(___)返回用于生成仿体断层图像的矩阵E。

~ 一个最经典的仿体断层图像模型就是Shepp-Logan模型，matlab里也给了相应的示例：

2. 用Radon变换对仿体断层图像进行平行束扫描，获得投影数据。

以下是Radon变换的matlab文档：（注：Radon变换具体原理见成像系统课件-“2021-04-06 CT成像”）