3D myocardial tissue tracking with slice followed cine DENSE MRI

Analysis of left ventricle regional myocardial motion for cardiac radioablation: Left ventricular motion analysis

PURPOSE

Left ventricle (LV) regional myocardial displacement due to cardiac motion was assessed using cardiovascular magnetic resonance (CMR) cine images to establish region-specific margins for cardiac radioablation treatments.

METHODS

CMR breath-hold cine images and LV myocardial tissue contour points were analyzed for 200 subjects, including controls (n = 50) and heart failure (HF) patients with preserved ejection fraction (HFpEF, n = 50), mid-range ejection fraction (HFmrEF, n = 50), and reduced ejection fraction (HFrEF, n = 50). Contour points were divided into segments according to the 17-segment model. For each patient, contour point displacements were determined for the long-axis (all 17 segments) and short-axis (segments 1-12) directions. Mean overall, tangential (longitudinal or circumferential), and normal (radial) displacements were calculated for the 17 segments and for each segment level.

RESULTS

The greatest overall motion was observed in the control group-long axis: 4.5 ± 1.2 mm (segment 13 [apical anterior] epicardium) to 13.8 ± 3.0 mm (segment 6 [basal anterolateral] endocardium), short axis: 4.3 ± 0.8 mm (segment 9 [mid inferoseptal] epicardium) to 11.5 ± 2.3 mm (segment 1 [basal anterior] endocardium). HF patients exhibited lesser motion, with the smallest overall displacements observed in the HFrEF group-long axis: 4.3 ± 1.7 mm (segment 13 [apical anterior] epicardium) to 10.6 ± 3.4 mm (segment 6 [basal anterolateral] endocardium), short axis: 3.9 ± 1.3 mm (segment 8 [mid anteroseptal] epicardium) to 7.4 ± 2.8 mm (segment 1 [basal anterior] endocardium).

CONCLUSIONS

This analysis provides an estimate of epicardial and endocardial displacement for the 17 segments of the LV for patients with normal and impaired LV function. This reference data can be used to establish treatment planning margin guidelines for cardiac radioablation. Smaller margins may be used for patients with higher degree of impaired heart function, depending on the LV segment.

A comprehensive comparison between shortest-path HARP refinement, SinMod, and DENSEanalysis processing tools applied to CSPAMM and DENSE images

We addressed comprehensively the performance of Shortest-Path HARP Refinement (SP-HR), SinMod, and DENSEanalysis using 2D slices of synthetic CSPAMM and DENSE images with realistic contrasts obtained from 3D phantoms. The three motion estimation techniques were interrogated under ideal and no-ideal conditions (with MR induced artifacts, noise, and through-plane motion), considering several resolutions and noise levels. Under noisy conditions, and for isotropic pixel sizes of 1.5 mm and 3.0 mm in CSPAMM and DENSE images respectively, the nRMSE obtained for the circumferential and radial strain components were 10.7 ± 10.8% and 25.5 ± 14.8% using SP-HR, 11.9 ± 2.5% and 29.3 ± 6.5% using SinMod, and 6.4 ± 2.0% and 18.2 ± 4.6% using DENSEanalysis. Overall, the results showed that SP-HR tends to fail for large tissue motions, whereas SinMod and DENSEanalysis gave accurate displacement and strain field estimations, being the last which performed the best.

3D Motion Estimation of Left Ventricular Dynamics Using MRI and Track-to-Track Fusion

Objective: This study investigates the estimation of three dimensional (3D) left ventricular (LV) motion using the fusion of different two dimensional (2D) cine magnetic resonance (CMR) sequences acquired during routine imaging sessions. Although standard clinical cine CMR data is inherently 2D, the actual underlying LV dynamics lies in 3D space and cannot be captured entirely using single 2D CMR image sequences. By utilizing the image information from various short-axis and long-axis image sequences, the proposed method intends to estimate the dynamic state vectors consisting of the position and velocity information of the myocardial borders in 3D space. Method: The proposed method comprises two main components: tracking myocardial points in 2D CMR sequences and fusion of multiple trajectories correspond to the tracked points. The tracking which yields the set of corresponding temporal points representing the myocardial points is performed using a diffeomorphic nonrigid image registration approach. The trajectories obtained from each cine CMR sequence is then fused with the corresponding trajectories from other CMR views using an unscented Kalman smoother (UKS) and a track-to-track fusion algorithm. Results: We evaluated the proposed method by comparing the results against CMR imaging with myocardial tagging. We report a quantitative performance analysis by projecting the state vector estimates we obtained onto 2D tagged CMR images acquired from the same subjects and comparing them against harmonic phase estimates. The proposed algorithm yielded a competitive performance with a mean root mean square error of 1.3±0.5 pixels (1.8±0.6 mm) evaluated over 118 image sequences acquired from 30 subjects. Conclusion: This study demonstrates that fusing the information from short and long-axis views of CMR improves the accuracy of cardiac tissue motion estimation. Clinical Impact: The proposed method demonstrates that the fusion of tissue tracking information from long and short-axis views improves the binary classification of the automated regional function assessment.

Myocardial motion analysis based on an optical flow method using tagged MR images

