High-resolution magic-angle-spinning 1H NMR spectroscopy (HR-MAS NMR) is a well-established technique for assessing the biochemical composition of intact tissue samples. In this study, we utilized a method based on HR-MAS NMR spectroscopy with slice localization (SLS) to achieve spatial resolution of metabolites. The obtained 7 slice spectra from each of the model samples (i.e., chicken thigh muscle with skin and murine renal biopsy including medulla (M) and cortex (C)) showed distinct metabolite compositions. Furthermore, we analyzed previously acquired 1H HR-MAS NMR spectra of separated cortex and medulla samples using multivariate statistical methods. Concentrations of glycerophosphocholine (GPC) were found to be significantly higher in the renal medulla compared to the cortex. Using GPC as a biomarker, we identified the tissue slices that were predominantly the cortex or medulla. This study demonstrates that HR-MAS SLS combined with multivariate statistics has the potential for identifying tissue heterogeneity and detailed biochemical characterization of complex tissue samples.