Specific Ultrasound Data Acquisition for Tissue Motion and Strain Estimation: Initial Results

3D normalized cross-correlation for estimation of the displacement field in ultrasound elastography

This paper introduces a novel technique to estimate tissue displacement in quasi-static elastography. A major challenge in elastography is estimation of displacement (also referred to time-delay estimation) between pre-compressed and post-compressed ultrasound data. Maximizing normalized cross correlation (NCC) of ultrasound radio-frequency (RF) data of the pre- and post-compressed images is a popular technique for strain estimation due to its simplicity and computational efficiency. Several papers have been published to increase the accuracy and quality of displacement estimation based on NCC. All of these methods use 2D spatial windows in RF data to estimate NCC, wherein displacement is assumed to be constant within each window. In this work, we extend this assumption along the third dimension. Two approaches are proposed to get third dimension. In the first approach, we use temporal domain to exploit neighboring samples in both spatial and temporal directions. Considering temporal information is important since traditional and ultrafast ultrasound machines are, respectively, capable of imaging at more than 30 frame per second (fps) and 1000 fps. Another approach is to use time-delayed pre-beam formed data (channel data) instead of RF data. In this method information of all channels that are recorded as pre-beam formed data of each RF line will be considered as 3^rd dimension. We call these methods as spatial temporal normalized cross correlation (STNCC) and channel data normalized cross correlation (CNCC) and show that they substantially outperforms NCC using simulation, phantom and in-vivo experiments. Given substantial improvements of results in addition to the relative simplicity of implementing STNCC and CNCC, the proposed approaches can potentially have a large impact in both academic and commercial work on ultrasound elastography.

Investigation of out-of-plane motion artifacts in 2D noninvasive vascular ultrasound elastography

Ultrasound noninvasive vascular elastography (NIVE) has shown its potential to measure strains of carotid arteries to predict plaque instability. When two-dimensional (2D) strain estimation is performed, either in longitudinal or cross-sectional view, only in-plane motions are considered. The motions in elevation direction (i.e. perpendicular to the imaging plane), can induce estimation artifacts affecting the accuracy of 2D NIVE. The influence of such out-of-plane motions on the performance of axial strain and axial shear strain estimations has been evaluated in this study. For this purpose, we designed a diseased carotid bifurcation phantom with a 70% stenosis and an in vitro experimental setup to simulate orthogonal out-of-plane motions of 1 mm, 2 mm and 3 mm. The Lagrangian speckle model estimator (LSME) was used to estimate axial strains and shears under pulsatile conditions. As anticipated, in vitro results showed more strain estimation artifacts with increasing magnitudes of motions in elevation. However, even with an out-of-plane motion of 2.0 mm, strain and shear estimations having inter-frame correlation coefficients higher than 0.85 were obtained. To verify findings of in vitro experiments, a clinical LSME dataset obtained from 18 participants with carotid artery stenosis was used. Deduced out-of-plane motions (ranging from 0.25 mm to 1.04 mm) of the clinical dataset were classified into three groups: small, moderate and large elevational motions. Clinical results showed that pulsatile time-varying strains and shears remained reproducible for all motion categories since inter-frame correlation coefficients were higher than 0.70, and normalized cross-correlations (NCC) between radiofrequency (RF) images were above 0.93. In summary, the performance of LSME axial strain and shear estimations appeared robust in the presence of out-of-plane motions (<2 mm) as encountered during clinical ultrasound imaging.

Electromagnetic tracking-based freehand 3D quasi-static elastography with 1D linear array: a phantom study

Recent developments in hardware and scanning protocols have advanced conventional 2D quasi-static elastography to 3D level, which provides an intuitive visualization of lesions. A 2D linear array or scanning mechanism is typically required for 3D quasi-static elastography, requiring expensive and specifically designed hardware. In this study, we propose a novel method based on a commercial electromagnetic tracking system for freehand 3D quasi-static elastography with 1D linear array. Phantom experiments are performed to validate the feasibility of the proposed method. During data acquisition, the probe contacts the surface of an elasticity phantom and moves in the elevational direction, while applying sinusoidal-like axial compression to the phantom. For each frame of ultrasound data, the 3D coordinates and orientations of the probe are obtained from an electromagnetic tracking system. A correlation-based algorithm is adopted to obtain a series of axial strain images. Volumetric strain data are reconstructed by using the recorded 3D coordinates and orientations of the probe corresponding to each strain image. The diameters of inclusions are then obtained from the slice plots of the volumetric strain data. The volumes of inclusions are estimated from the isosurface plots. The experimental result shows that the volume estimation of the inclusions has good accuracy, with errors within 2%, while the diameters of the inclusions estimated from three orthogonal planes have larger errors up to 18%. In conclusion, the present framework would promise a reliable and effective solution for freehand 3D quasi-static elastography with 1D linear array.