ct扇束等角滤波反投影重建算法matlab源码

资 源 简 介

ct扇束等角滤波反投影重建算法matlab源码

详 情 说 明

CT扇束等角滤波反投影重建是医学成像领域的经典算法，其实现通常可分为三个核心模块：源图像创建、投影计算以及滤波反投影重建。三个文件协同完成从数字体模到重建图像的完整流程。

源图像创建模块负责构建测试用的数字体模。该模块通过定义不同形状和密度的区域来模拟人体组织的衰减系数分布，常见做法包括椭圆模型或自定义灰度矩阵。图像矩阵的每个像素值代表该位置的线性衰减系数，这是后续投影计算的基础数据。

投影计算模块实现扇束几何下的射线追踪过程。该模块需要处理X射线源围绕物体旋转时各个角度的衰减积分。核心是计算每条射线经过图像矩阵时的路径长度加权求和，输出为各个角度下的投影数据（正弦图）。扇束几何需特别注意探测器单元与射线源之间的角度关系。

滤波反投影模块包含两个关键步骤：首先对投影数据进行斜坡滤波（常用Ram-Lak滤波器），消除简单反投影导致的模糊效应；随后将滤波后的投影数据按对应角度反向投射到图像空间并累加。等角扇束重建需要特殊的加权处理来校正射线密度不均匀性。重建质量受投影角度数量、滤波函数选择等因素影响。