Principal Strain Vascular Elastography: Simulation and Preliminary Clinical Evaluation

Characterizing atherosclerotic tissues: in silico analysis of mechanical properties using intravascular ultrasound and inverse finite element methods

Atherosclerosis is a prevalent cause of acute coronary syndromes that consists of lipid deposition inside the artery wall, creating an atherosclerotic plaque. Early detection may prevent the risk of plaque rupture. Nowadays, intravascular ultrasound (IVUS) is the most common medical imaging technology for atherosclerotic plaque detection. It provides an image of the section of the coronary wall and, in combination with new techniques, can estimate the displacement or strain fields. From these magnitudes and by inverse analysis, it is possible to estimate the mechanical properties of the plaque tissues and their stress distribution. In this paper, we presented a methodology based on two approaches to characterize the mechanical properties of atherosclerotic tissues. The first approach estimated the linear behavior under particular pressure. In contrast, the second technique yielded the non-linear hyperelastic material curves for the fibrotic tissues across the complete physiological pressure range. To establish and validate this method, the theoretical framework employed in silico models to simulate atherosclerotic plaques and their IVUS data. We analyzed different materials and real geometries with finite element (FE) models. After the segmentation of the fibrotic, calcification, and lipid tissues, an inverse FE analysis was performed to estimate the mechanical response of the tissues. Both approaches employed an optimization process to obtain the mechanical properties by minimizing the error between the radial strains obtained from the simulated IVUS and those achieved in each iteration. The second methodology was successfully applied to five distinct real geometries and four different fibrotic tissues, getting median R 2 of 0.97 and 0.92, respectively, when comparing the real and estimated behavior curves. In addition, the last technique reduced errors in the estimated plaque strain field by more than 20% during the optimization process, compared to the former approach. The findings enabled the estimation of the stress field over the hyperelastic plaque tissues, providing valuable insights into its risk of rupture.

Reconstructing the Spatial Distribution of the Relative Shear Modulus in Quasi-static Ultrasound Elastography: Plane Stress Analysis

Quasi-static ultrasound elastography (QSUE) is an imaging technique that mainly provides axial strain maps of tissues when the latter are subjected to compression. In this article, a method for reconstructing the relative shear modulus distribution within a linear elastic and isotropic medium, in QSUE, is introduced. More specifically, the plane stress inverse problem is considered. The proposed method is based on the variational formulation of the equilibrium equations and on the choice of adapted discretization spaces, and only requires displacement fields in the analyzed media to be determined. Results from plane stress and 3-D numerical simulations, as well as from phantom experiments, showed that the method is able to reconstruct the different regions within a medium, with shear modulus contrasts that unambiguously reveal whether inclusions are stiffer or softer than the surrounding material. More specifically, for the plane stress simulations, inclusion-to-background modulus ratios were found to be very accurately estimated, with an error lower than 3%. For the 3-D simulations, for which the plane stress conditions are no longer satisfied, these ratios were, as expected, less accurate, with an error that remained lower than 10% for two of the three cases analyzed but was around 34% for the last case. Concerning the phantom experiments, a comparison with a shear wave elastography technique from a clinical ultrasound scanner was also made. Overall, the inclusion-to-background shear modulus ratios obtained with our approach were found to be closer to those given by the phantom manufacturer than the ratios provided by the clinical system.

Ultrafast Myocardial Principal Strain Ultrasound Elastography During Stress Tests: In Vitro Validation and In Vivo Feasibility

