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

Unsupervised deep learning-based displacement estimation for vascular elasticity imaging applications

Objective. Arterial wall stiffness can provide valuable information on the proper function of the cardiovascular system. Ultrasound elasticity imaging techniques have shown great promise as a low-cost and non-invasive tool to enable localized maps of arterial wall stiffness. Such techniques rely upon motion detection algorithms that provide arterial wall displacement estimation.Approach. In this study, we propose an unsupervised deep learning-based approach, originally proposed for image registration, in order to enable improved quality arterial wall displacement estimation at high temporal and spatial resolutions. The performance of the proposed network was assessed through phantom experiments, where various models were trained by using ultrasound RF signals, or B-mode images, as well as different loss functions.Main results. Using the mean square error (MSE) for the training process provided the highest signal-to-noise ratio when training on the B-modes images (30.36 ± 1.14 dB) and highest contrast-to-noise ratio when training on the RF signals (32.84 ± 1.89 dB). In addition, training the model on RF signals demonstrated the capability of providing accurate localized pulse wave velocity (PWV) maps, with a mean relative error (MREPWV) of 3.32 ± 1.80% and anR2 of 0.97 ± 0.03. Finally, the developed model was tested in human common carotid arteriesin vivo, providing accurate tracking of the distension pulse wave propagation, with an MREPWV= 3.86 ± 2.69% andR2 = 0.95 ± 0.03.Significance. In conclusion, a novel displacement estimation approach was presented, showing promise in improving vascular elasticity imaging techniques.

A New Plane Wave Compounding Scheme Using Phase Compensation for Motion Detection

Plane wave (PW) transmission has enabled multiple new applications, such as shear wave elastography, ultrafast Doppler imaging, and functional ultrasound imaging. PW compounding (PWC), which coherently sums the echo signals from multiple PW transmits with different angles, is widely used to improve B-mode image quality. When the motion between two speckle images is estimated, PWC suffers from an inherent displacement estimation error. This is derived theoretically and experimentally demonstrated in this work. We show that the phase difference between the acquired data with PW emissions with different angles is related to this error. When the absolute value of the phase difference is larger than π /2, the displacement estimation error occurs. A new scheme, named initial-phase-compensated PWC (IPCPWC), is proposed, which compensates the phase of echo signals from each PW transmit and maintains the absolute value of the phase difference smaller than π /2. The increased signal-to-noise ratio and reduced jitter of IPCPWC in motion data are demonstrated using tissue mimicking phantoms compared with PWC.

Feasibility of longitudinal monitoring of atherosclerosis with pulse wave imaging in a swine model

Objective.Atherosclerosis is a vascular disease characterized by compositional and mechanical changes in the arterial walls that lead to a plaque buildup. Depending on its geometry and composition, a plaque can ruptured and cause stroke, ischemia or infarction. Pulse wave imaging (PWI) is an ultrasound-based technique developed to locally quantify the stiffness of arteries. This technique has shown promising results when applied to patients. The objective of this study is to assess the capability of PWI to monitor the disease progression in a swine model that mimics human pathology.Approach.The left common carotid of three hypercholesterolemic Wisconsin miniature swines, fed an atherogenic diet, was ligated. Ligated and contralateral carotids were imaged once a month over 9 months, at a high-frame-rate, with a 5-plane wave compounding sequence and a 5 MHz linear array. Each acquisition was repeated after probe repositioning to evaluate the reproducibility. Wall displacements were estimated from the beamformed RF-data and were arranged as spatiotemporal maps depicting the wave propagation. The pulse wave velocity (PWV) estimated by tracking the 50% upstroke of the wave was converted in compliance using the Bramwell-Hill model. At the termination of the experiment, the carotids were extracted for histology analysis.Main results.PWI was able to monitor the evolution of compliance in both carotids of the animals. Reproducibility was demonstrated as the difference of PWV between cardiac cycles was similar to the difference between acquisitions (9.04% versus 9.91%). The plaque components were similar to the ones usually observed in patients. Each animal presented a unique pattern of compliance progression, which was confirmed by the plaque composition observed histologically.Significance.This study provides important insights on the vascular wall stiffness progression in an atherosclerotic swine model. It therefore paves the way for a thorough longitudinal study that examines the role of stiffness in both the plaque formation and plaque progression.

Feasibility of Bilinear Mechanical Characterization of the Abdominal Aorta in a Hypertensive Mouse Model

