Noninvasive carotid artery elastography using multielement synthetic aperture imaging: Phantom and in vivo evaluation

Comparing Focused-Tracked and Plane Wave-Tracked ARFI Log(VoA) In Silico and in Application to Human Carotid Atherosclerotic Plaque, Ex Vivo

A significant risk factor for ischemic stroke is carotid atherosclerotic plaque that is susceptible to rupture, with rupture potential conveyed by plaque morphology. Human carotid plaque composition and structure have been delineated noninvasively and in vivo by evaluating log(VoA), a parameter derived as the decadic log of the second time derivative of displacement induced by an acoustic radiation force impulse (ARFI). In prior work, ARFI-induced displacement was measured using conventional focused tracking; however, this requires a long data acquisition period, thereby reducing framerate. We herein evaluate if ARFI log(VoA) framerate can be increased without a reduction in plaque imaging performance using plane wave tracking instead. In silico, both focused- and plane wave-tracked log(VoA) decreased with increasing echobrightness, quantified as signal-to-noise ratio (SNR), but did not vary with material elasticity for SNRs below 40 dB. For SNRs of 40-60 dB, both focused- and plane wave-tracked log(VoA) varied with SNR and material elasticity. Above 60 dB SNR, both focused- and plane wave-tracked log(VoA) varied with material elasticity alone. This suggests that log(VoA) discriminates features according to a combination of their echobrightness and mechanical property. Further, while both focused- and plane-wave tracked log(VoA) values were artifactually inflated by mechanical reflections at inclusion boundaries, plane wave-tracked log(VoA) was more strongly impacted by off-axis scattering. Applied to three excised human cadaveric carotid plaques with spatially aligned histological validation, both log(VoA) methods detected regions of lipid, collagen, and calcium (CAL) deposits. These findings support that plane wave tracking performs comparably to focused tracking for log(VoA) imaging and that plane wave-tracked log(VoA) is a viable approach to discriminating clinically relevant atherosclerotic plaque features at a 30-fold higher framerate than by focused tracking.

Deformability of ascending thoracic aorta aneurysms assessed using ultrafast ultrasound and a principal strain estimator: In vitro evaluation and in vivo feasibility

BACKGROUND

Noninvasive vascular strain imaging under conventional line-by-line scanning has a low frame rate and lateral resolution, and depends on the coordinate system. It is thus affected by high deformations due to image decorrelation between frames.

PURPOSE

To develop an ultrafast time-ensemble regularized tissue-Doppler optical-flow principal strain estimator for aorta deformability assessment in a long-axis view.

METHODS

This approach alleviated the impact of lateral resolution using image compounding and that of the coordinate system dependency using principal strain. Accuracy and feasibility were evaluated in two aorta-mimicking phantoms first, and then in four age-matched individuals with either a normal aorta or a pathological ascending thoracic aorta aneurysm (TAA).

RESULTS

Instantaneous aortic maximum and minimum principal strain maps and regional accumulated strains during each cardiac cycle were estimated at systolic and diastolic phases to characterize the normal aorta and TAA. In vitro, principal strain results matched sonomicrometry measurements. In vivo, a significant decrease in maximum and minimum principal strains was observed in TAA cases, whose range was respectively 7.9 ± 6.4% and 8.2 ± 2.6% smaller than in normal aortas.

CONCLUSIONS

The proposed principal strain estimator showed an ability to potentially assess TAA deformability, which may provide an individualized and reliable evaluation method for TAA rupture risk assessment. This article is protected by copyright. All rights reserved.

Quantitative assessment of ensemble coherency in contrast-free ultrasound microvasculature imaging

Purpose

Contrast‐free visualization of microvascular blood flow (MBF) using ultrasound can play a valuable role in diagnosis and detection of diseases. In this study, we demonstrate the importance of quantifying ensemble coherence for robust MBF imaging. We propose a novel approach to quantify ensemble coherence by estimating the local spatiotemporal correlation (LSTC) image, and evaluate its efficacy through simulation and in vivo studies.

Methods