We developed a method of velocimetry based on an optical flow method using quantitative analyses of tagged magnetic resonance (MR) images (tagged MR-optical flow velocimetry, tMR-O velocimetry). The purpose of our study was to examine the accuracy of measurement of the proposed tMR-O velocimetry. We performed retrospective pseudo-electrocardiogram (ECG) gating tagged cine MR imaging on a rotating phantom. We optimized imaging parameters for tagged MR imaging, and validated the accuracy of tMR-O velocimetry. Our results indicated that the difference between the reference velocities and the computed velocities measured using optimal imaging parameters was less than 1%. In addition, we performed tMR-O velocimetry and echocardiography on 10 healthy volunteers, for four sections of the heart (apical, midventricular, and basal sections aligned with the short-axis, and a four-chamber section aligned with the long-axis), and obtained radial and longitudinal myocardial velocities in these sections. We compared the myocardial velocities obtained using tMR-O velocimetry with those obtained using echocardiography. Our results showed good agreement between tMR-O velocimetry and echocardiography in the radial myocardial velocities in three short-axial sections and longitudinal myocardial velocities on the midventricular portion of the four-chamber section in the long-axis. In the study conducted on the rotating phantom, tMR-O velocimetry showed high accuracy; moreover, in the healthy volunteers, the myocardial velocities obtained using tMR-O velocimetry were relatively similar to those obtained using echocardiography. In conclusion, tMR-O velocimetry is a potentially feasible method for analyzing myocardial motion in the human heart.

Segmental biventricular analysis of myocardial function using high temporal and spatial resolution tissue phase mapping

OBJECTIVE

Myocardial dysfunction of the right ventricle (RV) is an important indicator of RV diseases, e.g. RV infarction or pulmonary hypertension. Tissue phase mapping (TPM) has been widely used to determine function of the left ventricle (LV) by analyzing myocardial velocities. The analysis of RV motion is more complicated due to the different geometry and smaller wall thickness. The aim of this work was to adapt and optimize TPM to the demands of the RV.

MATERIALS AND METHODS

TPM measurements were acquired in 25 healthy volunteers using a velocity-encoded phase-contrast sequence and kt-accelerated parallel imaging in combination with optimized navigator strategy and blood saturation. Post processing was extended by a 10-segment RV model and a detailed biventricular analysis of myocardial velocities was performed.

RESULTS

High spatio-temporal resolution (1.0 × 1.0 × 6 mm³, 21.3 ms) and the optimized blood saturation enabled good delineation of the RV and its velocities. Global and segmental velocities, as well as time to peak velocities showed significant differences between the LV and RV. Furthermore, complex timing of the RV could be demonstrated by segmental time to peak analysis.

CONCLUSION

High spatio-temporal resolution TPM enables a detailed biventricular analysis of myocardial motion and might provide a reliable tool for description and detection of diseases affecting left and right ventricular function.

Assessment of reverse remodeling predicted by myocardial deformation on tissue tracking in patients with severe aortic stenosis: a cardiovascular magnetic resonance imaging study

BACKGROUND

The technique of tissue tracking with balanced steady-state free precession cine sequences was introduced, and allowed myocardial strain to be derived directly, offering advantages over traditional myocardial tagging. The aim of this study was to evaluate the correlation between reverse remodeling as an outcome and left ventricular strain using cardiovascular magnetic resonance imaging (CMR) tissue tracking, and to evaluate prediction of reverse remodeling by myocardial deformation in patients with severe aortic stenosis (AS).

METHODS

We enrolled 63 patients with severe AS and normal left ventricular (LV) systolic function (ejection fraction > 60%), who underwent both CMR and transthoracic echocardiography (Echo) before surgical aortic valve replacement (AVR). CMR at 1.5 T, including non and post-contrast T1 mapping for extracellular volume (ECV), was carried out to define the amount of myocardial fibrosis. Cardiac Performance Analysis software was used to derive myocardial deformation as strain parameters from three short-axis cine views (basal, mid and apical levels) and apical 2, 3, and 4 chamber views. The primary outcome was reverse remodeling, as evaluated by regression of left ventricular mass index (LVMI).

RESULTS

Median follow-up was 28.8 months (interquartile range 11.3-38.3 months). As evaluated by LVMI between baseline and follow-up, mass regression was significantly improved after AVR (baseline 145.9 ± 37.0 [g/m2] vs. follow-up 97.7 ± 22.2[g/m2], p < 0.001). Statistically significant Pearson's correlations with LVMI regression were observed for longitudinal global strain (r = -0.461, p < 0.001), radial strain (r = 0.391, p = 0.002), and circumferential strain (r = -0.334, p = 0.009). A simple linear regression analysis showed that all strain parameters could predict the amount of LVMI regression (P < 0.05), as well as non-contrast T1 value (beta = -0.314, p < 0.001) and ECV (beta = -2.546, p = 0.038). However, ECV had the lowest predictive power (multiple r2 = 0.071). Multiple regression analysis showed strain could independently predict the amount of LVMI regression and the longitudinal global strain (beta = -3.335, p < 0.001).

CONCLUSION

