Dynamic land use change simulation using cellular automata with spatially nonstationary transition rules

The dynamic relationships between land use change and its driving forces vary spatially and can be identified by geographically weighted regression (GWR). We present a novel cellular automata (GWR-CA) model that incorporates GWR-derived spatially varying relationships to simulate land use change. Our GWR-CA model is characterized by spatially nonstationary transition rules that fully address local interactions in land use change. More importantly, each driving factor in our GWR model contains effects that both promote and resist land use change. We applied GWR-CA to simulate rapid land use change in Suzhou City on the Yangtze River Delta from 2000 to 2015. The GWR coefficients were visualized to highlight their spatial patterns and local variation, which are closely associated with their effects on land use change. The transition rules indicate low land conversion potential in the city’s center and outer suburbs, but higher land conversion potential in the inner near suburbs along the belt expressway. Residual statistics show that GWR fits the input data better than logistic regression (LR). Compared with an LR-based CA model, GWR-CA improves overall accuracy by 4.1% and captures 5.5% more urban growth, suggesting that GWR-CA may be superior in modeling land use change. Our results demonstrate that the GWR-CA model is effective in capturing spatially varying land transition rules to produce more realistic results, and is suitable for simulating land use change and urban expansion in rapidly urbanizing regions.     

Spatiotemporal variability of urban growth factors: A global and local perspective on the megacity of Mumbai

The rapid growth of megacities requires special attention among urban planners worldwide, and particularly in Mumbai, India, where growth is very pronounced. To cope with the planning challenges this will bring, developing a retrospective understanding of urban land-use dynamics and the underlying driving-forces behind urban growth is a key prerequisite. This research uses regression-based land-use change models – and in particular non-spatial logistic regression models (LR) and auto-logistic regression models (ALR) – for the Mumbai region over the period 1973–2010, in order to determine the drivers behind spatiotemporal urban expansion. Both global models are complemented by a local, spatial model, the so-called geographically weighted logistic regression (GWLR) model, one that explicitly permits variations in driving-forces across space. The study comes to two main conclusions. First, both global models suggest similar driving-forces behind urban growth over time, revealing that LRs and ALRs result in estimated coefficients with comparable magnitudes. Second, all the local coefficients show distinctive temporal and spatial variations. It is therefore concluded that GWLR aids our understanding of urban growth processes, and so can assist context-related planning and policymaking activities when seeking to secure a sustainable urban future.     

The pyrogeography of sub-Saharan Africa: a study of the spatial non-stationarity of fire–environment relationships using GWR

This study analyses the relationship between fire incidence and some environmental factors, exploring the spatial non-stationarity of the phenomenon in sub-Saharan Africa. Geographically weighted regression (GWR) was used to study the above relationship. Environment covariates comprise land cover, anthropogenic and climatic variables. GWR was compared to ordinary least squares, and the hypothesis that GWR represents no improvement over the global model was tested. Local regression coefficients were mapped, interpreted and related with fire incidence. GWR revealed local patterns in parameter estimates and also reduced the spatial autocorrelation of model residuals. All the covariates were non-stationary and in terms of goodness of fit, the model replicates the data very well (R ² = 87%). Vegetation has the most significant relationship with fire incidence, with climate variables being more important than anthropogenic variables in explaining variability of the response. Some coefficient estimates exhibit locally different signs, which would have gone undetected by a global approach. This study provides an improved understanding of spatial fire–environment relationships and shows that GWR is a valuable complement to global spatial analysis methods. When studying fire regimes, effects of spatial non-stationarity need to be incorporated in vegetation-fire modules to have better estimates of burned areas and to improve continental estimates of biomass burning and atmospheric emissions derived from vegetation fires.     

Examining the influences of air quality in China's cities using multi‐scale geographically weighted regression

