Multicellular spheroids (hereinafter referred to as spheroids) are 3D biological models. The metabolomic profiles inside spheroids provide crucial information reflecting the molecular phenotypes and microenvironment of cells. To study the influence of an anticancer drug on the spatially resolved metabolites, spheroids were cultured using HCT-116 colorectal cancer cells, treated with the anticancer drug Irinotecan under a series of time- and concentration-dependent conditions. The Single-probe mass spectrometry imaging (MSI) technique was utilized to conduct the experiments. The MSI data were analyzed using advanced data analysis methods to efficiently extract metabolomic information. Multivariate curve resolution alternating least square (MCR-ALS) was used to decompose each MS image into different components with grouped species. To improve the efficiency of data analysis, both supervised (Random Forest) and unsupervised (cluster large applications (CLARA)) machine learning (ML) methods were employed to cluster MS images according to their metabolomic features. Our results indicate that anticancer drug significantly affected the abundances of a variety of metabolites in different regions of spheroids. This integrated experiment and data analysis approach can facilitate the studies of metabolites in different types of 3D tumor models and tissues and potentially benefit the drug discovery, therapeutic resistance, and other biological research fields.