Longitudinal global strain measured by CMR tissue tracking as a technique was correlated with reverse remodeling as LVMI regression and was predictive of this outcome. As a simple and practical method, tissue tracking is promising to assess strain and predict reverse remodeling in severe AS, especially in patients with suboptimal Echo image quality.

Phase Vector Incompressible Registration Algorithm for Motion Estimation From Tagged Magnetic Resonance Images

Tagged magnetic resonance imaging has been used for decades to observe and quantify motion and strain of deforming tissue. It is challenging to obtain 3-D motion estimates due to a tradeoff between image slice density and acquisition time. Typically, interpolation methods are used either to combine 2-D motion extracted from sparse slice acquisitions into 3-D motion or to construct a dense volume from sparse acquisitions before image registration methods are applied. This paper proposes a new phase-based 3-D motion estimation technique that first computes harmonic phase volumes from interpolated tagged slices and then matches them using an image registration framework. The approach uses several concepts from diffeomorphic image registration with a key novelty that defines a symmetric similarity metric on harmonic phase volumes from multiple orientations. The material property of harmonic phase solves the aperture problem of optical flow and intensity-based methods and is robust to tag fading. A harmonic magnitude volume is used in enforcing incompressibility in the tissue regions. The estimated motion fields are dense, incompressible, diffeomorphic, and inverse-consistent at a 3-D voxel level. The method was evaluated using simulated phantoms, human brain data in mild head accelerations, human tongue data during speech, and an open cardiac data set. The method shows comparable accuracy to three existing methods while demonstrating low computation time and robustness to tag fading and noise.

Clinical feasibility and validation of 3D principal strain analysis from cine MRI: comparison to 2D strain by MRI and 3D speckle tracking echocardiography

Two-dimensional (2D) strain analysis is constrained by geometry-dependent reference directions of deformation (i.e. radial, circumferential, and longitudinal) following the assumption of cylindrical chamber architecture. Three-dimensional (3D) principal strain analysis may overcome such limitations by referencing intrinsic (i.e. principal) directions of deformation. This study aimed to demonstrate clinical feasibility of 3D principal strain analysis from routine 2D cine MRI with validation to strain from 2D tagged cine analysis and 3D speckle tracking echocardiography. Thirty-one patients undergoing cardiac MRI were studied. 3D strain was measured from routine, multi-planar 2D cine SSFP images using custom software designed to apply 4D deformation fields to 3D cardiac models to derive principal strain. Comparisons of strain estimates versus those by 2D tagged cine, 2D non-tagged cine (feature tracking), and 3D speckle tracking echocardiography (STE) were performed. Mean age was 51 ± 14 (36% female). Mean LV ejection fraction was 66 ± 10% (range 37–80%). 3D principal strain analysis was feasible in all subjects and showed high inter- and intra-observer reproducibility (ICC range 0.83–0.97 and 0.83–0.98, respectively—p < 0.001 for all directions). Strong correlations of minimum and maximum principal strain were respectively observed versus the following: 3D STE estimates of longitudinal (r = 0.81 and r = −0.64), circumferential (r = 0.76 and r = −0.58) and radial (r = −0.80 and r = 0.63) strain (p < 0.001 for all); 2D tagged cine estimates of longitudinal (r = 0.81 and r = −0.81), circumferential (r = 0.87 and r = −0.85), and radial (r = −0.76 and r = 0.81) strain (p < 0.0001 for all); and 2D cine (feature tracking) estimates of longitudinal (r = 0.85 and −0.83), circumferential (r = 0.88 and r = −0.87), and radial strain (r = −0.79 and r = 0.84, p < 0.0001 for all). 3D principal strain analysis is feasible using routine, multi-planar 2D cine MRI and shows high reproducibility with strong correlations to 2D conventional strain analysis and 3D STE-based analysis. Given its independence from geometry-related directions of deformation this technique may offer unique benefit for the detection and prognostication of myocardial disease, and warrants expanded investigation.

Myocardial strain in healthy adults across a broad age range as revealed by cardiac magnetic resonance imaging at 1.5 and 3.0T: Associations of myocardial strain with myocardial region, age, and sex

Purpose

To assess myocardial strain using cine displacement encoding with stimulated echoes (DENSE) using 1.5T and 3.0T MRI in healthy adults.

Materials and Methods

Healthy adults without any history of cardiovascular disease underwent magnetic resonance imaging (MRI) at 1.5T and 3.0T within 2 days. The MRI protocol included balanced steady‐state free‐precession (b‐SSFP), 2D cine‐echo planar imaging (EPI)‐DENSE, and late gadolinium enhancement in subjects >45 years. Acquisitions were divided into six segments; global and segmental peak longitudinal and circumferential strain were derived and analyzed by field strength, age, and gender.

Results

