A high-frame-rate ultrasound system for the study of tissue motions

Multi-Frequency 3D Shear Wave Absolute Vibro-Elastography (S-WAVE) System for the Prostate

This article describes a novel system for quantitative and volumetric measurement of tissue elasticity in the prostate using simultaneous multi-frequency tissue excitation. Elasticity is computed by using a local frequency estimator to measure the three-dimensional local wavelengths of steady-state shear waves within the prostate gland. The shear wave is created using a mechanical voice coil shaker which transmits simultaneous multi-frequency vibrations transperineally. Radio frequency data is streamed directly from a BK Medical 8848 transrectal ultrasound transducer to an external computer where tissue displacement due to the excitation is measured using a speckle tracking algorithm. Bandpass sampling is used that eliminates the need for an ultra-fast frame rate to track the tissue motion and allows for accurate reconstruction at a sampling frequency that is below the Nyquist rate. A roll motor with computer control is used to rotate the transducer and obtain 3D data. Two commercially available phantoms were used to validate both the accuracy of the elasticity measurements as well as the functional feasibility of using the system for in vivo prostate imaging. The phantom measurements were compared with 3D Magnetic Resonance Elastography (MRE), where a high correlation of 96% was achieved. In addition, the system has been used in two separate clinical studies as a method for cancer identification. Qualitative and quantitative results of 11 patients from these clinical studies are presented here. Furthermore, an AUC of 0.87±0.12 was achieved for malignant vs. benign classification using a binary support vector machine classifier trained with data from the latest clinical study with leave one patient out cross-validation.

3-D Transducer Mounted Shear Wave Absolute Vibro-Elastography: Proof of Concept

Quantitative tissue stiffness characterization using ultrasound (US) has been shown to improve prostate cancer (PCa) detection in multiple studies. Shear wave absolute vibro-elastography (SWAVE) allows quantitative and volumetric assessment of tissue stiffness using external multifrequency excitation. This article presents a proof of concept of a first-of-a-kind 3-D hand-operated endorectal SWAVE system designed to be used during systematic prostate biopsy. The system is developed with a clinical US machine, requiring only an external exciter that can be mounted directly to the transducer. Subsector acquisition of radio frequency (RF) data allows imaging of shear waves with a high effective frame rate (up to 250 Hz). The system was characterized using eight different quality assurance phantoms. Due to the invasive nature of prostate imaging, at this early stage of development, validation of in vivo human tissue was instead carried out by intercostally scanning the livers of n = 7 healthy volunteers. The results are compared with 3-D magnetic resonance elastography (MRE) and an existing 3-D SWAVE system with a matrix array transducer (M-SWAVE). High correlations were found with MRE (99% in phantoms, 94% in liver data) and with M-SWAVE (99% in phantoms, 98% in liver data).

SWAVE 2.0 Imaging of Placental Elasticity and Viscosity: Potential Biomarkers for Placenta-Mediated Disease Detection

Pregnancy complications such as pre-eclampsia (PE) and intrauterine growth restriction (IUGR) are associated with structural and functional changes in the placenta. Different elastography techniques with an ability to assess the mechanical properties of tissue can identify and monitor the pathological state of the placenta. Currently available elastography techniques have been used with promising results to detect placenta abnormalities; however, limitations include inadequate measurement depth and safety concerns from high negative pressure pulses. Previously, we described a shear wave absolute vibro-elastography (SWAVE) method by applying external low-frequency mechanical vibrations to generate shear waves and studied 61 post-delivery clinically normal placentas to explore the feasibility of SWAVE for placental assessment and establish a measurement baseline. This next phase of the study, namely, SWAVE 2.0, improves the previous system and elasticity reconstruction by incorporating a multi-frequency acquisition system and using a 3-D local frequency estimation (LFE) method. Compared with its 2-D counterpart, the proposed system using 3-D LFE was found to reduce the bias and variance in elasticity measurements in tissue-mimicking phantoms. In the aim of investigating the potential of improved SWAVE 2.0 measurements to identify placental abnormalities, we studied 46 post-delivery placentas, including 26 diseased (16 IUGR and 10 PE) and 20 normal control placentas. By use of a 3.33-MHz motorized curved-array transducer, multi-frequency (80,100 and 120 Hz) elasticity measures were obtained with 3-D LFE, and both IUGR (15.30 ± 2.96 kPa, p = 3.35e^-5) and PE (12.33 ± 4.88 kPa, p = 0.017) placentas were found to be significantly stiffer compared with the control placentas (8.32 ± 3.67 kPa). A linear discriminant analysis (LDA) classifier was able to classify between healthy and diseased placentas with a sensitivity, specificity and accuracy of 87%, 78% and 83% and an area under the receiver operating curve of 0.90 (95% confidence interval: 0.8-0.99). Further, the pregnancy outcome in terms of neonatal intensive care unit admission was predicted with a sensitivity, specificity and accuracy of 70%, 71%, 71%, respectively, and area under the receiver operating curve of 0.78 (confidence interval: 0.62-0.93). A viscoelastic characterization of placentas using a fractional rheological model revealed that the viscosity measures in terms of viscosity parameter n were significantly higher in IUGR (2.3 ± 0.21) and PE (2.11 ± 0.52) placentas than in normal placentas (1.45 ± 0.65). This work illustrates the potential relevance of elasticity and viscosity imaging using SWAVE 2.0 as a non-invasive technology for detection of placental abnormalities and the prediction of pregnancy outcomes.

