Open-source Gauss-Newton-based methods for refraction-corrected ultrasound computed tomography

This work presents refraction-corrected sound speed reconstruction techniques for transmission-based ultrasound computed tomography using a circular transducer array. Pulse travel times between element pairs can be calculated from slowness (the reciprocal of sound speed) using the eikonal equation. Slowness reconstruction is posed as a nonlinear least squares problem where the objective is to minimize the error between measured and forward-modeled pulse travel times. The Gauss-Newton method is used to convert this problem into a sequence of linear least-squares problems, each of which can be efficiently solved using conjugate gradients. However, the sparsity of ray-pixel intersection leads to ill-conditioned linear systems and hinders stable convergence of the reconstruction. This work considers three approaches for resolving the ill-conditioning in this sequence of linear inverse problems: 1) Laplacian regularization, 2) Bayesian formulation, and 3) resolution-filling gradients. The goal of this work is to provide an open-source example and implementation of the algorithms used to perform sound speed reconstruction, which is currently being maintained on Github: https://github.com/ rehmanali1994/refractionCorrectedUSCT.github.io

[1]  T. V. Oughton,et al.  Breast imaging in coronal planes with simultaneous pulse echo and transmission ultrasound. , 1981, Science.

[2]  S. Haykin.  Kalman Filtering and Neural Networks , 2001 .

[3]  Michael Zapf,et al.  Realization of an optimized 3D USCT , 2011, Medical Imaging.

[4]  Haijiang Zhang,et al.  An improved automatic time-of-flight picker for medical ultrasound tomography. , 2009, Ultrasonics.

[5]  Lihong V. Wang,et al.  Enhancement of photoacoustic tomography by ultrasonic computed tomography based on optical excitation of elements of a full-ring transducer array. , 2013, Optics letters.

[6]  Martin Vetterli,et al.  Robust ultrasound travel-time tomography using the bent ray model , 2010, Medical Imaging.

[7]  A. Kak,et al.  Simultaneous Algebraic Reconstruction Technique (SART): A Superior Implementation of the Art Algorithm , 1984, Ultrasonic imaging.

[8]  Alex M. Andrew,et al.  Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science (2nd edition) , 2000 .

[9]  Ichiro Sakuma,et al.  Novel automatic first-arrival picking method for ultrasound sound-speed tomography , 2015 .

[10]  T. M. Kolb,et al.  Occult cancer in women with dense breasts: detection with screening US--diagnostic yield and tumor characteristics. , 1998, Radiology.

[11]  Cuiping Li,et al.  Travel time denoising in ultrasound tomography , 2012, Medical Imaging.

[12]  Clifford H. Thurber,et al.  Parameter estimation and inverse problems , 2005 .

[13]  N. Duric,et al.  Combining time of flight and diffraction tomography for high resolution breast imaging: initial in vivo results (L). , 2012, The Journal of the Acoustical Society of America.

[14]  Uwe Fischer,et al.  Contrast Enhancement on Cone-Beam Breast-CT for Discrimination of Breast Cancer Immunohistochemical Subtypes123 , 2017, Translational oncology.

[15]  S J Norton,et al.  Correcting for ray refraction in velocity and attenuation tomography: a perturbation approach. , 1982, Ultrasonic imaging.

[16]  Ming Jiang,et al.  Convergence of the simultaneous algebraic reconstruction technique (SART) , 2003, Conference Record of Thirty-Fifth Asilomar Conference on Signals, Systems and Computers (Cat.No.01CH37256).

[17]  Nicole V. Ruiter,et al.  USCT data challenge , 2017, Medical Imaging.

[18]  Nicole V. Ruiter,et al.  USCT reference data base: conclusions from the first SPIE USCT data challenge and future directions , 2018, Medical Imaging.

[19]  T. M. Kolb,et al.  Comparison of the performance of screening mammography, physical examination, and breast US and evaluation of factors that influence them: an analysis of 27,825 patient evaluations. , 2002, Radiology.

[20]  Cuiping Li,et al.  Waveform inversion with source encoding for breast sound speed reconstruction in ultrasound computed tomography , 2014, IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control.

[21]  J. Sethian,et al.  Fast-phase space computation of multiple arrivals , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[22]  B T Cox,et al.  k-Wave: MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields. , 2010, Journal of biomedical optics.

[23]  Youli Quan,et al.  Sound-speed tomography using first-arrival transmission ultrasound for a ring array , 2007, SPIE Medical Imaging.

[24]  J. J. Gisvold,et al.  Ultrasound transmission computed tomography of the breast. , 1984, Radiology.

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

[26]  Stephen J. Wright,et al.  Conjugate Gradient Methods , 1999 .

[27]  N. Duric,et al.  In vivo breast sound-speed imaging with ultrasound tomography. , 2009, Ultrasound in medicine & biology.

[28]  Richard Su,et al.  Clinical feasibility study of combined optoacoustic and ultrasonic imaging modality providing coregistered functional and anatomical maps of breast tumors , 2012, Photonics West - Biomedical Optics.

[29]  Aly A. Farag,et al.  MultiStencils Fast Marching Methods: A Highly Accurate Solution to the Eikonal Equation on Cartesian Domains , 2007, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[30]  N. Duric,et al.  Detection of breast cancer with ultrasound tomography: first results with the Computed Ultrasound Risk Evaluation (CURE) prototype. , 2007, Medical physics.

[31]  Neb Duric,et al.  Resolution limitation of travel time tomography: beyond the first Fresnel zone , 2013, Medical Imaging.

[32]  J. Baerentzen,et al.  On the Implementation of Fast Marching Methods for 3d Lattices , .

[33]  R. E. Kalman,et al.  A New Approach to Linear Filtering and Prediction Problems , 2002 .

[34]  D. Borup,et al.  Non-linear inverse scattering: high resolution quantitative breast tissue tomography. , 2012, The Journal of the Acoustical Society of America.

[35]  J. Krebs,et al.  Fast full-wavefield seismic inversion using encoded sources , 2009 .

[36]  J. Sethian,et al.  3-D traveltime computation using the fast marching method , 1999 .

[37]  T. Başar,et al.  A New Approach to Linear Filtering and Prediction Problems , 2001 .