- Open Access
A dual propagation contours technique for semi-automated assessment of systolic and diastolic cardiac function by CMR
© Feng et al; licensee BioMed Central Ltd. 2009
Received: 26 March 2009
Accepted: 13 August 2009
Published: 13 August 2009
Although cardiovascular magnetic resonance (CMR) is frequently performed to measure accurate LV volumes and ejection fractions, LV volume-time curves (VTC) derived ejection and filling rates are not routinely calculated due to lack of robust LV segmentation techniques. VTC derived peak filling rates can be used to accurately assess LV diastolic function, an important clinical parameter. We developed a novel geometry-independent dual-contour propagation technique, making use of LV endocardial contours manually drawn at end systole and end diastole, to compute VTC and measured LV ejection and filling rates in hypertensive patients and normal volunteers.
39 normal volunteers and 49 hypertensive patients underwent CMR. LV contours were manually drawn on all time frames in 18 normal volunteers. The dual-contour propagation algorithm was used to propagate contours throughout the cardiac cycle. The results were compared to those obtained with single-contour propagation (using either end-diastolic or end-systolic contours) and commercially available software. We then used the dual-contour propagation technique to measure peak ejection rate (PER) and peak early diastolic and late diastolic filling rates (ePFR and aPFR) in all normal volunteers and hypertensive patients.
Compared to single-contour propagation methods and the commercial method, VTC by dual-contour propagation showed significantly better agreement with manually-derived VTC. Ejection and filling rates by dual-contour propagation agreed with manual (dual-contour – manual PER: -0.12 ± 0.08; ePFR: -0.07 ± 0.07; aPFR: 0.06 ± 0.03 EDV/s, all P = NS). However, the time for the manual method was ~4 hours per study versus ~7 minutes for dual-contour propagation. LV systolic function measured by LVEF and PER did not differ between normal volunteers and hypertensive patients. However, ePFR was lower in hypertensive patients vs. normal volunteers, while aPFR was higher, indicative of altered diastolic filling rates in hypertensive patients.
Dual-propagated contours can accurately measure both systolic and diastolic volumetric indices that can be applied in a routine clinical CMR environment. With dual-contour propagation, the user interaction that is routinely performed to measure LVEF is leveraged to obtain additional clinically relevant parameters.
Temporal changes in left ventricular (LV) volume over the cardiac cycle provides fundamental information regarding systolic and diastolic function of the heart but is difficult to measure by standard clinical techniques. Cine cardiovascular magnetic resonance (CMR) using serial short axis slices is well accepted as a gold standard for measuring geometry-independent ventricular volumes [1, 2]. Measurement of LV ED and ES volumes is based on drawing contours at the ED and ES time points. If contours could be reliably identified in all acquired time frames, ventricular volume-time curves (VTC) could be constructed, from which important parameters of ventricular function such as peak ejection rates (PER) and peak filling rates (PFR) [3, 4] can be derived. These parameters may complement flow indices  in the assessment of diastolic function.
Fully automated contouring techniques have been a research topic for many years [6–15], and, more recently, techniques have been developed for propagating contours drawn at single time frame to the remaining time frames [16–22]. While the accuracy of these methods continues to improve, contour review and editing by a trained expert is still mandatory. A common problem encountered in myocardial contour identification is the presence of papillary muscles; following the echocardiographic convention, papillary muscles are often excluded from the endocardial contour. At ED in the short axis CMR image, papillary muscles are usually not a problem because they are separated from the LV wall. During systole, however, the papillary muscles move close to the LV wall, and it can be difficult to distinguish papillary muscles from the heart wall without carefully examining the images. For this reason, fully automatic contouring routines often have difficulty detecting papillary muscles, and may include papillary muscle volume as part of the LV cavity volume in ED and as outside the LV cavity (i.e., in the myocardial muscle volume) in ES. This potentially affects the derived volumes and masses.
Consequently, we propose a semi-automated method which leverages the user interaction in drawing ED and ES contours by automatically propagating them to all other time frames in a typical cardiac scan. This dual-contour propagation technique has the potential to more accurately exclude papillary muscles from the LV wall than single-contour propagation techniques or fully automated techniques. The proposed dual-contour propagation technique will not require additional work by the user, because at most institutions contours are already routinely drawn at ED and ES to compute standard volumes, myocardial mass, and ejection fraction. The purpose of this study was to develop a novel semi-automated technique using dual-contour propagation to measure ventricular volumes throughout the cardiac cycle, and compare this method to manual and single-contour techniques in normal volunteers and hypertensive patients.
The study was approved by the appropriate institutional review boards and informed consent was obtained from all the participants. 39 normal human volunteers (NLs) and 49 hypertensive (HTN) patients consecutively enrolled in a study of resistant hypertension (defined as requiring 3 or more anti-hypertensive medications to achieve blood pressure < 140/90 mmHg) participated in this study. All patients were in sinus rhythm at the time of CMR.
CMR was performed on a 1.5-T scanner (CV/i, GE Healthcare, Milwaukee, WI) optimized for cardiac application. ECG-gated, breath-hold steady state free precision technique was used to obtain standard (2, 3 and 4 Chamber, Short Axis) views using the following parameters – slice thickness 8 mm with no gap between short-axis slices, field-of-view 44 × 44 cm, scan matrix 256 × 128, flip angle 45°, typical TR/TE = 3.8/1.6 ms; typical acquired temporal resolution approximately 40 ms); data reconstructed to 20 cardiac phases.
In all scans, LVED and LVES endocardial contours were manually drawn on all short axis slices between the mitral annulus and apex  with exclusion of the papillary muscles. These contours were then automatically propagated to all the other frames in the acquisition using the dual-contour propagation algorithm described below. For validation, LV contours were manually drawn on all time frames in 18 randomly-selected normal scans by a Level 3 trained CMR specialist. These contours were used as a gold standard for evaluating and validating the dual-contour propagation algorithm.
Non-rigid registration (NRR)  was used to propagate the contours manually drawn at end-diastole and end-systole to all other time frames in the acquisition. The NRR algorithm computed a deformation field that warped a source image to fit a template image. The deformation field was then used to propagate contours defined on the template image to the source image. Details of this algorithm are provided in the Appendix. All algorithms were implemented in MATLAB (The Mathworks, Natick, MA).
The LV volume at each time frame was computed by summing the volumes defined by the contours in each slice. The contour propagation procedure, however, propagated contours in all slices that were contoured at ED, and, near the base, the LV margin may have moved through the image plane in systole. To address this problem, the NRR algorithm was used in a long-axis slice to track a user-selected point near the mitral annulus through the image sequence. The displacement of this point was used to determine how much each short-axis slice should be included in the volume computation [4, 25]. For example, if the mitral annulus displaced 12 mm between ED and the current time frame and the slice thickness was 8 mm, the most basal slice would not be used in the volume computation and 50% of the second-most-basal slice volume would be used in the total volume. This mitral annulus tracking procedure was used to compute volumes for all the contouring techniques discussed in this paper, including manual contouring and automatic contour propagation algorithms.
Once the volumes were computed in each time frame, a VTC curve was constructed and differentiated with respect to time. End-diastole was defined as the maximum-volume time frame, and end-systole was defined as the minimum-volume time frame. Early diastole and late diastole were defined as the first and second halves respectively of the diastolic interval. The peak ejection rate (PER) was defined as the maximum negative time derivative during the systolic interval. The early diastolic and late diastolic peak filling rates (ePFR and aPFR) were defined as the maximum derivative during the early and late diastole.
Comparison Between Single and Dual-Contour Propagation
Most existing contour propagation techniques propagate contours from either ED or ES time frames [8, 16, 20–22, 26]. Volumes computed from dual-propagated contours, single-propagated contours from ED and ES using the NRR method described above, and single-propagated contours from ED and ES using CAAS MRV for Windows, version 3.2 (Pie Medical Imaging, Maastricht, the Netherlands), software were compared to volumes computed from manual contours on nine randomly-selected normal studies (identified by NV1-NV9). A VTC was computed for each type of contours for each study. To compare VTCs computed from different types of contours, differences were computed at each time point by subtracting the manual volume from the propagated volumes.
To assess inter-user variability in volumes computed from propagated contours, nine normal studies (NV10-NV18) were randomly selected. For each study, two sets of contours were manually drawn at ED and ES by different users. Each user was a Level 3 trained CMR specialist or equivalent. Each set of ED and ES contours was propagated using the dual contour technique, and VTCs and ejection/filling rates were computed.
LV Mass Evaluation
To evaluate stability of dual-propagated contours throughout the cardiac cycle, LV mass measurements were computed in each time frame for studies NV10-NV18 from dual-propagated contours and compared to that from one set of manually-drawn contours.
Comparison of PER and PFR Values in Normals and Hypertensives
The dual-contour propagation algorithm was used to propagate contours to all time frames and compute VTCs and ejection/filling rates in all 39 normals and 49 hypertensives.
Comparisons of LV volumes computed from different propagation schemes and comparisons of LV masses computed from dual-propagated and manual contours were performed using mixed modeling via PROC MIXED (SAS version 9.1). To account for the repeated measures within a subject, a compound symmetry correlation structure was assumed. In the LV volume study, confidence intervals on the differences based on the fitted mixed model were constructed each at 99% level to achieve a joint confidence level of at least 95% for this set of confidence intervals using Bonferroni adjustment .
Comparisons of PER and PFR values computed from dual-propagated and manually-drawn contours and comparisons of contours propagated by two different users were performed using two-tailed paired t-tests, correlation analysis, and Bland-Altman analysis. PER and PFR values derived from dual-propagated contours in hypertensive patients were compared to normals using unpaired t-tests. In all these statistical tests, a P-value less than 5% was considered statistically significant.
Comparison Between Single and Dual-Contour Propagation
Differences between LV volumes (expressed as fraction of EDV) computed from propagated contours and manually-drawn contours.
Volume Difference (EDV)
Mean ± SE
99% Confidence Interval
Dual NRR – Manual
-0.19 ± 0.56
ED NRR – Manual
1.61 ± 0.56
ES NRR – Manual
-3.50 ± 0.56
ED CAAS – Manual
-6.54 ± 0.55
ES CAAS – Manual
-11.05 ± 0.55
The average computation time for dual-contour propagation (not including manually contouring the ED and ES contours) was 7.3 minutes for a single study on a 2.6 GHz dual-core personal computer with 4 GB of RAM. Automated contour propagation using CAAS MRV required less than 1 minute per study. Manual contouring of all slices and phases (typically 12 to 14 short axis slices × 20 cardiac phases) required approximately 4 hours per study.
Validation of Functional Parameters
Differences between peak ejection and filling rates computed from manually-drawn and dual-propagated contours.
Rate Difference (EDV/s)
Mean ± SE
95% Confidence Interval
-0.12 ± 0.08
-0.07 ± 0.07
0.06 ± 0.03
LV Mass Evaluation
Peak Ejection and Filling Rates in HTN
In this paper, we described a novel dual-contour propagation technique for measuring volume-time curves (VTCs), validated it against manually drawn contours, and demonstrated its utility in a clinically-relevant patient population. This method requires nothing more than standard short-axis and long-axis CMR acquisitions and routinely drawn ED and ES contours. We show that the dual-propagated contours can be used to accurately measure peak ejection and filling rates compared to the reference standard of manually-drawn contours. The dual-contour propagation technique provides a fast, practical means of measuring volume-based indices of systolic and diastolic ventricular function from routine clinical CMR.
Several techniques have been proposed for propagating contours in CMR [16–22] and other modalities [28–33], but these techniques only propagate contours from a single time frame. While propagating contours from only one time frame requires less user interaction, we found that the resulting volumes are less accurate compared to dual propagation. If contours are only defined at ED, propagated contours with NRR may not be able to separate papillary muscles from the LV wall at ES resulting in the volume differences at ES (Figure 2). Propagating only ES contours with NRR may solve this problem, but volume differences occur in late diastole as demonstrated in Figure 2.
van Guens, et al.  proposed an automated method for drawing contours at ED and ES. The required user input was minimal – only four manually-drawn epicardial contours on two and four-chamber views at ED and ES – but volumes were only validated at ED and ES. In addition, for registration purposes, this method requires that both long axis and multiple short-axis acquisitions be performed with reproducible breath-hold positions, which can sometimes be difficult to obtain under clinical conditions. The contour propagation method proposed in this paper, however, does not have this limitation since contours are propagated in each slice independently.
Investigators have previously described use of volume time indices for measuring systolic and diastolic function [34, 35]. CMR allows measurement of ventricular volumes throughout the cardiac cycle independent of geometric assumptions. The excellent spatial resolution and image contrast make it potentially the most accurate clinically-applicable non-invasive technique for assessment of systolic and diastolic function. To provide an illustration of the utility of our propagation method in clinical assessment of patients with risk factors for heart failure, the dual-contour propagation technique was employed to assess the physiology of LV systolic and diastolic function in 49 patients consecutively enrolled in a study of resistant hypertension. The images in this study contain the normal range of image quality and presence of artifacts encountered under routine clinical conditions. The concentrically hypertrophied LV in the HTN patients had a normal LV ejection fraction and LV peak ejection rate; however, early peak filling was decreased and late filling rate was increased, consistent with diastolic dysfunction.
The inter-user and intra-user variability in the propagated contours depend on the inter-user and intra-user variability of the semi-automatically-drawn contours at ED and ES. This variability has been studied in [15, 36, 37]. In our study, no significant difference was found between PER, ePFR and aPFR values computed from contours propagated from ED and ES contours drawn by two different users.
A limitation of contour propagation algorithms in general is that any errors in the seed contours get propagated to all other time frames. Consequently, it is especially important to ensure accurate seed contours before propagation. Also in this paper, papillary muscle volume was considered part of the LV blood volume. Since the papillary muscle volume is relatively constant throughout the cardiac cycle, subtracting the papillary muscle volume would reduce the blood volume by the same amount in all phases and would not significantly affect filling or ejection rates, which are the key parameters determined in this work. Although not evaluated in the present study, we believe that another potential advantage of the dual-contour propagation approach is that defining contours at two time points provides increased robustness to imaging artifacts.
In conclusion, the dual-contour propagation technique provides a fast, accurate and practical means of measuring volume-based indices of systolic and diastolic ventricular function from routine clinical CMR.
First, a square region of interest (ROI) was defined that enclosed epicardial contours in the template image. Since the NRR was applied sequentially starting from either ED or ES, an epicardial contour (either manually-drawn or propagated) was always available in the template image.
where p is a point in the template ROI, β(p) is a B-spline basis function, μ i is a control point, C is the number of control points and v i is the knot location associated with the i- th control point. Eight control points were used in each dimension. The B-spline order and number of control points were determined empirically.
where Ω is the set of pixels in the template ROI.
The cost function in Eq. (1) was optimized using a multi-resolution strategy to speed up computation and avoid local minima in the cost function. First, the template ROI was resampled from 64 × 64 to 16 × 16 pixels and the number of control points was reduced to four in each dimension. The reduced resolution control points were then computed using Levenberg optimization algorithm with an analytical gradient and Hessian. This process was then repeated using template ROI resampled from 64 × 64 to 32 × 32 pixels. The number of control points was kept at four in each dimension, and the result of the coarser resolution optimization was used as the starting point for the optimization. Next, the displacement field was interpolated to eight control points in each dimension, and Eq. (1) was optimized on the full-resolution 64 × 64 template ROI. The source image was resampled accordingly in each multi-resolution layer to match the resolution of the template ROI. The result was a spatially continuous mapping of points from the template ROI to the source image. Finally, each contour point in the template ROI was mapped to the source image using the final displacement field.
Combining Contours Propagated From ED and ES
Control points for generating ED-propagated contour weights.
% Systolic Interval
ED – Propagated Contour Weight
% Diastolic Interval
This work was supported by NIH grant P50-HL077100.
- Bellenger NG, Davies LC, Francis JM, Coats AJ, Pennell DJ: Reduction in sample size for studies of remodeling in heart failure by the use of cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2000, 2: 271-278. 10.3109/10976640009148691.View ArticlePubMedGoogle Scholar
- Schalla S, Nagel E, Lehmkuhl H, Klein C, Bornstedt A, Schnackenburg B, Schneider U, Fleck E: Comparison of magnetic resonance real-time imaging of left ventricular function with conventional magnetic resonance imaging and echocardiography. Am J Cardiol. 2001, 87: 95-99. 10.1016/S0002-9149(00)01279-0.View ArticlePubMedGoogle Scholar
- Zeidan Z, Erbel R, Barkhausen J, Hunold P, Bartel T, Buck T: Analysis of global systolic and diastolic left ventricular performance using volume-time curves by real-time three-dimensional echocardiography. Journal of the American Society of Echocardiography. 2003, 16: 29-37. 10.1067/mje.2003.40.View ArticlePubMedGoogle Scholar
- Maceira AM, Prasad SK, Khan M, Pennell DJ: Normalized left ventricular systolic and diastolic function by steady state free precession cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2006, 8: 417-426. 10.1080/10976640600572889.View ArticlePubMedGoogle Scholar
- Rathi VK, Doyle M, Yamrozik J, Williams RB, Caruppannan K, Truman C, Vido D, Biederman RW: Routine evaluation of left ventricular diastolic function by cardiovascular magnetic resonance: a practical approach. J Cardiovasc Magn Reson. 2008, 10: 36-10.1186/1532-429X-10-36.PubMed CentralView ArticlePubMedGoogle Scholar
- Kaus MR, von Berg J, Weese J, Niessen W, Pekar V: Automated segmentation of the left ventricle in cardiac MRI. Med Image Anal. 2004, 8: 245-254. 10.1016/j.media.2004.06.015.View ArticlePubMedGoogle Scholar
- Kirschbaum SW, Baks T, Gronenschild EH, Aben JP, Weustink AC, Wielopolski PA, Krestin GP, de Feyter PJ, van Geuns RJ: Addition of the long-axis information to short-axis contours reduces interstudy variability of left-ventricular analysis in cardiac magnetic resonance studies. Invest Radiol. 2008, 43: 1-6. 10.1097/RLI.0b013e318154b1dc.View ArticlePubMedGoogle Scholar
- Lorenzo-Valdes M, Sanchez-Ortiz GI, Elkington AG, Mohiaddin RH, Rueckert D: Segmentation of 4D cardiac MR images using a probabilistic atlas and the EM algorithm. Med Image Anal. 2004, 8: 255-265. 10.1016/j.media.2004.06.005.View ArticlePubMedGoogle Scholar
- Lynch M, Ghita O, Whelan PF: Automatic segmentation of the left ventricle cavity and myocardium in MRI data. Comput Biol Med. 2006, 36: 389-407. 10.1016/j.compbiomed.2005.01.005.View ArticlePubMedGoogle Scholar
- Merrifield R, Keegan J, Firmin D, Yang GZ: Dual contrast TrueFISP imaging for left ventricular segmentation. Magn Reson Med. 2001, 46: 939-945. 10.1002/mrm.1280.View ArticlePubMedGoogle Scholar
- Mitchell SC, Bosch JG, Lelieveldt BP, van der Geest RJ, Reiber JH, Sonka M: 3-D active appearance models: segmentation of cardiac MR and ultrasound images. IEEE Trans Med Imaging. 2002, 21: 1167-1178. 10.1109/TMI.2002.804425.View ArticlePubMedGoogle Scholar
- Pluempitiwiriyawej C, Moura JM, Fellow , Wu YJ, Ho C: STACS: new active contour scheme for cardiac MR image segmentation. IEEE Trans Med Imaging. 2005, 24: 593-603. 10.1109/TMI.2005.843740.View ArticlePubMedGoogle Scholar
- Uzumcu M, van der Geest RJ, Sonka M, Lamb HJ, Reiber JH, Lelieveldt BP: Multiview active appearance models for simultaneous segmentation of cardiac 2- and 4-chamber long-axis magnetic resonance images. Invest Radiol. 2005, 40: 195-203. 10.1097/01.rli.0000154216.94316.f7.View ArticlePubMedGoogle Scholar
- van Assen HC, Danilouchkine MG, Frangi AF, Ordas S, Westenberg JJ, Reiber JH, Lelieveldt BP: SPASM: a 3D-ASM for segmentation of sparse and arbitrarily oriented cardiac MRI data. Med Image Anal. 2006, 10: 286-303. 10.1016/j.media.2005.12.001.View ArticlePubMedGoogle Scholar
- van Geuns RJ, Baks T, Gronenschild EH, Aben JP, Wielopolski PA, Cademartiri F, de Feyter PJ: Automatic quantitative left ventricular analysis of cine MR images by using three-dimensional information for contour detection. Radiology. 2006, 240: 215-221. 10.1148/radiol.2401050471.View ArticlePubMedGoogle Scholar
- Hautvast G, Lobregt S, Breeuwer M, Gerritsen F: Automatic contour propagation in cine cardiac magnetic resonance images. IEEE Trans Med Imaging. 2006, 25: 1472-1482. 10.1109/TMI.2006.882124.View ArticlePubMedGoogle Scholar
- Noble NM, Hill DL, Breeuwer M, Schnabel JA, Hawkes DJ, Gerritsen FA, Razavi R: Myocardial delineation via registration in a polar coordinate system. Acad Radiol. 2003, 10: 1349-1358. 10.1016/S1076-6332(03)00537-3.View ArticlePubMedGoogle Scholar
- Ranganath S: Contour extraction from cardiac MRI studies using snakes. IEEE Trans Med Imaging. 1995, 14: 328-338. 10.1109/42.387714.View ArticlePubMedGoogle Scholar
- Uzumcu M, van der Geest RJ, Swingen C, Reiber JH, Lelieveldt BP: Time continuous tracking and segmentation of cardiovascular magnetic resonance images using multidimensional dynamic programming. Invest Radiol. 2006, 41: 52-62. 10.1097/01.rli.0000194070.88432.24.View ArticlePubMedGoogle Scholar
- van der Geest RJ, Buller VG, Jansen E, Lamb HJ, Baur LH, van der Wall EE, de Roos A, Reiber JH: Comparison between manual and semiautomated analysis of left ventricular volume parameters from short-axis MR images. J Comput Assist Tomogr. 1997, 21: 756-765. 10.1097/00004728-199709000-00019.View ArticlePubMedGoogle Scholar
- van der Geest RJ, Lelieveldt BP, Angelie E, Danilouchkine M, Swingen C, Sonka M, Reiber JH: Evaluation of a new method for automated detection of left ventricular boundaries in time series of magnetic resonance images using an Active Appearance Motion Model. J Cardiovasc Magn Reson. 2004, 6: 609-617. 10.1081/JCMR-120038082.View ArticlePubMedGoogle Scholar
- Wierzbicki M, Drangova M, Guiraudon G, Peters T: Validation of dynamic heart models obtained using non-linear registration for virtual reality training, planning, and guidance of minimally invasive cardiac surgeries. Med Image Anal. 2004, 8: 387-401. 10.1016/j.media.2004.06.014.View ArticlePubMedGoogle Scholar
- Marcus JT, Gotte MJ, DeWaal LK, Stam MR, van der Geest RJ, Heethaar RM, Van Rossum AC: The influence of through-plane motion on left ventricular volumes measured by magnetic resonance imaging: implications for image acquisition and analysis. J Cardiovasc Magn Reson. 1999, 1: 1-6. 10.3109/10976649909080828.View ArticlePubMedGoogle Scholar
- Kybic J, Unser M: Fast parametric elastic image registration. IEEE Transactions on Image Processing. 2003, 12: 1427-1442. 10.1109/TIP.2003.813139.View ArticlePubMedGoogle Scholar
- Heiberg E, Wigström L, Carlsson M, Bolger AF, Karlsson M: Time Resolved Three-dimensional Automated Segmentation of the Left Ventricle. Proceedings of IEEE Computers in Cardiology. 2005, 599-602. 32Google Scholar
- von Berg J, Lorenz C: A geometric model of the beating heart. Methods Inf Med. 2007, 46: 282-286.PubMedGoogle Scholar
- Neter J, Kutner MH, Wasserman W, Nachtsheim CJ: Applied Linear Statistical Models. 1996, McGraw-Hill/Irwin, 4Google Scholar
- Chandrashekara R, Mohiaddin RH, Rueckert D: Analysis of 3-D myocardial motion in tagged MR images using nonrigid image registration. IEEE Trans Med Imaging. 2004, 23: 1245-1250. 10.1109/TMI.2004.834607.View ArticlePubMedGoogle Scholar
- Huang J, Abendschein D, Davila-Roman VG, Amini AA: Spatio-temporal tracking of myocardial deformations with a 4-D B-spline model from tagged MRI. IEEE Trans Med Imaging. 1999, 18: 957-972. 10.1109/42.811299.View ArticlePubMedGoogle Scholar
- Lorenz C, von Berg J: A comprehensive shape model of the heart. Med Image Anal. 2006, 10: 657-670. 10.1016/j.media.2006.03.004.View ArticlePubMedGoogle 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.847508.View ArticlePubMedGoogle Scholar
- Spottiswoode BS, Zhong X, Hess AT, Kramer CM, Meintjes EM, Mayosi BM, Epstein FH: Tracking myocardial motion from cine DENSE images using spatiotemporal phase unwrapping and temporal fitting. IEEE Trans Med Imaging. 2007, 26: 15-30. 10.1109/TMI.2006.884215.View ArticlePubMedGoogle Scholar
- Stegmann MB, Olafsdottir H, Larsson HB: Unsupervised motion-compensation of multi-slice cardiac perfusion MRI. Med Image Anal. 2005, 9: 394-410. 10.1016/j.media.2004.10.002.View ArticlePubMedGoogle Scholar
- Bonow RO, Bacharach SL, Green MV, Kent KM, Rosing DR, Lipson LC, Leon MB, Epstein SE: Impaired left ventricular diastolic filling in patients with coronary artery disease: assessment with radionuclide angiography. Circulation. 1981, 64: 315-323.View ArticlePubMedGoogle Scholar
- Harizi RC, Bianco JA, Alpert JS: Diastolic function of the heart in clinical cardiology. Arch Intern Med. 1988, 148: 99-109. 10.1001/archinte.148.1.99.View ArticlePubMedGoogle Scholar
- Karamitsos TD, Hudsmith LE, Selvanayagam JB, Neubauer S, Francis JM: Operator induced variability in left ventricular measurements with cardiovascular magnetic resonance is improved after training. J Cardiovasc Magn Reson. 2007, 9: 777-783. 10.1080/10976640701545073.View ArticlePubMedGoogle Scholar
- Steen H, Nasir K, Flynn E, El-Shehaby I, Lai S, Katus HA, Bluemcke D, Lima JA: Is magnetic resonance imaging the 'reference standard' for cardiac functional assessment? Factors influencing measurement of left ventricular mass and volumes. Clin Res Cardiol. 2007, 96: 743-751. 10.1007/s00392-007-0556-2.View ArticlePubMedGoogle 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.