An automated quantitation of short echo time MRS spectra in an open source software environment: AQSES.

This paper describes a new quantitation method called AQSES for short echo time magnetic resonance spectra. This method is embedded in a software package available online from www.esat.kuleuven.be/sista/members/biomed/new/ with a graphical user interface, under an open source license, which means that the source code is freely available and easy to adapt to specific needs of the user. The quantitation problem is mathematically formulated as a separable nonlinear least-squares fitting problem, which is numerically solved using a modified variable-projection procedure. A macromolecular baseline is incorporated into the fit via nonparametric modelling, efficiently implemented using penalized splines. Unwanted components such as residual water are removed with a maximum-phase FIR filter. Constraints on the phases, dampings and frequencies of the metabolites can be imposed. AQSES has been tested on simulated MR spectra with several types of disturbance and on short echo time in vivo proton MR spectra. Results show that AQSES is robust, easy to use and very flexible.

Quantification of prostate MRSI data by model-based time domain fitting and frequency domain analysis.

This paper compares two spectral processing methods for obtaining quantitative measures from in vivo prostate spectra, evaluates their effectiveness, and discusses the necessary modifications for accurate results. A frequency domain analysis (FDA) method based on peak integration was compared with a time domain fitting (TDF) method, a model-based nonlinear least squares fitting algorithm. The accuracy of both methods at estimating the choline + creatine + polyamines to citrate ratio (CCP:C) was tested using Monte Carlo simulations, empirical phantom MRSI data and in vivo MRSI data. The paper discusses the different approaches employed to achieve the quantification of the overlapping choline, creatine and polyamine resonances. Monte Carlo simulations showed induced biases on the estimated CCP:C ratios. Both methods were successful in identifying tumor tissue, provided that the CCP:C ratio was greater than a given (normal) threshold. Both methods predicted the same voxel condition in 94% of the in vivo voxels (68 out of 72). Both TDF and FDA methods had the ability to identify malignant voxels in an artifact-free case study using the estimated CCP:C ratio. Comparing the ratios estimated by the TDF and the FDA, the methods predicted the same spectrum type in 17 out of 18 voxels of the in vivo case study (94.4%).

Frequency-domain method based on the singular value decomposition for frequency-selective NMR spectroscopy.

In several applications of NMR spectroscopy the user is interested only in the components lying in a small frequency band of the spectrum. A frequency selective analysis deals precisely with this kind of NMR spectroscopy: parameter estimation of only those spectroscopic components that lie in a preselected frequency band of the NMR data spectrum, with as little interference as possible from the out-of-band components and in a computationally efficient way. In this paper we introduce a frequency-domain singular value decomposition (SVD)-based method for frequency selective spectroscopy that is computationally simple, statistically accurate, and which has a firm theoretical basis. To illustrate the good performance of the proposed method we present a number of numerical examples for both simulated and in vitro NMR data.

Automatic frequency alignment and quantitation of single resonances in multiple magnetic resonance spectra via complex principal component analysis.

Several algorithms for automatic frequency alignment and quantitation of single resonances in multiple magnetic resonance (MR) spectra are investigated. First, a careful comparison between the complex principal component analysis (PCA) and the Hankel total least squares-based methods for quantifying the resonances in the spectral sets of magnetic resonance spectroscopy imaging (MRSI) spectra is presented. Afterward, we discuss a method based on complex PCA plus linear regression and a method based on cross correlation of the magnitude spectra for correcting frequency shifts of resonances in sets of MR spectra. Their advantages and limitations are demonstrated on simulated MR data sets as well as on an in vivo MRSI data set of the human brain.

A subspace time-domain algorithm for automated NMR spectral normalization.

Recently, two methods have been proposed for quantitatively comparing NMR spectra of control and treated samples, in order to examine the possible occurring variations in cell metabolism and/or structure in response to numerous physical, chemical, and biological agents. These methods are the maximum superposition normalization algorithm (MaSNAl) and the minimum rank normalization algorithm (MiRaNAl). In this paper a new subspace-based time-domain normalization algorithm, denoted by SuTdNAl (subspace time-domain normalization algorithm), is presented. By the determination of the intersection of the column spaces of two Hankel matrices, the common signal poles and further on the components having proportionally varying amplitudes are detected. The method has the advantage that it is computationally less intensive than the MaSNAl and the MiRaNAl. Furthermore, no approximate estimate of the normalization factor is required. The algorithm was tested by Monte Carlo simulations on a set of simulation signals. It was shown that the SuTdNAl has a statistical performance similar to that of the MiRaNAl, which itself is an improvement over the MaSNAl. Furthermore, two samples of known contents are compared with the MiRaNAl, the SuTdNAl, and an older method using a standard. Finally, the SuTdNAl is tested on a realistic simulation example derived from an in vitro measurement on cells.