直接反投影重建算法（matlab）

Hannah的踩坑记录

已于 2022-08-04 14:36:08 修改

阅读量4.5k

点赞数 8

分类专栏：医学图像成像学习文章标签： matlab 算法图像处理

于 2022-05-10 22:23:49 首次发布

本文链接：https://blog.csdn.net/weixin_42118607/article/details/124698610

版权

医学图像成像学习专栏收录该内容

5 篇文章 3 订阅

订阅专栏

本文介绍了医学图像重建的基本原理，通过观看国立台湾大学的成像课程视频理解投影数据的反投影过程。实验部分使用MATLAB演示了两种重建方法：iradon函数和解析法。代码示例展示了如何使用iradon函数以及自定义的medfuncBackprojection函数进行图像重建，并给出了重建效果。

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

1.原理推导

将仿真得到的投影数据，或者扫描得到的真实投影数据，沿着扫描路径按照原路回抹到路径经过的像素上，然后累加；数学公式的推导可参考YouTube上国立台湾大学成像课程https://www.youtube.com/watch?v=0255KLvu75g&t=1295s

2.实验

方法一："iradon"函数

>> clear all;
>> close all;
>> I=phantom(256);
>> theta=0:179;
>> P=radon(I,theta);
>> rec=iradon(P, theta, 'none');
>> figure;
>> imshow(I, []), title("raw image");
>> figure;
>> imshow(rec, []), title("imaging");

重建效果

方法二：解析法

>> clear all;
>> close all;
>> 
>> 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);

>> figure;
>> imshow(I, []);
>> rec=medfuncBackprojection(theta_num, N, P1, delta);
>> figure;
>> imshow(rec, []);


>> edit medfuncBackprojection;
function rec= medfuncBackprojection(theta_num, N, R1, delta)

rec=zeros(N);
for m=1:theta_num
    pm=R1(:,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