Gross ecosystem product (GEP) accounting can be used to reveal the contributions of ecosystems to economic and social development and human well-being. However, most current studies have only focused on the accounting of GEP values within an administrative unit,this type of data cannot reflect the spatial heterogeneity within the statistical unit even though it can reflect the differences between the statistical units. In this paper, we construct a spatialization method using the net primary productivity (NPP) and tourism suitability index to spread the total GEP in administrative areas over geographic grids of certain sizes, and then, we take the Zhongshan District, Guizhou Province, China, as a case study for analysis. The results show that (1) the constructed spatialization method can reflect the spatial heterogeneity of the GEP within administrative units and ecosystem types, objectively representing the characteristics of the spatial distribution of the GEP. (2) The spatial distribution patterns of the value of product provisioning, the value of regulating services, and the total value of the GEP are consistent, i.e., high in the northwestern part of the study area and low in the southeastern part, while the value of cultural services exhibits the opposite spatial pattern. (3) The GEPs of the different ecosystem types differ significantly: cropland > forest land > wetlands > shrubland > grassland > urban land. Spatialization of the GEP can realize overlaying analysis with multiple driving factors, providing technical support for further analysis of the change patterns of the GEP.