Emmanuel Montagnon

Stolt's f-k migration for plane wave ultrasound imaging

Ultrafast ultrasound is an emerging modality that offers new perspectives and opportunities in medical imaging. Plane wave imaging (PWI) allows one to attain very high frame rates by transmission of planar ultrasound wavefronts. As a plane wave reaches a given scatterer, the latter becomes a secondary source emitting upward spherical waves and creating a diffraction hyperbola in the received RF signals. To produce an image of the scatterers, all the hyperbolas must be migrated back to their apexes. To perform beamforming of plane wave echo RFs and return high-quality images at high frame rates, we propose a new migration method carried out in the frequency-wavenumber (f-k) domain. The f-k migration for PWI has been adapted from the Stolt migration for seismic imaging. This migration technique is based on the exploding reflector model (ERM), which consists in assuming that all the scatterers explode in concert and become acoustic sources. The classical ERM model, however, is not appropriate for PWI. We showed that the ERM can be made suitable for PWI by a spatial transformation of the hyperbolic traces present in the RF data. In vitro experiments were performed to outline the advantages of PWI with Stolts f-k migration over the conventional delay-and-sum (DAS) approach. The Stolts f-k migration was also compared with the Fourier-based method developed by J.-Y. Lu. Our findings show that multi-angle compounded f-k migrated images are of quality similar to those obtained with a stateof- the-art dynamic focusing mode. This remained true even with a very small number of steering angles, thus ensuring a highly competitive frame rate. In addition, the new FFT-based f-k migration provides comparable or better contrast-to-noise ratio and lateral resolution than the Lus and DAS migration schemes. Matlab codes for the Stolts f-k migration for PWI are provided.

Shear Wave Induced Resonance Elastography of Venous Thrombi: A Proof-of-Concept

Shear wave induced resonance elastography (SWIRE) is proposed for deep venous thrombosis (DVT) elasticity assessment. This new imaging technique takes advantage of properly polarized shear waves to induce resonance of a confined mechanical heterogeneity. Realistic phantoms (n = 9) of DVT total and partial clot occlusions with elasticities from 406 to 3561 Pa were built for in vitro experiments. An ex vivo study was also performed to evaluate the elasticity of two fresh porcine venous thrombi in a pig model. Transient shear waves at 45-205 Hz were generated by the vibration of a rigid plate (plane wavefront) or by a needle to simulate a radiation pressure on a line segment (cylindrical wavefront). Induced propagation of shear waves was imaged with an ultrafast ultrasound scanner and a finite element method was developed to simulate tested experimental conditions. An inverse problem was then formulated considering the first resonance frequency of the DVT inclusion. Elasticity agreements between SWIRE and a reference spectroscopy instrument (RheoSpectris) were found in vitro for total clots either in plane (r2 = 0.989) or cylindrical (r2 = 0.986) wavefront configurations. For total and partial clots, elasticity estimation errors were 9.0 ±4.6% and 9.3 ±11.3%, respectively. Ex vivo, the blood clot elasticity was 498 ±58 Pa within the inferior vena cava and 436 ±45 Pa in the right common iliac vein (p = 0.22). To conclude, the SWIRE technique seems feasible to quantitatively assess blood clot elasticity in the context of DVT ultrasound imaging.

Real-time processing in dynamic ultrasound elastography: A GPU-based implementation using CUDA

This paper addresses the computational cost of the normalized cross-correlation (NCC) algorithm in ultrasound elastography. Parallel implementations of the NCC algorithm based on multicore architectures and a graphical processor unit (GPU) are formulated and applied to radio-frequency (RF) data from dynamic elastography experiments. Compared to single computer processor unit (CPU) performances, results show that parallel implementation of the NCC algorithm allows speedups of less than 5 for multi-threaded execution on CPU and up to 85 using a GPU. Processing frame rates from 80 to 173 sec-1 have been achieved for large fields of view with good spatial resolution. The trade-off between accuracy, spatial resolution and computational cost in displacement estimation using the NCC algorithm therefore appears obsolete. Open source codes for implementing the NCC algorithm on GPU are made available at www.lbum-crchum.com.