This study evaluates the influences of air pollution in China using a recently proposed model—multi‐scale geographically weighted regression (MGWR). First, we review previous research on the determinants of air quality. Then, we explain the MGWR model, together with two global models: ordinary least squares (OLS) and OLS containing a spatial lag variable (OLSL) and a commonly used local model: geographically weighted regression (GWR). To detect and account for any variation of the spatial autocorrelation of air pollution over space, we construct two extra local models which we call GWR with lagged dependent variable (GWRL) and MGWR with lagged dependent variable (MGWRL) by including the lagged form of the dependent variable in the GWR model and the MGWR model, respectively. The performances of these six models are comprehensively examined and the MGWR and MGWRL models outperform the two global models as well as the GWR and GWRL models. MGWRL is the most accurate model in terms of replicating the observed air quality index (AQI) values and removing residual dependency. The superiority of the MGWR framework over the GWR framework is demonstrated—GWR can only produce a single optimized bandwidth, while MGWR provides covariate‐specific optimized bandwidths which indicate the different spatial scales that different processes operate.     

Analysing regional industrialisation in Jiangsu province using geographically weighted regression

 Industry is the most important sector in the Chinese economy. To identify the spatial interaction between the level of regional industrialisation and various factors, this paper takes Jiangsu province of China as a case study. To unravel the existence of spatial nonstationarity, geographically weighted regression (GWR) is employed in this article. Conventional regression analysis can only produce `average' and `global' parameter estimates rather than `local' parameter estimates which vary over space in some spatial systems. Geographically weighted regression (GWR), on the other hand, is a relatively simple, but useful new technique for the analysis of spatial nonstationarity. Using the GWR technique to study regional industrialisation in Jiangsu province, it is found that there is a significant difference between the ordinary linear regression (OLR) and GWR models. The relationships between the level of regional industrialisation and various factors show considerable spatial variability. Received: 4 April 2001 / Accepted: 17 November 2001     

Multicollinearity and correlation among local regression coefficients in geographically weighted regression

Present methodological research on geographically weighted regression (GWR) focuses primarily on extensions of the basic GWR model, while ignoring well-established diagnostics tests commonly used in standard global regression analysis. This paper investigates multicollinearity issues surrounding the local GWR coefficients at a single location and the overall correlation between GWR coefficients associated with two different exogenous variables. Results indicate that the local regression coefficients are potentially collinear even if the underlying exogenous variables in the data generating process are uncorrelated. Based on these findings, applied GWR research should practice caution in substantively interpreting the spatial patterns of local GWR coefficients. An empirical disease-mapping example is used to motivate the GWR multicollinearity problem. Controlled experiments are performed to systematically explore coefficient dependency issues in GWR. These experiments specify global models that use eigenvectors from a spatial link matrix as exogenous variables.This study was supported by grant number 1 R1 CA95982-01, Geographic-Based Research in Cancer Control and Epidermiology, from the National Cancer Institute. The author thank the anonymous reviewers and the editor for their helpful comments.     

Mapping local regression for spatial object-pairs

Local regression methods such as geographically weighted regression (GWR) can provide specific information about individual locations (or places) in spatial analysis that is useful for mapping nonstationary covariate relationships. However, the distance-based weighting schemes used in GWR are only adaptable for spatial objects that are point or area features. In particular, spatial object-pairs pose a challenge for local analysis because they have a linear dimensionality rather than a point dimensionality. This paper proposes using an alternative local regression model – quantile regression (QR) – for investigating the stationarity of regression parameters with respect to these linear features as well as facilitating the visualization of the results. An empirical example of a gravity model analysis of trade patterns within Europe is used to illustrate the utility of the proposed method.     

Do Factors Causing Wildfires Vary in Space? Evidence from Geographically Weighted Regression

This paper describes the results of a geo-statistical analysis carried out at the provincial level in Southern Europe to model wildfire occurrence from socio-economic and demographic indicators together with land cover and agricultural statistics. We applied a classical ordinary least squares (OLS) linear regression together with a geographically weighted regression (GWR) to explain long-term wild-fire occurrence patterns (mean annual density of >1 ha fires). The explanatory power of the OLS model increased from 52% to 78% as a result of the non-constant relationships between fire occurrence and the underlying explanatory variables throughout the Mediterranean Basin. The global model we developed (i.e., OLS regression) was not sufficient to fully describe the underlying causal factors in wildfire occurrence modeling. Indeed, local approaches (i.e., GWR) can complement the global model in overcoming the problem of non-stationarity or missing variables. Our results confirm the importance of agrarian activities, land abandonment, and development processes as underlying factors of fire occurrence. The identification of regions with spatially varying relationships can contribute to the better understanding of the fire problem, especially over large geographic areas, while at the same time recognizing its local character. This can be very important for fire management and policy.     

