Method and apparatus for accelerating magnetic resonance temperature imaging

• US 7,521,930 B2
• Filed: 03/27/2008
• Issued: 04/21/2009
• Est. Priority Date: 03/29/2007
• Status: Expired due to Fees

• Current Assignee: Siemens Healthcare GMBH (Siemens AG)
• Original Assignee: Siemens AG
• Inventors: Chung, Yiu-Cho, Zhang, Qiang, Li, Guo Bin, Zhou, Xiao Dong
• Primary Examiner(s): Vargas; Dixomara
• Application Number:
  US12/056,607
• Publication Number:
  US 20080238423A1
• Time in Patent Office: 390 Days
• Field of Search: 324300-322, 600410-435
• US Class Current: 324/315
• CPC Class Codes: G01K 2213/00 Temperature mapping
• CPC Class Codes: G01K 7/38 the variations of temperature influencing the magnetic permeability
• CPC Class Codes: G01R 33/4804 Spatially selective measurement of temperature or pH
• CPC Class Codes: G01R 33/4814 MR combined with ultrasound

Abstract

In a method and apparatus for accelerating MR temperature imaging used in MR-monitored high intensity focused ultrasound (HIFU) therapy, temperature changes are determined at the focus of the ultrasound during MR temperature imaging; determining the ideal acceleration rate needed for data sampling according to the temperature changes at said focus is determined, the variable-density (VD) data sampling in k-space is adjusted according to the determined ideal acceleration rate, and the data obtained from sampling are reconstructed to form an image. The capability of accelerating MR temperature imaging with both good temporal resolution and good spatial resolution is improved by determining the acceleration rate according to temperature changes at the ultrasound focus and by adjusting the VD data sampling of k-space and thereby the benefits of good flexibility, feasibility and stability are achieved.

Claims

1. A method for accelerating magnetic resonance (MR) temperature imaging for monitoring high intensity focused ultrasound (HIFU) therapy, comprising the steps of:

conducting a HIFU therapy procedure on a subject situated in an MR imaging apparatus, by emitting high intensity focused ultrasound into the subject from a focus;

operating said MR imaging apparatus in an MR temperature imaging sequence to determine temperature changes at said focus during said HIFU therapy procedure, said MR temperature imaging sequence including acquiring MR temperature imaging data and entering said data into k-space by sampling k-space with a variable density;

automatically electronically determining an ideal acceleration rate for said sampling of k-space dependent on said temperature changes at said focus;

adjusting said variable density sampling of k-space according to the determined ideal acceleration rate; and

reconstructing a temperature image of the subject from the MR temperature imaging data entered into k-space.

2. A method as claimed in claim 1 comprising automatically electronically determining said ideal acceleration rate by determining a higher rate for use as said ideal acceleration rate when said temperature changes quickly at said focus, and determining a lower rate for use as said ideal acceleration rate when the temperature is changing slowly at said focus.

3. A method as claimed in claim 1 comprising executing said MR temperature imaging sequence with phase coding in a phase coding direction, and electronically dividing k-space into a plurality of segments in a same direction as said phase coding direction, and adjusting the variable density sampling of k-space according to the ideal acceleration rate by adjusting the segmentation of k-space for variable density sampling and also adjusting respective acceleration factors for sampling the respective segments, dependent on said ideal acceleration rate.

4. A method as claimed in claim 3 comprising adjusting said segmentation and said acceleration factor by, when the temperature at said focus is changing quickly, making an adjustment selected from the group consisting of reducing respective sizes of segments having lower acceleration factors and increasing the respective acceleration factors for outer segments of k-space and, when the temperature at said focus is changing slowly, making an adjustment selected from the selected from the group consisting of increasing respective sizes of segments having lower acceleration factors and reducing the respective acceleration factors for outer segments of k-space.

5. A method as claimed in claim 3 comprising implementing said MR temperature imaging sequence with phase coding in a phase coding direction, and adjusting the variable density sampling of k-space dependent on the ideal acceleration rate by, if a same ideal acceleration rate is applied for multiple repetitions of said MR temperature imaging sequence, alternatingly sampling phase coding lines in k-space in said multiple repetitions.

6. A method as claimed in claim 1 comprising implementing said MR temperature imaging sequence with frequency coding in a frequency coding direction and phase coding in a phase coding direction, and comprising applying an inverse discrete Fourier transform in the frequency coding direction of k-space and select phase coding lines in k-space containing an area of interest for image reconstruction before reconstructing said temperature image, and reconstructing the temperature image from data in said phase coding lines in k-space in said area of interest.

7. A method as claimed in claim 1 comprising reconstructing said temperature image using an image reconstruction technique selected from the group consisting of echo sharing technique and iterative K-t techniques.

8. A method as claimed in claim 7 comprising using said echo sharing technique to reconstruct said temperature image when real-time updating of said temperature image is required, and using iterative K-t techniques to reconstruct said temperature image when thermal dose estimation is required, and using both said echo sharing technique and an iterative K-t technique to reconstruct said temperature image when both real time updating and thermal dose estimation are required.

9. A method as claimed in claim 1 comprising automatically determining said temperature changes at said focus by re-estimating said temperature changes at said focus, before implementing said HIFU therapy procedure.

10. A method as claimed in claim 1 comprising automatically electronically determining said temperature changes at said focus by periodically feeding back temperature values at said focus obtained from said reconstructed temperature image, and automatically predicting subsequent changes of said temperature at said focus based on said temperature values over a period of time.

11. A method as claimed in claim 10 comprising reducing a period of feeding back said temperature values when said temperature values at said focus exhibit larger changes, and increasing the period of feeding back said temperature values at said focus when said temperature values at said focus exhibit smaller changes.

12. An apparatus for accelerating magnetic resonance (MR) temperature imaging during MR-monitored HIFU therapy wherein an HIFU therapy head emits high intensity focused ultrasound into a subject from a focus, and wherein MR temperature image data are acquired for entry into k-space by variable-density sampling of k-space, said apparatus comprising:

an acceleration rate determining unit that determines temperature changes at said focus during said HIFU therapy and that determines an ideal acceleration rate for sampling k-space dependent on said temperature changes; and

a signal processor that adjusts the variable density sampling of k-space dependent on said ideal acceleration rate, and that reconstructs a temperature image from the MR temperature data entered into k-space.

13. An apparatus as claimed in claim 12 wherein said acceleration rate determining unit comprises a temperature change determining module that determines said temperature changes at said focus, and an acceleration rate determining module that determines said ideal acceleration rate dependent on the temperature changes determined by the temperature change determining module.

14. An apparatus as claimed in claim 12 wherein said signal processing unit comprises a data sampling module that adjusts said variable density sampling of k-space dependent on said ideal acceleration rate, and an image reconstruction module that reconstructs said temperature image from said data in k-space.

15. An apparatus as claimed in claim 12 comprising a temperature feedback unit that periodically feeds back said temperature values at said focus, obtained from said reconstructive temperature image, to said acceleration rate determining unit, and wherein said acceleration rate determining unit determines said temperature change from said temperature values.

16. An apparatus as claimed in claim 15 comprising a temperature feedback control unit that controls a period during which said temperature feedback unit feeds back said temperature values to said acceleration rate determining unit.