Scientifically identifying the territorially ecological restoration zoning is a vital prerequisite for implementing ecological restoration projects and enhancing environmental quality. Based on remote sensing data, we syste-matically assessed supply and demand for ecosystem service and their relationship in Guyuan City, China by using the InVEST model, coordination degree model, and spatial autocorrelation analysis. We carried out territorially ecological restoration zoning by coupling the ecosystem service supply and demand. Furthermore, the corresponding optimization strategies were put forward according to the regional characteristics of natural resources and socio-economic development. The results showed that the areas with high water yield, carbon storage, soil conservation and habitat quality were mainly located in the southern part of Guyuan City. The areas with high population density, economic development level and high ecological demand were mainly located in the central urban area and the location towns of the county government. The spatial mismatch of supply and demand for ecosystem services in Guyuan City was dominant, as indicated by the fact that the areas with high ecological supply having low ecological demand and low ecological supply having high ecological demand. The average coordination index was 0.5, indicating that the relationship between ecological supply and ecological demand was basically coordinated. Based on the diffe-rences between supply and demand of ecosystem services and regional natural geographical pattern, the administrative township units in Guyuan could be classified into six ecological restoration zonings: key ecological restoration area, potential ecological restoration area, ecological economic reconstruction area, characteristic agricultural development area, ecological core protection area, and ecological industry construction area. Different management strategies were proposed to provide scientific support for ecological restoration.