A joint reconstruction and lambda tomography regularization technique for energy-resolved x-ray imaging

Here we present new joint reconstruction and regularization techniques inspired by ideas in microlocal analysis and lambda tomography, for the simultaneous reconstruction of the attenuation coefficient and electron density from X-ray transmission (i.e., X-ray CT) and backscattered data (assumed to be primarily Compton scattered). To demonstrate our theory and reconstruction methods, we consider the "parallel line segment" acquisition geometry of Webber and Miller ("Compton scattering tomography in translational geometries." Inverse Problems 36, no. 2 (2020): 025007), which is motivated by system architectures currently under development for airport security screening. We first present a novel microlocal analysis of the parallel line geometry which explains the nature of image artefacts when the attenuation coefficient and electron density are reconstructed separately. We next introduce a new joint reconstruction scheme for low effective $Z$ (atomic number) imaging ($Z<20$) characterized by a regularization strategy whose structure is derived from lambda tomography principles and motivated directly by the microlocal analytic results. Finally we show the effectiveness of our method in combating noise and image artefacts on simulated phantoms.

[1]  Shlomo Sternberg,et al.  Geometric Asymptotics, Revised edition , 1977 .

[2]  M. Meju,et al.  Joint two-dimensional DC resistivity and seismic travel time inversion with cross-gradients constraints , 2004 .

[3]  Per Christian Hansen,et al.  IR Tools: a MATLAB package of iterative regularization methods and large-scale test problems , 2017, Numerical Algorithms.

[4]  Ron Kimmel,et al.  A general framework for low level vision , 1998, IEEE Trans. Image Process..

[5]  Eric Todd Quinto,et al.  Microlocal Analysis in Tomography , 2015, Handbook of Mathematical Methods in Imaging.

[6]  Per Christian Hansen,et al.  Regularization Tools version 4.0 for Matlab 7.3 , 2007, Numerical Algorithms.

[7]  Eric Todd Quinto,et al.  Singularities of the X-ray transform and limited data tomography , 1993 .

[8]  Eric L. Miller,et al.  Stabilizing dual-energy x-ray computed tomography reconstructions using patch-based regularization , 2014, ArXiv.

[9]  J. H. Hubbell,et al.  Tables of X-Ray Mass Attenuation Coefficients and Mass Energy-Absorption Coefficients 1 keV to 20 MeV for Elements Z = 1 to 92 and 48 Additional Substances of Dosimetric Interest , 1995 .

[10]  Gaël Rigaud,et al.  Compton Scattering Tomography: Feature Reconstruction and Rotation-Free Modality , 2017, SIAM J. Imaging Sci..

[11]  L. Hörmander The analysis of linear partial differential operators , 1990 .

[12]  David Atkinson,et al.  Joint reconstruction of PET-MRI by exploiting structural similarity , 2014, Inverse Problems.

[13]  Robert E. Huffman,et al.  Photon Cross Sections , 1992 .

[14]  E. Miller,et al.  A geometric approach to joint inversion with applications to contaminant source zone characterization , 2013, 1303.5109.

[15]  Wojciech Samek,et al.  Learning the invisible: a hybrid deep learning-shearlet framework for limited angle computed tomography , 2018, Inverse Problems.

[16]  V. Palamodov A uniform reconstruction formula in integral geometry , 2011, 1111.6514.

[17]  Eldad Haber,et al.  Model Fusion and Joint Inversion , 2013, Surveys in Geophysics.

[18]  E. T. Quinto,et al.  Local Tomography in Electron Microscopy , 2008, SIAM J. Appl. Math..

[19]  W. Marsden I and J , 2012 .

[20]  Eric Todd Quinto,et al.  Artifacts in Incomplete Data Tomography with Applications to Photoacoustic Tomography and Sonar , 2014, SIAM J. Appl. Math..

[21]  Allan Hanbury,et al.  Metrics for evaluating 3D medical image segmentation: analysis, selection, and tool , 2015, BMC Medical Imaging.

[22]  L. Hörmander,et al.  The Analysis of Linear Partial Differential Operators IV , 1985 .

[23]  Gaël Rigaud,et al.  3D Compton scattering imaging: study of the spectrum and contour reconstruction , 2019, ArXiv.

[24]  V. Palamodov,et al.  An analytic reconstruction for the Compton scattering tomography in a plane , 2011 .