In all, 89 volunteers (mean age 44.8 ± 18.0 years, range: 18–87 years) underwent MRI at 1.5T, and 88 of these subjects underwent MRI at 3.0T (1.4 ± 1.4 days between the scans). Compared with 3.0T, the magnitudes of global circumferential (–19.5 ± 2.6% vs. –18.47 ± 2.6%; P = 0.001) and longitudinal (–12.47 ± 3.2% vs. –10.53 ± 3.1%; P = 0.004) strain were greater at 1.5T. At 1.5T, longitudinal strain was greater in females than in males: –10.17 ± 3.4% vs. –13.67 ± 2.4%; P = 0.001. Similar observations occurred for circumferential strain at 1.5T (–18.72 ± 2.2% vs. –20.10 ± 2.7%; P = 0.014) and at 3.0T (–17.92 ± 1.8% vs. –19.1 ± 3.1%; P = 0.047). At 1.5T, longitudinal and circumferential strain were not associated with age after accounting for sex (longitudinal strain P = 0.178, circumferential strain P = 0.733). At 3.0T, longitudinal and circumferential strain were associated with age (P < 0.05). Longitudinal strain values were greater in the apico‐septal, basal‐lateral, and mid‐lateral segments and circumferential strain in the inferior, infero‐lateral, and antero‐lateral LV segments.

Conclusion

Myocardial strain parameters as revealed by cine‐DENSE at different MRI field strengths were associated with myocardial region, age, and sex. J. Magn. Reson. Imaging 2016;44:1197–1205.

Clinical experience of strain imaging using DENSE for detecting infarcted cardiac segments

BACKGROUND

We hypothesised that myocardial deformation determined with magnetic resonance imaging (MRI) will detect myocardial scar.

METHODS

Displacement Encoding with Stimulated Echoes (DENSE) was used to calculate left ventricular strain in 125 patients (29 women and 96 men) with suspected coronary artery disease. The patients also underwent cine imaging and late gadolinium enhancement. 57 patients had a scar area >1% in at least one segment, 23 were considered free from coronary artery disease (control group) and 45 had pathological findings but no scar (mixed group). Peak strain was calculated in eight combinations: radial and circumferential strain in transmural, subendocardial and epicardial layers derived from short axis acquisition, and transmural longitudinal and radial strain derived from long axis acquisitions. In addition, the difference between strain in affected segments and reference segments, "differential strain", from the control group was analysed.

RESULTS

In receiver-operator-characteristic analysis for the detection of 50% transmurality, circumferential strain performed best with area-under-curve (AUC) of 0.94. Using a cut-off value of -17%, sensitivity was 95% at a specificity of 80%. AUC did not further improve with differential strain. There were significant differences between the control group and global strain circumferential direction (-17% versus -12%) and in the longitudinal direction (-13% versus -10%). Interobserver and scan-rescan reproducibility was high with an intraclass correlation coefficient (ICC) >0.93.

CONCLUSIONS

DENSE-derived circumferential strain may be used for the detection of myocardial segments with >50 % scar area. The repeatability of strain is satisfactory. DENSE-derived global strain agrees with other global measures of left ventricular ejection fraction.

Assessment of the relationships between myocardial contractility and infarct tissue revealed by serial magnetic resonance imaging in patients with acute myocardial infarction

Imaging changes in left ventricular (LV) volumes during the cardiac cycle and LV ejection fraction do not provide information on regional contractility. Displacement ENcoding with Stimulated Echoes (DENSE) is a strain-encoded cardiac magnetic resonance (CMR) technique that measures strain directly. We investigated the relationships between strain revealed by DENSE and the presence and extent of infarction in patients with recent myocardial infarction (MI). 50 male subjects were invited to undergo serial CMR within 7 days of MI (baseline) and after 6 months (follow-up; n = 47). DENSE and late gadolinium enhancement (LGE) images were acquired to enable localised regional quantification of peak circumferential strain (Ecc) and the extent of infarction, respectively. We assessed: (1) receiver operating characteristic (ROC) analysis for the classification of LGE, (2) strain differences according to LGE status (remote, adjacent, infarcted) and (3) changes in strain revealed between baseline and follow-up. 300 and 258 myocardial segments were available for analysis at baseline and follow-up respectively. LGE was present in 130/300 (43 %) and 97/258 (38 %) segments, respectively. ROC analysis revealed moderately high values for peak Ecc at baseline [threshold 12.8 %; area-under-curve (AUC) 0.88, sensitivity 84 %, specificity 78 %] and at follow-up (threshold 15.8 %; AUC 0.76, sensitivity 85 %, specificity 64 %). Differences were observed between remote, adjacent and infarcted segments. Between baseline and follow-up, increases in peak Ecc were observed in infarcted segments (median difference of 5.6 %) and in adjacent segments (1.5 %). Peak Ecc at baseline was indicative of the change in LGE status between baseline and follow-up. Strain-encoded CMR with DENSE has the potential to provide clinically useful information on contractility and its recovery over time in patients with MI.

Accurate high-resolution measurements of 3-D tissue dynamics with registration-enhanced displacement encoded MRI

Displacement fields are important to analyze deformation, which is associated with functional and material tissue properties often used as indicators of health. Magnetic resonance imaging (MRI) techniques like DENSE and image registration methods like Hyperelastic Warping have been used to produce pixel-level deformation fields that are desirable in high-resolution analysis. However, DENSE can be complicated by challenges associated with image phase unwrapping, in particular offset determination. On the other hand, Hyperelastic Warping can be hampered by low local image contrast. The current work proposes a novel approach for measuring tissue displacement with both DENSE and Hyperelastic Warping, incorporating physically accurate displacements obtained by the latter to improve phase characterization in DENSE. The validity of the proposed technique is demonstrated using numerical and physical phantoms, and in vivo small animal cardiac MRI.

