matlab卷积反投影,[转载]X-CT图像重建的卷积反投影图解法

摘要：X-CT(即X射线计算机断层成像)图像重建的卷积反投影原理对医学影像专业学生较为抽象难懂，本文在简化条件下采用图解法较为系统地描述该原理，达到直观易懂的目的。

1　X-CT图像重建的卷积反投射原理

X-CT成像原理在不少刊物上有详细介绍[1-3]。所谓CT图像重建是按照采集后的数据，求解图像矩阵中像素，然后重新构造图像的过程；而图像矩阵的求解由计算机完成，在X-CT图像重建的解析法中，当前最常用的是采用卷积运算的滤波反投影法，也称卷积反射影法。

卷积反投影重建图像时，先把由检测器上获得的原始数据与一个滤波函数进行了卷积运算，得到各方向卷积的投影函数；然后再把它们从各方向进行反投影，即按其原路径平均分配到每一矩阵元上，进行叠加后得到每一矩阵元的CT值；再经过适当处理后就可以得到被扫描物体的断层图像，卷积反投影可消除单纯的反投影产生的边缘失锐效应，补偿投影中的高频成分和降低投影中心密度，并保证重建图像边缘清晰和内部分布均匀。

卷积反投影法的数学表达式可由雷当投影定理推导出[2]；下面直接给出表达式。

如图1所示，取物体某断层面为x-y坐标系，X射线管A发出的单能窄束X射线束穿过衰减系数为μ(x,y)的物体断层，沿S路径X线束的投影可写为：

式中I0是X线管输出强度，Im是检测器测得的X线束强度。