The in vivo patient studies included three volunteers with a suspicious breast tumor, 15 volunteers with a suspicious thyroid tumor, and two healthy volunteers for renal MBF imaging. The breast data displayed negligible prior motion and were used for simulation analysis involving synthetically induced motion, to assess its impact on ensemble coherency and motion artifacts in MBF images. The in vivo thyroid data involved complex physiological motion due to its proximity to the pulsating carotid artery, which was used to assess the in vivo efficacy of the proposed technique. Further, in vivo renal MBF images demonstrated the feasibility of using the proposed ensemble coherence metric for curved array‐based MBF imaging involving phase conversion. All ultrasound data were acquired at high imaging frame rates and the tissue signal was suppressed using spatiotemporal clutter filtering. Thyroid tissue motion was estimated using two‐dimensional normalized cross correlation‐based speckle tracking, which was subsequently used for ensemble motion correction. The coherence of the MBF image was quantified based on Casorati correlation of the Doppler ensemble.

Results

The simulation results demonstrated that an increase in ensemble motion corresponded with a decrease in ensemble coherency, which reciprocally degraded the MBF images. Further the data acquired from breast tumors demonstrated higher ensemble coherency than that from thyroid tumors. Motion correction improved the coherence of the thyroid MBF images, which substantially improved its visualization. The proposed coherence metrics were also useful in assessing the ensemble coherence for renal MBF imaging. The results also demonstrated that the proposed coherence metric can be reliably estimated from downsampled ensembles (by up to 90%), thus allowing improved computational efficiency for potential applications in real‐time MBF imaging.

Conclusions

This pilot study demonstrates the importance of assessing ensemble coherency in contrast‐free MBF imaging. The proposed LSTC image quantified coherence of the Doppler ensemble for robust MBF imaging. The results obtained from this pilot study are promising, and warrant further development and in vivo validation.

Shear Induced Non-Linear Elasticity Imaging: Elastography for Compound Deformations

The goal of non-linear ultrasound elastography is to characterize tissue mechanical properties under finite deformations. Existing methods produce high contrast non-linear elastograms under conditions of pure uni-axial compression, but exhibit bias errors of 10-50% when the applied deformation deviates from the uni-axial condition. Since freehand transducer motion generally does not produce pure uniaxial compression, a motion-agnostic non-linearity estimator is desirable for clinical translation. Here we derive an expression for measurement of the Non-Linear Shear Modulus (NLSM) of tissue subject to combined shear and axial deformations. This method gives consistent nonlinear elasticity estimates irrespective of the type of applied deformation, with a reduced bias in NLSM values to 6-13%. The method combines quasi-static strain imaging with Single-Track Location-Shear Wave Elastography (STL-SWEI) to generate local estimates of axial strain, shear strain, and Shear Wave Speed (SWS). These local values were registered and non-linear elastograms reconstructed with a novel nonlinear shear modulus estimation scheme for general deformations. Results on tissue mimicking phantoms were validated with mechanical measurements and multiphysics simulations for all deformation types with an error in NLSM of 6-13%. Quantitative performance metrics with the new compound-motion tracking strategy reveal a 10-15 dB improvement in Signal-to-Noise Ratio (SNR) for simple shear versus pure compressive deformation for NLSM elastograms of homogeneous phantoms. Similarly, the Contrast-to-Noise Ratio (CNR) of NLSM elastograms of inclusion phantoms improved by 25-30% for simple shear over pure uni-axial compression. Our results show that high fidelity NLSM estimates may be obtained at ~30% lower strain under conditions of shear deformation as opposed axial compression. The reduction in strain required could reduce sonographer effort and improve scan safety.

A comparison between unfocused and focused transmit strategies in cardiac strain imaging

Unfocused ultrasound imaging, particularly coherent compounding with diverging waves, is a commonly employed high-frame rate transmit strategy in cardiac strain imaging. However, the accuracy and precision of diverging wave imaging compared to focused-beam transmit approaches in human subjects is unknown. Three transmit strategies-coherent compounding imaging, composite focused imaging with ECG gating and narrow-beams, and focused imaging with wide-beams-were compared in simulation and in transthoracic imaging of healthy human subjects (n  =  7). The focused narrow-beam sequence estimated radial end-systolic cumulative strains of a simulated left ventricular deformation with 26%  ±  1.5% and 34%  ±  1.5% greater accuracy compared with compounding and wide-beam imaging, respectively. Strain estimation precision in transthoracic imaging was then assessed with the Strain Filter on cumulative end-systolic radial strains. Within the strain values where statistically significant differences in precision (E(SNR_e|ε)) were found between transmit strategies, the narrow-beam sequence estimated radial strain 13%  ±  0.71% and 34%  ±  8.9% more precisely on average compared to compounding or wide-beam imaging, respectively.

