Spline-based cardiac motion tracking using velocity-encoded magnetic resonance imaging

Validation of cardiac magnetic-resonance-derived left ventricular strain measurements from free-breathing motion-corrected cine imaging

BACKGROUND

Myocardial strain is an important measure of cardiac function and can be assessed on cardiac magnetic resonance (MR) through the current gold standard of breath-held segmented steady-state free precession (SSFP) cine imaging. Novel free-breathing techniques have been validated for volumetry and systolic function, allowing for evaluation of sicker and younger children who cannot reliably hold their breath. It is unclear whether strain measurements can be reliably performed on free-breathing, motion-corrected, re-binning cine images.

OBJECTIVE

To compare strain analysis from motion-corrected retrospective re-binning images to the breath-held SSFP cine images to explore their validity.

MATERIALS AND METHODS

Twenty-five children and young adults, ages (2.1-18.6 years) underwent breath-held and motion-corrected retrospective re-binning cine techniques during the same MR examination on a 1.5-tesla magnet. We measured endocardial end-systolic global circumferential strain and endocardial averaged segmental strain using commercial software (MEDIS QStrain 2.1). We used Pearson correlation coefficients to test agreement across techniques.

RESULTS

Analysis was possible in all 25 breath-held and motion-corrected retrospective re-binning studies. Global circumferential strain and endocardial averaged segmental strain obtained by motion-corrected retrospective re-binning compared favorably to breath-held studies. Global circumferential strain linear regression models demonstrated acceptable agreement, with coefficients of determination of 0.75 for breath-held compared to motion-corrected retrospective re-binning (P<0.001) and for endocardial averaged segmental strain comparisons yielded 0.77 for breath-held vs. motion-corrected retrospective re-binning (P<0.001). Bland-Altman assessment demonstrated minimal bias for breath-held compared to motion-corrected retrospective re-binning (mean 2.4 and 1.9, respectively, for global circumferential strain and endocardial averaged segmental strain).

CONCLUSION

Free-breathing imaging by motion-corrected retrospective re-binning cine imaging provides adequate spatial and temporal resolution to measure myocardial deformation when compared to the gold-standard breath-held SSFP cine imaging in children with normal or borderline systolic function.

Spatiotemporal Strategies for Joint Segmentation and Motion Tracking From Cardiac Image Sequences

Although accurate and robust estimations of the deforming cardiac geometry and kinematics from cine tomographic medical image sequences remain a technical challenge, they have significant clinical value. Traditionally, boundary or volumetric segmentation and motion estimation problems are considered as two sequential steps, even though the order of these processes can be different. In this paper, we present an integrated, spatiotemporal strategy for the simultaneous joint recovery of these two ill-posed problems. We use a mesh-free Galerkin formulation as the representation and computation platform, and adopt iterative procedures to solve the governing equations. Specifically, for each nodal point, the external driving forces are individually constructed through the integration of data-driven edginess measures, prior spatial distributions of myocardial tissues, temporal coherence of image-derived salient features, imaging/image-derived Eulerian velocity information, and cyclic motion model of myocardial behavior. The proposed strategy is accurate and very promising application results are shown from synthetic data, magnetic resonance (MR) phase contrast, tagging image sequences, and gradient echo cine MR image sequences.

Novel insight into the detailed myocardial motion and deformation of the rodent heart using high-resolution phase contrast cardiovascular magnetic resonance

BACKGROUND

Phase contrast velocimetry cardiovascular magnetic resonance (PC-CMR) is a powerful and versatile tool allowing assessment of in vivo motion of the myocardium. However, PC-CMR is sensitive to motion related artifacts causing errors that are geometrically systematic, rendering regional analysis of myocardial function challenging. The objective of this study was to establish an optimized PC-CMR method able to provide novel insight in the complex regional motion and strain of the rodent myocardium, and provide a proof-of-concept in normal and diseased rat hearts with higher temporal and spatial resolution than previously reported.

METHODS

A PC-CMR protocol optimized for assessing the motion and deformation of the myocardium in rats with high spatiotemporal resolution was established, and ten animals with different degree of cardiac dysfunction underwent examination and served as proof-of-concept. Global and regional myocardial velocities and circumferential strain were calculated, and the results were compared to five control animals. Furthermore, the global strain measurements were validated against speckle-tracking echocardiography, and inter- and intrastudy variability of the protocol were evaluated.

RESULTS

The presented method allows assessment of regional myocardial function in rats with high level of detail; temporal resolution was 3.2 ms, and analysis was done using 32 circumferential segments. In the dysfunctional hearts, global and regional function were distinctly altered, including reduced global peak values, increased regional heterogeneity and increased index of dyssynchrony. Strain derived from the PC-CMR data was in excellent agreement with echocardiography (r = 0.95, p < 0.001; limits-of-agreement -0.02 ± 3.92%strain), and intra- and interstudy variability were low for both velocity and strain (limits-of-agreement, radial motion: 0.01 ± 0.32 cm/s and -0.06 ± 0.75 cm/s; circumferential strain: -0.16 ± 0.89%strain and -0.71 ± 1.67%strain, for intra- and interstudy, respectively).

