Magnetic resonance imaging phase correction from pixel positions method and apparatus

• US 6,466,016 B1
• Filed: 10/10/2000
• Issued: 10/15/2002
• Est. Priority Date: 10/22/1999
• Status: Expired due to Fees

• Current Assignee: GE Medical Systems Global Technology Company LLC A Delaware Limited Liability Company
• Original Assignee: GE Medical Systems Global Technology Company LLC (GE Healthcare Technologies, Inc.)
• Inventors: Miyoshi, Mitsuharu
• Primary Examiner(s): Lefkowitz, Edward
• Assistant Examiner(s): VARGAS, DIXOMARA
• Application Number:
  US09/685,816
• Time in Patent Office: 735 Days
• Field of Search: 324/312, 324/307, 324/309
• US Class Current: 324/312
• CPC Class Codes: G01R 33/565 Correction of image distortions, e.g. due to magnetic field inhomogeneities
• CPC Class Codes: G01R 33/56563 caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0 G01R33/56509, G01R33/56518, G01R33/56536 take precedence

Abstract

To calculate a precise phase map even when there is a local phase disturbance and perform this phase correction, an image acquired by magnetic resonance imaging is subjected to lowpass filtering (702), pixel positions for which the ratio of a value after filtering to a value before filtering does not exceed a predetermined ratio are detected (704), a phase distribution excluding pixel data at these positions is calculated, and the phases at the excluded pixel positions are estimated from the phases of neighboring pixel positions to provide compensation so as to complete a phase map.

Claims

1. A phase correction method comprising the steps of:

lowpass filtering an image obtaind by magnetic resonance imaging;

detecting a range of pixel positions defined by values “

e”

at one end and “

E”

at another end, wherein $\frac{\mathrm{signal}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{intensity}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{after}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{filtering}}{\mathrm{signal}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{intensity}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{before}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{filtering}}\le <br><br>e<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{and}$$\frac{\mathrm{signal}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{intensity}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{before}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{filtering}}{\mathrm{signal}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{intensity}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{after}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{filtering}}\ge <br><br>E$

so as to provide a range of detected pixel positions;

calculating a phase distribution based on the image before said lowpass filtering or the image after said lowpass filtering excluding pixel data within said range of detected pixel positions so as to provide a calculated phase distribution;

obtaining a compensated phase distribution by estimating the phase at said excluded range of detected pixel positions in said calculated phase distribution from the phase at neighboring pixel positions; and

performing phase correction of the image using said compensated phase distribution.

2. A phase correcting apparatus comprising:

filtering means for lowpass filtering an image obtained by magnetic resonance imageing;

detecting means for detecting a range of pixel positions defined by values “

e”

at one end and “

E”

at another end, wherein $\frac{\mathrm{signal}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{intensity}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{after}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{filtering}}{\mathrm{signal}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{intensity}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{before}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{filtering}}\le <br><br>e<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{and}$$\frac{\mathrm{signal}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{intensity}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{before}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{filtering}}{\mathrm{signal}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{intensity}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{after}<br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{filtering}}\ge <br><br>E$

so as to provide a range of detected pixel positions;

calculating means for calculating a phase distribution based on the image before lowpass filtering or the image after said low-pass filtering excluding pixel data within said range of detected pixel positions so as to provide a calculated phase distribution;

compensating means for obtaining a compensated phase distribution by estimating the phase at said excluded range of detected pixel positions in said calculated phase distribution from the phase at neighboring pixel positions; and

correcting means for phase correcting said image using said compensated phase distribution.

3. A phase correction method comprising the steps of:

lowpass filtering an image obtained by magnetic resonance imaging to produce a lowpass filtered image;

calculating a phase distribution based on said lowpass filtered image to provide a calculated phase distribution;

phase correcting the image before said lowpass filtering based on said calculated phase distribution to provide a phase corrected image;

lowpass filtering said phase corrected image to provide a new lowpass filtered image;

calculating a new phase distribution based on said new lowpass filtered image to provide a new calculated phase distribution; and

phase correcting said phase corrected image based on said new calculated phase distribution.

4. The method of claim 3, wherein all of the foregoing steps are repeated cyclically.

5. The method of claim 4, wherein said lowpass filtering is progressively weakened each time said lowpass filtering is provided.

