Regional heart wall motion: two-dimensional analysis and functional imaging with MR imaging

Two-dimensional spatial and temporal displacement and deformation field fitting from cardiac magnetic resonance tagging

Tagged magnetic resonance imaging is a specially developed technique to noninvasively assess contractile function of the heart. Several methods have been developed to estimate myocardial deformation from tagged image data. Most of these methods do not explicitly impose a continuity constraint through time although myocardial motion is a continuous physical phenomenon. In this paper, we propose to model the spatio-temporal myocardial displacement field by a cosine series model fitted to the entire tagged dataset. The method has been implemented in two dimensions (2D)+time. Its accuracy was successively evaluated on actual tagged data and on a simulated two-dimensional (2D) moving heart model. The simulations show that an overall theoretical mean accuracy of 0.1 mm can be attained with adequate model orders. The influence of the tagging pattern was evaluated and computing time is provided as a function of the model complexity and data size. This method provides an analytical and hierarchical model of the 2D+time deformation inside the myocardium. It was applied to actual tagged data from a healthy subject and from a patient with ischemia. The results demonstrate the adequacy of the proposed model for this evaluation.

Three-dimensional myocardial strain analysis based on short- and long-axis magnetic resonance tagged images using a 1D displacement field

A robust algorithm to estimate three-dimensional strain in the left-ventricular heart wall, based on magnetic resonance (MR) grid-tagging in two sets of orthogonal image planes, is presented. Starting-point of this study was to minimize global interpolation and smoothing. Only the longitudinal displacement was interpolated between long-axis images. Homogeneous strain analysis was performed using small tetrahedrons. The method was tested using a stack of short-axis images and three long-axis images in six healthy volunteers. In addition, the method was subjected to an analytical test case, in which the effect of noise in tag point position on the observed strains was explored for normally distributed noise (0.5 mm RMS). In volunteers, the error in the longitudinal displacement due to interpolation between the long-axis image planes was -0.10 +/- 0. 48 mm (mean +/- SD). The resulting error in the longitudinal strain epsilon(l) was -0.003 +/- 0.02. The analytical test case was used to quantify the effects of three sources of errors on the observed strain. The SD of the difference between homogeneous strain and true strain was 0.06 for epsilon(r.) The error due to the 3-D reconstruction was 0.004 for epsilon(r.) The error in epsilon(r) resulting from simulated noise in the tag point position was 0.10. Equivalent results were obtained for all other strain parameters; thus, the error resulting from noise in the tag point position dominates the error introduced by approximations in the method. Because the proposed method uses a minimum of global interpolation and smoothing, it offers the prospect to detect small regions of aberrant contraction.

Cardiac phase contrast gradient echo MRI: measurement of myocardial wall motion in healthy volunteers and patients

A number of methods have been proposed for the noninvasive measurement of myocardial wall motion. The paper describes a strategy for assessing myocardial motion based on the sensitivity of the phase of the MR-signal to motion using a breath-hold phase contrast technique. A motion-sensitized and a motion-compensated MR-signal are measured during successive scans. The difference between the two MR-signals is used to calculate myocardial velocity in all three spatial dimensions. Postprocessing includes the transformation of the measured velocities into an internal coordinate system of the left ventricle. Also various presentation modes and further processing of the received velocity information are provided including calculation of global motion parameters. We examined 20 patients suffering from myocardial infarction. The overall left ventricular motion can be characterized by appropriate parameters describing the rotation and contraction or expansion, respectively. Regional motional disturbances are visualized using parametric images. Contrary to the highly consistent interindividual data in normal volunteers, patients showed significant localized motion deficits.

Recognition of infarct localization by specific changes in intramural myocardial mechanics

BACKGROUND

After transmural myocardial infarction (MI), changes occur in intramural myocardial function. This has been described in anterior MI only. The aim of this study was to determine the relation between variable infarct locations and intramural deformation in patients with a first MI.

