%使用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)