2-D Versus 3-D Cross-Correlation-Based Radial and Circumferential Strain Estimation Using Multiplane 2-D Ultrafast Ultrasound in a 3-D Atherosclerotic Carotid Artery Model

Longitudinal intravital imaging of mouse placenta

Studying placental functions is crucial for understanding pregnancy complications. However, imaging placenta is challenging due to its depth, volume, and motion distortions. In this study, we have developed an implantable placenta window in mice that enables high-resolution photoacoustic and fluorescence imaging of placental development throughout the pregnancy. The placenta window exhibits excellent transparency for light and sound. By combining the placenta window with ultrafast functional photoacoustic microscopy, we were able to investigate the placental development during the entire mouse pregnancy, providing unprecedented spatiotemporal details. Consequently, we examined the acute responses of the placenta to alcohol consumption and cardiac arrest, as well as chronic abnormalities in an inflammation model. We have also observed viral gene delivery at the single-cell level and chemical diffusion through the placenta by using fluorescence imaging. Our results demonstrate that intravital imaging through the placenta window can be a powerful tool for studying placenta functions and understanding the placental origins of adverse pregnancy outcomes.

Influence of key parameters on motion artifacts in lateral strain estimation with spatial angular compounding

Strain imaging can reveal the changes in tissue mechanical properties related to pathological alterations by estimating tissue strains in the lateral and axial directions of ultrasound imaging. The estimation performance in the lateral direction is usually worse than that in the axial direction. Spatial angular compounding (SAC) has been demonstrated to improve the quality of lateral estimation by deriving the lateral displacements using axial displacements obtained from multi-angle transmissions. However, motion and deformation of tissues during multiple transmissions may cause motion artifacts, and thus deteriorate the quality of strain estimation. These artifacts can be reduced by choosing appropriate imaging parameters. However, few studies have been conducted to evaluate the influences of key parameters in strain estimation, such as the pulse repetition frequency (PRF), the number of steering angles (NSA), and the maximum steering angles (MSA), in terms of performance optimization. Therefore, this study aims to investigate the effects of these parameters through simulations and phantom experiments. The performance of strain estimation is evaluated by measuring the root-mean-square error (RMSE) and the standard deviation (SD) in the simulations and phantom experiments, respectively. The contrast-to-noise ratio (CNR) of strain images is calculated in both the simulations and phantom experiments. The results show that motion artifacts in strain estimation can be reduced by increasing the PRF to 1 kHz. When the PRF reaches 1 kHz, further increase of the PRF shows little obvious improvement in strain estimation. An increase in the NSA can cause larger motion artifacts and deteriorate the quality of strain images, and the improvement of strain estimation is limited when the NSA is increased from 3 to 7. An NSA of 3 is thus recommended to balance the influences of motion artifacts and the improvement for strain estimation. The MSA has little influence on the motion artifacts, while increased MSA can achieve improved lateral estimation performance at the cost of a smaller imaging region. In light of the lateral strain estimation performance and imaging region, an MSA of 15° is recommended. The influences of these key parameters obtained from this study may provide insights for parameter optimization in strain estimation with SAC to minimize the effects of motion artifacts.

Multicomponent material property characterization of atherosclerotic human carotid arteries through a Bayesian Optimization based inverse finite element approach

OBJECTIVE

Plaque rupture in atherosclerotic carotid arteries is a main cause of ischemic stroke and it is correlated with high plaque stresses. Hence, analyzing stress patterns is essential for plaque specific rupture risk assessment. However, the critical information of the multicomponent material properties of atherosclerotic carotid arteries is still lacking greatly. This work aims to characterize component-wise material properties of atherosclerotic human carotid arteries under (almost) physiological loading conditions.

METHODS

An inverse finite element modeling (iFEM) framework was developed to characterize fibrous intima and vessel wall material properties of 13 cross sections from five carotids. The novel pipeline comprised ex-vivo inflation testing, pre-clinical high frequency ultrasound for deriving plaque deformations, pre-clinical high-magnetic field magnetic resonance imaging, finite element modeling, and a sample efficient machine learning based Bayesian Optimization.

RESULTS

