Physical motion information capturing of a subject during magnetic resonce imaging automatically motion corrected by the magnetic resonance system

• US 8,653,816 B2
• Filed: 11/04/2009
• Issued: 02/18/2014
• Est. Priority Date: 11/04/2009
• Status: Expired due to Fees

• Current Assignee: International Business Machines Corporation
• Original Assignee: International Business Machines Corporation
• Inventors: Lake, David S., Manduca, Armando, McAllister, Jeffrey S., Ramirez, Nelson
• Primary Examiner(s): Koval, Melissa
• Assistant Examiner(s): FETZNER, TIFFANY A
• Application Number:
  US12/612,172
• Publication Number:
  US 20110101978A1
• Time in Patent Office: 1,567 Days
• Field of Search: 324300-322, 600407-435, 382128-131
• US Class Current: 324/309
• CPC Class Codes: G01R 33/56509 due to motion, displacement or flow, e.g. gradient moment nulling G01R33/567 takes precedence

Abstract

Systems, methods and articles of manufacture are disclosed for compensating for motion of a subject during an MRI scan of the subject. k-space data may be received from the MRI scan of the subject. Motion information may be received for the subject. Based on the received motion information, a translational motion of the subject may be determined between a first point in time and a second point in time. A search space for motion correction may be reduced using the determined change and an error margin of the capturing technique. A motion-compensated, graphical image of the subject may be generated using the reduced search space.

Claims

1. A computer-implemented method in order to compensate for physical motion of a subject during a magnetic resonance imaging (MRI) scan of the subject, based on frequency domain modification and reduction of a search space, the search space itself including a set of potential physical motions of the subject in three-dimensional space of the physical world between a first point in time and a second point in time, the reduction of the search space being based on physical motion information that has been captured by a capturing technique during the MRI scan of the subject, the computer-implemented method comprising:

receiving, from a transmitting element and by operation of one or more computer processors of a system compatible with MRI, k-space data as a result of the MRI scan being performed on the subject,

wherein the k-space data comprises frequency domain data values that are used in generating a graphical image of the subject, and

wherein the frequency domain data values include a plurality of rows including;

a first row of data acquired at the first point in time during the MRI scan and

a second row of data acquired at the second point in time during the MRI scan;

receiving the physical motion information that has been captured by the capturing technique during the MRI scan of the subject, wherein the received physical motion information characterizes physical movement of the subject during the MRI scan;

determining a first physical translational motion of the subject between the first point in time at which the first row of data is acquired during the MRI scan and the second point in time at which the second row of data is acquired during the MRI scan, based on the received physical motion information that characterizes the physical movement of the subject;

programmatically reducing the search space to a set of candidate physical motions falling within a predefined error margin of the determined first physical translational motion, wherein the set of candidate physical motions is a subset of the set of potential physical motions; and

generating a motion-compensated, graphical image of the subject by operation of the one or more computer processors and by modifying one or more of the frequency domain data values of the received k-space data based on the reduced search space, wherein the motion-compensated, graphical image is output for display.

2. The computer-implemented method of claim 1, wherein modifying one or more of the frequency domain data values of the received k-space data based on the reduced search space comprises, for each physical motion falling within the predefined error margin of the determined first physical translational motion:

modifying at least the second row of the frequency domain data values of the k-space data in order to compensate for the respective physical motion and in order to generate modified k-space data;

generating a graphical image of the subject based on the modified k-space data;

programmatically evaluating sharpness of the generated graphical image based on a predefined image sharpness criterion; and

determining whether to accept the modified k-space data based on the evaluation.

3. The computer-implemented method of claim 2, wherein modifying at least the second row of the frequency domain data values of the k-space data comprises:

determining a corresponding change to a value in the k-space data in order to compensate for the respective physical motion, wherein the value is a phase of a sinusoid in the frequency domain representation; and

modifying at least the second row of the frequency domain data values of the k-space data based on the determined change in order to shift an image corresponding to the first row of the frequency domain data values of the k-space data towards an image corresponding to the second row of the frequency domain data values of the k-space data.

4. The computer-implemented method of claim 2, wherein generating a graphical image of the subject comprises performing a two-dimensional fast Fourier transform on the k-space data.

5. The computer-implemented method of claim 3, wherein the corresponding change is determined using the Fourier shift theorem.

6. The computer-implemented method of claim 3, wherein modifying at least the second row of the frequency domain data values of the k-space data comprises multiplying a complex number of the k-space data by another complex number.

7. The computer-implemented method of claim 1, wherein the search space, prior to being reduced, includes a set of all potentially occurring motions of the subject between the first point in time at which the first row of data is acquired during the MRI scan and the second point in time at which the second row of data is acquired during the MRI scan, wherein the predefined error margin is specific to the capturing technique, and wherein at least one physical motion in the set of all potentially occurring physical motions is not included in the set of candidate physical motions.

8. The computer-implemented method of claim 7, wherein the plurality of rows of the k-space data further includes a third row of data acquired at a third point in time during the MRI scan, wherein the first point in time at which the first row of data is acquired during the MRI scan, the second point in time at which the second row of data is acquired during the MRI scan, and the third point in time at which the third row of data is acquired during the MRI scan are distinct, and wherein the method further comprises:

determining a second physical translational motion of the subject between the second point in time at which the second row of data is acquired during the MRI scan and the third point in time at which the third row of data is acquired during the MRI scan, based on the received physical motion information, wherein the determined second physical translational motion is different from the determined first physical translational motion, wherein at least the third row of the frequency domain data values of the k-space data is modified in order to compensate for the determined second physical translational motion, wherein the third row of the frequency domain data values of the k-space data is modified differently than the second row of the frequency domain data values of the k-space data.

9. The computer-implemented method of claim 8, wherein the physical motion is compensated for by a motion correction application, wherein the motion correction application includes:

(i) a k-space manager;

(ii) a motion manager component;

(iii) an interval manager component;

(iv) a k-space modifier component;

(v) an image generator component; and

(vi) an image evaluator component;

wherein the k-space manager component is configured to receive the k-space data from the MRI scan of the subject;

wherein the motion manager component is configured to;

(i) receive the physical motion information captured by the capturing technique during the MRI scan of the subject and

(ii) determine the first physical translational motion of the subject between the first point in time at which the first row of data is acquired during the MRI scan and the second point in time at which the second row of data is acquired during the MRI scan;

wherein the interval manager component is configured to synchronize the physical motion information and the k-space data;

wherein the k-space modifier component is configured to multiply at least one data element in the k-space data by a complex number in order to alter at least one of a phase and an amplitude of the at least one data element;

wherein the image generator component is configured to generate graphical images based on the k-space data; and

wherein the image evaluator component is configured to programmatically evaluate generated graphical images based on a set of predefined criteria, wherein the set of predefined criteria includes the predefined image sharpness criterion.

10. The computer-implemented method of claim 9, wherein modifying the received k-space data based on the reduced search space comprises, for each physical motion falling within the predefined error margin of the determined first physical translational motion:

modifying at least the second row of the frequency domain data values of the k-space data in order to compensate for the respective physical motion;

generating a candidate graphical image of the subject based on the modified k-space data;

programmatically evaluating sharpness of the candidate graphical image based on the predefined image sharpness criterion; and

determining whether to accept the modified k-space data based on the evaluation; and

wherein the modified k-space data is accepted for only one physical motion falling within the predefined error margin of the determined first physical translational motion, the one physical motion having a corresponding candidate graphical image having a highest measure of sharpness among the candidate graphical images for the candidate motions, and wherein only one modification is accepted for each of the determined first physical translational motion and the determined second physical translational motion.

11. The computer-implemented method of claim 10, wherein modifying at least the second row of the frequency domain data values of the k-space data comprises:

determining a corresponding change to a value in the k-space data in order to compensate for the respective physical motion, wherein the value is a phase of a sinusoid in the frequency domain representation; and

modifying at least the second row of the frequency domain data values of the k-space data based on the determined change in order to shift an image corresponding to the first row of the frequency domain data values of the k-space data towards an image corresponding to the second row of the frequency domain data values of the k-space data; and

wherein generating a graphical image of the subject comprises performing a two-dimensional fast Fourier transform on the k-space data, wherein the corresponding change is determined using the Fourier shift theorem, and wherein modifying at least the second row of the frequency domain data values of the k-space data comprises multiplying a data element in the k-space data, comprising a first complex number, by a second complex number in order to alter at least one of a phase and an amplitude of the data element in the k-space data.

12. The computer-implemented method of claim 11, wherein the motion correction application is configured to reduce the search space for motion correction using each of a plurality of distinct capturing techniques including:

(i) a camera-based approach;

(ii) a periodically rotated overlapping parallel lines with enhanced reconstruction (PROPELLER) approach; and

(iii) an optical-fiber-based approach;

wherein the camera-based approach is associated with a first predefined error margin, wherein the PROPELLER approach is associated with a second predefined error margin, wherein the camera-based approach is associated with a third predefined error margin, wherein each of the first predefined error margin, the second predefined error margin, and the third predefined error margin is distinct in magnitude, wherein the motion correction application is configured to reduce the search space for motion correction using the respective error margin specific to each of the plurality of distinct capturing techniques.

13. A non-transitory computer-readable medium containing a program which, when executed by operation of one or more computer processors of a system compatible with magnetic resonance imaging (MRI), performs an operation in order to compensate for physical motion of a subject during an MRI scan of the subject, based on frequency domain modification and reduction of a search space, the search space itself including a set of potential physical motions of the subject in three-dimensional space of the physical world between a first point in time and a second point in time, the reduction of the search space being based on physical motion information that has been captured by a capturing technique during the MRI scan of the subject, the operation comprising:

receiving k-space data as a result of the MRI scan being performed on the subject, wherein the k-space data comprises frequency domain data values that are used in generating a graphical image of the subject, and wherein the frequency domain data values include a plurality of rows including a first row of data acquired at the first point in time during the MRI scan and a second row of data acquired at the second point in time during the MRI scan;

receiving the physical motion information that has been captured by the capturing technique during the MRI scan of the subject, wherein the received physical motion information characterizes physical movement of the subject during the MRI scan;