6. The method of claim 4, wherein said lowpass filtering is progressively strengthened each time said lowpass filtering is provided.

7. A phase correction apparatus comprising:

filtering means for lowpass filtering an image obtained by magnetic resonance imaging;

calculating means for calculating a phase distribution based on said lowpass filtered image;

correcting means for phase correcting the image before said lowpass filtering based on said calculated phase distribution; and

control means for causing lowpass filtering of said phase corrected image again by said filtering means, calculating of a new phase distribution by said calculating means based on said new lowpass filtered image, and phase correcting of said phase corrected image again by said correcting means based on said new calculated phase distribution.

8. The apparatus of claim 7, wherein said control means causes cyclical repeated operation of said filtering means, said calculating means, and said correcting means.

9. The apparatus of claim 8, wherein said control means causes said filtering means to provide lowpass filtering which is progressively weakened each time of filtering.

10. The apparatus of claim 8, wherein said control means causes said filtering means to provide lowpass filtering which is progressively strengthened each time of filtering.

11. A magnetic resonance imaging apparatus comprising:

imaging means for acquiring an image of an imaging object using magnetic resonance;

filtering means for lowpass filtering said image obtained by magnetic resonance imaging;

detecting means for detecting pixel positions for which ratio of signal intensity after lowpass filtering to signal intensity before lowpass filtering does not exceed a predetermined value and for providing detected pixel positions;

calculating means for calculating a phase distribution based on said image before said lowpass filtering or said image after said lowpass filtering excluding pixel data at said detected pixel positions;

compensation means for estimating phase at said excluded pixel positions in said calculated phase distribution from phase at neighboring pixel positions to provide compensated phase distribution;

correcting means for phase correcting said image using said compensated phase distribution; and

generating means for generating a water image and a lipid image separately using a phase difference of pixel data of phase corrected image.

12. A magnetic resonance imaging method comprising the steps of:

taking an image of an imaging object using magnetic resonance;

lowpass filtering said image obtained by magnetic resonance imaging;

detecting pixel positions for which ratio of signal intensity before lowpass filtering to signal intensity after lowpass filtering exceeds . a predetermined value and providing detected pixel positions;

calculating a phase distribution based on said image before said lowpass filtering or said image after said lowpass filtering excluding pixel data at said detected pixel positions;

estimating phase at said excluded pixel positions in said calculated phase distribution from pase at neighboring pixel positions to provided compensated phase distribution;

correcting said image using said compensated phase distribution to obtain a phase corrected image; and

generating a water image and a lipid image separately using a phase difference of pixel data of said phase corrected image.

13. A magnetic resoance imaging apparatus comprising:

imaging means for acquiring an image of an imaging object using magnetic resonance imaging;

filtering means for lowpass filtering said image obtained by magnetic resonance imaging to provide a low pass filtered image;

calculating means for calculating a phase distribution based on said lowpass filtered image to provide a calculated phase distribution;

correcting means for phase correcting said image before lowpass filtering based on said calculated phase distribution to provided a phase corrected image;

control means for causing lowpass filtering of said phase corrected image again by said filtering means, calculating of a new phase distribution by said calculating means based on a new lowpass filtered image, and correcting said phase corrected image again by said correcting means based on a new calculated phase distribution; and

generating means for generating a water image and a lipid image separately using a phase difference of pixel data of the phase corrected image.

14. The apparatus of claim 13, wherein said control means causes cyclical operation of said filtering means, said calculating means and said correcting means.

15. The apparatus of claim 14, further comprising means for adjusting said filtering means so that lowpass filtering is progressively weakened for each filtering step.

16. The apparatus of claim 14, further comprising means for adjusting said filtering means so that lowpass filtering is progressively strengthened with each filtering step.

17. A magnetic resonance imaging method comprising the steps of:

taking an image of an imaging object using magnetic resonance;

lowpass filtering said image obtained by magnetic resonance imaging;

calculating a phase distribution based on the lowpass filtered image;

phase correcting the image before lowpass filtering based on the calculated phase distribution obtained in the prior step;

lowpass filtering the phase corrected image obtained in the prior step;

calculating a new phase distribution based on the new lowpass filtered image obtained in the prior step;

phase correcting the phase corrected image again based on the new calculated phase distribution obtained in the prior step; and

generating a water image and a lipid image separately using a phase difference of pixel data of the phase corrected image obtained in the prior step.