Displacement Encoding With Stimulated Echoes MRI Research Articles

Correction of phase offset errors and quantification of background noise, signal-to-noise ratio, and encoded-displacement uncertainty on DENSE MRI for kinematics of the descending thoracic and abdominal aorta

Displacement encoding with stimulated echoes (DENSE) MRI is a phase contrast technique that allows the encoding of tissue displacement into the phase of the magnetic resonance signal. Recent developments in this technique allow the imaging of relatively thin structures such as the aortic wall. Quantifying background noise associated to DENSE MRI is required to assess the uncertainty of derived displacement measurements and for the design and implementation of adequate noise-reduction techniques. Although noise and error management of cardiac DENSE MRI has been previously studied, developments for aortic applications are scarce. Herein, we evaluate the noise and uncertainty of DENSE MRI scans at three different locations along the descending aorta: the distal aortic arch (DAA), the descending thoracic aorta (DTA), and infrarenal abdominal aorta (IAA). Additionally, we analyze three datasets from in vitro validation experiments with polyvinyl alcohol phantoms. We implement and evaluate the effectiveness of an offset-error correction algorithm and noise filtering techniques on DENSE MRI for aortic motion applications. Our results show that the phase signal of pixels composing the static background was normally distributed, centered on average at 0.003 ± 0.02 rad and − 0.02 ± 0.024 rad for each phase directions, suggesting that background noise is random, isotropic, and DENSE MRI has little offset errors. However, background signal noise significantly increased with elapsed time of the cardiac cycle; and was spatially heterogeneous consistently increased towards the anterior space. Background noise showed no significant differences between the 3 aortic locations and the in vitro experiments. However, SNR depended on the displacement of the region of interest, in consequence it was found significantly larger at DAA (16.7 ± 8.5, p = 0.003) and DTA (15.4 ± 7.6, p = 0.008) than at the IAA (8.0 ± 4.1), but not significantly different than the SNR of in vitro experiments (8.0 ± 3.7), and had an overall average of 13 ± 7. The applied methods significantly reduced the offset error and effect of noise on the estimation of encoded displacements. Finally, this analysis suggests that the implemented DENSE MRI protocol is adequate to assess the motion of healthy human aortas. However, the relative effect of noise increased considerably on the analysis of an ageing or diseased aortas with impaired mobility, calling for further analyses on pathologically stiffened aortas.

Multi-frame biomechanical and relaxometry analysis during in vivo loading of the human knee by spiral dualMRI and compressed sensing.

Knee cartilage experiences repetitive loading during physical activities, which is altered during the pathogenesis of diseases like osteoarthritis. Analyzing the biomechanics during motion provides a clear understanding of the dynamics of cartilage deformation and may establish essential imaging biomarkers of early-stage disease. However, in vivo biomechanical analysis of cartilage during rapid motion is not well established. We used spiral displacement encoding with stimulated echoes (DENSE) MRI on in vivo human tibiofemoral cartilage during cyclic varus loading (0.5 Hz) and used compressed sensing on the k-space data. The applied compressive load was set for each participant at 0.5 times body weight on the medial condyle. Relaxometry methods were measured on the cartilage before (T1ρ , T2 ) and after (T1ρ ) varus load. Displacement and strain maps showed a gradual shift of displacement and strain in time. Compressive strain was observed in the medial condyle cartilage and shear strain was roughly half of the compressive strain. Male participants had more displacement in the loading direction compared to females, and T1ρ values did not change after cyclic varus load. Compressed sensing reduced the scanning time up to 25% to 40% when comparing the displacement maps and substantially lowered the noise levels. These results demonstrated the ease of which spiral DENSE MRI could be applied to clinical studies because of the shortened imaging time, while quantifying realistic cartilage deformations that occur through daily activities and that could serve as biomarkers of early osteoarthritis.

Improved computation of Lagrangian tissue displacement and strain for cine DENSE MRI using a regularized spatiotemporal least squares method.

