% Demonstration of radon.m
%
% Samuli Siltanen Jan 2015

% Construct a phantom
f = phantom(32);

% Compute the sinogram
sinog = radon(f,1:180);

% Compute a reconstruction
rec = iradon(sinog,1:180)

% Take a look
figure(1)
clf
subplot(1,3,1)
imagesc(f)
colormap gray
axis equal
axis off
subplot(1,3,2)
imagesc(sinog)
colormap gray
axis equal
axis off
subplot(1,3,3)
imagesc(rec)
colormap gray
axis equal
axis off
