radon1.asv

Matlab code to reconstruct using back projection

clear;
P=phantom('Modified Shepp-Logan',512);
output_size=max(size(P));

figure,imshow(P),colormap(gray),colorbar,title('Original Image');

theta1=0:10:170;
[R1,xp]=radon(P,theta1);
num_angle_R1=size(R1,2)
figure,imshow(theta1,xp,R1,[],'notruesize'),colormap(gray),colorbar,title('Radon Transform 1');

dtheta1=theta1(2) - theta1(1);
I1=iradon(R1,dtheta1,output_size);
figure,imshow(I1),colormap(gray),colorbar,title('Reconstructed Image 1');

theta2=0:5:175;
[R2,xp]=radon(P,theta2);
num_angle_R2=size(R2,2)
figure,imshow(theta2,xp,R2,[],'notruesize'),colormap(gray),colorbar,title('Radon Transform 2');

dtheta2=theta2(2) - theta2(1);
I2=iradon(R2,dtheta2,output_size);
figure,imshow(I2),colormap(gray),colorbar,title('Reconstructed Image 2');

theta3=0:2:178;
[R3,xp]=radon(P,theta3);
num_angle_R3=size(R3,2)
figure,imshow(theta3,xp,R3,[],'notruesize'),colormap(gray),colorbar,title(' Radon transform 3');

dtheta3=theta3(2) - theta3(1);
I3=iradon(R3,dtheta3,output_size);
figure,imshow(I3),colormap(gray),colorbar,title('Reconstructed Image 3');