Home | About | Journals | Submit | Contact Us | Français |

**|**HHS Author Manuscripts**|**PMC2868913

Formats

Article sections

Authors

Related links

Magn Reson Imaging. Author manuscript; available in PMC 2010 May 12.

Published in final edited form as:

Published online 2008 May 16. doi: 10.1016/j.mri.2008.01.050

PMCID: PMC2868913

NIHMSID: NIHMS58854

Magnetic Resonance Research, Department of Radiology, University of California San Diego, 200 West Arbor Drive, San Diego CA 92103-8226, Phone: 619 471 0520, Fax: 619 471 0503

Mark Bydder: ude.dscu@reddybm; Gavin Hamilton: ude.dscu@notlimahg; Takeshi Yokoo: ude.dscu@ookoyt; Claude B Sirlin: ude.dscu@nilrisc

The publisher's final edited version of this article is available at Magn Reson Imaging

See other articles in PMC that cite the published article.

A method is described for making a weighted linear combination of the spectra acquired by a phased array coil. Unlike most previous combination methods, no special reference points in the data are chosen to represent the coil weights. Instead all the data points are used, which results in more reliable estimation. The method uses singular value decomposition to identify the coils weights and extract the principal component of variation in the signal. Subsequent processing of the combined signal (e.g. Fourier transform, baseline correction, phasing) may proceed as per a single coil acquisition.

Several ways to combine multiple spectra acquired by phased array coils have been proposed (1–9). The basic approach is to form a weighted linear combination of the spectra or free induction decays (fid) using weights that ensure constructive addition of signals with maximized signal to noise ratio (SNR). The various methods differ in how the weights are obtained. These may involve a short calibration scan or identification of special reference points that have high signal.

