Free Software

Tomography Software and Example

Software Download: The zip file extracts into one MATLAB m-file and 62 jpg files. The latter images are showing different perspectives of a typical wave pattern in the Belousov-Zhabotinsky (BZ) reaction. The observation angle increases by 5.8 degrees from picture to picture. In the data provided here, the wave pattern essentially did not move ( less than 5 pixels) during data aquisition.
Running the Software: Once extracted you should change the path information in the m-file and excecute it with MATLAB. We have tested the program with MATLAB versions 6.5 and 7.4. The script requires the function IRADON from the image processing toolbox. The program will compute the three-dimensional tomographic reconstruction based on the 62 images. It will open three windows illustrating the reconstruction result in terms of (i) a reconstructed isosurface, (ii) the posterior half of it, and (iii) a slice-wise visualization. On our PCs the run time of the script is approximately 2-3 minutes.

How it Works: 3D image reconstruction is carried out slice-wise with the individual slices being normal to the rotation axis of the sample. The crucial step is to compute the inverse Radon transform of the data using "filtered backprojection". A concise description of the method is given in: Winfree et al., Chaos 6, 617-626 (1996).

About the Example: We are studying three-dimensional concentration waves in the BZ reaction. These waves are traveling zones of oxidation in a reduced medium. They are detected via local changes in the system’s absorption of white light. In simple terms, we see blue waves in a red system. During the experiment, a cylindrical cuvette with the reaction system is rotated around its symmetry axis at 12 rpm using a small electric motor and a simple belt drive. This rotation allows us to collect image data from various perspectives. Image sequences of the rotating sample are recorded with a monochrome video camera connected to PC-based frame grabber. There is no fluid motion. The wave patterns arise from autocatalytic reactions "infecting" neighboring regions by diffusion. This type of behavior is called excitation wave and shares many similarities with action potentials propagating in neuronal and cardiac tissue. In this specific example, the wave patterns is a three-dimensional spiral wave. This "scroll ring" is not static but subsequent analyses (not provided here) show that it rotates around a circular loop known as a "filament".

The figures show data reconstructed from experiments with the BZ reaction.

Chiral Walkers

The MatLab script CW.m allows you to run basic examples of a chiral walker on a square lattice. Background: The model is based on rotating walkers with trailing tails spanning L lattice bonds. These tail segments cannot overlap and their leading A segments cannot be crossed. As prescribed by their chirality, walkers must turn if possible, or go straight, or else correct earlier steps recursively. The resulting motion traces unbound trajectories and complex periodic orbits with various symmetries. Periods tend to decrease with increasing L and vary between L and L². Interacting walkers can form intricate pair states. Some orbits match pinned spiral tip trajectories observed experimentally in excitable systems.

Acknowledgment: This material is based upon work supported by the National Science Foundation under Grant No. CHM-0513912. If you have comments regarding this software, please contact us at
Selected References:

1. B. Marts, T. Bánsági, Jr. and O. Steinbock, "Evidence for Burgers' Equation Describing the Untwisting of Scroll Rings", Europhys. Lett. - accepted
2. T. Bánsági, Jr. and O. Steinbock, "Three-Dimensional Spiral Waves in an Excitable Reaction System: Initiation and Dynamics of Scroll Rings and Scroll Ring Pairs", Chaos 18, 026102, 1-8, 2008 pdf
3. T. Bánsági, Jr., K. J. Meyer and O. Steinbock, "Wave-Pinned Filaments of Scroll Rings", J. Chem. Phys. 128, 094503, 1-7, 2008 pdf
4. T. Bánsági, Jr. and O. Steinbock, "Negative Filament Tension of Scroll Rings in an Excitable System", Phys. Rev. E 76, 045202, 1-4, 2007 pdf
5. T. Bánsági, Jr., C. Palczewski and O. Steinbock, "Scroll Wave Filaments Terminate in the Back of Traveling Fronts", J. Phys. Chem. A 111, 2492-2497, 2007 pdf
6. T. Bánsági Jr. and O. Steinbock, "Nucleation and Collapse of Scroll Rings in Excitable Media", Phys. Rev. Lett. 97, 198301, 1-4, 2006 pdf