System and method for phase offset and time delay correction in magnetic resonance spectroscopy data

• US 8,598,871 B2
• Filed: 05/26/2010
• Issued: 12/03/2013
• Est. Priority Date: 11/26/2007
• Status: Active Grant

• Current Assignee: Mayo Foundation For Medical Education And Research (Mayo Clinic)
• Original Assignee: Mayo Foundation For Medical Education And Research (Mayo Clinic)
• Inventors: Bahn, Mark M
• Primary Examiner(s): Koval, Melissa
• Assistant Examiner(s): CHAN, EMILY Y
• Application Number:
  US12/788,057
• Publication Number:
  US 20110066025A1
• Time in Patent Office: 1,287 Days
• Field of Search: 324300-322, 600407-465
• US Class Current: 324/307
• CPC Class Codes: G01N 24/08 by using nuclear magnetic resonance G01N24/12 takes precedence
• CPC Class Codes: G01R 33/4625 Processing of acquired signals, e.g. elimination of phase errors, baseline fitting, chemometric analysis

Abstract

A method is provided for producing, with a magnetic resonance (MR) system, a nuclear magnetic resonance spectrum that has been corrected for errors arising from phase offsets and time shifts in the acquired spectroscopic data. A spectroscopic data set includes temporal information indicative of the underlying phase offsets and time shifts. In general, this temporal information is utilized to correct for errors in the acquired data. As a result, a plurality of acquired spectroscopic data sets are more accurately combined by first individually correcting each spectroscopic data set for such errors. Exemplary sources of the phase offsets and time shifts include the physical separation between a volume of interest and the detectors in the MRI system. Using the temporal information, T₂ decay in the produced spectrum is also corrected.

Claims

1. A method for producing a spectrum with a magnetic resonance imaging (MRI) system, the steps comprising:

a) exciting, with the MRI system, a volume of interest;

b) acquiring, with the MRI system, a spectroscopic data set from the volume of interest;

c) transforming the spectroscopic data set to produce a spectrum;

d) determining a correction factor using temporal information from the spectroscopic data set that includes a zero time indicative of a time at which a signal is first emitted from the volume of interest following the excitation in step a) and a time delay between the zero time and a time at which the signal is detected by the MRI system in step b); and

e) applying the correction factor to the spectrum produced in step c) to correct for errors attributable to phase and time shifts in the spectroscopic data set acquired in step b).

2. The method as recited in claim 1 in which the errors attributable to phase and time shifts in the acquired spectroscopic data set result from at least one of a time delay occurring between steps a) and b), and a location of the volume of interest with respect to the MRI system.

3. The method as recited in claim 1 in which step d) further comprises determining a separation correction using information about a physical separation between the volume of interest and the MRI system.

4. The method as recited in claim 3 in which step e) includes multiplying each data point in the spectrum by the correction factor and the separation correction factor.

5. The method as recited in claim 1 in which steps a)-c) are repeated to produce a plurality of spectroscopic data sets and corresponding spectra.

6. The method as recited in claim 5 in which step d) includes averaging the plurality of spectra to produce an average spectrum.

7. The method as recited in claim 6 in which step d) includes producing a phase plot by plotting a phase of each data point in the average spectrum as a function of frequency.

8. The method as recited in claim 7 in which step d) includes:

determining the slope of the phase plot;

estimating a shifted slope; and

calculating the correction factor that shifts the determined slope to the estimated slope.

9. The method as recited in claim 5 in which step d) includes:

selecting a subset of frequencies; and

calculating an initial correction factor by performing a regression on a subset of the average spectrum corresponding to the selected subset of frequencies.

10. The method as recited in claim 9 in which step d) includes calculating a correction factor for each spectrum produced in step c) by performing a regression, using the calculated initial correction factor, on a subset of each spectrum corresponding to the selected subset of frequencies.

11. The method as recited in claim 1 further including:

f) producing, from the acquired spectroscopic data set, a plurality of time-shifted spectroscopic data sets; and

g) transforming the plurality of time-shifted spectroscopic data sets to produce a plurality of time-shifted spectra.

12. The method as recited in claim 11 in which step d) includes:

d)i) calculating a zero time from the spectroscopic data set acquired in step b) and the time-shifted spectroscopic data sets produced in step f); and

d)ii) calculating a plurality of phase offsets from the spectrum produced in step c) and the time-shifted spectra produced in step g).

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

producing a plurality of difference spectra, wherein each difference spectrum is produced by subtracting two time-shifted spectra; and

determining the zero time using the difference spectra.

14. The method as recited in claim 13 in which step d)i) includes:

calculating a slope value for each difference spectra using a regression method; and

producing a phase difference line using the calculated slope values.

15. The method as recited in claim 12 in which step d)ii) includes producing a plurality of circle structures using the spectrum produced in step c) and the time-shifted spectra produced in step g), wherein one circle structure is produced for each spectral frequency represented in the spectrum produced in step c).

16. The method as recited in claim 1 further including:

f) correcting the spectroscopic data set acquired in step b) with the correction factor determined in step d) in order to produce a corrected spectroscopic data set;

g) producing, from the corrected spectroscopic data set, a plurality of corrected time-shifted spectroscopic data sets;

h) transforming the plurality of corrected time-shifted spectroscopic data sets to produce a corresponding plurality of corrected time-shifted spectra; and

i) producing an improved spectrum using the plurality of corrected time-shifted spectra.

17. The method as recited in claim 16 in which step i) includes estimating a complex magnitude value for each spectral frequency in the corrected time-shifted spectra.

18. The method as recited in claim 17 in which each complex magnitude value is estimated in step i) by fitting each value in the plurality of corrected time-shifted spectra for a given spectral frequency to a sinusoidal curve.