METHOD AND SYSTEM FOR CONSTRAINED RECONSTRUCTION APPLIED TO MAGNETIC RESONANCE TEMPERATURE MAPPING

• US 20080292167A1
• Filed: 01/17/2008
• Published: 11/27/2008
• Est. Priority Date: 05/24/2007
• Status: Abandoned Application

• Current Assignee: University of Utah Research Foundation (State of Utah)
• Original Assignee: University of Utah Research Foundation (State of Utah)
• Inventors: Parker, Dennis L., DiBella, Edward V.R., Raja, Ganesh Sharma Adluru Venkata, Todd, Nick
• Application Number:
  US12/016,071
• Publication Number:
  US 20080292167A1
• Time in Patent Office: Days
• US Class Current: 382/131
• CPC Class Codes: G01R 33/4804 Spatially selective measurement of temperature or pH
• CPC Class Codes: G01R 33/4814 MR combined with ultrasound
• 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/5619 by temporal sharing of data, e.g. keyhole, block regional interpolation scheme for k-Space [BRISK]
• CPC Class Codes: G06T 2207/10088 Magnetic resonance imaging [MRI]
• CPC Class Codes: G06T 7/11 Region-based segmentation

Abstract

A method, system, and computer-readable medium are provided which perform reconstruction of an image from undersampled k-space data. Imaging data of an image area is received. The imaging data is thermal magnetic resonance imaging data in k-space generated at less than the Nyquist rate. A cost function is minimized based on an image estimate and the received imaging data. An image of the image area is defined based on the minimized cost function. The received imaging data may include current k-space data and a summation of k-space data from previous time frames. Additionally, the image may be defined before imaging data is received for a next timeframe.

Claims

1. A system for performing image reconstruction of undersampled image data, the system comprising:

a processor; and

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

receive imaging data of an image area, wherein the imaging data is thermal magnetic resonance imaging data in k-space generated at less than the Nyquist rate;

minimize a cost function based on an image estimate and the received imaging data; and

define an image of the image area based on the minimized cost function.

2. The system of claim 1, further comprising a magnetic resonance imaging machine configured to generate the imaging data of the image area.

3. A computer-readable medium comprising computer-readable instructions therein that, upon execution by a processor, cause the processor to perform image reconstruction of undersampled image data, the instructions configured to cause a computing device to:

receive imaging data of an image area, wherein the imaging data is thermal magnetic resonance imaging data in k-space generated at less than the Nyquist rate;

minimize a cost function based on an image estimate and the received imaging data; and

construct an image of the image area based on the minimized cost function.

4. A method of performing Image reconstruction of undersampled image data, the method comprising:

(a) receiving imaging data of an image area, wherein the imaging data is thermal magnetic resonance imaging data in k-space generated at less than the Nyquist rate;

(b) minimizing a cost function based on an image estimate and the received imaging data; and

(c) defining an image of the image area based on the minimized cost function.

5. The method of claim 4, further comprising:

(d) receiving second imaging data of the image area after receiving the imaging data, wherein the second imaging data is thermal magnetic resonance imaging data in k-space generated at less than the Nyquist rate;

(e) minimizing a second cost function based on an image estimate, the received imaging data, and the received second imaging data; and

(f) defining a second image of the image area based on the minimized second cost function.

6. The method of claim 5, wherein the image is defined before receiving the second imaging data.

7. The method of claim 4, wherein the received imaging data defined for a first time comprises the imaging data obtained for the first time and a summation of the imaging data received at one or more previous times to form a full k-space of the image area.

8. The method of claim 7, repeating (a)-(c) for a second time to support monitoring of a temperature of at least a portion of the image area, wherein the second time occurs after the first time.

9. The method of claim 7, wherein the received imaging data defined for the first time further comprises the imaging data obtained for a second time, wherein the second time occurs after the first time.

10. The method of claim 4, further comprising, before (a);

receiving second imaging data of the image area, wherein the second imaging data is thermal magnetic resonance imaging data in k-space generated at the Nyquist rate; and

initializing the cost function based on the received second imaging data.

11. The method of claim 4, further comprising, before (b), initializing the cost function using a sliding window technique.

12. The method of claim 4, wherein minimizing the cost function comprises applying an iterative gradient descent method.

13. The method of claim 12, wherein the iterative gradient descent method comprises iteratively updating the image data according to {tilde over (m)}^m+1={tilde over (m)}^o−

λ

C⁺( m), where n represents an iteration number, λ

is a step size of a gradient descent method, and C⁺({tilde over (m)}) is an Euler-Lag range derivative of the cost function.

14. The method of claim 4, wherein the cost function comprises a data fidelity term and a constraint term.

15. The method of claim 14, wherein the data fidelity term comprises ∥

WF{tilde over (m)}−

d′

∥

₂² wherein W represents a binary sparsifying pattern used to obtain the received imaging data, F represents a Fourier transform, {tilde over (m)} represents the image estimate, d′

represents the received imaging data, and represents the L2 norm.

16. The method of claim 14, wherein the constraint term comprises α

₁Ψ

({tilde over (m)}), wherein α

₁ is a weighting factor, {tilde over (m)} represents the image estimate, and Ψ

({tilde over (m)}) is a constraint.

17. The method of claim 16, wherein the constraint comprises

${\sum }_{}^{}$



∇

? 

m 

? 

2 2 , 

? 

indicates text missing or illegible when filed

wherein N is the number of pixels in a time frame, ∇

, is a temporal gradient operator, and {tilde over (m)}_i is the i^th pixel of the image estimate over time, and represents the L2 norm.

18. The method of claim 16, wherein the constraint comprises

$\sum$

? 

( ∇

r 

m ~ 

? ) 2 + β

2 

? 

, 

? 

indicates text missing or illegible when filed

wherein N is the number of pixels in a time frame, ∇

, is a temporal gradient operator, and {tilde over (m)}_i is the i^th pixel of the image estimate over time, β

is a small positive constant, and represents the L1 norm.

19. The method of claim 4, further comprising calculating a thermal dose of at least a portion of the image area and presenting the calculated thermal dose to support monitoring of a temperature of at least a portion of the image area.

20. The method of claim 4, further comprising calculating a temperature change of at least a portion of the image area and presenting the calculated temperature change to support monitoring of a temperature of at least a portion of the image area.