Cardiac motion and deformation recovery from MRI: a review

Non-invasive Model-Based Assessment of Passive Left-Ventricular Myocardial Stiffness in Healthy Subjects and in Patients with Non-ischemic Dilated Cardiomyopathy

Patient-specific modelling has emerged as a tool for studying heart function, demonstrating the potential to provide non-invasive estimates of tissue passive stiffness. However, reliable use of model-derived stiffness requires sufficient model accuracy and unique estimation of model parameters. In this paper we present personalised models of cardiac mechanics, focusing on improving model accuracy, while ensuring unique parametrisation. The influence of principal model uncertainties on accuracy and parameter identifiability was systematically assessed in a group of patients with dilated cardiomyopathy ([Formula: see text]) and healthy volunteers ([Formula: see text]). For all cases, we examined three circumferentially symmetric fibre distributions and two epicardial boundary conditions. Our results demonstrated the ability of data-derived boundary conditions to improve model accuracy and highlighted the influence of the assumed fibre distribution on both model fidelity and stiffness estimates. The model personalisation pipeline-based strictly on non-invasive data-produced unique parameter estimates and satisfactory model errors for all cases, supporting the selected model assumptions. The thorough analysis performed enabled the comparison of passive parameters between volunteers and dilated cardiomyopathy patients, illustrating elevated stiffness in diseased hearts.

Non-invasive Model-Based Assessment of Passive Left-Ventricular Myocardial Stiffness in Healthy Subjects and in Patients with Non-ischemic Dilated Cardiomyopathy

Patient-specific modelling has emerged as a tool for studying heart function, demonstrating the potential to provide non-invasive estimates of tissue passive stiffness. However, reliable use of model-derived stiffness requires sufficient model accuracy and unique estimation of model parameters. In this paper we present personalised models of cardiac mechanics, focusing on improving model accuracy, while ensuring unique parametrisation. The influence of principal model uncertainties on accuracy and parameter identifiability was systematically assessed in a group of patients with dilated cardiomyopathy (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n=3$$\end{document}n=3) and healthy volunteers (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n=5$$\end{document}n=5). For all cases, we examined three circumferentially symmetric fibre distributions and two epicardial boundary conditions. Our results demonstrated the ability of data-derived boundary conditions to improve model accuracy and highlighted the influence of the assumed fibre distribution on both model fidelity and stiffness estimates. The model personalisation pipeline—based strictly on non-invasive data—produced unique parameter estimates and satisfactory model errors for all cases, supporting the selected model assumptions. The thorough analysis performed enabled the comparison of passive parameters between volunteers and dilated cardiomyopathy patients, illustrating elevated stiffness in diseased hearts.

Automated tracking, segmentation and trajectory classification of pelvic organs on dynamic MRI

Pelvic organ prolapse is a major health problem in women where pelvic floor organs (bladder, uterus, small bowel, and rectum) fall from their normal position and bulge into the vagina. Dynamic Magnetic Resonance Imaging (DMRI) is presently used to analyze the organs' movements from rest to maximum strain providing complementary support for diagnosis. However, there is currently no automated or quantitative approach to measure the movement of the pelvic organs and their correlation with the severity of prolapse. In this paper, a two-stage method is presented to automatically track and segment pelvic organs on DMRI followed by a multiple-object trajectory classification method to improve the diagnosis of pelvic organ prolapse. Organs are first tracked using particle filters and K-means clustering with prior information. Then, they are segmented using the convex hull of the cluster of particles. Finally, the trajectories of the pelvic organs are modeled using a new Coupled Switched Hidden Markov Model (CSHMM) to classify the severity of pelvic organ prolapse. The tracking and segmentation results are validated using Dice Similarity Index (DSI) whereas the classification results are compared with two manual clinical measurements. Results demonstrate that the presented method is able to automatically track and segment pelvic organs with a DSI above 82% for 26 out of 46 cases and DSI above 75% for all 46 tested cases. The accuracy of the trajectory classification model is also better than current manual measurements.

Regional infarction identification from cardiac CT images: a computer-aided biomechanical approach