The nonlinear Yeoh constants for the fibrous intima and wall layers were successfully obtained. The optimization scheme of the iFEM reached the global minimum with a mean error of 3.8% in 133 iterations on average. The uniqueness of the results were confirmed with the inverted Gaussian Process (GP) model trained during the iFEM protocol.

CONCLUSION

The developed iFEM approach combined with the inverted GP model successfully predicted component-wise material properties of intact atherosclerotic human carotids ex-vivo under physiological-like loading conditions.

SIGNIFICANCE

We developed a novel iFEM framework for the nonlinear, component-wise material characterization of atherosclerotic arteries and utilized it to obtain human atherosclerotic carotid material properties. The developed iFEM framework has great potential to be advanced for patient-specific in-vivo application.

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.

4D cardiac electromechanical activation imaging

Cardiac abnormalities, a major cause of morbidity and mortality, affect millions of people worldwide. Despite the urgent clinical need for early diagnosis, there is currently no noninvasive technique that can infer to the electrical function of the whole heart in 3D and thereby localize abnormalities at the point of care. Here we present a new method for noninvasive 4D mapping of the cardiac electromechanical activity in a single heartbeat for heart disease characterization such as arrhythmia and infarction. Our novel technique captures the 3D activation wave of the heart in vivo using high volume-rate (500 volumes per second) ultrasound with a 32 × 32 matrix array. Electromechanical activation maps are first presented in a normal and infarcted cardiac model in silico and in canine heart during pacing and re-entrant ventricular tachycardia in vivo. Noninvasive 4D electromechanical activation mapping in a healthy volunteer and a heart failure patient are also determined. The technique described herein allows for direct, simultaneous and noninvasive visualization of electromechanical activation in 3D, which provides complementary information on myocardial viability and/or abnormality to clinical imaging.

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 Blood Velocity Vector Imaging Using Adaptive Velocity Compounding in the Carotid Artery Bifurcation

Visualization and quantification of blood flow are considered important for early detection of atherosclerosis and patient-specific diagnosis and intervention. As conventional Doppler imaging is limited to 1-D velocity estimates, 2-D and 3-D techniques are being developed. We introduce an adaptive velocity compounding technique that estimates the 2-D velocity vector field using predominantly axial displacements estimated by speckle tracking from dual-angle plane wave acquisitions. Straight-vessel experiments with a 7.8-MHz linear array transducer connected to a Verasonics Vantage ultrasound system revealed that the technique performed with a maximum velocity magnitude bias and angle bias of -3.7% (2.8% standard deviation) and -0.16° (0.41° standard deviation), respectively. In vivo, complex flow patterns were visualized in two healthy and three diseased carotid arteries and quantified using a vector complexity measure that increased with increasing wall irregularity. This measure could potentially be a relevant clinical parameter which might aid in early detection of atherosclerosis.

Interoperator Reproducibility of Carotid Elastography for Identification of Vulnerable Atherosclerotic Plaques

Ultrasound-based carotid elastography has been developed to evaluate the vulnerability of carotid atherosclerotic plaques. The aim of this study was to investigate the in vivo interoperator reproducibility of carotid elastography for the identification of vulnerable plaques, with high-resolution magnetic resonance imaging (MRI) as reference. Ultrasound radio-frequency data of 45 carotid arteries (including 53 plaques) from 32 volunteers were acquired separately by two experienced operators in the longitudinal view and then were used to estimate the interframe axial strain rate (ASR) with a two-step optical flow method. The maximum 99th percentile of absolute ASR of all plaques in a carotid artery was used as the elastographic index. MRI scanning was also performed on each volunteer to identify the vulnerable plaque. The results showed no systematic bias in the Bland-Altman plot and an intraclass correlation coefficient of 0.66 between the two operators. In addition, no statistical significance was found between the receiver operating characteristic (ROC) curves from the two operators ( ), and their areas under the ROC curves were 0.83 and 0.77, respectively. Using the mean measurements of the two operators as the classification criterion, a sensitivity of 71.4%, a specificity of 87.1%, and an accuracy of 82.2% were obtained with a cutoff value of 1.37 [Formula: see text]. This study validates the interoperator reproducibility of ultrasound-based carotid elastography for identifying vulnerable carotid plaques.

Investigation of out-of-plane motion artifacts in 2D noninvasive vascular ultrasound elastography

