基于DEM的甘孜县土地利用格局空间自相关分析

陈家良 夏建国 雷苑

摘要 基于甘孜县土地利用、DEM数据对甘孜县土地利用格局进行了空间自相关分析，全县共划分346个网格单元，以每个单元内各地类面积为观测量分别计算了全局和局部空间自相关莫兰指数（Morans I），绘制LISA分布图，结果表明，甘孜县耕地、林地、草地、水域和其他用地均呈现显著的空间聚集性，且为HH型聚集，林地还存在LL聚集的现象，结合DEM提取出的海拔、坡度因子，分析地形因素对空间自相关聚集特征的影响。甘孜县各地类的聚集特征均与相应区域的地形情况密切相关，地形是土地利用格局的重要影响因素。

关键词 土地利用格局，空间自相关，DEM，Morans I指数，地形分析，甘孜县

中图分类号 F301.2文献标识码 A

文章编号 0517-6611（2020）06-0059-05

Abstract We analyzed the spatial autocorrelation of land use patterns in Ganzi County based on land use and DEM data from Ganzi County.346 grid units were divided， and the global and local spatial autocorrelation Moran index was calculated based on the observations of the area of each type within each cell.The LISA distribution map was drawn，and the cultivated land， forest land， grassland， water area and other land in Ganzi County showed significant spatial aggregation， and it was HH clustering， and there was also LL clustering phenomenon in forest lands.Combined with the elevation and slope factors extracted from DEM， the influence of topographic factors on spatial autocorrelation aggregation characteristics was analyzed.The results show that the agglomeration characteristics of various types of Ganzi County are closely related to the topographical conditions of the corresponding regions， and the topography is an important factor influencing the land use pattern.

Key words Land use patterns，Spatial autocorrelation，DEM，Morans I index，Terrain analysis，Ganzi County

在同一地理空间中，一些变量所呈现出的与其他变量的潜在互相依赖和关联性即为空间自相关性，空间自相关能够度量变量的空间分布特点及其对相邻单元的作用情况，反映空间分布的聚集情况和程度[1-5]。在土地利用的研究领域中，土地利用格局的空间自相关为在同一研究区中，各土地利用观测量的空间分布特点和对相邻土地单元的影响，表现形式为空间聚集。相关研究中，国内学者构造用地比率[6]、土地利用程度综合指数[7]等作为空间自相关分析观测值，该研究构造渔网，以网格单元内某地类占网格面积为变量，分析区域土地利用格局的空间特点。

自然条件是人类活动的基础，是土地利用的主要影响因素和限制条件，而地形要素在自然条件中起主要作用，常作为国内外土地利用格局和景观特征研究[8-10]中的重要因子进行分析。为获取甘孜县土地利用格局的空间聚集特点，笔者应用全局和局部空间自相关分析甘孜县土地利用的自相关格局，并利用分区统计地形因子相关值，链接地形与土地利用显著聚集特征的关系，以期为地方政府未来土地规划提供依据，也进一步探索自然条件与土地利用关系，丰富土地利用相关研究。

1 资料与方法

1.1 研究区概况

甘孜县位于四川甘孜藏族自治州（99°08′～100°25′E、31°24′～32°54′N），面积6 739.39 km2，位于横断山脉东北边缘，海拔在3 325～5 688 m，县城海拔3 410 m，是丘原向山原的过渡区，属青藏高原的一部分，高山寒温带气候，县域地形起伏大，地貌丰富，具有较好的典型性和一定的代表性。甘孜县辖1镇、21个乡，境内多民族聚居。甘孜县土地利用类型主要为耕地、林地、草地、建设用地、水域和其他用地等，其中其他用地为以荒地为主的自然保留地，草地为主要地类，面积占全县面积比达79%，其次为林地、耕地，分别占13.0%和3.7%。

1.2 数据来源及处理

甘孜县土地利用现状图来源于甘孜县人民政府，数据类型为栅格数据，比例尺1∶100 000，利用ArcGIS軟件对其进行矢量化处理，得到甘孜县土地利用图（图1），构造渔网，将研究区划分为 346个5 000 m×5 000 m的网格单元，每个网格单元内各类用地的面积为空间自相关分析的观测量，DEM 数字高程来自美国国家航空航天局，分辨率为30 m，数据类型为栅格数据，坐标系为UTM/WGS84，利用ArcGIS软件，采用自然间断点分级法（Jenks）分5段显示甘孜县域海拔情况，提取坡度因子，也使用自然间断点分级法分6段显示（图2）。

