# Tomographic reconstruction

Tomographic imaging is applied in Computed Tomography to obtain cross-sectional images of patients. This article applies in general to tomographic reconstruction for all kinds of tomography, but some of the terms and physical descriptions refer directly to X-ray computed tomography. The mathematical basis for tomographic imaging was laid down by Johann Radon.

## Description

Figure 1: Parallel beam geometry utilized in tomography and tomographic reconstruction. Each projection, resulting from tomography under a specific angle, is made up of the set of line integrals through the object.
Resulting tomographic image from a plastic skull phantom. Projected X-rays are clearly visible on this slice taken with a CT-scan as image artifacts, due to limited amount of projection slices over angles.

The projection of an object, resulting from the tomographic measurement process at a given angle ${\displaystyle \theta }$, is made up of a set of line integrals (see Fig. 1). A set of many such projections under different angles organized in 2D is called sinogram (see Fig. 3). In X-ray CT, the line integral represents the total attenuation of the beam of x-rays as it travels in a straight line through the object. As mentioned above, the resulting image is a 2D (or 3D) model of the attenuation coefficient. That is, we wish to find the image ${\displaystyle \mu (x,y)}$. The simplest and easiest way to visualise the method of scanning is the system of parallel projection, as used in the first scanners. For this discussion we consider the data to be collected as a series of parallel rays, at position ${\displaystyle r}$, across a projection at angle ${\displaystyle \theta }$. This is repeated for various angles. Attenuation occurs exponentially in tissue:

${\displaystyle I=I_{0}\exp \left({-\int \mu (x,y)\,ds}\right)}$

where ${\displaystyle \mu (x)}$ is the attenuation coefficient at position ${\displaystyle x}$ along the ray path. Therefore generally the total attenuation ${\displaystyle p}$ of a ray at position ${\displaystyle r}$, on the projection at angle ${\displaystyle \theta }$, is given by the line integral:

${\displaystyle p(r,\theta )=\ln(I/I_{0})=-\int \mu (x,y)\,ds}$

Using the coordinate system of Figure 1, the value of ${\displaystyle r}$ onto which the point ${\displaystyle (x,y)}$ will be projected at angle ${\displaystyle \theta }$ is given by:

${\displaystyle x\cos \theta +y\sin \theta =r\ }$

So the equation above can be rewritten as

${\displaystyle p(r,\theta )=\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }f(x,y)\delta (x\cos \theta +y\sin \theta -r)\,dx\,dy}$

where ${\displaystyle f(x,y)}$ represents ${\displaystyle \mu (x,y)}$. This function is known as the Radon transform (or sinogram) of the 2D object. The projection-slice theorem tells us that if we had an infinite number of one-dimensional projections of an object taken at an infinite number of angles, we could perfectly reconstruct the original object, ${\displaystyle f(x,y)}$. So to get ${\displaystyle f(x,y)}$ back, from the above equation means finding the inverse Radon transform. It is possible to find an explicit formula for the inverse Radon transform. However, the inverse Radon transform proves to be extremely unstable with respect to noisy data.

## Usage and reconstruction algorithms

In practice of tomographic image reconstruction, often a stabilized and discretized version of the inverse Radon transform is used, known as the filtered back projection algorithm. Recent developments have seen the Radon transform and its inverse used for tasks related to realistic object insertion required for testing and evaluating Computed Tomography use in Airport Security.[1] An alternative family of tomographic reconstruction algorithms are the Algebraic Reconstruction Technique ones.

## Tomographic reconstruction software

For flexible tomographic reconstruction, open source toolboxes are available, such as the ASTRA toolbox .[2] The ASTRA Toolbox is a MATLAB toolbox of high-performance GPU primitives for 2D and 3D tomography, from 2009-2014 developed by iMinds-Vision Lab, University of Antwerp and since 2014 jointly developed by iMinds-VisionLab, UAntwerpen and CWI, Amsterdam. The toolbox supports parallel, fan, and cone beam, with highly flexible source/detector positioning. A large number of reconstruction algorithms are available, including FBP, SIRT, SART, CGLS.

## Gallery

Shown in the gallery is the complete process for a simple object tomography and the following tomographic reconstruction based on ART.

## References

1. ^ Megherbi, N., Breckon, T.P., Flitton, G.T., Mouton, A. (October 2013). "Radon Transform based Metal Artefacts Generation in 3D Threat Image Projection". Proc. SPIE Optics and Photonics for Counterterrorism, Crime Fighting and Defence (PDF). 8901. SPIE. pp. 1–7. doi:10.1117/12.2028506. Retrieved 5 November 2013.
2. ^ Van Aarle, W., Palenstijn, W.J., De Beenhouwer, J., Altantzis T., Bals S., Batenburg K. J., and J. Sijbers (October 2015). "The ASTRA Toolbox: a platform for advanced algorithm development in electron tomography". Ultramicroscopy. 157: 35–47. doi:10.1016/j.ultramic.2015.05.002.