Specular Beamforming

Quantifying the effect of fiber pennation angle on shear wave viscoelastography estimates: In silico and phantom studies

Ultrasound shear wave elastography can be useful for assessing muscle pathology. The effect of anisotropy on shear wave elasticity estimates of skeletal muscle has been reported previously. However, muscle is inherently viscoelastic, and hence, tissue viscosity is also an important material parameter to assess. The goal of this study was to systematically quantify the effect of fiber pennation angle on shear wave viscoelasticity imaging estimates. Numerical phantom simulations of skeletal muscle-mimicking phantoms were analyzed. Anisotropic polyvinyl alcohol phantoms embedded with polysulfone fibers were developed to mimic the viscoelasticity and appearance of muscle in B-mode images. Shear wave dispersion analysis, assuming a Kelvin-Voigt model, was performed to estimate the shear modulus and viscosity of the phantoms along the fibers (in-plane) and across the fibers (cross-plane) with varying pennation angles (0°-30°). A decreasing trend was observed in shear modulus estimates with increasing fiber pennation angle in the in-plane orientation for all phantoms. Notably, simulations showed that viscosity estimates decreased with increasing angle. These results provide a systematic quantification of the effect of fiber pennation angle on viscoelastic estimates under controlled conditions, which will be useful for assessing the performance of shear wave viscoelasticity imaging approaches for muscle assessment.

Enhanced Needle Visualization With Reflection Tuned Apodization Based on the Radon Transform for Ultrasound Imaging

In ultrasound (US)-guided interventions, accurately tracking and visualizing needles during in-plane insertions are significant challenges due to strong directional specular reflections. These reflections violate the geometrical delay and apodization estimations in the conventional delay and sum beamforming (DASB) degrading the visualization of needles. This study proposes a novel reflection tuned apodization (RTA) to address this issue and facilitate needle enhancement through DASB. The method leverages both temporal and angular information derived from the Radon transforms of the radio frequency (RF) data from plane-wave imaging to filter the specular reflections from the needle and their directivity. The directivity information is translated into apodization center maps through time-to-space mapping in the Radon domain, which is subsequently integrated into DASB. We assess the influence of needle angulations, projection angles in the Radon transform, needle gauge sizes, and the presence of multiple specular interfaces on the approach. The analysis shows that the method surpasses conventional DASB in enhancing the image quality of needle interfaces while preserving the diffuse scattering from the surrounding tissues without significant computational overhead. The work offers promising prospects for improved outcomes in US-guided interventions and better insights into characterizing US reflections with Radon transforms.

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.

On the physics of ultrasound transmission for in-plane needle tracking in guided interventions

Objective.In ultrasound (US) guided interventions, the accurate visualization and tracking of needles is a critical challenge, particularly during in-plane insertions. An inaccurate identification and localization of needles lead to severe inadvertent complications and increased procedure times. This is due to the inherent specular reflections from the needle with directivity depending on the angle of incidence of the US beam, and the needle inclination.Approach.Though several methods have been proposed for improved needle visualization, a detailed study emphasizing the physics of specular reflections resulting from the interaction of transmitted US beam with the needle remains to be explored. In this work, we discuss the properties of specular reflections from planar and spherical wave US transmissions respectively through multi-angle plane wave (PW) and synthetic transmit aperture (STA) techniques for in-plane needle insertion angles between 15°-50°.Main Results.The qualitative and quantitative results from simulations and experiments reveal that the spherical waves enable better visualization and characterization of needles than planar wavefronts. The needle visibility in PW transmissions is severely degraded by the receive aperture weighting during image reconstruction than STA due to greater deviation in reflection directivity. It is also observed that the spherical wave characteristics starts to alter to planar characteristics due to wave divergence at large needle insertion depths.Significance.The study highlights that synergistic transmit-receive imaging schemes addressing the physical properties of reflections from the transmit wavefronts are imperative for the precise imaging of needle interfaces and hence have strong potential in elevating the quality of outcomes from US guided interventional practices.

The influence of intra-cortical microstructure on the contrast in ultrasound images of the cortex of long bones: A 2D simulation study