Rheological assessment of a polymeric spherical structure using a three-dimensional shear wave scattering model in dynamic spectroscopy elastography

With the purpose of assessing localized rheological behavior of pathological tissues using ultrasound dynamic elastography, an analytical shear wave scattering model was used in an inverse problem framework. The proposed method was adopted to estimate the complex shear modulus of viscoelastic spheres from 200 to 450 Hz. The inverse problem was formulated and solved in the frequency domain, allowing assessment of the complex viscoelastic shear modulus at discrete frequencies. A representative rheological model of the spherical obstacle was determined by comparing storage and loss modulus behaviors with Kelvin¿Voigt, Maxwell, Zener, and Jeffrey models. The proposed inversion method was validated by using an external vibrating source and acoustic radiation force. The estimation of viscoelastic properties of three-dimensional spheres made softer or harder than surrounding tissues did not require a priori rheological assumptions. The proposed method is intended to be applied in the context of breast cancer imaging.

Frequency adaptation for enhanced radiation force amplitude in dynamic elastography

In remote dynamic elastography, the amplitude of the generated displacement field is directly related to the amplitude of the radiation force. Therefore, displacement improvement for better tissue characterization requires the optimization of the radiation force amplitude by increasing the push duration and/or the excitation amplitude applied on the transducer. The main problem of these approaches is that the Food and Drug Administration (FDA) thresholds for medical applications and transducer limitations may be easily exceeded. In the present study, the effect of the frequency used for the generation of the radiation force on the amplitude of the displacement field was investigated. We found that amplitudes of displacements generated by adapted radiation force sequences were greater than those generated by standard nonadapted ones (i.e., single push acoustic radiation force impulse and supersonic shear imaging). Gains in magnitude were between 20 to 158% for in vitro measurements on agar-gelatin phantoms, and 170 to 336% for ex vivo measurements on a human breast sample, depending on focus depths and attenuations of tested samples. The signal-to-noise ratio was also improved more than 4-fold with adapted sequences. We conclude that frequency adaptation is a complementary technique that is efficient for the optimization of displacement amplitudes. This technique can be used safely to optimize the deposited local acoustic energy without increasing the risk of damaging tissues and transducer elements.

Viscoelastic characterization of elliptical mechanical heterogeneities using a semi-analytical shear-wave scattering model for elastometry measures.

This paper presents a semi-analytical model of shear wave scattering by a viscoelastic elliptical structure embedded in a viscoelastic medium, and its application in the context of dynamic elastography imaging. The commonly used assumption of mechanical homogeneity in the inversion process is removed introducing a priori geometrical information to model physical interactions of plane shear waves with the confined mechanical heterogeneity. Theoretical results are first validated using the finite element method for various mechanical configurations and incidence angles. Secondly, an inverse problem is formulated to assess viscoelastic parameters of both the elliptic inclusion and its surrounding medium, and applied in vitro to characterize mechanical properties of agar-gelatin phantoms. The robustness of the proposed inversion method is then assessed under various noise conditions, biased geometrical parameters and compared to direct inversion, phase gradient and time-of-flight methods. The proposed elastometry method appears reliable in the context of estimating confined lesion viscoelastic parameters.

Generation of remote adaptive torsional shear waves with an octagonal phased array to enhance displacements and reduce variability of shear wave speeds: comparison with quasi-plane shear wavefronts.

