trial1.m

Matlab code to reconstruct using simple back projection

% This program Reconstruct the image from its projections                  %
% The following command will genrate and display "Phantom                  %
% head model" of size 256, 256                                             %
P = phantom(256);
imshow(P)
% To generate projections,"theta" is varied from 0 to 180                  %
% with different incremental steps. and "Radon" function is called         %
theta1 = 0:5:175; [R1,xp] = radon(P,theta1);
figure, imagesc(theta1,xp,R1); colormap(gray); colorbar
xlabel('\theta'); ylabel('x\prime');
% Inverse Radon function is called to Reconstruct the image                %
I1 = iradon(R1,theta1);
figure, imshow(I1)