Method For Reducing Maximum Local Specific Absorption Rate In Magnetic Resonance Imaging

• US 20100134105A1
• Filed: 10/15/2009
• Published: 06/03/2010
• Est. Priority Date: 10/15/2008
• Status: Active Grant

• Current Assignee: Massachusetts Institute of Technology
• Original Assignee: Massachusetts Institute of Technology
• Inventors: Goyal, Vivek, Zelinski, Adam C., Setsompop, Kawin, Adalsteinsson, Elfar
• Granted Patent: US 8,148,985 B2
• Time in Patent Office: Days
• US Class Current: 324/309
• CPC Class Codes: G01R 33/288 Provisions within MR facilities for enhancing safety during MR, e.g. reduction of the specific absorption rate [SAR], detection of ferromagnetic objects in the scanner room
• CPC Class Codes: G01R 33/5612 Parallel RF transmission, i.e. RF pulse transmission using a plurality of independent transmission channels
• CPC Class Codes: G01R 33/583 Calibration of signal excitation or detection systems, e.g. for optimal RF excitation power or frequency G01R33/246 takes precedence

Abstract

A method for reducing maximum local specific absorption rate (“SAR”) in a magnetic resonance imaging (“MRI”) system is disclosed. More specifically, a plurality of candidate radio frequency (“RF”) pulses are designed and the manner in which they are applied to a subject is determined such that the maximum local SAR is substantially reduced relative to applying the candidate RF pulse that produces the lowest maximum local SAR alone. Put another way, this “time-multiplexing” of a set of RF pulses that each produce approximately the same excitation pattern yields a lower maximum local SAR than does transmitting the individual RF pulse having the lowest local SAR over many repetition times (“TRs”). A convex optimization method is utilized to determine the manner in which the RF pulses are multiplexed in time such that a substantially lower maximum local SAR is achieved.

Claims

1. A method for producing a radio frequency excitation field with an RF coil that forms a part of an MRI system, the steps comprising:

a) estimating an electric and a magnetic field produced by the RF coil;

b) estimating an RF transmission profile of the RF coil using the estimated magnetic field;

c) selecting locations in k-space to which RF energy is to be deposited;

d) determining, using the estimated RF transmission profile and selected locations in k-space, a plurality of candidate RF pulse waveforms corresponding to a plurality of candidate RF pulses that have substantially similar excitation profiles but different specific absorption ratio (SAR) profiles;

e) calculating a local SAR value at each spatial location within a selected portion of the subject for each candidate RF pulse waveform;

f) calculating a plurality of weighting factors, using the determined plurality of candidate RF pulse waveforms and calculated local SAR values for each candidate RF pulse waveform, the plurality of weighting factors defining a number of times to apply RF pulses corresponding to selected ones of the plurality of candidate RF pulse waveforms in order to substantially minimize maximum local SAR in the selected portion of the subject; and

g) directing the MRI system, using the calculated plurality of weighting factors and the associated selected ones of the plurality of candidate RF pulse waveforms, to produce a plurality of RF excitation pulses during a selected time period in order to produce an RF excitation field.

2. The method as recited in claim 1 in which the RF excitation field produced in step g) has a maximum local SAR that is lower than the maximum local SAR produced by applying any one of the plurality of RF excitation pulses alone over the selected time period.

3. The method as recited in claim 1 in which step e) includes calculating an SAR value at each spatial location in a field-of-excitation that contains the subject using the estimated electric field.

4. The method as recited in claim 3 in which step e) includes identifying a cubical region of space that contains substantially N-grams of mass associated with the subject for each spatial location in the selected portion of the subject.

5. The method as recited in claim 4 in which step e) includes calculating an N-gram SAR at each spatial location in the selected portion of the subject using the calculated SAR value and identified cubical region of space corresponding to the associated spatial location.

6. The method as recited in claim 5 in which the N-gram SAR is calculated in step e) by averaging the calculated SAR values for each spatial location in the identified cubical region of space.

7. The method as recited in claim 1 in which step f) includes forming an objective function that relates a maximum local SAR in the selected portion of the subject to the calculated local SAR for each candidate RF pulse waveform and the plurality of weighting factors.

8. The method as recited in claim 7 in which step f) includes minimizing the formed objective function such that the maximum local SAR in the selected portion of the subject is substantially minimized.

9. The method as recited in claim 1 in which the RF coil is an RF coil array that includes a plurality of RF transmission channels.

10. The method as recited in claim 9 in which step a) includes estimating an electric and a magnetic field produced by each of the plurality of RF transmission channels in the RF coil array.

11. The method as recited in claim 10 in which step b) includes estimating an RF transmission profile for each of the plurality of RF transmission channels in the RF coil array using the associated estimated magnetic field.

12. The method as recited in claim 10 in which step e) includes calculating an SAR value at each spatial location in a field-of-excitation that contains the subject using the estimated electric field produced by each of the plurality of RF transmission channels in the RF coil array.

13. The method as recited in claim 12 in which step e) includes:

i) identifying a cubical region of space that contains substantially N-grams of mass associated with the subject for each spatial location in the selected portion of the subject; and

ii) calculating an N-gram SAR at each spatial location in the selected portion of the subject by averaging the calculated SAR values for each spatial location in the identified cubical region of space.

14. The method as recited in claim 1 in which step g) includes producing a timing schedule that defines the timing and order in which each of the selected ones of the plurality of candidate RF pulse waveforms are to be used by the MRI system to produce an RF excitation pulse.

15. The method as recited in claim 1 in which step c) includes forming an objective function using the selected locations in k-space and the estimated RF transmission profile.

16. The method as recited in claim 15 in which step c) includes minimizing the formed objective function using a magnitude least squares estimation.

17. The method as recited in claim 15 in which the objective function in step c) includes a term that penalizes candidate RF pulse waveforms that have an l₂ energy above a selected threshold.