A Parametric Level-Set Approach to Simultaneous Object Identification and Background Reconstruction for Dual-Energy Computed Tomography

Dual-energy computerized tomography has gained great interest because of its ability to characterize the chemical composition of a material rather than simply providing relative attenuation images as in conventional tomography. The purpose of this paper is to introduce a novel polychromatic dual-energy processing algorithm, with an emphasis on detection and characterization of piecewise constant objects embedded in an unknown cluttered background. Physical properties of the objects, particularly the Compton scattering and photoelectric absorption coefficients, are assumed to be known with some level of uncertainty. Our approach is based on a level-set representation of the characteristic function of the object and encompasses a number of regularization techniques for addressing both the prior information we have concerning the physical properties of the object and the fundamental physics-based limitations associated with our ability to jointly recover the Compton scattering and photoelectric absorption properties of the scene. In the absence of an object with appropriate physical properties, our approach returns a null characteristic function and, thus, can be viewed as simultaneously solving the detection and characterization problems. Unlike the vast majority of methods that define the level-set function nonparametrically, i.e., as a dense set of pixel values, we define our level set parametrically via radial basis functions and employ a Gauss-Newton-type algorithm for cost minimization. Numerical results show that the algorithm successfully detects objects of interest, finds their shape and location, and gives an adequate reconstruction of the background.

[1]  E. Miller,et al.  Efficient determination of multiple regularization parameters in a generalized L-curve framework , 2002 .

[2]  S. Kay Fundamentals of statistical signal processing: estimation theory , 1993 .

[3]  Eric L. Miller,et al.  A projection-based level-set approach to enhance conductivity anomaly reconstruction in electrical resistance tomography , 2007 .

[4]  Simon R. Arridge,et al.  Image reconstruction in optical tomography using local basis functions , 2003, J. Electronic Imaging.

[5]  Ronald Fedkiw,et al.  Level set methods and dynamic implicit surfaces , 2002, Applied mathematical sciences.

[6]  V. A. Morozov,et al.  Methods for Solving Incorrectly Posed Problems , 1984 .

[7]  P. Joseph,et al.  A Method for Correcting Bone Induced Artifacts in Computed Tomography Scanners , 1978, Journal of computer assisted tomography.

[8]  Guowei Zhang,et al.  A practical reconstruction method for dual energy computed tomography , 2008 .

[9]  O. Dorn,et al.  Level set methods for inverse scattering , 2006 .

[10]  A. Macovski,et al.  Energy-selective reconstructions in X-ray computerised tomography , 1976, Physics in medicine and biology.

[11]  Ken D. Sauer,et al.  A unified approach to statistical tomography using coordinate descent optimization , 1996, IEEE Trans. Image Process..

[12]  Faming Liang,et al.  Statistical and Computational Inverse Problems , 2006, Technometrics.

[13]  J. Sethian,et al.  FRONTS PROPAGATING WITH CURVATURE DEPENDENT SPEED: ALGORITHMS BASED ON HAMILTON-JACOB1 FORMULATIONS , 2003 .

[14]  Willi A. Kalender,et al.  Dual-source cardiac computed tomography: image quality and dose considerations , 2008, European Radiology.

[15]  Max A. Meju,et al.  A quadratic programming approach for joint image reconstruction: mathematical and geophysical examples , 2005 .

[16]  Patrick Dupont,et al.  An iterative maximum-likelihood polychromatic algorithm for CT , 2001, IEEE Transactions on Medical Imaging.

[17]  Jeffrey A. Fessler,et al.  Edge-preserving tomographic reconstruction with nonlocal regularization , 1998, IEEE Transactions on Medical Imaging.

[18]  E Somersalo,et al.  Statistical inversion for medical x-ray tomography with few radiographs: II. Application to dental radiology. , 2003, Physics in medicine and biology.

[19]  Xiaoming Wang,et al.  A level set method for structural topology optimization , 2003 .

[20]  B.J. Heismann Atomic number measurement precision of spectral decomposition methods for CT , 2005, IEEE Nuclear Science Symposium Conference Record, 2005.

[21]  Isaac N. Bankman,et al.  Handbook of medical imaging , 2000 .

[22]  Fridrik Larusson,et al.  Hyperspectral imaging for diffuse optical tomography , 2009 .

[23]  Michael Unser,et al.  Variational B-Spline Level-Set: A Linear Filtering Approach for Fast Deformable Model Evolution , 2009, IEEE Transactions on Image Processing.

