% Shows Sinogram explicity
N = 256;
x0 = 192; 
y0 = 192;
    im = zeros(N);
    im(x0,y0) = 1;
figure(1);clf;
subplot(1,2,1); imagesc(im);title('Original Image');
%%

ang = [0:511];

g = radon(im,ang);
for t = 1:length(ang)
figure(1);clf;
subplot(1,2,1); imagesc(im);title('Original Image');
subplot(1,2,2); imagesc(g(:,1:t)); title('Sinogram')
%pause(0.1);
end

%% now show moving pixels and sinograms

for x0 = 16:8:240
    for y0 = 16:8:240
        im = zeros(N);
        im(x0,y0) = 1;      
        g = radon(im,ang);
        
        figure(1);clf;
        subplot(1,2,1); imagesc(im);title('Original Image');
        subplot(1,2,2); imagesc(g); title('Sinogram')
        pause(0.1);
    end
end