Objective myocardial contractility assessment during stress tests aims to improve the diagnosis of myocardial ischemia. Tissue Doppler imaging (TDI) or optical flow (OF) speckle tracking echocardiography (STE) has been used to quantify myocardial contractility at rest. However, this is more challenging during stress tests due to image decorrelation at high heart rates. Moreover, stress tests imply a high frame rate which leads to a limited lateral field of view. Therefore, a large lateral field-of-view robust ultrafast myocardial regularized OF-TDI principal strain estimator has been developed for high-frame-rate echocardiography of coherently compounded transmitted diverging waves. The feasibility and accuracy of the proposed estimator were validated in vitro (using sonomicrometry as the gold standard) and in vivo stress experiments. Compared with OF strain imaging, the proposed estimator improved the accuracy of principal major and minor strains during stress tests, with an average contrast-to-noise ratio improvement of 4.4 ± 2.7 dB ( p -value < 0.01). Moreover, there was a significant correlation and a very close agreement between the proposed estimator and sonomicrometry for tested heart rates between 60 and 180 beats per minute (bpm). The averages ± standard deviations (STD) of R² and biases ± STD between them were 0.96 ± 0.04 ( p -value < 0.01) and 0.01 ± 0.03% in the axial direction, respectively; and 0.94 ± 0.02 ( p -value < 0.01) and 0.04 ± 0.06% in the lateral direction, respectively. These results suggest that the proposed estimator could be useful clinically to provide an accurate and quantitative 2-D large lateral field-of-view myocardial strain assessment at high heart rates during stress echocardiography.

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.

Virtual Source Synthetic Aperture for Accurate Lateral Displacement Estimation in Ultrasound Elastography

Ultrasound elastography (USE) is an emerging noninvasive imaging technique in which pathological alterations can be visualized by revealing the mechanical properties of the tissue. Estimating tissue displacement in all directions is required to accurately estimate the mechanical properties. Despite capabilities of elastography techniques in estimating displacement in both axial and lateral directions, estimation of axial displacement is more accurate than lateral direction due to higher sampling frequency, higher resolution, and having a carrier signal propagating in the axial direction. Among different ultrasound imaging techniques, synthetic aperture (SA) has better lateral resolution than others, but it is not commonly used for USE due to its limitation in imaging depth of field. Virtual source synthetic aperture (VSSA) imaging is a technique to implement SA beamforming on the focused transmitted data to overcome the limitation of SA in depth of field while maintaining the same lateral resolution as SA. Besides lateral resolution, VSSA has the capability of increasing sampling frequency in the lateral direction without interpolation. In this article, we utilize VSSA to perform beamforming to enable higher resolution and sampling frequency in the lateral direction. The beamformed data are then processed using our recently published elastography technique, OVERWIND. Simulation and experimental results show substantial improvement in the estimation of lateral displacements.

Parameterized Strain Estimation for Vascular Ultrasound Elastography With Sparse Representation

Ultrasound vascular strain imaging has shown its potential to interrogate the motion of the vessel wall induced by the cardiac pulsation for predicting plaque instability. In this study, a sparse model strain estimator (SMSE) is proposed to reconstruct a dense strain field at a high resolution, with no spatial derivatives, and a high computation efficiency. This sparse model utilizes the highly-compacted property of discrete cosine transform (DCT) coefficients, thereby allowing to parameterize displacement and strain fields with truncated DCT coefficients. The derivation of affine strain components (axial and lateral strains and shears) was reformulated into solving truncated DCT coefficients and then reconstructed with them. Moreover, an analytical solution was derived to reduce estimation time. With simulations, the SMSE reduced estimation errors by up to 50% compared with the state-of-the-art window-based Lagrangian speckle model estimator (LSME). The SMSE was also proven to be more robust than the LSME against global and local noise. For in vitro and in vivo tests, residual strains assessing cumulated errors with the SMSE were 2 to 3 times lower than with the LSME. Regarding computation efficiency, the processing time of the SMSE was reduced by 4 to 25 times compared with the LSME, according to simulations, in vitro and in vivo results. Finally, phantom studies demonstrated the enhanced spatial resolution of the proposed SMSE algorithm against LSME.

Deep Learning for Carotid Plaque Segmentation using a Dilated U-Net Architecture

Carotid plaque segmentation in ultrasound longitudinal B-mode images using deep learning is presented in this work. We report on 101 severely stenotic carotid plaque patients. A standard U-Net is compared with a dilated U-Net architecture in which the dilated convolution layers were used in the bottleneck. Both a fully automatic and a semi-automatic approach with a bounding box was implemented. The performance degradation in plaque segmentation due to errors in the bounding box is quantified. We found that the bounding box significantly improved the performance of the networks with U-Net Dice coefficients of 0.48 for automatic and 0.83 for semi-automatic segmentation of plaque. Similar results were also obtained for the dilated U-Net with Dice coefficients of 0.55 for automatic and 0.84 for semi-automatic when compared to manual segmentations of the same plaque by an experienced sonographer. A 5% error in the bounding box in both dimensions reduced the Dice coefficient to 0.79 and 0.80 for U-Net and dilated U-Net respectively.

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.

