Method and apparatus for accelerating magnetic resonance imaging

• US 9,476,957 B2
• Filed: 07/23/2013
• Issued: 10/25/2016
• Est. Priority Date: 08/06/2012
• Status: Active Grant

• Current Assignee: Universitaetsklinikum Freiburg (Albert-Ludwigs-Universitaet Freiburg)
• Original Assignee: Universitaetsklinikum Freiburg (Albert-Ludwigs-Universitaet Freiburg)
• Inventors: Li, Guobin, Lee, Kuan Jin, Zaitsev, Maxim, Hennig, Juergen, Dragonu, Julius
• Primary Examiner(s): Gray, David
• Assistant Examiner(s): ROTH, LAURA K
• Application Number:
  US13/948,191
• Publication Number:
  US 20140035576A1
• Time in Patent Office: 1,190 Days
• US Class Current: 1/1
• CPC Class Codes: G01R 33/5608 Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels image data processing or generation, in general G06T
• CPC Class Codes: G01R 33/5611 Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE structural details of arrays of sub-coils G01R33/3415
• CPC Class Codes: G01R 33/5619 by temporal sharing of data, e.g. keyhole, block regional interpolation scheme for k-Space [BRISK]
• CPC Class Codes: G01R 33/56509 due to motion, displacement or flow, e.g. gradient moment nulling G01R33/567 takes precedence
• 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
• CPC Class Codes: G01R 33/5659 caused by a distortion of the RF magnetic field, e.g. spatial inhomogeneities of the RF magnetic field G01R33/56509, G01R33/56518, G01R33/56536 take precedence

Abstract

A method of MRI entails recording and storing MR signals from an object to produce an old set of data. A further measurement of the object is then initiated at a later time to record new MR data, whereby k-space is undersampled in the further measurement. The old data are corrected for changes in the new position of the object, for changes in the sensitivity and exact spatial positioning of the receiver coils as well as for changes in the actual field shimming. The old and new data are then combined to create a new, high resolution image of the object.

Claims

1. A method of magnetic resonance (MR) imaging with spatial encoding for generating an image of an object of investigation inside an MR apparatus, the method comprising the steps of:

a) selecting a volume of interest within the object of investigation;

b) applying a radio frequency (RF) pulse to generate a transverse magnetization within the object of investigation;

c) effecting spatial encoding through application of spatially encoding magnetic fields;

d) recording and storing MR signals originating from the object of investigation and assigning those signals to old MR data;

e) physically removing the object of investigation from the MR apparatus following step d);

f) repositioning the object inside the MR apparatus following step e);

g) repeating steps a) through d) for recording further MR signals originating from the same object of investigation to record and store new MR data, whereby step c) is performed by undersampling the k-space for shortening a measurement time;

h) correcting the old MR data recorded in step d) for changes in a new position of the object of investigation relative to the MR apparatus in step g);

i) correcting the old MR data recorded in step d) for changes in a sensitivity and exact spatial positioning of RF receiver coils of the MR apparatus in step g);

j) correcting old MR data recorded in step d) for changes in an actual B0 magnetic field of the MR apparatus in step g); and

k) reconstructing a high resolution MR image of the object of investigation from a synthesis of the old MR data recorded and stored in step d) and corrected in steps h), i) and j) with the new MR data recorded in step g).

2. The method of claim 1, wherein a coarse new image is reconstructed based solely on the new MR data recorded in step g), and the coarse new image is used to perform the corrections of steps h), i) and j).

3. The method of claim 2, wherein the corrections in step h) for compensating changes in a new position of the object of investigation relative to the MR apparatus comprise aligning a prior image reconstructed in an old examination to the coarse new image using an image registration method and optimizing an image alignment locally in image space using an iterative procedure.

4. The method of claim 3, wherein information acquired in the old examination and information acquired in a new examination are adaptively synthesized by extracting an image profile of an aligned prior image using a low-pass filter or by extracting the image profile of the aligned prior image using a fitting algorithm and extracting an image profile of the coarse new image using that fitting algorithm.

5. The method of claim 4, wherein the aligned prior image is replaced by multiplying the aligned prior image by the image profile of the coarse new image and dividing by the image profile of the aligned prior image, to obtain an estimated image.

6. The method of claim 4, wherein a similarity map is calculated in coarse scale by subtracting the estimated image from the coarse new image to obtain a difference image, applying a low-pass filter to the difference image to obtain a filtered difference image, inversing a signal intensity of the filtered difference image to obtain the similarity map and using this similarity map to guide a data synthesis step to obtain a coarsely synthesized image, wherein more information is taken from the estimated image in regions with high similarity.

7. The method of claim 6, wherein the coarsely synthesized image is transformed into k-space and undersampled with an identical sampling pattern applied to a new examination to obtain an estimated k-space.

8. The method of claim 7, where the estimated k-space is subtracted from a new acquired k-space to obtain a residual difference k-space and the residual difference k-space is reconstructed with a reconstruction algorithm to obtain a weighting map, wherein the coarsely synthesized image is combined with the coarse new image using the weighting map to obtain a final reconstructed image.

9. The method of claim 1, wherein, in step g), a same MR apparatus is used as in step d), however with different RF coils or wherein, in step g), an MR apparatus different from that in step d) is used.

10. The method of claim 1, wherein one of a conventional MRI sequence, a gradient echo, a spin echo, a fast low angle shot (FLASH), a fast spin echo (FSE) or an echo planar imaging (EPI) is used for generating image data.

11. The method of claim 1, wherein acceleration techniques or acceleration techniques with k-space (GRAPPA) or image space (SENSE) reconstructions are used.

12. The method of claim 1, wherein a compressed sensing technique is used.

13. An apparatus for accelerated recording of MR images using the method of claim 1, the apparatus comprising:

a database unit built to store patient information including acquired MR data, reconstructed images and used scanning parameters in historical examination;

an MR scanner unit reading information from said database unit for configuring scanning parameters in a new examination; and

an image reconstruction unit loading information from said database unit and said MR scanner and integrating prior information into reconstruction of a new image, wherein any new image reconstruction updates said database unit after each examination.