In displacement encoding with stimulated echoes (DENSE), tissue displacement is encoded in the signal phase such that the phase of each pixel in space and time provides an independent measurement of absolute tissue displacement. Previously for DENSE, estimation of Lagrangian displacement used two steps: first a spatial interpolation and, second, least squares fitting through time to a Fourier or polynomial model. However, there is no strong rationale for such a through-time model. To compute the Lagrangian displacement field from DENSE phase data, a minimization problem is introduced to enforce fidelity with the acquired Eulerian displacement data while simultaneously providing model-independent regularization in space and time, enforcing only spatiotemporal smoothness. A regularized spatiotemporal least squares (RSTLS) method is used to solve the minimization problem, and RSTLS was tested using two-dimensional DENSE data from 71 healthy volunteers. The mean absolute percent error (MAPE) between the Lagrangian displacements and the corresponding Eulerian displacements was significantly lower for the RSTLS method vs. the two-step method for both x- and y-directions (0.73±0.59 vs 0.83 ±0.1, p < 0.05) and (0.75±0.66 vs 0.82 ±0.1, p < 0.05), respectively. Also, peak early diastolic strain rate (PEDSR) was higher (1.81±0.58 (s-1) vs. 1.56±0. 63 (s-1), p<0.05) and the strain rate during diastasis was lower (0.14±0.18 (s-1) vs 0.35±0.2 (s-1), p < 0.05) for the RSTLS vs. the two-step method, with the former suggesting that the two-step method was over-regularized. The proposed RSTLS method provides more realistic measurements of Lagrangian displacement and strain from DENSE images without imposing arbitrary motion models.

Quantification of the heterogeneous effect of static and dynamic perivascular structures on patient-specific local aortic wall mechanics using inverse finite element modeling and DENSE MRI

Recent studies have highlighted the relevance of perivascular interactions on aortic wall mechanics. Most of the approaches assume static perivascular structures; however, the beating heart dynamically displaces the neighboring aorta. We develop a model to account for the effect of periaortic interactions due to static and dynamic structures by prescribing a moving elastic foundation boundary condition (EFBC) embedded into an inverse finite element algorithm using in vivo displacements from 2D displacement encoding with stimulated echoes (DENSE) MRI as target data. We applied this method at three different locations of interest, the distal aortic arch (DAA), descending thoracic aorta (DTA), and infrarenal abdominal aorta (IAA) for a total of 27 cases in healthy humans. The model reproduces the target diastole-to-systole deformation and bulk displacement of the aortic wall with median displacement errors below 0.5mm. The EFBC showed good agreement with the location of anatomical features and was consistent among individuals of similar characteristics. Results show that an energy source acting on the adventitia is required to reproduce the displacements measured at the vicinity of the heart, but not at the abdomen. The average adventitial load as a percentage of the luminal pulse-pressure was found to increase with age and to decrease along the descending aorta, from 61% at the DAA to 37% at the DTA, and 30% at the IAA. This approach offers a patient-specific method to estimate in vivo adventitial loads and aortic wall stiffness, which can bring a better understanding of normal and pathological in vivo aortic function.

A 3D model of the soleus reveals effects of aponeuroses morphology and material properties on complex muscle fascicle behavior

The soleus is an important plantarflexor muscle with complex fascicle and connective tissue arrangement. In this study we created an image-based finite element model representing the 3D structure of the soleus muscle and its aponeurosis connective tissue, including distinct fascicle architecture of the posterior and anterior compartments. The model was used to simulate passive and active soleus lengthening during ankle motion to predict tissue displacements and fascicle architecture changes. Both the model’s initial architecture and changes incurred during passive lengthening were consistent with prior in vivo data from diffusion tensor imaging. Model predictions of active lengthening were consistent with axial plane muscle displacements that we measured in eight subjects’ lower legs using cine DENSE (Displacement Encoding with Stimulated Echoes) MRI during eccentric dorsiflexion. Regional strains were variable and nonuniform in the model, but average fascicle strains were similar between the compartments for both passive (anterior: 0.18 ± 0.06, posterior: 0.19 ± 0.05) and active (anterior: 0.12 ± 0.05, posterior: 0.13 ± 0.06) lengthening and were two- to three-times greater than muscle belly strain (0.06). We used additional model simulations to investigate the effects of aponeurosis material properties on muscle deformation, by independently varying the longitudinal or transverse stiffness of the posterior or anterior aponeurosis. Results of model variations elucidate how properties of soleus aponeuroses contribute to fascicle architecture changes. Greater longitudinal stiffness of posterior compared to anterior aponeurosis promoted more uniform spatial distribution of muscle tissue deformation. Reduced transverse stiffness in both aponeuroses resulted in larger differences between passive and active soleus lengthening.