Functional and molecular imaging techniques in aortic aneurysm disease

PURPOSE OF REVIEW

Functional and molecular aortic imaging has shown great promise for evaluation of aortic disease, and may soon augment conventional assessment of aortic dimensions for the clinical management of patients.

RECENT FINDINGS

A range of imaging techniques is available for evaluation of patients with aortic disease. Magnetic resonance blood flow imaging can identify atherosclerosis prone aortic regions and may be useful for predicting aneurysm growth. Computational modeling can demonstrate significant differences in wall stress between abdominal aortic aneurysms of similar size and may better predict rupture than diameter alone. Metabolic imaging with fluorodeoxyglucose-PET [(FDG)-PET] can identify focal aortic wall inflammation that may portend rapid progression of disease. Molecular imaging with probes that target collagen and elastin can directly exhibit changes in the vessel wall associated with disease.

SUMMARY

The complexity of aortic disease is more fully revealed with new functional imaging techniques than with conventional anatomic analysis alone. This may better inform surveillance imaging regimens, medical management and decisions regarding early intervention for aortic disease.

Three-dimensional regional strain computation method with displacement encoding with stimulated echoes (DENSE) in non-ischemic, non-valvular dilated cardiomyopathy patients and healthy subjects validated by tagged MRI

PURPOSE

Fast cine displacement encoding with stimulated echoes (DENSE) MR has higher spatial resolution and enables rapid postprocessing. Thus we compared the accuracy of regional strains computation by DENSE with tagged MR in healthy and non-ischemic, non-valvular dilated cardiomyopathy (DCM) subjects.

MATERIALS AND METHODS

Validation of three-dimensional regional strains computed with DENSE was conducted in reference to standard tagged MRI (TMRI) in healthy subjects and patients with DCM. Additional repeatability studies in healthy subjects were conducted to increase confidence in DENSE. A meshfree multiquadrics radial point interpolation method (RPIM) was used for computing Lagrange strains in sixteen left ventricular segments. Bland-Altman analysis and Student's t-tests were conducted to observe similarities in regional strains between sequences and in DENSE repeatability studies.

RESULTS

Regional circumferential strains ranged from -0.21 ± 0.07 (Lateral-Apex) to -0.11 ± 0.05 (Posterorseptal-Base) in healthy subjects and -0.15 ± 0.04 (Anterior-Apex) to -0.02 ± 0.08 (Posterorseptal-Base) in DCM patients. Computed mean differences in regional circumferential strain from the DENSE-TMRI comparison study was 0.01 ± 0.03 (95% limits of agreement) in normal subjects, -0.01 ± 0.06 in DCM patients and 0.0 ± 0.02 in repeatability studies, with similar agreements in longitudinal and radial strains.

CONCLUSION

We found agreement between DENSE and tagged MR in patients and volunteers in terms of evaluation of regional strains.

Regional heart motion abnormality detection via multiview fusion

This study investigates regional heart motion abnormality detection via multiview fusion in cine cardiac MR images. In contrast to previous methods which rely only on short-axis image sequences, the proposed approach exploits the information from several other long-axis image sequences, namely, 2-chamber, 3-chamber and 4-chamber MR images. Our analysis follows the standard issued by American Heart Association to identify 17 standardized left ventricular segments. The proposed method first computes an initial sequence of corresponding myocardial points using a nonrigid image registration algorithm within each sequence. Then, these points were mapped to 3D space and tracked using Unscented Kalman Filter (UKS). We propose a maximum likelihood based track-to-track fusion approach to combine UKS tracks from multiple image views. Finally, we use a Shannon's differential entropy of distributions of potential classifiers obtained from multiview fusion estimates, and a naive Bayes classifier algorithm to automatically detect abnormal functional regions of the myocardium. We proved the benefits of the proposed method by comparing the classification results with and without fusion over 480 regional myocardial segments obtained from 30 subjects. The evaluations in comparisons to the ground truth classifications by radiologists showed that the proposed fusion yielded an area-under-the-curve (AUC) of 95.9%, bringing a significant improvement of 3.8% in comparisons to previous methods that use only short-axis images.

Undersampled cine 3D tagging for rapid assessment of cardiac motion

BACKGROUND

CMR allows investigating cardiac contraction, rotation and torsion non-invasively by the use of tagging sequences. Three-dimensional tagging has been proposed to cover the whole-heart but data acquisition requires three consecutive breath holds and hence demands considerable patient cooperation. In this study we have implemented and studied k-t undersampled cine 3D tagging in conjunction with k-t PCA reconstruction to potentially permit for single breath-hold acquisitions.

METHODS

The performance of undersampled cine 3D tagging was investigated using computer simulations and in-vivo measurements in 8 healthy subjects and 5 patients with myocardial infarction. Fully sampled data was obtained and compared to retrospectively and prospectively undersampled acquisitions. Fully sampled data was acquired in three consecutive breath holds. Prospectively undersampled data was obtained within a single breath hold. Based on harmonic phase (HARP) analysis, circumferential shortening, rotation and torsion were compared between fully sampled and undersampled data using Bland-Altman and linear regression analysis.

