Trajectory correction method and apparatus for k-space data in magnetic resonance imaging

Patent metadata

Field	Value
Publication number	US-9594146-B2
Application number	US-201313927376-A
Country	US
Kind code	B2
Filing date	Jun 26, 2013
Priority date	Jun 29, 2012
Publication date	Mar 14, 2017
Grant date	Mar 14, 2017

Abstract

In a trajectory correction method and apparatus for k-space data points in magnetic resonance imaging, a magnetic resonance data acquisition unit is operated to execute a sampling sequence to obtain k-space to be corrected. Empirical points are selected used to divide k-space to be corrected into a central region and a peripheral region. The trajectories of the data points in the central region and the peripheral region are corrected and the sampling sequence is executed again, with the corrected trajectories, to obtain corrected k-space.

First claim

We claim as our invention:

1. A method for acquiring magnetic resonance (MR) data comprising:
from a processor, operating an MR data acquisition scanner to execute an MR data acquisition pulse sequence to read out raw MR data from a subject, and to enter the raw MR data into a memory as k-space data;
from said processor, operating said MR data acquisition scanner in said data acquisition pulse sequence to read out said raw MR data by actuating, in each of three orthogonal directions of a Cartesian coordinate system, a readout gradient magnetic field that comprises a first gradient portion in which a magnitude of said gradient magnetic field changes with respect to time followed by a second portion in which said magnitude of said gradient magnetic field is constant, and entering said raw MR data into said memory as said k-space data with respect to axes in said memory respectively corresponding to said orthogonal directions, along respective actual trajectories that fill actual data entry locations in said memory with said MR raw data;
in said processor, dividing said data entry locations in said memory into a central region of said axes, filled with said raw MR data read out during said first portion of each readout gradient magnetic field, and a peripheral region of said axes, filled with said raw MR data read out during said second portion of each readout gradient magnetic field;
in said processor, automatically calculating, for each of said actual trajectories, a deviation of said data entry locations only in said peripheral region from respective filling locations along predetermined ideal trajectories for filling said memory and, from said deviation, automatically calculating a delay time for activating said readout gradient magnetic fields, relative to respective times at which said readout gradient magnetic fields were activated in said MR data acquisition pulse sequence, that will conform each actual trajectory in said peripheral region to each ideal trajectory;
from said processor, operating said MR data acquisition scanner to acquire further raw MR data with said readout gradient magnetic fields respectively activated with the calculated delay times, and thereby obtaining acquired corrected k-space data only in said peripheral region of said memory and re-acquired MR raw data in said central region of said memory;
in said processor, executing a correction calculation algorithm only for said re-acquired raw MR data in said central region of said memory, thereby obtaining calculated corrected k-space data; and
making said acquired corrected k-space data and said calculated corrected k-space data available from said processor in electronic form as a datafile.

2. The method as claimed in claim 1 , comprising:
from said processor, operating said MR data acquisition scanner in said MR data acquisition pulse sequence to radiate radio-frequency pulses at respectively excite nuclear spins in said subject in an imaging range FOV that comprises first and second layers of said subject that are symmetrical with respect to an isocenter of the MR data acquisition scanner, with the MR raw data at each of said data entry locations having a magnitude and a phase; and
in said processor, calculating the actual data entry locations in each of said trajectories by
identifying the phase ΔΦ 1 of the MR raw data at the respective data entry location of the first layer, and identifying the phase ΔΦ 2 of the raw MR data at a data entry location of the second excitation layer that is symmetric to said respective data entry location of said first excitation layer with said isocenter, and
calculating each actual data entry location k′(t) for all axes according to:

K
′

⁡

(
t
)

=

(

Δ
⁢

⁢

Φ
1

-

ΔΦ
2

)

*
FOV

D
r

*
2
⁢

⁢
π

wherein D r is a distance between the respective layer in which the data entry location is situated and the isocenter.

3. The method as claimed in claim 2 , comprising calculating the delay time by:
calculating, for each of said axes, the data entry locations L desired in each of said ideal trajectories according to:

L
desired

=

γ

2
⁢
π

⁢

∫

G
⁢

ⅆ
t

wherein γ is the gyromagnetic ratio of said nuclear spins, and G is the magnitude of the second portion of the gradient magnetic field activated along the respective axis;
calculating the deviation K shift for each axis according to:
K shift =K′−L desired

wherein K′ is the actual data entry location on the respective axis;
calculating the gradient delay time ΔT delay of said empirical points on the x axis, y axis and z axis according to:
K shift =G·ΔT delay .

4. The method as claimed in claim 1 , comprising, in said correction algorithm for the MR data in said central region
calculating a final data entry location K real (t) for said re-acquired raw MR data on all axes in said central region according to:
K real =K x sin θ cos φ· {right arrow over (x)}+K y sin θ sin φ· {right arrow over (y)}+K z cos θ·{right arrow over ( z )}

wherein K x , K y , K z are respectively corrected data entry locations of the re-acquired raw MR data on the x axis, y axis and z axis in said central region, θ is an angle between a re-acquired trajectory of the re-acquired raw MR data in said central region and the z axis, and φ is an angle between the re-acquired trajectory of the re-acquired raw MR data and the x axis.

5. A magnetic resonance (MR) apparatus comprising:
an MR data acquisition scanner;
a processor in communication with a memory;
said processor being configured to operate said MR data acquisition scanner to execute an MR data acquisition pulse sequence to read out raw MR data from a subject, and to enter the raw MR data into a memory as k-space data;
said processor being configured to operate said MR data acquisition scanner in said data acquisition pulse sequence to read out said raw MR data by

Assignees

• Weng De He
• Zhang Qiong
• Siemens Ag

Inventors

• Weng De He
• Zhang Qiong

Classifications

• G01R35/005 Primary

  Calibrating; Standards or reference devices, e.g. voltage or resistance standards, "golden" references (G01R33/0035, G01R35/002 take precedence) · CPC title

• G01R33/4824

  using a non-Cartesian trajectory · CPC title

• G01R33/4822

  in three dimensions · CPC title

• G01R33/482

  using a Cartesian trajectory · CPC title

• G01R33/565 Primary

  Correction of image distortions, e.g. due to magnetic field inhomogeneities · CPC title