Impact of imaging cross-section on visualization of thyroid microvessels using ultrasound: Pilot study

Non-invasive, contrast-free microvascular imaging of human thyroids can be potentially beneficial in reducing the large number of benign biopsies of suspicious nodules. However, motion incurred by thyroid due to its proximity to the pulsating carotid artery significantly impacts the visualization of blood flow in small vessels. Singular value based spatiotemporal clutter filtering (SVD-STF) improves the performance of tissue rejection in the presence of motion. However, despite effective clutter filtering, motion in thyroid imaging can impact coherent integration of the Doppler ensemble and degrade the visualization of the underlying vasculature. Recently studies have demonstrated that motion correction using 2D normalized cross-correlation based speckle tracking can address this issue, however, only in-plane motion can be tracked and corrected. Given the natural anatomical orientation of the rigid trachea, thyroid and the pulsating carotid artery, we hypothesize that imaging of thyroid microvessels may be more reliable in the longitudinal view than in the transverse. Specifically, distal presence of rigid trachea can limit out-of-plane motion in the longitudinal view. We tested this hypothesis on 48 acquisitions obtained from 24 thyroid patients having at least one suspicious nodule. In each patient, ultrasound images of the thyroid were acquired in both longitudinal and transverse views. Compounded plane-wave imaging was used to acquire the ultrasound images at high frame-rate, which is important for contrast-free small vessel blood flow imaging. Thyroid motion was tracked using 2D normalized cross-correlation based speckle tracking. Tissue clutter was rejected using singular value decomposition based spatiotemporal clutter filtering. The clutter-filtered Doppler ensemble was motion corrected prior to slow-time power Doppler integration. Signal-to-noise and contrast-to-noise ratios were computed to assess the improvement in quality of the power Doppler images. Out-of-plane motion was detected by estimating normalized ensemble cross-correlation coefficient. The results demonstrated that motion associated with the thyroid due to the carotid artery was primarily in the lateral direction, which could be estimated and corrected using 2D speckle tracking. However, the motion in the transverse view displayed increased speckle decorrelation. The average ensemble cross-correlation coefficient of the thyroid ultrasound images were significantly higher (p < 0.05) in the longitudinal view than in the transverse view. The largest improvement in SNR and CNR of the estimated PD images upon motion correction was observed in the longitudinal view (12.95 ± 3.76 dB and 16.48 ± 4.6 dB) than in the transverse view (3.72 ± 0.894 dB and 6.217 ± 1.689 dB). These preliminary results show that motion encountered by the thyroid due to carotid pulsations can be effectively tracked and corrected in the longitudinal view relative to transverse, which is important for reliably visualizing the underlying blood flow.

Adaptive background noise bias suppression in contrast-free ultrasound microvascular imaging

Non-invasive, contrast-free imaging of small vessel blood flow is diagnostically invaluable for detection, diagnosis and monitoring of disease. Recent advances in ultrafast imaging and tissue clutter-filtering have considerably improved the sensitivity of power Doppler (PD) imaging in detecting small vessel blood flow. However, suppression of tissue clutter exposes the depth-dependent time-gain compensated noise bias that noticeably degrades the PD image. We hypothesized that background suppression of PD images based on noise bias estimated from the entire clutter-filtered singular value spectrum can considerably improve flow signal visualization compared to currently existing techniques. To test our hypothesis, in vivo experiments were conducted on suspicious breast lesions in 10 subjects and deep-seated hepatic and renal microvasculatures in four healthy volunteers. Ultrasound PD images were acquired using a clinical ultrasound scanner, implemented with compounded plane wave imaging. The time gain compensated noise field was computed from the clutter-filtered Doppler ensemble (CFDE) based on its local spatio-temporal correlation, combined with low-rank signal estimation. Subsequently, the background bias in the PD images was suppressed by subtracting the estimated noise field. Background-suppressed PD images obtained using the proposed technique substantially improved visualization of the blood flow signal. The background bias in the noise suppressed PD images varied  <0.6 dB, independent of depth, which otherwise increased up to 13.8 dB. Further, the results demonstrated that the proposed technique efficaciously suppressed the background noise bias associated with smaller Doppler ensembles, which are challenging due to increased overlap between blood flow and noise components in the singular value spectrum. These preliminary results demonstrate the utility of the proposed technique to improve the visualization of small vessel blood flow in contrast-free PD images. The results of this feasibility study were encouraging, and warrant further development and additional in vivo validation.

