Fork me on GitHub

TomograPy by nbarbey

Siddon tomography algorithm in Python

What is TomograPy ?

This is a fast parallelized tomography projector / backprojector. It originates from solar tomography application but could be use for other applications. It uses the fast Siddon algorithm as its core. The parallelization is done with OpenMP on the C part of the code.

Application in solar tomography

For solar astrophysics, TomograPy allows to perform solar tomographic reconstructions of the corona. It can take SOHO or STEREO data and output a 3-dimensional map of the corona. For now only Extreme-Ultraviolet data is handled but handling of coronographic data is planed along with handling of data from other spacecrafts. For more information, you can take a look at the Solar Physics TomograPy paper [arxiv].

Requirements

You need numpy, pyfits, fitsarray for TomograPy to run.

Documentation

The documentation for TomograPy is embedded into the code as docstrings. To get started, take a look at TomograPy/__init__.py or in IPython, just do:
>>> import TomograPy
>>> TomograPy?

Exemple

For an example on how to use TomograPy, see exemple/test_siddon_simu.py For solar tomography, you need first to process data using http://www.lmsal.com/solarsoft/. Then you can take a look at exemple/siddon_secchi_mask.py or use directly the TomograPy/srt script.

Dependencies

numpy, pyfits, fitsarray, lo

Install

You should have distutils installed. Then, just do:
easy_install TomograPy
Alternatively, you can install using the development sources with:
git clone git://github.com/nbarbey/TomograPy.git
cd TomograPy
python setup.py install

License

CeCILL

Authors

Nicolas Barbey (nicolas.a.barbey@gmail.com)

Contact

nicolas.a.barbey@gmail.com

Download

You can download this project in either zip or tar formats.

You can also clone the project with Git by running:

$ git clone git://github.com/nbarbey/TomograPy