This Article Abstract Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Request Permissions Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Garot, J. Articles by Lima, J. A. C. Search for Related Content PubMed PubMed Citation Articles by Garot, J. Articles by Lima, J. A. C. Pubmed/NCBI databases Compound via MeSH Substance via MeSH Medline Plus Health Information MRI Scans Related Collections Cardiovascular imaging agents/Techniques CT and MRI

(Circulation. 2000;101:981.)
© 2000 American Heart Association, Inc.

Clinical Investigation and Reports

Fast Determination of Regional Myocardial Strain Fields From Tagged Cardiac Images Using Harmonic Phase MRI

Jérôme Garot, MD; David A. Bluemke, MD, PhD; Nael F. Osman, PhD; Carlos E. Rochitte, MD; Elliot R. McVeigh, PhD; Elias A. Zerhouni, MD; Jerry L. Prince, PhD; João A. C. Lima, MD

From the Cardiology Division of the Department of Medicine (J.G., C.E.R., J.A.C.L.) and the Departments of Radiology (D.A.B., E.A.Z.) and Biomedical Engineering (E.R.M.) of the Johns Hopkins University School of Medicine, and the Department of Electrical and Computer Engineering of the Whiting School of Engineering (N.F.O., J.L.P.), Johns Hopkins University, Baltimore, Md.

Correspondence to João A.C. Lima, MD, Cardiology Division, Blalock 569, Johns Hopkins Hospital, 600 N Wolfe St, Baltimore, MD 21287-6568. E-mail jlima{at}welchlink.welch.jhu.edu

	Abstract

Top Abstract Introduction Methods Results Discussion Appendix A Appendix B References

 
Background—Tagged MRI of the heart is difficult to implement clinically because of the lack of fast analytical techniques. We investigated the accuracy of harmonic phase (HARP) imaging for rapid quantification of myocardial strains and for detailed analysis of left ventricular (LV) function during dobutamine stimulation.

Methods and Results—Tagged MRI was performed in 10 volunteers at rest and during 5 to 20 µg^-1 · kg^-1 · min^-1 dobutamine and in 9 postinfarct patients at rest. We compared 2D myocardial strains (circumferential shortening, Ecc; maximal shortening, E₂; and E₂, direction) as assessed by a conventional technique and by HARP. Full quantitative analysis of the data was 10 times faster with HARP. For pooled data, the regression coefficient was r=0.93 for each strain (P<0.001). In volunteers, Ecc and E₂ were greater in the free wall than in the septum (P<0.01), but recruitable myocardial strain at peak dobutamine was greater in the LV septum (P<0.01). E₂ orientation shifted away from the circumferential direction at peak dobutamine (P<0.01). HARP accurately detected subtle changes in myocardial strain fields under increasing doses of dobutamine. In patients, HARP-determined Ecc and E₂ values were dramatically reduced in the asynergic segments as compared with remote (P<0.001), and E₂ direction shifted away from the circumferential direction (P<0.001).

Conclusions—HARP MRI provides fast, accurate assessment of myocardial strains from tagged MR images in normal subjects and in patients with coronary artery disease with wall motion abnormalities. HARP correctly indexes dobutamine-induced changes in strains and has the potential for on-line quantitative monitoring of LV function during stress testing.

Key Words: magnetic resonance imaging • ventricles • myocardium • contractility

	Introduction

Top Abstract Introduction Methods Results Discussion Appendix A Appendix B References

 
The analyses of myocardial wall motion abnormalities with dobutamine stress echocardiography¹ and MRI² are established methods for the detection of myocardial ischemia. However, the assessment of regional left ventricular (LV) function by either method is semiquantitative and subjective.^{2 3} Tagged MRI of the heart is a valuable technique for the quantitative noninvasive assessment of regional myocardial contractile performance.^{4 5 6} The spatial modulation of magnetization (SPAMM) technique generates 2 orthogonal sets of parallel planes of magnetic saturation by a sequence of nonselective radiofrequency pulses.^{7 8} The myocardium appears with a spatially encoded pattern that moves with the tissue and can be analyzed to reconstruct myocardial motion. Several methods for reconstructing LV strain fields have been proposed.^{5 6 9 10} The position of the tags must be measured with the use of a tag detection algorithm¹¹ and myocardial motion computed from displacement information from each tag.^{12 13} However, the process is time consuming because it requires the tracking of tag lines and myocardial contours with an interactive computer-aided technique.^{11 13} The lack of fast analytical techniques of myocardial strain quantification represents the main limitation to routine clinical utilization of tagged cardiac MRI.

SPAMM-tagged MR images correspond to a collection of spectral peaks in the Fourier domain.^{7 8} The inverse Fourier transform of one of these peaks is a complex image whose phase is linearly related to a directional component of the tissue displacement. A harmonic phase (HARP) image is the calculated phase of this complex image, which can be used to synthesize conventional tag lines and calculate 2D myocardial strain.¹⁴ This HARP imaging approach¹⁴ allows fast visualization and automated analysis of tagged cardiac MR images. Its potential in clinical cardiology depends on the demonstration of its sensitivity to small changes in myocardial strain during pharmacological stress testing and on its ability to accurately index regional wall motion abnormalities. We investigated the accuracy of HARP for quantitative assessment of 2D myocardial strain fields in normal individuals who underwent dobutamine stress–tagged MRI and in postinfarct patients with underlying abnormal regional LV function.

	Methods

