Method and system for motion correction in imaging systems

• US 8,155,389 B2
• Filed: 08/12/2008
• Issued: 04/10/2012
• Est. Priority Date: 10/02/2007
• Status: Active Grant

• Current Assignee: University of Utah Research Foundation (State of Utah)
• Original Assignee: University of Utah Research Foundation (State of Utah)
• Inventors: Parker, Dennis Lee, Kholmovski, Evgueni G.
• Primary Examiner(s): Lu, Tom Y
• Application Number:
  US12/190,287
• Publication Number:
  US 20090087057A1
• Time in Patent Office: 1,337 Days
• Field of Search: 382/107, 382/128, 382/131, 382/132, 324/307, 324/309
• US Class Current: 382/107
• CPC Class Codes: G01R 33/561 by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
• 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/56509 due to motion, displacement or flow, e.g. gradient moment nulling G01R33/567 takes precedence
• CPC Class Codes: G06T 2207/10088 Magnetic resonance imaging [MRI]
• CPC Class Codes: G06T 2207/20056 Discrete and fast Fourier transform, [DFT, FFT]
• CPC Class Codes: G06T 2207/20201 Motion blur correction
• CPC Class Codes: G06T 2207/30004 Biomedical image processing
• CPC Class Codes: G06T 5/10 using non-spatial domain filtering
• CPC Class Codes: G06T 5/50 using two or more images, e.g. averaging or subtraction
• CPC Class Codes: G06T 5/73 Deblurring; Sharpening
• CPC Class Codes: G06T 7/20 Analysis of motion motion estimation for coding, decoding, compressing or decompressing digital video signals H04N19/43, H04N19/51

Abstract

A method, a system, and a computer-readable medium are provided which perform motion correction of image data. A first set of data and a second set of data of k-space data of an object to be imaged are received. The first set of data and the second set of data include a plurality of phase encoded lines that encompass the object to be imaged. The first set of data correlates with the second set of data. A cross-correlation is calculated by multiplying the first set of data and the second set of data in k-space. A value of the motion of the object to be imaged that occurred between acquisition of the first set of data and acquisition of the second set of data is calculated using the cross-correlation. The second set of data is corrected using the calculated value to remove the motion. The correction process is repeated until the k-space data is completely processed.

Claims

1. A system for performing motion correction of image data, the system comprising:

a processor; and

a computer-readable medium operably coupled to the processor, the computer-readable medium comprising instructions that, when executed by the processor, perform operations comprising

receiving a first set of data that is a first subset of k-space data of an object to be imaged, the first set of data including a first plurality of phase encoded lines that encompass the object to be imaged;

receiving a second set of data that is a second subset of the k-space data of the object to be imaged, the second set of data including a second plurality of phase encoded lines that encompass the object to be imaged, wherein the first set of data correlates with the second set of data and the second set of data is shifted in a phase-encode direction relative to the first set of data;

calculating a cross-correlation by multiplying the received first set of data and the received second set of data in k-space;

calculating a value of the motion of the object to be imaged that occurred between acquisition of the first set of data and acquisition of the second set of data using the calculated cross-correlation; and

correcting the second set of data to remove the motion using the calculated value.

2. The system of claim 1, further comprising an imaging apparatus configured to generate the first set of data and to generate the second set of data.

3. The system of claim 2, further comprising a communication interface operably coupled to the processor to receive the first set of data and the second set of data from the imaging apparatus.

4. A non-transitory computer-readable medium comprising computer-readable instructions therein that, when executed by a processor, cause a computing device to:

receive a first set of data that is a first subset of k-space data of an object to be imaged, the first set of data including a first plurality of phase encoded lines that encompass the object to be imaged;

receive a second set of data that is a second subset of the k-space data of the object to be imaged, the second set of data including a second plurality of phase encoded lines that encompass the object to be imaged, wherein the first set of data correlates with the second set of data and the second set of data is shifted in a phase-encode direction relative to the first set of data;

calculate a cross-correlation by multiplying the received first set of data and the received second set of data in k-space;

calculate a value of the motion of the object to be imaged that occurred between acquisition of the first set of data and acquisition of the second set of data using the calculated cross-correlation; and

correct the second set of data to remove the motion using the calculated value.

6. The non-transitory computer-readable medium of claim 4, wherein the computer-readable instructions further cause the computing device to repeat (a)-(e) with the second set of data as the first set of data and with a third set of data that is a third subset of the k-space data as the second set of data until each subset of the k-space data is processed.

7. The non-transitory computer-readable medium of claim 4, wherein the computer-readable instructions further cause the computing device to repeat (a)-(e) with a third set of data that is a third subset of the k-space data as the second set of data until each subset of the k-space data is processed.

8. The non-transitory computer-readable medium of claim 4, wherein the k-space data provides a field of view that is at least twice as large as the object to be imaged in a phase-encode direction.

9. The non-transitory computer-readable medium of claim 4, wherein the computer-readable instructions further cause the computing device to define an image of the object using the first set of data and the corrected second set of data by calculating an inverse Fourier transform with the first set of data in combination with the corrected second set of data.

10. The non-transitory computer-readable medium of claim 9, wherein the computer-readable instructions further cause the computing device to output the defined image of the object wherein outputting the defined image comprises at least one of printing the defined image, displaying the defined image, and storing data associated with the defined image in a computer-readable medium.

11. The non-transitory computer-readable medium of claim 4, wherein the k-space data provides a field of view that is larger than the object to be imaged.

12. The non-transitory computer-readable medium of claim 11, wherein the k-space data provides a field of view that is at least twice as large as the object to be imaged in a phase-encode direction.

13. The non-transitory computer-readable medium of claim 4, wherein the first set of data is uniformly sampled in a phase-encode direction.

14. The non-transitory computer-readable medium of claim 4, wherein the first set of data is randomly sampled in a phase-encode direction.

15. The non-transitory computer-readable medium of claim 4, wherein the first set of data does not overlap with the second set of data in a phase-encode direction.

16. The non-transitory computer-readable medium of claim 4, wherein the first set of data is obtained in a first echotrain using a spin echo sequence.

17. The non-transitory computer-readable medium of claim 16, wherein the first echotrain includes a plurality of echoes.

18. The non-transitory computer-readable medium of claim 16, wherein the second set of data is obtained in a second echotrain using the spin echo sequence.

19. The non-transitory computer-readable medium of claim 4, wherein the calculated value is a linear phase shift and correcting the second set of data comprises applying the linear phase shift to the received second set of data.

5. A method of performing motion correction of image data, the method comprising:

(a) receiving a first set of data that is a first subset of k-space data of an object to be imaged, the first set of data including a first plurality of phase encoded lines that encompass the object to be imaged;

(b) receiving a second set of data that is a second subset of the k-space data of the object to be imaged, the second set of data including a second plurality of phase encoded lines that encompass the object to be imaged, wherein the first set of data correlates with the second set of data and the second set of data is shifted in a phase-encode direction relative to the first set of data;

(c) calculating a cross-correlation by multiplying the received first set of data and the received second set of data in k-space;

(d) calculating a value of the motion of the object to be imaged that occurred between acquisition of the first set of data and acquisition of the second set of data using the calculated cross-correlation; and

(e) correcting the second set of data to remove the motion using the calculated value.