Refraction-corrected ray-based inversion for three-dimensional ultrasound tomography of the breast

Ultrasound Tomography has seen a revival of interest in the past decade, especially for breast imaging, due to improvements in both ultrasound and computing hardware. In particular, three-dimensional ultrasound tomography, a fully tomographic method in which the medium to be imaged is surrounded by ultrasound transducers, has become feasible. In this paper, a comprehensive derivation and study of a robust framework for large-scale bent-ray ultrasound tomography in 3D for a hemispherical detector array is presented. Two ray-tracing approaches are derived and compared. More significantly, the problem of linking the rays between emitters and receivers, which is challenging in 3D due to the high number of degrees of freedom for the trajectory of rays, is analysed both as a minimisation and as a root-finding problem. The ray-linking problem is parameterised for a convex detection surface and three robust, accurate, and efficient ray-linking algorithms are formulated and demonstrated. To stabilise these methods, novel adaptive-smoothing approaches are proposed that control the conditioning of the update matrices to ensure accurate linking. The nonlinear UST problem of estimating the sound speed was recast as a series of linearised subproblems, each solved using the above algorithms and within a steepest descent scheme. The whole imaging algorithm was demonstrated to be robust and accurate on realistic data simulated using a full-wave acoustic model and an anatomical breast phantom, and incorporating the errors due to time-of-flight picking that would be present with measured data. This method can used to provide a low-artefact, quantitatively accurate, 3D sound speed maps. In addition to being useful in their own right, such 3D sound speed maps can be used to initialise full-wave inversion methods, or as an input to photoacoustic tomography reconstructions

Refraction-corrected ray-based inversion for three-dimensional ultrasound tomography of the breast

Ultrasound tomography (UST) has seen a revival of interest in the past decade, especially for breast imaging, due to improvements in both ultrasound and computing hardware. In particular, three-dimensional UST, a fully tomographic method in which the medium to be imaged is surrounded by ultrasound transducers, has become feasible. This has led to renewed attention on UST image reconstruction algorithms. In this paper, a comprehensive derivation and study of a robust framework for large-scale bent-ray UST in 3D for a hemispherical detector array is presented. Two ray

Limited Angle Ultrasound Tomography of the Compressed Breast

X-ray mammography is widely accepted as the clinical standard for breast cancer screening and diagnosis. However, reflection mode ultrasound has been known to outperform x-ray in screening performance in dense breasts. With newer modes of ultrasound, acoustic properties of breast tissue, such as the speed of sound and attenuation coefficient distributions, can be extracted from captured ultrasound signals and used to characterize breast tissue types and contribute to detection and diagnosis of malignancy. The same is possibly true for optical absorption via photoacoustic imaging. Recently, we have developed a dual-sided ultrasound scanner that can be integrated with existing x-ray mammographic systems and acquire images in the mammographic view and compression. Transmission imaging for speed of sound and attenuation coefficient in this geometry is termed limited angle tomography, as the beams at frequencies yielding high resolution cannot transit the long axis of the compressed breast. This approach, ideally, should facilitate the co-registration and comparisons between images from three modalities discussed here (x-ray, ultrasound and photoacoustic) and increase diagnostic detection confidence. However, potential limitations inherent in limited angle tomography have received minimal exploration up to this study, and existing imaging techniques developed for this approach are based on overly optimistic assumptions that hinder achievement of the desired image quality. This investigation of these problems should contribute valuable information to the validation and translation of the mammographically-configured, dual-sided ultrasound, or ultrasound and photoacoustic, scanner to the clinic. This dissertation first aims to extensively identify possible sources of error resulting from imaging in the limited angle tomography approach. Simulation findings mapping parametric conditions reveal that image artifacts arising in reflection mode (B-mode) can be modulated or mitigated by ultrasound gels with adequate acoustic properties. In addition, sound speed imaging was performed determining the level of significance for several key sources of error. Results suggest that imaging in transmission mode is the most sensitive to transducer misplacement in the signal propagation direction. This misplacement, however, could be minimized easily by routinely calibrating transducer positions. Next, this dissertation aims to advance speed of sound, attenuation, and photoacoustic image reconstruction algorithms for the limited angle tomography approach. This was done by utilizing both structural information of the imaged objects/tissues by means of the corresponding reflection mode images taken from the same imaging location, and a full acoustic modeling framework to account for complex acoustic interactions within the field of view. We have shown through simulations that both a priori information from reflection mode images and full acoustic modeling contribute to a noticeable improvement in the reconstructed images. Work done throughout the course of this dissertation should provide a foundation and insight necessary for improvements upon the existing dual-sided ultrasound scanner towards breast imaging in the clinic.PHDBiomedical EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttps://deepblue.lib.umich.edu/bitstream/2027.42/143944/1/rungroj_1.pd

