由投影重建图像
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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 模型的投影不同角度的图行有了直观的了解。