Nonlinear Shear Modulus Estimation With Bi-Axial Motion Registered Local Strain

Nonlinear elasticity imaging provides additional information about tissue behavior that is potentially diagnostic and avoids errors inherent in applying a linear elastic model to tissue under large strains. Nonlinear elasticity imaging is challenging to perform due to the large deformations required to obtain sufficient tissue strain to elicit nonlinear behavior. This work uses a method of axial and lateral displacement tracking to estimate local axial strain with simultaneous measurement of shear modulus at multiple compression levels. By following the change in apparent shear modulus and the stress deduced from the strain maps, we are able to accurately quantify nonlinear shear modulus (NLSM). We have validated our technique with a mechanical NLSM measurement system. Our results demonstrate that 2-D tracking provides more consistent NLSM estimates than those obtained by 1-D (axial) tracking alone, especially where lateral motion is significant. The elastographic contrast-to-noise ratio in heterogeneous phantoms was 12.5%-60% higher using our method than that of 1-D tracking. Our method is less susceptible to mechanical variations, with deviations in mean elastic values of 2%-4% versus 5%-37% for 1-D tracking.

Distributing Synthetic Focusing Over Multiple Push-Detect Events Enhances Shear Wave Elasticity Imaging Performance

Plane wave (PW) imaging is a commonly used method for tracking waves during shear wave elasticity imaging (SWEI), but its unfocused transmission beam reduces tracking accuracy and precision. Coherent compounding minimizes this problem, but SWEI's stringent frame rate requirement and the coarse pitch of most clinical transducers limit its effectiveness. Synthetic aperture imaging (SAI) is an alternate ultrasound imaging approach with a much tighter focus than PW imaging, but its lower transmission power has deterred researchers from using SAI in SWEI. Hadamard-encoded multielement SAI can overcome this limitation. However, only a limited number of subapertures (3-5) can be transmitted in a single push-detect event. We have developed methods to distribute more subapertures or more compounding angles over multiple push-detect events. In this paper, we report the results of experiments conducted on phantoms to assess SWEI's performance when using Hadamard-encoded distributed-multielement synthetic aperture (HDMSA) imaging or distributed plane wave compounding (DPWC) to track shear waves. Tracking shear waves with HDMSA improved the elastographic signal-to-noise ratio (SNR_e) by 61.6%-89.5% depending on the phantom employed. Similarly, DPWC tracking improved SNR_e by 56.2%-93.3% for the same group of phantoms. Compared to focused ultrasound tracking (at the focus), SNR_e improved by 28.6% and 32.5% when tracking shear waves with HDMSA and DPWC, respectively. Long acquisitions could introduce decoding errors that decrease the performance when performing HDMSA tracking within the clinical setting. Nevertheless, the results of studies performed on the bicep muscle of three healthy volunteers demonstrate that for stationary organs, tracking shear waves with HDMSA yielded repeatable elastograms that offer better elastographic performance than those produced with current tracking methods.

Non-invasive Small Vessel Imaging of Human Thyroid Using Motion-Corrected Spatiotemporal Clutter Filtering

