实验二扇形束投影数据采集与滤波反投影重建实验.pdfVIP

南昌大学实验报告 学生姓名： 胡文松 学 号： 6103413007 专业班级： 生医 131班 实验类型：□验证□综合■设计□创新 实验日期： 实验成绩： 实验二 、扇形束投影数据采集与滤波反投影重建实验 一、实验目的及要求： 1、 用 MA TLAB 中的 fanbeam 函数获取 Shepp-Logan 头模型的扇形束投影数据； 2、 显示扇形束投影数据并和平行束投影数据进行比较； 3、 用 MA TLAB 中的 ifanbeam 函数实现图像重建； 4、 改变投影参数，观察对重建图像的影响。 二、实验基本原理： fanbeam 函数，计算扇形束投影。 句法 F = fanbeam(I,D) F = fanbeam(...,param1,val1,param1,val2,...) [F,sensor_positions,fan_rotation_angles] = fanbeam(...) F = fanbeam(I,D) 计算图像 I 的扇形束投影 (正弦图 ) ， D 是扇形束顶点到旋转中心的距 离。 F 的每一列表示一个角度下的扇形采样数据。探测器圆弧状排布，探测器间相隔一度， 旋转 360 度 [0:359] 。 F = fanbeam(...,param1,val1,param1,val2,...) 指定参数控制扇形投影 FanRotationIncrement 指定旋转角度增量，单位为度，默认值为 1 度。 FanSensorGeometry 指定探测器的排布方式 arc—探测器圆弧状排布（默认值） line -- 探测器等距离线性排布 FanSensorSpacing 指定探测器间的间隔，探测器圆弧状排布的情况，单位为度，默认值 为 1。探测器等距离线性排布，单位为像素。 [F,sensor_positions,fan_rotation_angles] = fanbeam(...) 返回探测器位置和旋转角度信息 ifanbeam 函数，由扇形束投影数据重建图像 句法 I = ifanbeam(F,D) I = ifambeam(...,param1,val1,param2,val2,...) [I,H] = ifanbeam(...) I = ifanbeam(F,D) 由投影数据 F 重建图像 I 。 D 是扇形束顶点到旋转中心的距离 I = ifanbeam(...,param1,val1,param2,val2,...) 指定参数控制图像重建。 FanCoverage 指定旋转角度范围 cycle} — 旋转 360 度 [0,360). minimal –重建图像所需的最小旋转角度。 FanRotationIncrement 同 fanbeam FanSensorGeometry 同 fanbeam FanSensorSpacing 同 fanbeam OutputSize 指定重建图像矩阵的大小， 假如不指定， ifanbeam 自动决定矩阵大小 （扇形投影重建图像和原始图像矩阵大小可能不同） 。 三、主要仪器设备及实验耗材： 具有 XP 或 2000 系统，并装有 MATLAB系统的 PC机。 四、实验内容或步骤以及显示结果： 1．用 MATLAB图像处理工具箱的 phantom 生成 Shep-Logan 头模型； P=phantom(256); imshow(P) 2．用 MATLAB中的 fanbeam 函数获得头模型的扇形束和平行束投影数据； D = 200; F = fanbeam(P,D,FanRotationIncrement,0.5,FanSensorSpacing,0.5); theta=0: