Noninvasive carotid strain imaging using angular compounding at large beam steered angles: validation in vessel phantoms

Layer-specific radiofrequency ultrasound-based strain analysis in a porcine model of ischemic cardiomyopathy validated by a geometric model

Local layer-specific myocardial deformation after myocardial infarction (MI) has not been studied extensively although the sub-endocardium is more vulnerable to ischemia and interstitial fibrosis deposition. Radiofrequency (RF) ultrasound-based analysis could provide superior layer-specific radial strain estimation compared with clinically available deformation imaging techniques. In this study, we used RF-based myocardial deformation measurements to investigate layer-specific differences between healthy and damaged myocardium in a porcine model of chronic MI. RF data were acquired epicardially in healthy (n = 21) and infarcted (n = 5) regions of a porcine chronic MI model 12 wk post-MI. Radial and longitudinal strains were estimated in the sub-endocardial, mid-wall and sub-epicardial layers of the left ventricle. Collagen content was quantified in three layers of healthy and infarcted regions in five pigs. An analytical geometric model of the left ventricle was used to theoretically underpin the radial deformation estimated in different myocardial layers. Means ± standard errors of the peak radial and longitudinal strain estimates of the sub-endocardial, mid-wall and sub-epicardial layers of the healthy and infarcted tissue were: 82.7 ± 5.2% versus 39.9 ± 10.8% (p = 0.002), 63.6 ± 3.3% versus 38.8 ± 7.7% (p = 0.004) and 34.3 ± 3.0% versus 35.1 ± 5.2% (p = 0.9), respectively. The radial strain gradient between the sub-endocardium and the sub-epicardium had decreased 12 wk after MI, and histologic examination revealed the greatest increases in collagen in the sub-endocardial and mid-wall layers. Comparable normal peak radial strain values were found by geometric modeling when input values were derived from the in vivo measurements and literature. In conclusion, the estimated strain values are realistic and indicate that sub-endocardial radial strain in healthy tissue can amount to 80%. This high value can be explained by the cardiac geometry, as was illustrated by geometric modeling. After MI, strain values were decreased and collagen content was increased in the sub-endocardial and mid-wall layers. Layer-specific peak radial strain can be assessed by RF strain estimation and clearly differs between healthy and infarcted tissue. Although the relationship between tissue stiffness and tissue strain is not strictly local, this novel technique provides a valuable way to assess layer-specific regional cardiac function in a variety of myocardial diseases.

MRI strain imaging of the carotid artery: present limitations and future challenges

Rupture of atherosclerotic plaques in the carotid artery is a main cause of stroke. Current diagnostics are not sufficient to identify all rupture-prone plaques, and studies have shown that biomechanical factors improve current plaque risk assessment. Strain imaging may be a valuable contribution to this risk assessment. MRI is a versatile imaging technique that offers various methods that are capable of measuring tissue strain. In this review, MR imaging techniques with displacement (DENSE), velocity (PC MRI), or strain (SENC) encoding protocols are discussed, together with post-processing techniques based on time-resolved MRI data. Although several MRI techniques are being developed to improve time-resolved MR imaging, current technical limitations related to spatial and temporal resolutions render MRI strain imaging currently unfit for carotid plaque strain evaluation. A novel approach using non-rigid image registration of MR images to determine strain in carotid arteries based on black blood cine MRI is proposed in this review. This and other post-processing techniques based on time-resolved MRI data may provide a good estimate of plaque strain, but are also dependent on the spatial and temporal resolution of the MR images. However, they seem to be the most promising approach for MRI based plaque strain analysis in the near future.

Correlation of cognitive function with ultrasound strain indices in carotid plaque

Instability in carotid vulnerable plaque can generate cerebral micro-emboli, which may be related to both stroke and eventual cognitive abnormality. Strain imaging to detect plaque vulnerability based on regions with large strain fluctuations, with arterial pulsation, may be able to determine the risk of cognitive impairment. Plaque instability may be characterized by increased strain variations over a cardiac cycle. Radiofrequency signals for ultrasound strain imaging were acquired from the carotid arteries of 24 human patients using a Siemens Antares with a VFX 13-5 linear array transducer. These patients underwent standardized cognitive assessment (Repeatable Battery for the Assessment of Neuropsychological Status [RBANS]). Plaque regions were segmented by a radiologist at end-diastole using the Medical Imaging Interaction Toolkit. A hierarchical block-matching motion tracking algorithm was used to estimate the cumulated axial, lateral and shear strains within the imaging plane. The maximum, minimum and peak-to-peak strain indices in the plaque computed from the mean cumulated strain over a small region of interest in the plaque with large deformations were obtained. The maximum and peak-to-peak mean cumulated strain indices over the entire plaque region were also computed. All strain indices were then correlated with RBANS Total performance. Overall cognitive performance (RBANS Total) was negatively associated with values of the maximum strain and the peak-to-peak for axial and lateral strains, respectively. There was no significant correlation between the RBANS Total score and shear strain and strain indices averaged over the entire identified plaque for this group of patients. However, correlation of maximum lateral strain was higher for symptomatic patients (r = -0.650, p = 0.006) than for asymptomatic patients (r = -0.115, p = 0.803). On the other hand, correlation of maximum axial strain averaged over the entire plaque region was significantly higher for asymptomatic patients (r = -0.817, p = 0.016) than for symptomatic patients (r = -0.224, p = 0.402). The results reveal a direct relationship between the maximum axial and lateral strain indices in carotid plaque and cognitive impairment.

Ultrasound elastography using multiple images

Displacement estimation is an essential step for ultrasound elastography and numerous techniques have been proposed to improve its quality using two frames of ultrasound RF data. This paper introduces a technique for calculating a displacement field from three (or multiple) frames of ultrasound RF data. To calculate a displacement field using three images, we first derive constraints on variations of the displacement field with time using mechanics of materials. These constraints are then used to generate a regularized cost function that incorporates amplitude similarity of three ultrasound images and displacement continuity. We optimize the cost function in an expectation maximization (EM) framework. Iteratively reweighted least squares (IRLS) is used to minimize the effect of outliers. An alternative approach for utilizing multiple images is to only consider two frames at any time and sequentially calculate the strains, which are then accumulated. We formally show that, compared to using two images or accumulating strains, the new algorithm reduces the noise and eliminates ambiguities in displacement estimation. The displacement field is used to generate strain images for quasi-static elastography. Simulation, phantom experiments and in vivo patient trials of imaging liver tumors and monitoring ablation therapy of liver cancer are presented for validation. We show that even with the challenging patient data, where it is likely to have one frame among the three that is not optimal for strain estimation, the introduction of physics-based prior as well as the simultaneous consideration of three images significantly improves the quality of strain images. Average values for strain images of two frames versus ElastMI are: 43 versus 73 for SNR (signal to noise ratio) in simulation data, 11 versus 15 for CNR (contrast to noise ratio) in phantom data, and 5.7 versus 7.3 for CNR in patient data. In addition, the improvement of ElastMI over both utilizing two images and accumulating strains is statistically significant in the patient data, with p-values of respectively 0.006 and 0.012.

A novel experimental approach for three-dimensional geometry assessment of calcified human stenotic arteries in vitro

To improve diagnosis and understanding of the risk of rupture of atherosclerotic plaque, new strategies to realistically determine mechanical properties of atherosclerotic plaque need to be developed. In this study, an in vitro experimental method is proposed for accurate 3-D assessment of (diseased) vessel geometry using ultrasound. The method was applied to a vascular phantom, a healthy porcine carotid artery and human carotid endarterectomy specimens (n = 6). Vessel segments were pressure fixed and rotated in 10 ° steps. Longitudinal cross sections were imaged over 360 °. Findings were validated using micro-computed tomography (μCT). Results show good agreement between ultrasound and μCT-based geometries of the different segment types (ISI phantom = 0.94, ISI healthy = 0.79, ISI diseased = 0.75-0.80). The method does not suffer from acoustic shadowing effects present when imaging stenotic segments and allows future dynamic measurements to determine mechanical properties of atherosclerotic plaque in an in vitro setting.

Motion Artifact Reduction in Ultrasound Based Thermal Strain Imaging of Atherosclerotic Plaques Using Time Series Analysis

Large lipid pools in vulnerable plaques, in principle, can be detected using US based thermal strain imaging (US-TSI). One practical challenge for in vivo cardiovascular application of US-TSI is that the thermal strain is masked by the mechanical strain caused by cardiac pulsation. ECG gating is a widely adopted method for cardiac motion compensation, but it is often susceptible to electrical and physiological noise. In this paper, we present an alternative time series analysis approach to separate thermal strain from the mechanical strain without using ECG. The performance and feasibility of the time-series analysis technique was tested via numerical simulation as well as in vitro water tank experiments using a vessel mimicking phantom and an excised human atherosclerotic artery where the cardiac pulsation is simulated by a pulsatile pump.

Suppression of grating lobe artifacts in ultrasound images formed from diverging transmitting beams by modulation of receiving beams

In linear array transducers, owing to regular spacing of the array elements, grating lobes exist in transmission and reception. In ultrasonic imaging involving the use of diverging (unfocused) transmitting beams and steered receiving beams by linear transducer arrays, aperture apodization and spatial combination of steered receiving beams from multiple transmissions are not sufficient to suppress receive-grating lobe artifacts. To further suppress receive-grating lobe artifacts in reconstructed B-mode images, we propose a technique of modulating the receiving beams by a factor that is governed by the envelope of a corresponding signal, which is formed by filtering the receiving beam with a zero-phase low-pass filter with a cut-off frequency that is determined by the receiving beam steering angle. This technique suppressed receive-grating lobe artifacts without significant loss in spatial resolution in offline reconstructed B-mode images from simulation, phantom and in vivo imaging of the carotid artery. In a simulation of point scatterers, a relative reduction in grating lobe artifacts of 40 dB was realized in images from diverging beam scanning.

Noninvasive vascular displacement estimation for relative elastic modulus reconstruction in transversal imaging planes

Atherosclerotic plaque rupture can initiate stroke or myocardial infarction. Lipid-rich plaques with thin fibrous caps have a higher risk to rupture than fibrotic plaques. Elastic moduli differ for lipid-rich and fibrous tissue and can be reconstructed using tissue displacements estimated from intravascular ultrasound radiofrequency (RF) data acquisitions. This study investigated if modulus reconstruction is possible for noninvasive RF acquisitions of vessels in transverse imaging planes using an iterative 2D cross-correlation based displacement estimation algorithm. Furthermore, since it is known that displacements can be improved by compounding of displacements estimated at various beam steering angles, we compared the performance of the modulus reconstruction with and without compounding. For the comparison, simulated and experimental RF data were generated of various vessel-mimicking phantoms. Reconstruction errors were less than 10%, which seems adequate for distinguishing lipid-rich from fibrous tissue. Compounding outperformed single-angle reconstruction: the interquartile range of the reconstructed moduli for the various homogeneous phantom layers was approximately two times smaller. Additionally, the estimated lateral displacements were a factor of 2–3 better matched to the displacements corresponding to the reconstructed modulus distribution. Thus, noninvasive elastic modulus reconstruction is possible for transverse vessel cross sections using this cross-correlation method and is more accurate with compounding.

Noninvasive vascular elastography using plane-wave and sparse-array imaging

Stroke may occur when an atherosclerotic plaque ruptures in the carotid artery. Noninvasive vascular elastography (NIVE) visualizes the strain distribution within the carotid artery, which is related to its mechanical properties that govern plaque rupture. Strain elastograms obtained from the transverse plane of the carotid artery are difficult to interpret, because strain is estimated in Cartesian coordinates. Sparsearray (SA) elastography overcomes this problem by transforming shear and normal strain to polar coordinates. However, the SA's transmit power may be too weak to produce useful elastograms in the clinical setting. Consequently, we are exploring other imaging methods to solve this potential problem. This study evaluated the quality of elastograms produced with SA imaging, plane-wave (PW) imaging, and compounded-plane-wave (CPW) imaging. We performed studies on simulated and physical vessel phantoms, and the carotid artery of a healthy volunteer. All echo imaging was performed with a linear transducer array that contained 128 elements, operating at 5 MHz. In SA imaging, 7 elements were fired during transmission, but all 128 elements were active during reception. In PW imaging, all 128 elements were active during both transmission and reception. We created CPW images by steering the acoustic beam within the range of -15° to 15° in increments of 5°. SA radial and circumferential strain elastograms were comparable to those produced using PW and CPW imaging. Additionally, side-lobe levels incurred during SA imaging were 20 dB lower than those produced during PW imaging, and 10 dB lower than those computed using CPW imaging. Overall, SA imaging performs well in vivo; therefore, we plan to improve the technique and perform preclinical studies.

Noninvasive estimation of the blood pressure waveform in the carotid artery using continuous finger blood pressure monitoring