顾及空间非平稳特征的遥感干旱监测

遥感技术具备实时快速、时空连续、广覆盖尺度等独特优势,在全球气候恶化大背景下,利用遥感干旱监测方法相比于传统地面监测手段,能够提供实时、准确、稳定的旱情信息,辅助科学决策。目前常用遥感旱情监测方法大多依赖全域性数学模型建模,假定了旱情模式的空间平稳特性,因而难以准确反映旱情模式的局部差异特征。本文提出利用地理加权回归模型GWR (Geographically Weighted Regression),考虑旱情模式的空间非平稳特性,综合多种遥感地面旱情监测指数,以实现传统全域旱情监测模型的局部优化。以美国大陆为研究区,监测2002年—2011年共10年的旱情状态。研究表明,GWR模型能够提供空间变化的局部最佳估计模型参数,监测结果更加吻合标准美国旱情监测USDM (U.S Drought Monitor)验证数据,且与地面实测值的最高相关系数R达到0.8552,均方根误差RMSE达到0.972,显著优于其他遥感旱情监测模型。GWR模型具备空间非平稳探测优势,实现了旱情模式的局部精细探测,能够显著提升遥感旱情监测精度,具备较好的应用前景。     

Drivers of land cover and land use changes in St. Louis metropolitan area over the past 40 years characterized by remote sensing and census population data

In this study, we explored the spatial and temporal patterns of land cover and land use (LCLU) and population change dynamics in the St. Louis Metropolitan Statistical Area. The goal of this paper was to quantify the drivers of LCLU using long-term Landsat data from 1972 to 2010. First, we produced LCLU maps by using Landsat images from 1972, 1982, 1990, 2000, and 2010. Next, tract level population data of 1970, 1980, 1990, 2000, and 2010 were converted to 1-km square grid cells. Then, the LCLU maps were integrated with basic grid cell data to represent the proportion of each land cover category within a grid cell area. Finally, the proportional land cover maps and population census data were combined to investigate the relationship between land cover and population change based on grid cells using Pearson's correlation coefficient, ordinary least square (OLS), and local level geographically weighted regression (GWR). Land cover changes in terms of the percentage of area affected and rates of change were compared with population census data with a focus on the analysis of the spatial-temporal dynamics of urban growth patterns. The correlation coefficients of land cover categories and population changes were calculated for two decadal intervals between 1970 and 2010. Our results showed a causal relationship between LCLU changes and population dynamics over the last 40 years. Urban sprawl was positively correlated with population change. However, the relationship was not linear over space and time. Spatial heterogeneity and variations in the relationship demonstrate that urban sprawl was positively correlated with population changes in suburban area and negatively correlated in urban core and inner suburban area of the St. Louis Metropolitan Statistical Area. These results suggest that the imagery reflects processes of urban growth, inner-city decline, population migration, and social spatial inequality. The implications provide guidance for sustainable urban planning and development. We also demonstrate that grid cells allow robust synthesis of remote sensing and socioeconomic data to advance our knowledge of urban growth dynamics from both spatial and temporal scales and its association with population change.     

Using Landsat Imagery and Census Data for Urban Population Density Modeling in Port-au-Prince,Haiti

This research uses the most recent (2003) census data and a Landsat ETM+ image to build a population estimation model for Port-au-Prince, Haiti. The purpose of the study is to establish the linkage of population density with remotely sensed surface reflectance signals of an urban area, and use that to estimate population when census data are not available in a timely fashion. The research begins with deriving subpixel vegetation-impervious surface-soil (VIS) fractions derived from the Landsat ETM+ multispectral bands, and then uses the geographically weighted regression (GWR) model to examine how the variation of population density can be explained by the VIS variables and their derivatives. With comparison to the ordinary least square (OLS) model, the GWR model accounts for spatial non-stationarity in the relationship between population patterns and land characteristics in the study area. The study reveals that three VIS variables are significant in explaining population density: the mean value of houses fraction image, the mean value of vegetation fraction image, and the standard deviation of vegetation fraction image.     