Ultrasound noninvasive vascular elastography (NIVE) has shown its potential to measure strains of carotid arteries to predict plaque instability. When two-dimensional (2D) strain estimation is performed, either in longitudinal or cross-sectional view, only in-plane motions are considered. The motions in elevation direction (i.e. perpendicular to the imaging plane), can induce estimation artifacts affecting the accuracy of 2D NIVE. The influence of such out-of-plane motions on the performance of axial strain and axial shear strain estimations has been evaluated in this study. For this purpose, we designed a diseased carotid bifurcation phantom with a 70% stenosis and an in vitro experimental setup to simulate orthogonal out-of-plane motions of 1 mm, 2 mm and 3 mm. The Lagrangian speckle model estimator (LSME) was used to estimate axial strains and shears under pulsatile conditions. As anticipated, in vitro results showed more strain estimation artifacts with increasing magnitudes of motions in elevation. However, even with an out-of-plane motion of 2.0 mm, strain and shear estimations having inter-frame correlation coefficients higher than 0.85 were obtained. To verify findings of in vitro experiments, a clinical LSME dataset obtained from 18 participants with carotid artery stenosis was used. Deduced out-of-plane motions (ranging from 0.25 mm to 1.04 mm) of the clinical dataset were classified into three groups: small, moderate and large elevational motions. Clinical results showed that pulsatile time-varying strains and shears remained reproducible for all motion categories since inter-frame correlation coefficients were higher than 0.70, and normalized cross-correlations (NCC) between radiofrequency (RF) images were above 0.93. In summary, the performance of LSME axial strain and shear estimations appeared robust in the presence of out-of-plane motions (<2 mm) as encountered during clinical ultrasound imaging.

3-D Single Breath-Hold Shear Strain Estimation for Improved Breast Lesion Detection and Classification in Automated Volumetric Ultrasound Scanners

Automated breast volume scanner (ABVS) is an ultrasound imaging modality used in breast cancer screening. It has high sensitivity but limited specificity as it is hard to discriminate between benign and malignant lesions by echogenic properties. Specificity might be improved by shear strain imaging as malignant lesions, firmly bonded to its host tissue, show different shear patterns compared to benign lesions, often loosely bonded. Therefore, 3-D quasi-static elastography was implemented in an ABVS-like system. Plane wave instead of conventional focused transmissions were used to reduce scan times within a single breath hold. A 3-D strain tensor was obtained and shear strains were reconstructed in phantoms containing firmly and loosely bonded lesions. Experiments were also simulated in finite-element models (FEMs). Experimental results, confirmed by FEM-results, indicated that loosely bonded lesions showed increased maximal shear strains (~2.5%) and different shear patterns compared to firmly bonded lesions (~0.9%). To conclude, we successfully implemented 3-D elastography in an ABVS-like system to assess lesion bonding by shear strain imaging.

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.

Strain imaging of the lateral collateral ligament using high frequency and conventional ultrasound imaging: An ex-vivo comparison

Recent first attempts of in situ ultrasound strain imaging in collateral ligaments encountered a number of challenges and illustrated a clear need for additional studies and more thorough validation of the available strain imaging methods. Therefore, in this study we experimentally validated ultrasound strain measurements of ex vivo human lateral collateral ligaments in an axial loading condition. Moreover, the use of high frequency ultrasound (>20 MHz) for strain measurement was explored and its performance compared to conventional ultrasound. The ligaments were stretched up to 5% strain and ultrasound measurements were compared to surface strain measurements from optical digital image correlation (DIC) techniques. The results show good correlations between ultrasound based and DIC based strain measures with R² values of 0.71 and 0.93 for high frequency and conventional ultrasound, subsequently. The performance of conventional ultrasound was significantly higher compared to high frequency ultrasound strain imaging, as the high frequency based method seemed more prone to errors. This study demonstrates that ultrasound strain imaging is feasible in ex vivo lateral collateral ligaments, which are relatively small structures. Additional studies should be designed for a more informed assessment of optimal in vivo strain measurements in collateral knee ligaments.

A dynamic ultrasound simulation of a pulsating three-layered CCA for validation of two-dimensional wall motion and blood velocity estimation algorithms

PURPOSE