Ray-based inversion accounting for scattering for biomedical ultrasound computed tomography

An efficient and accurate image reconstruction algorithm for ultrasound tomography (UST) is described and demonstrated, which can recover accurate sound speed distribution from acoustic time series measurements made in soft tissue. The approach is based on a second-order iterative minimisation of the difference between the measurements and a model based on a ray-approximation to the heterogeneous Green's function. It overcomes the computational burden of full-wave solvers while avoiding the drawbacks of time-of-flight methods. Through the use of a second-order iterative minimisation scheme, applied stepwise from low to high frequencies, the effects of scattering are incorporated into the inversion

Ray-based inversion accounting for scattering for biomedical ultrasound computed tomography

An efficient and accurate image reconstruction algorithm for ultrasound tomography in soft tissue is described and demonstrated, which can recover accurate sound speed distribution from acoustic time series measurements. The approach is based on a second-order iterative minimisation of the difference between the measurements and a model based on a ray-approximation to the heterogeneous Green's function. It overcomes the computational burden of full-wave solvers while avoiding the drawbacks of time-of-flight methods. Through the use of a second-order iterative minimisation scheme, applied stepwise from low to high frequencies, the effects of scattering are incorporated into the inversion

Métodos de reconstrucción en dominio temporal para tomografía por transmisión de ultrasonidos

Tesis inédita de la Universidad Complutense de Madrid, Facultad de Ciencias Físicas, Departamento de Física Atómica, Molecular y Nuclear, leída el 06-06-2017Breast cancer (BC) is the leading cause of cancer-related death for women in Europe, and the second one after lung cancer in the US [World Cancer Report, 2008]. Early detection is very important for the survival rate of BC, because the smaller the local extension of the neoplasia, the better the output of the surgical treatments employed. Besides, early detection increases the possibility of preserving the breast and decreases the probability of needing more invasive treatments [Secretaría de Salud, 2007, Alteri et al., 2011]. Mammography is currently the standard procedure employed for breast screening programs around the world. Nevertheless, its efficiency has been questioned lately because: (i) it generates many abnormal findings not related to cancer, (ii) it requires irradiating the patient and (iii) it has low specificity with dense breasts [Santen and Mansel, 2005]. Consequently, complementary techniques to mammography are being proposed to improve the detection and characterization of BC. Among these techniques, is the Ultrasound Computed Tomography (USCT), in reflection mode (which provides qualitative maps with the concentration of scatterers in the tissue), and transmission mode (which provides quantitative maps of the sound speed (SS) and the acoustic attenuation (AA) of the tissues). The images provided by the transmission modality have been proposed for BC detection as they can improve the detectability of malignancies in the breast [Mast, 2000, Duric et al., 2009]...El cáncer de mama (CM) es el cáncer más mortal entre las mujeres europeas, y el segundo más común en Estados Unidos [World Cancer Report, 2008]. La detección temprana es un factor que condiciona en gran medida la tasa de supervivencia a esta enfermedad, ya que a menor tamaño de la neoplasia detectada, mejores resultados pueden esperarse para los tratamientos quirúrgicos que se realicen. Además, la detección temprana aumenta la posibilidad de conservar la mama después de la cirugía y disminuye la necesidad de emplear otros tratamientos más invasivos[Secretaría de Salud, 2007, Alteri et al., 2011]. La mamografía es actualmente el procedimiento estándar que se emplea para el cribado del CM. Sin embargo, en los últimas años su eficiencia está siendo muy cuestionada por varios factores: (i) alta tasa de falsos positivos, (ii) requiere la irradiación del paciente y (iii) baja especificidad en mamas densas 2. Debido a lo anterior, para mejorar la detección y caracterización del CM se han propuesto varias técnicas complementarias. Entre ellas está la tomografía ultrasónica (TU), que es una técnica en desarrollo que presenta dos modalidades principales: la reflexión (proporciona mapas cualitativos de la concentración de dispersores en el tejido) y la transmisión (proporciona mapas cuantitativos de la velocidad y atenuación del sonido en el tejido). Los mapas del modo transmisión han sido propuestos como una eficiente alternativa, libre de radiación, para la detección del CM, ya que proporcionan alto contraste y especificidad [Mast, 2000, Duric et al., 2009]...Depto. de Estructura de la Materia, Física Térmica y ElectrónicaFac. de Ciencias FísicasTRUEunpu

