A method for measuring spatially resolved high‐resolution nuclear magnetic resonance spectra for a spatially heterogeneous complex object has been examined by making images of several phantoms. The factors which determine the image quality are discussed from both theoretical and experimental viewpoints. Especially, the effects of field inhomogeneity, originating from magnet inhomogeneity and nonuniform magnetic susceptibility of the object itself, are examined in detail. The indirect spin‐spin coupling effect on imaging results is theoretically discussed. Although mathematical treatments are given for three‐dimensional objects, actual measurements were performed for spatially limited two‐dimensional regions by slicing the three‐dimensional objects using a selective excitation technique. The results show that the present method is feasible for practical chemical shift imaging purposes.