滤波反投影法:
滤波反投影法根据附件三所给接收信息,采用先修正、后投影重建图像的做法,可得到原始图像的吸收率信息。
其原理为:在得到某一角度下的投影函数(一维函数)后,对此函数做滤波处理,得一修正后的滤波函数,再对修正后的滤波函数做反投影运算,得待检测介质吸收率在正方形托盘中的每一点的分布密度函数。
图1给出了滤波反投影法重建原始图像的流程图。
图1滤波反投影法流程图
反投影法重建原始图像的步骤:
(1)在对应于投影函数的角度下对投影函数做一维Fourier变换;
(2)对(1)得到的变换结果乘以权重因子;
(3)对(2)加权后得到的结果做一维傅立叶;
(4)对(3)所得函数做直接反投影;
(5)改变投影角度,得到180个不同的投影角度,对每一角度,重复上述步骤(1)~(4)。
R-L(Ram-Lak)滤波函数:
此函数的基本条件是二维图像函数的频率是有界的,显然,此题所得附件五的所有数据满足此条件。
故频域中的滤波函数可表示为:
其函数图像如图1.
图1R-L滤波函数图像
连续的R-L卷积函数所得结果为:
离散的R-L卷积函数所得结果为:
根据上述滤波原理,在本题中,对附件五中数据的具体滤波过程可用Matlab内置的Ram-Lak命令实现。