Four dimensional chemical shift imaging was used to map the relative peak heights of phosphorus metabolites of the liver and overlying skeletal muscle of a normal subject and two patients. The technique provides 31P spectra localised on a voxel-by-voxel basis and may be valuable in mapping heterogeneous structural and metabolic changes in disease.