A change in elastin and collagen content is indicative of damage caused by hypertension, which changes the non-linear behavior of the vessel wall. This study was aimed at investigating the feasibility of monitoring the non-linear material behavior in an angiotensin II hypertensive mice model. Aortas from 13 hypertensive mice were imaged with pulse wave imaging (PWI) over 4 wk using a 40-MHz linear array. The pulse wave velocity was estimated using two wave features: (i) the maximum axial acceleration of the foot (PWV_dia) and (ii) the maximum axial acceleration of the dicrotic notch (PWV_end-sys). The Bramwell-Hill equation was used to derive the compliance at diastolic and end-systolic pressure. This study determined the potential of PWI in a hypertensive mouse model to image and quantify the non-linear material behavior in vivo. End-systolic compliance could differentiate between the sham and angiotensin II groups, whereas diastolic compliance could not, indicating that PWI can detect early collagen-dominated remodeling.

Time-Aligned Plane Wave Compounding Methods for High-Frame-Rate Shear Wave Elastography: Experimental Validation and Performance Assessment on Tissue Phantoms

Shear wave elastography (SWE) is an ultrasonic technique able to quantitatively assess the mechanical properties of tissues by combining acoustic radiation force and ultrafast imaging. While utilizing coherent plane wave compounding enhances echo and shear wave motion signal-to-noise ratio (SNR), it also reduces the effective pulse repetition frequency (PRF_e), affecting the accuracy of the measurements of motion and, consequently, of material properties. It is important to maintain both high-motion SNR and PRF_e, particularly for the characterization of (material and/or geometrical) dispersive tissues such as arteries. This work proposes a method for SWE measurements with high SNR, while maintaining a high PRF_e, using conventional clinical ultrasound scanners. A time alignment process is applied after acquiring data from plane wave transmissions at different angles. The time alignment uses interpolation to obtain data points at higher frame rates, and the time-aligned data are compounded to increase the SNR. The method is used for SWE in tissue-mimicking phantoms of different stiffness and is compared with traditional plane wave compounding. Increases of 58% and 36% in spatial and temporal bandwidth compared with conventional plane wave compounding, respectively, can be achieved for SWE measurements of representative arterial stiffness values. Improvements in phase velocity accuracy and bandwidth in an arterial phantom are also described, to emphasize the beneficial advantage in dispersive cases.

Pulse Wave Imaging Coupled With Vector Flow Mapping: A Phantom, Simulation, and In Vivo Study

Pulse wave imaging (PWI) is an ultrasound imaging modality that estimates the wall stiffness of an imaged arterial segment by tracking the pulse wave propagation. The aim of the present study is to integrate PWI with vector flow imaging, enabling simultaneous and co-localized mapping of vessel wall mechanical properties and 2-D flow patterns. Two vector flow imaging techniques were implemented using the PWI acquisition sequence: 1) multiangle vector Doppler and 2) a cross-correlation-based vector flow imaging (CC VFI) method. The two vector flow imaging techniques were evaluated in vitro using a vessel phantom with an embedded plaque, along with spatially registered fluid structure interaction (FSI) simulations with the same geometry and inlet flow as the phantom setup. The flow magnitude and vector direction obtained through simulations and phantom experiments were compared in a prestenotic and stenotic segment of the phantom and at five different time frames. In most comparisons, CC VFI provided significantly lower bias or precision than the vector Doppler method ( ) indicating better performance. In addition, the proposed technique was applied to the carotid arteries of nonatherosclerotic subjects of different ages to investigate the relationship between PWI-derived compliance of the arterial wall and flow velocity in vivo. Spearman's rank-order test revealed positive correlation between compliance and peak flow velocity magnitude ( rs = 0.90 and ), while significantly lower compliance ( ) and lower peak flow velocity magnitude ( ) were determined in older (54-73 y.o.) compared with young (24-32 y.o.) subjects. Finally, initial feasibility was shown in an atherosclerotic common carotid artery in vivo. The proposed imaging modality successfully provided information on blood flow patterns and arterial wall stiffness and is expected to provide additional insight in studying carotid artery biomechanics, as well as aid in carotid artery disease diagnosis and monitoring.

Helical toroid phantom for 3D flow imaging investigations

