The following codes are written in Scilab, using the methods explained in chapter 3 of Kak and Slaney (1988):

phantom_2a.sci - generates a Shepp Logan phantom and the corresponding sinogram.
recon_3a.sci - reconstructs the image using the Convolution Back Projection method.
test_script.sce - a script to test the above codes

Explanations on how to use these are are found in the codes.

The function phantom_2.sci outputs two data files:   phantom.txt,   sinogram.txt
and the function recon_3a.sci outputs one:               reconstruction.txt

These can respectively be shown as images using the following Gnuplot scripts:

show_phantom.plt
show_sinogram.plt
show_reconstruction.plt

The results are as shown here:


A Shepp Logan phantom


The sinogram, using 110 angles and 127 rays for each projection.


Reconstructed using the Convolution Back Projection method.

A Matlab version of the codes can be found here: Link. The background of the reconstructed image here looks brighter than the image from Matlab. The reason is that Gnuplot uses a slightly difference grayscale.

Last modified 10 June 2010.