Breast Cancer Detection using Multimodal Time Series Features from Ultrasound Shear Wave Absolute Vibro-Elastography

In shear wave absolute vibro-elastography (S-WAVE), a steady-state multi-frequency external mechanical excitation is applied to tissue, while a time-series of ultrasound radio-frequency (RF) data are acquired. Our objective is to determine the potential of S-WAVE to classify breast tissue lesions as malignant or benign. We present a new processing pipeline for feature-based classification of breast cancer using S-WAVE data, and we evaluate it on a new data set collected from 40 patients. Novel bi-spectral and Wigner spectrum features are computed directly from the RF time series and are combined with textural and spectral features from B-mode and elasticity images. The Random Forest permutation importance ranking and the Quadratic Mutual Information methods are used to reduce the number of features from 377 to 20. Support Vector Machines and Random Forest classifiers are used with leave-one-patient-out and Monte Carlo cross-validations. Classification results obtained for different feature sets are presented. Our best results (95% confidence interval, Area Under Curve = 95%1.45%, sensitivity = 95%, and specificity = 93%) outperform the state-of-the-art reported S-WAVE breast cancer classification performance. The effect of feature selection and the sensitivity of the above classification results to changes in breast lesion contours is also studied. We demonstrate that time-series analysis of externally vibrated tissue as an elastography technique, even if the elasticity is not explicitly computed, has promise and should be pursued with larger patient datasets. Our study proposes novel directions in the field of elasticity imaging for tissue classification.

Echocardiography image enhancement using texture-cartoon separation

Due to the speckled nature of cardiac ultrasound imaging, it is not easy to process and extract useful information directly from the acquired image. In this work, we have proposed a method to reduce the effect of speckle artifacts through the decomposition of echocardiography images into cartoon and texture components. The first component (i.e., cartoon image) contains image structures containing smooth areas and sharp edges, and the texture component is mainly composed of highly oscillating and repetitive patterns. To decompose the image into these two subcomponents, convolutional sparse coding has been utilized as a solid tool for solving the decomposition optimization function. The significant advantage of using convolutional sparse coding, compared to classical sparse coding methods, is image quality enhancement due to not using the block coding, making the classic solutions computationally feasible. The original image has been masked with the cartoon part leading to suppress speckle artifacts which result in image quality enhancement. Besides, it has been shown that using this speckle reduction scenario, considerable accuracy enhancement of the segmentation task can be achieved, compared to segmentation of the original image. Numerical results provide acceptable reasons to prove the efficiency of the proposed algorithm. Resulting echocardiography videos show a mean segmentation enhancement of 15.98 for Hausdorff distance (in pixels) and 0.0632 for the Dice similarity coefficient.

Three-Dimensional Multi-Frequency Shear Wave Absolute Vibro-Elastography (3D S-WAVE) With a Matrix Array Transducer: Implementation and Preliminary In Vivo Study of the Liver

Magnetic resonance elastography (MRE) is commonly regarded as the imaging-based gold-standard for liver fibrosis staging, comparable to biopsy. While ultrasound-based elastography methods for liver fibrosis staging have been developed, they are confined to a 1D or a 2D region of interest and to a limited depth. 3D Shear Wave Absolute Vibro-Elastography (S-WAVE) is a steady-state, external excitation, volumetric elastography technique that is similar to MRE, but has the additional advantage of multi-frequency excitation. We present a novel ultrasound matrix array implementation of S-WAVE that takes advantage of 3D imaging. We use a matrix array transducer to sample axial multi-frequency steady-state tissue motion over a volume, using a Color Power Angiography sequence. Tissue motion with the frequency components {40,50,60} and {45,55,65} Hz are acquired over a (90° lateral) × (40° elevational) × (16 cm depth) sector with an acquisition time of 12 seconds. We compute the elasticity map in 3D using local spatial frequency estimation. We characterize this new approach in tissue phantoms against measurements obtained with transient elastography and MRE. Six healthy volunteers and eight patients with chronic liver disease were imaged. Their MRE and S-WAVE volumes were aligned using T1 to B-mode registration for direct comparison in common regions of interest. S-WAVE and MRE results are correlated with R² = 0.92, while MRE and TE results are correlated with R² = 0.71. Our findings show that S-WAVE with matrix array has the potential to deliver a similar assessment of liver fibrosis as MRE in a more accessible, inexpensive way, to a broader set of patients.