METHODS

Forty patients (33 men and 7 women aged 57 +/- 11 years) with different infarct-related coronary arteries (25 left anterior descending, 7 circumflex, and 8 right coronary) were studied 6 +/- 3 days after infarction with magnetic resonance tissue tagging and 2-dimensional finite element analysis of myocardial deformation. Short-axis tagged images were acquired at base, mid, and apical level. Intramural deformation was measured in 6 circumferential segments per level. Results were compared with 9 age-matched healthy controls.

RESULTS

Each infarct area demonstrated a significant reduction of intramural deformation. At mid-ventricular level, segments with maximum impaired intramural function were the anteroseptal segment for left anterior descending-related MI (stretch: 16% vs 33% for controls, P <.001), the posterolateral segment for related MI (stretch: 20% vs 34%, P <. 01); and the inferior segment for right coronary artery related MI (stretch: 18% vs 25%, P =.082). In these infarct segments, the intramural regional systolic stretch was more circumferentially oriented compared with radially oriented stretch in the same segments in controls (P <.05).

CONCLUSION

The infarct area can be recognized by a specific spatial pattern of intramural deformation. In infarcted compared with noninfarcted myocardium, deformation is significantly reduced and systolic stretch deviates from the radial direction. Left anterior descending related infarcts were found to have larger regional differences in intramural deformation than circumflex or right coronary artery related MI of enzymatically the same size.

Cascaded MRI-SPAMM for LV motion analysis during a whole cardiac cycle

We present a new paradigm which incorporates multiple sets of tagged MRI data (MRI-SPAMM) acquired in a cascaded fashion in order to estimate the full 3-D motion of the left ventricle (LV) during its entire cardiac cycle. Our technique is based on an extension of our volumetric physics-based deformable models, whose parameters are functions. Using these parameters, we can characterize the local shape variation of an object with a small number of intuitive parameters. By integrating a cascaded sequence of SPAMM data sets into our modeling technique, we have extended the capability of the MRI-SPAMM technique and have provided an accurate representation of the LV motion during the full cardiac cycle (from end-diastole to end-diastole) to better understand cardiac mechanics.

A spatiotemporal model of cyclic kinematics and its application to analyzing nonrigid motion with MR velocity images

We present a method (DMESH) for nonrigid cyclic motion analysis using a series of velocity images covering the cycle acquired, for example, from phase-contrast magnetic resonance imaging. The method is based on fitting a dynamic finite-element mesh model to velocity samples of an extended region, at all time frames. The model offers a flexible tradeoff between accuracy and reproducibility with controllable built-in spatiotemporal smoothing, which is determined by the fineness of the initially defined mesh and the richness of included Fourier harmonics. The method can further provide a prediction of the analysis reproducibility, along with the estimated motion and deformation quantities. Experiments have been conducted to validate the method and to verify the reproducibility prediction. Use of the method for motion analysis using displacement information (e.g., from magnetic resonance tagging) has also been explored.

Kinematic analysis of left ventricular deformation in myocardial infarction using magnetic resonance cardiac tagging

The Magnetic Resonance (MR) tagging technique provides detailed information about 2D motion in the plane of observation. Interpretation of this information as a reflection of the 3D motion of the entire cardiac wall is a major problem. In finite element models of the mechanics of the infarcted heart, an infarcted region causes motional asymmetry, extending far beyond the infarct boundary. Here we present a method to quantify such asymmetry in amplitude and orientation. For this purpose images of a short-axis cross-section of the ejecting left ventricle were acquired from 9 healthy volunteers and 5 patients with myocardial infarction. MR-tags were applied in a 5 mm grid at end-diastole. The tags were tracked by video-image analysis. Tag motion was fitted to a kinematic model of cardiac motion. For the volunteers and the patients the center of the cavity displaced by about the same amount (p = 0.11) during the ejection phase: 3.8 +/- 1.4 and 3.0 +/- 0.9 mm (mean +/- sd), respectively. Cross-sectional rotation and the decrease in cross-sectional area of the cavity were both greater in the volunteers than in the patients: 6.4 +/- 1.5 vs. 3.0 +/- 0.8 degrees (p < 0.001), and 945 +/- 71 vs. 700 +/- 176 mm2 (p = 0.02), respectively. In the patients, asymmetry of wall motion, as expressed by a sine wave dependency of contraction around the circumference, was significantly enlarged (p = 0.02). The proposed method of kinematic analysis can be used to assess cardiac deformation in humans. We expect that by analyzing images of more cross-sections simultaneously, the 3D location and the degree of infarction can be assessed efficiently.