Top Abstract Introduction Methods Results Discussion Appendix A Appendix B References

 
Study Population
SPAMM-tagged MR images were obtained in the short-axis view under a dobutamine challenge in 10 healthy male volunteers (age 37±10 years). All volunteers had no history of chest pain, prior cardiovascular disease, or cardiovascular risk factors. We also studied 9 patients at baseline (8 men, age 43±11 years) 3±2 days after a first acute myocardial infarction (AMI, 4 anterior, 5 inferior). Informed written consent was obtained according to the standards of the Joint Committee for Clinical Investigation of the Johns Hopkins Hospital.

MRI Protocol
Images were acquired during multiple breath-holds on a 1.5-T whole-body magnet (Signa, General Electric Medical Systems). Anterior and posterior phased-array coils were used for signal acquisition. We used an ECG-triggered segmented k-space fast gradient-echo imaging pulse sequence.¹⁵ The tagging pulse sequence consisted of nonselective radiofrequency pulses separated by spatial modulation of magnetization-encoding gradients to achieve tag separation of 7 mm. After scout images were completed, contiguous stacks of 4 base-to-apex short-axis cross sections were prescribed. Two sets of identical short-axis views were acquired (the second set rotated by 90°). A slab saturation band was applied to presaturate the blood in the LV, resulting in "black blood" in the LV cavity. This imaging sequence allowed us to image 4 slices within 4 breath-holds (14 to 20 seconds each). The number of views per phase was decreased as heart rate increased to maximize temporal resolution. Scanner settings were field of view 36 cm, tag separation 7 mm, slice thickness 8 mm, TR 6.5 ms, TE 2.3 ms, tip angle 15°, and image matrix 256x160, with 5 to 7 phase-encoded views per movie frame and cardiac cycle.

Dobutamine Stress–MR study
A single-lead ECG was continuously monitored and blood pressure was recorded at baseline and every 3 minutes throughout the procedure. After baseline acquisitions, dobutamine was infused through a digital infusion pump at 5 and 20 µg^-1 · kg^-1 · min^-1. Imaging began 2 minutes after each dose increase and required 3 minutes for each of the 4 levels. Criteria for terminating the study were (1) acceleration of heart rate >100% of age-predicted maximal heart rate, (2) fall of systolic blood pressure >30 mm Hg, (3) chest pain compatible with angina, (4) frequent ventricular or supraventricular ectopic beats, (5) intolerable side effects of dobutamine.

Harmonic Phase Imaging
The SPAMM technique uses a special pulse sequence to spatially modulate the longitudinal magnetization of the myocardium before acquiring image data.^{7 8} SPAMM-tagged images have regularly distributed spectral peaks in k-space, and each peak contains information about tag motion in a given direction. HARP imaging is based on the use of isolated peaks extracted with a bandpass filter (Figure 1) (Appendix A). One spectral peak was extracted for each direction of tag lines.¹⁴ The inverse Fourier transform of one of these peaks is a complex image whose phase is linearly related to a directional component of the true motion. The principal value of the phase was used to construct a HARP image (Figure 1) (Appendix A), which is linearly related to a component of the 3D motion except that it is constrained to lie in the range (-,+). Slopes of phases reflect the frequency of the tag pattern and phase images reflect motion of the heart. HARP images can be used to measure 2D strains^{14 16} (Appendix B), described as normalized myocardial deformation in a specific direction given by a unit vector. Figure 2 shows examples of 2D strain maps of the LV short-axis view throughout systole in a healthy volunteer at 5 and 20 µg^-1 · kg^-1 · min^-1 dobutamine.

View larger version (58K): [in this window] [in a new window]   Figure 1. SPAMM-tagged images of the LV short-axis view (top), the magnitudes of the corresponding Fourier transform (middle), and the corresponding HARP images (bottom). Circles in the Fourier domain correspond to the size of the bandpass filter.

View larger version (79K): [in this window] [in a new window]   Figure 2. SPAMM-tagged images of the LV short-axis view throughout systole in a healthy volunteer at baseline (top) and the corresponding 2D strain maps by HARP (circumferential shortening) with 5 (middle) and 20 µg-1 · kg-1 · min-1 (bottom) dobutamine. Black indicates shortening; yellow indicates stretching. HARP allowed strain to be calculated at any arbitrary pixel within the myocardium.