[25]  M. Nguyen,et al.  Compton scatter tomography in annular domains , 2019, Inverse Problems.

[26]  Guillermo Sapiro,et al.  Anisotropic diffusion of multivalued images with applications to color filtering , 1996, IEEE Trans. Image Process..

[27]  Eric L. Miller,et al.  On the Fusion of Compton Scatter and Attenuation Data for Limited-view X-ray Tomographic Applications , 2017, ArXiv.

[28]  Y. Egorov,et al.  Fourier Integral Operators , 1994 .

[29]  Hyunjoong Kim,et al.  Functional Analysis I , 2017 .

[30]  Marta M. Betcke,et al.  Multicontrast MRI Reconstruction with Structure-Guided Total Variation , 2015, SIAM J. Imaging Sci..

[31]  L. Hörmander The Analysis of Linear Partial Differential Operators III , 2007 .

[32]  Simon R. Arridge,et al.  Vector-Valued Image Processing by Parallel Level Sets , 2014, IEEE Transactions on Image Processing.

[33]  O. Öktem,et al.  Shearlets as feature extractor for semantic edge detection: the model-based and data-driven realm , 2019, Proceedings of the Royal Society A.

[34]  Nicola Wadeson Modelling and Correction of Scatter in a Switched Source Multi-Ring Detector X-ray CT machine , 2011 .

[35]  Stephen J. Norton,et al.  Compton scattering tomography , 1994 .

[36]  Bernadette N. Hahn,et al.  3D Compton scattering imaging and contour reconstruction for a class of Radon transforms , 2018 .

[37]  E. Miller,et al.  Multi-parameter acoustic imaging of uniform objects in inhomogeneous soft tissue , 2012, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[38]  Eric L. Miller,et al.  A Parametric Level-Set Approach to Simultaneous Object Identification and Background Reconstruction for Dual-Energy Computed Tomography , 2011, IEEE Transactions on Image Processing.

[39]  James W. Webber,et al.  Compton scattering tomography in translational geometries , 2019, Inverse Problems.

[40]  Erik L. Ritman,et al.  Local Tomography II , 1997, SIAM J. Appl. Math..

[41]  P. Hansen Rank-Deficient and Discrete Ill-Posed Problems: Numerical Aspects of Linear Inversion , 1987 .

[42]  Eric Todd Quinto,et al.  Analyzing Reconstruction Artifacts from Arbitrary Incomplete X-ray CT Data , 2017, SIAM J. Imaging Sci..

[43]  Eric Todd Quinto,et al.  The dependence of the generalized Radon transform on defining measures , 1980 .

[44]  W. Kalender X-ray computed tomography , 2006, Physics in medicine and biology.

[45]  Gregory Beylkin,et al.  Inversion Of The Generalized Radon Transform , 1983, Other Conferences.

[46]  Alfred K. Louis,et al.  Combining Image Reconstruction and Image Analysis with an Application to Two-Dimensional Tomography , 2008, SIAM J. Imaging Sci..

[47]  James Webber,et al.  Microlocal Analysis of a Compton Tomography Problem , 2019, SIAM J. Imaging Sci..

[48]  Per Christian Hansen,et al.  AIR Tools II: algebraic iterative reconstruction methods, improved implementation , 2017, Numerical Algorithms.

[49]  Eric L. Miller,et al.  Tensor-Based Formulation and Nuclear Norm Regularization for Multienergy Computed Tomography , 2013, IEEE Transactions on Image Processing.

[50]  Erik L. Ritman,et al.  Local tomography , 1992 .

[51]  J. H. Hubbell,et al.  Photon cross sections, attenuation coefficients, and energy absorption coefficients from 10 keV to 100 GeV , 1969 .

[52]  Max A. Meju,et al.  Joint two‐dimensional cross‐gradient imaging of magnetotelluric and seismic traveltime data for structural and lithological classification , 2007 .

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

[54]  Tony F. Chan,et al.  Color TV: total variation methods for restoration of vector-valued images , 1998, IEEE Trans. Image Process..

[55]  김도헌,et al.  Calcium-sulfate , 1998, Reactions Weekly.

[56]  Simon R. Arridge,et al.  Preface to special issue on joint reconstruction and multi-modality/multi-spectral imaging , 2020, Inverse Problems.

[57]  L. Hörmander Fourier integral operators. I , 1995 .

[58]  J. Winkler Titanium Dioxide , 2014 .