Mapping propagation of mechanical activation in the paced heart with MRI tagging

The temporal evolution of three-dimensional (3-D) strain maps derived from magnetic resonance imaging (MRI) tagging were used to noninvasively evaluate mechanical activation in the left ventricle (LV) while seven canine hearts were paced in situ from three different sites: the base of the LV free wall (LVb), the right ventricular apex (RVa), and the right atrium (RA). Strain maps plotted against time showed the evolution of shortening over the entire LV midwall and were used to generate mechanical activation maps showing the onset of circumferential shortening. RA pacing showed rapid synchronous shortening; LVb pacing showed a wave front of mechanical activation propagating slowly and steadily from the pacing site, whereas RVa pacing showed regions of rapid and slower propagation. The mechanical (M) activation times correlated linearly with the electrical (E) activation (M = 1.06E + 8.4 ms, R = 0.95). The time for 90% activation of the LV was 63.1 +/- 24.3 ms for RA pacing, 130.2 +/- 9.8 ms for LVb pacing, and 121.3 +/- 17.9 ms for RVa pacing. The velocity of mechanical activation was calculated for LVb and RVa pacing and was similar to values reported for electrical conduction in myocardium. The propagation of mechanical activation for RVa pacing showed regional variations, whereas LVb pacing did not.

Simultaneous MRI tagging and through-plane velocity quantification: a three-dimensional myocardial motion tracking algorithm

A tracking algorithm was developed for calculation of three-dimensional point-specific myocardial motion. The algorithm was designed for images acquired with simultaneous magnetic resonance imaging (MRI) grid tagging and through-plane velocity quantification. The tagging grid provided the in-plane motion while the velocity quantification measured the through-plane motion. In four healthy volunteers, the in vivo performance was evaluated by comparing the systolic through-plane displacement with the displacement of tagging-grid intersections in long-axis images. The correlation coefficient was 0.93 (P < 0.001, N = 183). A t-test for paired samples revealed a small underestimation of the through-plane displacement by 0.04 +/- 0.09 cm (mean +/- SD, P < 0.001) on an average displacement of 0.77 +/- 0.23 cm toward the apex. The authors conclude that three-dimensional point-specific motion tracking based on simultaneous tagging and velocity quantification is competitive with other methods such as tagging in mutually orthogonal image planes or quantification of three orthogonal velocity components.

Toward high-resolution myocardial tagging

Magnetic resonance imaging with preceding tissue tagging is a robust method for assessing cardiac motion of the entire heartbeat cycle with a high degree of accuracy. One limitation of this technique, however, is the low resolution of the obtained displacement map of the labeled points within the myocardium. By a new tagging technique, which is based on the combination of two or more measurements of the same slice but with different grid positions, a highly improved resolution of cardiac motion data can be achieved. In combination with a multi-heart-phase echo-planar imaging sequence, such images with doubled grid frequency can be acquired in two short breath-hold periods.

Frequency-domain simulation of MR tagging