Cerebellar and Brainstem Displacement Measured with DENSE MRI in Chiari Malformation Following Posterior Fossa Decompression Surgery

Background Posterior fossa decompression (PFD) surgery is a treatment for Chiari malformation type I (CMI). The goals of surgery are to reduce cerebellar tonsillar crowding and restore posterior cerebral spinal fluid flow, but regional tissue biomechanics may also change. MRI-based displacement encoding with stimulated echoes (DENSE) can be used to assess neural tissue displacement. Purpose To assess neural tissue displacement by using DENSE MRI in participants with CMI before and after PFD surgery and examine associations between tissue displacement and symptoms. Materials and Methods In a prospective, HIPAA-compliant study of patients with CMI, midsagittal DENSE MRI was performed before and after PFD surgery between January 2017 and June 2020. Peak tissue displacement over the cardiac cycle was quantified in the cerebellum and brainstem, averaged over each structure, and compared before and after surgery. Paired t tests and nonparametric Wilcoxon signed-rank tests were used to identify surgical changes in displacement, and Spearman correlations were determined between tissue displacement and presurgery symptoms. Results Twenty-three participants were included (mean age ± standard deviation, 37 years ± 10; 19 women). Spatially averaged (mean) peak tissue displacement demonstrated reductions of 46% (79/171 µm) within the cerebellum and 22% (46/210 µm) within the brainstem after surgery (P < .001). Maximum peak displacement, calculated within a circular 30-mm2 area, decreased by 64% (274/427 µm) in the cerebellum and 33% (100/300 µm) in the brainstem (P < .001). No significant associations were identified between tissue displacement and CMI symptoms (r < .74 and P > .012 for all; Bonferroni-corrected P = .0002). Conclusion Neural tissue displacement was reduced after posterior fossa decompression surgery, indicating that surgical intervention changes brain tissue biomechanics. For participants with Chiari malformation type I, no relationship was identified between presurgery tissue displacement and presurgical symptoms. © RSNA, 2021 Online supplemental material is available for this article.

Suppression of artifact-generating echoes in cine DENSE using deep learning.

To use deep learning for suppression of the artifact-generating T1 -relaxation echo in cine displacement encoding with stimulated echoes (DENSE) for the purpose of reducing the scan time. A U-Net was trained to suppress the artifact-generating T1 -relaxation echo using complementary phase-cycled data as the ground truth. A data-augmentation method was developed that generates synthetic DENSE images with arbitrary displacement-encoding frequencies to suppress the T1 -relaxation echo modulated for a range of frequencies. The resulting U-Net (DAS-Net) was compared with k-space zero-filling as an alternative method. Non-phase-cycled DENSE images acquired in shorter breath-holds were processed by DAS-Net and compared with DENSE images acquired with phase cycling for the quantification of myocardial strain. The DAS-Net method effectively suppressed the T1 -relaxation echo and its artifacts, and achieved root Mean Square(RMS) error = 5.5 ± 0.8 and structural similarity index = 0.85 ± 0.02 for DENSE images acquired with a displacement encoding frequency of 0.10 cycles/mm. The DAS-Net method outperformed zero-filling (root Mean Square error = 5.8 ± 1.5 vs 13.5 ± 1.5, DAS-Net vs zero-filling, P < .01; and structural similarity index = 0.83 ± 0.04 vs 0.66 ± 0.03, DAS-Net vs zero-filling, P < .01). Strain data for non-phase-cycled DENSE images with DAS-Net showed close agreement with strain from phase-cycled DENSE. The DAS-Net method provides an effective alternative approach for suppression of the artifact-generating T1 -relaxation echo in DENSE MRI, enabling a 42% reduction in scan time compared to DENSE with phase-cycling.

Demonstration of circumferential heterogeneity in displacement and strain in the abdominal aortic wall by spiral cine DENSE MRI.