1.3 研究方法

1.3.1 空间全局自相关。

全局空间自相关是对数据集的空间自相关的总体估计[11]，经典空间自相关统计包括Morans I、Gearys C、Getiss G和标准偏差椭圆等，该研究的研究对象甘孜县土地利用为连续分布，选取全局Morans I进行格局分析，其计算公式[12-15]如下：

1.3.2 空间局部自相关。

全局空间自相关分析主要描述观测量在整个研究区内的空间特征，而无法表现其聚集或异常的具体分布，该研究进一步分析甘孜县局部空间自相关情况，采用Moran散点图和空间关联局域指标（local indicators of spatial association，LISA） [16-18]来刻画聚集或异常的具体情况，以进一步分析甘孜县土地利用空间分布特征。局部Morans I的计算公式如下：

1.3.3 分区统计。

分区统计由ArcGIS软件[19]空间分析工具实现，可通过叠加矢量图层和栅格图层对空间要素进行描述性统计。该研究将显著空间聚集或异常的网格单元作为分区统计的基本单位，汇总各聚集或异常的单元的海拔、坡度的最小值、最大值和平均值并分析其驱动因素。

2 结果与分析

2.1 甘孜县土地利用全局自相关分析

将划分的每个网格单元内各类用地的面积作为观测变量，通过全局自相关方法，利用ArcGIS软件计算出甘孜县耕地、林地、草地、建设用地和其他用地的Morans I值、Z score和P value（表1），以分析其土地利用的空间分布特征。

由表1可知，在置信度95%的水平上，耕地、林地、草地、水域和其他用地的Morans I值均为正值，由此可拒绝零假设，说明甘孜县土地利用中耕地、林地、草地、水域和其他用地均表现出了显著的空间自相关特征即空间聚集性，在各网格单元呈现正相关性。建设用地Morans I值为负，且P值高达0.944 394，未通过自相关显著性检验，即甘孜县建设用地的空间分布无空间聚集现象，为随机分布，故后续局部空间自相关分析及区域综合统计分析忽略此地类。

2.2 甘孜县土地利用局部空间自相关分析

2.2.1 Moran散点图。

土地利用局部空间自相关主要有4类情况：HH聚集、LL聚集、HL异常和LH异常。利用Geoda软件绘制Moran散点图（图3），其中，横坐标为每一个网格单元的各地类用地面积，纵坐标为每一个网格单元的各地类用地面积空间滞后，即相邻单元的地类用地面积，斜率为两者的相关关系，即Morans I指数。

从图3可以看出，数据集中分布在第一象限即为HH聚集，如耕地、林地、水域和其他用地，表明这4种地类在空间中的聚集特征为大面积的用地单元与大面积的用地单元集聚，数据分布在第二象限为LH异常，即低观测值与高观测值之间的集聚，总体而言，甘孜县土地利用LH异常较少，数据分布在第3象限为LL聚集，如草地，表明甘孜县草地的空间分布也呈现出小面积用地单元之间的集聚，数据分布在第四象限为HL聚集，甘孜县各地类此类聚集也相对较少。

2.2.2 LISA分布图。Moran散点图能够简明地表现研究区各地类聚集或异常的类型，但并不能指出其发生聚集或异常的具体区域，该研究绘制LISA分布图进一步透视和剖析甘孜县土地利用空间分布特征。从计算汇总每个网格单元内各类用地的面积并绘制的LISA分布图 （图4）可以看出，甘孜县耕地、林地、草地、水域和其他用地均有表现出显著的HH聚集，即5种地类各占网格单元的面积之间存在显著正相关，且各地类所占面积较高的单元之间发生空间聚集，相邻单元空间差异度小。草地所占面积较低的网格单元与同样低的发生空间聚集，即为LL聚集，呈现空间聚集格局。Z值<-1.96表示观测值之间存在显著负相关，草地网格单元之间高观测值与低观测值聚集在一起，即HL异常，而水域中一个网格单元中，低观测值倾向于与高观测值聚集在一起，即LH异常，呈现空间异常格局。

