Three dimensional magnetic resonance motion estimation on a single image plane

• US 20070016000A1
• Filed: 05/04/2006
• Published: 01/18/2007
• Est. Priority Date: 05/04/2005
• Status: Active Grant

• Current Assignee: Johns Hopkins University
• Original Assignee: Johns Hopkins University
• Inventors: Stuber, Matthias, Prince, Jerry, Osman, Nael, Abd-Elmoniem, Khaled
• Granted Patent: US 7,495,438 B2
• Time in Patent Office: Days
• US Class Current: 600/410
• CPC Class Codes: A61B 5/055 involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
• CPC Class Codes: G01R 33/56316 involving phase contrast techniques
• CPC Class Codes: G01R 33/56333 Involving spatial modulation of the magnetization within an imaged region, e.g. spatial modulation of magnetization [SPAMM] tagging perfusion imaging based on arterial spin tagging G01R33/56366

Abstract

Three-dimensional MR motion estimation on a single image plane based on tagged MRI and HARP processing. Tagged magnetic resonance imaging technique encodes and automatically tracks displacement of spatially modulated object in three dimensions, encoding both in plane and through-plane motion in a single image plane without affecting acquisition speed. Post-processing unravels encoding in order to directly track 3-D displacement of points within the image plane throughout image sequence. The invention is particularly suited to use on a heart for tracking and determining myocardial displacement. In one embodiment, an MR pulse sequence extends a slice following complementary spatial modulation of magnetization (CSPAMM) pulse sequence with two small z-encoding gradients immediately before the readouts in successive CSPAMM acquisitions, thereby adding a through-plane encoding from which through-plane motion can be computed from acquired images. HARP processing is used to determine in-plane motion, after which through-plane position can be determined using phase encodings. Use of balanced encodings and horizontal and vertical tags permits cancellation of systematic phase artifacts present in CSPAMM acquisitions.

Claims

1. A method of measuring three-dimensional motion of an object by magnetic resonance imaging from one acquired image orientation, comprising:

applying a magnetic resonance imaging pulse sequence to spatially modulate a region of interest of said object, the pulse sequence comprising a slice-following CSPAMM (complementary spatial modulation of magnetization) pulse sequence with added z-phase encode;

acquiring a set of one or more vertically tagged images and a set of one or more horizontally tagged images of said spatially modulated object; and

determining in-plane and through-plane motions of said spatially modulated object from the vertically and horizontally tagged images.

2. The method of claim 1, wherein the pulse sequence comprises a vertical tagging sequence having a first z-encoding gradient and a horizontal tagging sequence having a second z-encoding gradient of equal magnitude but opposite polarity, thereby adding through-plane encoding from which through-plane motion can be computed from the acquired images.

3. The method of claim 2, the determining including computing a first two harmonic phase images for a vertically tagged image and a second two harmonic phase images for a horizontally tagged image.

4. The method of claim 3, the determining further including solving a system of linear equations at a set of image coordinates for phase images indicating in-plane and through-plane motion, the system of linear equations formed by the four harmonic phase images computed for the vertically and horizontally tagged images.

5. The method of claim 4, the determining further including applying harmonic phase (HARP) processing to compute in-plane motion of said spatially modulated object.

6. The method of claim 5, the determining further including iterative application of a wrapping operator to compute through-plane motion of said spatially modulated object.

7. The method of claim 1, including employing said method to obtain a real-time image of cardiac motion.

8. The method of claim 1, including employing said method to obtain a real-time image of myocardial strain.

9. The method of claim 1, including employing said method to obtain a dense three-dimensional motion map of points in a plane for computing myocardial surface strain.

10. The method of claim 1, including employing said method with a plurality of image planes to track a set of points within the myocardium for computing a three-dimensional measure of a regional function and global function.

11. The method of claim 10, wherein the regional function comprises radial, circumferential, or longitudinal strain, principal strains, or the direction of maximal thickening.

12. The method of claim 10, wherein global function comprises longitudinal shortening, rotation, or torsion.

13. The method of claim 1, including employing said method to differentiate between in-plane strain and three-dimensional rotation, further comprising:

using a 3×

3 displacement gradient tensor to remove an apparent false strain due to through-plane rotation; and

computing a planar strain within the imaged slice.

14. A method of measuring three-dimensional motion of an object by magnetic resonance imaging from one acquired image orientation, comprising:

applying a magnetic resonance imaging pulse sequence to spatially modulate a region of interest of said object, the pulse sequence comprising a slice-following CSPAMM (complementary spatial modulation of magnetization) pulse sequence with added z-phase encode;

acquiring a set of one or more of grid tagged images of said spatially modulated object; and

determining in-plane and through-plane motions of said spatially modulated object from the grid tagged images.