由投影重建图像

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');

% 36 个角度

figure,imagesc(theta2,xp,R2);xlabel('\theta');ylabel('x\prime'); % 90 个角度

figure,imagesc(theta3,xp,R3);xlabel('\theta');ylabel('x\prime');

3．用MATLAB 中的iradon 函数对获得的投影数据进行滤波反投影重建，获得Shepp‐Logan 模型的重建图像：

I1=iradon(R1,10);I2=iradon(R2,5);I3=iradon(R3,2);

%显示重建图像：

figure,imshow(I1);

figure,imshow(I2);

figure,imshow(I3);

四、 实验总结

本次实验内容较少，通过本次实验我了解了反投影重建算法的方法.了解了利用radon 和iradon 函数实现平行束投影和反投影重建算法。了解了利用fanbeam 和ifanbeam 函数实现扇形束投影和反投影重建算法。对Shepp‐Logan 模型的投影不同角度的图行有了直观的了解。