LISA分布图不仅能表现空间聚集或异常的情况，还可直观反映出聚集或异常发生的具体区域以及网格单元数目。结合图1和图4可以看出，耕地HH聚集主要发生在县域南部的雅砻江河谷沿岸，且聚集网格单元数目较多，聚集规模较大，河流为农业生产提供灌溉水源，流水沉积带来较为肥沃的土壤，致使此区域耕地大面积集聚，林地在县域东部谷地呈现出显著的HH聚集，聚集网格单元数多，规模大，其聚集区域的山谷中有河流发育，较为湿润，且温度较高，适宜林木生长，草地在县域北部表现出了显著的HH聚集特征，而东部及西部少数区域呈显著LL聚集，北部为海拔较高的山区，温度较低，无法生长乔木和灌木，为草类适宜环境，水域主要在中北部和南部呈现HH聚集，主要为山谷河流易发育地带，聚集区域在县域内分布较广，而非集中在某一处，其他用地主要在西部和南部呈现HH聚集，因其他用地主要为自然保留地，其分布与西南部自然条件、河谷地带人类农业活动和城镇建设都密切相关。

2.3 甘孜县土地利用自相关格局 DEM 分析

地形在土地利用中起主导作用，与土地利用的空间分布特征密切相关，也使得土地利用自相关情况有规律可循，且人类在地理空间中的活动也是以地形为基础的[19]。因此，研究地形因素和土地利用空间聚集或异常特征的相互关系有助于更好地理解区域土地利用格局。通过ArcGIS软件分析甘孜县DEM，提取海拔、坡度數据，如图2。

基于网格单元各地类所占面积，结合土地利用类型LISA分布图（图4）和海拔、坡度因子分析各土地利用类型显著性聚集或异常区域内的地形特点。利用ArcGIS软件的空间区域分析工具进行分区统计，得到各地类海拔和坡度因子的数值变化幅度和平均值，以进一步分析土地利用中呈现空间聚集或异常区域的地形特征，也可依托此分析透视某土地利用类型在未来发生聚集或异常的可能性。计算得出显著性HH、LL聚集区下的地形因子统计值如表2和表3所示。

结合表2、表3及LISA分布图（图4）可以看出，耕地HH聚集主要发生在区域南部，相较其他几种土地利用类型而言，耕地聚集区海拔跨度更小，甘孜县主要农作物为青稞、豌豆、小麦，作物生长过程对光、水、热均有一定的要求，在县域整体海拔较高的基础上，作物分布在相对低海拔区能获得较好的生长发育条件，甘孜县位于横断山脉的东北部，以山地为主且地形起伏较大，县域总体而言较不适宜农业种植业发展，但县域南部为雅砻江河谷，坡度较小且地形起伏较县域其他地区更小，地勢平坦开阔，光照条件好，便于农业生产，故西南部耕地大面积集聚。林地在区域东部呈现出显著的HH聚集，东部主要为谷地，海拔较低，温度较高，且较为湿润，适宜林木生长，坡度条件对林木生长的限制性较小，故甘孜县林木大面积聚集区坡度变化幅度较大，且总体而言林木聚集区平均坡度最大。草地在区域北部表现出了显著的HH聚集特征，较其他植物而言，草类环境适应

性较强，甘孜县草地聚集区海拔跨度较高但坡度和缓，适宜生长，甘孜县主要为山寒温带气候，气温条件致使县域内草地大面积分布，尤其在北部海拔较高地区，东部及西部少数区域呈显著LL聚集，此区域为林地聚集区，故草地面积较小，但由于林木与草类可共生，且该区生长条件较好，故仍呈现小面积的聚集。水域主要在县域中北部和南部呈现HH聚集，结合高程图可知，中北部和南部多山谷，为河流易发育地带，故水域聚集区平均海拔较低且海拔跨度较低，甘孜县域内主要为雅砻江、达曲、泥曲3条河流顺其山脉走向流经县境，为常年性河流，水文特征较为稳定，故其聚集区坡度较低。其他用地在西部和南部呈现HH聚集，由于其他用地主要是自然保留地，是目前还未利用的土地，包括难利用地、撂荒地等，西部和南部为甘孜县城镇建设区域，该区不合理的土地利用会造成土地退化，且当地为生态脆弱区，土地利用难度较大，一旦出现土地退化现象，恢复难度较大，故其他用地大面积聚集。

3 结论

