In recent years, attention has been given to the construction and development of new educational centers, but their spatial distribution across the cities has received less attention. In this study, the Average Nearest Neighbor (ANN) and the optimized hot spot analysis methods have been used to determine the general spatial distribution of the schools. Also, in order to investigate the spatial distribution of the schools based on the substructure variables, which include the school building area, the results of the general and local Moran and Getis Ord analyses have been investigated. A differential Moran index was also used to study the spatial-temporal variations of the schools’ distribution patterns based on the net per capita variable, which is the amount of school building area per student. The results of the Average Nearest Neighbor (ANN) analysis indicated that the general spatial patterns of the primary schools, the first high schools, and the secondary high schools in the years 2011, 2016, 2018, and 2021 are clustered. Applying the optimized hot spot analysis method also identified the southern areas and the suburbs as cold polygons with less-density. Also, the results of the differential Moran analysis showed the positive trend of the net per capita changes for the primary schools and first high schools. However, the result is different for the secondary high schools.