fanbeamreconst.m

Matlab code to reconstruct using back projection

% Image Reconstruction from Fanbeam Projections algorithm
% Varie the D from 1 to 360
clear;
ph = phantom(256);
figure,subplot(121),imshow(ph),title('Original Image');
[F,Fpos,Fangles] = fanbeam(ph,250);
 
subplot(122),imshow(F,[],'XData',Fangles,'YData',Fpos,'InitialMagnification','fit');
axis normal
xlabel('Rotation Angles (degrees)')
ylabel('Sensor Positions (degrees)')

D=180;
I = ifanbeam(F,D);
figure,imshow(I),title('Reconstructed Image');