Reliable assessment of small vessel blood flow in the thyroid, without using any contrast agents, can be challenging because of increased physiological motion resulting from its proximity to the pulsating carotid artery. In this study, we hypothesized that correction of tissue motion prior to singular value decomposition (SVD)-based clutter filtering can improve the coherency of the tissue components and, thus, may allow better clutter suppression and visualization of small vessels in the thyroid. We corroborated this hypothesis by conducting phantom and in vivo studies using a clinical ultrasound scanner implemented with compounded plane wave imaging. The phantom studies were conducted using a homogeneous tissue-mimicking phantom to study the impact of motion on the covariance of the spatiotemporal Doppler data, in the absence of blood activity. The non-invasive in vivo study was conducted on a 74-y-old woman with a thyroid nodule suspicious of malignancy. A rigid body-based motion correction was performed using tissue displacements obtained from 2-D normalized cross-correlation-based speckle tracking. Subsequently, the power Doppler images were computed using SVD-based spatiotemporal clutter filtering. The results from the phantom study revealed that motion can considerably reduce the covariance of the spatiotemporal data and, thus, increase the rank of the tissue components. When the phantom was subjected to a total translation displacement of 6 pixels over the entire ensemble, in each direction (axial and lateral), the covariance dropped by more than 25%. The results obtained from the non-invasive in vivo study indicated that visualization of small vessel blood flow improved with motion correction of the power Doppler ensemble. The contrast-to-noise ratio of the blood signal in motion-corrected power Doppler images was considerably higher (8.17 and 8.32 dB), compared with that obtained using the standard SVD approach at an optimal threshold (0.87 and 4.33 dB) and a lower singular value threshold (1.92 and 3.05 dB). Further, the covariance of the in vivo thyroid spatiotemporal data increased by approximately 10% with motion correction. These preliminary results indicate that motion correction can be used to improve the visualization of small vessel blood flow in the thyroid, without using any contrast agents. The results of this feasibility study were encouraging, and warrant further development and more in vivo validation in moving tissues and organs.

Non-contrast agent based small vessel imaging of human thyroid using motion corrected power Doppler imaging

Singular value based spatiotemporal clutter filtering (SVD-STF) can significantly improve the sensitivity of blood flow imaging in small vessels without using contrast agents. However, despite effective clutter filtering, large physiological motion in thyroid imaging can impact coherent integration of the Doppler signal and degrade the visualization of the underlying vasculature. In this study, we hypothesize that motion correction of the clutter filtered Doppler ensemble, prior to the power Doppler estimation, can considerably improve the visualization of smalls vessels in suspicious thyroid nodules. We corroborated this hypothesis by conducting in vivo experiments on 10 female patients in the age group 44-82 yrs, with at least one thyroid nodule suspicious of malignancy, with recommendation for fine needle aspiration biopsy. Ultrasound images were acquired using a clinical ultrasound scanner, implemented with compounded plane wave imaging. Axial and lateral displacements associated with the thyroid nodules were estimated using 2D normalized cross-correlation. Subsequently, the tissue clutter associated with the Doppler ensemble was suppressed using SVD-STF. Motion correction of the clutter-filtered Doppler ensemble was achieved using a spline based sub-pixel interpolation. The results demonstrated that power Doppler images of thyroid nodules were noticeably degraded due to large physiological motion of the pulsating carotid artery in the proximity. The resultant power Doppler images were corrupted with signal distortion, motion blurring and occurrence of artificial shadow vessels and displayed visibly low signal-to-background contrast. In contrast, the power Doppler images obtained from the motion corrected ultrasound data addressed the issue and considerabley improved the visualization of blood flow. The signal-to-noise ratio and the contrast-to-noise ratio increased by up to 15.2 dB and 12.1 dB, respectively. Across the ten subjects, the highest improvement was observed for the nodule with the largest motion. These preliminary results show the ability of using motion correction to improve the visualization of small vessel blood flow in thyroid, without using any contrast agents. The results of this feasibility study were encouraging, and warrant further development and more in vivo validation in moving tissues and organs.

Diverging beam with synthetic aperture technique for rotation elastography: preliminary experimental results

Rotation Elastogram (RE) is a 2D spatial distribution map of the estimated local rigid-body rotation undergone by a target when subjected to an external compression, which is one of the recent variants in elastographic imaging. A recent study has shown that inclusion-contrast in RE is independent of inclusion-background modulus contrast and thus may be helpful in distinguishing between barely-stiff benign and malignant lesions. However, estimation of quality RE requires not only precise axial displacement estimates but also lateral displacement estimates. The widely used conventional focused beamforming technique using linear array (CFB-LA) provides better lateral resolution only over the depth of focus, which still results in poorer quality lateral displacement estimates compared to the axial displacement estimates. As an alternative to overcome this depth-dependent lateral resolution and obtain an improved lateral resolution, synthetic aperture-based approaches have been proposed in literature. Recently, we developed a synthetic aperture-based method, diverging beam with synthetic aperture technique (DB-SAT) that was aimed to not only reduce the ultrasound system complexity, but also provide improved lateral resolution throughout the depth of imaging and at higher frame-rate than that is possible in CFB-LA. In this paper, we report the preliminary experimental findings on the use of DB-SAT on RE and compare the resultant image quality against that obtained using often-employed CFB-LA and the synthetic transmit aperture (STA) technique. The investigation was done on tissue-mimicking phantoms and using contrast-to-noise ratio (CNR) as the metric for performance evaluation. The estimated CNR values from the REs obtained using CFB-LA, STA, and DB-SAT were 2.69  ±  0.81, 1.35  ±  0.22, and 14.71  ±  9.83, respectively, for inclusion present at 55 mm depth. The obtained results clearly demonstrated that the quality of RE can be improved significantly, especially at larger depth, using DB-SAT compared to that obtained using CFB-LA and STA technique.

