MRI reconstruction with motion-dependent regularization

• US 9,482,732 B2
• Filed: 10/29/2013
• Issued: 11/01/2016
• Est. Priority Date: 11/08/2012
• Status: Active Grant

• Current Assignee: Siemens Healthineers AG (Siemens AG)
• Original Assignee: Siemens Corp. (Siemens AG)
• Inventors: Chesneau, Nicolas, Janardhanan, Nirmal, Liu, Jun, Nadar, Mariappan S., Wang, Qiu, Yang, Zhili
• Primary Examiner(s): Kholdebarin, Iman K
• Application Number:
  US14/065,498
• Publication Number:
  US 20140126796A1
• Time in Patent Office: 1,099 Days
• Field of Search: None
• US Class Current: 1/1
• CPC Class Codes: G01R 33/4824 using a non-Cartesian trajectory
• 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/56308 Characterization of motion or flow; Dynamic imaging
• CPC Class Codes: G01R 33/56509 due to motion, displacement or flow, e.g. gradient moment nulling G01R33/567 takes precedence

Abstract

A method of image reconstruction for a magnetic resonance imaging (MRI) system includes obtaining k-space scan data captured by the MRI system, the k-space scan data being representative of an undersampled region over time, iteratively reconstructing preliminary dynamic images for the undersampled region from the k-space scan data via optimization of a first instance of a minimization problem, the minimization problem including a regularization term weighted by a weighting parameter array, generating a motion determination indicative of an extent to which each location of the undersampled region exhibits motion over time based on the preliminary dynamic images, and iteratively reconstructing motion-compensated dynamic images for the region from the k-space scan data via optimization of a second instance of the minimization problem, the second instance having the weighting parameter array altered as a function of the motion determination.

Claims

1. A method of image reconstruction for a magnetic resonance imaging (MRI) system comprising a plurality of coils, the method comprising:

obtaining k-space scan data captured by the MRI system via the plurality of coils, the k-space scan data being representative of an undersampled region over time;

iteratively reconstructing, with a processor, preliminary dynamic images for the undersampled region from the k-space scan data via optimization of a first instance of a minimization problem, the minimization problem comprising a regularization term weighted by a weighting parameter array;

generating a motion determination indicative of an extent to which each location of the undersampled region exhibits motion over time based on the preliminary dynamic images; and

iteratively reconstructing, with the processor, motion-compensated dynamic images for the region from the k-space scan data via optimization of a second instance of the minimization problem, the second instance having the weighting parameter array altered as a function of the motion determination.

2. The method of claim 1, wherein generating the motion determination comprises computing a difference metric for each location of the undersampled region based on the temporal comparison of the preliminary dynamic images with one another.

3. The method of claim 2, wherein values of the weighting parameter array are specified in accordance with whether the difference metric exceeds a threshold.

4. The method of claim 2, wherein computing the difference metric comprises computing a temporal standard deviation for each location of the undersampled region.

5. The method of claim 1, wherein the motion determination is a binary map in which each location is indicated as a motion pixel or as a static pixel if the location exhibits motion or does not exhibit motion, respectively.

6. The method of claim 5, wherein iteratively reconstructing the motion-compensated dynamic images comprises assigning weight elements in the weighting parameter array such that the static pixels have higher weights than the motion pixels.

7. The method of claim 6, wherein assigning the weight elements comprises providing high temporal frequencies with a higher weight relative to low temporal frequencies.

8. The method of claim 1, wherein the weighting parameter array is configured to emphasize temporal correlations in the dynamic images.

9. The method of claim 1, further comprising determining a respective coil sensitivity profile for the region for each coil of the plurality of coils, wherein the minimization problem is a function of the determined coil sensitivity profiles.

10. The method of claim 1, wherein the regularization term is a function of redundant Haar wavelet transforms of the dynamic images.

11. The method of claim 1, wherein the minimization problem is configured to minimize an objective function comprising a data fidelity term configured in accordance with an L2 norm.

12. A computer program product for implementing a method of image reconstruction for a magnetic resonance imaging (MRI) system comprising a plurality of coils, the computer program product comprising one or more non-transitory computer-readable storage media having stored thereon instructions executable by one or more processors of a computing system to cause the computing system to perform operations comprising:

obtaining k-space scan data captured by the MRI system, the k-space scan data being representative of an undersampled region over time;

determining a respective coil sensitivity profile for the region for each coil of the plurality of coils;

iteratively reconstructing preliminary dynamic images for the region from the k-space scan data via optimization of a first instance of a minimization problem, the minimization problem comprising a data fidelity term based on the coil sensitivity profiles and a regularization term weighted relative to the data fidelity term by a weighting parameter array;

comparing the preliminary dynamic images with one another to determine motion and static portions of the undersampled region based on an extent to which each location of the undersampled region exhibits motion over time; and

iteratively reconstructing motion-compensated dynamic images for the region from the k-space scan data via optimization of a second instance of the minimization problem, the second instance having the weighting parameter array adjusted in accordance with the motion and non-motion portions of the undersampled region.

13. The computer program product of claim 12, wherein comparing the preliminary dynamic images comprises computing a difference metric for each location of the undersampled region, the difference metric being indicative of the extent to which each location of the undersampled region exhibits motion over time.

14. The computer program product of claim 13, wherein values of the weighting parameter array are specified in accordance with whether the difference metric exceeds a threshold.

15. The computer program product of claim 12, wherein computing the difference metric comprises computing a temporal standard deviation for each location of the undersampled region.

16. The computer program product of claim 12, wherein comparing the preliminary dynamic images comprises determining a binary map in which each location is indicated as part of the motion portion or as part of the static portion.

17. The computer program product of claim 12, wherein iteratively reconstructing the motion-compensated dynamic images comprises assigning weight elements in the weighting parameter array such that pixels in the static portion have higher weights than pixels in the motion portion.

18. A data processing system for a magnetic resonance imaging (MRI) system comprising a plurality of coils, the data processing system comprising:

a data store in which k-space scan data captured by the plurality of coils is stored, the k-space scan data being representative of an undersampled region over time; and

a processor coupled to the data store, and configured to reconstruct preliminary dynamic images for the undersampled region via optimization of a first instance of a minimization problem, the minimization problem comprising a regularization term weighted by a weighting parameter array;

wherein the processor is further configured to determine a map of the undersampled region indicative of an extent to which each location of the undersampled region exhibits motion over time based on a temporal comparison of the preliminary dynamic images with one another;

wherein the processor is further configured to change the weighting parameter array as a function of the map, and to reconstruct motion-compensated dynamic images for the undersampled region from the k-space scan data via optimization of a second instance of the minimization problem, the second instance using the changed weighting parameter array.

19. The data processing system of claim 18, wherein determining the map comprises computing a temporal standard deviation for each location of the undersampled region based on the temporal comparison of the preliminary dynamic images with one another.

20. The data processing system of claim 19, wherein the map is a binary map in which each location is indicated as a motion pixel if the temporal standard deviation exceeds a threshold, or indicated as a static pixel if the temporal standard deviation does not exceed a threshold.