Strain Computation by Operator-Controlled Segmentation
Myocardial strains were assessed off-line by means of an established technique¹¹ and HARP by 2 independent observers. Coordinates of the posterior right ventricular–LV insertion point were calculated on the most basal slice and were used as reference landmarks for segmentation of the LV. The myocardium was divided into 5 segments as follows: segment 1, inferior; 2, inferolateral; 3, lateral; 4, anterior; and 5, septal wall. For the conventional analysis, images were processed with the use of an in-house–developed software program (Findtags). This method requires interactive and time-consuming detection of myocardial contours and tag lines and generates a detailed motion map through the use of interpolation.¹¹ Detection of myocardial contours or tag lines was not necessary with HARP. Strain was defined as the deformation gradient and was related to the derivative of the phase, that is, the local frequency. Eulerian strain was calculated by means of the single-shot harmonic phase (SHARP) approach.¹⁴ Three circles then were superimposed on the first image from the subendocardium, throughout the myocardium, to the subepicardium. The inner circle was placed as close as possible to the endocardium while avoiding papillary muscles. The outer circle was located as close as possible to the epicardium and the mid-circle, equally spaced between the inner and outer ones (Figure 3). For a single slice, 45 points (5 segmentsx3 layersx3 points per segment) were automatically tracked throughout systole by the cine-HARP approach.¹⁶ Strain changes were assessed between the reference and the deformed state (end-systole) by the fractional changes in length in the circumferential direction (Ecc) in each myocardial layer. A negative value stands for compression of a line segment between 2 material points. In addition, we calculated maximal shortening (E₂) as a principal strain given by the magnitude as well as the direction of the associated eigenvector.

View larger version (62K): [in this window] [in a new window]   Figure 3. Segmentation of the myocardium by the use of 3 equally spaced circles (each divided into 5 segments) throughout the wall. Top, a sequence of LV tagged images throughout systole is displayed in a patient with a recent anterior AMI (arrows). The corresponding 2D map of circumferential shortening is shown at the bottom. RV indicates right ventricle; LV, left ventricle; s, septum; i, inferior; il, inferolateral; l, lateral; and a, anterior wall.

Myocardial wall thickening (MWT) measured from cine-MR images was used as a gold standard for the detection of regional wall motion abnormality in patients with coronary artery disease (CAD). It was calculated at baseline in postinfarct patients and volunteers (control group), in the entire LV, by use of the automatic detection of endocardial and epicardial boundaries (Findtags),¹¹ as

where ESWT is end-systolic wall thickness and EDWT end-diastolic wall thickness. Wall motion abnormality (dysfunctional segments) was defined as a value of MWT <2 SD of the mean value in the control population and akinetic segments as segments with a MWT <5%.

Statistical Analysis
Mean values are expressed as mean±SD. Comparisons between both methods were assessed by linear regression analysis and Bland-Altman plots. Mean values of myocardial strains in each myocardial layer were compared by paired Student’s t test. Comparisons of myocardial strains under different doses of dobutamine were assessed by repeated-measures ANOVA. Percent increase in strain under dobutamine was compared by ² analysis. Interobserver reproducibility of measurements by HARP was assessed in 4 randomly selected volunteers and 4 patients by the use of a linear regression analysis and by calculating the mean difference between both series of measurements. All tests were 2-tailed, and a value of P<0.05 was considered statistically significant.

	Results

Top Abstract Introduction Methods Results Discussion Appendix A Appendix B References

 
Myocardial Strain by HARP and Conventional Method
Once the filter specifications are defined for each series of images, full quantitative assessment of myocardial strains in a single slice took <3 minutes by HARP. Quantitative analysis of a complete data set (baseline and 2 steps of dobutamine) required 10 hours by the conventional approach and 60 minutes by HARP. For pooled data in normal and dysfunctional myocardium, HARP led to reproducible results between 2 independent observers for Ecc, E₂, and E₂ direction. Regression coefficients between both analyses were r=0.98, 0.99, and 0.97, respectively (P<0.0001). Mean differences between both analyses were 5.4x10^-4±3.5x10^-3 (<2%), 1.8x10^-3±9.6x10^-3 (<2%), and 1.2±2.7° (<4%), respectively (NS).

For pooled data from volunteers and patients, the 2 methods showed good correlation. For Ecc, E₂, and E₂ direction, regression coefficients were r=0.93 (P<0.001) (Figure 4). Comparisons between data are displayed in Bland-Altman plots in Figure 5. In volunteers, the mean differences between both techniques were 1.1x10^-3±2.6x10^-2 (0.5%, NS), 3.6x10^-3±2.1x10^-2 (1.6%, P<0.01), and 1.1±4.1° (7%, P<0.01), for each strain, respectively. HARP led to a slight underestimation of each myocardial strain in the subendocardium (-0.231±0.042 vs -0.236±0.050, P=0.008, -0.256±0.049 vs -0.262±0.056, P<0.0001, and 15.6±6.4 vs 17.3±6.6°, P<0.0001, respectively). In patients with CAD, the mean differences between both methods were also very small and statistically not significant (5.1x10^-3±3.4x10^-2, 4.1x10^-3±3.5x10^-2, and 0.2±6.2°, respectively).

View larger version (19K): [in this window] [in a new window]   Figure 4. Correlation between myocardial strain as assessed by HARP (y-axis) and the conventional method (x-axis) for Ecc (A), E2 (B), and E2 direction (C). Black symbols indicate volunteers; white symbols, patients with CAD.

View larger version (25K): [in this window] [in a new window]   Figure 5. Bland-Altman plots for the comparison of myocardial circumferential shortening (A), maximal shortening (B), and the direction of the eigenvector (C). Black symbols indicate volunteers; white symbols, patients with CAD.

