基于MATLAB的CT直接反投影算法的实验研究

看星河的兔子

已于 2023-04-30 14:35:27 修改

阅读量628

点赞数 2

文章标签： matlab 算法开发语言

于 2023-04-30 14:34:45 首次发布

本文链接：https://blog.csdn.net/fcxgfdjy/article/details/130448385

版权

该篇博客介绍了如何使用MATLAB的iradon函数进行CT平行束投影数据的直接反投影重建，并对比了使用内置函数与自编算法的重建结果。作者首先展示了iradon函数的简单应用，然后详细推导并实现了自定义的直接反投影算法，最后通过图像显示比较了两种方法的重建图像。

摘要由CSDN通过智能技术生成

推导参考我的上一篇博客基于MATLAB的CT平行束投影数据的仿真
实验任务
1、利用“iradon”函数进行直接反投影重建。

I=phantom(256);
theta=0:179;
P=radon(I,theta);
rec=iradon(P,theta,'None');
figure;
imshow(I,[]),title('原始图像');
figure;
imshow(I,[]),title('直接反投影图像');

在这里插入图片描述
2、采用本节推导出来的直接反投影算法进行重建。

clear;
N=256;
I=phantom(N);
delta=pi/180;
theta=0:1:179;
theta_num=length(theta);
%%产生投影数据
P=radon(I,theta);
[mm,nn]=size(P);
e=floor((mm-N-1)/2+1)+1;
P=P(e:N+e-1,:);
P1=reshape(P,N,theta_num);
%%投影重建
rec=zeros(N);
for m=1:theta_num
    pm=P1(:,m);
    Cm=(N/2)*(1-cos((m-1)*delta)-sin((m-1)*delta));
    for k1=1:N
        for k2=1:N
            Xrm=Cm+(k2-1)*cos((m-1)*delta)+(k1-1)*sin((m-1)*delta);
            n=floor(Xrm);
            t=Xrm-floor(Xrm);
            n=max(1,n);
            n=min(n,N-1);
            p=(1-t)*pm(n)+t*pm(n+1);
            rec(N+1-k1,k2)=rec(N+1-k1,k2)+p;
        end
    end
end
figure;
imshow(I,[]),title('原始图像');
figure;
imshow(rec,[]),title('解析法直接反投影图像');