Lagrangian carotid strain imaging indices normalized to blood pressure for vulnerable plaque

OBJECTIVE

Ultrasound Lagrangian carotid strain imaging (LCSI) utilizes physiological deformation caused by arterial pressure variations to generate strain tensor maps of the vessel walls and plaques. LCSI has been criticized for the lack of normalization of magnitude-based strain indices to physiological stimuli, namely blood pressure. We evaluated the impact of normalization of magnitude-based strain indices to blood pressure measured immediately after the acquisition of radiofrequency (RF) data loops for LCSI.

MATERIALS AND METHODS

A complete clinical ultrasound examination along with RF data loops for LCSI was performed on 50 patients (30 males and 20 females) who presented with >60% carotid stenosis and were scheduled for carotid endarterectomy. Cognition was assessed using the 60-minute neuropsychological test protocol.

RESULTS

For axial strains correlation of maximum accumulated strain indices (MASI), cognition scores were -0.46 for non-normalized and -0.45, -0.49, -0.37, and -0.48 for systolic, diastolic, pulse pressure, and mean arterial pressure normalized data, respectively. The corresponding area under the curve (AUC) values for classifiers designed using maximum likelihood estimation of a binormal distribution with a median-split of the executive function cognition scores were 0.73, 0.70, 0.71, 0.70, and 0.71, respectively.

CONCLUSIONS

No significant differences in the AUC estimates were obtained between normalized and non-normalized magnitude-based strain indices.

A novel method for describing biomechanical properties of the aortic wall based on the three-dimensional fluid-structure interaction model

OBJECTIVES

Our goal was to present a novel non-invasive approach for assessment of aortic wall displacement to describe its biomechanical properties during the cardiac cycle.

METHODS

The fluid-structure interaction (FSI) technique was used to reconstruct aortic wall displacement based on computed tomography angiography and 2-dimensional speckle-tracking technique (2DSTT) data collected from 20 patients [10 with healthy aortas (AA) and 10 with abdominal aortic aneurysms (AAAs)]. The mechanical properties of the wall of the aorta were described by the Yeoh hyperelastic materials model with α and β parameters, and wall displacement was determined with 2DSTT. The mechanical parameters of the wall of the aorta in the FSI model were automatically updated in the calculation loop until the calculated and clinically measured wall movements were the same.

RESULTS

Results showed 98% accuracy of FSI compared to 2DSTT for AA and AAA (P > 0.05). The mean wall deformation for AA was 2.45 ± 0.12 mm and 2.49 ± 0.10 mm for FSI and 2DSTT, respectively (P = 0.40), whereas that for AAA was 2.84 ± 0.44 mm and 2.88 ± 0.45 mm, respectively (P = 0.83). The FSI analysis indicated that the α and β parameters for AA were equal to 14.35 ± 1.30 N⋅cm-2 and 9.33 ± 1.08 N⋅cm-2, respectively; and for AAA, α was 11.00 ± 0.49 N⋅cm-2 and β was 79.46 ± 4.32 N⋅cm-2.

CONCLUSIONS

The FSI technique may be successfully applied to assess the mechanical parameters of patient-specific aortic walls using computed tomography angiographic and 2DSTT measurements.

3-D Strain Imaging of the Carotid Bifurcation: Methods and in-Human Feasibility

Atherosclerotic plaque development in the carotid artery bifurcation elevates the risk for stroke, which is often initiated by plaque rupture. The risk-to-rupture of a plaque is related to its composition. Two-dimensional non-invasive carotid elastography studies have found a correlation between wall strain and plaque composition. This study introduces a technique to perform non-invasive volumetric elastography in vivo. Three-dimensional ultrasound data of carotid artery bifurcations were acquired in four asymptomatic individuals using an electrocardiogram-triggered multislice acquisition device that scanned over a length of 35 mm (350 slices) using a linear transducer (L11-3, f_c = 9 MHz). For each slice, three-angle ultrasound plane wave data were acquired and beamformed. A correction for breathing-induced motion was applied to spatially align the slices, enabling 3-D cross-correlation-based compound displacement, distensibility and strain estimation. Distensibility values matched with previously published values, while the corresponding volumetric principal strain maps revealed locally elevated compressive and tensile strains. This study presents for the first time 3-D elastography of carotid arteries in vivo.