Myocardial Strain Changes Induced by Dobutamine
The mean time in the magnet was 30 minutes. We did not observe any significant side effect during dobutamine infusion. The evolution of 2D strains under dobutamine is displayed in Figure 6. For each strain, both techniques were able to detect a transmural strain gradient. Recruitable myocardial strain (% increase at peak dobutamine) as assessed by HARP increased similarly in subendocardium and subepicardium (22.4% vs 21.6% for Ecc and 27.2% vs 26.5% for E₂, P=NS). The 2 methods were able to depict subtle changes in myocardial strains during dobutamine stimulation (5 µg^-1 · kg^-1 · min^-1 vs baseline, P<0.01, and 20 vs 5 µg^-1 · kg^-1 · min^-1, P<0.001). E₂ orientation shifted further from the circumferential direction at 20 µg^-1 · kg^-1 · min^-1 as compared with baseline (18±6° vs 12±5°, P<0.01) (Figure 6). At baseline, Ecc and E₂ were greater in LV free wall (ie, lateral and posterolateral wall) than in the septum (P<0.01) (Table). However, recruitable strain at peak dobutamine was greater in LV septum than in free wall (24.1% vs 14.6% and 27.6% vs 16.1%, for Ecc and E₂, respectively, P<0.01). At baseline, Ecc and E₂ were greater at the apex compared with the base (P<0.01) (Table), but recruitable deformation was not significantly different at the apex versus the base (21.2% vs 19.8% and 21.2% vs 19.7%, respectively, NS).

View larger version (32K): [in this window] [in a new window]   Figure 6. Evolution of circumferential shortening (A), maximal shortening (B), and the direction of the eigenvector (C) as assessed by the conventional technique (gray bars) and by HARP (white bars) in the subendocardium (solid bars), midwall (striped bars), and subepicardium (hatched bars) at baseline and under increasing doses of dobutamine. *P<0.01 vs baseline, P<0.001 vs 5 µg-1 · kg-1 · min-1 dobutamine. P<0.05 and §P<0.01 vs conventional technique.

View this table: [in this window] [in a new window]   Table 1. Dobutamine-Induced Modifications of Circumferential (Ecc) and Maximal Shortening (E2) Analyzed Segment by Segment and at Each Level in the Left Ventricle

Myocardial Strain in Patients With CAD
Myocardial strain was measured in 9 patients (8 men, 43±11 years) at baseline, 3±2 days after a first AMI (4 anterior, 5 inferior), and compared with MWT obtained from cine-MR images. Of the 180 segments analyzed (5 segmentsx4 slicesx9 patients), 98 were classified as dysfunctional according to wall thickening analysis (ie, MWT <2 SD). Of these, 46 were classified as akinetic (MWT <5%). Figure 3 shows an example of a strain map obtained by HARP in a patient with anterior AMI. In each layer, Ecc and E₂ were decreased in dysfunctional myocardial segments when compared with remote (P<0.001) and further decreased in akinetic segments (P<0.01 vs dysfunctional) (Figure 7). E₂ direction was increased in dysfunctional segments when compared with remote (34±19 vs 19±10°, in the subendocardium, P<0.001) and was further augmented in akinetic segments (41±21°, P<0.01 vs dysfunctional).

View larger version (19K): [in this window] [in a new window]   Figure 7. Ecc (A), E2 (B), and E2 direction (C) in dysfunctional (hypo+akinetic; gray bars) and akinetic segments (white bars) in patients with CAD. *P<0.001 vs remote (black bars), P<0.01 vs dysfunctional.

	Discussion

Top Abstract Introduction Methods Results Discussion Appendix A Appendix B References

 
We describe and validate a new image processing method (HARP) based on the use of isolated spectral peaks in SPAMM-tagged MR images, which allows for rapid analysis of myocardial strain from tagged cardiac images. We also use MRI tagging to measure the effect of graded dobutamine stimulation on regional myocardial mechanics in the normal human LV. HARP is based on the fact that the inverse Fourier transform of a spectral peak in the Fourier domain is a complex image with phase linearly related to a directional component of tissue displacement. This permits the isolation of tag motion components, making the analysis of regional function completely automatic. We demonstrate that measurements of 2D myocardial strains by HARP are reproducible and similar to those obtained by a conventional tag motion tracking technique, both in patients with CAD with wall motion abnormalities and in normal volunteers at rest and during inotropic stimulation.

Dobutamine-Induced Myocardial Strain Alterations in the Normal Human Heart
Previous studies using tagged MRI have reported normal values and regional variations in 2D systolic strains at rest.^{17 18} We compared myocardial deformation at 20 µg^-1 · kg^-1 · min^-1 versus baseline and found a significant increase in both circumferential and maximal shortening. Other studies focusing on the effect of dobutamine on regional LV function as assessed by tagged MRI in normal subjects have shown an increase in myocardial strain from baseline to 10 µg^-1 · kg^-1 · min^-1 and no dobutamine-induced wall motion heterogeneity.^{19 20} Our data confirm that dobutamine does not further exacerbate variations in regional wall motion contractility documented at rest. Although Ecc and E₂ were less in the anteroseptal than in the LV free wall at baseline, percent increase in myocardial strain under dobutamine was greater in the anteroseptal than in the free wall, resulting in a relative homogeneity of myocardial contraction under inotropic stimulation. When analyzed by cross-sectional level along the LV long axis, dobutamine infusion resulted in a uniform increase in myocardial strains, in agreement with recently reported data.²⁰

