Method for correcting image distortion and system, and magnetic resonance imaging equipment

Patent metadata

Field	Value
Publication number	US-9404987-B2
Application number	US-201313872292-A
Country	US
Kind code	B2
Filing date	Apr 29, 2013
Priority date	Apr 27, 2012
Publication date	Aug 2, 2016
Grant date	Aug 2, 2016

Abstract

In a method for an apparatus correcting image distortion in diffusion-weighted echo-planar magnetic resonance imaging, a marker sequence is applied before a diffusion-weighted echo planar imaging sequence, to form a combined sequence. The combined sequence is used to obtain marked images with different preset b values and different preset diffusion directions. The diffusion-weighted echo planar imaging sequence is used to obtain diffusion-weighted echo planar images with the same b values and diffusion directions as the marked images. A stretching coefficient and a displacement coefficient are calculated for each image data column of the diffusion-weighted echo planar image. The stretching coefficient and displacement coefficient are used to correct the diffusion-weighted echo planar images.

First claim

I claim as my invention:

1. A method for correcting image distortion in diffusion-weighted echo planar imaging, comprising:
utilizing a combined sequence for acquisition to obtain a plurality of marked images with different preset b values and different preset diffusion directions, wherein the combined sequence comprises a diffusion-weighted echo planar imaging sequence and a marker sequence applied before the same; and utilizing the diffusion-weighted echo planar imaging sequence for acquisition to obtain a plurality of diffusion-weighted echo planar images with the same b values and diffusion directions as the marked images;
subtracting each of the marked images from the diffusion-weighted echo planar image having the same b values and the same diffusion directions to obtain a marker image corresponding to each of the diffusion-weighted echo planar images, wherein a marker image with the b value being zero is a reference marker image, and a marker image with the b value not being zero is a distortion marker image; and
for each diffusion-weighted echo planar image with the b value not being zero, utilizing the distortion marker image having the same b value and the same diffusion direction as the diffusion-weighted echo planar image and a reference marker image corresponding to the diffusion-weighted echo planar image to calculate a stretching coefficient and a displacement coefficient corresponding to each image data column of the diffusion-weighted echo planar image, and utilizing the stretching coefficient and displacement coefficient to correct the diffusion-weighted echo planar image.

2. The method as claimed in claim 1 , wherein calculating a stretching coefficient and a displacement coefficient corresponding to each image data column of the diffusion-weighted echo planar image comprises:
selecting a plurality of distorted image data columns from the distortion marker image, and for each of the distorted image data columns selected, selecting a reference image data column in a corresponding position from the reference marker image; or, selecting a plurality of reference image data columns from the reference marker image, and for each of the selected reference image data columns, selecting a distorted image data column in a corresponding position from the distortion marker image;
for the currently selected reference image data column and the corresponding distorted image data column thereof, determining the position of a wave peak immediately adjacent to one side of the highest peak in a Fourier transform spectrum curve of the reference image data column, performing position fitting on the wave peak, and taking the peak value position of the fitted wave peak as a reference marker frequency; and determining the position of a wave peak immediately adjacent to one side of the highest peak in the Fourier transform spectrum curve of the distorted image data column, performing position fitting on the wave peak, and taking the peak value position of the fitted wave peak as a distortion marker frequency;
taking the ratio between the distortion marker frequency and the reference marker frequency as the stretching coefficient of the current distorted image data column relative to the current reference image data column;
utilizing the stretching coefficient to correct the distorted image data column, successively displacing the corrected distorted image data column in a preset direction according to a preset displacement interval, calculating a relevance value between the distorted image data column and the reference image data column after each displacement, and taking the displacement amount corresponding to the maximum relevance value as the displacement coefficient of the distorted image data column relative to the reference image data column; and
after calculation of the stretching coefficients and displacement coefficients corresponding to a plurality of distorted image data columns, respectively performing linear fitting and extrapolation on a plurality of stretching coefficients and a plurality of displacement coefficients to obtain the stretching coefficient and the displacement coefficient corresponding to each distorted image data column of the entire distortion marker image, as the stretching coefficient and the displacement coefficient corresponding to each image data column of the diffusion-weighted echo planar image.

3. The method as claimed in claim 2 , wherein:
selecting a plurality of distorted image data columns from the distortion marker image comprises: successively extracting each image data column in the distortion marker image, for the currently extracted image data column, calculating a Fourier transform spectrum curve of the image data column, and if the ratio between the peak value of the highest peak and the peak value of a wave peak on one side of the highest peak in the spectrum curve is greater than a preset selection threshold value, then selecting the currently extracted image data column as a distorted image data column; and
selecting a plurality of reference image data columns from the reference marker image comprises: successively extracting each image data column in the reference marker image, for the currently extracted image data column, calculating a Fourier transform spectrum curve of the image data column, and if the ratio between the peak value of the highest peak and the peak value of a wave peak on one side of the highest peak in the spectrum curve is greater than a preset selection threshold value, then selecting the currently extracted image data column as a reference image data column.

4. A magnetic resonance system for correcting image distortion in diffusion-weighted echo planar imaging, comprising:
a magnetic resonance data image acquisition unit for utilizing a combined sequence for acquisition to obtain a plurality of marked images with different preset b values and different preset diffusion directions, wherein the combined sequence comprises a diffusion-weighted echo planar imaging sequence and a marker sequence applied before the same; and utilizing the diffusion-weighted echo planar imaging sequence for acquisition to obtain a plurality of diffusion-weighted echo planar images with the same b values and diffusion directions as the marked images;
a marker image generation unit for subtracting each of the marked images from the diffusion-weighted echo planar images having the same b values and the same diffusion directions to obtain a marker image corresponding to each of the diffusion-weighted echo planar images, wherein a marker image with the b value being zero is a reference marker image, and a marked image with the b value not being zero is a distortion marker image;
a distortion coefficient calculation unit for, regarding a current diffusion-weighted echo planar image in the diffusion-weighted echo planar images with the b value not being zero, utilizing a distortion marker image having the same b value and the same diffusion direction as the current diffusion-weighted echo planar image and a reference marker image corresponding to the current diffusion-weighted echo planar image to calculate a stretching coefficient and a displacement coefficient corresponding to each image data column of the current diffusion-weighted echo planar image; and
an image correction unit for utilizing the stretching coefficient and displacement coefficient corresponding to each image data column of the current diffusion-weighted echo planar image obtained by the distortion coefficient calculation unit to correct the current diffusion-weighted echo planar image.

5. The system as claimed in claim 4 , wherein the distortion coefficient calculation unit comprises:
an image data column selection subunit for selecting a plurality of distorted image data columns from the distortion marker imag

Assignees

• Zhou Kun
• Siemens Ag

Inventors

• Zhou Kun

Classifications

• G01R33/56518

  due to eddy currents, e.g. caused by switching of the gradient magnetic field · CPC title

• G01R33/56572 Primary

  caused by a distortion of a gradient magnetic field, e.g. non-linearity of a gradient magnetic field (G01R33/56509, G01R33/56518, G01R33/56536 take precedence) · CPC title

• G01R33/56341 Primary

  Diffusion imaging · CPC title