In vivo carotid strain imaging using principal strains in longitudinal view

Carotid plaque rupture can result in stroke or transient ischemic attack that can be devastating for patients. Ultrasound strain imaging provides a noninvasive method to identify unstable plaque likely to rupture. Axial, lateral and shear strains in carotid plaque have been shown to be linked to carotid plaque instability. Recently, there has been interest in using principal strains, which do not depend on angle of insonification of the carotid artery for quantifying instability in plaque along the longitudinal view. In this work relationships between angle dependent axial, lateral and shear strain along with axis independent principal strains are compared. Three strain indices were defined, 1) Average Mean Strain (AMS), 2) Maximum Mean Strain (MMS) and 3) Mean Standard Deviation (MSD) to identify relationships between these five strain image types in a group of 76 in vivo patients. The maximum principal strain demonstrated the highest strain values when compared to axial strain for all patients with a linear regression slope of 1.6 and a y intercept of 2.4 percent strain for AMS. The maximum shear strain when compared to shear strain had a slope of 1.15 and a y intercept of 0.21 percent for AMS. Next, the effect of insonification angle, which is the angle subtended by the artery at the location of plaque was studied. Patients were divided into three sub groups, i.e. less than 5 degrees (n = 31), between 5 and 10 degrees (n = 24) and above 10 degrees (n = 21). The angle of insonification did not make a significant difference between the three angle groups when comparing the relationship between the angle dependent and independent strain values.

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.

Spatial Angular Compounding With Affine-Model-Based Optical Flow for Improvement of Motion Estimation

Tissue motion estimation is an essential step for ultrasound elastography. Our previous study has shown that the affine-model-based optical flow (OF) method outperforms the normalized cross-correlation-based block matching (BM) method in motion estimation. However, the quality of lateral estimation using OF is still low due to inherent limitation of ultrasound imaging. BM-based spatial angular compounding (SAC) has been developed to obtain better motion estimation. In this paper, OF-based SAC (OF-SAC) is proposed to further improve the performance of lateral (and axial) estimation, and it is compared with BM-based SAC (BM-SAC). Plane wave as well as focused wave is transmitted in both simulations and phantom experiments on a linear array. In order to compare the performance quantitatively, the root-mean-square error (RMSE) of axial/lateral displacement and strain, and signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) of axial/lateral strain are used as the evaluation criteria in the simulations. In the phantom experiments, the SNR and CNR are used to assess the quality of axial/lateral strain. The results show that for both OF and BM, SAC improves the performance of motion estimation, regardless of using plane or focused wave transmission. More importantly, OF-SAC is shown to outperform BM-SAC with lower RMSE, higher SNR, and higher CNR. In addition, preliminary in vivo experiments on the carotid artery of a healthy human subject also prove the superiority of OF-SAC. These results suggest that OF-SAC is preferred for both axial and lateral motion estimation to BM-SAC.

Intraluminal Ultrasonic Palpation Imaging Technique Revisited for Anisotropic Characterization of Healthy and Atherosclerotic Coronary Arteries: A Feasibility Study

Accurate mechanical characterization of coronary atherosclerotic lesions remains essential for the in vivo detection of vulnerable plaques. Using intravascular ultrasound strain measurements and based on the mechanical response of a circular and concentric vascular model, E. I. Céspedes, C. L. de Korte and A. F. van der Steen developed an elasticity-palpography technique in 2000 to estimate the apparent stress-strain modulus palpogram of the thick subendoluminal arterial wall layer. More recently, this approach was improved by our group to consider the real anatomic shape of the vulnerable plaque. Even though these two studies highlighted original and promising approaches for improving the detection of vulnerable plaques, they did not overcome a main limitation related to the anisotropic mechanical behavior of the vascular tissue. The present study was therefore designed to extend these previous approaches by considering the orthotropic mechanical properties of the arterial wall and lesion constituents. Based on the continuum mechanics theory prescribing the strain field, an elastic anisotropy index was defined. This new anisotropic elasticity-palpography technique was successfully applied to characterize ten coronary plaque and one healthy vessel geometries of patients imaged in vivo with intravascular ultrasound. The results revealed that the anisotropy index-palpograms were estimated with a good accuracy (with a mean relative error of 26.8 ± 48.8%) compared with ground true solutions.

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.