RESULTS

In computer simulations, the error for circumferential shortening was 2.8 ± 2.3% and 2.7 ± 2.1% for undersampling rates of R = 3 and 4 respectively. Errors in ventricular rotation were 2.5 ± 1.9% and 3.0 ± 2.2% for R = 3 and 4. Comparison of results from fully sampled in-vivo data acquired with prospectively undersampled acquisitions showed a mean difference in circumferential shortening of -0.14 ± 5.18% and 0.71 ± 6.16% for R = 3 and 4. The mean differences in rotation were 0.44 ± 1.8° and 0.73 ± 1.67° for R = 3 and 4, respectively. In patients peak, circumferential shortening was significantly reduced (p < 0.002 for all patients) in regions with late gadolinium enhancement.

CONCLUSION

Undersampled cine 3D tagging enables significant reduction in scan time of whole-heart tagging and facilitates quantification of shortening, rotation and torsion of the left ventricle without adding significant errors compared to previous 3D tagging approaches.

MR assessment of regional myocardial mechanics

Regional myocardial function can be measured by several MR techniques including tissue tagging, phase velocity mapping, and more recently, displacement encoding with stimulated echoes (DENSE) and strain encoding (SENC). Each of these techniques was developed separately and has undergone significant change since its original implementation. As a result, in the current literature, the common features and the differences between the techniques and what they measure are often unclear and confusing. This review article delivers an extensively referenced introductory text which clarifies the current methodology from the starting point of the Bloch equations. By doing this in a consistent way for each method, the similarities and differences between them are highlighted. In addition, their capabilities and limitations are discussed, together with their relative advantages and disadvantages. While the focus is on sequence design and development, the principal parameters measured by each technique are also summarized, together with brief results, with the reader being directed to the extensive literature on data processing and clinical applications for more detail.

Temporal 3D Lagrangian strain from 2D slice-followed cine DENSE MRI

A quantitative analysis of myocardial mechanics is fundamental to the understanding of cardiac function, diagnosis of heart disease and assessment of therapeutic intervention. In clinical practice, most cardiac imaging analyses are performed in 2D because of the limited scan time available. However, the obtained information from a 2D measurement is limited. This study presents a method to obtain temporal evolutions of transmural 3D Lagrangian strains from two intersecting 2D planes of slice-followed cine displacement encoding with stimulated echoes (DENSE) data using a bilinear-cubic polynomial element to resolve strain from the displaced myocardial positions. The method was validated against an analytical standard and has been applied to in vivo data acquired on a 3 T magnetic resonance system from a healthy volunteer to quantify systolic strains at the anterior-basal region of left ventricular myocardium. The method demonstrates accurate results when validated in the analytical model, and the in vivo results agree within experimental accuracy with values reported in the literature. Even with a short scan time, this method provides the full 3D Lagrangian strain tensor from two 2D DENSE measurements.

Cardiac motion and deformation recovery from MRI: a review

Magnetic resonance imaging (MRI) is a highly advanced and sophisticated imaging modality for cardiac motion tracking and analysis, capable of providing 3D analysis of global and regional cardiac function with great accuracy and reproducibility. In the past few years, numerous efforts have been devoted to cardiac motion recovery and deformation analysis from MR image sequences. Many approaches have been proposed for tracking cardiac motion and for computing deformation parameters and mechanical properties of the heart from a variety of cardiac MR imaging techniques. In this paper, an updated and critical review of cardiac motion tracking methods including major references and those proposed in the past ten years is provided. The MR imaging and analysis techniques surveyed are based on cine MRI, tagged MRI, phase contrast MRI, DENSE, and SENC. This paper can serve as a tutorial for new researchers entering the field.

Incompressible deformation estimation algorithm (IDEA) from tagged MR images

Measuring the 3D motion of muscular tissues, e.g., the heart or the tongue, using magnetic resonance (MR) tagging is typically carried out by interpolating the 2D motion information measured on orthogonal stacks of images. The incompressibility of muscle tissue is an important constraint on the reconstructed motion field and can significantly help to counter the sparsity and incompleteness of the available motion information. Previous methods utilizing this fact produced incompressible motions with limited accuracy. In this paper, we present an incompressible deformation estimation algorithm (IDEA) that reconstructs a dense representation of the 3D displacement field from tagged MR images and the estimated motion field is incompressible to high precision. At each imaged time frame, the tagged images are first processed to determine components of the displacement vector at each pixel relative to the reference time. IDEA then applies a smoothing, divergence-free, vector spline to interpolate velocity fields at intermediate discrete times such that the collection of velocity fields integrate over time to match the observed displacement components. Through this process, IDEA yields a dense estimate of a 3D displacement field that matches our observations and also corresponds to an incompressible motion. The method was validated with both numerical simulation and in vivo human experiments on the heart and the tongue.

Myocardial tagging by cardiovascular magnetic resonance: evolution of techniques--pulse sequences, analysis algorithms, and applications