Visualizing Angle-Independent Principal Strains in the Longitudinal View of the Carotid Artery: Phantom and In Vivo Evaluation

Non-invasive vascular elastography can evaluate the stiffness of the carotid artery by visualizing the vascular strain distribution. Axial strain estimates of the longitudinal cross section of the carotid artery are sensitive to the angle between the artery and the transducer. Anatomical variations in branching and arching of the carotid artery can affect the assessment of arterial stiffness. In this study, we hypothesized that principal strain elastograms computed using compounded plane wave imaging can reliably visualize the strain distribution in the carotid artery, independent of the transducer angle. We corroborated this hypothesis by conducting phantom and in vivo studies using a commercial ultrasound scanner (Sonix RP, Ultrasonix Medical Corp., Richmond, BC, Canada). The phantom studies were conducted using a homogeneous cryogel vessel phantom. The goal of the phantom study was to assess the feasibility of visualizing the radial deformation in the longitudinal plane of the vessel phantom, independent of the transducer angle (±30°, ±20°, ±10° and 0°). The in vivo studies were conducted on 20 healthy human volunteers in the age group 50-60 y. All echo imaging was performed at a transmit frequency of 5 MHz and sampling frequency of 40 MHz. The elastograms obtained from the phantom study revealed that for straight vessels, which had their lumen parallel to the transducer, principal strains were similar to axial strains. At non-parallel configurations (angles ±30°, ±20° and ±10°), the magnitudes of the mean principal strains were within 2.5% of the parallel configuration (0° angle) estimates and, thus, were observed to be relatively unaffected by change in angle. However, in comparison, the magnitude of the axial strain decreased with increase in angle because of coordinate dependency. Further, the pilot in vivo study indicated that the principal and axial strain elastograms were similar for subjects with relatively straight arteries. However, for arteries with arched geometry, axial strains were significantly lower (p <0.01) than the corresponding principal vascular strains, which was consistent with the results obtained from the phantom study. In conclusion, the results of the phantom and in vivo studies revealed that principal strain elastograms computed using CPW imaging could reliably visualize angle-independent vascular strains in the longitudinal plane of the carotid artery.

Optimization of Transmit Parameters in Cardiac Strain Imaging With Full and Partial Aperture Coherent Compounding

Coherent compounding methods using the full or partial transmit aperture have been investigated as a possible means of increasing strain measurement accuracy in cardiac strain imaging; however, the optimal transmit parameters in either compounding approach have yet to be determined. The relationship between strain estimation accuracy and transmit parameters-specifically the subaperture, angular aperture, tilt angle, number of virtual sources, and frame rate-in partial aperture (subaperture compounding) and full aperture (steered compounding) fundamental mode cardiac imaging was thus investigated and compared. Field II simulation of a 3-D cylindrical annulus undergoing deformation and twist was developed to evaluate accuracy of 2-D strain estimation in cross-sectional views. The tradeoff between frame rate and number of virtual sources was then investigated via transthoracic imaging in the parasternal short-axis view of five healthy human subjects, using the strain filter to quantify estimation precision. Finally, the optimized subaperture compounding sequence (25-element subperture, 90° angular aperture, 10 virtual sources, 300-Hz frame rate) was compared to the optimized steered compounding sequence (60° angular aperture, 15° tilt, 10 virtual sources, 300-Hz frame rate) via transthoracic imaging of five healthy subjects. Both approaches were determined to estimate cumulative radial strain with statistically equivalent precision (subaperture compounding E(SNRe %) = 3.56, and steered compounding E(SNRe %) = 4.26).