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.