基于FBP算法的CT图像重建仿真系统

资 源 简 介

本系统旨在实现高质量的CT图像重建，其核心采用经典的滤波反投影（Filtered Back Projection, FBP）算法。项目主要针对标准头部仿真模型（如Shepp-Logan模型）进行投影数据的模拟采集与解析重建。具体功能涵盖了从原始投影数据（Sinogram）的生成、R-L（Ram-Lak）滤波器或Shepp-Logan滤波器的设计与频域应用、到反投影算子的离散数值实现全过程。系统支持不同探测器通道数量和旋转步长的参数化配置，能够精确还原模型内部的细微组织结构。除了基础重建功能外，该项目还集成了对重建质量的定量分析模块，通过计算均方误差（MSE）、峰值信噪比（PSNR）以及结构相似性（SSIM）等关键指标评估算法表现。本程序旨在为医学影像处理领域的学生和研究人员提供一个可观测、可扩展的学习研究平台，支持用户根据需求自定义滤波器参数以探讨伪影消除及图像对比度增强的优化策略。通过该系统，用户可以深入理解Radon变换原理、傅里叶切片定理及其在实际计算机断层扫描（CT）成像中的工程实现，是一个集算法仿真模拟、可视化结果展示与性能科学评价于一体的综合性研究工具。

详 情 说 明

基于滤波反投影(FBP)算法的头部仿真模型CT图像重建系统

项目介绍

本系统是一个专注于医学影像重建技术的仿真平台，核心旨在通过数值计算模拟计算机断层扫描（CT）的完整物理与数学过程。系统基于经典的滤波反投影（Filtered Back Projection, FBP）算法，实现了从解析头部仿真模型设计、投影数据采集（Sinogram生成）到图像恢复重建的全流程。通过该系统，用户可以观察不同滤波器对图像质量的影响，并利用量化的客观指标评估重建效果，是理解Radon变换、傅里叶切片定理及离散反投影算子的理想研究工具。

功能特性

1. 参数化模型生成：内置标准Shepp-Logan头模生成算子，支持自定义组织形状、亮度及几何参数。
2. 投影数据仿真：模拟平行束CT扫描，支持自定义旋转步长与探测器通道数量。
3. 频域滤波引擎：集成Ram-Lak（坡度滤波器）与Shepp-Logan滤波器，有效消除反投影带来的星状伪影。
4. 高精度反投影：基于几何映射与线性插值的离散反投影实现，确保投影值准确回填至像素空间。
5. 质量定量评估：集成了均方误差(MSE)、峰值信噪比(PSNR)以及结构相似性(SSIM)三大关键评价指标。
6. 多维度可视化：一键生成包含原始图像、正弦图、滤波器响应、重建图像、残差图及评价报表的综合看板。

使用方法

1. 启动环境：在MATLAB开发环境中打开项目。
2. 配置参数：根据需求修改初始配置设置，如图像分辨率（默认256）、扫描角度（默认0-179度）或滤波器类型。
3. 执行计算：运行主程序，系统将依次执行幻影生成、正弦图仿真、频域滤波、图像重建及指标计算。
4. 结果观测：程序自动弹出可视化窗口，展示重建后的解剖结构及对应性能数据报表。

系统要求

1. 软件环境：MATLAB R2016b 或更高版本。
2. 硬件要求：建议内存4GB以上以确保在处理高分辨率（512x512及以上）投影时的运算速度。
3. 依赖库：主要使用MATLAB内置数值计算库，包括信号处理相关函数及图像展示模块。

详细实现逻辑与算法细节

1. 仿真幻影生成

系统通过定义的椭圆参数矩阵构建Shepp-Logan头部模型。每个椭圆由亮度、长短轴长度、中心坐标及旋转角度六个参数驱动。在笛卡尔坐标系下，通过判断像素点是否在旋转平移后的椭圆方程范围内，动态累加亮度值，最终形成反映人体头部组织的灰度图像。

2. 正弦图(Sinogram)仿真采集

为了模拟CT扫描过程，系统对原始图像进行离散角度的投影运算。

• 旋转映射：利用双线性插值算法对图像进行旋转。
• 积分投影：对旋转后的图像在垂直方向进行累加，模拟探测器接收到的X射线能量衰减。
• 坐标对齐：根据探测器通道数量，对投影数据进行几何中心对齐，确保旋转中心与图像物理中心重合。

3. 频域滤波器设计

滤波是解决反投影“边缘模糊”现象的核心步骤。系统执行以下过程：

• 补零扩充：为防止时域卷积中的环绕混叠，将投影数据长度通过补零扩展至2的幂次。
• 滤波器构造：生成归一化的Ram-Lak滤波器（绝对值函数）；若选择Shepp-Logan滤波器，则在此基础上叠加sinc函数窗口，以抑制高频噪声。
• 频域处理：将投影信号经由FFT转换至频域，与所选滤波器响应相乘后，通过IFFT还原至时域。

4. 离散反投影重建

反投影过程旨在将滤波后的投影值按原扫描路径回填至图像矩阵。

• 坐标映射：利用像素点的坐标(X, Y)与当前旋转角度，根据傅里叶切片定理计算该点在探测器上的投影位置。
• 线性插值：由于计算出的投影位置通常非整数，系统采用离散插值技术在探测器通道间提取精确的投影强度。
• 幅度修正：对所有角度的反投影图像进行累加，并叠加权重系数修正由于极坐标变换带来的幅度增益。

5. 图像质量评估算法

• 均方误差(MSE)：计算原始图像与重建图像之间像素强度的平均平方差异。
• 峰值信噪比(PSNR)：基于最大像素值与MSE的关系，衡量图像的保真度，单位为dB。
• 结构相似性(SSIM)：通过对比原始图像与重建图像的亮度平均值、方差及协方差，综合评价图像在亮度、对比度和结构上的相似程度，结果范围在-1到1之间。

6. 结果展示逻辑

系统利用图形化界面展示重建流水线的关键环节。除直观的图像对比外，残差图（Error Map）能清晰暴露重建过程中的几何畸变或由于滤波不足留下的伪影。生成的报表模块提供了自动化的结论建议，若SSIM指标超过0.9，则判定为高质量重建。