High Resolution 3D Ultrasonic Breast Imaging by Time-Domain Full Waveform Inversion

Ultrasound tomography (UST) scanners allow quantitative images of the human breast's acoustic properties to be derived with potential applications in screening, diagnosis and therapy planning. Time domain full waveform inversion (TD-FWI) is a promising UST image formation technique that fits the parameter fields of a wave physics model by gradient-based optimization. For high resolution 3D UST, it holds three key challenges: Firstly, its central building block, the computation of the gradient for a single US measurement, has a restrictively large memory footprint. Secondly, this building block needs to be computed for each of the $10^3-10^4$ measurements, resulting in a massive parallel computation usually performed on large computational clusters for days. Lastly, the structure of the underlying optimization problem may result in slow progression of the solver and convergence to a local minimum. In this work, we design and evaluate a comprehensive computational strategy to overcome these challenges: Firstly, we introduce a novel gradient computation based on time reversal that dramatically reduces the memory footprint at the expense of one additional wave simulation per source. Secondly, we break the dependence on the number of measurements by using source encoding (SE) to compute stochastic gradient estimates. Also we describe a more accurate, TD-specific SE technique with a finer variance control and use a state-of-the-art stochastic LBFGS method. Lastly, we design an efficient TD multi-grid scheme together with preconditioning to speed up the convergence while avoiding local minima. All components are evaluated in extensive numerical proof-of-concept studies simulating a bowl-shaped 3D UST breast scanner prototype. Finally, we demonstrate that their combination allows us to obtain an accurate 442x442x222 voxel image with a resolution of 0.5mm using Matlab on a single GPU within 24h

Fully three-dimensional sound speed-corrected multi-wavelength photoacoustic breast tomography

Photoacoustic tomography is a contrast agent-free imaging technique capable of visualizing blood vessels and tumor-associated vascularization in breast tissue. While sophisticated breast imaging systems have been recently developed, there is yet much to be gained in imaging depth, image quality and tissue characterization capability before clinical translation is possible. In response, we have developed a hybrid photoacoustic and ultrasound-transmission tomographic system PAM3. The photoacoustic component has for the first time three-dimensional multi-wavelength imaging capability, and implements substantial technical advancements in critical hardware and software sub-systems. The ultrasound component enables for the first time, a three-dimensional sound speed map of the breast to be incorporated in photoacoustic reconstruction to correct for inhomogeneities, enabling accurate target recovery. The results demonstrate the deepest photoacoustic breast imaging to date namely 48 mm, with a more uniform field of view than hitherto, and an isotropic spatial resolution that rivals that of Magnetic Resonance Imaging. The in vivo performance achieved, and the diagnostic value of interrogating angiogenesis-driven optical contrast as well as tumor mass sound speed contrast, gives confidence in the system's clinical potential.Comment: 33 pages Main Body, 9 pages Supplementary Materia