A dual propagation contours technique for semi-automated assessment of systolic and diastolic cardiac function by CMR
Journal of Cardiovascular Magnetic Resonance volume 11, Article number: 30 (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 dual-contour propagation scheme shown in Figure 1 was used to propagate both ED and ES contours to all other time frames in the sequence. First, the NRR algorithm was used to propagate ED contours forward in time through systole and backward in time through diastole (white arrows in Figure 1). Next, ES contours were propagated forward in time through diastole and backward in time through systole (gray arrows in Figure 1). The contours propagated from ED and ES were then combined, as described in the Appendix, into a single set of endocardial and epicardial contours.
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 propagated contours and manual contours resulted in differences in VTCs. Figure 2 shows VTCs from a normal volunteer. ED-propagated contours with NRR resulted in volume overestimation near ES, and ES-propagated contours with NRR produced volume underestimation in early systole and late diastole. Dual-propagated contours showed excellent agreement throughout the entire cardiac cycle. Both ED and ES propagated contours using CAAS MRV underestimated the volumes as compared to the manually drawn, gold standard volumes, more than NRR propagated contours throughout the cardiac cycle; the CAAS MRV propagation method also changes the manually-drawn ED and ES contours slightly, so the volume difference is not zero at ED or ES in these curves.
Table 1 shows confidence intervals of the volume differences between each propagation method and manual. ED-propagated contours with NRR overestimate LV volume, whereas ES-propagated contours with NRR underestimate LV volume. In comparison, both ED-propagated and ES-propagated contours with CAAS MRV underestimate LV volume by a larger margin. However, the dual-propagated volumes were not statistically different from manually-contoured volumes.
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
Figure 3 shows the VTCs in normal volunteers derived from the manually-drawn contours and dual-propagated contours. The manual and dual-propagated VTCs were quite close to each other in all nine studies – particularly during systole and early diastole. This similarity between manual and propagated VTCs means that the contours manually drawn at ED and ES were consistently propagated to the other time frames in the cine sequence.
No statistically-significant differences were found between PER, ePFR and aPFR rates computed from manually-drawn contours and dual-propagated contours (Table 2). The correlation coefficients between the PER, ePFR and aPFR values were 0.92, 0.95 and 0.96 respectively (all P < 0.001). Figure 4 shows scatter and Bland-Altman plots comparing the manual and dual-propagated measurements of ejection and filling rates.
No significant difference was found between PER, ePFR and aPFR values computed from contours propagated with ED and ES contours drawn by two different users (User1 and User2). The differences (User2-User1) between PER, ePFR and aPFR values were 0.07 ± 0.16 EDV/s (P = 0.24), -0.03 ± 0.05 EDV/s (P = 0.11) and -0.01 ± 0.05 EDV/s (P = 0.50) respectively. The correlation coefficients for the PER, ePFR and aPFR values were 0.95 (P < 0.0001), 0.99 (P < 0.0001) and 0.99 (P < 0.0001) respectively. Figure 5 shows scatter and Bland-Altman plots comparing User1 and User2 measurements of PER, ePFR and aPFR.
LV Mass Evaluation
No significant difference was found between normalized LV masses (LV mass/ED LV mass) computed from dual-propagated contours and manual contours. The normalized LV mass difference (propagated-manual) was -0.015 ± 0.0077 (P = 0.08). Figure 6 shows the mean normalized LV mass throughout the cardiac cycle averaged over studies NV10-NV18. The normalized LV mass computed from dual-propagated contours was close to that from the manual contours and both remain stable throughout the cardiac cycle, indicating that the propagated contours were as stable as the manual contours.
Peak Ejection and Filling Rates in HTN
Figure 7 shows typical VTCs for a normal volunteer and a hypertensive patient measured from dual-propagated contours. Peak ejection rates are similar in both curves, but the early diastolic filling rate is lower in the hypertensive patient than in the normal while the late filling rate is higher.
Figure 8 shows the mean peak ejection and filling rates measured from dual-propagated contours in all 39 normals and 49 patients with hypertension. In hypertensives, PER was not different from normal (3.4 ± 0.1 vs. 3.2 ± 0.1 EDV/sec, P = NS). Diastolic filling rates, however, were altered compared to normals, demonstrating diastolic dysfunction in hypertension that is common in this patient group: ePFR was lower than normal (2.6 ± 0.1 vs. 3.2 ± 0.1 EDV/sec, P < 0.0001), but aPFR was higher than normal (2.4 ± 0.1 vs. 1.6 ± 0.1 EDV/sec, P < 0.0001).
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.
Next, a two-dimensional displacement field was defined on the template ROI. The deformation field was parameterized by a tensor product of quadratic, uniformly-spaced B-splines:
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.
The B-spline control points, μ, were computed to minimize the sum of squared differences between the pixel intensities in the template image, I t , and the source image, I s :
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
As described above, the NRR was used to propagate contours from both ED and ES to all other time frames in a sequence. These propagations resulted in two contours for each time frame (except at ED and ES). The two contours were combined into a single B-spline contour using a weighted-least-squares fit. The ED-propagated contour weight for a given time was computed using cubic-spline interpolation from the empirically-determined control points in Table 3. The end-systolic ES-propagated contour weight is one minus the ED contour weight. The ED-propagated and ES-propagated weights at a given frame are based on their distances from the ED and ES frames. For example, as the distance of a frame from ED increases, its ED-propagated weights decreases and its ES-propagated weights increases.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Ranganath S: Contour extraction from cardiac MRI studies using snakes. IEEE Trans Med Imaging. 1995, 14: 328-338. 10.1109/42.387714.
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.
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.
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.
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.
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.
Kybic J, Unser M: Fast parametric elastic image registration. IEEE Transactions on Image Processing. 2003, 12: 1427-1442. 10.1109/TIP.2003.813139.
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. 32
von Berg J, Lorenz C: A geometric model of the beating heart. Methods Inf Med. 2007, 46: 282-286.
Neter J, Kutner MH, Wasserman W, Nachtsheim CJ: Applied Linear Statistical Models. 1996, McGraw-Hill/Irwin, 4
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
This work was supported by NIH grant P50-HL077100.
The authors declare that they have no competing interests.
WF developed the dual-contour propagation algorithm, drew one set of contours for the inter-user variability study, and drafted the manuscript. HN performed manual contouring of normal datasets including one set of contours for the inter-user variability analysis. HG participated in the study design and helped draft the manuscript. SGL participated in the study design and helped draft the manuscript. IA participated in the study design and performed the statistical analysis. GJP helped draft the manuscript. DAC participated in study coordination. LJD participated in study coordination and helped draft the manuscript. TSD conceived of the study, participated in study design and coordination and helped draft the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Feng, W., Nagaraj, H., Gupta, H. et al. A dual propagation contours technique for semi-automated assessment of systolic and diastolic cardiac function by CMR. J Cardiovasc Magn Reson 11, 30 (2009). https://doi.org/10.1186/1532-429X-11-30