IRIS-HSVD algorithm for automatic quantitation of in vivo 31P MRS.

The rapid development of (31)P magnetic resonance spectroscopy (MRS) has enhanced non-invasive measurement of brain metabolites, which is important for biomedical research. The accuracy and efficiency of data post processing and quantification is paramount for MRS applications. One of the difficulties with in vivo(31)P MRS data quantification is the separation of broad line-width resonances from chemical compounds' resonances under a low signal-to-noise ratio condition. Furthermore, the chemical shift of some compounds caused by pH and Mg(2+) concentration can be troublesome. This work aims to develop an automatic algorithm using a state-space based quantification approach to solve the above mentioned problems. To achieve this aim, we utilized an HSVD based adaptive optimizing prior knowledge algorithm, which uses so called "interference" signals to optimize prior knowledge iteratively for parameter optimization. We termed this algorithm IRIS-HSVD, which stands for Iterative Reduction of Interference Signal HSVD. The Monte Carlo evaluations of the algorithm were conducted with simulated data using in vivo parameters commonly obtained from a 4T scanner. The performance of this algorithm using simulated data was compared to those of other automatic methods including HSVD and HTLS-PK. Examples of in vivo(31)P data obtained from brains of healthy subjects on a 4T MRI scanner were also presented, which demonstrated the superiority of the new method. The results were compared with those using AMARES.

[1]  Pietro Luigi Indovina,et al.  A new time-domain frequency-selective quantification algorithm. , 2002, Journal of magnetic resonance.

[2]  S Van Huffel,et al.  Subspace-based MRS data quantitation of multiplets using prior knowledge. , 2004, Journal of magnetic resonance.

[3]  K. Arun,et al.  State-space and singular-value decomposition-based approximation methods for the harmonic retrieval problem , 1983 .

[4]  J. Ruíz-Cabello,et al.  A history of biological applications of NMR spectroscopy , 1995 .

[5]  H. Barkhuijsen,et al.  Aspects of the computational efficiency of LPSVD , 1985 .

[6]  Philipp Birken,et al.  Numerical Linear Algebra , 2011, Encyclopedia of Parallel Computing.

[7]  S Mierisová,et al.  MR spectroscopy quantitation: a review of frequency domain methods , 2001, NMR in biomedicine.

[8]  R. de Beer,et al.  Java-based graphical user interface for MRUI, a software package for quantitation of in vivo/medical magnetic resonance spectroscopy signals , 2001, Comput. Biol. Medicine.

[9]  D. Graveron-Demilly,et al.  Time-domain quantitation of 1H short echo-time signals: background accommodation , 2004, Magnetic Resonance Materials in Physics, Biology and Medicine.

[10]  J. Pettegrew,et al.  Postprocessing method to segregate and quantify the broad components underlying the phosphodiester spectral region of in vivo 31P brain spectra , 2001, Magnetic resonance in medicine.

[12]  D. van Ormondt,et al.  Retrieval of frequencies, amplitudes, damping factors, and phases from time-domain signals using a linear least-squares procedure , 1985 .

[13]  I. Smith,et al.  Magnetic resonance spectroscopy in medicine: clinical impact , 2002 .

[14]  Sabine Van Huffel,et al.  Using prior knowledge in SVD-based parameter estimation for magnetic resonance spectroscopy-the ATP example , 2004, IEEE Transactions on Biomedical Engineering.

[15]  W Gough The graphical analysis of a Lorentzian function and a differentiated Lorentzian function , 1968 .

[16]  Gene H. Golub,et al.  Matrix computations , 1983 .

[17]  Vanhamme,et al.  Improved method for accurate and efficient quantification of MRS data with use of prior knowledge , 1997, Journal of magnetic resonance.

[18]  Jullie W Pan,et al.  Quantitative 31P spectroscopic imaging of human brain at 4 Tesla: Assessment of gray and white matter differences of phosphocreatine and ATP , 2001, Magnetic resonance in medicine.

[19]  Sabine Van Huffel,et al.  Improved methods for exponential parameter estimation in the presence of known poles and noise , 1997, IEEE Trans. Signal Process..

[20]  L. Trefethen,et al.  Numerical linear algebra , 1997 .

[21]  A Heerschap,et al.  Common processing of in vivo MR spectra , 2001, NMR in biomedicine.

[22]  D. van Ormondt,et al.  SVD-based quantification of magnetic resonance signals , 1992 .

[23]  D van Ormondt,et al.  Cramér–Rao bounds: an evaluation tool for quantitation , 2001, NMR in biomedicine.

[24]  D. van Ormondt,et al.  Improved algorithm for noniterative time-domain model fitting to exponentially damped magnetic resonance signals , 1987 .

[25]  Z. S. Liu,et al.  Spectral fitting of NMR spectra using an alternating optimization method with a priori knowledge. , 1999, Journal of magnetic resonance.

[26]  D. van Ormondt,et al.  Analysis of NMR Data Using Time Domain Fitting Procedures , 1992 .

[27]  Sabine Van Huffel,et al.  Parameter Estimation with Prior Knowledge of Known Signal Poles for the Quantification of NMR Spectroscopy Data in the Time Domain , 1996 .

[28]  P ? ? ? ? ? ? ? % ? ? ? ? , 1991 .

[29]  J Higinbotham,et al.  Use of voigt lineshape for quantification of in vivo 1H spectra , 1997, Magnetic resonance in medicine.

[30]  Gene H. Golub,et al.  Singular value decomposition and least squares solutions , 1970, Milestones in Matrix Computation.

[31]  van Ormondt D,et al.  Cramer-Rao bound expressions for parametric estimation of overlapping peaks: influence of prior knowledge , 2000, Journal of magnetic resonance (San Diego, Calif. 1997 : Print).

[32]  S. Vanhuffel,et al.  Algorithm for time-domain NMR data fitting based on total least squares , 1994 .

[33]  S. Huffel,et al.  MR spectroscopy quantitation: a review of time‐domain methods , 2001, NMR in biomedicine.

[34]  K. Uğurbil,et al.  High contrast and fast three‐dimensional magnetic resonance imaging at high fields , 1995, Magnetic resonance in medicine.

[35]  D. Graveron-Demilly,et al.  Java-based graphical user interface for the MRUI quantitation package , 2001, Magnetic Resonance Materials in Physics, Biology and Medicine.

[36]  B. De Moor,et al.  The fit of a sum of exponentials to noisy data , 1987 .

[37]  S Van Huffel,et al.  Improved Lanczos algorithms for blackbox MRS data quantitation. , 2002, Journal of magnetic resonance.

[38]  Petre Stoica,et al.  Frequency-domain method based on the singular value decomposition for frequency-selective NMR spectroscopy. , 2003, Journal of magnetic resonance.