Magnetic resonance spectroscopic imaging has been recognized for a long time as a powerful tool for biochemical imaging. However, its practical utility is still rather limited due to poor spatial resolution, low signal-to-noise ratio, and long data acquisition times. In this work, we propose a new technique that enables reconstruction of metabolite maps with high spatial resolution. This technique uses a statistical model to incorporate known anatomical boundaries for edge-preserving noise filtering. This statistical reconstruction scheme makes it possible to use very noisy data, thereby enabling the collection of high-resolution data in a reasonable amount of time. We illustrate the performance of this method with images of the N-acetyl-L-aspartate distribution from an in vivo mouse brain.