This repository contains different iterative reconstruction algorithms for parallel beam tomography.
The following algorithms are included:
-
the Simultaneous Iterative Reconstruction Technique (SIRT);
-
the Maximum-Likelihood Expectation-Maximization (MLEM);
-
the Separable Paraboloidal Surrogate (SPS);
-
the Alternate Directions Method of Multipliers (ADMM).
Basic compilers like gcc and g++ and the FFTW library are required. The simplest way to use the code is with an Anaconda environment equipped with python-2.7, scipy, scikit-image and Cython.
Procedure:
-
Create the Anaconda environment (if not created yet):
conda create -n iter-rec python=2.7 anaconda
. -
Install required Python packages:
conda install -n iter-rec scipy scikit-image Cython
. -
Activate the environment:
source activate iter-rec
. -
git clone git@github.com:arcaduf/iterative_tomographic_reconstruction_algorithms.git
. -
Install routines in C:
python setup.py
.
If setup.py
runs without giving any error all subroutines in C have been installed and
your python version meets all dependencies.
If you run python setup.py 1
(you can use any other character than 1), the
all executables, temporary and build folders are deleted, the test data are
placed in .zip files. In this way, the repository is restored to its original
status, right after the download.
Go inside the folder "scripts/" and run: python run_all_tests.py