Exploring the Impact of Non-normality on Spatial Non-stationarity in Geographically Weighted Regression Analyses: Tobacco Outlet Density in New Jersey

The principal rationale for applying geographically weighted regression (GWR) techniques is to investigate the potential spatial non-stationarity of the relationship between the dependent and independent variables—i.e., that the same stimulus would provoke different responses in different locations. The calibration of GWR employs a geographically weighted local least squares regression approach. To obtain meaningful inference, it assumes that the regression residual follows a normal or asymptotically normal distribution. In many classical econometric analyses, the assumption of normality is often readily relaxed, although it has been observed that such relaxation might lead to unreliable inference of the estimated coefficients' statistical significance. No studies, however, have examined the behavior of residual non-normality and its consequences for the modeled relationships in GWR. This study attempts to address this issue for the first time by examining a set of tobacco-outlet-density and demographic variables (i.e., percent African American residents, percent Hispanic residents, and median household income) at the census tract level in New Jersey in a GWR analysis. The regression residual using the raw data is apparently non-normal. When GWR is estimated using the raw data, we find that there is no significant spatial variation of the coefficients between tobacco outlet density and percentage of African American and Hispanics. After transforming the dependent variable and making the residual asymptotically normal, all coefficients exhibit significant variation across space. This finding suggests that relaxation of the normality assumption could potentially conceal the spatial non-stationarity of the modeled relationships in GWR. The empirical evidence of the current study implies that researchers should verify the normality assumption prior to applying GWR techniques in analyses of spatial non-stationarity.     

Mapping the Results of Geographically Weighted Regression

Abstract

Geographically weighted regression (GWR) is a local spatial statistical technique for exploring spatial nonstationarity. Previous approaches to mapping the results of GWR have primarily employed an equal step classification and sequential no-hue colour scheme for choropleth mapping of parameter estimates. This cartographic approach may hinder the exploration of spatial nonstationarity by inadequately illustrating the spatial distribution of the sign, magnitude, and significance of the influence of each explanatory variable on the dependent variable. Approaches for improving mapping of the results of GWR are illustrated using a case study analysis of population density–median home value relationships in Philadelphia, Pennsylvania, USA. These approaches employ data classification schemes informed by the (nonspatial) data distribution, diverging colour schemes, and bivariate choropleth mapping.     

利用搜索引擎数据模拟疾病空间分布

互联网记录了人们的日常生活,对带有位置信息的搜索引擎数据进行分析和挖掘可以获得隐藏于其中的地理信息。本文通过分析中国各省流感月度发病数与相关关键词百度搜索指数之间的相关性,选取相关性较高关键词的百度指数作为解释变量,发病数作为因变量,在采用主成分分析法消除变量共线性后,分别使用普通最小二乘回归（OLS）、地理加权回归（GWR）及时空地理加权回归（GTWR）构建流感发病数的空间分布模型。模型的拟合度能够从OLS的0.737、GWR的0.915提高到GTWR的0.959,赤池信息准则（AIC）也表明,GTWR模型明显优于OLS与GWR模型。验证结果显示,GTWR模型能准确识别流感高发地区,将该方法与搜索引擎数据结合能较好地模拟流感空间分布,为空间流行病学的研究提供预测模型和统计解释。     

Geographically weighted regression‐based determinants of malaria incidences in northern China

Geographically weighted regression (GWR) is an important local method to explore spatial non‐stationarity in data relationships. It has been repeatedly used to examine spatially varying relationships between epidemic diseases and predictors. Malaria, a serious parasitic disease around the world, shows spatial clustering in areas at risk. In this article, we used GWR to explore the local determinants of malaria incidences over a 7‐year period in northern China, a typical mid‐latitude, high‐risk malaria area. Normalized difference vegetation index (NDVI), land surface temperature (LST), temperature difference, elevation, water density index (WDI) and gross domestic product (GDP) were selected as predictors. Results showed that both positively and negatively local effects on malaria incidences appeared for all predictors except for WDI and GDP. The GWR model calibrations successfully depicted spatial variations in the effect sizes and levels of parameters, and also showed substantially improvements in terms of goodness of fits in contrast to the corresponding non‐spatial ordinary least squares (OLS) model fits. For example, the diagnostic information of the OLS fit for the 7‐year average case is R² = 0.243 and AICc = 837.99, while significant improvement has been made by the GWR calibration with R² = 0.800 and AICc = 618.54.     