To noninvasively estimate the blood pressure continuously in the common carotid artery (CCA), we obtained the distension waveform of the CCA from seven healthy volunteers and 20 hypertensive patients using radio-frequency ultrasound. Consequently, it was calibrated by the mean and diastolic pressure measured in the finger artery and compared with applanation tonometry, calibrated using the systolic and diastolic pressure in the brachial artery. The mean difference in estimating the mean blood pressure was 0.3 mm Hg (limits of agreement: -11.7 to 12.3 mm Hg). In estimating the systolic blood pressure, the mean difference was 8.0 mm Hg (limits of agreement: -29.8 to 45.8 mm Hg) and showed increasing variation with blood pressure. The systolic blood pressure values can be expected between 0.83 and 1.35 times the control method. In this study, we obtained proof-of-principle for noninvasively measuring blood pressure in the CCA using continuous finger blood pressure monitoring. This opens the way to estimating location specific arterial stiffness and intra-plaque elasticity.

Methods for robust in vivo strain estimation in the carotid artery

A hierarchical block-matching motion tracking algorithm for strain imaging is presented. Displacements are estimated with improved robustness and precision by utilizing a Bayesian regularization algorithm and an unbiased subsample interpolation technique. A modified least-squares strain estimator is proposed to estimate strain images from a noisy displacement input while addressing the motion discontinuity at the wall-lumen boundary. Methods to track deformation over the cardiac cycle incorporate a dynamic frame skip criterion to process data frames with sufficient deformation to produce high signal-to-noise displacement and strain images. Algorithms to accumulate displacement and/or strain on particles in a region of interest over the cardiac cycle are described. New methods to visualize and characterize the deformation measured with the full 2D strain tensor are presented. Initial results from patients imaged prior to carotid endarterectomy suggest that strain imaging detects conditions that are traditionally considered high risk including soft plaque composition, unstable morphology, abnormal hemodynamics and shear of plaque against tethering tissue can be exacerbated by neoangiogenesis. For example, a maximum absolute principal strain exceeding 0.2 is observed near calcified regions adjacent to turbulent flow, protrusion of the plaque into the arterial lumen and regions of low echogenicity associated with soft plaques. Non-invasive carotid strain imaging is therefore a potentially useful tool for detecting unstable carotid plaque.

Visualizing the radial and circumferential strain distribution within vessel phantoms using synthetic-aperture ultrasound elastography

Noninvasive elastography (NIVE) produces elastograms that are difficult to interpret because NIVE visualizes strain in the transducer coordinate system. In this paper, we hypothesized that transforming normal and shear strain elastograms to the vessel coordinate system will produce better strain elastograms. To corroborate this hypothesis, we acquired synthetic-aperture (SA) ultrasound data from simulated and physical vessel phantoms. In both studies, SA echo frames were reconstructed from data acquired with a sparse transducer array. The simulation study was performed with homogeneous and heterogenous phantoms, but in the experimental study we used a modified ultrasound scanner to acquire SA data from homogeneous (n = 1) and heterogeneous (n = 3) vessel phantoms. Axial and lateral displacements were estimated by performing two-dimensional cross-correlation analysis on the beamformed RF echo frames. We generated radial and circumferential strain elastograms by transforming normal and shear strain elastograms to the vessel coordinate system. The results revealed: 1) radial and circumferential strain elastograms acquired from simulated data had a relative root mean squared error on the order of 0.1%; 2) experimentally acquired radial and circumferential strain elastograms had elastographic contrast-to-noise ratio (CNRe) between 10 and 40 dB, and elastographic signal-to-noise ratio (SNRe) between 10 and 35 dB, depending on the number of active transmission elements employed during imaging; 3) radial and circumferential strain elastograms produced with fewer than 8 active transmission elements were inferior to those computed with a greater number of active elements; and 4) plaques were evident in the strain elastograms, except in those obtained with the sparsest transducer array. This study demonstrated that a syntheticaperture ultrasound system could visualize radial and circumferential strain noninvasively.

Filter-based compounded delay estimation with application to strain imaging

