//------------------------
//
//  Tomography test codes:
//
//  To generate phantom image, sinogram, and
//  to do reconstruction using the
//  Convolution Back Projection method.
//
//  Kai Hock
//  Cockcroft Institute
//  10 June 2010
//------------------------
clear

//getf('heaviside.sci') ;
getf('phantom_2a.sci') ;
getf('recon_3a.sci') ;

rho1  = -1;     // lower limit of position in each projection
rho2  = 1;      // upper limit of position in each projection
nproj = 110;    // number of projections (or angles)
npos  = 127;    // number of rays (or positions)

[sinogram, angles, positions, centre, weights] = phantom_2a(rho1, rho2, nproj, npos) ;

x0 = rho1:.01:rho2;
y0 = x0;

recon_3a(x0, y0, sinogram, angles, positions, centre, weights, nproj, npos);