基于不同权重的土地利用空间自相关模型对比分析

摘要 ［目的］研究不同空间权重矩阵对土地利用空间自相关模型的影响。［方法］以元谋县土地利用变量为基础，首先分析不同权重矩阵下土地利用类型及其驱动因子的空间自相关性，建立基于queen权重矩阵、rook权重矩阵和距离阈值权重矩阵的耕地时空格局演变的空间自回归模型，选取拟合优度、最大似然对数值、赤池信息准则、施瓦茨信息准则、影响因子数量和模型残差空间自相关性等参数作为模型评价指标，对比分析经典线性回归模型、空间滞后模型和空间误差模型的差异。［结果］2018年元谋县土地利用类型和土地利用驱动因子在不同空间权重矩阵下均呈现出较强的空间正相关关系。通过对比3种权重矩阵的经典线性回归模型、空间滞后模型和空间误差模型，发现在同一空间自相关模型中，基于距离阈值权重矩阵的空间自回归模型的拟合效果较好，基于同一个空间权重矩阵，空间误差模型的拟合效果较好。［结论］空间相关性与空间权重矩阵有关，基于距离阈值矩阵的空间误差模型的拟合度最高，解释能力最强，能更好表征云南高原山区的土地利用格局的空间演变。

关键词 土地利用；空间权重矩阵；空间自相关；空间自回归模型；元谋县

中图分类号 F301.2  文献标识码 A  文章编号 0517-6611（2024）11-0054-06

doi：10.3969/j.issn.0517-6611.2024.11.012

Comparative Analysis of Spatial Autocorrelation Models of Land Use Based on Different Weights—A Case Study of Yuanmou County

HU Rui1，2， LIU Zhi-wei1， ZHANG Xiong-hao3 et al

（1.Southwest Survey and Planning Institute of National Forestry and Grassland Administration， Kunming， Yunnan 650031;2.Faculty of Land Resource Engineering， Kunming University of Science and Technology， Kunming， Yunnan 650093；3.Surveying and Mapping Engineering Institute of Yunnan Province， Kunming， Yunnan 650033）

Abstract ［Objective］To study the impact of different spatial weight matrices on the spatial correlation model of land use. ［Method］Based on the land-use variables of the administrative village of Yuanmou County， this paper first analyzed the spatial autocorrelation of land-use types and their driving factors under different weight matrices， and established an algorithm based on the queen weight matrix， rook weight matrix and distance threshold weight matrix. Spatial autoregressive model for the evolution of the spatio-temporal pattern of cultivated land. Parameters such as goodness of fit， maximum likelihood logarithm， Akaike information criterion， Schwartz information criterion， number of influence factors and spatial autocorrelation of model residuals were selected as model evaluation indicators.The differences between the classic linear regression model， the spatial lag model and the spatial error model were compared and analyzed. ［Result］In 2018， land use types and land use drivers in Yuanmou County showed a strong positive spatial correlation under different spatial weight matrices. By comparing the classic linear regression model， spatial lag model and spatial error model of the three weight matrices， it was found that in the same spatial autocorrelation model， the fitting effect of the spatial autoregressive model based on the distance threshold weight matrix was better; based on the same space，the weight matrix and the spatial error model fit better. ［Conclusion］The spatial correlation was related to the spatial weight matrix. The spatial error model based on the distance threshold matrix had the best fitting degree and the strongest interpretation ability， which could better represent the spatial evolution of the land use pattern in the mountainous plateau of Yunnan.

Key words Land use;Spatial weight matrix;Spatial autocorrelation;Spatial autoregression model;Yuanmou County

基金项目 国家自然科学基金项目（41761081）。

作者简介 胡锐（1995—），男，云南文山人，工程师，硕士，从事地理信息系统及自然保护地规划研究。

*通信作者，教授，博士，博士生导师，从事GIS和国土资源信息化研究。

收稿日期 2023-08-21

