Reflection Matrix Imaging for Wave Velocity Tomography

Besides controlling wave trajectory inside complex media, wave velocity constitutes a relevant bio-marker for medical imaging. In a transmission configuration, wave-front distortions can be unscrambled to provide a map of the wave velocity landscape $c(\mathbf{r})$. However, most in-vivo applications correspond to a reflection configuration for which only back-scattered echoes generated by short-scale fluctuations of $c(\mathbf{r})$ can be harvested. Under a single scattering assumption, this speckle wave-field cannot provide the long-scale variations of $c(\mathbf{r})$. In this paper, we go beyond the first Born approximation and show how a map of $c(\mathbf{r})$ can be retrieved in epi-detection. To that aim, a reflection matrix approach of wave imaging is adopted. While standard reflection imaging methods generally rely on confocal focusing operations, matrix imaging consists in decoupling the location of the incident and received focal spots. Following this principle, a self-portrait of the focusing process can be obtained around each point of the medium. The Gouy phase shift exhibited by each focal spot is leveraged to finely monitor the wave velocity distribution $c(\mathbf{r})$ inside the medium. Experiment in a tissue-mimicking phantom and numerical simulations are first presented to validate our method. Speed-of-sound tomography is then applied to ultrasound data collected on the liver of a difficult-to-image patient. Beyond paving the way towards quantitative ultrasound, our approach can also be extremely rewarding for standard imaging. Indeed, each echo can be assigned to the actual position of a scatterer. It allows an absolute measurement of distance, an observable often used for diagnosis but generally extremely sensitive to wave velocity fluctuations.

Three-Dimensional Ultrasound Matrix Imaging

Matrix imaging paves the way towards a next revolution in wave physics. Based on the response matrix recorded between a set of sensors, it enables an optimized compensation of aberration phenomena and multiple scattering events that usually drastically hinder the focusing process in heterogeneous media. Although it gave rise to spectacular results in optical microscopy or seismic imaging, the success of matrix imaging has been so far relatively limited with ultrasonic waves because wave control is generally only performed with a linear array of transducers. In this paper, we extend ultrasound matrix imaging to a 3D geometry. Switching from a 1D to a 2D probe enables a much sharper estimation of the transmission matrix that links each transducer and each medium voxel. Here, we first present an experimental proof of concept on a tissue-mimicking phantom through ex-vivo tissues and then, show the potential of 3D matrix imaging for transcranial applications.

Ultrasound Matrix Imaging. II. The distortion matrix for aberration correction over multiple isoplanatic patches

This is the second article in a series of two which report on a matrix approach for ultrasound imaging in heterogeneous media. This article describes the quantification and correction of aberration, i.e. the distortion of an image caused by spatial variations in the medium speed-of-sound. Adaptive focusing can compensate for aberration, but is only effective over a restricted area called the isoplanatic patch. Here, we use an experimentally-recorded matrix of reflected acoustic signals to synthesize a set of virtual transducers. We then examine wave propagation between these virtual transducers and an arbitrary correction plane. Such wave-fronts consist of two components: (i) An ideal geometric wave-front linked to diffraction and the input focusing point, and; (ii) Phase distortions induced by the speed-of-sound variations. These distortions are stored in a so-called distortion matrix, the singular value decomposition of which gives access to an optimized focusing law at any point. We show that, by decoupling the aberrations undergone by the outgoing and incoming waves and applying an iterative strategy, compensation for even high-order and spatially-distributed aberrations can be achieved. As a proof-of-concept, ultrasound matrix imaging (UMI) is applied to the in-vivo imaging of a human calf. A map of isoplanatic patches is retrieved and is shown to be strongly correlated with the arrangement of tissues constituting the medium. The corresponding focusing laws yield an ultrasound image with an optimal contrast and a transverse resolution close to the ideal value predicted by diffraction theory. UMI thus provides a flexible and powerful route towards computational ultrasound.

Ultrasound Matrix Imaging. I. The focused reflection matrix and the F-factor

This is the first article in a series of two dealing with a matrix approach \alex{for} aberration quantification and correction in ultrasound imaging. Advanced synthetic beamforming relies on a double focusing operation at transmission and reception on each point of the medium. Ultrasound matrix imaging (UMI) consists in decoupling the location of these transmitted and received focal spots. The response between those virtual transducers form the so-called focused reflection matrix that actually contains much more information than a raw ultrasound image. In this paper, a time-frequency analysis of this matrix is performed, which highlights the single and multiple scattering contributions as well as the impact of aberrations in the monochromatic and broadband regimes. Interestingly, this analysis enables the measurement of the incoherent input-output point spread function at any pixel of this image. A focusing criterion can then be built, and its evolution used to quantify the amount of aberration throughout the ultrasound image. In contrast to the standard coherence factor used in the literature, this new indicator is robust to multiple scattering and electronic noise, thereby providing a highly contrasted map of the focusing quality. As a proof-of-concept, UMI is applied here to the in-vivo study of a human calf, but it can be extended to any kind of ultrasound diagnosis or non-destructive evaluation.