Cardiovascular magnetic resonance (CMR) tagging has been established as an essential technique for measuring regional myocardial function. It allows quantification of local intramyocardial motion measures, e.g. strain and strain rate. The invention of CMR tagging came in the late eighties, where the technique allowed for the first time for visualizing transmural myocardial movement without having to implant physical markers. This new idea opened the door for a series of developments and improvements that continue up to the present time. Different tagging techniques are currently available that are more extensive, improved, and sophisticated than they were twenty years ago. Each of these techniques has different versions for improved resolution, signal-to-noise ratio (SNR), scan time, anatomical coverage, three-dimensional capability, and image quality. The tagging techniques covered in this article can be broadly divided into two main categories: 1) Basic techniques, which include magnetization saturation, spatial modulation of magnetization (SPAMM), delay alternating with nutations for tailored excitation (DANTE), and complementary SPAMM (CSPAMM); and 2) Advanced techniques, which include harmonic phase (HARP), displacement encoding with stimulated echoes (DENSE), and strain encoding (SENC). Although most of these techniques were developed by separate groups and evolved from different backgrounds, they are in fact closely related to each other, and they can be interpreted from more than one perspective. Some of these techniques even followed parallel paths of developments, as illustrated in the article. As each technique has its own advantages, some efforts have been made to combine different techniques together for improved image quality or composite information acquisition. In this review, different developments in pulse sequences and related image processing techniques are described along with the necessities that led to their invention, which makes this article easy to read and the covered techniques easy to follow. Major studies that applied CMR tagging for studying myocardial mechanics are also summarized. Finally, the current article includes a plethora of ideas and techniques with over 300 references that motivate the reader to think about the future of CMR tagging.

Phase II trials in heart failure: the role of cardiovascular imaging

The development of new therapies for heart failure (HF), especially acute HF, has proven to be quite challenging; and therapies evaluated in HF have greatly outnumbered treatments that are eventually successful in obtaining regulatory approval. Thus, the development of therapies for HF remains a vexing problem for pharmaceutical and device companies, clinical trialists, and health care professionals. Nowhere is this more apparent than in the phase II HF clinical trial, in which the goal is to determine whether an investigational agent should move forward to a phase III trial. Recent advancements in noninvasive cardiovascular imaging have allowed a new era of comprehensive phenotyping of cardiac structure and function in phase II HF trials. Besides using imaging parameters to predict success of subsequent phase III outcome studies, it is essential to also use imaging in phase II HF trials in a way that increases understanding of drug or device mechanism. Determination of the patients who would benefit most from a particular drug or device could decrease heterogeneity of phase III trial participants and lead to more successful HF clinical trials. In this review, we outline advantages and disadvantages of imaging various aspects of cardiac structure and function that are potential targets for therapy in HF, compare and contrast imaging modalities, provide practical advice for the use of cardiovascular imaging in drug development, and conclude with some novel uses of cardiac imaging in phase II HF trials.

Emerging MRI methods in translational cardiovascular research

Cardiac magnetic resonance imaging (CMR) has become a reference standard modality for imaging of left ventricular (LV) structure and function and, using late gadolinium enhancement, for imaging myocardial infarction. Emerging CMR techniques enable a more comprehensive examination of the heart, making CMR an excellent tool for use in translational cardiovascular research. Specifically, emerging CMR methods have been developed to measure the extent of myocardial edema, changes in ventricular mechanics, changes in tissue composition as a result of fibrosis, and changes in myocardial perfusion as a function of both disease and infarct healing. New CMR techniques also enable the tracking of labeled cells, molecular imaging of biomarkers of disease, and changes in calcium flux in cardiomyocytes. In addition, MRI can quantify blood flow velocity and wall shear stress in large blood vessels. Almost all of these techniques can be applied in both pre-clinical and clinical settings, enabling both the techniques themselves and the knowledge gained using such techniques in pre-clinical research to be translated from the lab bench to the patient bedside.

Acceleration of tissue phase mapping by k-t BLAST: a detailed analysis of the influence of k-t-BLAST for the quantification of myocardial motion at 3T

BACKGROUND

The assessment of myocardial motion with tissue phase mapping (TPM) provides high spatiotemporal resolution and quantitative motion information in three directions. Today, whole volume coverage of the heart by TPM encoding at high spatial and temporal resolution is limited by long data acquisition times. Therefore, a significant increase in imaging speed without deterioration of the quantitative motion information is required. For this purpose, the k-t BLAST acceleration technique was combined with TPM black-blood functional imaging of the heart. Different k-t factors were evaluated with respect to their impact on the quantitative assessment of cardiac motion.

RESULTS

It is demonstrated that a k-t BLAST factor of two can be used with a marginal, but statistically significant deterioration of the quantitative motion data. Further increasing the k-t acceleration causes substantial alteration of the peak velocities and the motion pattern, but the temporal behavior of the contraction is well maintained up to an acceleration factor of six.

CONCLUSIONS

The application of k-t BLAST for the acceleration of TPM appears feasible. A reduction of the acquisition time of almost 45% could be achieved without substantial loss of quantitative motion information.

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.

Strain and torsion quantification in mouse hearts under dobutamine stimulation using 2D multiphase MR DENSE

