Construction of a two-parameter empirical model of left ventricle wall motion using cardiac tagged magnetic resonance imaging data
© Shi et al.; licensee BioMed Central Ltd. 2012
Received: 5 April 2012
Accepted: 1 October 2012
Published: 24 October 2012
A one-parameter model was previously proposed to characterize the short axis motion of the LV wall at the mid-ventricle level. The single parameter of this model was associated with the radial contraction of myocardium, but more comprehensive model was needed to account for the rotation at the apex and base levels. The current study developed such model and demonstrated its merits and limitations with examples.
Materials and methods
The hearts of five healthy individuals were visualized using cardiac tagged magnetic resonance imaging (tMRI) covering the contraction and relaxation phases. Based on the characteristics of the overall dynamics of the LV wall, its motion was represented by a combination of two components - radial and rotational. Each component was represented by a transformation matrix with a time-dependent variable α or β.
Image preprocessing step and model fitting algorithm were described and applied to estimate the temporal profiles of α and β within a cardiac cycle at the apex, mid-ventricle and base levels. During this process, the tagged lines of the acquired images served as landmark reference for comparing against the model prediction of the motion. Qualitative and quantitative analyses were performed for testing the performance of the model and thus its validation.
The α and β estimates exhibited similarities in values and temporal trends once they were scaled by the radius of the epicardium (r epi )and plotted against the time scaled by the period of the cardiac cycle (T cardiac ) of each heart measured during the data acquisition. α/r epi peaked at about Δt/T cardiac =0.4 and with values 0.34, 0.4 and 0.3 for the apex, mid-ventricle and base level, respectively. β/r epi similarly maximized in amplitude at about Δt/T cardiac =0.4, but read 0.2 for the apex and - 0.08 for the base level. The difference indicated that the apex twisted more than the base.
It is feasible to empirically model the spatial and temporal evolution of the LV wall motion using a two-parameter formulation in conjunction with tMRI-based visualization of the LV wall in the transverse planes of the apex, mid-ventricle and base. In healthy hearts, the analytical model will potentially allow deriving biomechanical entities, such as strain, strain rate or torsion, which are typically used as diagnostic, prognostic or predictive markers of cardiovascular diseases including diabetes.
KeywordsCardiac modeling Myocardial deformation Left ventricle wall motion Tagged MRI
Research in biomechanical modeling of left ventricle (LV) aims to determine sensitive and specific entities that can serve as clinical biomarkers for evaluating cardiac function or diagnosing an abnormality. With this in mind, rigorous computational models have been developed for mimicking the LV wall motion and determining global and regional parameters for assessing cardiac performance [1–6]. Any measured deviation in a model parameter of interest from its normative value is considered as an indication of an abnormality and offers key diagnostics for detecting dysfunction associated with heart disease.
In general, the LV wall motion of a healthy heart in short-axis view is a complex event involving actions of radial and rotational components during both systolic contraction and diastolic relaxation phases of the cardiac cycle . However, the LV wall at mid-ventricle moves mostly in radial direction, meaning it contracts isotropically with negligible amount of rotation. Based on these unique characteristics, a simple one-parameter model was developed earlier to mathematically represent the motion field at the mid-ventricle . The current study was under taken to increase the capacity of this empirical model further to account for the twist-type motions of the apex and base. Specifically, a new parameter was introduced into the previous model. As before, the new proposed model with now two parameters was tested if it could sufficiently accomplish this task. The validation work was performed with the cardiac data acquired from healthy human hearts using tagged magnetic resonance imaging (tMRI) modality .
In the following, we first describe the assumptions and considerations that led to the construction of the two-parameter model. Next, we summarize the tMRI data acquisition protocol for visualizing the heart under clinical situation and introduce image preprocessing and registration algorithms for fitting the model to the data. Then, we give details of the construction of the model and demonstrate its capabilities with examples. Finally, we discuss the results and state our conclusions.
Two-parameter model of the LV wall motion
Justification of the model
LV shape is governed by a cylindrically symmetric geometry along its long-axis.
Origin of the coordinates is positioned at the center of the LV wall in the short-axis view.
Beginning of the systolic contraction corresponds to the end-diastolic phase denoted by time t=0. At this moment, the LV wall appears at its most dilated state and experiences no deformation or torsion.
Vectors r=(x,y,t=0) and r ′=(x ′,y′ t) denote Cartesian coordinates of the myocardial tissue in the LV wall in the short-axis plane before (at the end-diastole) and after going through the radial motion and/or rotation during either the contraction or relaxation phases of the LV. The difference r ′ r denotes displacement. We note that this time-dependent vectorial representation is consistent with Lagrangian description of the motion .
Magnitudes and are distances to the origin before and after the motion of the myocardial tissue, respectively.
The new spatial coordinates of a moving myocardial tissue are written in terms of its position at the end-diastole. This requires breaking up the motion into its radial and rotational components. Each component is independently represented by a combination of an isotropic radial transformation and a non-uniform rotational transformation in two-dimensional space, respectively.
The shortening of the LV wall along its long axis during the contraction or relaxation produces negligible amount of out-of-plane motion.
Formulation of the isotropic radial transformation
α(t) is a time-dependent parameter. The minus sign in the diagonal elements of T indicates contraction. Rearranging Eq. (1) yields radial tissue displacement |r ′ − r| = α2/|r| that is inversely proportional to the distance to the origin. 1/r-dependence in the displacement reflects that the motion is more significant near the inner boundary (endocardium) than the outer boundary (epicardium) of the LV wall. Larger α yields larger displacement. Similarly, |r ' − r|/|r| = α2/|r|2 indicates that the relative displacement is inversely proportional to the square of the distance to the origin.
Formulation of the non-uniform rotational transformation
where , and .
However, close examination of the in vivo images of heart indicates that the myocardium goes through a lesser degree of rotation near the endocardium than the epicardium . Our empirical analysis of the images in polar coordinates suggested that the radially non-uniform behavior of the rotation can be absorbed into the formulation by the substitution Θ=β(t)/|r|. In this notation, a positive (or negative) β corresponds to a clockwise (or counterclockwise) rotation. Larger β corresponds to an increase in rotation, but β=0 represents no rotation as in the case experienced by the myocardium of the mid-ventricle.
The time-dependent variables α and β together constitute the two parameters of the model whose values at time t to be estimated by following the algorithm described below.
It should be noted that in order to avoid nonlinear entanglement of the cascaded transformations in Eq. (3), Eq. (1) is first applied orderly before Eq. (2) since 1/r is invariant under rotation but the angle Θ changes with the application of the radial transformation in Eq. (2). This ordering reduces the complexity of the computational process in estimating the parameters α and β to describe the underlying real motion of the LV wall. The inverse transformation associated with the time-reversal operation was given in Appendix.
In practice, to successfully characterize the spatial and temporal profile of the LV wall motion or the other features of the motion dynamics, the analytical model in Eq. (3) needs to be validated by comparing the model based motion predictions against the real motion of the LV wall during the different phases of the cardiac cycle - while contracting in the systolic phase or relaxing in the diastolic phase. The steps involved in the validation process are described next.
Materials and methods
tMRI data acquisition
Preprocessing of the tMRI data
An algorithm was given earlier for the processing of the acquired tMRI data sets prior to performing the motion analysis on the LV wall . But, the current study applied a slightly different preprocessing approach to reduce the level of complexity. This was achieved by normalizing the variables in the formulation by the key parameters describing the shape and size of the LV wall. The preprocessing step involved locating the center (long-axis) of the LV wall and its epicardial and endocardial radii in every 1st end-diastolic image from different slices. For example, to accomplish this task at the slice corresponding to the mid-ventricle level, a series of stroboscopic images were first selected from the underlying slice location. The 1st image frame in the series was displayed to depict the end-diastolic LV on the computer screen. The LV wall on this image appeared at its largest cross sectional dimension compared to the remaining images in the series. Then, a circle was aligned with the epicardium of the LV wall. The center and radius of the circle were taken as the center and epicardial wall radius (r epi ) of the LV at the end-diastole. r endo was similarly estimated from a smaller circle by aligning its circumference with endocardium. These parameters for the LV wall in the other slices at the apex and base levels were also measured similarly by repeating the above process on the corresponding image frames.
Validation of the model
where (x′ i ,y′ i ) and N denote the coordinate and number of pixels of the calculated mesh, respectively. Note that F(α,β,t) ranged from 0 for a perfect match to 1 for a completely mismatch between the calculated mesh and the tags on the LV wall. A minimization of F(α,β,t) yielded a set of “best-fit” for the parameters α and β. In the model fitting code, the minimization of F(α,β,t) was accomplished by scanning the two-dimensional α–β space in relatively narrow range as we were only seeking α and β with two or three-digit accuracy.
The use of the tMRI modality was instrumental in building and validating the two-parameter model since the LV wall was conveniently visualized during its contraction and relaxation. Since the tag lines in the acquired images preserved the material coordinates, the lines provided a true system of spatial reference for tracking the motion and matching the simulated grid mesh for determining the model parameters α and β. These parameters could have also been estimated from images acquired with other imaging modalities such as cardiac cine imaging or echography using standard block matching algorithms that work on image intensity features. However, alterations in the intensity pattern of the cine images during the cardiac cycle and speckle decorrelation of the myocardial echogeneity in these modalities would yield estimates with lower precision.
Extensive visual assessments of the images such as those in Figure 4 indicated that the calculated meshes coded by red color mostly overlapped with the tag cells in the LV walls. These qualitative evaluations, combined with the quantitative estimates of the model parameters in Figures 5 and 6, collectively confirmed the validity of the proposed two-parameter model and its capability of representing the underlying LV wall motion reasonably well. The optimization algorithm produced robust and accurate estimates of the parameters α and β at different time points of the cardiac cycle. This might be a generic property of the model that may be extended to cover the LV wall motion in other species. From this aspect alone, the model can potentially be a valuable tool in preclinical translational studies aimed at better understanding and comparing the LV wall motion in experimental studies. The scaling of α and β by r epi introduced above helps these efforts as it minimizes the intra and inter species variations between the LV wall sizes and the cardiac cycle periods.
Active research focuses on determining specific and sensitive biomarkers at global or regional level for quantitatively evaluating the LV performance under normal or pathological conditions. The two-parameters α and β governing the model can serve this purpose. The model formula in Eq. (3) allows calculating the position of the myocardial tissue in the LV wall at any given spatial point and time in the short-axis plane during the cardiac cycle after the substitution of the α and β estimates promptly. The resulting analytical expressions governing the model can be evaluated to derive additional parameters as surrogate biomarkers for assessing the ongoing cardiac events. For example, from Eq. (1), the ratio α2/r end represents the change in the radius (r-r′) and the reduction of the cross-sectional area of the LV cavity in the short-axis plane is π(r2-r′ 2)| r=end . Using the approximation (r2-r′ 2) ~ (r+r′) (r-r′) ~ 2r(r-r′), α2 can be seen as proportional to the reduction of the LV cavity area due to the contraction. The characteristics of such markers could potentially be developed into critical differential diagnostics between the healthy and diseased hearts. One possibility of potential application is in evaluating diabetic heart where the myocardial tissue contraction is modified due to fibrotic tissue deposition. Other possibilities are cases where the heart goes through remodeling after infarct or myocardial hypertrophy. The current preliminary study was, however, performed on a limited number of subjects. By increasing the number of subjects further and performing statistical analysis on the α and β estimates can yield group averages to better embody the mean biomechanics of the LV wall motion as normative measures in healthy individuals. Therefore, evaluations of how well the model would perform in the circumstances of cardiac disease and further extension of the model to represent the overall LV wall motion in every transverse plane remain to future studies.
We note that the model is valid for representing the global cardiac motion only. This presents a limitation since it does not distinguish the differences in regional motions occuring at different segments of the short-axis of the LV wall due to local myocardial abnormalities. Hence, all the regional cardiac function measurements, such as displacement, strain and strain rate, using the model would remain the same for each segment.
Another limitation of the model is that it fails to account for the counter-clockwise rotation of the apex that occurs early on during systole. Insufficient temporal resolution between the image frames may explain such aspect of the apex motion was not present in our data.
Our model fitting procedure involved converting tagged images to binary images. The tag lines typically fade towards the end of the cardiac cycle in the image series. This is a common event and the tag line fading gets worse with a shorter T1 relaxation time of the underlying myocardial tissue or longer period of cardiac cycle. Our study was conducted on a 3 T clinical scanner. T1 measures longer with increased magnetic field, but still this effect was visible in our data. Nevertheless, the model fitting code was able to perform well since the code operated on binary images not the original tagged images. This enabled matching the faded tag lines to the calculated meshes. The ability to perform the fitting process seamlessly using images with relatively poor quality of tag lines is another advantage of our approach.
This study demonstrated that it is feasible to empirically model the spatial and temporal evolution of the LV wall motion using a two-parameter formulation in conjunction with tMRI-based visualization of the LV wall in the transverse planes at the apex, mid-ventricle and base levels. The algorithm developed and implemented in a simple fashion provides robust and accurate estimates of the two-parameters of the model when tMRI data was from normal human hearts. The performance of the estimator does not degrade with the fading of the tag lines. Once a model is constructed, a measure of interest characterizing certain aspect of the LV wall biomechanics, such as strain, strain rate or torsion, can be derived and expressed analytically. Any deviation in such measure from its normative value would indicate cardiac abnormality, and thereby serve as a surrogate biomarker for detecting cardiac dysfunction and evaluating its severity in clinical investigations or experimental studies with translational focus.
where r ′ = (x ′ 2 + y ′ 2)1/2. When the radial displacement is small (4α2/r ' 2 < < 1), the approximation (1 + √ (1 + 4α2/r ' 2))/2 ≃ 1 + α2/r ' 2 holds. But, the radial displacement remains large during most cardiac cycle especially near the end systole and 4α2/r ' 2 ∼ 1 holds.
This work was partially supported by the NIH Grant number P20 RR016475 from the INBRE Program of the National Center for Research Resources (IVS), the American Heart Association Scientist Development award (IVS) and the University of Malaya Research Grant award RG323-11HTM (MB).
- Frangi AF, Niessen WJ, Viergever MA: Three-dimensional modeling for functional analysis of cardiac images: a review. IEEE Trans Med Imaging 2001, 20: 2–25. 10.1109/42.906421View ArticleGoogle Scholar
- Luo G, Heng PA: LV shape and motion: B-spline-based deformable model and sequential motion decomposition. IEEE Trans Inf Technol Biomed 2005, 9: 430–446. 10.1109/TITB.2005.847508View ArticleGoogle Scholar
- Hashim S, Richens D: Finite element method in cardiac surgery. Interact Cardiovasc Thorac Surg 2006, 5: 5–8.View ArticleGoogle Scholar
- Kermani S, Moradi MH, Abrishami-Moghaddam H, Saneei H, Marashi MJ, Shahbazi-Gahrouei D: Quantitative analysis of left ventricular performance from sequences of cardiac magnetic resonance imaging using active mesh model. Comput Med Imaging Graph 2009, 33: 222–234. 10.1016/j.compmedimag.2008.12.005View ArticleGoogle Scholar
- Moroni D, Colantonio S, Salvetti O, Salvetti M: Heart deformation pattern analysis through shape modelling. Pattern Recognition and Image Analysis 2009, 19: 262–270. 10.1134/S1054661809020084View ArticleGoogle Scholar
- Bilgen M: Harmonic phase interference for the detection of tag line crossings and beyond in homogeneous strain analysis of cardiac tagged MRI data. Australas Phys Eng Sci Med 2010, 33: 357–366. 10.1007/s13246-010-0039-zView ArticleGoogle Scholar
- Sengupta PP, Tajik AJ, Chandrasekaran K, Khandheria BK: Twist mechanics of the left ventricle: principles and application. JACC Cardiovasc Imaging 2008, 1: 366–376. 10.1016/j.jcmg.2008.02.006View ArticleGoogle Scholar
- Alrefae T, Smirnova IV, Cook LT, Bilgen M: A model-based time-reversal of left ventricular motion improves cardiac motion analysis using tagged MRI data. Biomed Eng Online 2008, 7: 15. 10.1186/1475-925X-7-15View ArticleGoogle Scholar
- Zerhouni EA, Parish DM, Rogers WJ, Yang A, Shapiro EP: Human heart: tagging with MR imaging–a method for noninvasive assessment of myocardial motion. Radiology 1988, 169: 59–63.View ArticleGoogle Scholar
- Russel IK, Gotte MJ, Bronzwaer JG, Knaapen P, Paulus WJ, van Rossum AC: Left ventricular torsion: an expanding role in the analysis of myocardial dysfunction. JACC Cardiovasc Imaging 2009, 2: 648–655. 10.1016/j.jcmg.2009.03.001View ArticleGoogle Scholar
- Bower AF: Applied Mechanics of Solids. CRS Press, Boca Raton, FL; 2010.Google Scholar
- Akagawa E, Murata K, Tanaka N, Yamada H, Miura T, Kunichika H, Wada Y, Hadano Y, Tanaka T, Nose Y, et al.: Augmentation of left ventricular apical endocardial rotation with inotropic stimulation contributes to increased left ventricular torsion and radial strain in normal subjects: quantitative assessment utilizing a novel automated tissue tracking technique. Circ J 2007, 71: 661–668. 10.1253/circj.71.661View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.