The medical physics community has hitherto lacked an effective calibration phantom to holistically evaluate the performance of three-dimensional (3D) flow imaging techniques. Here, we present the design of a new omnidirectional, three-component (3-C) flow phantom whose lumen is consisted of a helical toroid structure (4 mm lumen diameter; helically winded for 5 revolutions over a torus with 10 mm radius; 5 mm helix radius). This phantom's intraluminal flow trajectory embraces all combinations of x, y, and z directional components, as confirmed using computational fluid dynamics (CFD) simulations. The phantom was physically fabricated via lost-core casting with polyvinyl alcohol cryogel (PVA) as the tissue mimic. 3D ultrasound confirmed that the phantom lumen expectedly resembled a helical toroid geometry. Pulsed Doppler measurements showed that the phantom, when operating under steady flow conditions (3 ml s^-1 flow rate), yielded flow velocity magnitudes that agreed well with those derived from CFD at both the inner torus (-47.6 ± 5.7 versus -52.0 ± 2.2 cm s^-1; mean ± 1 S.D.) and the outer torus (49.5 ± 4.2 versus 48.0 ± 1.7 cm s^-1). Additionally, 3-C velocity vectors acquired from multi-angle pulsed Doppler showed good agreement with CFD-derived velocity vectors (<7% and 10° difference in magnitude and flow angle, respectively). Ultrasound color flow imaging further revealed that the phantom's axial flow pattern was aligned with the CFD-derived flow profile. Overall, the helical toroid phantom has strong potential as an investigative tool in 3D flow imaging innovation endeavors, such as the development of flow vector estimators and visualization algorithms.

Receive/Transmit Aperture Selection for 3D Ultrasound Imaging with a 2D Matrix Transducer

Recently, we realized a prototype matrix transducer consisting of 48 rows of 80 elements on top of a tiled set of Application Specific Integrated Circuits (ASICs) implementing a row-level control connecting one transmit/receive channel to an arbitrary subset of elements per row. A fully sampled array data acquisition is implemented by a column-by-column (CBC) imaging scheme (80 transmit-receive shots) which achieves 250 volumes/second (V/s) at a pulse repetition frequency of 20 kHz. However, for several clinical applications such as carotid pulse wave imaging (CPWI), a volume rate of 1000 per second is needed. This allows only 20 transmit-receive shots per 3D image. In this study, we propose a shifting aperture scheme and investigate the effects of receive/transmit aperture size and aperture shifting step in the elevation direction. The row-level circuit is used to interconnect elements of a receive aperture in the elevation (row) direction. An angular weighting method is used to suppress the grating lobes caused by the enlargement of the effective elevation pitch of the array, as a result of element interconnection in the elevation direction. The effective aperture size, level of grating lobes, and resolution/sidelobes are used to select suitable reception/transmission parameters. Based on our assessment, the proposed imaging sequence is a full transmission (all 80 elements excited at the same time), a receive aperture size of 5 and an aperture shifting step of 3. Numerical results obtained at depths of 10, 15, and 20 mm show that, compared to the fully sampled array, the 1000 V/s is achieved at the expense of, on average, about two times wider point spread function and 4 dB higher clutter level. The resulting grating lobes were at −27 dB. The proposed imaging sequence can be used for carotid pulse wave imaging to generate an informative 3D arterial stiffness map, for cardiovascular disease assessment.

3-D Super-Resolution Ultrasound Imaging With a 2-D Sparse Array

High-frame-rate 3-D ultrasound imaging technology combined with super-resolution processing method can visualize 3-D microvascular structures by overcoming the diffraction-limited resolution in every spatial direction. However, 3-D super-resolution ultrasound imaging using a full 2-D array requires a system with a large number of independent channels, the design of which might be impractical due to the high cost, complexity, and volume of data produced. In this study, a 2-D sparse array was designed and fabricated with 512 elements chosen from a density-tapered 2-D spiral layout. High-frame-rate volumetric imaging was performed using two synchronized ULA-OP 256 research scanners. Volumetric images were constructed by coherently compounding nine-angle plane waves acquired at a pulse repetition frequency of 4500 Hz. Localization-based 3-D super-resolution images of two touching subwavelength tubes were generated from 6000 volumes acquired in 12 s. Finally, this work demonstrates the feasibility of 3-D super-resolution imaging and super-resolved velocity mapping using a customized 2-D sparse array transducer.

Arterial wall mechanical inhomogeneity detection and atherosclerotic plaque characterization using high frame rate pulse wave imaging in carotid artery disease patients in vivo

