matlab 图像合成_MATLAB图像处理：14：从投影数据重建图像

此示例示出了如何使用radon，iradon，fanbeam，和ifanbeam，从样本图像建立投影，屏从投影重建图像。同时radon，iradon对投影使用平行光束几何，fanbeam并ifanbeam使用扇形光束几何。为了比较平行光束几何和扇形光束几何，以下示例为每种几何创建了合成投影，然后使用这些合成投影来重建原始图像。

需要重建图像的实际应用，是X射线吸收层析成像，其中通过测量以不同角度穿过物理样本的辐射的衰减来形成投影。原始图像可以认为是穿过样品的横截面，其中强度值代表样品的密度。通过特殊的医学成像设备收集投影，然后使用iradon或ifanbeam重建标本的内部图像。

函数iradon从平行光束投影重建图像。在平行光束几何中，每个投影都是通过以特定角度组合通过图像的一组线积分来形成的。函数ifanbeam从具有一个发射器和多个传感器的扇形束投影中重建图像。

请参阅《 Image Processing Toolbox™用户指南》以获取说明两种几何形状的图表。

创建头部影像

测试图像是可以使用函数生成的Shepp-Logan头部模型phantom。图像说明了在现实世界的人体头部层析成像中发现的许多质量。沿外部的明亮椭圆形外壳类似于头骨，内部的许多椭圆类似于脑部特征或肿瘤。

 P = phantom(256);
 imshow(P)

平行光束-计算合成投影

使用平行光束几何计算合成投影并改变投影角度的数量。对于每个radon调用，输出都是一个矩阵，其中每一列是对应theta角度中