滤波反投影算法matlab

滤波反投影算法（Filtered Back Projection, FBP）是计算机断层扫描（Computed Tomography, CT）图像重建中最基本的算法之一。在Matlab中，可以使用以下代码实现滤波反投影算法：

% 读入投影数据
data = load('projection_data.txt'); 

% 设置投影角度
theta = 0 : 1 : 179;

% 设置重建图像大小
N = 256;

% 设置滤波器类型
filter_type = 'Ram-Lak';

% 生成滤波器
h = fspecial('unsharp');

% 进行滤波和反投影重建
recon = iradon(data, theta, filter_type, N, 'none', 1, h);

% 显示重建图像
imshow(recon, []);

其中，`projection_data.txt`是投影数据文件，`theta`是投影角度，`N`是重建图像大小，`filter_type`是滤波器类型，可以选择`Ram-Lak`、`Shepp-Logan`或`Cosine`等，`h`是滤波器，可以根据需要自己生成。最后，使用`iradon`函数进行滤波和反投影重建，并使用`imshow`函数显示重建图像。