PURPOSE

Regional infarction identification is important for heart disease diagnosis and management, and myocardial deformation has been shown to be effective for this purpose. Although tagged and strain-encoded MR images can provide such measurements, they are uncommon in clinical routine. On the contrary, cardiac CT images are more available with lower costs, but they only provide motion of cardiac boundaries and additional constraints are required to obtain the myocardial strains. The goal of this study is to verify the potential of contrast-enhanced CT images on computer-aided regional infarction identification.

METHODS

We propose a biomechanical approach combined with machine learning algorithms. A hyperelastic biomechanical model is used with deformable image registration to estimate 3D myocardial strains from CT images. The regional strains and CT image intensities are input to a classifier for regional infarction identification. Cross-validations on ten canine image sequences with artificially induced infarctions were used to study the performances of using different feature combinations and machine learning algorithms.

RESULTS

Radial strain, circumferential strain, first principal strain, and image intensity were shown to be discriminative features. The highest identification accuracy ([Formula: see text] %) was achieved when combining radial strain with image intensity. Random forests gave better results than support vector machines on less discriminative features. Random forests also performed better when all strains were used together.

CONCLUSION

Although CT images cannot directly measure myocardial deformation, with the use of a biomechanical model, the estimated strains can provide promising identification results especially when combined with CT image intensity.

Left ventricular diastolic dysfunction in type 2 diabetes patients: a novel 2D strain analysis based on cardiac magnetic resonance imaging

OBJECTIVE

This study was to develop a strain analysis method to evaluate the left ventricular (LV) functions in type 2 diabetic patients with an asymptomatic LV diastolic dysfunction.

METHODS

Two groups (10 asymptomatic type 2 diabetic subjects and 10 control ones) were considered. All of the subjects had normal ejection fraction values but impaired diastolic functions assessed by the transmitral blood flow velocity. For each subject, based on cardiac MRI, global indexes including LV volume, LV myocardial mass, cardiac index (CI), and transmitral peak velocity, were measured, and regional indexes (i.e., LV deformation, strain and strain rate) were calculated through an image-registration technology.

RESULTS

Most of the global indexes did not differentiate between the two groups, except for the CI, LV myocardial mass and transmitral peak velocity. While for the regional indexes, the global LV diastolic dysfunction of the diabetic indicated an increased strain (0.08 ± 0.044 vs. -0.031 ± 0.077, p = 0.001) and a reduced strain rate (1.834 ± 0.909 vs. 3.791 ± 2.394, p = 0.033) compared to the controls, moreover, the local LV diastolic dysfunction reflected by the strain and strain rate varied, and the degree of dysfunction gradually decreased from the basal level to the apical level.

CONCLUSIONS

The results showed that the strain and strain rates are effective to capture the subtle alterations of the LV functions, and the proposed method can be used to estimate the LV myocardial function based on cardiac MRI.

Estimation of Myocardial Strain and Contraction Phase From Cine MRI Using Variational Data Assimilation

This paper presents a new method to estimate left ventricle deformations using variational data assimilation that combines image observations from cine MRI and a dynamic evolution model of the heart. The main contribution of the model is that it embeds parameters modeling the contraction / relaxation process. It estimates myocardial motion and contraction parameters simultaneously, providing accurate complementary information for diagnosis. The method was applied to synthetic datasets with known ground truth motion and to 47 patients MRI datasets acquired at three slice locations (base, mid-ventricle and apex). Radial and circumferential strain components were compared to those obtained with a reference tag tracking software, exhibiting good agreement with intraclass correlation coefficients (ICC) above 0.8. Results were also evaluated against wall motion score indices used to assess cardiac kinetics in clinical practice. The assimilation process overcame issues caused by temporal artifacts as a result of the dynamic model, compared to using the observation term alone. Moreover we found that the new dynamic model, consisting of a piecewise transport model acting independently on systole and diastole performed better than the standard continuous transport model, which oversmooths temporal variations. Estimated strain and contraction parameters significantly correlated to clinical scores, making them promising features for diagnosing not only hypokinesia but also dyskinesia.