Simultaneous Vascular Strain and Blood Vector Velocity Imaging Using High-Frequency Versus Conventional-Frequency Plane Wave Ultrasound: A Phantom Study

Plaque strain and blood vector velocity imaging of stenosed arteries are expected to aid in diagnosis and prevention of cerebrovascular disease. Ultrafast plane wave imaging enables simultaneous strain and velocity estimation. Multiple ultrasound vendors are introducing high-frequency ultrasound probes and systems. This paper investigates whether the use of high-frequency ultrafast ultrasound is beneficial for assessing blood velocities and strain in arteries. The performance of strain and blood flow velocity estimation was compared between a high-frequency transducer (MS250, fc = 21 MHz) and a clinically utilized transducer (L12-5, fc = 9 MHz). Quantitative analysis based on straight tube phantom experiments revealed that the MS250 outperformed the L12-5 in the superficial region: low velocities near the wall were more accurately estimated and wall strains were better resolved. At greater than 2-cm echo depth, the L12-5 performed better due to the high attenuation of the MS250 probe. Qualitative comparison using a perfused patient-specific carotid bifurcation phantom confirmed these findings. Thus, in conclusion, for strain and blood velocity estimation for depths up to ~2 cm, a high-frequency probe is recommended.

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).

Model-based vascular elastography improves the detection of flow-induced carotid artery remodeling in mice

Increased arterial thickness measured with ultrasound correlates with future cardiovascular events, but conventional ultrasound imaging techniques cannot distinguish between intima, media, or atherosclerotic plaque in the carotid artery. In this work, we evaluated how well vascular elastography can detect intimal changes in a mouse model of carotid remodeling. We ligated the left external and internal branches of the carotid artery of male FVB mice and performed sham operations for 2 weeks. High-resolution ultrasound imaging accurately detected lower blood velocities and low blood volume flow in the carotid arteries after ligation in FVB mice. However, ultrasound could not detect differences in the carotid wall even at 2 weeks post-surgery. The Young’s modulus was measured based on displacements of the carotid artery wall, and Young’s modulus was 2-fold greater in shams at 1 week post ligation, and 3-fold greater 2 weeks after ligation. Finally, the higher Young’s modulus was most associated with higher intimal thickness but not medial or adventitial thickness as measured by histology. In conclusion, we developed a robust ultrasound-based elastography method for early detection of intimal changes in small animals.

Revisiting the Cramér Rao Lower Bound for Elastography: Predicting the Performance of Axial, Lateral and Polar Strain Elastograms

We derived the Cramér Rao lower bound for 2-D estimators employed in quasi-static elastography. To illustrate the theory, we modeled the 2-D point spread function as a sinc-modulated sine pulse in the axial direction and as a sinc function in the lateral direction. We compared theoretical predictions of the variance incurred in displacements and strains when quasi-static elastography was performed under varying conditions (different scanning methods, different configuration of conventional linear array imaging and different-size kernels) with those measured from simulated or experimentally acquired data. We performed studies to illustrate the application of the derived expressions when performing vascular elastography with plane wave and compounded plane wave imaging. Standard deviations in lateral displacements were an order higher than those in axial. Additionally, the derived expressions predicted that peak performance should occur when 2% strain is applied, the same order of magnitude as observed in simulations (1%) and experiments (1%-2%). We assessed how different configurations of conventional linear array imaging (number of active reception and transmission elements) influenced the quality of axial and lateral strain elastograms. The theoretical expressions predicted that 2-D echo tracking should be performed with wide kernels, but the length of the kernels should be selected using knowledge of the magnitude of the applied strain: specifically, longer kernels for small strains (<5%) and shorter kernels for larger strains. Although the general trends of theoretical predictions and experimental observations were similar, biases incurred during beamforming and subsample displacement estimation produced noticeable differences.