[24]  Guowei Zhang,et al.  Exact Reconstruction for Dual Energy Computed Tomography Using an H-L Curve Method , 2006, 2006 IEEE Nuclear Science Symposium Conference Record.

[25]  Yi Pan,et al.  Fast and Processor Efficient Parallel Matrix Multiplication Algorithms on a Linear Array With a Reconfigurable Pipelined Bus System , 1998, IEEE Trans. Parallel Distributed Syst..

[26]  R. Lewitt Alternatives to voxels for image representation in iterative reconstruction algorithms , 1992, Physics in medicine and biology.

[27]  D. Marquardt An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[28]  Arnold Neumaier,et al.  Introduction to Numerical Analysis , 2001 .

[29]  Anthony J. Yezzi,et al.  Curve evolution implementation of the Mumford-Shah functional for image segmentation, denoising, interpolation, and magnification , 2001, IEEE Trans. Image Process..

[30]  M. Meju,et al.  Characterization of heterogeneous near‐surface materials by joint 2D inversion of dc resistivity and seismic data , 2003 .

[31]  A. Macovski,et al.  Generalized image combinations in dual KVP digital radiography. , 1981, Medical physics.

[32]  Jeffrey A. Fessler,et al.  Maximum-likelihood dual-energy tomographic image reconstruction , 2002, SPIE Medical Imaging.

[33]  E. Alhava,et al.  Bone densitometry of the spine and femur in children by dual-energy x-ray absorptiometry. , 1992, Bone and mineral.

[34]  Josiane Zerubia,et al.  A Level Set Model for Image Classification , 1999, International Journal of Computer Vision.

[35]  Ville Kolehmainen,et al.  Limited Data X-Ray Tomography Using Nonlinear Evolution Equations , 2008, SIAM J. Sci. Comput..

[36]  C. Crawford,et al.  Dual energy computed tomography for explosive detection , 2006 .

[37]  Jorge Nocedal,et al.  Steering exact penalty methods for nonlinear programming , 2008, Optim. Methods Softw..

[38]  Jie Tang,et al.  Performance comparison between total variation (TV)-based compressed sensing and statistical iterative reconstruction algorithms , 2009, Physics in medicine and biology.

[39]  Joseph Fryar,et al.  Element specific imaging in computerised tomography using a tube source of X-rays and low energy-resolution detector system , 1993 .

[40]  W. D. Friedman,et al.  Review of dual-energy computed tomography techniques , 1990 .

[41]  P M Joseph,et al.  The effects of scatter in x-ray computed tomography. , 1982, Medical physics.

[42]  Patrick Dupont,et al.  Metal streak artifacts in X-ray computed tomography: a simulation study , 1998 .

[43]  G J Michael Tissue analysis using dual energy CT. , 1992, Australasian physical & engineering sciences in medicine.

[44]  M. Reiser,et al.  Material differentiation by dual energy CT: initial experience , 2007, European Radiology.

[45]  Jeffrey A. Fessler,et al.  Regularized Image Reconstruction Algorithms for Dual-Isotope Myocardial Perfusion SPECT (MPS) Imaging Using a Cross-Tracer Prior , 2011, IEEE Transactions on Medical Imaging.

[46]  Dimitri P. Bertsekas,et al.  Nonlinear Programming , 1997 .

[47]  Bruce R. Whiting,et al.  Noise simulation in x-ray CT , 2005, SPIE Medical Imaging.

[48]  J. O’Sullivan,et al.  Properties of preprocessed sinogram data in x-ray computed tomography. , 2006, Medical physics.

[49]  Maneesha Singh,et al.  Explosives detection systems (EDS) for aviation security , 2003, Signal Process..

[50]  G. Poludniowski Calculation of x-ray spectra emerging from an x-ray tube. Part II. X-ray production and filtration in x-ray targets. , 2007, Medical physics.

[51]  J. H. Hubbell,et al.  XCOM: Photon Cross Section Database (version 1.2) , 1999 .

[52]  B. Whiting,et al.  On two-parameter models of photon cross sections: application to dual-energy CT imaging. , 2006, Medical physics.

[53]  Hiroshi Yamashita,et al.  A Primal-Dual Exterior Point Method for Nonlinear Optimization , 2010, SIAM J. Optim..

