British Inverse Problems Society
One-day Workshop on Inverse Problems supported by the LMS
University of Manchester
Tuesday 12th June 2007
See the programme
Truncated Projections and Region-of-Interest Reconstruction
in Classical X-Ray Computed Tomography
Rolf Clackdoyle
Abstract
The theory behind classical (2-dimensional) Computed
Tomography has been well-understood for several decades now,
with applications in many industrial and medical arenas. Direct
(non-iterative) image reconstruction methods are usually based
on a Fourier transform approach, and the well-known filtered
backprojection (FBP) algorithm has been a standard for many
years. One important requirement is that to reconstruct a
cross-section of the sample (the patient, in the medical
context), x-rays must be passed through every point in the
sample from every direction. If the sample is too large for
the detector, the projection is 'truncated' and it has long
been known that it is not possible to reconstruct the entire
sample without strong supplemental a priori information. (In
many cases a reasonable image can be obtained, but quantitative
recovery of local densities is in principle impossible.)
The particular case of some internal region being illuminated
from all directions has also been thoroughly explored - in
the mathematical literature it is known as the 'interior
problem.' Perhaps surprisingly, no part of this internal
region can be accurately reconstructed in general. However,
when examining the FBP algorithm closely one notes that the
reconstruction at a specific point uses _all_ the x-ray data,
even from rays that do not pass through that point. So, until
the turn of the century, the understanding was that classical
computed tomography was 'all or nothing' - either all rays
through the sample are measured, in which case the sample
can be reconstructed; or some rays are missing (typically,
some projections are truncated) in which case nothing can
be recovered.
In this talk, it will be shown that even with truncated
projections, certain regions of interest (ROIs) can be
accurately recovered. Over the last 5 years, new mathematical
results have shown that (nearly) any ROI that is illuminated
from all angles can be reconstructed provided it is not
internal. For some cases, relatively simple new algorithms
have been devised that verify the theory and can be readily
implemented for direct image reconstruction. It turns out that
the Hilbert transform is the relevant operator, and not the
Fourier transform.
.
The speaker
Rolf Clackdoyle is an internationally recognised expert in medical
image reconstruction from projections, for the cases of fully 3D
reconstruction for positron emission tomography (PET), single photon
emission computed tomography (SPECT) as well as fan-beam and
cone-beam tomography (CT). His current research interests include
reconstructing regions of interest from truncated fan-beam projection
data. Rolf was based at the University of Utah for a
number of years, but recently moved to France and is now research director at the University of St Etienne.
References
[1] R.Clackdoyle, F.Noo, G.Guo, J.Roberts. 2004. "Quantitative reconstruction
from truncated projections in classical tomography." IEEE Trans. Nucl. Sci.
51(3), 2570-2578. DOI 10.1109/TNS.2004.835781
[2] F.Noo, M.Defrise, R.Clackdoyle, H.Kudo. 2002. "Image reconstruction from fan-beam
projections on less than a short-scan." Phys. Med. Biol. 47, 2525-2546.
[3] R.Clackdoyle, M.Defrise, F.Noo, H.Kudo. 2006. "Two-Dimensional Region-of-Interest
Tomography." Oberwolfach Reports, G.-M.Greuel, S.Klaus (eds.) ISSN 1660-8933 Vol 3
No 3 Pages: 2070-2074.
[4] M.Defrise, F.Noo, R.Clackdoyle, H.Kudo. 2006. "Truncated Hilbert transform and
image reconstruction from limited tomographic data." Inverse Problems 22, 1037-1053.