Knowledge of tissue properties of the abdominal aorta can improve understanding of vascular disease and guide interventional approaches. Existing MRI methods to quantify aortic wall displacement and strain are unable to discern circumferential heterogeneity. To assess regional variation in abdominal aortic wall displacement and strain as a function of circumferential position using spiral cine displacement encoding with stimulated echoes (DENSE). Prospective. Cardiovascular disease-free men (n = 8) and women (n = 9) ages 30-42. Prospective electrocardiogram (ECG)-gated and navigator echo-gated spiral, cine 2D DENSE and retrospective ECG-gated phase contrast MR (PCMR) sequences at 3T. In-plane displacement values of the aortic wall acquired with DENSE were used to determine radial and circumferential aortic wall motion. A quadrilateral-based 2D strain calculation method was implemented to determine strain from the displacement field. Peak displacement and its radial and circumferential contributions as well as peak circumferential strain were compared among eight circumferential wall segments. Distensibility was calculated using PCMR and compared with homogenized circumferential strain. To account for repeated measurements in volunteers, linear mixed models for mean sector values were created for displacement magnitude, circumferential displacement, radial displacement, and circumferential strain. Comparisons were made between sectors. Calculated distensibility and homogenized circumferential strain were compared using Bland-Altman analysis. Statistical significance was defined as P < 0.05. Displacement was highest in the anterior wall (1.5 ± 0.7 mm) and was primarily in the radial as compared with circumferential direction (1.04 ± 0.05 mm vs. 0.81 ± 0.42 mm). Circumferential strain was highest in the lateral walls (left 0.16 ± 0.05 and right 0.21 ± 0.12) with homogenized circumferential strain of 0.14 ± 0.05. DENSE imaging in the abdominal aortic wall demonstrated that the anterior aortic wall exhibits the greatest displacement, while the lateral wall experiences the largest circumferential strain. 3 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;49:731-743.

Optimal configuration of respiratory navigator gating for the quantification of left ventricular strain using spiral cine displacement encoding with stimulated echoes (DENSE) MRI.

To determine the optimal respiratory navigator gating configuration for the quantification of left ventricular strain using spiral cine displacement encoding with stimulated echoes (DENSE) MRI. Two-dimensional spiral cine DENSE was performed on a 3 Tesla MRI using two single-navigator configurations (retrospective, prospective) and a combined "dual-navigator" configuration in 10 healthy adults and 20 healthy children. The adults also underwent breathhold DENSE as a reference standard for comparisons. Peak left ventricular strains, signal-to-noise ratio (SNR), and navigator efficiency were compared. Subjects also underwent dual-navigator gating with and without visual feedback to determine the effect on navigator efficiency. There were no differences in circumferential, radial, and longitudinal strains between navigator-gated and breathhold DENSE (P = 0.09-0.95) (as confidence intervals, retrospective: [-1.0%-1.1%], [-7.4%-2.0%], [-1.0%-1.2%]; prospective: [-0.6%-2.7%], [-2.8%-8.3%], [-0.3%-2.9%]; dual: [-1.6%-0.5%], [-8.3%-3.2%], [-0.8%-1.9%], respectively). The dual configuration maintained SNR compared with breathhold acquisitions (16 versus 18, P = 0.06). SNR for the prospective configuration was lower than for the dual navigator in adults (P = 0.004) and children (P < 0.001). Navigator efficiency was higher (P < 0.001) for both retrospective (54%) and prospective (56%) configurations compared with the dual configuration (35%). Visual feedback improved the dual configuration navigator efficiency to 55% (P < 0.001). When quantifying left ventricular strains using spiral cine DENSE MRI, a dual navigator configuration results in the highest SNR in adults and children. In adults, a retrospective configuration has good navigator efficiency without a substantial drop in SNR. Prospective gating should be avoided because it has the lowest SNR. Visual feedback represents an effective option to maintain navigator efficiency while using a dual navigator configuration. 2 J. Magn. Reson. Imaging 2017;45:786-794.

Accelerated cine DENSE MRI using compressed sensing and parallel imaging

Background Cine DENSE (Displacement Encoding with Stimulated Echoes) provides accurate and high-resolution displacement and strain imaging of the heart; however, image acquisition times are relatively long and, due to properties inherent to stimulated echoes, signal-to-noise ratio (SNR) is relatively low. Accelerated cine DENSE could substantially shorten scan times and/or provide improved spatial and/or temporal resolution; however, the low SNR, requirement to preserve phase information, and cardiac motion present technical challenges. We aimed to develop acceleration methods that overcome these challenges.