determining a first physical translational motion of the subject between the first point in time at which the first row of data is acquired during the MRI scan and the second point in time at which the second row of data is acquired during the MRI scan, based on the received physical motion information that characterizes the physical movement of the subject;

programmatically reducing the search space to a set of candidate physical motions falling within a predefined error margin of the determined first physical translational motion, wherein the set of candidate physical motions is a subset of the set of potential physical motions; and

generating a motion-compensated, graphical image of the subject by operation of the one or more computer processors when executing the program and by modifying one or more of the frequency domain data values of the received k-space data based on the reduced search space, wherein the motion-compensated, graphical image is output for display.

14. The non-transitory computer-readable medium of claim 13, wherein modifying one or more of the frequency domain data values of the received k-space data based on the reduced search space comprises, for each physical motion falling within the predefined error margin of the determined first physical translational motion:

modifying at least the second row of the frequency domain data values of the k-space data in order to compensate for the respective physical motion and in order to generate modified k-space data;

generating a graphical image of the subject based on the modified k-space data;

programmatically evaluating sharpness of the generated graphical image based on a predefined image sharpness criterion; and

determining whether to accept the modified k-space data based on the evaluation.

15. The non-transitory computer-readable medium of claim 14, wherein modifying at least the second row of the frequency domain data values of the k-space data comprises:

determining a corresponding change to a value in the k-space data in order to compensate for the respective physical motion, wherein the value is a phase of a sinusoid in the frequency domain representation; and

modifying at least the second row of the frequency domain data values of the k-space data based on the determined change in order to shift an image corresponding to the first row of the frequency domain data values of the k-space data towards an image corresponding to the second row of the frequency domain data values of the k-space data.

16. The non-transitory computer-readable medium of claim 14, wherein generating a graphical image of the subject comprises performing a two-dimensional fast Fourier transform on the k-space data.

17. The non-transitory computer-readable medium of claim 15, wherein the corresponding change is determined using the Fourier shift theorem.

18. The non-transitory computer-readable medium of claim 15, wherein modifying at least the second row of the frequency domain data values of the k-space data comprises multiplying a complex number of the k-space data by another complex number.

19. A system, compatible with magnetic resonance imaging (MRI), that compensates for physical motion of a subject during an MRI scan of the subject, based on frequency domain modification and reduction of a search space, the search space itself including a set of potential physical motions of the subject in three-dimensional space of the physical world between a first point in time and a second point in time, the reduction of the search space being based on physical motion information that has been captured by a capturing technique during the MRI scan of the subject, the system comprising:

a processor; and

a memory containing a program, which when executed by the processor is configured to perform an operation comprising;

receiving k-space data as a result of the MRI scan being performed on the subject, wherein the k-space data comprises frequency domain data values that are used in generating a graphical image of the subject, and wherein the frequency domain data values include a plurality of rows including a first row of data acquired at the first point in time during the MRI scan and a second row of data acquired at the second point in time during the MRI scan;

receiving the physical motion information that has been captured by the capturing technique during the MRI scan of the subject, wherein the received physical motion information characterizes physical movement of the subject during the MRI scan;

determining a first physical translational motion of the subject between the first point in time at which the first row of data is acquired during the MRI scan and the second point in time at which the second row of data is acquired during the MRI scan, based on the received physical motion information that characterizes the physical movement of the subject;

programmatically reducing the search space to a set of candidate physical motions falling within a predefined error margin of the determined first physical translational motion, wherein the set of candidate physical motions is a subset of the set of potential physical motions; and

generating a motion-compensated, graphical image of the subject by modifying one or more of the frequency domain data values of the received k-space data based on the reduced search space, wherein the motion-compensated, graphical image is output for display.

20. The system of claim 19, wherein modifying one or more of the frequency domain data values of the received k-space data based on the reduced search space comprises, for each physical motion falling within the predefined error margin of the determined first physical translational motion:

modifying at least the second row of the frequency domain data values of the k-space data in order to compensate for the respective physical motion and in order to generate modified k-space data;

generating a graphical image of the subject based on the modified k-space data;

programmatically evaluating sharpness of the generated graphical image based on a predefined image sharpness criterion; and

determining whether to accept the modified k-space data based on the evaluation.

21. The system of claim 20, wherein modifying at least the second row of the frequency domain data values of the k-space data comprises:

determining a corresponding change to a value in the k-space data in order to compensate for the respective physical motion, wherein the value is a phase of a sinusoid in the frequency domain representation; and

modifying at least the second row of the frequency domain data values of the k-space data based on the determined change in order to shift an image corresponding to the first row of the frequency domain data values of the k-space data towards an image corresponding to the second row of the frequency domain data values of the k-space data.

22. The system of claim 20, wherein generating a graphical image of the subject comprises performing a two-dimensional fast Fourier transform on the k-space data.

23. The system of claim 21, wherein the corresponding change is determined using the Fourier shift theorem.

24. The system of claim 21, wherein modifying at least the second row of the frequency domain data values of the k-space data comprises multiplying a complex number of the k-space data by another complex number.