10-8.m

matlab图像处理源码

%使用radon函数和iradon函数计算图像的投影并从投影中重构图像
P = phantom(256);
%应用在X线断层摄影术里广泛使用一个测试图像:：Shepp-Logan Head影像
imshow(P)
%三种不同角度投影模式
theta1 = 0:10:170; [R1,xp] = radon(P,theta1);
%共存在18个角度投影
theta2 = 0:5:175;  [R2,xp] = radon(P,theta2);
%共存在36个角度投影
theta3 = 0:2:178;  [R3,xp] = radon(P,theta3);
%共存在90个角度投影
figure, imagesc(theta3,xp,R3); colormap(hot); colorbar
%显示Shepp-Logan Head影像的radon变换
xlabel('\theta'); ylabel('x\prime');
I1 = iradon(R1,10);
I2 = iradon(R2,5);
I3 = iradon(R3,2);
%三种情况的逆radon变换，重构图像
figure;subplot(131)
imshow(I1)
subplot(132), imshow(I2)
subplot(133), imshow(I3)