Spatial autocorrelation plays an important role in geographical analysis; however, there is still room for improvement of this method. The formula for Moran’s index is complicated, and several basic problems remain to be solved. Therefore, I will reconstruct its mathematical framework using mathematical derivation based on linear algebra and present four simple approaches to calculating Moran’s index. Moran’s scatterplot will be ameliorated, and new test methods will be proposed. The relationship between the global Moran’s index and Geary’s coefficient will be discussed from two different vantage points: spatial population and spatial sample. The sphere of applications for both Moran’s index and Geary’s coefficient will be clarified and defined. One of theoretical findings is that Moran’s index is a characteristic parameter of spatial weight matrices, so the selection of weight functions is very significant for autocorrelation analysis of geographical systems. A case study of 29 Chinese cities in 2000 will be employed to validate the innovatory models and methods. This work is a methodological study, which will simplify the process of autocorrelation analysis. The results of this study will lay the foundation for the scaling analysis of spatial autocorrelation.

A number of spatial statistic measurements such as Moran’sIand Geary’sCcan be used for spatial autocorrelation analysis. Spatial autocorrelation modeling proceeded from the 1-dimension autocorrelation of time series analysis, with time lag replaced by spatial weights so that the autocorrelation functions degenerated to autocorrelation coefficients. This paper develops 2-dimensional spatial autocorrelation functions based on the Moran index using the relative staircase function as a weight function to yield a spatial weight matrix with a displacement parameter. The displacement bears analogy with the time lag in time series analysis. Based on the spatial displacement parameter, two types of spatial autocorrelation functions are constructed for 2-dimensional spatial analysis. Then the partial spatial autocorrelation functions are derived by using the Yule-Walker recursive equation. The spatial autocorrelation functions are generalized to the autocorrelation functions based on Geary’s coefficient and Getis’ index. As an example, the new analytical framework was applied to the spatial autocorrelation modeling of Chinese cities. A conclusion can be reached that it is an effective method to build an autocorrelation function based on the relative step function. The spatial autocorrelation functions can be employed to reveal deep geographical information and perform spatial dynamic analysis, and lay the foundation for the scaling analysis of spatial correlation.

Moran’s index is an important spatial statistical measure used to determine the presence or absence of spatial autocorrelation, thereby determining the selection orientation of spatial statistical methods. However, Moran’s index is chiefly a statistical measurement rather than a mathematical model. This paper is devoted to establishing spatial autocorrelation models by means of linear regression analysis. Using standardized vector as independent variable, and spatial weighted vector as dependent variable, we can obtain a set of normalized linear autocorrelation equations based on quadratic form and vector inner product. The inherent structure of the models’ parameters are revealed by mathematical derivation. The slope of the equation gives Moran’s index, while the intercept indicates the average value of standardized spatial weight variable. The square of the intercept is negatively correlated with the square of Moran’s index, but omitting the intercept does not affect the estimation of the slope value. The datasets of a real urban system are taken as an example to verify the reasoning results. A conclusion can be reached that the inner product equation of spatial autocorrelation based on Moran’s index is effective. The models extend the function of spatial analysis, and help to understand the boundary values of Moran’s index.

In geo-statistics, the Durbin-Watson test is frequently employed to detect the presence of residual serial correlation from least squares regression analyses. However, the Durbin-Watson statistic is only suitable for ordered time or spatial series. If the variables comprise cross-sectional data coming from spatial random sampling, the test will be ineffectual because the value of Durbin-Watson’s statistic depends on the sequence of data points. This paper develops two new statistics for testing serial correlation of residuals from least squares regression based on spatial samples. By analogy with the new form of Moran’s index, an autocorrelation coefficient is defined with a standardized residual vector and a normalized spatial weight matrix. Then by analogy with the Durbin-Watson statistic, two types of new serial correlation indices are constructed. As a case study, the two newly presented statistics are applied to a spatial sample of 29 China’s regions. These results show that the new spatial autocorrelation models can be used to test the serial correlation of residuals from regression analysis. In practice, the new statistics can make up for the deficiencies of the Durbin-Watson test.

Urban form takes on properties similar to random growing fractals and can be described in terms of fractal geometry. However, a model of simple fractals is not effectual enough to characterize both the global and local features of urban patterns. In this paper multifractal measurements are employed to model urban form and analyze urban growth. The capacity dimension D ₀ , information dimension D ₁ , and correlation dimension D ₂ of a city's pattern can be estimated utilizing the box-counting method. If D ₀ > D ₁ > D ₂ significantly, the city can be treated as a system of multifractals, and two sets of fractal parameters, including global and local parameters, can be used to spatially analyze urban growth. In this case study, multifractal geometry was applied to Beijing city, China. The results based on the remote-sensing images taken in 1988, 1992, 1999, 2006, and 2009 show that the urban landscape of Beijing bears multiscaling fractal attributes. The dimension spectrum curves show several abnormal aspects, especially the upper limit of the global dimension breaks through the Euclidean dimension of embedding space and the local dimension fails to converge in a proper way. The geographical features of Beijing's spatiotemporal evolution are discussed, and the conclusions may be instructive for spatial optimization and city planning in the future.

Spatial correlation modeling comprises both spatial autocorrelation and spatial cross-correlation processes. The spatial autocorrelation theory has been well-developed. It is necessary to advance the method of spatial cross-correlation analysis to supplement the autocorrelation analysis. This paper presents a set of models and analytical procedures for spatial cross-correlation analysis. By analogy with Moran’s index newly expressed in a spatial quadratic form, a theoretical framework is derived for geographical cross-correlation modeling. First, two sets of spatial cross-correlation coefficients are defined, including a global spatial cross-correlation coefficient and local spatial cross-correlation coefficients. Second, a pair of scatterplots of spatial cross-correlation is proposed, and the plots can be used to visually reveal the causality behind spatial systems. Based on the global cross-correlation coefficient, Pearson’s correlation coefficient can be decomposed into two parts: direct correlation (partial correlation) and indirect correlation (spatial cross-correlation). As an example, the methodology is applied to the relationships between China’s urbanization and economic development to illustrate how to model spatial cross-correlation phenomena. This study is an introduction to developing the theory of spatial cross-correlation, and future geographical spatial analysis might benefit from these models and indexes.