SWAVE Imaging of Placental Elasticity and Viscosity: Proof of Concept

The placenta is the interface between the fetus and the mother and is vital for fetal development. Ultrasound elastography provides a non-invasive way to examine in vivo the stiffness of the placenta; increased stiffness has previously been linked to fetal growth restriction. This study used a previously developed dynamic elastography method, called shear wave absolute vibro-elastography, to study 61 post-delivery clinically normal placentas. The shear wave speeds in the placenta were recorded under five different low-frequency mechanical excitations. The elasticity and viscosity were estimated through rheological modeling. The shear wave speeds at excitation frequencies of 60, 80, 90, 100 and 120 Hz were measured to be 1.23 ± 0.44, 1.67 ± 0.76, 1.74 ± 0.72, 1.80 ± 0.78 and 2.25 ± 0.80 m/s. The shear wave speed values we obtained are consistent with previous studies. In addition, our multi-frequency acquisition approach enables us to provide viscosity estimates that have not been previously reported.

Shear Wave Imaging of Breast Tissue by Color Doppler Shear Wave Elastography

Shear wave elastography is a distinctive method to access the viscoelastic characteristic of the soft tissue that is difficult to obtain by other imaging modalities. This paper proposes a novel shear wave elastography [color Doppler shear wave imaging (CD SWI)] for breast tissue. Continuous shear wave is produced by a small lightweight actuator, which is attached to the tissue surface. Shear wave wavefront that propagates in tissue is reconstructed as a binary pattern that consists of zero and the maximum flow velocities on color flow image (CFI). Neither any modifications of the ultrasound color flow imaging instrument nor a high frame rate ultrasound imaging instrument is required to obtain the shear wave wavefront map. However, two conditions of shear wave displacement amplitude and shear wave frequency are needed to obtain the map. However, these conditions are not severe restrictions in breast imaging. This is because the minimum displacement amplitude is [Formula: see text] for an ultrasonic wave frequency of 12 MHz and the shear wave frequency is available from several frequencies suited for breast imaging. Fourier analysis along time axis suppresses clutter noise in CFI. A directional filter extracts shear wave, which propagates in the forward direction. Several maps, such as shear wave phase, velocity, and propagation maps, are reconstructed by CD SWI. The accuracy of shear wave velocity measurement is evaluated for homogeneous agar gel phantom by comparing with the acoustic radiation force impulse method. The experimental results for breast tissue are shown for a shear wave frequency of 296.6 Hz.

Two-dimensional shear-wave elastography on conventional ultrasound scanners with time-aligned sequential tracking (TAST) and comb-push ultrasound shear elastography (CUSE)

Two-dimensional shear-wave elastography presents 2-D quantitative shear elasticity maps of tissue, which are clinically useful for both focal lesion detection and diffuse disease diagnosis. Realization of 2-D shear-wave elastography on conventional ultrasound scanners, however, is challenging because of the low tracking pulse-repetition-frequency (PRF) of these systems. Although some clinical and research platforms support software beamforming and plane-wave imaging with high PRF, the majority of current clinical ultrasound systems do not have the software beamforming capability, which presents a critical challenge for translating the 2-D shear-wave elastography technique from laboratory to clinical scanners. To address this challenge, this paper presents a time-aligned sequential tracking (TAST) method for shear-wave tracking on conventional ultrasound scanners. TAST takes advantage of the parallel beamforming capability of conventional systems and realizes high-PRF shear-wave tracking by sequentially firing tracking vectors and aligning shear wave data in the temporal direction. The comb-push ultrasound shear elastography (CUSE) technique was used to simultaneously produce multiple shear wave sources within the field-of-view (FOV) to enhance shear wave SNR and facilitate robust reconstructions of 2-D elasticity maps. TAST and CUSE were realized on a conventional ultrasound scanner. A phantom study showed that the shear-wave speed measurements from the conventional ultrasound scanner were in good agreement with the values measured from other 2-D shear wave imaging technologies. An inclusion phantom study showed that the conventional ultrasound scanner had comparable performance to a state-of-the-art shear-wave imaging system in terms of bias and precision in measuring different sized inclusions. Finally, in vivo case analysis of a breast with a malignant mass, and a liver from a healthy subject demonstrated the feasibility of using the conventional ultrasound scanner for in vivo 2-D shear-wave elastography. These promising results indicate that the proposed technique can enable the implementation of 2-D shear-wave elastography on conventional ultrasound scanners and potentially facilitate wider clinical applications with shear-wave elastography.

