Nael F. Osman

Cardiac Motion Tracking Using CINE Harmonic Phase (HARP) Magnetic Resonance Imaging

This article introduces a new image processing technique for rapid analysis of tagged cardiac magnetic resonance image sequences. The method uses isolated spectral peaks in SPAMM‐tagged magnetic resonance images, which contain information about cardiac motion. The inverse Fourier transform of a spectral peak is a complex image whose calculated angle is called a harmonic phase (HARP) image. It is shown how two HARP image sequences can be used to automatically and accurately track material points through time. A rapid, semiautomated procedure to calculate circumferential and radial Lagrangian strain from tracked points is described. This new computational approach permits rapid analysis and visualization of myocardial strain within 5–10 min after the scan is complete. Its performance is demonstrated on MR image sequences reflecting both normal and abnormal cardiac motion. Results from the new method are shown to compare very well with a previously validated tracking algorithm. Magn Reson Med 42:1048–1060, 1999.

Imaging heart motion using harmonic phase MRI

Describes a new image processing technique for rapid analysis and visualization of tagged cardiac magnetic resonance (MR) images. The method is based on the use of isolated spectral peaks in spatial modulation of magnetization (SPAMM)-tagged magnetic resonance images. The authors call the calculated angle of the complex image corresponding to one of these peaks a harmonic phase (HARP) image and show that HARP images can be used to synthesize conventional tag lines, reconstruct displacement fields for small motions, and calculate two-dimensional (2-D) strain. The performance of this new approach is demonstrated using both real and simulated tagged MR images. Potential for use of HARP images in fast imaging techniques and three-dimensional (3-D) analyses are discussed.

Fast determination of regional myocardial strain fields from tagged cardiac images using harmonic phase MRI

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 microg(-1). kg(-1). min(-1) dobutamine and in 9 postinfarct patients at rest. We compared 2D myocardial strains (circumferential shortening, Ecc; maximal shortening, E(2); and E(2), 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(2) 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(2) 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(2) values were dramatically reduced in the asynergic segments as compared with remote (P<0.001), and E(2) 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.

Imaging longitudinal cardiac strain on short-axis images using strain-encoded MRI.

This article presents a new method for measuring longitudinal strain in a short‐axis section of the heart using harmonic phase magnetic resonance imaging (HARP‐MRI). The heart is tagged using 1‐1 SPAMM at end‐diastole with tag surfaces parallel to a short‐axis imaging plane. Two or more images are acquired such that the images have different phase encodings in a direction orthogonal to the image plane. A dense map of the longitudinal strain can be computed from these images using a simple, fast computation. Simulations are conducted to study the effect of noise and the choice of out‐of‐plane phase encoding values. Longitudinal strains acquired from a normal human male are shown. Magn Reson Med 46:324‐334, 20001.

Three-Dimensional Mapping of Optimal Left Ventricular Pacing Site for Cardiac Resynchronization

Background— The efficacy of cardiac resynchronization therapy (CRT) depends on placement of the left ventricular lead within the late-activated territory. The geographic extent and 3-dimensional distribution of left ventricular (LV) locations yielding optimal CRT remain unknown. Methods and Results— Normal or tachypacing-induced failing canine hearts made dyssynchronous by right ventricular free wall pacing or chronic left bundle-branch ablation were acutely instrumented with a nonconstraining epicardial elastic sock containing 128 electrodes interfaced with a computer-controlled stimulation/recording system. Biventricular CRT was performed using a fixed right ventricular site and randomly selected LV sites covering the entire free wall. For each LV site, global cardiac function (conductance catheter) and mechanical synchrony (magnetic resonance imaging tagging) were determined to yield 3-dimensional maps reflecting CRT impact. Optimal CRT was achieved from LV lateral wall sites, slightly more anterior than posterior and more apical than basal. LV sites yielding ≥70% of the maximal dP/dtmax increase covered ≈43% of the LV free wall. This distribution and size were similar in both normal and failing hearts. The region was similar for various systolic and diastolic parameters and correlated with 3-dimensional maps based on mechanical synchrony from magnetic resonance imaging strain analysis. Conclusions— In hearts with delayed lateral contraction, optimized CRT is achieved over a fairly broad area of LV lateral wall in both nonfailing and failing hearts, with modest anterior or posterior deviation still capable of providing effective CRT. Sites selected to achieve the most mechanical synchrony are generally similar to those that most improve global function, confirming a key assumption underlying the use of wall motion analysis to optimize CRT.

Visualizing myocardial function using HARP MRI

Harmonic phase magnetic resonance imaging (HARP) is a new technique for measuring the motion of the left ventricle of the heart. HARP uses magnetic resonance tagging, Fourier filtering and special processing algorithms to calculate key indices of myocardial motion including Eulerian and Lagrangian strain. This paper presents several new methods for visualizing myocardial motion based on HARP. Quantities that are computed and visualized include motion grids, velocity fields, strain rates, pathlines, tracked Eulerian strain, and contraction angle. The computations are fast and fully automated and have the potential for clinical application.

Myocardial tissue tagging with cardiovascular magnetic resonance

Cardiovascular magnetic resonance (CMR) is currently the gold standard for assessing both global and regional myocardial function. New tools for quantifying regional function have been recently developed to characterize early myocardial dysfunction in order to improve the identification and management of individuals at risk for heart failure. Of particular interest is CMR myocardial tagging, a non-invasive technique for assessing regional function that provides a detailed and comprehensive examination of intra-myocardial motion and deformation. Given the current advances in gradient technology, image reconstruction techniques, and data analysis algorithms, CMR myocardial tagging has become the reference modality for evaluating multidimensional strain evolution in the human heart. This review presents an in depth discussion on the current clinical applications of CMR myocardial tagging and the increasingly important role of this technique for assessing subclinical myocardial dysfunction in the setting of a wide variety of myocardial disease processes.

Quantitative assessment of regional myocardial function with MR-tagging in a multi-center study: interobserver and intraobserver agreement of fast strain analysis with Harmonic Phase (HARP) MRI

PURPOSE To assess the reproducibility of Harmonic Phase (HARP) analysis of myocardial MR tagged images acquired in the Multi-Center Study of Atherosclerosis (MESA). METHODS Using the HARP method, three independent observers performed two separate quantitative strain analyses of myocardial cine MR-tagging images blindly in 24 participants. The images were obtained in four different centers and analyzed at a single core lab. Each study comprised 3 short-axis slices subdivided in 12 segments (24 x 3 x 12 = 864 segments), each with three layers. Normal strains (circumferential [Ecc] and radial [Err]), principal strains (Lambda1, Lambda2), and the angle alpha (between Ecc-Lambda2) were calculated. Intraclass correlation coefficient (R) for peak systolic strains, and all pooled systolic and diastolic strain data were used to determine inter- and intraobserver agreement. Two observers also visually graded study quality. R values were related to the image quality in different myocardial regions and layers. RESULTS Overall, HARP yielded an excellent inter- and intraobserver agreement for peak systolic strain data (for Ecc, R = 0.84 and 0.89, respectively) and all systolic pooled data (for Ecc, interobserver R = 0.82, intraobserver R = 0.69-0.76). Both inter and intraobserver agreement were lower for diastolic pooled data (R = 0.69 and 0.58-0.62, respectively). There was a direct relationship between image quality and performance of the HARP analysis, with increasing inter- and intraobserver R values in studies with longer tag persistence. Both inter- and intraobserver agreement were better in the anterior and septal myocardial regions, and in the midwall layer. The intraobserver agreement was similar among the three observers. CONCLUSION Employing the HARP method for quantitative strain analysis of myocardial MR tagged images provides a high inter- and intraobserver agreement. These good results are obtained in case of good to excellent MR image quality.

Fast tracking of cardiac motion using 3D-HARP

Magnetic resonance (MR) tagging is capable of accurate, noninvasive quantification of regional myocardial function. Routine clinical use, however, is hindered by cumbersome and time-consuming postprocessing procedures. We propose a fast, semiautomatic method for tracking three-dimensional (3-D) cardiac motion from a temporal sequence of short- and long-axis tagged MR images. The new method, called 3-D-HARmonic Phase (3D-HARP), extends the HARP approach, previously described for two-dimensional (2-D) tag analysis, to 3-D. A 3-D material mesh model is built to represent a collection of material points inside the left ventricle (LV) wall at a reference time. Harmonic phase, a material property that is time-invariant, is used to track the motion of the mesh through a cardiac cycle. Various motion-related functional properties of the myocardium, such as circumferential strain and left ventricular twist, are computed from the tracked mesh. The correlation analysis of 3D-HARP and FINDTAGS + Tag Strain(E) Analysis (TEA), which are well-established tag analysis techniques, shows that the regression coefficients of circumferential strain (E/sub CC/) and twist angle are r/sup 2/=0.8605 and r/sup 2/=0.8645, respectively. The total time required for tracking 3-D cardiac motion is approximately 10 min in a 9 timeframe tagged MRI dataset and has the potential to be much faster.

Quantitative Ischemia Detection During Cardiac Magnetic Resonance Stress Testing by Use of FastHARP

Background—Because ECG alterations caused by ischemia cannot be reliably detected in the high-field MRI environment, detection of wall motion abnormalities is often used to ensure patient safety during stress testing. However, an experienced observer is needed to detect these abnormalities. In this study, we investigate the use of fast harmonic phase (FastHARP) MRI for the quantitative, operator-independent detection of the onset of ischemia during acute coronary occlusion. Methods and Results—Eight mongrel dogs underwent an acute 2-minute closed-chest coronary artery occlusion while continuous FastHARP images were acquired. Full regional wall strain was determined every other heartbeat in a single short-axis imaging slice. After 5 minutes of reperfusion, a second 2-minute ischemic episode was induced during the acquisition of conventional cine wall-motion images. The time at which ECG alterations were observed during the first ischemic period was recorded. The time from occlusion to the detection of ischemia, based on a consensus of 2 blinded observers, was determined for MRI. No significant ischemia was present in 2 animals. In the remaining animals, the onset of ischemia was detected significantly earlier by FastHARP than by cine MRI (9.5±5 versus 33±14 seconds, P <0.01). HARP ischemia detection preceded ECG changes, on average, by 54 seconds. Conclusions—The rapid acquisition and detection of induced ischemia with FastHARP MRI shows promise as a nonsubjective method to diagnose significant coronary lesions during MR stress testing.