A method based on adaptive torsional shear waves (ATSW) is proposed to overcome the strong attenuation of shear waves generated by a radiation force in dynamic elastography. During the inward propagation of ATSW, the magnitude of displacements is enhanced due to the convergence of shear waves and constructive interferences. The proposed method consists in generating ATSW fields from the combination of quasi-plane shear wavefronts by considering a linear superposition of displacement maps. Adaptive torsional shear waves were experimentally generated in homogeneous and heterogeneous tissue mimicking phantoms, and compared to quasi-plane shear wave propagations. Results demonstrated that displacement magnitudes by ATSW could be up to 3 times higher than those obtained with quasi-plane shear waves, that the variability of shear wave speeds was reduced, and that the signal-to-noise ratio of displacements was improved. It was also observed that ATSW could cause mechanical inclusions to resonate in heterogeneous phantoms, which further increased the displacement contrast between the inclusion and the surrounding medium. This method opens a way for the development of new noninvasive tissue characterization strategies based on ATSW in the framework of our previously reported shear wave induced resonance elastography (SWIRE) method proposed for breast cancer diagnosis.

Acoustic radiation force induced resonance elastography of coagulating blood: theoretical viscoelasticity modeling and ex vivo experimentation

Deep vein thrombosis is a common vascular disease that can lead to pulmonary embolism and death. The early diagnosis and clot age staging are important parameters for reliable therapy planning. This article presents an acoustic radiation force induced resonance elastography method for the viscoelastic characterization of clotting blood. The physical concept of this method relies on the mechanical resonance of the blood clot occurring at specific frequencies. Resonances are induced by focusing ultrasound beams inside the sample under investigation. Coupled to an analytical model of wave scattering, the ability of the proposed method to characterize the viscoelasticity of a mimicked venous thrombosis in the acute phase is demonstrated. Experiments with a gelatin-agar inclusion sample of known viscoelasticity are performed for validation and establishment of the proof of concept. In addition, an inversion method is applied in-vitro for the kinetic monitoring of the blood coagulation process of six human blood samples obtained from two volunteers. The computed elasticity and viscosity values of blood samples at the end of the 90 min kinetics were estimated at 411 ± 71 Pa and 0.25 ± 0.03 Pa.s for volunteer #1, and 387 ± 35 Pa and 0.23 ± 0.02 Pa.s for volunteer #2, respectively. The proposed method allowed reproducible time-varying thrombus viscoelastic measurements from samples having physiological dimensions.

Validation and application of a nondestructive and contactless method for rheological evaluation of biomaterials

Hydrogels are extensively used for tissue engineering, cell therapy or controlled release of bioactive factors. Nondestructive techniques that can follow their viscoelastic properties during polymerization, remodeling, and degradation are needed, since these properties are determinant for their in vivo efficiency. In this work, we proposed the viscoelastic testing of bilayered materials (VeTBiM) as a new method for nondestructive and contact-less mechanical characterization of soft materials. The VeTBiM method measures the dynamic displacement response of a material, to a low amplitude vibration in order to characterize its viscoelastic properties. We validated VeTBiM by comparing data obtained on various agar and chitosan hydrogels with data from rotational rheometry, and compression tests. We then investigated its potential to follow the mechanical properties of chitosan hydrogels during gelation and in the presence of papain and lysozyme that induce fast or slow enzymatic degradation. Due to this nondestructive and contactless approach, samples can be removed from the instrument and stored in different conditions between measurements. VeTBiM is well adapted to follow biomaterials alone or with cells, over long periods of time. This new method will help in the fine tuning of the mechanical properties of biomaterials used for cell therapy and tissue engineering.

A Robotic Ultrasound Scanner for Automatic Vessel Tracking and Three-Dimensional Reconstruction of B-Mode Images

Locating and evaluating the length and severity of a stenosis is very important for planning adequate treatment of peripheral arterial disease (PAD). Conventional ultrasound (US) examination cannot provide maps of entire lower limb arteries in 3-D. We propose a prototype 3D-US robotic system with B-mode images, which is nonionizing, noninvasive, and is able to track and reconstruct a continuous segment of the lower limb arterial tree between the groin and the knee. From an initialized cross-sectional view of the vessel, automatic tracking was conducted followed by 3D-US reconstructions evaluated using Hausdorff distance, cross-sectional area, and stenosis severity in comparison with 3-D reconstructions with computed tomography angiography (CTA). A mean Hausdorff distance of