In this study, a 2D multiphase magnetic resonance displacement encoding with stimulated echoes (DENSE) imaging and analysis method was developed for direct quantification of Lagrangian strain in the mouse heart. Using the proposed method, <10 ms temporal resolution and 0.56 mm in-plane resolution were achieved. A validation study that compared strain calculation by displacement encoding with stimulated echoes and by magnetic resonance tagging showed high correlation between the two methods (R(2) > 0.80). Regional ventricular wall strain and twist were characterized in mouse hearts at baseline and under dobutamine stimulation. Dobutamine stimulation induced significant increase in radial and circumferential strains and torsion at peak systole. A rapid untwisting was also observed during early diastole. This work demonstrates the capability of characterizing cardiac functional response to dobutamine stimulation in the mouse heart using 2D multiphase magnetic resonance displacement encoding with stimulated echoes.

Single-breath-hold multiple-slice DENSE MRI

A method to acquire multiple displacement encoded slices within a single breath hold is presented. Efficiency is improved over conventional DENSE without compromising image quality by readout of multiple slices in the same cardiac cycle, thus utilizing the position-encoded stimulated echo available in the whole heart. The method was evaluated by comparing strain values obtained using the proposed method to strain values obtained by conventional separate breath-hold single-slice DENSE acquisitions. Good agreement (Lagrangian E(2) strain bias = 0.000, 95% limits of agreement +/- 0.04, root-mean-square-difference 0.02 [9.4% of mean end-systolic E(2)]) was found between the methods, indicating that the proposed method can replace a multiple breath-hold acquisition. Eliminating the need for multiple breath holds reduces the risk of changes in breath-hold positions or heart rate, results in higher patient comfort, and facilitates inclusion of DENSE in a clinical routine protocol.

Analytical method to measure three-dimensional strain patterns in the left ventricle from single slice displacement data

BACKGROUND

Displacement encoded Cardiovascular MR (CMR) can provide high spatial resolution measurements of three-dimensional (3D) Lagrangian displacement. Spatial gradients of the Lagrangian displacement field are used to measure regional myocardial strain. In general, adjacent parallel slices are needed in order to calculate the spatial gradient in the through-slice direction. This necessitates the acquisition of additional data and prolongs the scan time. The goal of this study is to define an analytic solution that supports the reconstruction of the out-of-plane components of the Lagrangian strain tensor in addition to the in-plane components from a single-slice displacement CMR dataset with high spatio-temporal resolution. The technique assumes incompressibility of the myocardium as a physical constraint.

RESULTS

The feasibility of the method is demonstrated in a healthy human subject and the results are compared to those of other studies. The proposed method was validated with simulated data and strain estimates from experimentally measured DENSE data, which were compared to the strain calculation from a conventional two-slice acquisition.

CONCLUSION

This analytical method reduces the need to acquire data from adjacent slices when calculating regional Lagrangian strains and can effectively reduce the long scan time by a factor of two.

Deformation analysis of 3D tagged cardiac images using an optical flow method

BACKGROUND

This study proposes and validates a method of measuring 3D strain in myocardium using a 3D Cardiovascular Magnetic Resonance (CMR) tissue-tagging sequence and a 3D optical flow method (OFM).

METHODS

Initially, a 3D tag MR sequence was developed and the parameters of the sequence and 3D OFM were optimized using phantom images with simulated deformation. This method then was validated in-vivo and utilized to quantify normal sheep left ventricular functions.

RESULTS

Optimizing imaging and OFM parameters in the phantom study produced sub-pixel root-mean square error (RMS) between the estimated and known displacements in the x (RMSx = 0.62 pixels (0.43 mm)), y (RMSy = 0.64 pixels (0.45 mm)) and z (RMSz = 0.68 pixels (1 mm)) direction, respectively. In-vivo validation demonstrated excellent correlation between the displacement measured by manually tracking tag intersections and that generated by 3D OFM (R >or= 0.98). Technique performance was maintained even with 20% Gaussian noise added to the phantom images. Furthermore, 3D tracking of 3D cardiac motions resulted in a 51% decrease in in-plane tracking error as compared to 2D tracking. The in-vivo function studies showed that maximum wall thickening was greatest in the lateral wall, and increased from both apex and base towards the mid-ventricular region. Regional deformation patterns are in agreement with previous studies on LV function.

CONCLUSION

A novel method was developed to measure 3D LV wall deformation rapidly with high in-plane and through-plane resolution from one 3D cine acquisition.

TRUHARP: SINGLE BREATH-HOLD MRI FOR HIGH RESOLUTION CARDIAC MOTION AND STRAIN QUANTIFICATION

MRI techniques for tissue motion and strain quantifications have limited resolution because of interference from the conjugate echo or spectral peak in Fourier space. Methods have been proposed to remove this interference in order to improve resolution; however, these methods are clinically impractical due to long image acquisition time. In this paper, we propose TruHARP, an MRI motion and strain quantification methodology that involves a novel single breath-hold imaging protocol. In post-processing, TruHARP separates the spectral peaks in the acquired datasets, enabling high resolution motion and strain quantification. The impact of high resolution on circumferential and radial strain is studied using a realistic simulation and the improvement in strain maps is demonstrated in an in-vivo human study.

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.