% test_script.m
%
%  Testing of phantom and reconstruction codes
%
%
%  Kai Hock
%  Cockcroft Institute
%  10 June 2010
%---------------------------------------
clear

close all

rho1 = -1;
rho2 =  1;
nproj = 110;
npos  = 127;

%----- phantom and sinogram ---------------------

PNG = 1;   % to save figure as PNG image
EPS = 0;   % to save figure as EPS image
MAT = 1;   % to save data as MAT file

ProjectFile = 'sino';   %  MAT file name

phantom_2(ProjectFile, rho1, rho2, nproj, npos, PNG, EPS, MAT);

%r----- reconstruction ------------------------

x0 = -1:.01:1;
y0 = -1:.01:1;

PNG = 1;   % to save figure as PNG image
EPS = 0;   % to save figure as EPS image
MAT = 0;   % to save data as MAT file

ReconFile = 'recon';   %  MAT file name

recon_3a(ProjectFile, ReconFile, x0, y0, PNG, EPS, MAT);