Method and procedure for signal estimation and data harmonization for magnetic resonance spectroscopy (MRS)
A method and a system for analysis of raw MRS data, in the form of signal strength versus chemical shift (ppm), from multiple scanners, includes “signal estimation” from each raw data set, followed by cross-scanner “data harmonization” of results. The final resulting MRS signals are consistent from one scanner to another, and are used for analysis by radiologists and other physicians.
1. A magnetic resonance spectroscopy (MRS) system, comprising:
target and reference MRS scanners producing MRS data of subjects; and
an analysis system that performs signal estimation to eliminate outlier spectra and a data harmonization process in which the spectra from control subjects from the target scanner are adjusted and normalized scans to ensure that the amplitudes and frequencies of the peaks correspond to scans from the reference scanner to facilitate comparison of the relative abundance of specific metabolites as measured by the target and reference scanners;
wherein the data harmonization includes dividing signals from the target and reference MRS scanners for the control subjects by respective standard deviations to get two normalized signals for the target and reference scanners, then a ratio of normalized signals of the reference scanner to normalized signals of the target scanner are used to adjust the signals from the target scanner, then the adjusted signals of the target scanner are multiplied by the standard deviation of the signals of the target scanner.
2. A system as claimed in claim 1 , wherein the control subjects are scanned in each of the scanners.
3. A system as claimed in claim 1 , further comprising calculating a mean as a function of parts per million for each of the control subjects.
4. A system as claimed in claim 1 , further comprising repeating the signal estimation until a threshold condition is met.
5. A system as claimed in claim 4 , wherein the threshold condition is a minimum percentage of raw data that must be included to compute a mean.
6. A system as claimed in claim 4 , wherein the threshold condition is a specified signal-to-noise (SNR) ratio.
7. A system as claimed in claim 1 , wherein the multiple scanners are pair-wise harmonized with the reference scanner.
8. A method for harmonizing magnetic resonance spectroscopy (MRS) data, comprising:
producing MRS data of control subjects with a target scanner and a reference scanner;
performing signal estimation to eliminate outlier spectra and data harmonization in which the spectra from the control subjects from the target scanner and the reference scanner are adjusted and normalized to ensure that the amplitudes and frequencies of the peaks correspond to scans from a reference scanner to facilitate comparison of the relative abundance of specific metabolites as measured by the multiple scanners, wherein the data harmonization includes dividing signals from the target and reference MRS scanners for the control subjects by respective standard deviations to get two normalized signals for the target and reference scanners, then a ratio of normalized signals of the reference scanner to normalized signals of the target scanner are used to adjust the signals from the target scanner, then the adjusted signals of the target scanner are multiplied by the standard deviation of the signals of the target scanner.
9. A method as claimed in claim 8 , further comprising using a resulting mapping to harmonize further MRS data.
10. A method as claimed in claim 8 , wherein the control subjects are scanned in each of the scanners.
11. A method as claimed in claim 8 , further comprising calculating a mean as a function of parts per million for each of the control subjects.
12. A method as claimed in claim 8 , further comprising repeating the signal estimation until a threshold condition is met.
13. A method as claimed in claim 12 , wherein the threshold condition is a minimum percentage of raw data that must be included to compute a mean.
14. A method as claimed in claim 12 , wherein the threshold condition is a specified signal-to-noise (SNR) ratio.
15. A method as claimed in claim 8 , wherein the multiple scanners are pair-wise harmonized with the reference scanner.