Tomographic reconstruction
You don't need to be Editor-In-Chief to add or edit content to WikiDoc. You can begin to add to or edit text on this WikiDoc page by clicking on the edit button at the top of this page. Next enter or edit the information that you would like to appear here. Once you are done editing, scroll down and click the Save page button at the bottom of the page.
The mathematical basis for tomographic imaging was laid down by Johann Radon. It 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 projection of an object at a given angle, θ is made up of a set of line integrals. 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 μ(x,y). The simplest and easiest way to visualise 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 r, across a projection at angle θ. This is repeated for various angles. Attenuation occurs exponentially in tissue:
where μ(x) is the attenuation coefficient at position x along the ray path. Therefore generally the total attenuation of a ray at position r, on the projection at angle θ, is given by the line integral:
Using the coordinate system of Figure 1, the value of r onto which the point (x,y) will be projected at angle θ is given by:
So the equation above can be rewritten as
where f(x,y) represents μ(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, f(x,y). So to get 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. In practice, a stabilized and discretized version of the inverse Radon transform is used, known as the filtered back projection algorithm.
Further reading
- Avinash Kak & Malcolm Slaney (1988), Principles of Computerized Tomographic Imaging, IEEE Press, ISBN 0-87942-198-3.
External links
Acknowledgement and Attribution Regarding Sources of Content
Some of the initial content on this page may be incorporated in part from copyleft sources in the public domain including wikis such as Wikipedia and AskDrWiki. Drug information for patients came from the The National Library of Medicine. Infectious disease information may have come from the Centers for Disease Control (CDC). Differential Diagnoses are drawn from clinicians as well as an amalgamation of 3 sources: 1.The Disease Database; 2. Kahan, Scott, Smith, Ellen G. In A Page: Signs and Symptoms. Malden, Massachusetts: Blackwell Publishing, 2004:3; 3. Sailer, Christian, Wasner, Susanne. Differential Diagnosis Pocket. Hermosa Beach, CA: Borm Bruckmeir Publishing LLC, 2002:7 .