[54]  W. Kalender,et al.  Evaluation of a prototype dual-energy computed tomographic apparatus. I. Phantom studies. , 1986, Medical physics.

[55]  K. Maute,et al.  A parametric level-set approach for topology optimization of flow domains , 2010 .

[56]  Eric L. Miller,et al.  Parametric Level Set Methods for Inverse Problems , 2010, SIAM J. Imaging Sci..

[57]  P M Evans,et al.  Extraction of primary signal from EPIDs using only forward convolution. , 1997, Medical physics.

[58]  Michael Gj Tissue analysis using dual energy CT. , 1992 .

[59]  Denis Friboulet,et al.  Compactly Supported Radial Basis Functions Based Collocation Method for Level-Set Evolution in Image Segmentation , 2007, IEEE Transactions on Image Processing.

[60]  Soung Chang Liew,et al.  Correction of energy-dependent systematic errors in dual-energy X-ray CT using a basis material coefficients transformation method , 1996 .

[61]  P. Suetens,et al.  Metal streak artifacts in X-ray computed tomography: a simulation study , 1998, 1998 IEEE Nuclear Science Symposium Conference Record. 1998 IEEE Nuclear Science Symposium and Medical Imaging Conference (Cat. No.98CH36255).

[62]  Eric L. Miller,et al.  Cortical constraint method for diffuse optical brain imaging , 2004, SPIE Optics + Photonics.

[63]  Eric L. Miller,et al.  Hyperspectral image reconstruction for diffuse optical tomography , 2011, Biomedical optics express.

[64]  Kaj Madsen,et al.  Methods for Non-Linear Least Squares Problems , 1999 .

[65]  F. Santosa A Level-set Approach Inverse Problems Involving Obstacles , 1995 .

[66]  W. Kalender,et al.  An algorithm for noise suppression in dual energy CT material density images. , 1988, IEEE transactions on medical imaging.

[67]  E. Miller,et al.  A shape reconstruction method for electromagnetic tomography using adjoint fields and level sets , 2000 .

[68]  Gene H. Golub,et al.  Generalized cross-validation as a method for choosing a good ridge parameter , 1979, Milestones in Matrix Computation.

[69]  E Somersalo,et al.  Statistical inversion for medical x-ray tomography with few radiographs: I. General theory. , 2003, Physics in medicine and biology.

[70]  F Verhaegen,et al.  SpekCalc: a program to calculate photon spectra from tungsten anode x-ray tubes , 2009, Physics in medicine and biology.

[71]  Charles R. Johnson,et al.  Topics in Matrix Analysis , 1991 .

[72]  Patrick J. La Rivière,et al.  Penalized-likelihood sinogram restoration for computed tomography , 2006, IEEE Transactions on Medical Imaging.

[73]  Jeffrey A. Fessler,et al.  Segmentation-free statistical image reconstruction for polyenergetic x-ray computed tomography with experimental validation , 2003 .

[74]  Philip M Evans,et al.  Calculation of x-ray spectra emerging from an x-ray tube. Part I. electron penetration characteristics in x-ray targets. , 2007, Medical physics.

[75]  C. Bouman,et al.  A recursive filter for noise reduction in statistical iterative tomographic imaging , 2006, Electronic Imaging.

[76]  Predrag Sukovic,et al.  Penalized weighted least-squares image reconstruction for dual energy X-ray transmission tomography , 2000, IEEE Transactions on Medical Imaging.

[77]  W. Clem Karl,et al.  A curve evolution approach to object-based tomographic reconstruction , 2003, IEEE Trans. Image Process..

[78]  T. M. Williams,et al.  Practical Methods of Optimization. Vol. 1: Unconstrained Optimization , 1980 .

[79]  Til Aach,et al.  Scattered radiation in cone beam computed tomography : analysis, quantification and compensation , 2007 .

[80]  Butler,et al.  RE-THINKING CHECKED-BAGGAGE SCREENING , 2002 .

[81]  Tony F. Chan,et al.  Active contours without edges , 2001, IEEE Trans. Image Process..

[82]  T. M. Williams Practical Methods of Optimization. Vol. 2 — Constrained Optimization , 1982 .

[83]  J. H. Hubbell,et al.  XCOM : Photon Cross Sections Database , 2005 .

[84]  Zhengrong Ying,et al.  Dual Energy Volumetric X-ray Tomographic Sensor for Luggage Screening , 2007, 2007 IEEE Sensors Applications Symposium.