B1 mapping in MRI system using k-space spatial frequency domain filtering

• US 20100239142A1
• Filed: 03/19/2009
• Published: 09/23/2010
• Est. Priority Date: 03/19/2009
• Status: Active Grant

• Current Assignee: Toshiba Medical Systems Corporation (Canon Inc.)
• Original Assignee: Toshiba Medical Systems Corporation (Canon Inc.), Kabushiki Kaisha Toshiba (Toshiba Corporation)
• Inventors: Dannels, Wayne R., Wheaton, Andrew J.
• Granted Patent: US 8,077,955 B2
• Time in Patent Office: Days
• US Class Current: 382/131
• CPC Class Codes: G01R 33/246 Spatial mapping of the RF magnetic field B1
• CPC Class Codes: G01R 33/443 Assessment of an electric or a magnetic field, e.g. spatial mapping, determination of a B0 drift or dosimetry
• CPC Class Codes: G01R 33/56333 Involving spatial modulation of the magnetization within an imaged region, e.g. spatial modulation of magnetization [SPAMM] tagging perfusion imaging based on arterial spin tagging G01R33/56366
• CPC Class Codes: G06T 11/003 Reconstruction from projections, e.g. tomography

Abstract

Frequency filtering of spatially modulated or “tagged” MRI data in the spatial frequency k-space domain with subsequent 2DFT to the spatial domain and pixel-by-pixel arithmetic calculations provide robust ratio values that can be subjected to inverse trigonometric functions to derive B1 maps for an MRI system.

Claims

1. A method for generating a B1 map in an MRI system, said method comprising using an MRI system to:

(a) acquire at least one set of B1 amplitude tagged MRI data signals in spatial frequency domain k-space from MR nuclei within an imaged volume of the MRI system;

(b) process said k-space data to produce at least two sub-sets of frequency-filtered k-space data;

(i) a baseline low-frequency sub-set, and (ii) at least one of a positive higher frequency sub-set including a harmonic version of said baseline sub-set and a negative higher frequency sub-set including a harmonic version of said baseline sub-set;

(c) separately transform each of said sub-sets to respectively corresponding spatial domain data sets;

(i) a baseline magnitude set, and (ii) at least one of a positive harmonic magnitude set and a negative harmonic magnitude set;

(d) arithmetically combine said baseline magnitude set with said at least one of the harmonic magnitude sets on a pixel-by-pixel basis to provide upper and lower magnitude data sets;

(e) process said upper and lower magnitude data set values on a pixel-by-pixel basis to generate a B1 map value based on value related to a ratio of said upper and lower pixel data values at a given pixel location in the spatial domain; and

(f) store said B1 map of values in a computer accessible and readable memory device for subsequent use by said MRI system.

2. The method of claim 1, wherein step (a) comprises:

applying at least two RF excitation pulses with at least one interleaved magnetic gradient pulse pre-sequence before applying an MR imaging pulse sequence to elicit successive NMR responses in the time domain which are mapped to respectively corresponding portions of a two-dimensional set of k-space data in a spatial frequency domain.

3. The method of claim 2, wherein step (a) further comprises:

applying a 2DFT (two-dimensional Fourier transform) to an originally-acquired set of two-dimensional k-space data to provide a two-dimensional set of spatial domain data; and

applying a 2DFT to said two-dimensional set of spatial domain data to provide another two-dimensional set of k-space spatial frequency domain data for use in subsequent step (b).

4. The method of claim 2, wherein said pre-sequence is a SPAMM (spatial modulation of magnetization) sequence and step (e) uses an inverse cosine function of said ratio to generate B1 map values.

5. The method of claim 2, wherein said pre-sequence is a DANTE (delays alternating with nutations for tailored excitation) sequence.

6. The method of claim 1, wherein step (a) comprises:

acquiring a first set of k-space data using a first tagging pre-sequence of at least two RF excitation pulses having a same nutation angle sense;

acquiring a second set of k-space data using a second tagging sequence of at least two RF excitation pulses having different alternating nutation angle senses; and

arithmetically combining said first and second sets of k-space data on a pixel-by-pixel basis to produce a combined set of k-space data to be used in subsequent steps.

7. The method of claim 6, wherein step (a) further comprises:

applying a 2DFT (two-dimensional Fourier transform) to said combined set of two-dimensional k-space data to provide a further two-dimensional set of spatial domain data; and

applying a 2DFT to said further two-dimensional set of spatial domain data to provide yet another two-dimensional set of k-space spatial frequency domain data for use in subsequent step (b).

8. The method of claim 1, wherein step (b) comprises:

extracting k-space data from a mid-portion of k-space to produce said baseline low-frequency sub-set;

extracting k-space data from a different portion of k-space disposed to one side of said mid-portion to produce said positive higher frequency sub-set; and

extracting k-space data from another different portion of k-space disposed to another side of said mid-portion to produce said negative higher frequency sub-set.

9. The method of claim 8, wherein said extracting steps multiply respectively corresponding shaped window functions by k-space data.

10. The method of claim 8, wherein said extracting steps extract horizontal or vertical strips across k-space data.

11. The method of claim 8, wherein said extracting steps extract rectangular or square areas within k-space data.

12. The method of claim 1, wherein step (c) comprises:

separately applying a 2DFT to each of said sub-sets of frequency-filtered k-space data to convert such into the spatial domain and then using only the magnitude value for each resulting pixel data point to respectively provide said magnitude sets of spatial domain data.

13. The method of claim 1, wherein step (d) comprises:

adding all three magnitude sets of spatial domain data on a pixel-by-pixel basis to provide a map of peak or upper signal magnitude envelope data set; and

subtracting from the baseline magnitude set the sum of the harmonic magnitude sets on a pixel-by-pixel basis to provide a map of trough or lower signal magnitude envelope data set.

14. The method of claim 1, wherein step (e) comprises:

creating a tag depth ratio (TD) data set on a pixel-by-pixel basis where each pixel of the tag depth ratio data set corresponds to TD=min/max where min is the lower magnitude data set value and max is the upper magnitude data set value; and

calculating the arccos of TD on a pixel-by-pixel basis to provide a B1 map of effective tag nutation flip angles.

15. The method of claim 14, wherein step (e) further comprises:

calculating a normalized value for effective/intended tag nutation flip angles on a pixel-by-pixel basis to provide a B1 map of normalized tag nutation flip angles.

16. The method of claim 14, wherein step (e) further comprises:

calculating B1 strength based on said tag nutation flip angles on a pixel-by-pixel basis to provide a B1 map of effective B1 field strengths.

17. The method of claim 16, wherein step (e) further comprises:

calculating a normalized value for effective/intended B1 field strength on a pixel-by-pixel basis to provide a B1 map of normalized B1 field strengths.

18. The method of claim 1, further comprising:

(g) displaying, on a display screen or printer, said B1 map of data values for visual observation by a user of the MRI system.

19. The method of claim 1, further comprising:

(g) using said stored B1 map of data values to provide image compensation to a diagnostic image acquired by said MRI system, thereby providing a corrected diagnostic image of improved quality.

20. The method of claim 1, wherein step (b) processes said k-space data to produce more than three sub-sets of frequency filtered k-space data by including further sub-sets encompassing respectively associated higher harmonic components of the spatial modulation.

21. The method of claim 20, wherein determination of a B1 map value at each pixel includes a data fitting process using the associated higher harmonic components.

22. The method of claim 1, wherein step (a) comprises constraining the time shift between pre-sequence tagging RF pulses so that the shift of the tag lines for off-resonance difference between fat and water species nuclei is substantially a multiple of the time period needed for one cycle of phase to accrue between these species.

23. The method of claim 1, wherein step (a) utilizes plural RF transmit channels and an acquisition pulse sequence is repeated plural times using different ones or combinations of the RF transmit channels for respectively corresponding different acquisition pulse sequences.

24. A computer-readable storage medium containing computer program code which, when executed, effects the method of claim 1.

25. An MRI system comprising:

(a) an MRI scanner configured to acquire at least one set of B1 amplitude tagged MRI data signals in spatial frequency domain k-space from MR nuclei within an imaged volume of the MRI system;

at least one processor configured to;

(b) process said k-space data to produce at least two sub-sets of frequency-filtered k-space data;

(i) a baseline low-frequency sub-set, and (ii) at least one of a positive higher frequency sub-set including a harmonic version of said baseline sub-set and a negative higher frequency sub-set including a harmonic version of said baseline sub-set;

separately transform each of said sub-sets to respectively corresponding spatial domain data sets;

(i) a baseline magnitude set, and (ii) at least one of a positive harmonic magnitude set and a negative harmonic magnitude set;

arithmetically combine said baseline magnitude set with said at least one of the harmonic magnitude sets on a pixel-by-pixel basis to provide upper and lower magnitude data sets;

process said upper and lower magnitude data set values on a pixel-by-pixel basis to generate a B1 map value based on a value related to a ratio of said upper and lower pixel data values at a given pixel location in the spatial domain; and

(c) a data memory configured to store said B1 map of values in a computer accessible and readable memory device for subsequent use by said MRI system.