湖北省县域粮食单产影响因素的空间异质性研究

在城镇化进程快速推进耕地保护形势严峻的背景下,粮食单产的区域差异研究对地区粮食安全具有重要意义。本文以湖北省粮食单产数据为基础,采用探索性数据分析方法和地理加权回归模型揭示省内县域粮食单产的空间关系和影响因素的空间异质性。结果表明:湖北省县域粮食单产具有显著的空间自相关特征,整体水平稳中有增。农村劳动力、化肥施用量、农村机械总动力和有效灌溉面积比对粮食单产具有正向促进作用和一定的空间分异规律,对农村用电量呈现出先正后负的影响,各因素的空间异质性显著。结合县域现状和因素的区域特质采取对应的有效措施应对粮食安全问题具有深远的现实意义。     

Using Contextualized Geographically Weighted Regression to Model the Spatial Heterogeneity of Land Prices in Beijing,China

Geographically Weighted Regression (GWR) is a method of spatial statistical analysis allowing the modeled relationship between a response variable and a set of covariates to vary geographically across a study region. Its use of geographical weighting arises from the expectation that observations close together by distance are likely to share similar characteristics. In practice, however, two points can be geographically close but socially distant because the contexts (or neighborhoods) within which they are situated are not alike. Drawing on a previous study of geographically and temporally weighted regression, in this article we develop what we describe as contextualized Geographically Weighted Regression (CGWR), applying it to the field of hedonic house price modeling to examine spatial heterogeneity in the land parcel prices of Beijing, China. Contextual variables are incorporated into the analysis by adjusting the geographical weights matrix to measure proximity not only by distance but also with respect to an attribute space defined by measures of each observation's neighborhood. Comparing CGWR with GWR suggests that adding the contextual information improves the model fit.     

基于IGGⅢ的地理加权回归模型研究

针对离群值存在时地理加权回归模型拟合效果较差的问题,本文提出了基于IGGⅢ的地理加权回归方法。核心是采用IGGⅢ方案中的权函数计算权重矩阵,将权因子用于地理加权回归参数估计模型。利用模拟数据和真实数据与GWR、ACV-GWR进行对比试验,以MSE、MAE和R²作为指标对结果进行评价。模拟试验结果显示,IGGⅢ-GWR比GWR性能分别提升了51.14%、23.77%、28.4%,比ACV-GWR分别提升了49.96%、22.57%、27.1%;真实试验结果显示,IGGⅢ-GWR比GWR性能分别提升了12.65%、7.44%、0.37%,比ACV-GWR分别提升了11.85%、6.96%、0.34%。试验结果表明,基于IGGⅢ的地理加权回归可提高模型的抗差能力,拟合效果更好。     

南昌市二手房价格影响因子的空间异质性

以安居客网站爬取的2018年10月894个南昌市住宅小区二手房价格为研究对象,利用地理加权回归模型探讨了建筑特征、邻里特征、区位特征等方面各影响因子对住宅价格的作用差异.研究结果表明:1)地理加权回归(GWR)模型的拟合结果优于OLS模型,将回归系数结果空间可视化发现南昌市二手房价格影响因子具有空间异质性.2)不同因子...     

融合地理探测器和地理加权回归的太原市建设用地变化因素研究

准确地识别城市化进程中建设用地的变化情况及其背后的驱动力,对城市后续的可持续发展具有重要意义。本文首先以2000—2020年遥感影像为基础,对太原市建设用地空间分布变化进行研究,然后结合地理探测器模型和地理加权回归模型,对研究区建设用地的空间分布影响驱动力因子进行研究,得到以下结论:除政策因素外,现有的城市建设用地空间分布变化还受到高程、交通、GDP、人口等因素的显著作用。太原市城市建设用地变化的布局不单是GDP变化、人口变化、海拔高度、公路网密度4个显著性因子均匀、独立、直接作用的结果,而是具有空间异质性的各因子两两交互作用后增效的产物。本文成果有望为城市建设用地驱动力研究提供一种新思路。