Because of wall shear (differences between subepicardial and subendocardial rotational deformation), maximal shortening does not occur in the same direction as circumferential shortening (Figure 8). We report a slight increase in the angle between Ecc and E₂ orientations as a result of dobutamine stimulation, indicating that dobutamine slightly amplifies myocardial wall shear. In agreement with previous studies,²⁰ the angle remained <20°, indicating that maximal shortening still remained circumferentially oriented. Previous work has shown that dobutamine leads to increased rotation of both subepicardial and subendocardial layers in the normal left ventricle.²¹ Because rotation increases in both layers, changes in the direction of the eigenvector are blunted, maintaining LV efficiency.

View larger version (15K): [in this window] [in a new window]   Figure 8. Schematic representation of LV short-axis slice. Systolic deformation of a triangle and a circle located within the myocardium is shown. The direction of the normal circumferential strain (Ecc) is parallel to the tangent of the myocardium in respect to the heart surface. The direction of the normal radial strain (Err) is perpendicular to the circumferential direction, toward the ventricular centroid. Because of wall shear, a circle would be deformed into an ellipse during systole. Principal strains are given by the magnitude and direction of the eigenvector. The major and minor axes of the ellipse are representations of the minimal principal strain (maximal shortening, E2) and the greatest principal strain (maximal elongation, E1). In addition to the magnitude of E2 and Ecc, we assessed angle between the direction of Ecc and that of E2.

Potential Advantages of HARP
HARP can be used with any tagging technique provided the tag pattern is planar and tag lines are uniformly apart from each other. We assessed myocardial strain off-line from the tagged-image data set. We defined the filter specifications for each series of images, which represents the most time-consuming part before the HARP analysis (30 to 40 minutes). Once the filter is set, a full quantitative analysis of the data typically takes <3 minutes. Standardized acquisitions of tagged images allow subsequent presetting of the filter and can provide very fast myocardial strain mapping. Furthermore, HARP images can be extracted directly from the raw data, allowing very fast display of 2D strain fields. In other words, HARP may provide on-line detailed quantitative assessment of 2D myocardial strains.

Future Clinical Applications
HARP may have important clinical applications by overcoming the main limitation to routine clinical utilization of tagged cardiac MRI. When used in combination with dobutamine-tagged MRI, HARP might be of great value for the detection of myocardial ischemia during stress testing. Similarly, it might provide fast and accurate quantification of functional recovery in stunned or hibernating but viable myocardium.^{22 23} It also might be useful for studying dynamic changes in regional LV function after acute infarction by allowing serial quantitative examinations over time.^{24 25} Alterations in LV torsional deformation may be important in several pathological states.²⁶ Because HARP has the potential for other applications of any tracking motion technique, it allows for rapid noninvasive assessment of twist mechanics in the human heart, in different myocardial segments, and at each myocardial layer.

Study Limitations
We measured 2D myocardial strain fields and did not compensate for through-plane translation of the heart.^{27 28} It is known that 3D deformation allows more accurate evaluation of cardiac mechanics.²⁹ However, 2D analysis of regional myocardial function by tissue tagging is sufficiently powerful surements of small planar displacements by HARP are readily extendable to 3D motion by use of the out-of-plane tag direction for quantitative assessment of a 3D strain tensor.^{14 16}

Conclusions
HARP MRI provides fast, accurate measurements of 2D strain fields from tagged MR images in normal individuals and in patients with wall motion abnormalities caused by CAD. It allows for detailed analysis of the effect of graded dobutamine stimulation on regional myocardial mechanics in the normal human LV. Finally, with standardized acquisitions, HARP has the potential for providing on-line quantitative monitoring of LV function during stress testing in humans.

	Acknowledgments

 
This work was supported by NIH/NHLBI grants R29HL-47405 and HL-45090 and by grant NO1-HC95162. Dr Garot was supported by the Institut Roche Cardiovasculaire (Paris, France). The authors thank T. Faranesh for the development of the pulse sequence and conventional strain analysis software.

	Appendix A

Top Abstract Introduction Methods Results Discussion Appendix A Appendix B References

 
Harmonic Images
Let yR² represent the coordinates of a point in the imaging plane. The intensity of this point at time t is given by the scalar quantity I(y, t). The tagged image can be written as

where I^o represents the image without tag lines and f is a periodic function describing the tag pattern. We can use the Fourier series to expand the SPAMM-tag pattern into a harmonic summation, as

where K depends on the tag pattern specifications. Each term I_k is a complex image corresponding to a harmonic image. A harmonic image is given by a magnitude D_k and a phase _k as

with

where _k is the tag frequency vector of the corresponding harmonic peak. The phase is related to the motion according to the reference mapping function q(y, t) that maps any point at (y, t) into its reference location q when the tag pattern was imposed.

Wrapping Artifact
We obtained a wrapped version of the phase corresponding to the harmonic phase image a_k. HARP images are related to the actual phase _k

