Imaging of turbulence with magnetic resonance

• US 5,277,192 A
• Filed: 09/18/1992
• Issued: 01/11/1994
• Est. Priority Date: 09/18/1992
• Status: Expired due to Fees

• Current Assignee: General Electric Company
• Original Assignee: General Electric Company
• Inventors: Dumoulin, Charles L
• Primary Examiner(s): Smith, Ruth S.
• Application Number:
  US07/946,702
• Time in Patent Office: 480 Days
• Field of Search: 128/653.2, 128/653.3, 128/708, 324/306, 324/309, 324/312
• US Class Current: 600/413
• CPC Class Codes: A61B 5/0263 using NMR
• CPC Class Codes: A61B 5/7257 using Fourier transforms
• CPC Class Codes: G01R 33/56308 Characterization of motion or flow; Dynamic imaging
• CPC Class Codes: G01R 33/56316 involving phase contrast techniques

Abstract

Turbulent fluid flow within a subject is imaged using a magnetic resonance imaging method which is sensitive to non-coherent changes in the magnetic resonance signal intensity and phase. Flow-induced phase shifts are generated in the magnetic resonance signal with a bi-polar magnetic field gradient pulse applied in an arbitrary orientation. Data is acquired in the presence of a uniform flow encoding field gradient applied in an selected direction independent of the bi-polar field gradient pulse. The flow encoding gradient induces a frequency shift in the data which is proportional to position along the direction of the gradient. The spatial distribution of magnetic resonance spins within the flowing fluids along this applied frequency-encoding magnetic field gradient is determined by Fourier transformation. Data acquisition and Fourier transformation are repeated for a plurality of cardiac cycles of the subject and the standard deviation of the magnetic resonance signal at each position along the applied magnetic field gradient is calculated to produce a one-dimensional projection of the random component of the flow-induced phase shift. The steps of data acquisition, Fourier transformation and standard deviation calculation are then repeated for multiple directions of the frequency-encoding magnetic field gradient. A two-dimensional spatial distribution of the non-coherent component of the magnetic resonance signal is then determined by projection reconstruction of the standard deviation projection data sets.

Claims

1. A method for obtaining magnetic resonance (MR) images of turbulence of fluid within a subject situated in a substantially homogeneous magnetic field comprising the steps of:

a) applying a slice-selective MR excitation pulse to excite a slice of said subject;

b) applying a flow-encoding magnetic field gradient pulse to the excited slice;

c) applying a readout magnetic field gradient pulse to the excited slice at a selected projection angle θ

_m m=1,2,3, . . . M where M represents a total number of projection angles;

d) collecting an MR signal from said subject responsive to the MR excitation and magnetic field gradient pulses;

e) Fourier transforming the MR signal to result in a set of flow velocity values for each position x_i along an axis in a direction of the selected projection angle;

f) repeating steps "a"-"e" during a single cardiac cycle of said subject to obtain a plurality of sets of flow velocity values;

g) repeating steps "a"-"f" for a plurality of N cardiac cycles of said subject to obtain N sets of sets of flow velocity values; and

h) producing a projection p_m being a turbulence vs. position image of fluid within said subject by determining the standard deviation of the sets of flow velocity values for each position x_i of the projection angle θ

_m.

2. The method for obtaining MR images as recited in claim 1 further comprising the steps of:

a) repeating the steps "a"-"h" a plurality M repetitions to arrive at a plurality of M projections p_m, each repetition having the readout magnetic field gradient pulse oriented at a different projection angle θ

_m, m=1,2,3, . . . M; and

b) back projecting the M projections p_m to arrive at a two-dimensional image of fluid turbulence.

3. The method for obtaining MR images as recited in claim 1 in which the flow-encoding magnetic field gradient pulse is a bi-polar pulse.

4. The method for obtaining MR images as recited in claim 1 wherein the flow-encoding magnetic field gradient pulse has more than two lobes and images fluid having motion described by:

##EQU3## where x_i is the displacement along the readout magnetic field gradient pulse direction.

5. The method for obtaining MR images as recited in claim 1 wherein steps "a"-"h" are initiated after detection of a cardiac trigger event in said subject.

6. A method for obtaining magnetic resonance (MR) images of turbulence of fluid within a subject situated in a substantially homogeneous magnetic field comprising the steps of:

a) applying a slice-selective MR excitation pulse to excite a slice of said subject;

b) applying a flow-encoding magnetic field gradient pulse to the excited slice;

c) applying a readout magnetic field gradient pulse to the excited slice at a selected projection angle θ

_m ;

d) collecting a first MR signal from said subject responsive to the rf pulse and magnetic field gradients;

e) Fourier transforming the first MR signal to result in a first set of flow velocity values for each position x_i along an axis in a direction of the selected projection angle;

f) repeating steps "a"-"c";

g) collecting a second MR signal from said subject responsive to the rf pulse and magnetic field gradients;

h) Fourier transforming the second MR signal to result in a second set of flow velocity values;

i) subtracting the first set of flow velocity values from the second set of flow velocity values to obtain a motion encoded signal;

j) repeating steps "a"-"i" during a single cardiac cycle of said subject to obtain a plurality of motion encoded signals; and

k) repeating steps "a"-"j" for a plurality of N cardiac cycles of said subject to obtain N sets of motion encoded signals; and

l) producing a projection p_m being a turbulence vs. position image of fluid within said subject having a zero average magnitude by determining the standard deviation of the motion encoded signals for each position x_i of the projection angle θ

_m.

7. The method for obtaining MR images as recited in claim 6 further comprising the steps of:

a) repeating the steps "a"-"l" a plurality M repetitions to arrive at a plurality of M projections p_m, each repetition having the readout magnetic field gradient pulse oriented at a different projection angle θ

_m, m=1,2,3, . . . M; and

b) back projecting the M projections p_m to arrive at a two-dimensional image of fluid turbulence.