该研究基于甘孜县土地利用数据及 DEM 数据，划分346个网格单元，以每个单元内各地类面积为观测量分别计算了全局和局部空间自相关莫兰指数（Morans Ⅰ），绘制LISA分布图，结合海拔和坡度因子进一步分析了甘孜县各地类呈显著集聚区的地形特征，主要结论如下：

（1）甘孜县域土地利用类型空间聚集性耕地、林地、草地、水域及其他用地4种地类在置信度95%的水平的Morans Ⅰ 指数均为正值，表明甘孜县这4种地类的空间分布均呈现了显著的空间聚集性，建设用地未通过显著性检验，无空间聚集，为随机分布。

（2）甘孜县耕地HH聚集主要发生在县域南部的雅砻江河谷沿岸，林地在县域东部谷地呈现出显著的HH聚集，草地在县域北部表现出了显著的HH聚集特征，而耕地、林地和草地在县域东部及西部少数区域呈显著LL聚集，水域主要在中北部和南部呈现HH聚集，其他用地主要在西部和南部呈现HH聚集。

（3）结合由DEM提取的海拔、坡度2个地形因子和LISA分布图，分析各地类HH和LL聚集的地形特征，地形条件是甘孜县域各土地利用类型聚集或异常的空间特征的重要影响因素。

参考文献

[1] 林琳， 马飞.广州市人口老龄化的空间分布及趋势[J].地理研究， 2007， 26（5）：1043-1054.

[2] 黎夏.珠江三角洲发展走廊1988-1997年土地利用变化特征的空间分析[J].自然资源学报， 2004， 19（3）：307-315.

[3] 赵作权.地理空间分布整体统计研究进展[J].地理科学进展，2009， 28（1）： 1-8.

[4] 田红霞， 郑海霞.太原市中心城区建设用地集约利用空间相关性分析[J].山西农经， 2018（1）：35-37.

[5] 曾豪.土地利用格局及其空间自相关动态变化分析：以四川中江县为例[D].成都： 成都理工大学， 2017.

[6] 谷建立， 张海涛， 陈家赢， 等.基于DEM的县域土地利用空间自相关格局分析[J].农业工程学报， 2012， 28（23）：216-224.

[7] 汪煜， 朱振宇， 刘玥.基于DEM的土地利用空间自相关格局分析[J].安徽农业科学，2014，42（23）：8009-8013.

[8] 赵东娟， 齐伟， 赵胜亭， 等.基于GIS的山区县域土地利用格局优化研究[J].农业工程学报， 2008， 24（2）：101-106.

[9] 贾宁凤， 段建南， 乔志敏.土地利用空间分布与地形因子相关性分析方法[J].经济地理， 2007， 27（2）：310-312.

[10] 孟斌， 王劲峰， 张文忠， 等.基于空间分析方法的中国区域差异研究[J].地理科学， 2005， 25（4）：393-400.

[11] 付金霞， 郑粉莉， 李媛媛.小理河流域土地利用空间自相关格局与影响因素分析[J].农业机械学报， 2017， 48（1）：128-138.

[12] 庞宇凤， 周德全.基于空间自相关的县域土地覆盖空间特征分析[J].贵州师范学院学报， 2015， 31（3）：41-45.

[13] 焦利民， 刘耀林， 刘艳芳.区域城镇基准地价水平的空间自相关格局分析[J].武汉大学学报（信息科学版）， 2009， 34（7）：873-877.

[14] 陈彦光.基于Moran 统计量的空间自相关理论发展和方法改进[J].地理研究， 2009， 28（6）：1449-1463.

[15] ANSELIN L.Local indicators of spatial association—LISA[J].Geographical analysis， 1995， 27（2）：93-115.

[16] ANSELIN L，SYABRI I，KHO Y.GeoDa：An introduction to spatial data analysis[J].Geographical analysis， 2006， 38（1）：5-22.

[17] LI H F，CALDER C A，CRESSIE N.Beyond Morans I：Testing for spatial dependence based on the spatial autoregressive model[J].Geographical analysis， 2007，39（4）：357-375.

[18] 吴静， 何必， 李海涛.ArcGIS 9.3 Desktop地理信息系统应用教程[M].北京：清华大学出版社，2011.

[19] 杨昕， 汤国安， 刘学军， 等.数字地形分析的理论、方法与应用[J].地理学报， 2009， 64（9）：1058-1070.