fbp算法 matlab

FBP算法是一种在计算机断层扫描（CT）成像中广泛使用的重建算法，它基于滤波反投影（FBP）原理。在MATLAB中，可以使用Radon变换和IRadon变换来实现FBP算法。

具体实现步骤如下：
1. 对原始图像进行Radon变换，得到投影数据。
2. 对投影数据进行滤波操作，通常使用Ram-Lak滤波器或Butterworth滤波器。
3. 对滤波后的数据进行反投影，得到重建图像。

以下是MATLAB代码示例：

% 读取原始图像
img = imread('lena.png');
% 进行Radon变换
theta = 0:179;
[R,xp] = radon(img,theta);
% 进行滤波操作
filter = 'ram-lak';
filtered_R = imfilter(R,fspecial('unsharp'),filter,'replicate');
% 进行反投影
reconstructed_img = iradon(filtered_R,theta,'linear','none',1,size(img,1));
% 显示重建图像
imshow(reconstructed_img,[]);