An automated method for analysis of in vivo proton magnetic resonance (MR) spectra and reconstruction of metabolite distributions from MR spectroscopic imaging (MRSI) data is described. A parametric spectral model using acquisition specific, a priori information is combined with a wavelet-based, nonparametric characterization of baseline signals. For image reconstruction, the initial fit estimates were additionally modified according to a priori spatial constraints. The automated fitting procedure was applied to four different examples of MRS data obtained at 1.5 T and 4.1 T. For analysis of major metabolites at medium TE values, the method was shown to perform reliably even in the presence of large baseline signals and relatively poor signal-to-noise ratios typical of in vivo proton MRSI. Identification of additional metabolites was also demonstrated for short TE data. Automated formation of metabolite images will greatly facilitate and expand the clinical applications of MR spectroscopic imaging.