Quantitative phase-contrast tomography using polychromatic radiation

We discuss theoretical, experimental and numerical aspects of several new techniques for quantitative phase-contrast tomography using, for example, unfiltered radiation from a polychromatic X-ray microfocus source. The proposed algorithms allow one to reconstruct the three-dimensional distribution of complex refractive index in a sample consisting of one or more constituent materials, given one or more projection images per view angle. If the sample is weakly absorbing or consists predominantly of a single material, these reconstruction algorithms can be simplified and fewer projections may be required for an unambiguous quantitative reconstruction of the spatial distribution of the refractive index. In the case of weakly absorbing samples, the reconstruction algorithm is shown to be achromatic and stable with respect to high-spatial-frequency noise, in contrast to conventional tomography. A variation of the algorithm exploits the natural combination of binary tomography with a phase-retrieval method that makes explicit use of the single-material nature of the sample. Such consistent use of a priori knowledge dramatically reduces the number of required projections, implying significantly reduced dose and scanning time when compared to most alternative phase-contrast tomography methods. Experimental demonstrations are also given, using data from a point-projection X-ray microscope. The refractive index distribution, in test samples of both a polymer fibre scaffold and an adult mouse, is accurately reconstructed from polychromatic phase-contrast data. Applications of the new techniques to rapid non-destructive testing in materials science and biomedical imaging are considered.

[1]  S. Wilkins,et al.  Phase-and-amplitude computer tomography , 2006 .

[2]  O. Bunk,et al.  Hard x-ray phase tomography with low-brilliance sources. , 2007, Physical review letters.

[3]  Gabor T. Herman,et al.  A Discrete Tomography Algorithm for Improving the Quality of 3DXRD Grain Maps , 2006 .

[4]  Françoise Peyrin,et al.  Observation of microstructure and damage in materials by phase sensitive radiography and tomography , 1997 .

[5]  Atsushi Momose,et al.  Phase-contrast x-ray computed tomography for observing biological specimens and organic materials , 1995 .

[6]  T E Gureyev,et al.  Phase-contrast tomography of single-material objects from few projections. , 2008, Optics express.

[7]  V. G. Kohn,et al.  Phase-contrast microtomography with coherent high-energy synchrotron x rays , 1996 .

[8]  Timur E. Gureyev,et al.  Polychromatic cone-beam phase-contrast tomography , 2007 .

[9]  Atsushi Momose,et al.  Phase Tomography by X-ray Talbot Interferometry for Biological Imaging , 2006 .

[10]  Lajos Hajdu,et al.  An algorithm for discrete tomography , 2001 .

[11]  Peter Gritzmann,et al.  Discrete Tomography: Determination of Finite Sets by X-Rays , 1995, Universität Trier, Mathematik/Informatik, Forschungsbericht.

[12]  Kentaro Uesugi,et al.  High‐resolution visualization of airspace structures in intact mice via synchrotron phase‐contrast X‐ray imaging (PCXI) , 2008, Journal of anatomy.

[13]  T. Gureyev,et al.  Quasi-one-dimensional tomography , 1996 .

[14]  Robert J. Vanderbei,et al.  The Discrete Radon Transform and Its Approximate Inversion Via Linear Programming , 1997, Discret. Appl. Math..

[15]  Kees Joost Batenburg A Network Flow Algorithm for Binary Image Reconstruction from Few Projections , 2006, DGCI.

[16]  Paul J. Besl,et al.  Active, optical range imaging sensors , 1988, Machine Vision and Applications.

[17]  Christoph Schnörr,et al.  Discrete tomography by convex-concave regularization and D.C. programming , 2005, Discret. Appl. Math..

[18]  S. Wilkins,et al.  Simultaneous phase and amplitude extraction from a single defocused image of a homogeneous object , 2002, Journal of microscopy.

[19]  T. Gureyev,et al.  Quantitative quasi-local tomography using absorption and phase contrast , 2007 .

[20]  T. Gureyev,et al.  The binary dissector: phase contrast tomography of two- and three-material objects from few projections. , 2008, Optics express.

[21]  Sven de Vries,et al.  Approximating Binary Images from Discrete X-Rays , 2000, SIAM J. Optim..

[22]  Andrei V. Bronnikov,et al.  Reconstruction formulas in phase-contrast tomography , 1999 .

[23]  S. Wilkins,et al.  Phase-contrast imaging using polychromatic hard X-rays , 1996, Nature.

[24]  S. Wilkins,et al.  Linear algorithms for phase retrieval in the Fresnel region. 2. Partially coherent illumination , 2006 .

[25]  R M Leahy,et al.  Derivation and analysis of a filtered backprojection algorithm for cone beam projection data. , 1991, IEEE transactions on medical imaging.

[26]  L. Feldkamp,et al.  Practical cone-beam algorithm , 1984 .

[27]  F. Natterer The Mathematics of Computerized Tomography , 1986 .

[28]  Christoph Schnörr,et al.  Binary Tomography by Iterating Linear Programs from Noisy Projections , 2004, IWCIA.

[29]  Timur E. Gureyev,et al.  Tomography of objects with a priori known internal geometry , 1998 .

[30]  S. Wilkins,et al.  Generalized eikonal of partially coherent beams and its use in quantitative imaging. , 2004, Physical review letters.

[31]  S. Wilkins,et al.  X-ray phase-contrast microscopy and microtomography. , 2003, Optics express.

[32]  Ulrich Bonse,et al.  X-ray computed microtomography (μCT) using synchrotron radiation (SR) , 1996 .

[33]  I. Orion,et al.  Computed tomography of x-ray index of refraction using the diffraction enhanced imaging method. , 2000, Physics in medicine and biology.

[34]  Xizeng Wu,et al.  X-Ray cone-beam phase tomography formulas based on phase-attenuation duality. , 2005, Optics Express.