Decreased thickness of the bone cortex due to bone loss in the course of ageing and osteoporosis is associated with reduced bone strength. Cortical thickness measurement from ultrasound images was recently demonstrated in young adults. This requires the identification of both the outer (periosteum) and inner (endosteum) surfaces of the bone cortex. However, with bone loss, the cortical porosity and the size of the vascular pores increase resulting in enhanced ultrasound scattering which may prevent the detection of the endosteum. The aim of this work was to study the influence of cortical bone microstructure variables, such as porosity and pore size, on the contrast of the endosteum in ultrasound images. We wanted to estimate the range of these variables for which ultrasound imaging of the endosteum is feasible. We generated synthetic data using a two-dimensional time-domain code to simulate the propagation of elastodynamic waves. A synthetic aperture imaging sequence with an array transducer operating at a center frequency of 2.5 MHz was used. The numerical simulations were conducted for 105 cortical microstructures obtained from high resolution X-ray computed tomography images of ex vivo bone samples with a porosity ranging from 2% to 24 %. Images were reconstructed using a delay-and-sum (DAS) algorithm with optimized f-number, correction of refraction at the periosteum, and sample-specific wave-speed. We observed a range variation of 18 dB of endosteum contrast in our data set depending on the bone microstructure. We found that as porosity increases, speckle intensity inside the bone cortex increases whereas the intensity of the signal from the endosteum decreases. Also, a microstructure with large pores (diameter >250 μm) was associated with poor endosteum visibility, compared with a microstructure with equal porosity but a more narrow distribution of pore sizes. These findings suggest that ultrasound imaging of the bone cortex with a probe operating at a central frequency of 2.5 MHz using refraction-corrected DAS is capable of detecting the endosteum of a cortex with moderate porosity (less than about 10%) if the largest pores remain smaller than about 200 μm.

Ultrasound Matrix Imaging-Part 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. After a numerical validation of the process, ultrasound matrix imaging (UMI) is applied to the in-vivo imaging of a gallbladder. A map of isoplanatic modes 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 drastically improved contrast and transverse resolution. UMI thus provides a flexible and powerful route towards computational ultrasound.

Correlation-based ultrasound imaging of strong reflectors with phase coherence filtering

Two main metrics are usually employed to assess the quality of medical ultrasound (US) images, namely the contrast and the spatial resolution. A number of imaging algorithms have been proposed to improve one of those metrics, often at the expense of the other one. This paper presents the application of a correlation-based ultrasound imaging method, called Excitelet, to medical US imaging applications and the inclusion of a new Phase Coherence (PC) metric within its formalism. The main idea behind this algorithm, originally developed and validated for Non-Destructive Testing (NDT) applications, is to correlate a reference signal database with the measured signals acquired from a transducer array. In this paper, it is shown that improved lateral resolutions and a reduction of imaging artifacts are obtained over the Synthetic Aperture Focusing Technique (SAFT) when using Excitelet in conjunction with a PC filter. This novel method shows potential for the imaging of specular reflectors, such as invasive surgical tools. Numerical and experimental results presented in this paper demonstrate the benefit, in terms of contrast and resolution, of using the Excitelet method combined with PC for the imaging of strong reflectors.

Pixel Intensity Vector Field: An Inside Out Approach of Looking at Ultrasound Reflections from the Lung at High Frame Rates

Ultrasound (US) imaging is becoming the routine modality for the diagnosis and prognosis of lung pathologies. Lung US imaging relies on artifacts from acoustic impedance (Z) mismatches to distinguish and interpret the normal and pathological lung conditions. The air-pleura interface of the normal lung displays specularity due to the huge Z mismatches. However, in the presence of pathologies, the interface alters exhibiting a diffuse behavior due to increased density and reduced spatial distribution of air in the sub-pleural space. The specular or the diffuse behavior influences the reflected acoustic intensity distribution. This study aims to understand the reflection pattern in a normal and pathological lung through a novel approach of determining pixel-level acoustic intensity vector field (IVF) at high frame rates. Detailed lung modeling procedures using k-Wave US toolbox under normal, edematous, and consolidated conditions are illustrated. The analysis of the IVF maps on the three lung models clearly shows the drifting of the air-pleura interface from specular to diffuse with the severity of the pathology.Clinical Relevance- The presented acoustic simulation lung models are an aid to teaching and research by providing a quick visual and intuitive understanding of lung ultrasound physics. The proposed intensity vector field maps are supplementary information to aid diagnostics and characterization of any tissue composed of specular and diffuse components.

Sub-wavelength lateral detection of tissue-approximating masses using an ultrasonic metamaterial lens

Practically applied techniques for ultrasonic biomedical imaging employ delay-and-sum (DAS) beamforming which can resolve two objects down to 2.1λ within the acoustic Fresnel zone. Here, we demonstrate a phononic metamaterial lens (ML) for detection of laterally subwavelength object features in tissue-like phantoms beyond the phononic crystal evanescent zone and Fresnel zone of the emitter. The ML produces metamaterial collimation that spreads 8x less than the emitting transducer. Utilizing collimation, 3.6x greater lateral resolution beyond the Fresnel zone limit was achieved. Both hard objects and tissue approximating masses were examined in gelatin tissue phantoms near the Fresnel zone limit. Lateral dimensions and separation were resolved down to 0.50λ for hard objects, with tissue approximating masses slightly higher at 0.73λ. The work represents the application of a metamaterial for spatial characterization, and subwavelength resolution in a biosystem beyond the Fresnel zone limit.

Traditional methods for ultrasound detection in biomedical application suffer from limited lateral resolution. Here, the authors show that a phononic metamaterial lens can be used for spatial characterisation of subwavelength objects, even beyond the Fresnel zone of the emitting transducer.