Simulation of MR images is a useful tool for offline sequence development and as an aid to understanding image formation. One particular application of simulation is MR tagging, which is used for tracking myocardial motion. Simple spatial-domain methods cannot adequately represent effects common in these images, such as motion artifact and signal wrap. An existing frequency-domain model is shown to be inappropriate for tagged images, and an extension based on the Bloch equations and Fourier shift theorem is described to correct this. Software incorporating the new model is used to generate ideal tag intensity profiles and to accurately simulate tagged images. The shifted k-space patterns associated with tagged images, and their dependence on the order of the binomial tagging sequence, are explained. An application of the Fourier shift theorem is suggested that allows more rapid simulation of static tagged images.

Analysis of myocardial motion based on velocity measurements with a black blood prepared segmented gradient-echo sequence: methodology and applications to normal volunteers and patients

The paper describes a strategy for measuring and characterizing myocardial motion in terms of velocity parameters derived from measurements with a segmented black blood prepared phase contrast gradient echo sequence. The characteristic parameters are calculated by transforming the velocities measured on a pixel-by-pixel basis across the left ventricle from the laboratory frame of reference into a cylindrical coordinate system, in which the motion velocities within the short axis plane are represented in polar coordinates and which is located at the center of the myocardium and moving with it over the ECG cycle. First results in a study with 12 healthy volunteers gave highly consistent values for the radial (expansion/compression) as well as the rotational velocities. Except for one volunteer, motion at the R wave of the ECG starts with clockwise rotation, followed by contraction and expansion accompanied by counterclockwise rotation. First examinations of patients with global and focal disease demonstrate the potential to detect disturbances in the local as well as the overall motion patterns.

Estimation of deformation gradient and strain from cine-PC velocity data

Phase contrast magnetic resonance imaging (MRI) can provide in vivo myocardial velocity field measurements. These data allow densely spaced material points to be tracked throughout the whole heart cycle using, for example, the Fourier tracking algorithm. To process the tracking results for myocardial deformation and strain quantification, we developed a method that is based on fitting the tracking results to an appropriate local deformation model. We further analyzed the accuracy and precision of the method and provided performance predictions for several local models. In order to validate the method and the theoretical performance analysis, we conducted controlled computer simulations and a phantom study. The results agreed well with expectations. Human heart data were also acquired and analyzed, and provided encouraging results. At the signal-to-noise ratio (SNR) level and spatial resolution expected in clinical settings, the study predicts strain quantification accuracy and precision that may allow the technique to become a practical and powerful noninvasive approach for the study of cardiac function, although clinically acceptable data acquisition strategies for three-dimensional (3-D) data are still a challenge.

Myocardial function in infarcted and remote regions early after infarction in man: assessment by magnetic resonance tagging and strain analysis

Early after infarction in the perfusion bed of the left anterior descending coronary artery, cine MRI with spatial modulation of magnetization (SPAMM) tagging (7-mm grid) was used for short- and long-axis cardiac imaging. Two-dimensional strain analysis of triangular finite elements was performed between end-diastole and end-systole. Patients (n = 10) were compared with age-matched healthy subjects (n = 8). The anteroseptal region at midventricular level was considered representative for "infarcted" and the posterolateral region at basal level was considered "remote". The left ventricular end-diastolic volume index was larger in the patients (69 +/- 15 ml/m2 versus 56 +/- 4 ml/m2, P < 0.05). Short-axis images showed in the infarcted region a decrease of first principal strain (greatest systolic lengthening: 1.10 +/- .06 versus 1.27 +/- 0.04, P < 0.0001), and in the remote region an increase (1.48 +/- 0.11 versus 1.36 +/- 0.07, P < 0.025). The lateral and inferior ventricular regions at mid- and basal levels were found to function normally. Long-axis images yielded similar results. Early after infarction, regions with dysfunction, normal function, and hyperfunction can be delineated with MR tagging. The compensatory increased contraction in the remote region is possibly triggered by the Frank-Starling mechanism.

Noninvasive measurement of cardiac strain with MRI