Dictionary-Driven Ischemia Detection From Cardiac Phase-Resolved Myocardial BOLD MRI at Rest

Cardiac Phase-resolved Blood-Oxygen-Level Dependent (CP-BOLD) MRI provides a unique opportunity to image an ongoing ischemia at rest. However, it requires post-processing to evaluate the extent of ischemia. To address this, here we propose an unsupervised ischemia detection (UID) method which relies on the inherent spatio-temporal correlation between oxygenation and wall motion to formalize a joint learning and detection problem based on dictionary decomposition. Considering input data of a single subject, it treats ischemia as an anomaly and iteratively learns dictionaries to represent only normal observations (corresponding to myocardial territories remote to ischemia). Anomaly detection is based on a modified version of One-class Support Vector Machines (OCSVM) to regulate directly the margins by incorporating the dictionary-based representation errors. A measure of ischemic extent (IE) is estimated, reflecting the relative portion of the myocardium affected by ischemia. For visualization purposes an ischemia likelihood map is created by estimating posterior probabilities from the OCSVM outputs, thus obtaining how likely the classification is correct. UID is evaluated on synthetic data and in a 2D CP-BOLD data set from a canine experimental model emulating acute coronary syndromes. Comparing early ischemic territories identified with UID against infarct territories (after several hours of ischemia), we find that IE, as measured by UID, is highly correlated (Pearson's r=0.84) with respect to infarct size. When advances in automated registration and segmentation of CP-BOLD images and full coverage 3D acquisitions become available, we hope that this method can enable pixel-level assessment of ischemia with this truly non-invasive imaging technique.

Myocardial strain assessment by cine cardiac magnetic resonance imaging using non-rigid registration

AIMS

To evaluate a novel post-processing method for assessment of longitudinal mid-myocardial strain in standard cine cardiac magnetic resonance (CMR) imaging sequences.

METHODS AND RESULTS

Cine CMR imaging and tagged cardiac magnetic resonance imaging (TMRI) were performed in 15 patients with acute myocardial infarction (AMI) and 15 healthy volunteers served as control group. A second group of 37 post-AMI patients underwent both cine CMR and late gadolinium enhancement (LGE) CMR exams. Speckle tracking echocardiography (STE) was performed in 36 of these patients. Cine CMR, TMRI and STE were analyzed to obtain longitudinal strain. LGE-CMR datasets were analyzed to evaluate scar extent. Comparison of peak systolic strain (PSS) measured from CMR and TMRI yielded a strong correlation (r=0.86, p<0.001). PSS measured from CMR and STE correlated well (r=0.75, p<0.001). A cutoff longitudinal PSS value of -13.14% differentiated non-infarction from any infarcted myocardium, with a sensitivity of 93% and a specificity of 89% (area under curve (AUC) 0.95). PSS value of -9.39% differentiated non-transmural from transmural infarcted myocardium, with a sensitivity of 75% and a specificity of 67% (AUC 0.78).

CONCLUSION

The present study showed a novel off-line post-processing method for segmental longitudinal strain analysis in mid-myocardium layer based on cine CMR data. The method was found to be highly correlated with strain measurements obtained by TMRI and STE. This tool allows accurate discrimination between different transmurality states of myocardial infarction.

Spatio-Temporal Tensor Decomposition of a Polyaffine Motion Model for a Better Analysis of Pathological Left Ventricular Dynamics

Given that heart disease can cause abnormal motion dynamics over the cardiac cycle, understanding and quantifying cardiac motion can provide insight for clinicians to aid with diagnosis, therapy planning, and determining prognosis. The goal of this paper is to extract population-specific motion patterns from 3D displacements in order to identify the mean motion in a population, and to describe pathology-specific motion patterns in terms of the spatial and temporal components. Since there are common motion patterns observed in patients with the same condition, extracting these can lead towards a better understanding of the disease. Quantifying cardiac motion at a population level is not a simple task since images can vary widely in terms of image quality, size, resolution, and pose. To overcome this, we analyze the parameters obtained from a cardiac-specific Polyaffine motion-tracking algorithm, which are aligned both spatially and temporally to a common reference space. Once all parameters are aligned, different subjects can be compared and analyzed in the space of Polyaffine transformations by projecting the transformations to a reduced order subspace in which dominant motion patterns in each population can be extracted. Using tensor decomposition, the spatial and temporal aspects can be decoupled in order to study the components individually. The proposed method was validated on healthy volunteers and Tetralogy of Fallot patients according to known spatial and temporal behavior for each population. A key advantage of this method is the ability to regenerate motion sequences from the models, which can be visualized in terms of the full motion.

MRI-based strain and strain rate analysis of left ventricle: a modified hierarchical transformation model

BACKGROUND

Different from other indicators of cardiac function, such as ejection fraction and transmitral early diastolic velocity, myocardial strain is promising to capture subtle alterations that result from early diseases of the myocardium. In order to extract the left ventricle (LV) myocardial strain and strain rate from cardiac cine-MRI, a modified hierarchical transformation model was proposed.

METHODS

A hierarchical transformation model including the global and local LV deformations was employed to analyze the strain and strain rate of the left ventricle by cine-MRI image registration. The endocardial and epicardial contour information was introduced to enhance the registration accuracy by combining the original hierarchical algorithm with an Iterative Closest Points using Invariant Features algorithm. The hierarchical model was validated by a normal volunteer first and then applied to two clinical cases (i.e., the normal volunteer and a diabetic patient) to evaluate their respective function.

RESULTS

Based on the two clinical cases, by comparing the displacement fields of two selected landmarks in the normal volunteer, the proposed method showed a better performance than the original or unmodified model. Meanwhile, the comparison of the radial strain between the volunteer and patient demonstrated their apparent functional difference.

CONCLUSIONS

The present method could be used to estimate the LV myocardial strain and strain rate during a cardiac cycle and thus to quantify the analysis of the LV motion function.

Strain analysis from 4-D cardiac CT image data

Strain is a discriminative parameter of regional myocardial dysfunction. Despite the large body of research on myocardial strain analysis in echocardiography and MR images, such techniques have not often been applied to cardiac CT data. Reasons for this include the challenges of sparse image deformation clues and the low temporal resolution. In the current study, we propose an algorithm that uses cardiac CT data to evaluate the mechanical function of the left ventricle. The algorithm is based on a deformable LV model that contains both the myocardium and the blood pool regions and that accounts for the elasticity and incompressibility of the myocardium with the rapid contraction of the blood pool. Our algorithm uses the image intensities of the trabeculle and papillary muscles as well as the border edges in an optical flow manner to extract the 3-D velocities. The resulting strains and rotational values derived from a set of normal patients correlate highly with values from the research literature. We validated our algorithm against 2-D speckle tracking analysis and against visual scores obtained by an expert. Our study shows that strain analysis using CT data can be used in clinical practice.

Images as drivers of progress in cardiac computational modelling

Computational models have become a fundamental tool in cardiac research. Models are evolving to cover multiple scales and physical mechanisms. They are moving towards mechanistic descriptions of personalised structure and function, including effects of natural variability. These developments are underpinned to a large extent by advances in imaging technologies. This article reviews how novel imaging technologies, or the innovative use and extension of established ones, integrate with computational models and drive novel insights into cardiac biophysics. In terms of structural characterization, we discuss how imaging is allowing a wide range of scales to be considered, from cellular levels to whole organs. We analyse how the evolution from structural to functional imaging is opening new avenues for computational models, and in this respect we review methods for measurement of electrical activity, mechanics and flow. Finally, we consider ways in which combined imaging and modelling research is likely to continue advancing cardiac research, and identify some of the main challenges that remain to be solved.

A robust and accurate center-frequency estimation (RACE) algorithm for improving motion estimation performance of SinMod on tagged cardiac MR images without known tagging parameters

A robust and accurate center-frequency (CF) estimation (RACE) algorithm for improving the performance of the local sine-wave modeling (SinMod) method, which is a good motion estimation method for tagged cardiac magnetic resonance (MR) images, is proposed in this study. The RACE algorithm can automatically, effectively and efficiently produce a very appropriate CF estimate for the SinMod method, under the circumstance that the specified tagging parameters are unknown, on account of the following two key techniques: (1) the well-known mean-shift algorithm, which can provide accurate and rapid CF estimation; and (2) an original two-direction-combination strategy, which can further enhance the accuracy and robustness of CF estimation. Some other available CF estimation algorithms are brought out for comparison. Several validation approaches that can work on the real data without ground truths are specially designed. Experimental results on human body in vivo cardiac data demonstrate the significance of accurate CF estimation for SinMod, and validate the effectiveness of RACE in facilitating the motion estimation performance of SinMod.

Deformable models with sparsity constraints for cardiac motion analysis

Deformable models integrate bottom-up information derived from image appearance cues and top-down priori knowledge of the shape. They have been widely used with success in medical image analysis. One limitation of traditional deformable models is that the information extracted from the image data may contain gross errors, which adversely affect the deformation accuracy. To alleviate this issue, we introduce a new family of deformable models that are inspired from the compressed sensing, a technique for accurate signal reconstruction by harnessing some sparseness priors. In this paper, we employ sparsity constraints to handle the outliers or gross errors, and integrate them seamlessly with deformable models. The proposed new formulation is applied to the analysis of cardiac motion using tagged magnetic resonance imaging (tMRI), where the automated tagging line tracking results are very noisy due to the poor image quality. Our new deformable models track the heart motion robustly, and the resulting strains are consistent with those calculated from manual labels.

Sparse deformable models with application to cardiac motion analysis

Deformable models have been widely used with success in medical image analysis. They combine bottom-up information derived from image appearance cues, with top-down shape-based constraints within a physics-based formulation. However, in many real world problems the observations extracted from the image data often contain gross errors, which adversely affect the deformation accuracy. To alleviate this issue, we introduce a new family of deformable models that are inspired from compressed sensing, a technique for efficiently reconstructing a signal based on its sparseness in some domain. In this problem, we employ sparsity to represent the outliers or gross errors, and combine it seamlessly with deformable models. The proposed new formulation is applied to the analysis of cardiac motion, using tagged magnetic resonance imaging (tMRI), where the automated tagging line tracking results are very noisy due to the poor image quality. Our new deformable models track the heart motion robustly, and the resulting strains are consistent with those calculated from manual labels.

MRI strain imaging of the carotid artery: present limitations and future challenges

Rupture of atherosclerotic plaques in the carotid artery is a main cause of stroke. Current diagnostics are not sufficient to identify all rupture-prone plaques, and studies have shown that biomechanical factors improve current plaque risk assessment. Strain imaging may be a valuable contribution to this risk assessment. MRI is a versatile imaging technique that offers various methods that are capable of measuring tissue strain. In this review, MR imaging techniques with displacement (DENSE), velocity (PC MRI), or strain (SENC) encoding protocols are discussed, together with post-processing techniques based on time-resolved MRI data. Although several MRI techniques are being developed to improve time-resolved MR imaging, current technical limitations related to spatial and temporal resolutions render MRI strain imaging currently unfit for carotid plaque strain evaluation. A novel approach using non-rigid image registration of MR images to determine strain in carotid arteries based on black blood cine MRI is proposed in this review. This and other post-processing techniques based on time-resolved MRI data may provide a good estimate of plaque strain, but are also dependent on the spatial and temporal resolution of the MR images. However, they seem to be the most promising approach for MRI based plaque strain analysis in the near future.

Groupwise elastic registration by a new sparsity-promoting metric: application to the alignment of cardiac magnetic resonance perfusion images

This paper proposes a methodology for the joint alignment of a sequence of images based on a groupwise registration procedure by using a new family of metrics that exploit the expected sparseness of the temporal intensity curves corresponding to the aligned points. Therefore, this methodology is able to tackle the alignment of temporal sequences of images in which the represented phenomenon varies in time. Specifically, we have applied it to the correction of motion in contrast-enhanced first-pass perfusion cardiac magnetic resonance images. The time sequence is elastically registered as a whole by using the aforementioned family of multi-image metrics and jointly optimizing the parameters of the transformations involved. The proposed metrics are able to cope with dynamic changes in the intensity content of corresponding points in the sequence guided by the assumption that these changes allow for a sparse representation in a properly selected frame. Results have shown the statistically significant improvement in the performance of the proposed metric with respect to previous groupwise registration metrics for the problem at hand, which is especially relevant to correct for elastic deformations.

Analysis of 3D cardiac deformations with 3D SinMod

In this paper, we propose a novel 3D sine wave modeling (3D SinMod) approach to automatic analysis of 3D cardiac deformations. An accelerated 3D complementary spatial modulation of magnetization (CSPAMM) tagging technique was used to modulate the myocardial tissue and to acquire 3D MR data sets of the whole-heart including three orthog- onal tags within three breath-holds. Each tag set is able to assess the motion along a direction perpendicular to the tag lines. With the application of CSPAMM, the effect of tag fading encountered in SPAMM tagging due to T1 relaxation is mitigated and tag deformations can be visualized for the entire cardiac cycle, including diastolic phases. In the proposed approach, the environment around each voxel in the 3D volume is modeled as a moving sine wavefront with local frequency and amplitude. The biggest advantage of the proposed technique is that the entire framework, from data acquisition to data analysis is in the 3D domain, which permits quantification of both the in-plane and through-plane motion components. The accuracy and the effectiveness of the proposed method has been validated using both simulated and in vivo tag data.

Imaging techniques for cardiac strain and deformation: comparison of echocardiography, cardiac magnetic resonance and cardiac computed tomography

Myocardial function assessment is essential for determining the health of the myocardium. Global assessment of myocardial function is widely performed (by estimating the ejection fraction), but many common cardiac diseases initially affect the myocardium on a regional, rather than global basis. Regional myocardial wall motion can be quantified using myocardial strain analysis (a normalized measure of deformation). Myocardial strain can be measured in terms of three normal strains (longitudinal strain, radial strain and circumferential) and six shear strains. Cardiac MRI (cMRI) is usually considered the reference standard for measurement of myocardial strain. The most common cMRI method, termed tagged cMRI, allows full, 3D assessment of regional strain. However, due to its complexity and lengthy times for analysis, tagged cMRI is not usually used outside of academic centers. Tagged cMRI is also primarily used only in research studies. Echocardiography combined with tissue Doppler imaging or a speckle tracking technique is now widely available in the clinical setting. Myocardial strain measurement by echocardiography shows reasonable agreement with cMRI. Limited standardization and differences between vendors represent current limitations of the technique. Cardiac computed tomography (CCT) is the newest and most rapidly growing modality for noninvasive imaging of the heart. While CCT studies are most commonly applied to assess the coronary arteries, CCT is easily adapted to provide functional information for both the left and right ventricles. New methods for CCT assessment of regional myocardial function are being developed. This review outlines the current literature on imaging techniques related to cardiac strain analysis and discusses the strengths and weaknesses of various methods for myocardial strain analysis.

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.

Impact of the Advisa MRI pacing system on the diagnostic quality of cardiac MR images and contraction patterns of cardiac muscle during scans: Advisa MRI randomized clinical multicenter study results

BACKGROUND

The Advisa MRI system is designed to safely undergo magnetic resonance imaging (MRI). Its influence on image quality is not well known.

OBJECTIVE

To evaluate cardiac magnetic resonance (CMR) image quality and to characterize myocardial contraction patterns by using the Advisa MRI system.

METHODS

In this international trial with 35 participating centers, an Advisa MRI system was implanted in 263 patients. Of those, 177 were randomized to the MRI group and 150 underwent MRI scans at the 9-12-week visit. Left ventricular (LV) and right ventricular (RV) cine long-axis steady-state free precession MR images were graded for quality. Signal loss along the implantable pulse generator and leads was measured. The tagging CMR data quality was assessed as the percentage of trackable tagging points on complementary spatial modulation of magnetization acquisitions (n=16) and segmental circumferential fiber shortening was quantified.

RESULTS

Of all cine long-axis steady-state free precession acquisitions, 95% of LV and 98% of RV acquisitions were of diagnostic quality, with 84% and 93%, respectively, being of good or excellent quality. Tagging points were trackable from systole into early diastole (360-648 ms after the R-wave) in all segments. During RV pacing, tagging demonstrated a dyssynchronous contraction pattern, which was not observed in nonpaced (n = 4) and right atrial-paced (n = 8) patients.

CONCLUSIONS

In the Advisa MRI study, high-quality CMR images for the assessment of cardiac anatomy and function were obtained in most patients with an implantable pacing system. In addition, this study demonstrated the feasibility of acquiring tagging data to study the LV function during pacing.

Generating anatomical models of the heart and the aorta from medical images for personalized physiological simulations

The anatomy and motion of the heart and the aorta are essential for patient-specific simulations of cardiac electrophysiology, wall mechanics and hemodynamics. Within the European integrated project euHeart, algorithms have been developed that allow to efficiently generate patient-specific anatomical models from medical images from multiple imaging modalities. These models, for instance, account for myocardial deformation, cardiac wall motion, and patient-specific tissue information like myocardial scar location. Furthermore, integration of algorithms for anatomy extraction and physiological simulations has been brought forward. Physiological simulations are linked closer to anatomical models by encoding tissue properties, like the muscle fibers, into segmentation meshes. Biophysical constraints are also utilized in combination with image analysis to assess tissue properties. Both examples show directions of how physiological simulations could provide new challenges and stimuli for image analysis research in the future.

Validation of continuously tagged MRI for the measurement of dynamic 3D skeletal muscle tissue deformation

PURPOSE

Typically spatial modulation of the magnetization (SPAMM) tagged magnetic resonance imaging (MRI) requires many repeated motion cycles limiting the applicability to highly repeatable tissue motions only. This paper describes the validation of a novel SPAMM tagged MRI and post-processing framework for the measurement of complex and dynamic 3D soft tissue deformation following just three motion cycles. Techniques are applied to indentation induced deformation measurement of the upper arm and a silicone gel phantom.

METHODS

A SPAMM tagged MRI methodology is presented allowing continuous (3.3-3.6 Hz) sampling of 3D dynamic soft tissue deformation using non segmented 3D acquisitions. The 3D deformation is reconstructed by the combination of three mutually orthogonal tagging directions, thus requiring only three repeated motion cycles. In addition a fully automatic post-processing framework is presented employing Gabor scale-space and filter-bank analysis for tag extrema segmentation and triangulated surface fitting aided by Gabor filter bank derived surface normals. Deformation is derived following tracking of tag surface triplet triangle intersections. The dynamic deformation measurements were validated using indentation tests (∼20 mm deep at 12 mm/s) on a silicone gel soft tissue phantom containing contrasting markers which provide a reference measure of deformation. In addition, the techniques were evaluated in vivo for dynamic skeletal muscle tissue deformation measurement during indentation of the biceps region of the upper arm in a volunteer.

RESULTS

For the phantom and volunteer tag point location precision were 44 and 92 μm, respectively resulting in individual displacements precisions of 61 and 91 μm, respectively. For both the phantom and volunteer data cumulative displacement measurement accuracy could be evaluated and the difference between initial and final locations showed a mean and standard deviation of 0.44 and 0.59 mm for the phantom and 0.40 and 0.73 mm for the human data. Finally accuracy of (cumulative) displacement was evaluated using marker tracking in the silicone gel phantom. Differences between true and predicted marker locations showed a mean of 0.35 mm and a standard deviation of 0.63 mm.

CONCLUSIONS

A novel SPAMM tagged MRI and fully automatic post-processing framework for the measurement of complex 3D dynamic soft tissue deformation following just three repeated motion cycles was presented. The techniques demonstrate dynamic measurement of complex 3D soft tissue deformation at subvoxel accuracy and precision and were validated for 3.3-3.6 Hz sampling of deformation speeds up to 12 mm/s.