用matlab对一个已知的矩阵进行后向投影成像的代码

以下是MATLAB代码，实现了对已知矩阵进行后向投影成像：

% 已知矩阵
img = imread('example.jpg');
theta = 0:179; % 投影角度范围
[R,xp] = radon(img,theta); % Radon变换

% 后向投影
backproj = zeros(size(img));
for i = 1:length(theta)
    angle = theta(i);
    proj = R(:,i);
    proj = repmat(proj,1,size(img,2));
    proj = imrotate(proj,-angle,'crop');
    backproj = backproj + proj;
end

% 显示结果
figure;
subplot(1,2,1);imshow(img);title('原图');
subplot(1,2,2);imshow(backproj,[]);title('重建图像');

首先，我们使用`imread`加载一个示例图像，然后定义投影角度范围`theta`。接下来，使用`radon`函数对图像进行Radon变换，得到投影矩阵`R`和对应的投影坐标`xp`。

然后，我们使用循环遍历每个投影角度，将投影矩阵重复到与原图像相同大小，然后使用`imrotate`函数将投影矩阵旋转回原来的角度，并将所有投影矩阵相加，得到后向投影图像。

最后，使用`imshow`函数显示原图像和重建图像的结果。注意，在显示重建图像时，我们将像素值限制在`[]`范围内，以使亮度比较均匀。