The motion sensitivity of cardiac magnetic resonance imaging (MRI) can be exploited to measure the motion patterns within the heart wall and thus to noninvasively calculate the intramyocardial strain. The resulting large data sets pose a challenge for visualization, but offer the potential of a greatly improved picture of cardiac dynamics. This may have both basic research and clinical applications.

Left ventricular myocardial tagging

Nuclear Magnetic Resonance myocardial tagging is a potent non-invasive technique which enables the quantification of myocardial deformation, globally but also regionally at different time points during the cardiac cycle. By the use of presaturating pulses prior to the actual imaging sequence non-invasive markers or tags can be placed on the myocardium at end diastole, which move and deform with the underlying myocardium on which they were inscribed. Through combination of perpendicular sets of short- and long-axis images with tags, a three-dimensionally reconstructed left ventricle is obtained, subdivided in 32 myocardial cuboids for which the 3D coordinates of the corners are known at different time points in the cardiac cycle. From these data global and regional strains and quantitative measures of shape can be computed.

Simulation of two-dimensional tagged MRI

MR tagging is a recent imaging development that, in cardiac applications, makes possible the tracking of points in the myocardium during the cardiac cycle. Researchers have developed semiautomated, computer-based methods for analyzing tagged images, but the images are complex and present a challenge to automated tracking systems. Simulation can provide an inexhaustible supply of images for testing and validation of tag tracking software and preview the effect of parameter changes in acquisition. SIMTAG is an interactive computer program that simulates two-dimensional tagged-MR experiments. The mathematic model used in the simulation and algorithms for simulating image noise and object deformation are described. Examples of the use of simulated images in SPAMM parameter selection, a comparison of tag contrast in signal-averaged SPAMM and CSPAMM, and simulated images as test sets for tag-tracking software are presented.

Identifying and tracking a guide wire in the coronary arteries during angioplasty from X-ray images

During angioplasty, a guide wire (GW) is routinely placed in the coronary artery. Balloon inflation during angioplasty causes transient occlusion of the coronary artery and regional dysfunction. Thus, it is of major importance to monitor myocardial function, which may be impaired during this period. Since the GW moves with the coronary arteries, information regarding myocardial function can potentially be extracted from the GW motion. An algorithm is suggested which is a step toward such monitoring. The algorithm presented is a semiautomatic procedure for identifying and tracking the GW using specific characteristics of the GW. This algorithm is based on working in limited active windows. A preprocessing stage which enhances the GW by the use of a modified Laplacian filter or a modified Marr-Hildreth filter is introduced. The second stage of the algorithm is the tracking of the GW, which is based on fitting a second-degree polynomial to the GW using the Hough transform in each window. To further improve the results further modifications of the basic algorithms that were taken. A single set of parameters, which enabled good tracking for a large number of images taken during angioplasty, was fitted to the final algorithm.

Polyvinyl alcohol cryogel: an ideal phantom material for MR studies of arterial flow and elasticity

The authors present a unique application of polyvinyl alcohol (PVA) cryogel as an anthropomorphic, elastic, vascular phantom material that can be used in MR imaging. The composition consists of two nontoxic ingredients: water and PVA. The biomechanical and MR properties can be adjusted to be similar to those of excised porcine aortas by varying the number of freeze-thaw cycles to which the PVA solution is exposed. The authors present the T1, T2, shrinkage, and tensile properties of PVA cryogel tubes as a function of freeze-thaw cycles. MR images of a dual elastic aortic phantom undergoing pulsatile motion are shown.

Dynamic magnetic resonance imaging assessment of the effect of ventricular wall curvature on regional function in hypertrophic cardiomyopathy