Ultrasonic wave interference produces local fluctuations in both the envelope, known as speckle, and phase of echoes. Furthermore, such fluctuations are correlated in space, and subsequent motion estimation from the envelope and/or phase signal produces patterned, correlated errors. Compounding, or combining information from multiple decorrelated looks, reduces such effects. We propose using a filter bank to create multiple looks to produce a compounded motion estimate. In particular, filtering in the lateral direction is shown to preserve delay estimation accuracy in the filtered sub-bands while creating decorrelation between sub-bands at the expense of some lateral resolution. For Gaussian apodization, we explicitly compute the induced signal decorrelation produced by Gabor filters. Furthermore, it is shown that lateral filtering is approximately equivalent to steering, in which filtered sub-bands correspond to signals extracted from shifted sub-apertures. Field II simulation of a point spread function verifies this claim. We use phase zero and its variants as displacement estimators for our compounded result. A simplified deformation model is used to provide computer simulations of deforming an elastic phantom. Simulations demonstrate root mean square error (RMSE) reduction in both displacement and strain of the compounded result over conventional and its laterally blurred versions. Then we apply the methods to experimental data using a commercial elastic phantom, demonstrating an improvement in strain SNR.

Ultrasound-based radial and longitudinal strain estimation of the carotid artery: a feasibility study

Ultrasound-based estimation of arterial wall elasticity is commonly used to assess arterial stiffness. However, previous elastography studies have mostly addressed radial strain measurements, and the longitudinal strain has been more or less ignored. This study shows the feasibility of a speckle-tracking-based algorithm for simultaneous estimation of radial and longitudinal strain in the carotid artery in silico. Additionally, these results were preliminarily confirmed in vivo.

Vascular ultrasound for atherosclerosis imaging

Cardiovascular disease is a leading cause of death in the Western world. Therefore, detection and quantification of atherosclerotic disease is of paramount importance to monitor treatment and possible prevention of acute events. Vascular ultrasound is an excellent technique to assess the geometry of vessel walls and plaques. The high temporal as well as spatial resolution allows quantification of luminal area and plaque size and volume. While carotid arteries can be imaged non-invasively, scanning of coronary arteries requires invasive intravascular catheters. Both techniques have already demonstrated their clinical applicability. Using linear array technology, detection of disease as well as monitoring of pharmaceutical treatment in carotid arteries are feasible. Data acquired with intravascular ultrasound catheters have proved to be especially beneficial in understanding the development of atherosclerotic disease in coronary arteries. With the introduction of vascular elastography not only the geometry of plaques but also the risk for rupture of plaques might be identified. These so-called vulnerable plaques are frequently not flow-limiting and rupture of these plaques is responsible for the majority of cerebral and cardiac ischaemic events. Intravascular ultrasound elastography studies have demonstrated a high correlation between high strain and vulnerable plaque features, both ex vivo and in vivo. Additionally, pharmaceutical intervention could be monitored using this technique. Non-invasive vascular elastography has recently been developed for carotid applications by using compound scanning. Validation and initial clinical evaluation is currently being performed. Since abundance of vasa vasorum (VV) is correlated with vulnerable plaque development, quantification of VV might be a unique tool to even prevent this from happening. Using ultrasound contrast agents, it has been demonstrated that VV can be identified and quantified. Although far from routine clinical application, non-invasive and intravascular ultrasound VV imaging might pave the road to prevent atherosclerotic disease in an early phase. This paper reviews the conventional vascular ultrasound techniques as well as vascular ultrasound strain and vascular ultrasound VV imaging.

An angular compounding technique using displacement projection for noninvasive ultrasound strain imaging of vessel cross-sections

Strain is considered to be a useful indicator of atherosclerotic plaque vulnerability. This study introduces an alternative for a recently introduced strain imaging method that combined beam steered ultrasound acquisitions to construct radial strain images of transverse cross-sections of superficial arteries. In that study, axial strains were projected in the radial direction. Using the alternative method introduced in this study, axial displacements are projected radially, followed by a least squares estimation of radial strains. This enables the use of a larger projection angle. Consequently, fewer acquisitions at smaller beam steering angles are required to construct radial strain images. Simulated and experimentally obtained radio-frequency data of radially expanding vessel phantoms were used to compare the two methods. Using only three beam steering angles (-30°, 0° and 30°), the new method outperformed the older method that used seven different angles and up to 45° of beam steering: the root mean squared error was reduced by 38% and the elastographic signal- and contrast-to-noise ratios increased by 1.8 dB and 4.9 dB, respectively. The new method was also superior for homogeneous and heterogeneous phantoms with eccentric lumens. To conclude, an improved noninvasive method was developed for radial strain imaging in transverse cross-sections of superficial arteries.