A completely noninvasive three-dimensional (3-D) static magnetic field magnitude spatially localized ³¹P spectroscopy technique has been developed and applied to study the in vivo canine myocardium at 9.4 T. The technique incorporates both Fourier series windows and selective Fourier transform methods utilizing all three orthogonal gradients for 3-D phase encoding. The number of data acquisitions for each phase-encoding step was weighted according to the Fourier coefficients to define cylindrical voxels. Spatially localized ³¹P spectra can be generated for voxels of desired location within the field of view as a postprocessing step. The quality of localization was first demonstrated by using a three-compartment phantom. The technique was then applied to in vivo canine models and yielded ³¹P cardiac spectra with an excellent signal-to-noise ratio. The in vivo validation experiments, using an implanted 2-phosphoenolpyruvate-containing marker, demonstrated that the technique is capable of measuring at least two transmural layers of left ventricular myocardium representing the subepicardium and subendocardium.