We hypothesized that contraction within the ventricular septum in hypertrophic cardiomyopathy (HC) may be related to its abnormal morphology because ventricular wall stress is related to wall curvature by the Laplace equation. To test this, we studied 17 HC patients with various septal morphologies using dynamic magnetic resonance imaging techniques. Short- and long-axis curvatures of the basal septal and basal lateral walls were determined on cine images as the reciprocal of the radius of the arc best fit to the endocardial contour, which was negative if the wall was convex to the cavity of the left ventricle. Endocardial and epicardial intramyocardial circumferential shortening (% circumferential shortening) was measured in the septal and lateral walls on basal short-axis myocardial tagging images. Septal walls were flatter in the short-axis plane and more convex toward the left ventricular cavity in the long-axis plane than lateral walls, as indicated by smaller short- and long-axis curvatures. Septal percent circumferential shortening was significantly lower than the lateral percent circumferential shortening, suggesting reduced septal contraction. Endocardial and epicardial percent circumferential shortening showed significant positive correlations with wall curvatures. Multiple stepwise linear regression analysis revealed that both short- and long-axis curvatures significantly contributed to percent circumferential shortening (r=0.87 for endocardial and r=0.70 for epicardial, both p<0.0001). In conclusion, wall curvature is related to wall function in HC; the more convex toward the left ventricular cavity the wall is, the less it contracts. Reduced contraction of the septum in HC may be partly due to its abnormal curvature.

Analysis of left ventricular wall motion based on volumetric deformable models and MRI-SPAMM

We present a new approach for the analysis of the left ventricular shape and motion based on the development of a new class of volumetric deformable models. We estimate the deformation and complex motion of the left ventricle (LV) in terms of a few parameters that are functions and whose values vary locally across the LV. These parameters capture the radial and longitudinal contraction, the axial twisting, and the long-axis deformation. Using Lagrangian dynamics and finite-element theory, we convert these volumetric primitives into dynamic models that deform due to forces exerted by the datapoints. We present experiments where we used magnetic tagging (MRI-SPAMM) to acquire datapoints from the LV during systole. By applying our method to MRI-SPAMM datapoints, we were able to characterize the 3-D shape and motion of the LV both locally and globally, in a clinically useful way. In addition, based on the model parameters we were able to extract quantitative differences between normal and abnormal hearts and visualize them in a way that is useful to physicians.

Cardiovascular applications of magnetic resonance flow and velocity measurements

With recent developments of MR techniques for blood flow measurements, qualitative and quantitative information on both flow volume and flow velocity in the major vessels can be obtained. MR flow quantitation uses the phase, rather than the amplitude of the MR signal, to reconstruct the images. Previous validation studies have demonstrated the accuracy of the phase shift techniques for measuring flow velocities. This technique is now being applied successfully in the cardiovascular system to quantify global and regional ventricular function, valvular heart disease, pulmonary artery disease, thoracic aortic disease, congenital heart disease, and ischemic heart disease.

Semi-automatic tracking of myocardial motion in MR tagged images

Tissue tagging using magnetic resonance (MR) imaging has enabled quantitative noninvasive analysis of motion and deformation in vivo. One method for MR tissue tagging is Spatial Modulation of Magnetization (SPAMM). Manual detection and tracking of tissue tags by visual inspection remains a time-consuming and tedious process. The authors have developed an interactively guided semi-automated method of detecting and tracking tag intersections in cardiac MR images. A template matching approach combined with a novel adaptation of active contour modeling permits rapid analysis of MR images. The authors have validated their technique using MR SPAMM images of a silicone gel phantom with controlled deformations. Average discrepancy between theoretically predicted and semi-automatically selected tag intersections was 0.30 mm+/-0.17 [mean+/-SD, NS (P<0.05)]. Cardiac SPAMM images of normal volunteers and diseased patients also have been evaluated using the authors' technique.

Hybrid DANTE and phase-contrast imaging technique for measurement of three-dimensional myocardial wall motion

Characterization of myocardial stress and strain is necessary for a complete understanding of myocardial function. The precise quantification of regional myocardial strain is complicated by its time-varying pattern and regional variation resulting from the anisotropy of the myocardium and by complex torsional and shortening motions of the heart during the cardiac cycle. The authors have developed a technique for point-specific tracking of myocardial motion along all three axes in a constant selected section of myocardium by combining prospective section selection with in-plane DANTE (delays alternating with nutations for tailored excitation) tissue tagging and phase-contrast detection of motion perpendicular to the image plane. With this technique, it is possible to determine point-specific myocardial strain values in vivo.