Three-dimensional sheaf of ultrasound planes reconstruction (SOUPR) of ablated volumes

This paper presents an algorithm for 3-D reconstruction of tumor ablations using ultrasound shear wave imaging with electrode vibration elastography. Radio-frequency ultrasound data frames are acquired over imaging planes that form a subset of a sheaf of planes sharing a common axis of intersection. Shear wave velocity is estimated separately on each imaging plane using a piecewise linear function fitting technique with a fast optimization routine. An interpolation algorithm then computes velocity maps on a fine grid over a set of C-planes that are perpendicular to the axis of the sheaf. A full 3-D rendering of the ablation can then be created from this stack of C-planes; hence the name "Sheaf Of Ultrasound Planes Reconstruction" or SOUPR. The algorithm is evaluated through numerical simulations and also using data acquired from a tissue mimicking phantom. Reconstruction quality is gauged using contrast and contrast-to-noise ratio measurements and changes in quality from using increasing number of planes in the sheaf are quantified. The highest contrast of 5 dB is seen between the stiffest and softest regions of the phantom. Under certain idealizing assumptions on the true shape of the ablation, good reconstruction quality while maintaining fast processing rate can be obtained with as few as six imaging planes suggesting that the method is suited for parsimonious data acquisitions with very few sparsely chosen imaging planes.

Real-time 1-D/2-D transient elastography on a standard ultrasound scanner using mechanically induced vibration

Transient elastography has been well established in the literature as a means of assessing the elasticity of soft tissue. In this technique, tissue elasticity is estimated from the study of the propagation of the transient shear waves induced by an external or internal source of vibration. Previous studies have focused mainly on custom single-element transducers and ultrafast scanners which are not available in a typical clinical setup. In this work, we report the design and implementation of a transient elastography system on a standard ultrasound scanner that enables quantitative assessment of tissue elasticity in real-time. Two new custom imaging modes are introduced that enable the system to image the axial component of the transient shear wave, in response to an externally induced vibration, in both 1-D and 2-D. Elasticity reconstruction algorithms that estimate the tissue elasticity from these transient waves are also presented. Simulation results are provided to show the advantages and limitations of the proposed system. The performance of the system is also validated experimentally using a commercial elasticity phantom.

Single-heartbeat electromechanical wave imaging with optimal strain estimation using temporally unequispaced acquisition sequences

Electromechanical Wave Imaging (EWI) is a non-invasive, ultrasound-based imaging method capable of mapping the electromechanical wave (EW) in vivo, i.e. the transient deformations occurring in response to the electrical activation of the heart. Optimal imaging frame rates, in terms of the elastographic signal-to-noise ratio, to capture the EW cannot be achieved due to the limitations of conventional imaging sequences, in which the frame rate is low and tied to the imaging parameters. To achieve higher frame rates, EWI is typically performed by combining sectors acquired during separate heartbeats, which are then combined into a single view. However, the frame rates achieved remain potentially sub-optimal and this approach precludes the study of non-periodic arrhythmias. This paper describes a temporally unequispaced acquisition sequence (TUAS) for which a wide range of frame rates are achievable independently of the imaging parameters, while maintaining a full view of the heart at high beam density. TUAS is first used to determine the optimal frame rate for EWI in a paced canine heart in vivo and then to image during ventricular fibrillation. These results indicate how EWI can be optimally performed within a single heartbeat, during free breathing and in real time, for both periodic and non-periodic cardiac events.

Remote ultrasound palpation for robotic interventions using absolute elastography

Although robotic surgery has addressed many of the challenges presented by minimally invasive surgery, haptic feedback and the lack of knowledge of tissue stiffness is an unsolved problem. This paper presents a system for finding the absolute elastic properties of tissue using a freehand ultrasound scanning technique, which utilizes the da Vinci Surgical robot and a custom 2D ultrasound transducer for intraoperative use. An external exciter creates shear waves in the tissue, and a local frequency estimation method computes the shear modulus. Results are reported for both phantom and in vivo models. This system can be extended to any 6 degree-of-freedom tracking method and any 2D transducer to provide real-time absolute elastic properties of tissue.

Electromechanical wave imaging for arrhythmias

