Quantitative sparse array vascular elastography: the impact of tissue attenuation and modulus contrast on performance

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.

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.

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.

Effect of Clot Stiffness on Recombinant Tissue Plasminogen Activator Lytic Susceptibility in Vitro

The lytic recombinant tissue plasminogen activator (rt-PA) is the only drug approved by the Food and Drug Administration for treating ischemic stroke. Less than 40% of patients with large vessel occlusions who are treated with rt-PA have improved blood flow. However, up to 6% of all patients receiving rt-PA develop intracerebral hemorrhage. Predicting the efficacy of rt-PA treatment a priori could help guide therapeutic decision making, such that rt-PA is administered only to those individuals who would benefit from this treatment. Clot composition and structure affect the lytic efficacy of rt-PA and have an impact on elasticity. However, the relationship between clot elasticity and rt-PA lytic susceptibility has not been adequately investigated. The goal of this study was to quantify the relationship between clot elasticity and rt-PA susceptibility in vitro. Human and porcine highly retracted and mildly retracted clots were fabricated in glass pipettes. The rt-PA lytic susceptibility was evaluated in vitro using the percent clot mass loss. The Young's moduli of the clots were estimated using ultrasound-based single-track-location shear wave elasticity imaging. The percent mass loss in mildly retracted porcine and human clots (28.9 ± 6.1% and 45.2 ± 7.1%, respectively) was significantly higher (p < 0.05) than in highly retracted porcine and human clots (10.9 ± 2.1% and 25.5 ± 10.0%, respectively). Furthermore, the Young's moduli of highly retracted porcine and human clots (5.33 ± 0.92 and 3.21 ± 1.97 kPa, respectively) were significantly higher (p < 0.05) than those of mildly retracted porcine and human clots (2.66 ± 0.55 and 0.79 ± 0.21 kPa, respectively). The results revealed an inverse relationship between the percent clot mass loss and Young's modulus. These findings motivate continued investigation of ultrasound-based methods to assess clot stiffness in order to predict rt-PA thrombolytic efficacy.

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.

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.

Frequency-sum beamforming for passive cavitation imaging

Beamforming includes a variety of spatial filtering techniques that may be used for determining sound source locations from near-field sensor array recordings. For this scenario, beamforming resolution depends on the acoustic frequency, array geometry, and target location. Random scattering in the medium between the source and the array may degrade beamforming resolution with higher frequencies being more susceptible to degradation. The performance of frequency-sum (FS) beamforming for reducing such sensitivity to mild scattering while increasing resolution is reported here. FS beamforming was used with a data-dependent [minimum variance (MV)] or data-independent (delay-and-sum, DAS) weight vector to produce higher frequency information from lower frequency signal components via a quadratic product of complex signal amplitudes. The current findings and comparisons are based on simulations and passive cavitation imaging experiments using 3 MHz and 6 MHz emissions recorded by a 128-element linear array. FS beamforming results are compared to conventional DAS and MV beamforming using four metrics: point spread function (PSF) size, axial and lateral contrast, and computation time. FS beamforming produces a smaller PSF than conventional DAS beamforming with less computation time than MV beamforming in free space and mild scattering environments. However, it may fail when multiple unknown sound sources are present.

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

PURPOSE

Vascular elastography can visualize the strain distribution in the carotid artery, which could be useful in assessing the propensity of advanced plaques to rupture. In our previous studies, we demonstrated that sparse synthetic aperture (SA) imaging can produce high quality vascular strain elastograms. However, the low output power of SA imaging may hamper its clinical utility. In this study, we hypothesize that multi-element defocused emissions can overcome this limitation and improve the quality of the vascular strain elastograms.

METHODS

To assess the impact of attenuation on the elastographic performance of SA and (multi-element synthetic aperture) MSA imaging, we conducted experiments using heterogeneous vessel phantoms with ideal (0.1 dB cm^-1 MHz^-1 ) and realistic (0.75 dB cm^-1 MHz^-1 ) attenuation. Further, we validated the results of the phantom study in vivo, on a healthy male volunteer. All echo imaging was performed at a transmit frequency of 5 MHz, using a commercially available ultrasound scanner (Sonix RP, Ultrasonix Medical Corp., Richmond, BC, Canada).

RESULTS

The results from the phantom results demonstrated that plaques were visible in all strain elastograms, but those produced using MSA imaging had less artifacts. MSA imaging improved the elastographic contrast to noise ratio (CNRe) of the vascular elastograms by 14.58 dB relative to SA imaging, and 9.1 dB relative to compounded plane wave (CPW) imaging. Further, the results demonstrated that the elastographic performance of MSA imaging improved with increase in (a) the number of transmit-receive events and (b) the size of the transmit sub-aperture, up to 13 elements. Using larger sub-apertures degraded the elastographic performance. The results from the in vivo study were in good agreement with the phantom results. Specifically, using a defocused multi-element transmit sub-aperture for SA imaging improved the performance of vascular elastography.

CONCLUSIONS

The results suggested that MSA imaging can produce reliable vascular stain elastograms. Future studies will involve using coded excitations to improve the CNRe and frame-rate of the proposed technique for vascular elastography.

Principal Strain Vascular Elastography: Simulation and Preliminary Clinical Evaluation

It is difficult to produce reliable polar strain elastograms (radial and circumferential) because the center of the carotid artery is typically unknown. Principal strain imaging can overcome this limitation, but suboptimal lateral displacement estimates make this an impractical approach for visualizing mechanical properties within the carotid artery. We hypothesized that compounded plane wave imaging can minimize this problem. To test this hypothesis, we performed (i) simulations with vessels of varying morphology and mechanical behavior (i.e., isotropic and transversely isotropic), and (ii) a pilot study with 10 healthy volunteers. The accuracy of principal and polar strain (computed using knowledge of the precise vessel center) elastograms varied between 7% and 17%. In both types of elastograms, strain concentrated at the junction between the fibrous cap and the vessel wall, and the strain magnitude decreased with increasing fibrous cap thickness. Elastograms of healthy volunteers were consistent with those of transversely isotropic homogeneous vessels; they were spatially asymmetric, a trend that was common to both principal and polar strains. No significant differences were observed in the mean strain recovered from principal and polar strains (p > 0.05). This investigation indicates that principal strain elastograms measured with compounding plane wave imaging overcome the problems incurred when polar strain elastograms are computed with imprecise estimates of the vessel center.

Contrast-Enhanced Quantitative Intravascular Elastography: The Impact of Microvasculature on Model-Based Elastography

Model-based intravascular ultrasound elastography visualizes the stress distribution within vascular tissue-information that clinicians could use to predict the propensity of atherosclerotic plaque rupture. However, there are concerns that clusters of microvessels may reduce the accuracy of the estimated stress distribution. Consequently, we have developed a contrast-enhanced intravascular ultrasound system to investigate how plaque microvasculature affects the performance of model-based elastography. In simulations, diameters of 200, 400 and 800 μm were used, where the latter diameter represented a cluster of microvessels. In phantoms, we used a microvessel with a diameter of 750 μm. Peak stress errors of 3% and 38% were incurred in the fibrous cap when stress recovery was performed with and without a priori information about microvessel geometry. The results indicate that incorporating geometric information about plaque microvasculature obtained with contrast-enhanced ultrasound imaging improves the accuracy of estimates of the stress distribution within the fibrous cap precisely.