Evaluation of cardiac function with magnetic resonance imaging

A large body of evidence has accumulated to substantiate the accuracy of functional MR measurements of both ventricles. Because of good accuracy and superior reproducibility, MR imaging may be considered the gold standard for in vivo quantification of left and right ventricular ejection fraction, myocardial mass, and wall stress. New prospects for functional MR imaging include determination of the end-systolic volume-pressure relation as an index of myocardial contractility. The ability of MR imaging to detect wall motion disturbances may be enhanced further by combining myocardial tagging techniques with finite element analysis. Conventional MR imaging is limited by long examination times, but recent ultrafast modifications of echo-planar imaging allow completion of a functional heart study within seconds. Implementation of ultrafast MR imaging will greatly increase the usefulness of MR imaging for routine evaluation of cardiac function.

Three-dimensional left ventricular deformation in hypertrophic cardiomyopathy

BACKGROUND

In hypertrophic cardiomyopathy, ejection fraction is normal or increased, and force-length relations are reduced. However, three-dimensional (3D) motion and deformation in vivo have not been assessed in this condition. We have reconstructed the 3D motion of the left ventricle (LV) during systole in 7 patients with hypertrophic cardiomyopathy (HCM) and 12 normal volunteers by use of magnetic resonance tagging.

METHODS AND RESULTS

Transmural tagging stripes were automatically tracked to subpixel resolution with an active contour model. A 3D finite-element model was used to interpolate displacement information between short- and long-axis slices and register data on a regional basis. Displacement and strain data were averaged into septal, posterior, lateral, and anterior regions at basal, midventricular, and apical levels. Radial motion (toward the central long axis) decreased slightly in patients with HCM, whereas longitudinal displacement (parallel to the long axis) of the base toward the apex was markedly reduced: 7.5 +/- 2.5mm (SD) versus 12.5 +/- 2.0 mm, P < .001. Circumferential and longitudinal shortening were both reduced in the septum (P < .01 at all levels). The principal strain associated with 3D maximal contraction was slightly depressed in many regions, significantly in the basal septum (-0.18 +/- 0.05 versus -0.22 +/- 0.02, P < .05) and anterior (-0.20 +/- 0.05 versus -0.23 +/- 0.02, P < .05) walls. In contrast, LV torsion (twist of the apex about the long axis relative to the base) was greater in HCM patients (19.9 +/- 2.4 degrees versus 14.6 +/- 2.7 degrees, P < .01).

CONCLUSIONS

HCM patients had reduced 3D myocardial shortening on a regional basis; however, LV torsion was increased.

Assessment of regional left ventricular long-axis motion with MR velocity mapping in healthy subjects

The pattern of left ventricular long-axis motion during early diastole was assessed with magnetic resonance (MR) velocity mapping in 31 healthy volunteers. Regional long-axis velocity varied with time and position around the ventricle. During systole, the base descended toward the apex. The greatest magnitude of long-axis velocity occurred during early diastole. The lateral wall had the highest velocity (140 mm/sec +/- 40 [mean +/- standard deviation]); the anterior and inferior walls had lower velocities (96 mm/sec +/- 27 and 92 mm/sec +/- 34, respectively). The inferoseptal area consistently had the lowest velocities (87 mm/sec +/- 40). Absolute values of peak early-diastolic velocity declined with age (r = -.64, P < .001). Peak early-diastolic velocity was not dependent on heart rate (r = .014, P = .94). Regional variations in left ventricular wall motion were seen. MR velocity mapping is a useful technique for assessing regional left ventricular long-axis heart function.

Two-dimensional left ventricular deformation during systole using magnetic resonance imaging with spatial modulation of magnetization

BACKGROUND

Myocardial tissue tagging with the use of magnetic resonance imaging allows noninvasive regional analysis of heart wall motion and deformation. However, any evaluation of the effect of disease or treatment requires a baseline reference of normal values and variation. We studied the two-dimensional motion of material points imaged within the left ventricular wall using spatial modulation of magnetization (SPAMM) in 12 normal human volunteers.

METHODS AND RESULTS

Five parallel short-axis and five parallel long-axis slices were acquired at five times during systole. SPAMM tags were generated at end diastole using a 7-mm grid. Intersection point data were analyzed for displacement, rotation, and torsion, and triangles of points were analyzed for local rotation and principal strains. Short-axis displacement was the least in the septum for all longitudinal levels (P < .001). Torsion about the long axis was uniform circumferentially because of the motion of the centroids used to reference the rotation. In the long-axis images, the base displaced longitudinally toward the apex, with the posterior wall moving farther than the anterior wall (13.4 +/- 2.2 versus 9.7 +/- 1.8 mm, P < .001) in this direction. The largest principal strain (maximum lengthening) was approximately radially oriented in both views. In the short-axis images, the minimum principal strain (maximum shortening) increased in magnitude toward the apex (P < .001) with little circumferential variation, except at midventricle, where the anterior wall showed greater contraction than the posterior wall (-0.21 +/- 0.03 versus -0.19 +/- 0.02, P < .02).

CONCLUSIONS

Consistent regional variations in deformation are seen in the normal human heart. Displacement and maximum shortening strains are well characterized with two-dimensional magnetic resonance tagging; however, higher-resolution images will be required to study transmural variations.

Automatic tracking of SPAMM grid and the estimation of deformation parameters from cardiac MR images

Presents a new approach for the automatic tracking of SPAMM (Spatial Modulation of Magnetization) grid in cardiac MR images and consequent estimation of deformation parameters. The tracking is utilized to extract grid points from MR images and to establish correspondences between grid points in images taken at consecutive frames. These correspondences are used with a thin plate spline model to establish a mapping from one image to the next. This mapping is then used for motion and deformation estimation. Spatio-temporal tracking of SPAMM grid is achieved by using snakes-active contour models with an associated energy functional. The authors present a minimizing strategy which is suitable for tracking the SPAMM grid. By continuously minimizing their energy functionals, the snakes lock on to and follow the in-slice motion and deformation of the SPAMM grid. The proposed algorithm was tested with excellent results on 123 images (three data sets each a multiple slice 2D, 16 phase Cine study, three data sets each a multiple slice 2D, 13 phase Cine study and three data sets each a multiple slice 2D, 12 phase Cine study).

Magnetic resonance measurement of velocity and flow: technique, validation, and cardiovascular applications

With a newly developed magnetic resonance (MR) technique for blood flow measurements, qualitative and quantitative information on both flow volume and flow velocity in the great vessels can be obtained. MR flow quantitation is performed with a gradient-echo MR sequence with high temporal resolution enabling measurements at frequent intervals throughout the cardiac cycle. MR flow quantitation uses the phase rather than the amplitude of the MR signal to reconstruct the images. These images, often referred to as MR velocity maps or velocity-encoded cine MR images, are two-dimensional displays of flow velocity. From these velocity maps, velocity and volume flow data can be obtained. Previous validation experiments have demonstrated the accuracy of MR velocity mapping, and this technique is now being applied successfully in several clinical fields. MR velocity mapping may be of considerable value when Doppler echocardiography results are unsatisfactory or equivocal, particularly because MR is suited for the analysis of volumetric flow and complex flow patterns. Among the vastly growing number of clinical cardiovascular applications that have been reported are the great arteries and veins, coronary vessels, valvular disease, and the abdominal and peripheral vessels. These items are reviewed, and some aspects of the technique that need improvement are discussed.