CONCLUSION

We demonstrate, for the first time, that PC-CMR enables high-resolution evaluation of in vivo circumferential strain in addition to myocardial motion of the rat heart. In combination with the superior geometric robustness of CMR, this ultimately provides a tool for longitudinal studies of regional function in rodents with high level of detail.

Assessing myocardial recovery following ST-segment elevation myocardial infarction: short- and long-term perspectives using cardiovascular magnetic resonance

Myocardial recovery after revascularization for ST-segment elevation myocardial infarction (STEMI) remains a significant diagnostic and, despite novel treatment strategies, a therapeutic challenge. Cardiovascular magnetic resonance (CMR) has emerged as a valuable clinical and research tool after acute STEMI. It represents the gold standard for functional and morphological evaluation of the left ventricle. Gadolinium-based perfusion and late-enhancement viability imaging has expanded our knowledge about the underlying pathologies of inadequate myocardial recovery. T2-weighted imaging of myocardial salvage after early reperfusion of the infarct-related artery underlines the effectiveness of current invasive treatment for STEMI. In the last decade, the number of publications on CMR after acute STEMI continued to rise, with no plateau in sight. Currently, CMR research is gathering robust prognostic data on standardized CMR protocols with the aim to substantially improve patient care and prognosis. Beyond established CMR protocols, more specific methods such as magnetic resonance relaxometry, myocardial tagging, 4D phase-contrast imaging and novel superparamagnetic contrast agents are emerging. This review will discuss the currently available data on the use of CMR after acute STEMI and take a brief look at developing new methods currently under investigation.

A magnetic resonance software simulator for the evaluation of myocardial deformation estimation

This paper proposes a methodology to design a physiologically realistic computer simulator of images of the left ventricle myocardium based on a patient-specific biomechanical model. The simulator takes a magnetic resonance image of a given patient at end diastole, uses a manual segmentation of that image to model the geometry of the myocardium and sets the parameters of the constitutive model used for biomechanical simulation according to a regional labeling of the contractility of the myocardium for that patient. The simulated deformations are used to warp the magnetic resonance dataset throughout the cardiac cycle to generate different image modalities. The simulator is validated by quantifying its ability to model actual deformations in a set of patients affected by an acute myocardial infarction. Specifically a high correlation has been encountered between the ejection fraction derived from the simulated end systolic deformation of the myocardium and the myocardium segmented from actual data. Additionally, most of the parameters that describe the simulated deformation compare well with reported values. Overall, the simulator is intended as a testbed for extensive comparisons of myocardial motion tracking methods due to its ability to relate the impaired myocardial function with the associated ventricular remodeling, a novel contribution in the literature of cardiac image simulators.

Regional heart motion abnormality detection: an information theoretic approach

Tracking regional heart motion and detecting the corresponding abnormalities play an essential role in the diagnosis of cardiovascular diseases. Based on functional images, which are subject to noise and segmentation/registration inaccuracies, regional heart motion analysis is acknowledged as a difficult problem and, therefore, incorporation of prior knowledge is desirable to enhance accuracy. Given noisy data and a nonlinear dynamic model to describe myocardial motion, an unscented Kalman smoother is proposed in this study to estimate the myocardial points. Due to the similarity between the statistical information of normal and abnormal heart motions, detecting and classifying abnormality is a challenging problem. We use the Shannon's differential entropy of the distributions of potential classifier features to detect and locate regional heart motion abnormality. A naive Bayes classifier algorithm is constructed from the Shannon's differential entropy of different features to automatically detect abnormal functional regions of the myocardium. Using 174 segmented short-axis magnetic resonance cines obtained from 58 subjects (21 normal and 37 abnormal), the proposed method is quantitatively evaluated by comparison with ground truth classifications by radiologists over 928 myocardial segments. The proposed method performed significantly better than other recent methods, and yielded an accuracy of 86.5% (base), 89.4% (mid-cavity) and 84.5% (apex). The overall classification accuracy was 87.1%. Furthermore, standard kappa statistic comparisons between the proposed method and visual wall motion scoring by radiologists showed that the proposed algorithm can yield a kappa measure of 0.73.

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.

Coronary microembolization causes long-term detrimental effects on regional left ventricular function

OBJECTIVES

To investigate whether coronary microemboli have long-term effects on left ventricular (LV) function in an experimental model. Furthermore, to determine if first-pass perfusion and late gadolinium enhancement (LGE) patterns differs between small- and large-sized microemboli.

DESIGN

Six pigs underwent left anterior descending (LAD)-coronary microembolization with small-sized (40-120 μm, n ∼ 250 000) microemboli using a combined x-ray and MRI-system. MR-images before, one hour after and 7-8 weeks after microembolization were obtained. Results were compared to MRI obtained by large-sized (100-300 μm, n ∼ 7200) microemboli.

RESULTS

