基于扇形波束投影数据的图像重建

从扇形光束投影数据重建图像，使用ifanbeam函数。使用此功能，您将投影数据指定为参数，并且在创建投影数据时风扇光束投影的顶点与旋转中心之间的距离。例如，此代码重新创建图像一世从投影数据P.和距离D.．

I = ifanbeam (P、D);

默认情况下，ifanbeam函数假设扇形光束投影数据是使用圆弧扇形传感器几何体创建的，光束以1度角度间隔，投影以1度增量在整个360度范围内进行。与粉丝功能，可以使用ifanbeam参数为指定投影数据的这些特征的其他值。使用相同的值对于创建投影数据时使用的这些参数。有关这些参数的更多信息，请参阅ifanbeam．

这ifanbeam函数将扇形光束投影数据转换为平行光束投影数据FAN2PARA.函数，然后调用阿拉顿函数执行图像重建。由于这个原因，IFANFEAM功能支持某些功能万博1manbetx阿拉顿参数，它传递给阿拉顿函数。看到氡反变换有关的更多信息阿拉顿函数。

使用反扇形光束投影重建图像

此示例显示了如何使用粉丝和ifanbeam从样本图像形成投影，然后根据投影重建图像。

生成测试图像并显示它。测试图像是谢普-洛根头部幻影，可以由幻影函数。幻影图像说明了人类头部现实世界断层摄影成像中发现的许多品质。

P =幻影(256);imshow (P)

计算测试图像的扇束投影数据，使用FanSensorSpacing参数来改变传感器间距。该示例使用扇形光束弧几何，因此您可以通过指定光束的角间距来指定传感器之间的间距。第一呼叫间隔梁2度;第二个是1度;第三个是0.25度。在每次调用中，旋转中心和投影顶点之间的距离恒定为250像素。此外,粉丝以1度增量旋转中心像素周围的投影。

d = 250;dsensor1 = 2;f1 = fanbeam（p，d，“FanSensorSpacing”, dsensor1);dsensor2 = 1;F2 = fanbeam (P D“FanSensorSpacing”，dsensor2）;dsensor3 = 0.25;[F3，Sensor_POS3，FAN_ROT_ANGLES3] = FANBEAM（P，D，......“FanSensorSpacing”, dsensor3);

绘制投影数据F3．因为粉丝以0到360度的旋转角度计算投影数据，相同的图案以180度的偏移发生。从两侧采样相同的特征。

图，ImageC（FAN_ROT_ANGLES3，SENSOR_POS3，F3）COLOROMAP（热）;彩色杆xlabel（“风扇旋转角度(度)”) ylabel (“风扇传感器位置(度)”的）

使用扇形光束投影数据重建图像ifanbeam．在每次重建中，将风扇传感器间距与之前创建投影数据时使用的间距进行匹配。示例使用输出参数将每个重建的输出大小约束为原始图像的大小相同P.．在输出中，注意重建的质量如何随着投影中的光束的数量而变得更好。第一张图片，Ifan1，使用梁的2度间距来创建;第二个图像，Ifan2，使用梁的1度间距来创建;第三个图像，Ifan3，使用梁的0.25间距来创建。

output_size = max（大小（p））;ifan1 = ifanbeam（f1，d，......“FanSensorSpacing”dsensor1,'输出', output_size);图，imshow（ifan1）标题（“Ifan1”的）

Ifan2 = ifanbeam (F2 D......“FanSensorSpacing”dsensor2,'输出', output_size);图，imshow（ifan2）标题（“Ifan2”的）

ifan3 = ifanbeam（f3，d，......“FanSensorSpacing”dsensor3,'输出', output_size);图中,imshow (Ifan3)标题(“Ifan3”的）