Prinicipal component analysis of myocardial strains to optimize cardiac resynchronization therapy patient selection

Background Cardiac resynchronization therapy (CRT) is effective for selected heart failure (HF) patients but improved methods for identification of mechanical dyssynchrony are needed. Cine displacement encoding with stimulated echoes (DENSE) MRI can be used to obtain high quality strain data for LV dyssynchrony quantification. An ideal dyssynchrony parameter should be time-independent and require minimal user interaction. We hypothesize that the time-independent, data-driven quantification of dyssynchrony using principal component analysis (PCA) will help improve patient selection for CRT. Methods 2D cine DENSE MRI was performed in multiple shortaxis views to obtain LV circumferential strains (Ecc) for HF patients referred for CRT. In order to minimize confounding factors with an influence on CRT response, all patients in this analysis (N = 25) were required to have absence of scar as determined from late gadoliniumenhanced (LGE) MRI, as well as LV lead placement in a late-activated segment with time from QRS onset to the LV lead electrogram at the final lead position at least 50% of the QRS width. CRT response was defined as at least a 15% reduction in the LV end systolic volume 6 months after CRT implantation. Typical cine DENSE imaging parameters for HF patients include field of view = 340 - 400 mm 2 , matrix = 128 × 128, slice thickness = 8 mm, no. of spiral interleaves = 6, fat suppression, temporal resolution = 17 ms and displacement encoding frequency = 0.1 cycles/mm 1 .T he LV Ecc curves were

Emulation of human and rodent cardiac motion with a computer‐controlled cardiac phantom using DENSE MRI