The model of the spectroscopy fid signal considered in the present study is comparable to those used in previous studies. However, for simplicity the coils are considered to be perfectly decoupled and subject to additive complex Gaussian noise with variance *σ*^{2} that is independently and identically distributed. These ideal conditions can be approximately satisfied by hardware calibrations or by pre-processing the data using measurements of the noise correlation, although this problem is not addressed in the present study. Thus the fid signal *S _{j}* (

$${S}_{j}(t)={A}_{j}exp\left(i{\phi}_{j}\right)s(t)+{\epsilon}_{j}(t)$$

[1]

The goal of combining the coils is to make a weighted linear sum of the *n* fids that takes maximum advantage of the SNR offered by the phased array, which can be expressed as in Eq 2 for complex weights *w _{j}*.

$$S(t)=\sum {w}_{j}{S}_{j}(t)$$

[2]

As discussed in Ref (7) and (10), the ideal choice should ensure that the signals are correctly aligned and that coils more sensitive to the spectroscopy voxel are weighted more heavily. This suggests using *w _{j}* =

Ref (7) proposed using the first time point of the fids as weights, i.e. *w _{j}* =

Some of the above approaches use a fixed reference point where the signal is assumed to be large. This introduces a risk of poor estimation, since there are instances when any given point in the fid or spectrum may be degraded; e.g. when a symmetric echo or saturation/inversion pulses are employed or cancellation occurs due to J-coupling. The use of a single reference point also means that estimates are vulnerable to individual noise fluctuations in *ε _{j}* (

A number of ways to do this have been proposed, notably in Ref (8) using a so-called multichannel nonparametric singular value decomposition (SVD) algorithm. Two other methods were described in Ref (8). For the case when the signals are uncorrelated but the noise variance is different on each channel, a minimization procedure was proposed for a specific error metric although it is known that this does not give a maximum likelihood estimate (11). For the case when the number of components in the spectrum is known, Ref (8) suggested a different approach however it is unclear what the value of such a method might be since it blurs the distinction between data acquisition and data analysis that might be better suited to dedicated tools (12).

By contrast, the SVD approach can be shown to produce the optimal SNR for the situation relevant to Eq 1; a theoretical proof is given in Theory. However the SVD method as presented in Ref (8) contains a scaling ambiguity, insofar as the sum of squares of the combined fid *S*(*t*) must be unity. The present article describes a way to overcome the scaling problem and compares combined spectra using SVD and the method of Ref (7) of the over a range of SNR.

The fid data contained in *S _{j}* (

$$\mathbf{H}=\overline{\mathbf{H}}+\mathbf{N}$$

[3]

In the general case, there are a number of ways to estimate from **H** in this situation (13–15). In the present case a simplification is possible as it is known that the rows of are scalar multiples of the signal vector *s*(*t*) and so is a rank-1 matrix. Ref (8) proposed making a rank-1 approximation of **H** using SVD, **H** = **USV′**, then the estimate of the signal vector from the principal column of **U**. The weights *w _{j}* associated with this combination method are contained in the principal column of

Identical weights can be obtained from the eigenvector corresponding to the largest eigenvalue of **H′ H**. As shown in Ref (16), this choice follows from a formal maximization of the SNR and so the SVD method of Ref (8) may be said to be optimal. A simplified proof is given below based on the comprehensive treatment in Ref (16). The optimization is posed a search for the *n* -vector **v** that maximizes the SNR metric defined in Eq 4.

$$\frac{\mathit{signal}}{\mathit{noise}}=\frac{\left|\right|\overline{\mathbf{H}}\mathbf{v}\left|\right|}{\left|\right|\mathbf{Nv}\left|\right|}=\frac{{\mathbf{v}}^{\prime}{\overline{\mathbf{H}}}^{\prime}\overline{\mathbf{H}}\mathbf{v}}{{\mathbf{v}}^{\prime}{\mathbf{N}}^{\prime}\mathbf{Nv}}$$

[4]

The assumption made in the Introduction that noise is uncorrelated and identically distributed means **N′ N** = *σ*^{2}**I** and so the denominator is constant. Therefore the optimization amounts to maximizing **v′ ′ v** subject to **v′ v** = 1, which is a Rayleigh quotient with known maximum when **v** is the eigenvector corresponding to the largest eigenvalue of **′ ** (17). In practice **′ ** is not available and **H′ H** must be used. For a statistically large sample (*m* *n*) the cross-terms vanish and **H′ H** → **′ ** + *σ*^{2}**I**; since the second term is constant with respect to **v** the location of the maximum is unchanged. This concludes the proof.

Note that the definition of SVD requires the columns of **U** and **V** to have unit norm, which leads to the scaling ambiguity mentioned in Ref (8). This is troublesome for comparing spectra acquired at different TE (say, to correct for T2 decay) because the signals are scaled into the same dynamic range. The scaling problem can be overcome by constructing

$${\mathbf{H}}_{\mathbf{1}}=\mathbf{U}{\mathbf{S}}_{\mathbf{1}}{\mathbf{V}}^{\prime}$$

[5]

where **S _{1}** has all but the largest singular value set to zero. The matrix

Spectra were acquired on a 3T Signa TwinSpeed scanner (General Electric Healthcare, Milwaukee, WI) using an 8-channel head coil and a spectroscopy head phantom. The phantom contained the main brain metabolites: NAA (2.02 ppm), choline (3.02 and 3.91 ppm), creatine (3.2 ppm) and myo-inositol (3.55 ppm). A 10×10×10 mm^{3} voxel was selected near the centre of the phantom and identical voxel positions and shim values were used for all acquisitions. A water-suppressed stimulated echo acquisition mode (STEAM) protocol was performed: TR 3000 ms, TE 10 ms, TM 13.7 ms, 5 kHz bandwidth, 2048 sampled points and with 2–64 averages (NEX) to vary the SNR. Raw fid signals from all coils were processed offline using MATLAB (The Mathworks, Natick, MA). Reconstruction times were less than one second on a 3GHz personal computer.

Figures 1 and and22 show plots of the spectra produced by the two methods for the 16 NEX and 64 NEX data sets; in both Figures panel A shows the method of Ref (7) and panel B shows the method of the present study. Visually it is difficult to assess the differences in noise level therefore two quantitative approaches were used to provide a more definitive comparison.

The construction of **H _{1}** facilitates comparisons with the method of Ref (7), which requires a specific linear combination of the columns of

A second way to compare the methods is to look at the consistency of the weights *w _{j}* in different data sets. Since there is no reason to expect any difference in the weights other than from noise, the vulnerability to individual noise fluctuations should show up as increased variability. Defining the

Tables II and III show examples of the weights (relative to the first coil) determined by the two methods for 16 and 64 NEX data sets. The amplitudes and angles for the method of Ref (7) exhibit large differences whereas those determined by the SVD method are more similar. These qualitative observations are confirmed by the *Error* given in Table IV. The *Error* is substantially lower for the SVD method indicating weights are more consistently estimated from experiment to experiment.

The method proposed in Ref (8) for spectroscopy coil combination using SVD has been investigated and compared with another previously described method (7). It is experimentally verified that the SNR is higher with the SVD method and that the estimates of the coil weights are more consistent from experiment to experiment.

Although only two of the many published methods for spectroscopy coil combination have been compared (1–9), it is meaningful to distinguish between those that use a small number of data points to estimate the weights and those that use all of the data points. Whereas the methods of Ref (1–7,9) are representative of the former, Ref (8) and the present study are representative of the latter. The reason for making this distinction is because the data points are subject to random noise and so the more points used to estimate the weights the greater the reliability of the estimates. A proof of the SNR optimality of using all the data points can be made, based on arguments described in (16,17).

A limitation of the SVD combination method is the requirement for noise to be uncorrelated and identically distributed for the optimality criterion to hold. There are standard ways of pre-processing the signals so that they conform to this model using noise correlation measurements, as reviewed in Refs (7,16) and elsewhere. Noise correlation measurements require separate calibrations or an assumption that some portion of the acquired data is representative of the noise only. A second limitation, shared with most MR acquisitions, is how to compensate for eddy currents, motion during and other sources of non-Gaussian error during the scan. The choice of weights that yield maximum SNR under a Gaussian noise model is unlikely to produce a maximum likelihood estimate with respect to other types of error, although to make any definitive statement would require the error distribution to be defined and the maximum likelihood estimator for that distribution to be calculated.

The authors thank General Electric Healthcare for research support.

**Publisher's Disclaimer: **This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

1. Hardy CJ, Bottomley PA, Rohling KW, Roemer PB. An NMR phased array for human cardiac 31P spectroscopy. Magn Reson Med. 1992;28:54. [PubMed]

2. Wald LL, Moyher SE, Day MR, Nelson SJ, Vigneron DB. Proton spectroscopic imaging of the human brain using phased array detectors. Magn Reson Med. 1995;34:440. [PubMed]

3. Maril N, Lenkinski RE. An Automated Algorithm for Combining Multivoxel MRS Data Acquired With Phased-Array Coils. J Magn Reson Imaging. 2005;21:317. [PubMed]

4. Wright S, Wald LL. Theory and application of array coils in MR spectroscopy. NMR Biomed. 1997;10:394. [PubMed]

5. Schaffter T, Bornert P, Leussler C, Carlsen IC, Leibfritz D. Fast 1H spectroscopic imaging using a multi-element head-coil array. Magn Reson Med. 1998;40:185. [PubMed]

6. Prock T, Collins DJ, Dzik-Jurasz ASK, Leach MO. An algorithm for the optimum combination of data from arbitrary magnetic resonance phased array probes. Phys Med Biol. 2002;47:N39. [PubMed]

7. Brown MA. Time-Domain Combination of MR Spectroscopy Data Acquired Using Phased-Array Coils. Magn Reson Med. 2004;52:1207. [PubMed]

8. Sandgren N, Stoica P, Frigo FJ, Selen Y. Spectral analysis of multichannel MRS data. J Magn Reson. 2005;175:79. [PubMed]

9. Dong Z, Peterson B. The rapid and automatic combination of proton MRSI data using multi-channel coils without water suppression. Magn Reson Imaging. 2007;25:1148. [PMC free article] [PubMed]

10. Roemer PB, Edelstein WA, Hayes CE, Souza SP, Mueller OM. The NMR Phased Array. Magn Reson Med. 1990;16:192. [PubMed]

11. Stoica P, Li J. On nonexistence of the maximum likelihood estimate in blind multichannel identification. IEEE Signal Process Mag. 2005;22:99.

12. Vanhamme L, van den Boogaart A, Van Huffel S. Improved Method for Accurate and Efficient Quantification of MRS Data with Use of Prior Knowledge. J Magn Res. 1997;129:35. [PubMed]

13. De Moor B. The Singular Value Decomposition and Long and Short Spaces of Noisy Matrices. IEEE Transactions on Signal Processing. 1993;41:2826.

14. Van Huffel S. Enhanced resolution based on minimum variance estimation and exponential data modeling. Signal Processing. 1993;33:333.

15. Bydder M, Du J. Noise Reduction in Multiple Echo Data Sets using Singular Value Decomposition. Magn Reson Imaging. 2006;24:849. [PubMed]

16. Walsh DO, Gmitro AF, Marcellin MW. Adaptive Reconstruction of Phased Array MR Imagery. Magn Reson Med. 2000;43:682. [PubMed]

17. Engeln-Müllges G, Uhlig F. In: Numerical Algorithms with C. Schon M, Uhlig F, translators. Springer–Verlag; Berlin, Heidelberg, New York: 1996. p. 166.

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |