n

theta = linspace(0,360,900);
%bwelow we have to multiply the projections by the voxel size or else each
%voxel is given a "unit" wieght. This means for the same image simply
%changing the matrix size would change the projections values which is
%wrong, so we must account for this.
projs = radon(imin,theta)*voxelSize;
%plot the projections
figure;imagesc(projs,[0 max(max(projs))]);colormap('gray');axis off;axis equal;title('Original Sinogram');