Reconstruction from projections

The problem of reconstructing a multidimensional signal from its projection is uniquely multidimensional, having no 1-D counterpart. It has applications that range from computer-aided tomography to geophysical signal processing. It is a problem which can be explored from several points of view—as a deconvolution problem, a modeling problem, an estimation problem, or an interpolation problem.[1][2]

Motivation and applications

edit

The problem of reconstruction from projections has arisen independently in a large number of scientific fields, since it is widely applied in areas such as medical imaging, geophysical tomography, industrial radiography and so on. For example, by using CT scanner, the lesion information of the patients can be presented in 3D on the computer, which offers a new and accurate approach in diagnosis and thus has vital clinical value.

Problem statement and basics

edit

A projection is a linear mapping of an   dimensional signal into an   dimensional one, where  . And the objective of reconstruction is to restore the   dimensional signal based on the   dimensional signal. The following case is a 2-D signal projected into 1D signal.[3] The signal in the original coordinate is denoted as  . Now consider a collimated beam of radiation coming from the opposite orientation of  , producing a projection along  .   and   are normal to each other, and the angle between   and   is theta. The signal obtained along   axis is defined to be  . The relationship between the original coordinate and the rotated coordinate is given by

 
Collimated beam scanning geometry

 

or inversely,

 

Then we have

 

By varying theta, a large number of projections can be obtained.

Given the projection-slice theorem,  ,the slice of the Fourier transform of   at angle theta, is equivalent to  , the Fourier Transform of the projection  . Therefore, the unknown   can be obtained from its Fourier transform by means of the Fourier transform inversion integral

         

By taking the inverse Fourier Transform and assuming  , we get  

Approaches

edit

In practice, there are a wide rarity of methods that are utilized, most of which are reconstruct 3-D information (volume) from 2-D signals (image). Typically used methods are CT, MRI, PET and SPECT. And the filtered back projection based on the principles introduced above are commonly applied.[4][5]

 
Computed tomography of human brain - large

Computed Tomography (CT)

edit

In CT, a volume is formed by stacking the axial slices. The software cuts the volume in a different plane (usually orthogonal). Commonly, slice data is generated using an X-ray source that rotates around the object. X-ray sensors are positioned on the opposite side of the circle from the X-ray source.

Magnetic resonance imaging (MRI)

edit

In MRI, energy from an oscillating magnetic field is temporarily applied to the patient at the appropriate resonance frequency. The protons (hydrogen atoms) emit a radio frequency signal which is measured by a receiving coil. The radio signal can be made to encode position information by varying the main magnetic field using gradient coils.[6]

Positron emission tomography (PET)

edit

The system detects pairs of gamma rays emitted indirectly by a positron-emitting radionuclide (tracer), which is introduced into the body on a biologically active molecule. Three-dimensional images of tracer concentration within the body are then constructed by computer analysis. In modern PET-CT scanners, three dimensional imaging is often accomplished with the aid of a CT X-ray scan performed on the patient during the same session, in the same machine.

Single-photon emission computed tomography (SPECT)

edit

SPECT imaging is performed by using a gamma camera to acquire multiple 2-D images (projections) from multiple angles. Multiple projections are used to yield a 3-D data set. This data set may then be manipulated to show thin slices along any chosen axis of the body. SPECT is similar to PET in its use of radioactive tracer material and detection of gamma rays, while the tracers used in SPECT emit gamma radiation that is measured more directly.[7]

See also

edit

References

edit
  1. ^ Dudgeon, Dan E., and Russell M. Mersereau. "Multidimensional digital signal processing." Prentice-Hall Signal Processing Series, Englewood Cliffs: Prentice-Hall, 1984 1 (1984).
  2. ^ Gonnor, Ricrrann, and Gabor T. Herman. "Three-dimensional reconstruction from projections: a review of algorithms." International Review of Cytology 38 (1974): 111.
  3. ^ Dudgeon, Dan E., and Russell M. Mersereau. "Multidimensional digital signal processing." Prentice-Hall Signal Processing Series, Englewood Cliffs: Prentice-Hall, 1984 1 (1984).
  4. ^ Doi, Kunio. "Computer-aided diagnosis in medical imaging: historical review, current status and future potential." Computerized medical imaging and graphics 31.4 (2007): 198-211.
  5. ^ Ledley, Robert S., and William R. Ayers. "Computerized Medical Imaging and Graphics." (1995).
  6. ^ Schultz, Gerrit, et al. "MR image reconstruction from generalized projections." Magnetic Resonance in Medicine 72.2 (2014): 546-557.
  7. ^ Budinger, T. F. "Single photon emission computed tomography." Diagnostic nuclear medicine. Second edition. 1988.
edit