Distortion matrix approach for ultrasound imaging of random scattering media

Focusing waves inside inhomogeneous media is a fundamental problem for imaging. Spatial variations of wave velocity can strongly distort propagating wave fronts and degrade image quality. Adaptive focusing can compensate for such aberration but is only effective over a restricted field of view. Here, we introduce a full-field approach to wave imaging based on the concept of the distortion matrix. This operator essentially connects any focal point inside the medium with the distortion that a wave front, emitted from that point, experiences due to heterogeneities. A time-reversal analysis of the distortion matrix enables the estimation of the transmission matrix that links each sensor and image voxel. Phase aberrations can then be unscrambled for any point, providing a full-field image of the medium with diffraction-limited resolution. Importantly, this process is particularly efficient in random scattering media, where traditional approaches such as adaptive focusing fail. Here, we first present an experimental proof of concept on a tissue-mimicking phantom and then, apply the method to in vivo imaging of human soft tissues. While introduced here in the context of acoustics, this approach can also be extended to optical microscopy, radar, or seismic imaging.

Towards A Pixel-Level Reconfigurable Digital Beamforming Core for Ultrasound Imaging

Ultrasound (US) imaging systems typically employ a single beamforming scheme which is the delay and sum (DAS) beamforming due to its reduced complexity. However, DAS results in images with limited resolution and contrast. The limitations of DAS have been overcome by, delay multiply and sum (DMAS) beamforming, making it especially preferable in cases where finer image details are required in larger depth of scans for an accurate diagnosis. But, DMAS is confined to transducer frequencies where the generated harmonics also fall in the processable frequency range of the US system. However, if US systems could provide the flexibility to reconfigure beamforming considering the restrictions of each beamforming scheme, it is possible to select the best beamforming according to the clinical requirement and system constraints. This work is a fundamental step towards enabling reconfigurable beamforming for on-the-fly selection among the DAS and DMAS beamforming schemes, with low reconfiguration overhead, specifically for each imaging scenario to aid better diagnosis. Two novel architectures are proposed, that reconfigures between DAS and DMAS beamforming as a function of transducer's center frequency with minimum additional computational overhead. The implementation results of the proposed architectures on xc7z010clg400-1 FPGA are reported. The possibilities of pixel-level beamforming reconfigurability, where the different tissue regions are beamformed with either DAS or DMAS are also shown through simulations.

The Generalized Contrast-to-Noise Ratio: A Formal Definition for Lesion Detectability

In the last 30 years, the contrast-to-noise ratio (CNR) has been used to estimate the contrast and lesion detectability in ultrasound images. Recent studies have shown that the CNR cannot be used with modern beamformers, as dynamic range alterations can produce arbitrarily high CNR values with no real effect on the probability of lesion detection. We generalize the definition of CNR based on the overlap area between two probability density functions. This generalized CNR (gCNR) is robust against dynamic range alterations; it can be applied to all kind of images, units, or scales; it provides a quantitative measure for contrast; and it has a simple statistical interpretation, i.e., the success rate that can be expected from an ideal observer at the task of separating pixels. We test gCNR on several state-of-the-art imaging algorithms and, in addition, on a trivial compression of the dynamic range. We observe that CNR varies greatly between the state-of-the-art methods, with improvements larger than 100%. We observe that trivial compression leads to a CNR improvement of over 200%. The proposed index, however, yields the same value for compressed and uncompressed images. The tested methods showed mismatched performance in terms of lesion detectability, with variations in gCNR ranging from -0.08 to +0.29. This new metric fixes a methodological flaw in the way we study contrast and allows us to assess the relevance of new imaging algorithms.

Ultrasonic Synthetic-Aperture Interface Imaging

Synthetic-aperture (SA) imaging is a popular method to visualize the reflectivity of an object from ultrasonic reflections. The method yields an image of the (volume) contrast in acoustic impedance with respect to the embedding. Typically, constant mass density is assumed in the underlying derivation. Due to the band-limited nature of the recorded data, the image is blurred in space, which is quantified by the associated point spread function. SA volume imaging is valid under the Born approximation, where it is assumed that the contrast is weak. When objects are large with respect to the wavelength, it is questionable whether SA volume imaging should be the method-of-choice. Herein, we propose an alternative solution that we refer to as SA interface imaging. This approach yields a vector image of the discontinuities of acoustic impedance at the tissue interfaces. Constant wave speed is assumed in the underlying derivation. The image is blurred in space by a tensor, which we refer to as the interface spread function. SA interface imaging is valid under the Kirchhoff approximation, where it is assumed that the wavelength is small compared to the spatial dimensions of the interfaces. We compare the performance of volume and interface imaging on synthetic data and on experimental data of a gelatin cylinder with a radius of 75 wavelengths, submerged in water. As expected, the interface image peaks at the gelatin-water interface, while the volume image exposes a peak and trough on opposing sides of the interface.