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