AbstractIn this work, human myocardial motion is studied in an electrical and in an MRI‐compatible (pneumatic) version of a specialized, commercially available, human cardiac phantom at 1.5 Tesla using Displacement Encoding with Stimulated Echoes (DENSE). Work is extended to a prototype rodent phantom, designed, manufactured, implemented, and imaged on a high field (7 T) scanner using urethane‐elastomers under dynamic conditions. Mechanical properties of phantom composition (including commercially available urethane‐elastomers of variable elasticity, poly(glycerol) sebacate (PGS), and polyvinyl acetate (PVA)‐elastomeric samples were evaluated with dynamic mechanical analysis (DMA) and atomic force microscopy (AFM). Both phantoms were cyclically driven at 90–100 bpm and were molded to match accurately human and rodent anatomy. Constructed waveforms attained human and rodent torsions that ranged between 0–40° and 0–20°, respectively. The human phantom used PVA elastomer‐materials, whereas solid, cylindrical urethane‐elastomeric, and PGS samples were tested in the rodent phantom. Measured bulk storage moduli of elasticity varied between 0.7 and 1.06 MPa using the DMA, whereas AFM measurements independently confirmed the graded surface stiffness of elastomers and PGS samples. Displacement maps were generated with DENSE at time intervals ending at 25% of the cardiac cycle (due to the reduced DENSE‐MRI SNR at subsequent intervals) and yielded basal, middle, and apical longitudinal displacements that ranged between −4–6, −8–8, and −10–8 mm, respectively. Elicited in‐plane strain results yielded LV phantom values during the first five cardiac phases that ranged between 0.03–3.7% for Exx and 0.03–4.8% for Eyy. © 2013 Wiley Periodicals, Inc. Concepts Magn Reson Part A 42A: 59–71, 2013.

In vivo SNR in DENSE MRI; temporal and regional effects of field strength, receiver coil sensitivity and flip angle strategies

Aim The influences on the signal-to-noise ratio (SNR) of Displacement ENcoding with Stimulated Echoes (DENSE) MRI of field strength, receiver coil sensitivity and choice of flip angle strategy have been previously investigated individually. In this study, all of these parameters have been investigated in the same setting, and a mutual comparison of their impact on SNR is presented. Materials and methods Ten healthy volunteers were imaged in a 1.5 T and a 3 T MRI system, using standard five- or six-channel cardiac coils as well as 32-channel coils, with four different excitation patterns. Variation of spatial coil sensitivity was assessed by regional SNR analysis. Results SNR ranging from 2.8 to 30.5 was found depending on the combination of excitation patterns, coil sensitivity and field strength. The SNR at 3 T was 53±26% higher than at 1.5 T ( P<.001), whereas spatial differences of 59±26% were found in the ventricle ( P<.001). Thirty-two-channel coils provided 52±29% higher SNR compared to standard five- or six-channel coils ( P<.001). A fixed flip angle strategy provided an excess of 50% higher SNR in half of the imaged cardiac cycle compared to a sweeping flip angle strategy, and a single-phase acquisition provided a sixfold increase of SNR compared to a cine acquisition. Conclusion The effect of field strength and receiver coil sensitivity influences the SNR with the same order of magnitude, whereas flip angle strategy can have a larger effect on SNR. Thus, careful choice of imaging hardware in combination with adaptation of the acquisition protocol is crucial in order to realize sufficient SNR in DENSE MRI.

Cine DENSE MRI for circumferential and radial dyssynchrony in patients referred for cardiac resynchronization therapy

Cine DENSE (Displacement Encoding with Stimulated Echoes) MRI offers high resolution circumferential, radial, and longitudinal strain assessment. Dyssynchrony assessment using midwall circumferential strain corresponds to the predominant physiologic orientation of myofibers, and its effectiveness for cardiac resynchronization therapy (CRT) patient selection has recently been demonstrated [1]. Current echocardiographic dyssynchrony measures are based on radial or longitudinal velocity/strain.

Motion-guided segmentation for cine DENSE MRI.

Defining myocardial contours is often the most time-consuming portion of dynamic cardiac MRI image analysis. Displacement encoding with stimulated echoes (DENSE) is a quantitative MRI technique that encodes tissue displacement into the phase of the complex MRI images. Cine DENSE provides a time series of these images, thus facilitating the non-invasive study of myocardial kinematics. Epicardial and endocardial contours need to be defined at each frame on cine DENSE images for the quantification of regional displacement and strain as a function of time. This work presents a reliable and effective two-dimensional semi-automated segmentation technique that uses the encoded motion to project a manually-defined region of interest through time. Contours can then easily be extracted for each cardiac phase. This method boasts several advantages, including, (1) parameters are based on practical physiological limits, (2) contours are calculated for the first few cardiac phases, where it is difficult to visually distinguish blood from myocardium, and (3) the method is independent of the shape of the tissue delineated and can be applied to short- or long-axis views, and on arbitrary regions of interest. Motion-guided contours were compared to manual contours for six conventional and six slice-followed mid-ventricular short-axis cine DENSE datasets. Using an area measure of segmentation error, the accuracy of the segmentation algorithm was shown to be similar to inter-observer variability. In addition, a radial segmentation error metric was introduced for short-axis data. The average radial epicardial segmentation error was 0.36+/-0.08 and 0.40+/-0.10 pixels for slice-followed and conventional cine DENSE, respectively, and the average radial endocardial segmentation error was 0.46+/-0.12 and 0.46+/-0.16 pixels for slice following and conventional cine DENSE, respectively. Motion-guided segmentation employs the displacement-encoded phase shifts intrinsic to DENSE MRI to accurately propagate a single set of pre-defined contours throughout the remaining cardiac phases.

Increasing the signal-to-noise ratio in DENSE MRI by combining displacement-encoded echoes.

A new technique was developed to increase the signal-to-noise ratio (SNR) in displacement encoding with stimulated echoes (DENSE) MRI. This signal-averaged DENSE (sav-DENSE) technique is based on the SNR advantage of extracting a pair of DENSE images with uncorrelated noise from the complex complementary spatial modulation of the magnetization image, and combining them during image reconstruction. Eleven healthy volunteers were imaged at three short-axis locations with the use of sav-DENSE, cine DENSE, and myocardial tagging pulse sequences. In this study, sav-DENSE increased the SNR by 15-34% as compared to cine DENSE. Circumferential strain values measured by sav-DENSE and myocardial tagging were strongly correlated (slope = 0.95, intercept = -0.02, R = 0.92) and within the 95% limits of agreement. The breath-hold sav-DENSE technique yielded relatively accurate and precise quantification of 2D intramyocardial function, with a 40.2-ms temporal resolution and a 3.5 x 3.5 mm2 spatial resolution.