where W is the nonlinear wrapping function given by

	Appendix B

Top Abstract Introduction Methods Results Discussion Appendix A Appendix B References

 
The Eulerian strain is related to the difference in displacement between adjacent parts of the myocardium. It can be computed from 2 HARP images (a₁ and a₂) having 2 linearly independent vectors ₁ and ₂. The apparent strain in the direction e is given by

where

the matrix R^2x2=[₁₂] and

for k=1,2.

Received July 28, 1999; revision received September 20, 1999; accepted October 6, 1999.

	References

Top Abstract Introduction Methods Results Discussion Appendix A Appendix B References

 

1. Marwick TH. Current status of stress echocardiography for diagnosis and prognostic assessment of coronary artery disease. Coron Artery Dis. 1998;9:411–426.[Medline] [Order article via Infotrieve]
   
2. Nagel E, Lehmkuhl HB, Bocksch W, Klein C, Vogel U, Frantz E, Ellmer A, Dreysse E, Fleck E. Noninvasive diagnosis of ischemia-induced wall motion abnormalities with the use of high-dose dobutamine stress magnetic resonance imaging: comparison with dobutamine stress echocardiography. Circulation. 1999;99:763–770.[Abstract/Free Full Text]
   
3. Picano E, Lattanzi F, Orlandini A, Marini C, L’Abbate A. Stress echocardiography and the human factor: importance of being expert. J Am Coll Cardiol. 1991;17:666–669.[Abstract]
   
4. 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.[Abstract/Free Full Text]
   
5. Young AA, Axel L. Three-dimensional motion and deformation of the heart wall: estimation with spatial modulation of magnetization: a model-based approach. Radiology. 1992;185:241–247.[Abstract/Free Full Text]
   
6. O’Dell WG, Moore CC, Hunter WC, Zerhouni EA, McVeigh ER. Three-dimensional myocardial deformations calculation with displacement field fitting to tagged MR images. Radiology. 1995;195:829–835.[Abstract/Free Full Text]
   
7. Axel L, Dougherty L. MR imaging of motion with spatial modulation of magnetization. Radiology. 1989;171:841–845.[Abstract/Free Full Text]
   
8. Axel L, Dougherty L. Heart wall motion: improved method of spatial modulation of magnetization for MR imaging. Radiology. 1989;172:349–350.[Abstract/Free Full Text]
   
9. Denney TS Jr, Prince JL. Reconstruction of 3D left ventricular motion from planar tagged cardiac MR images: an estimation theoretic approach. IEEE Trans Med Imaging. 1995;14:625–635.[Medline] [Order article via Infotrieve]
   
10. Young AA, Kraitchman DL, Dougherty L, Axel L. Tracking and finite element analysis of stripe deformation in magnetic resonance tagging. IEEE Trans Med Imaging. 1995;14:413–421.[Medline] [Order article via Infotrieve]
    
11. Guttman M, Prince JL, McVeigh ER. Tag and contour detection in tagged cardiac magnetic resonance images of the left ventricle. IEEE Trans Med Imaging. 1994;13:74–88.[Medline] [Order article via Infotrieve]
    
12. Moore CC, Reeder SB, McVeigh ER. Tagged MR imaging in a deforming phantom: photographic validation. Radiology. 1994;190:765–769.[Abstract/Free Full Text]
    
13. Guttman MA, Zerhouni EA, McVeigh ER. Analysis and visualization of cardiac function from MR images. IEEE Comp Graph Appl. 1997;17:30–38.
    
14. Osman NF, Prince JL. Angle images for measuring heart motion from tagged MRI. Proc IEEE Conf Image Proc. 1998;I:704–708.
    
15. McVeigh ER, Atalar E. Cardiac tagging with breath-hold cine MRI. Magn Reson Med.. 1992;28:318–327.[Medline] [Order article via Infotrieve]
    
16. Osman NF, Faranesh AZ, McVeigh ER, Prince JL. Tracking cardiac motion using cine harmonic phase (HARP) magnetic resonance imaging. Proc Int Soc Magn Reson Med. 1999;22–28.
    
17. Young AA, Imai H, Chang CN, Axel L. Two-dimensional left ventricular deformation during systole using magnetic resonance imaging with a spatial modulation of magnetization. Circulation. 1994;89:740–752.[Abstract/Free Full Text]
    
18. Clark NR, Reichek N, Bergey P, Hoffman EA, Brownson D, Palmon L, Axel L. Circumferential myocardial shortening in the normal human left ventricle: assessment by magnetic resonance imaging using spatial modulation of magnetization. Circulation. 1991;84:67–74.[Abstract/Free Full Text]
    
19. Power TP, Kramer CM, Shaffer AL, Theobald TM, Petruolo S, Reichek N, Rogers WJ Jr. Breath-hold dobutamine magnetic resonance myocardial tagging: normal left ventricular response. Am J Cardiol. 1997;80:1203–1207.[Medline] [Order article via Infotrieve]
    