A dynamic ultrasound simulation model for the common carotid artery (CCA) with three arterial layers for validation of two-dimensional wall motion and blood velocity estimation algorithms is proposed in the present study. This model describes layers with not only characteristics of echo distributions conforming to clinical ones but also varying thicknesses, axial, and radial displacements with pulsatile blood pressure during a cardiac cycle.

METHODS

The modeling process is as follows: first, a geometrical model according with the clinical structure size of a CCA is built based on the preset layer thicknesses and the diameter of lumen. Second, a three-dimensional scatterer model is constructed by a mapping with a Hilbert space-filling curve from the one-dimensional scatterer distribution with the position and amplitude following Gamma and Gaussian distributions, respectively. The characteristics of three layers and blood are depicted by smoothly adjusting the scatterer density, the scale, and shape parameters of the Gamma distribution as well as the mean and standard deviation of the Gaussian distribution. To obtain the values of parameters of scatterer distributions, including the shape parameter, density, and intensity, for arterial layers and blood, the envelope signals simulated from different configurations of scatterer distribution are compared with those from different kinds of tissue of CCAs in vivo through a statistic analysis. Finally, the dynamic scatterer model is realized based on the blood pressure, elasticity modulus of intima-media (IM) and adventitia, varying IM thickness, axial displacement of IM as well as blood flow velocity at central axis during a cardiac cycle. Then, the corresponding radiofrequency (RF) signals, envelope signals, and B-mode images of the pulsatile CCA are generated in a dynamic scanning mode using Field II platform.

RESULTS

The three arterial layers, blood, and surrounding tissue in simulated B-mode ultrasound images are clearly legible. The results based on a statistical analysis for the envelope signals from 30 simulations indicate that the echo characteristics of blood, intima, media, and adventitia are in accordant with clinical ones. The maximum relative errors between the preset geometrical sizes and the measured ones from the simulated images for the diameter of the lumen and the thicknesses of the intima, media, and adventitia are 0.13%, 3.89%, 1.35%, and 0.06%, respectively. For the dynamic parameters, the variation in IM thickness, the radial displacements of lumen and adventitia as well as the axial displacement of IM and blood flow velocity are measured with the mean relative errors of 68.03%, 9.27%, 2.10%, 4.93%, and 17.34%, respectively.

CONCLUSION

The simulated results present static sizes and dynamical variations according with preset values; echo distributions conforming to clinical versions. Therefore, the presented simulation model could be useful as a data source to evaluate the performance of studies on measurements of ultrasound-based tissue structures and dynamic parameters for the CCA layers.

Feasibility and Validation of 4-D Pulse Wave Imaging in Phantoms and In Vivo

Pulse wave imaging (PWI) is a noninvasive technique for tracking the propagation of the pulse wave along the arterial wall. The 3-D ultrasound imaging would aid in objectively estimating the pulse wave velocity (PWV) vector. This paper aims to introduce a novel PWV estimation method along the propagation direction, validate it in phantoms, and test its feasibility in vivo. A silicone vessel phantom consisting of a stiff and a soft segment along the longitudinal axis and a silicone vessel with a plaque were constructed. A 2-D array with a center frequency of 2.5 MHz was used. Propagation was successfully visualized in 3-D in each phantom and in vivo in six healthy subjects. In three of the healthy subjects, results were compared against conventional PWI using a linear array. PWVs were estimated in the stiff (3.42 ± 0.23 m [Formula: see text]) and soft (2.41 ± 0.07 m [Formula: see text]) phantom segments. Good agreement was found with the corresponding static testing values (stiff: 3.41 m [Formula: see text] and soft: 2.48 m [Formula: see text]). PWI-derived vessel compliance values were validated with dynamic testing. Comprehensive views of pulse propagation in the plaque phantom were generated and compared against conventional PWI acquisitions. Good agreement was found in vivo between the results of 4-D PWI (4.80 ± 1.32 m [Formula: see text]) and conventional PWI (4.28±1.20 m [Formula: see text]) ( n=3 ). PWVs derived for all of the healthy subjects ( n = 6 ) were within the physiological range. Thus, the 4-D PWI was successfully validated in phantoms and used to image the pulse wave propagation in normal human subjects in vivo.