从土地利用/覆被变化（land use/cover change，LUCC）计划提出至今，已经在世界各国得到了广泛开展和实施。经过20多年的发展，研究内容从全球气候变化效应扩展到不同尺度的LUCC变化过程、驱动机制、模拟预测、资源以及生态环境效应影响等领域，并逐渐成为地理学、遥感科学、生态学等学科的研究重点。在LUCC研究中，土地利用空间格局的空间自相关性成为众多学者研究的一个热点。土地空间自相关反映了空间中2个相邻空间单元土地利用的空间依赖关系，开展土地利用空间自相关的研究对区域土地利用发展具有重要意义。Hope等［1］在2005年对土地利用格局分析时，首先开展了土地利用空间自相关分析，引入了空间自回归模型，将空间自相关正式引入到土地科学研究中。随后，Cai等［2］深入剖析了地形指数间的空间自相关关系；谢花林等［3］、邱炳文等［4］、刘敏等［5］分别以内蒙古、福建省、贵州山区作为研究区域，将研究区进行不同尺度的划分，建立了多尺度的空间滞后模型；袁磊等［6］兼顾土地利用多尺度特征，构建了面向对象的土地利用多尺度时空数据模型，定量分析了土地利用变化在行政尺度层级之间的耦合关系；付金霞等［7］基于全局Moran’s I、Moran散点图和LISA图来分析小理河流域土地利用的空间自相关格局；宁秀红等［8］认为空间自回归模型能够对土地利用驱动因子进行全局参数估计；陈曦炜等［9］研究认为不同土地利用类型在不同尺度的最优空间自回归模型不同，最优空间自回归模型的选取具有尺度效应。总体而言，目前国内已有的研究不同空间权重矩阵的土地利用格局空间自相关性，基于不同空间权重矩阵建立空间自回归模型的文献较少，对土地利用格局的空间自相关性分析及其建立的空间自回归模型侧重于距离权重矩阵，忽略了不同空间权重矩阵下土地利用格局的空间异质性。笔者以元谋县为例，78个行政村作为研究单元，基于不同空间权重矩阵的空间自相关分析方法，分析兼有半山区、丘陵山区和坝区地形的县域土地利用类型及其驱动因子的空间相关性，建立不同空间权重矩阵的耕地时空格局演变的空间自回归模型，探讨不同权重矩阵的土地利用空间自回归模型的差异，以期为云南高原山区的土地利用政策改革、结构调整以及区域规划提供借鉴。

1 资料与方法

1.1 研究区概况

元谋县位于云南省中北部（101°35′～102°06′E、25°23′～26°06′N），是楚雄彝族自治州下辖的一个县，东与武定县接壤，南与禄丰县毗邻，西与大姚县连接，北与四川省会理县交界。地处金沙江中上游，地势为东南高、西北低，地形分为临河坝区、台地半山区和四周丘陵山区3个区域。属南亚热带干热季风气候。元谋县属于云南高原山区的一个典型代表。根据第二次全国土地调查一级类土地分类标准，全县土地一级类可分为6种，以林地和草地为主，分别占土地总面积的41.86%、33.90%，其次为耕地，占土地总面积的19.50%，建设用地、水域、其他用地分别占土地总面积的2.06%、2.08%和0.60%（图1）。

1.2 数据来源

1.2.1 土地利用驱动因子的选取。

土地利用变化是由自然因素和社会经济因素共同作用形成的自然经济综合体，不仅受气候、地形地貌、土壤、水文地质等自然因子的影响，还受社会制度、科学技术、交通条件、人口密度等人为社会因子影响。在云南高原山地区，地形条件决定了土地利用的难易程度，交通条件决定运输成本，总人口和人口密度反映土地的压力状态，人均GDP反映人们对土地利用的程度。由于研究区为高原山区，地表水缺乏，降水补给成为水源获取的重要方式，气象和水源因素也成为影响土地利用的重要因素。结合研究区现状和数据的可获取性、可定量化和空间化［10］，选取的土地利用驱动因子具体包括：①地形因子，包括高程、坡度、坡向、地形起伏度和地表粗糙度；②水资源因子，包括距最近水系的距离；③气象因子，包括降水量和年平均气温；④交通因子，包括距最近乡镇中心的距离、距铁路的距离、距公路的距离和距农村道路的距离；⑤社会经济因子，包括人均GDP、总人口和人口密度。

1.2.2 数据来源及处理。

该研究主要数据来源：①土地利用基础数据来源于2018年元谋县土地利用年度更新数据库，包括道路、水系和居民点等要素（矢量）。 ②中国科学院资源环境科学数据中心提供的2018年云南省的降水量、年平均气温数据及云南省SRTM 30 m分辨率的DEM数据，以元谋县行政区为掩膜提取元谋县降水量、年平均气温和DEM数据。基于元谋县DEM数据提取高程、地形起伏度、地表粗糙度、坡度和坡向5种地形因子。③研究区距离数据通过ArcGIS 10.6的距离分析功能得到。④元谋县行政村人口和社会经济数据来源于2019年元谋县年鉴，通过ArcGIS 10.6连接赋值于研究区矢量图层。

1.3 研究方法

1.3.1 空间权重矩阵。

为了揭示空间的相关性，需要探索空间对象之间的空间依赖关系。在空间统计和计量中，这种空间依赖性可以用空间权重矩阵表示［11］。全局与局部Moran’s I指数的计算均需要空间权重矩阵。因此，研究不同的空间权重矩阵对空间自相关性的影响具有重要意义。

目前有多种空间权重定义方式，如queen权重矩阵、rook权重矩阵、Delaunay权重矩阵、Dacey矩阵、Cliff-Ord矩阵、K-最近点权重矩阵和距离阈值权重矩阵等。在建立空间自回归模型时需要对称权重，故该研究主要基于一阶queen（图2A）和rook（图2B）邻近空间权重矩阵、距离阈值矩阵3种空间对称权重矩阵进行空间自相关分析。