Pulse wave imaging (PWI) is a non-invasive, ultrasound-based technique, which provides information on arterial wall stiffness by estimating the pulse wave velocity (PWV) along an imaged arterial wall segment. The aims of the present study were to: (1) utilize the PWI information to automatically and optimally divide the artery into the segments with most homogeneous properties and (2) assess the feasibility of this method to provide arterial wall mechanical characterization in normal and atherosclerotic carotid arteries in vivo. A silicone phantom consisting of a soft and stiff segment along its longitudinal axis was scanned at the stiffness transition, and the PWV in each segment was estimated through static testing. The proposed algorithm detected the stiffness interface with an average error of 0.98  ±  0.49 mm and 1.04  ±  0.27 mm in the soft-to-stiff and stiff-to-soft pulse wave transmission direction, respectively. Mean PWVs estimated in the case of the soft-to-stiff pulse wave transmission direction were 2.47 [Formula: see text] 0.04 m s^-1 and 3.43 [Formula: see text] 0.08 m s^-1 for the soft and stiff phantom segments, respectively, while in the case of stiff-to-soft transmission direction PWVs were 2.60 [Formula: see text] 0.18 m s^-1 and 3.72 [Formula: see text] 0.08 m s^-1 for the soft and stiff phantom segments, respectively, which were in good agreement with the PWVs obtained through static testing (soft segment: 2.41 m s^-1, stiff segment: 3.52 m s^-1). Furthermore, the carotid arteries of N  =  9 young subjects (22-32 y.o.) and N  =  9 elderly subjects (60-73 y.o.) with no prior history of carotid artery disease were scanned, in vivo, as well as the atherosclerotic carotid arteries of N  =  12 (59-85 y.o.) carotid artery disease patients. One-way ANOVA with Holm-Sidak correction showed that the number of most homogeneous segments in which the artery was divided was significantly higher in the case of carotid artery disease patients compared to young (3.25 [Formula: see text] 0.86 segments versus 1.00 [Formula: see text] 0.00 segments, p -value  <  0.0001) and elderly non-atherosclerotic subjects (3.25 [Formula: see text] 0.86 segments versus 1.44 [Formula: see text] 0.51 segments p -value  <  0.0001), indicating increased wall inhomogeneity in atherosclerotic arteries. The compliance provided by the proposed algorithm was significantly higher in non-calcified/high-lipid plaques as compared with calcified plaques (3.35 [Formula: see text] 2.45 *[Formula: see text] versus 0.22 [Formula: see text] 0.18 * [Formula: see text], p -value  <  0.01) and the compliance estimated in elderly subjects (3.35 [Formula: see text] 2.45 * [Formula: see text] versus 0.79 [Formula: see text] 0.30 * [Formula: see text], p -value  <  0.01). Moreover, lower compliance was estimated in cases where vulnerable plaque characteristics were present (i.e. necrotic lipid core, thrombus), compared to stable plaque components (calcification), as evaluated through plaque histological examination. The proposed algorithm was thus capable of evaluating arterial wall inhomogeneity and characterize wall mechanical properties, showing promise in vascular disease diagnosis and monitoring.

Pulse Wave Imaging in Carotid Artery Stenosis Human Patients in Vivo

Carotid stenosis involves narrowing of the lumen in the carotid artery potentially leading to a stroke, which is the third leading cause of death in the United States. Several recent investigations have found that plaque structure and composition may represent a more direct biomarker of plaque rupture risk compared with the degree of stenosis. In this study, pulse wave imaging was applied in 111 (n = 11, N = 13 plaques) patients diagnosed with moderate (>50%) to severe (>80%) carotid artery stenosis to investigate the feasibility of characterizing plaque properties based on the pulse wave-induced arterial wall dynamics captured by pulse wave imaging. Five (n = 5 patients, N = 20 measurements) healthy volunteers were also imaged as a control group. Both conventional and high-frame-rate plane wave radiofrequency imaging sequences were used to generate piecewise maps of the pulse wave velocity (PWV) at a single depth along stenotic carotid segments, as well as intra-plaque PWV mapping at multiple depths. Intra-plaque cumulative displacement and strain maps were also calculated for each plaque region. The Bramwell-Hill equation was used to estimate the compliance of the plaque regions based on the PWV and diameter. Qualitatively, wave convergence, elevated PWV and decreased cumulative displacement around and/or within regions of atherosclerotic plaque were observed and may serve as biomarkers for plaque characterization. Intra-plaque mapping revealed the potential to capture wave reflections between calcified inclusions and differentiate stable (i.e., calcified) from vulnerable (i.e., lipid) plaque components based on the intra-plaque PWV and cumulative strain. Quantitatively, one-way analysis of variance indicated that the pulse wave-induced cumulative strain was significantly lower (p < 0.01) in the moderately and severely calcified plaques compared with the normal controls. As expected, compliance was also significantly lower in the severely calcified plaques regions compared with the normal controls (p < 0.01). The results from this pilot study indicated the potential of pulse wave imaging coupled with strain imaging to differentiate plaques of varying stiffness, location and composition. Such findings may serve as valuable information to compensate for the limitations of currently used methods for the assessment of stroke risk.