20. Scott CH, Sutton MS, Gusani N, Fayad Z, Kraitchman D, Keane MG, Axel L, Ferrari VA. Effect of dobutamine on regional left ventricular function measured by tagged magnetic resonance imaging in normal subjects. Am J Cardiol. 1999;83:412–417.[Medline] [Order article via Infotrieve]
    
21. Buchalter MB, Weiss JL, Rogers WJ, Zerhouni EA, Weisfeldt ML, Beyar R, Shapiro EP. Noninvasive quantification of left ventricular rotational deformation in normal humans using magnetic resonance imaging myocardial tagging. Circulation. 1990;81:1236–1244.[Abstract/Free Full Text]
    
22. Geskin G, Kramer CM, Rogers WJ, Theobald TM, Pakstis D, Hu YL, Reichek N. Quantitative assessment of myocardial viability after infarction by dobutamine magnetic resonance tagging. Circulation. 1998;98:217–223.[Abstract/Free Full Text]
    
23. Baer FM, Theissen P, Schneider CA, Voth E, Sechtem U, Schicha H, Erdmann E. Dobutamine magnetic resonance imaging predicts contractile recovery of chronically dysfunctional myocardium after successful revascularization. J Am Coll Cardiol. 1998;31:1040–1048.[Abstract/Free Full Text]
    
24. Kramer CM, Lima JAC, Reichek N, Ferrari VA, Llaneras MR, Palmon LC, Yeh IT, Tallant B, Axel L. Regional differences in function within noninfarcted myocardium during left ventricular remodeling. Circulation. 1993;88:1279–1288.[Abstract/Free Full Text]
    
25. Lima JAC, Ferrari VA, Reichek N, Kramer CM, Palmon L, Llaneras MR, Tallant B, Young AA, Axel L. Segmental motion and deformation of transmurally infarcted myocardium in acute postinfarct period. Am J Physiol. 1995;268:H1304–H1312.[Abstract/Free Full Text]
    
26. Hansen DE, Daughters GT II, Alderman EL, Ingels NB Jr, Miller DC. Torsional deformation of the left ventricular midwall in human hearts with intramyocardial markers: regional heterogeneity and sensitivity to the inotropic effects of abrupt rate changes. Circ Res. 1988;62:941–952.[Abstract/Free Full Text]
    
27. Rogers WJ Jr, Shapiro EP, Weiss JL, Buchalter MB, Rademakers FE, Weisfeldt ML, Zerhouni EA. Quantification of and correction for left ventricular systolic long-axis shortening by magnetic resonance tissue tagging and slice isolation. Circulation. 1991;84:721–731.[Abstract/Free Full Text]
    
28. Stuber M, Scheidegger MB, Fischer SE, Nagel E, Steinemann F, Hess OM, Boesiger P. Alterations in the local myocardial motion pattern in patients suffering from pressure overload due to aortic stenosis. Circulation. 1999;100:361–368.[Abstract/Free Full Text]
    
29. Lima JAC, Jeremy R, Guier W, Bouton S, Zerhouni EA, McVeigh ER, Buchalter MB, Weisfeldt ML, Shapiro EP, Weiss JL. Accurate systolic wall thickening by nuclear magnetic resonance imaging with tissue tagging: correlation with sonomicrometers in normal and ischemic myocardium. J Am Coll Cardiol. 1993;21:1741–1751.[Abstract]
    