Electromechanical wave imaging (EWI) is a novel ultrasound-based imaging modality for mapping of the electromechanical wave (EW), i.e. the transient deformations occurring in immediate response to the electrical activation. The correlation between the EW and the electrical activation has been established in prior studies. However, the methods used previously to map the EW required the reconstruction of images over multiple cardiac cycles, precluding the application of EWI for non-periodic arrhythmias such as fibrillation. In this study, new imaging sequences are developed and applied based on flash- and wide-beam emissions to image the entire heart at very high frame rates (2000 fps) during free breathing in a single heartbeat. The methods are first validated by imaging the heart of an open-chest canine while simultaneously mapping the electrical activation using a 64-electrode basket catheter. Feasibility is then assessed by imaging the atria and ventricles of closed-chest, conscious canines during sinus rhythm and during right-ventricular pacing following atrio-ventricular dissociation, i.e., during a non-periodic rhythm. The EW was validated against electrode measurements in the open-chest case, and followed the expected electrical propagation pattern in the closed-chest setting. These results indicate that EWI can be used for the characterization of non-periodic arrhythmias in conditions similar to the clinical setting, in a single heartbeat, and during free breathing.

Bandpass sampling of high-frequency tissue motion

The characterization of tissue viscoelastic properties requires the measurement of tissue motion over a region of interest at frequencies that significantly exceed the frame rates of conventional ultrasound systems. In this paper, we propose that the bandpass sampling technique be applied to tissue motion sampling. With this approach, high-frequency signals limited to a frequency band can be sampled and reconstructed without aliasing at a sampling frequency that is lower than the Nyquist rate. We first review this approach and discuss the selection of the tissue excitation frequency band and of the feasible sampling frequencies that allow signal reconstruction without aliasing. We then demonstrate the approach using simulations based on the finite element method and ultrasound simulations. Finally, we perform experiments on tissue-mimicking materials and demonstrate accurate motion estimation using a lower sampling rate than that required by the conventional sampling theorem. The estimated displacements were used to measure the elasticity and viscosity in a phantom in which an inclusion has been correctly delineated. Thus, with bandpass sampling, it is feasible to use conventional beamforming on diagnostic ultrasound systems to perform high-frequency dynamic elastography. The method is simple to implement because it does not require beam interleaving, additional hardware, or synchronization.

Frequency-locked pulse sequencer for high-frame-rate monochromatic tissue motion imaging

To overcome the inherent low frame rate of conventional ultrasound, we have previously presented a system that can be implemented on conventional ultrasound scanners for high-frame-rate imaging of monochromatic tissue motion. The system employs a sector subdivision technique in the sequencer to increase the acquisition rate. To eliminate the delays introduced during data acquisition, a motion phase correction algorithm has also been introduced to create in-phase displacement images. Previous experimental results from tissue- mimicking phantoms showed that the system can achieve effective frame rates of up to a few kilohertz on conventional ultrasound systems. In this short communication, we present a new pulse sequencing strategy that facilitates high-frame-rate imaging of monochromatic motion such that the acquired echo signals are inherently in-phase. The sequencer uses the knowledge of the excitation frequency to synchronize the acquisition of the entire imaging plane to that of an external exciter. This sequencing approach eliminates any need for synchronization or phase correction and has applications in tissue elastography, which we demonstrate with tissue-mimicking phantoms.

2-D high-frame-rate dynamic elastography using delay compensated and angularly compounded motion vectors: preliminary results

This paper describes a new ultrasound-based system for high-frame-rate measurement of periodic motion in 2-D for tissue elasticity imaging. Similarly to conventional 2-D flow vector imaging, the system acquires the RF signals from the region of interest at multiple steering angles. A custom sector subdivision technique is used to increase the temporal resolution while keeping the total acquisition time within the range suitable for real-time applications. Within each sector, 1-D motion is estimated along the beam direction. The intra- and inter-sector delays are compensated using our recently introduced delay compensation algorithm. In-plane 2-D motion vectors are then reconstructed from these delay-compensated 1-D motions. We show that Young's modulus images can be reconstructed from these 2-D motion vectors using local inversion algorithms. The performance of the system is validated quantitatively using a commercial flow phantom and a commercial elasticity phantom. At the frame rate of 1667 Hz, the estimated flow velocities with the system are in agreement with the velocity measured with a pulsed-wave Doppler imaging mode of a commercial ultrasound machine with manual angle correction. At the frame rate of 1250 Hz, phantom Young's moduli of 29, 6, and 54 kPa for the background, the soft inclusion, and the hard inclusion, are estimated to be 30, 11, and 53 kPa, respectively.