matlab从投影重建图像,由投影重建图像

4-5 由投影重建图像

一、 实验目的

了解反投影重建算法的方法.

二、 实验内容

1.利用radon 和iradon 函数实现平行束投影和反投影重建算法

2. 利用fanbeam 和ifanbeam 函数实现扇形束投影和反投影重建算法

三、 实验步骤

1.用MATLAB图像处理工具箱的phantom 生成Shep‐Logan 头模型；

P=phantom(256);

imshow(P);

2.用MATLAB中的radon 函数获得Shepp‐Logan 模型的投影数据：

theta1=0:10:170;

[R1,xp]=radon(P,theta1); %计算Shep‐Logen头模型18 个角度

theta2=0:5:175; [R2,xp]=r adon(P,theta2); %36 个角度

theta3=0:2:178;[R3,xp]=radon(P,theta3); % 90 个角度的投应

%显示投影数据：

%18 个角度

figure,imagesc(theta1,xp,R1);xlabel('\theta');ylabel('x\prime');