MATLAB医学CT图像滤波反投影重建系统

资 源 简 介

该项目基于MATLAB实现了医学CT投影数据的滤波反投影重建算法，包含Radon变换、滤波函数设计和反投影重建三大模块。支持平行束和扇形束扫描几何，能够从原始投影数据重建出高质量的横截面医学图像，适用于医学影像处理研究。

详 情 说 明

基于MATLAB的医学CT图像滤波反投影重建系统

项目介绍

本项目实现医学CT投影数据的滤波反投影图像重建算法，主要用于从原始CT扫描投影数据重建出高质量的横截面图像。系统支持Radon变换、滤波函数设计和反投影重建三大核心模块，可处理平行束和扇形束扫描几何结构，并提供多种滤波函数选择以获得不同特性的重建图像。

功能特性

• 核心算法实现：完整实现Radon变换与逆Radon变换算法、频域滤波函数设计与卷积处理、反投影累加重建算法
• 多扫描几何支持：支持平行束和扇形束两种扫描几何结构
• 可配置滤波函数：提供Ram-Lak、Shepp-Logan、Cosine等多种滤波函数选择
• 重建质量评估：内置均方误差（MSE）和峰值信噪比（PSNR）量化评估指标
• 可视化功能：实时显示滤波后的投影数据和反投影累积过程动画
• 灵活参数配置：可自定义投影角度范围、探测器间距等扫描参数

使用方法

1. 准备输入数据：载入M×N二维投影数据矩阵（sinogram数据），其中M为投影角度数，N为探测器单元数

1. 配置重建参数：

- 设置投影角度范围（默认0:179度） - 配置探测器间距参数 - 选择扫描几何类型（平行束/扇形束） - 指定滤波函数类型

1. 执行图像重建：运行主程序启动重建流程，系统将自动完成滤波反投影计算

1. 查看重建结果：

- 获取K×K二维灰度重建图像 - 查看质量评估指标（MSE、PSNR） - 观察中间过程可视化效果

系统要求

• 软件环境：MATLAB R2018a或更高版本
• 必要工具箱：Image Processing Toolbox
• 内存要求：建议4GB以上RAM（根据数据规模调整）

文件说明

主程序文件集成了一整套完整的CT图像重建流程，其核心能力包括：系统参数初始化与用户交互界面管理、投影数据的预处理与规范化操作、扫描几何结构的自动识别与相应算法切换、频域滤波器的动态设计与卷积滤波处理、反投影重建算法的执行与图像生成、重建质量的量化评估计算，以及中间过程与最终结果的可视化展示功能。