Cine MRI showed an acute drop in ejection fraction (from 49.5 ± 2.6% to 32.5 ± 2.8) that substantially recovered at 7-8 weeks (47.5 ± 3.2%). Regional LV-function assessed as circumferential, longitudinal and radial strain declined in both microinfarcts and remote regions followed by partial recovery at 7-8 weeks. The decline in LV function and the severe perfusion deficit from the small microemboli was similar to the large microemboli at one hour. There was a significant recovery in perfusion at 7-8 weeks in the microinfarcts. LGE demonstrated the microinfarcts at 7-8 weeks but not at one hour and the microinfarcts were confirmed by histopathology.

CONCLUSION

Microembolization causes long-term, regional LV dysfunction and this study confirmed the need of a comprehensive MRI-protocol for the detection of microinfarcts. These findings suggest that even small microemboli (40-120 μm in diameter), which may escape the distal protective devices influence cardiac function.

Effects of gadolinium contrast agent on aortic blood flow and myocardial strain measurements by phase-contrast cardiovascular magnetic resonance

BACKGROUND

Quantitative blood flow and aspects of regional myocardial function such as myocardial displacement and strain can be measured using phase-contrast cardiovascular magnetic resonance (PC-CMR). Since a gadolinium-based contrast agent is often used to measure myocardial infarct size, we sought to determine whether the contrast agent affects measurements of aortic flow and myocardial displacement and strain. Phase-contrast data pre and post contrast agent was acquired during free breathing using 1.5T PC-CMR.

RESULTS

For aortic flow and regional myocardial function 12 and 17 patients were analysed, respectively. The difference pre and post contrast agent was 0.03±0.16 l/min for cardiac output, and 0.1±0.5 mm for myocardial displacement. Linear regression for myocardial displacement (MD) after and before contrast agent (CA) showed MDpostCA=0.95MDpreCA+0.05 (r=0.95, p<0.001). For regional myocardial function, the contrast-to-noise ratios for left ventricular myocardial wall versus left ventricular lumen were pre and post contrast agent administration 7.4±3.3 and 4.4±8.9, respectively (p<0.001). The contrast-to-noise ratios for left ventricular myocardial wall versus surrounding tissue were pre and post contrast agent administration -16.9±22 and -0.2±6.3, respectively (p<0.0001).

CONCLUSIONS

Quantitative measurements of aortic flow yield equal results both in the absence and presence of gadolinium contrast agent. The total examination time may thereby be reduced when assessing both viability and quantitative flow using PC-CMR, by assessing aortic flow post contrast agent administration. Phase-contrast information for myocardial displacement is also assessable both in the absence and presence of contrast agent. However, delineation of the myocardium may be difficult or impossible post contrast agent due to the lower image contrast. Acquisition of myocardial displacement should therefore be performed pre contrast agent using current PC-CMR sequences.

Motion estimation of tagged cardiac magnetic resonance images using variational techniques

This work presents a new method for motion estimation of tagged cardiac magnetic resonance sequences based on variational techniques. The variational method has been improved by adding a new term in the optical flow equation that incorporates tracking points with high stability of phase. Results were obtained through simulated and real data, and were validated by manual tracking and with respect to a reference state-of-the-art method: harmonic phase imaging (HARP). The error, measured in pixels per frame, obtained with the proposed variational method is one order of magnitude smaller than the one achieved by the reference method, and it requires a lower computational cost.

Persistent decline in longitudinal and radial strain after coronary microembolization detected on velocity encoded phase contrast magnetic resonance imaging

PURPOSE

To use velocity-encoded phase contrast (PC) MRI in assessing the effect of coronary microembolization on longitudinal and radial myocardial strain.

MATERIALS AND METHODS

A combined X-ray and MR system (XMR) was used for selective left anterior descending artery catheterization and microinfarct assessment in swine (n = 6). The embolized area at risk was defined on perfusion MRI followed by administration of a 7500 count (size = 100-300 microm) of the embolic agent. Quantification of strain and microinfarction was performed at 1 h and 1 week using PC-MRI and delayed enhancement (DE) MRI, respectively. At postmortem, sliced hearts were stained to define microinfarction.

RESULTS

Baseline longitudinal and radial strain did not differ between area-at-risk and remote myocardium. The embolized territory (area at risk) showed significant decline in longitudinal strain from -11.5 +/- 3.2% to 1.8 +/- 2.5% at 1 h (P < 0.05) and -3.9 +/- 1.1% at 1 week (P < 0.05). Similarly, regional radial strain progressively declined from 23.6 +/- 2.5% at baseline to 12.5 +/- 3.7% at 1 h (P < 0.05) and 4.8 +/- 5.0% at 1 week (P < 0.01). The size of microinfarction was not significantly different between DE-MRI and histochemical staining.

CONCLUSION

PC-MRI is sensitive in assessing changes in regional longitudinal and radial strain after coronary embolization. Longitudinal and radial strain of the hyperenhanced patchy microinfarction demonstrates persistent decline over the course of 1 week.