|
|
||||||||
Technical Developments |
1 From the Departments of Anatomy with Radiology (A.A.Y.), Medicine (B.R.C., S.F.T.), and Engineering Science (W.J.H.), University of Auckland, 85 Park Road, Auckland, New Zealand; and the Department of Medicine, University of Alabama at Birmingham, Ala (L.J.D.). Received August 13, 1999; revision requested October 5; revision received November 12; accepted December 20. Supported in part by the Auckland Medical Research Foundation and the Health Research Council of New Zealand. Address correspondence to A.A.Y. (e-mail: a.young@auckland.ac.nz).
| ABSTRACT |
|---|
|
|
|---|
Index terms: Heart, volume, 51.92 Magnetic resonance (MR), phase imaging, 51.12144 Magnetic resonance (MR), rapid imaging, 51.121416 Magnetic resonance (MR), volume measurement, 51.12144
| INTRODUCTION |
|---|
|
|
|---|
A major limitation of the MR imaging section summation method is the prohibitive time necessary to outline the inner and outer boundaries of the LV in each section at ED and ES. This image contouring bottleneck precludes the use of multiplanar MR imaging in routine clinical care and limits its application to small research trials. Many semiautomated image segmentation algorithms have been applied to this problem (810), but these algorithms are not sufficiently robust for routine clinical use. Image pixel intensities are insufficient to adequately constrain the segmentation problem, owing to the limited temporal and spatial resolution, presence of image artifacts, and lack of contrast between blood and muscle. The endocardial trabeculae and papillary muscles also make the inner boundary difficult to define. At present, the amount of time spent on manual editing and correction renders automated methods almost as slow as manual contouring in clinical practice.
In this study, we evaluated a fast, accurate method of calculating LV mass and volumes from multiple MR sections, without the need for image contouring. The method relies on an accurate mathematic model of the LV that is fitted to a relatively small number of data points (guide points) provided by the user. No image processing is required and the method is therefore independent of image quality.
| Materials and Methods |
|---|
|
|
|---|
Volunteers.A group of 15 healthy subjects (10 men and five women; mean age, 23 years; age range, 1832 years) underwent MR imaging with a 1.5-T imager (Vision; Siemens, Erlangen, Germany). Prospectively gated cardiac cine images were acquired in eight to nine short-axis sections and three long-axis sections with use of a segmented k-space pulse sequence (repetition time msec/echo time msec of 8/5, flip angle of 10°, field of view of 280350 mm) with view sharing (1119 frames per section) (3,4). Each section was acquired during a breath hold of 1519 cardiac cycles. The short-axis sections spanned the heart from apex to base with a section thickness of 8.0 mm and section gap of 0.03.0 mm.
Phase-contrast velocity images (11) were acquired in a section positioned perpendicular to the ascending aorta approximately 20 mm above the aortic valve. Since the entire imaging protocol performed with the volunteers consisted of several hours of imaging (including tagging and many flow-quantification sequences), all subjects had a break of 12 hours between cine anatomic imaging (including the short-axis cine sections used in the volume calculations) and the phase-contrast flow imaging. At the latter, prospectively gated cine gradient-echo imaging was performed (24/6, flip angle of 30°, field of view of 250 mm, section thickness of 8.0 mm, velocity encoding of 150 or 250 cm/sec, and 2545 frames through the cardiac cycle). Velocity sensitization was encoded in the through-plane (section-select) direction.
Patients.MR imaging data from 13 patients, who underwent MR imaging as part of a previous study (12) at 2 weeks (n = 6) and 3 months (n = 7) after their first acute Q wave myocardial infarction, were selected for analysis of global function. The 13 studies were selected consecutively in alphabetical order from a database of 70 studies, details of which are given by Johnson et al (12). All MR studies were performed with a 1.5-T imager (Gyroscan; Philips, Shelton, Conn). Respiration-compensated gradient-echo cine images (30/14, flip angle of 40°, field of view of 300400 mm) were acquired in nine to 10 short-axis sections, each 8.0-mm thick with a 1.0-mm section gap.
Animal studies.MR data from eight dogs imaged 56 months after induction of mitral regurgitation, as part of a previous study (13) approved by the institutional animal care committee, were analyzed by means of guide-point modeling. Mitral regurgitation was induced by means of percutaneous chordal rupture of the mitral valvular apparatus, as described in reference 13. MR imaging was performed with a 1.5-T imager (Gyroscan) with the same gradient-echo cine imaging sequences as were used with the patients (six to eight short-axis sections, 8-mm-thick sections, 1.64.0 section gap, field of view of 350 mm). After imaging, the hearts were arrested with KCl and removed from the chest. After the atria and right ventricular free wall were removed from the interventricular septum and LV, the portions were weighed.
Image Analysis
Manual contouring and section summation.The inner and outer boundaries of the LV were defined on each image as follows. First, a semiautomatic region-growing method was used to locate preliminary boundaries (10). These boundaries were then manually corrected for each image. To the best ability of the observer (S.F.T.), papillary muscles were excluded from the myocardium and included in the blood pool. Great care was taken to identify the exact boundary of the myocardium. The contours were then reviewed by another observer (B.R.C.), and a consensus was achieved in all cases where there was a difference in opinion. Figure 1 shows short-axis sections at ED and ES together with the final boundaries (called "manual contours" in this article since extensive manual correction of the segmentation result was required). The areas outlined in each image were then multiplied by the section thickness plus the section gap and the resultant volumes summed to produce estimates of ED and ES volumes. LV mass was calculated from the mean of the ED and ES myocardial volumes (as the difference between the volumes enclosed by the epicardial and endocardial contours) multiplied by 1.05 g/mL.
|
Model Definition
A finite element model was used to represent the geometry of the LV. This model was similar to those used previously to describe LV deformation (14,15). The model consisted of 16 elements, each with cubic interpolation in the circumferential and longitudinal directions. Linear interpolation was used to couple the inner and outer surfaces into a coherent 3D model. The model was defined in a polar coordinate system (Appendix). This simplified the solution process, allowing the radial coordinate (
) of the model to be fitted as a function of the two angular coordinates (
and µ in the circumferential and longitudinal directions, respectively). Use of a polar coordinate system imposes two constraints on the shape of the model. First, there must be a straight "central axis" passing within the LV cavity from apex to base. Second, the location of the LV inner and outer surfaces must be uniquely specified on the basis of a radial distance from the central axis at all points. These constraints are not as severe as those typically applied in the case of echocardiography (16), and the piecewise cubic interpolation allows a variety of normal and deformed LV shapes to be accurately modeled, including LV aneurysms and severe hypertrophy (as in reference 14).
Initially, the model was scaled to each heart according to the distance between the base and the apex of the LV. The initial shape of the LV model was a regular ellipsoid, which was obtained by setting the inner and outer surfaces to a constant radial value. The extent of the model in the longitudinal (µ) direction was set to correspond to the most basal extent of the LV in the long-axis images at each time point.
Guide-Point Modeling
Guide points were interactively placed on the images by the user and the two-dimensional image coordinates converted to 3D coordinates with the section position information encoded in the image header. The model shape was fitted to the 3D guide-point positions by minimizing an error function consisting of the sum of a smoothing term and a term penalizing the distance between each guide point and the corresponding model position (Appendix). The smoothing term penalized changes in slope and curvature around the LV, allowing the model to realistically interpolate the very sparse guide-point data. The error function was quadratic in the model parameters, resulting in a simple linear least squares solution process. As the user placed or modified guide points, the model fit was updated in real time. The intersections of the model with the image planes were recalculated after each fit by means of a subdivision algorithm (17) and redisplayed on each image. The user continued to interactively add or modify guide points until the model-image intersections provided a good representation of the boundaries of the LV. Figure 2 shows the intersections of the model with the images shown in Figure 1.
|
Model Mass and Volume Comparisons
LV mass and volume could be estimated from the model in a number of ways. To assess the accuracy of the technique, we needed a consistent method of volume calculation that could be used for both manual and guide-point modeling methods; therefore, we used the intersections of the model with each image section plane as model contours at ED and ES (Fig 2). The model contours were then input to the same custom-written software used to calculate the volumes from the manual contours. In this way, the ability of the modeling process to reproduce the manual contours for the purposes of volume calculation could be directly assessed.
Statistical Tests
Volumes and mass were compared between the two methods by means of paired t tests. Differences with a P value of .05 were considered statistically significant.
| Results |
|---|
|
|
|---|
|
|
|
|
|
|
| Discussion |
|---|
|
|
|---|
The guide-point modeling method relies on the user to successively refine placement of the model on the images. No image processing is required. Current methods that rely on image processing lack the robustness and accuracy necessary for successful clinical application; however, local image information can be included in the model optimization (18).
LV absolute volumes tended to be overestimated with the model contours of operator 1 relative to operator 2 and the manual contours in the healthy group. These differences were spread evenly over all sections and could be equated with a small bias in the perceived position of the boundary. It is vital in any institution that a consensus be achieved among all operators in the placement of ventricular boundaries on the images. Small biases in contour placement can easily lead to large errors in volumes and mass. In this study, the SDs of the differences between manual and modeling methods were small, and the bias is unlikely to be clinically important.
Results of analysis of covariance comparing phase-contrast aortic flow measurements with manual and model estimates of stroke volume showed that there was no effect (other than heart rate) due to the measurement method. The comparison of LV mass with postmortem measurements showed an apparent overestimation of mass with the guide-point method, since papillary muscles were included in the postmortem weighing but were excluded from the guide-point model. The eight dogs were studied after 56 months of mitral regurgitation caused by rupture of the chordae and thus exhibited severe LV remodeling due to chronic volume overload (13). The contribution of the papillary muscles to LV mass was therefore likely to be smaller than usual, and some papillary muscle may have been included in the guide-point model owing to failure to distinguish it from LV wall. Also, errors due to trabeculation of the inner surface of the LV, possible inclusion of epicardial fat (together with the specific gravity factor of 1.05 g/mL used for comparison), and the slight overestimation of model mass relative to manual contours (Tables 2, 3) may have combined to produce the apparent overestimation.
The guide-point modeling method was found to reduce the time required to estimate LV mass, ED volume, ES volume, stroke volume, and ejection fraction by a factor of 510 without loss of accuracy. In conjunction with recent advances in fast cardiac imaging (3,4), a comprehensive cardiac examination including imaging (scout images plus eight to 12 cine sections in short- and long-axis orientations) and analysis (global and possibly regional function) would typically require less than half an hour to complete. This time saving enables high-spatial-resolution multiplanar MR studies of cardiac function to be efficiently used in routine clinical practice.
| APPENDIX |
|---|
|
|
|---|
n are the nodal values,
n are basis functions that give the relative weighting of each nodal value, and (
1,
2,
3) are the element coordinates. The geometric field
was defined to be the radial coordinate in a prolate spheroidal coordinate system:
, µ,
) are the radial, longitudinal, and circumferential coordinates of the polar system and (x, y, z) are the corresponding rectangular cartesian coordinates. The focal length f of the prolate system was chosen so that the
= 1 surface gave a good initial approximation of the LV epicardial surface. This provided an overall scale factor for each case. In each element, element coordinates
1,
2, and
3 were in the circumferential, longitudinal, and transmural directions, respectively. The basis functions were chosen to be bicubic Hermite in
1 and
2 and linear in
3 (19). Nodal values were shared between neighboring elements to give continuity in both slope and position throughout the LV. The model consisted of 16 3D elements: four elements in the longitudinal direction by four elements in the circumferential direction (Fig 3).
The model was fitted to guide-point data by minimizing the following error function:
g are the positions of the guide-point data and
(
g) are the model positions at element coordinates
g corresponding to
g. The element coordinates were found by projecting the guide data onto the model along lines of constant µ and
. S(
) denotes a smoothing term that was included in the error function to constrain the model to smoothly interpolate between the sparse guide points. As done previously (20), we used a weighted Sobolev norm that penalized displacement of the model from a prior shape:
-
*, where
* is the prior shape (in this article,
* = 1 [ie, a regular ellipsoid] is used for the prior shape). The weights
1 and
2 penalized the slope of the displacement field in the circumferential and longitudinal directions, respectively, whereas the weights ß1 and ß2 penalized curvature, and the weights
1,
2, and
3 coupled slopes between directions. The smoothing weights used in this article are
1 =
2 = 0, ß1 = ß2 = 0.001,
1 =
2 =
3 = 0.01. The issue of which weights are "optimal" for guide-point fitting remains an active area of research. However, the optimality criterion is difficult to quantify since the user always places the guide points in such a way as to adequately model the ventricle. A wide range of weights were tried in our experiments, and these values provided good results in all cases.
The error function (Eq [A3]) was quadratic in the unknown nodal parameters (
n), allowing direct solution by means of the standard least squares methods.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
Author contributions: Guarantor of integrity of entire study, A.A.Y.; study concepts, A.A.Y., B.R.C.; study design, A.A.Y.; definition of intellectual content, A.A.Y., B.R.C.; literature research, A.A.Y.; clinical studies, B.R.C., L.J.D.; experimental studies, A.A.Y., S.F.T., B.R.C., L.J.D.; data acquisition, A.A.Y., S.F.T., B.R.C., L.J.D.; data analysis, A.A.Y., S.F.T., W.J.H.; statistical analysis, A.A.Y.; manuscript preparation, A.A.Y.; manuscript editing and review, A.A.Y., S.F.T., W.J.H., B.R.C.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. A. Young and A. F. Frangi Computational cardiac atlases: from patient to population and back Exp Physiol, May 1, 2009; 94(5): 578 - 596. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. A.H. Stewart, O. C. Raffel, A. J. Kerr, R. Gabriel, I. Zeng, A. A. Young, and B. R. Cowan Pilot Study to Assess the Influence of {beta}-Blockade on Mitral Regurgitant Volume and Left Ventricular Work in Degenerative Mitral Valve Disease Circulation, September 2, 2008; 118(10): 1041 - 1046. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Chan, F. Khafagi, A. A. Young, B. R. Cowan, C. Thompson, and T. H. Marwick Impact of coronary revascularization and transmural extent of scar on regional left ventricular remodelling Eur. Heart J., July 1, 2008; 29(13): 1608 - 1617. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. P. Finn, K. Nael, V. Deshpande, O. Ratib, and G. Laub Cardiac MR Imaging: State of the Technology. Radiology, November 1, 2006; 241(2): 338 - 354. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. J. M. van Geuns, T. Baks, E. H. B. M. Gronenschild, J.-P. M. M. Aben, P. A. Wielopolski, F. Cademartiri, and P. J. de Feyter Automatic Quantitative Left Ventricular Analysis of Cine MR Images by Using Three-dimensional Information for Contour Detection. Radiology, July 1, 2006; 240(1): 215 - 221. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. J.S. Cooper, A. R.J. Phillips, S. Y. Choong, B. L. Leonard, D. J. Crossman, D. H. Brunton, E. L. Saafi, A. M. Dissanayake, B. R. Cowan, A. A. Young, et al. Regeneration of the Heart in Diabetes by Selective Copper Chelation Diabetes, September 1, 2004; 53(9): 2501 - 2508. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. J. Francois, D. S. Fieno, S. M. Shors, and J. P. Finn Left Ventricular Mass: Manual and Automatic Segmentation of True FISP and FLASH Cine MR Images in Dogs and Pigs Radiology, February 1, 2004; 230(2): 389 - 395. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. G. Fonseca, H. C. Oxenham, B. R. Cowan, C. J. Occleshaw, and A. A. Young Aging alters patterns of regional nonuniformity in LV strain relaxation: a 3-D MR tissue tagging study Am J Physiol Heart Circ Physiol, July 11, 2003; 285(2): H621 - H630. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| RADIOLOGY | RADIOGRAPHICS | RSNA JOURNALS ONLINE |