This article has been cited by other articles:

				T. A. Dorfman, B. D. Rosen, M. A. Perhonen, T. Tillery, R. McColl, R. M. Peshock, and B. D. Levine Diastolic suction is impaired by bed rest: MRI tagging studies of diastolic untwisting J Appl Physiol, April 1, 2008; 104(4): 1037 - 1044. [Abstract] [Full Text] [PDF]

				A. M. Sironi, A. Pingitore, S. Ghione, D. De Marchi, B. Scattini, V. Positano, E. Muscelli, D. Ciociaro, M. Lombardi, E. Ferrannini, et al. Early Hypertension Is Associated With Reduced Regional Cardiac Function, Insulin Resistance, Epicardial, and Visceral Fat Hypertension, February 1, 2008; 51(2): 282 - 288. [Abstract] [Full Text] [PDF]

				M. Tumkosit, C. G. Martin, E. Bayram, T. M. Morgan, K. S. Lane, P. Rerkpattanapipat, C. A. Hamilton, K. M. Link, and W. G. Hundley Left Ventricular Spherical Remodeling and Apical Myocardial Relaxation: Cardiovascular MR Imaging Measurement of Myocardial Segments Radiology, August 1, 2007; 244(2): 411 - 418. [Abstract] [Full Text] [PDF]

				E. Pettersen, T. Helle-Valle, T. Edvardsen, H. Lindberg, H.-J. Smith, B. Smevik, O. A. Smiseth, and K. Andersen Contraction Pattern of the Systemic Right Ventricle: Shift From Longitudinal to Circumferential Shortening and Absent Global Ventricular Torsion J. Am. Coll. Cardiol., June 26, 2007; 49(25): 2450 - 2456. [Abstract] [Full Text] [PDF]

				K. Nasir, M. Tsai, B. D. Rosen, V. Fernandes, D. A. Bluemke, A. R. Folsom, and J. A.C. Lima Elevated Homocysteine Is Associated With Reduced Regional Left Ventricular Function: The Multi-Ethnic Study of Atherosclerosis Circulation, January 16, 2007; 115(2): 180 - 187. [Abstract] [Full Text] [PDF]

				L. C. Amado, K. H. Schuleri, A. P. Saliaris, A. J. Boyle, R. Helm, B. Oskouei, M. Centola, V. Eneboe, R. Young, J. A.C. Lima, et al. Multimodality Noninvasive Imaging Demonstrates In Vivo Cardiac Regeneration After Mesenchymal Stem Cell Therapy J. Am. Coll. Cardiol., November 21, 2006; 48(10): 2116 - 2124. [Abstract] [Full Text] [PDF]

				W. Liu, M. W. Ashford, J. Chen, M. P. Watkins, T. A. Williams, S. A. Wickline, and X. Yu MR tagging demonstrates quantitative differences in regional ventricular wall motion in mice, rats, and men Am J Physiol Heart Circ Physiol, November 1, 2006; 291(5): H2515 - H2521. [Abstract] [Full Text] [PDF]

				B. D. Rosen, J. A.C. Lima, K. Nasir, T. Edvardsen, A. R. Folsom, S. Lai, D. A. Bluemke, and M. Jerosch-Herold Lower Myocardial Perfusion Reserve Is Associated With Decreased Regional Left Ventricular Function in Asymptomatic Participants of the Multi-Ethnic Study of Atherosclerosis Circulation, July 25, 2006; 114(4): 289 - 297. [Abstract] [Full Text] [PDF]

				V. R.S. Fernandes, J. F. Polak, T. Edvardsen, B. Carvalho, A. Gomes, D. A. Bluemke, K. Nasir, D. H. O'Leary, and J. A.C. Lima Subclinical Atherosclerosis and Incipient Regional Myocardial Dysfunction in Asymptomatic Individuals: The Multi-Ethnic Study of Atherosclerosis (MESA) J. Am. Coll. Cardiol., June 20, 2006; 47(12): 2420 - 2428. [Abstract] [Full Text] [PDF]

				I. Paetsch, C. Jahnke, V. A. Ferrari, F. E. Rademakers, P. A. Pellikka, W. G. Hundley, D. Poldermans, J. J. Bax, K. Wegscheider, E. Fleck, et al. Determination of interobserver variability for identifying inducible left ventricular wall motion abnormalities during dobutamine stress magnetic resonance imaging Eur. Heart J., June 2, 2006; 27(12): 1459 - 1464. [Abstract] [Full Text] [PDF]

				B. D. Rosen, M. F. Saad, S. Shea, K. Nasir, T. Edvardsen, G. Burke, M. Jerosch-Herold, D. K. Arnett, S. Lai, D. A. Bluemke, et al. Hypertension and Smoking Are Associated With Reduced Regional Left Ventricular Function in Asymptomatic Individuals: The Multi-Ethnic Study of Atherosclerosis J. Am. Coll. Cardiol., March 21, 2006; 47(6): 1150 - 1158. [Abstract] [Full Text] [PDF]

				S. E. Petersen, B. A. Jung, F. Wiesmann, J. B. Selvanayagam, J. M. Francis, J. Hennig, S. Neubauer, and M. D. Robson Myocardial Tissue Phase Mapping with Cine Phase-Contrast MR Imaging: Regional Wall Motion Analysis in Healthy Volunteers Radiology, March 1, 2006; 238(3): 816 - 826. [Abstract] [Full Text] [PDF]

				T. Edvardsen, R. Detrano, B. D. Rosen, J. J. Carr, K. Liu, S. Lai, S. Shea, L. Pan, D. A. Bluemke, and J. A.C. Lima Coronary Artery Atherosclerosis Is Related to Reduced Regional Left Ventricular Function in Individuals Without History of Clinical Cardiovascular Disease: The Multiethnic Study of Atherosclerosis Arterioscler. Thromb. Vasc. Biol., January 1, 2006; 26(1): 206 - 211. [Abstract] [Full Text] [PDF]

				T. Helle-Valle, J. Crosby, T. Edvardsen, E. Lyseggen, B. H. Amundsen, H.-J. Smith, B. D. Rosen, J. A.C. Lima, H. Torp, H. Ihlen, et al. New Noninvasive Method for Assessment of Left Ventricular Rotation: Speckle Tracking Echocardiography Circulation, November 15, 2005; 112(20): 3149 - 3156. [Abstract] [Full Text] [PDF]

				I. Paetsch, C. Jahnke, E. Fleck, and E. Nagel Current clinical applications of stress wall motion analysis with cardiac magnetic resonance imaging Eur J Echocardiogr, October 1, 2005; 6(5): 317 - 326. [Abstract] [Full Text] [PDF]

B. D. Rosen, T. Edvardsen, S. Lai, E. Castillo, L. Pan, M. Jerosch-Herold, S. Sinha, R. Kronmal, D. Arnett, J. R. Crouse III, et al. Left Ventricular Concentric Remodeling Is Associated With Decreased Global and Regional Systolic Function: The Multi-Ethnic Study of Atherosclerosis Circulation, August 16, 2005; 112(7): 984 - 991. [Abstract]