摘要：以云南省尋甸縣為例，運(yùn)用空間自相關(guān)分析方法，研究8種土地利用類型在垂直方向上的分布特征。結(jié)果表明，8種土地利用類型在垂直尺度上都有較強(qiáng)的空間自相關(guān)性，在15個(gè)海拔梯度內(nèi)分布不均勻，基本都是隨著海拔的上升，數(shù)量先增加后逐漸減少，在15個(gè)海拔梯度內(nèi)的土地利用結(jié)構(gòu)存在差異。8種土地利用類型在垂直方向上呈現(xiàn)空間相關(guān)性，其中其他用地的空間自相關(guān)性最大，其次分別為園地、林地、工礦倉(cāng)儲(chǔ)用地、耕地、建設(shè)用地和水域及水利設(shè)施用地，草地空間自相關(guān)性最小。

關(guān)鍵詞：烏蒙山區(qū)；土地利用；空間自相關(guān)性；海拔梯度；尋甸縣

中圖分類號(hào)：F321.1" " " " "文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：0439-8114（2023）03-0088-08

DOI：10.14088/j.cnki.issn0439-8114.2023.03.016 開放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：

Vertical distribution characteristics of land use in Wumeng Mountain area：

A case study of Xundian County， Yunnan Province

FU Jing1a，1b， DUAN Qing-song1c， LIU Shi-xin1a，1b， XIANG Dong-lei2， LI Jian-hua1a，1c

（1a.Engineering Research Center of Science and Technology of Land and Resources；1b.College of Resources and Environment；1c.College of Water Resources， Yunnan Agricultural University， Kunming" 605201，China；2.Yunnan Land and Resources Planning and Design Institute， Kunming" 605201，China）

Abstract： Taking Xundian County of Yunnan Province as an example， and using the spatial autocorrelation analysis method， the vertical distribution characteristics of eight land use types were studied. The results showed that the eight land use types had strong spatial autocorrelation on the vertical scale， and the distribution was not uniform in the 15 elevation gradients. The number of land use types increased first and then decreased gradually with the elevation rising， and there were differences in the land use structure in the 15 elevation gradients. The eight land use types showed a spatial correlation in the vertical direction， and the spatial autocorrelation of other land use was the highest， followed by garden land， forest land， mining and storage land， farmland， construction land， water area and water conservancy facility land， and grassland had the lowest spatial autocorrelation.

Key words： Wumeng Mountain area； land use； spatial autocorrelation； elevation gradients； Xundian County

土地資源是承載人類發(fā)展的基石，是人類使用、保護(hù)和改造的對(duì)象［1，2］。在生態(tài)文明建設(shè)和鄉(xiāng)村振興戰(zhàn)略大背景下，研究土地利用分布特征可以為低碳農(nóng)業(yè)、特色農(nóng)業(yè)、區(qū)域經(jīng)濟(jì)發(fā)展、生態(tài)環(huán)境變化等研究提供重要依據(jù)。烏蒙山區(qū)位于云貴高原和四川盆地結(jié)合部，地形以山地高原為主，人地矛盾突出，生態(tài)環(huán)境非常脆弱，脫貧群眾對(duì)政策的依懶性大，內(nèi)生動(dòng)力比較差，返貧風(fēng)險(xiǎn)高，城鄉(xiāng)發(fā)展不平衡、不充分，鞏固拓展脫貧攻堅(jiān)成果任務(wù)繁重［3］。土地利用分布特征研究有利于烏蒙山區(qū)推動(dòng)土地資源有效利用，推動(dòng)特色產(chǎn)業(yè)發(fā)展、區(qū)域經(jīng)濟(jì)水平提高。

國(guó)內(nèi)外現(xiàn)有的研究主要圍繞著土地利用類型變化、土地利用類型及其他因子之間的影響研究［4-9］和建立土地利用/土地覆蓋模型系統(tǒng)的宏觀研究［10-15］，但涉及土地利用在較小尺度不同海拔垂直方向上分布規(guī)律的研究較為鮮見。烏蒙山區(qū)立體特征明顯，本研究以海拔為研究尺度，探討土地利用在垂直方向上的分布規(guī)律。選取云南省尋甸回族彝族自治縣（以下簡(jiǎn)稱尋甸縣）為典型縣域代表，一方面因?yàn)閷さ榭h自然、經(jīng)濟(jì)、社會(huì)情況與烏蒙山區(qū)高度吻合，另一方面是因?yàn)閷さ榭h地貌類型復(fù)雜、山多壩區(qū)面積少的特點(diǎn)突出，土地利用水平較低，土地利用結(jié)構(gòu)嚴(yán)重不合理，因此選取尋甸縣作為研究對(duì)象，對(duì)尋甸縣的土地利用垂直規(guī)律進(jìn)行研究，以期為烏蒙山區(qū)的土地利用提供參考。

1 研究區(qū)概況

研究區(qū)域選擇尋甸縣，昆明市轄縣，位于云南省東北部［16］。全境東西長(zhǎng)82.5 km，南北寬75.4 km，土地總面積3 588 km2。尋甸縣轄8個(gè)鎮(zhèn)，2個(gè)鄉(xiāng)，132個(gè)村民委員會(huì)（社區(qū)居民委員會(huì)），1 219個(gè)自然村。尋甸縣礦產(chǎn)資源豐富，海拔垂直差異大，氣候類型多樣，地形交錯(cuò)復(fù)雜、立體特征明顯，海拔西北高而東南低，呈東南傾斜階梯狀，海拔1 441～3 319 m，各海拔梯度內(nèi)不同土地利用數(shù)量占地面積比例差異明顯。

2 數(shù)據(jù)來(lái)源與研究方法

2.1 數(shù)據(jù)獲取與處理

尋甸縣土地利用現(xiàn)狀矢量數(shù)據(jù)來(lái)源于2015年尋甸縣土地利用變更數(shù)據(jù)庫(kù)（因倘甸、鳳合、聯(lián)合、金源4個(gè)鄉(xiāng)鎮(zhèn)自2011年1月1日起托管于轎子雪山開發(fā)區(qū)，故本研究所處理數(shù)據(jù)均不包含上述4個(gè)鄉(xiāng)鎮(zhèn)，且2015年并未建立塘子街道辦事處、金所街道辦事處），DEM柵格數(shù)據(jù)來(lái)源于地理空間數(shù)據(jù)云的GDEMV2 30 m分辨率數(shù)字高程數(shù)據(jù)。運(yùn)用ArcGIS、ENVI等軟件，對(duì)所獲取的數(shù)據(jù)資料進(jìn)行預(yù)處理。將基礎(chǔ)圖件進(jìn)行配準(zhǔn)校正，統(tǒng)一坐標(biāo)系，進(jìn)行裁剪、合并、空間分析等操作，利用ArcGIS軟件的重分類工具將DEM數(shù)字高程數(shù)據(jù)重分類為15個(gè)海拔梯度單元，進(jìn)而得出新的柵格數(shù)據(jù)。利用柵格轉(zhuǎn)面工具將該柵格數(shù)據(jù)轉(zhuǎn)化為具有15個(gè)海拔梯度單元的矢量數(shù)據(jù)。運(yùn)用GeoDa和SPSS軟件對(duì)預(yù)處理的數(shù)據(jù)進(jìn)行處理與分析［17］。

2.2 研究方法

2.2.1 全局自相關(guān) 全局自相關(guān)（Spatial autocorrelation anlysis）是用于描述某種現(xiàn)象在空間的分布情況，以判斷該現(xiàn)象是否在空間上有集聚現(xiàn)象存在。衡量全局空間自相關(guān)（Moran’s I）的指標(biāo)主要有全局Moran’s I指數(shù)、全局Getis’s G指數(shù)和全局Geary’s C指數(shù)，其中最常用的是全局Moran’s I指數(shù)。

式中，[S2=1nxi-x2]；[x=1ni=1nxi]；n為樣本數(shù)量；xi和xj為樣本觀測(cè)值；[x]為n個(gè)觀測(cè)值的平均值；wij為要素i和j的空間權(quán)重。空間權(quán)重的構(gòu)造一般有Contiguity Weight、Distance Weight和k-Nearest Neighbors 3種方法，本研究選取Distance Weigh方法構(gòu)造空間權(quán)重。

依照式（1）所得到的Moran’s I，介于［-1，1］，Moran’s Igt;0表示呈正相關(guān)，Moran’s Ilt;0則表示呈負(fù)相關(guān)，且Moran’s I越大代表相關(guān)性越高，即空間分布上呈集聚現(xiàn)象，Moran’s I=0表示空間呈隨機(jī)分布。

2.2.2 局部空間自相關(guān) 局部空間自相關(guān)區(qū)別于僅判斷某種現(xiàn)象在空間上是否有集聚現(xiàn)象存在的全局空間自相關(guān)，它不僅可以反映局部小區(qū)域單元與相鄰單元之間某種現(xiàn)象的相關(guān)程度，而且可以進(jìn)一步指出集聚發(fā)生的具體位置。

式中，[S2=1ni=1，j≠1n（xj-x）2]；n為樣本數(shù)量；xi和xj為樣本觀測(cè)值；[x]為n個(gè)觀測(cè)值的平均值；wij為要素i和j的空間權(quán)重。

3 結(jié)果與分析

3.1 尋甸縣海拔梯度下土地利用現(xiàn)狀分析

通過(guò)DEM數(shù)字高程數(shù)據(jù)得出，尋甸縣海拔為" "1 638～3 091 m，根據(jù)其氣候、植被、環(huán)境等現(xiàn)實(shí)狀況［15-17］，以100 m為間距，把尋甸縣海拔劃分為" " " "1 600～1 700 m、1 701～1 800 m、1 801～1 900 m、1 901～2 000 m、2 001～2 100 m、2 101～2 200 m、2 201～" " " "2 300 m、2 301～2 400 m、2 401～2 500 m、2 501～" " " "2 600 m、2 601～2 700 m、2 701～2 800 m、2 801～" " " "2 900 m、2 901～3 000 m、3 001～3 100 m共計(jì)15個(gè)垂直梯度帶（圖1），然后以每個(gè)梯度為研究單元對(duì)尋甸縣海拔垂直分布規(guī)律進(jìn)行研究。

尋甸縣土地?cái)?shù)量在每個(gè)海拔梯度內(nèi)不等，隨著海拔的逐漸上升，其數(shù)量呈先急劇上升后緩慢下降的趨勢(shì)，其中土地面積最多的海拔梯度為2 001～" " 2 100 m，土地面積為55 030.61 hm2，占尋甸縣總土地面積的19.60%，土地面積最少的海拔梯度是" " " "3 001～3 100 m，土地面積為50.08 hm2，占尋甸縣總土地面積的0.02%，兩者之間相差1 097.85倍。在15個(gè)海拔梯度內(nèi)面積占比大于10%的海拔梯度有4個(gè)，且集中在1 901～2 300 m，尋甸縣63.68%的土地都分布在這個(gè)范圍（表1）。

3.2 垂直尺度下尋甸縣土地利用面積分布特征

由表2可知，尋甸縣8種土地利用類型在海拔梯度上的分布峰值基本上都出現(xiàn)在2 001～2 100 m，包括耕地、建設(shè)用地、園地、草地、工礦倉(cāng)儲(chǔ)用地及其他用地，只有水域及水利設(shè)施用地分布峰值出現(xiàn)在" " "1 801～1 900 m海拔梯度內(nèi)，林地分布峰值出現(xiàn)在" " 2 101～2 200 m海拔梯度內(nèi)。

耕地、建設(shè)用地幾乎分布于所有海拔梯度內(nèi)，但在3 001～3 100 m海拔梯度內(nèi)，只存在草地、林地和其他用地。其中建設(shè)用地主要集中于1 801～2 200 m 4個(gè)海拔梯度內(nèi)，相比其他用地類型更加集中，4個(gè)海拔梯度內(nèi)建設(shè)用地面積之和占尋甸縣建設(shè)用地總面積的84.01%。與之相匹配的有耕地、園地、水域及水利設(shè)施用地，耕地4個(gè)海拔梯度內(nèi)面積之和占尋甸縣耕地總面積的62.27%，園地4個(gè)海拔梯度內(nèi)面積之和占尋甸縣園地總面積的89.08%，水域及水利設(shè)施用地4個(gè)海拔梯度內(nèi)面積之和占尋甸縣水域及水利設(shè)施用地總面積的93.44%。

尋甸縣耕地在海拔梯度上的分布與尋甸縣土地在海拔上的分布大致吻合，其面積均呈先逐漸上升后下降的趨勢(shì)。尋甸縣建設(shè)用地在海拔梯度上的分布與尋甸縣土地在海拔上的分布也是大致吻合，其面積也是呈先上升后下降的趨勢(shì)。尋甸縣園地在海拔梯度上的分布與其他用地有著較大的區(qū)別，園地只分布在1 600～2 700 m 11個(gè)海拔梯度內(nèi)，但同時(shí)也遵循著隨海拔的上升其面積呈先逐漸上升后逐漸下降的趨勢(shì)。尋甸縣8種土地利用類型中，只有林地和草地在尋甸縣所有海拔梯度上都有分布，林地也遵循著隨海拔上升其面積呈先上升后下降的趨勢(shì)，在剩余的13個(gè)海拔梯度內(nèi)，林地面積均較大。尋甸縣草地是在所有海拔梯度上都有分布的土地利用類型之一，草地在海拔梯度內(nèi)分布相對(duì)較平均而且其分布呈波浪式變化趨勢(shì)，即隨著海拔的上升，草地面積呈先上升后下降，然后又上升又下降的趨勢(shì)。尋甸縣水域及水利設(shè)施用地集中分布于1 801～2 200 m海拔梯度內(nèi)，隨著海拔上升其面積呈先增加后減少，然后又上升又下降的趨勢(shì)，相比其他地類，尋甸縣水域及水利設(shè)施用地在海拔分布上最為集中。尋甸縣工礦倉(cāng)儲(chǔ)用地只分布在1 701～2 800 m共計(jì)11個(gè)海拔梯度內(nèi)，隨著海拔上升其面積呈先上升后下降的趨勢(shì)。在所有用地類型中，其他用地在垂直方向上的分布規(guī)律與尋甸縣土地在垂直方向上的分布規(guī)律最為相似。

3.3 尋甸縣土地利用類型集聚區(qū)垂直分布規(guī)律

3.3.1 尋甸縣土地利用垂直全局相關(guān)性分析 利用GeoDa軟件對(duì)尋甸縣耕地、建設(shè)用地、園地、林地、草地、水域及水利設(shè)施用地、工礦倉(cāng)儲(chǔ)用地和其他用地8種土地利用類型的全局空間自相關(guān)程度進(jìn)行全局空間自相關(guān)分析。由表3可知，各土地利用類型在垂直分布上的全局 Moran’s I 結(jié)果均為正值，P等于0.001，因此說(shuō)明2015年尋甸縣各土地利用類型在垂直方向上的空間分布不是隨機(jī)的，與在水平方向上的空間分布相同，也表現(xiàn)出非常顯著的空間集聚性。即各土地利用類型在整個(gè)尋甸縣范圍內(nèi)的垂直分布上存在極顯著的空間正相關(guān)關(guān)系。Moran’s I 反映了土地利用類型在垂直方向上的空間分布集聚規(guī)律，數(shù)值的大小反映了其空間密集程度，越大反映其在垂直方向上的分布越密集，越小說(shuō)明其在垂直方向上的分布越分散，其中其他用地的空間自相關(guān)性最大，其次分別為園地、林地、工礦倉(cāng)儲(chǔ)用地、耕地、建設(shè)用地和水域及水利設(shè)施用地，草地空間自相關(guān)性最小。

3.3.2 尋甸縣土地利用垂直局部相關(guān)性分析 為了進(jìn)一步分析在垂直方向上發(fā)生空間相關(guān)的具體海拔梯度以及發(fā)生的空間集聚具體類型，對(duì)尋甸縣土地利用類型在海拔上進(jìn)行局部空間自相關(guān)性分析。其中，HH正相關(guān)單元表示該類型土地在海拔區(qū)域?yàn)榭臻g集聚類型，LL正相關(guān)單元表示該類型土地在海拔區(qū)域?yàn)榭臻g孤立類型，LH負(fù)相關(guān)單元表示該類型土地在海拔區(qū)域空間異常，其余則為不顯著相關(guān)單元。結(jié)果如圖2、表4所示，具體分析如下。

1）耕地局部垂直相關(guān)性分析。尋甸縣耕地不顯著相關(guān)單元類型共有5個(gè)，這5個(gè)海拔區(qū)域內(nèi)土地面積占總土地面積的11.38%，說(shuō)明尋甸縣耕地在垂直方向上的空間自相關(guān)性較強(qiáng)。其中HH單元所在海拔區(qū)域土地面積占總土地面積的78.81%，主要集中在1 801～2 400 m 6個(gè)梯度內(nèi)，為中下海拔梯度，周圍海拔梯度內(nèi)耕地面積也較大，說(shuō)明大部分耕地在海拔梯度上呈現(xiàn)空間集聚；LL單元有3個(gè)，所在海拔區(qū)域土地面積占總土地面積的9.39%，主要集中在海拔較高的梯度內(nèi)，該梯度內(nèi)耕地面積較小，周邊梯度耕地面積也較小，說(shuō)明耕地在海拔梯度上空間孤立相比空間集聚較弱；空間異常LH單元所在海拔區(qū)域土地面積占總土地面積的0.43%，為1 701～1 800 m海拔梯度，代表該梯度的耕地面積較小，其周邊的梯度如1 801～1 900 m則耕地面積較大。

2）建設(shè)用地局部垂直相關(guān)性分析。尋甸縣建設(shè)用地不顯著相關(guān)單元類型共有4個(gè)，這4個(gè)海拔區(qū)域內(nèi)土地面積占總土地面積的22.23%，說(shuō)明尋甸縣建設(shè)用地在垂直方向上的空間自相關(guān)性較強(qiáng)。其中HH單元所在海拔區(qū)域土地面積占總土地面積的43.61%，主要集中在1 801～1 900 m和2 001～2 200 m 3個(gè)梯度內(nèi)，為中下海拔梯度，周圍海拔梯度內(nèi)建設(shè)用地面積也較大，說(shuō)明大部分建設(shè)用地在海拔梯度上呈現(xiàn)空間集聚；LL單元所在海拔區(qū)域土地面積占總土地面積的12.79%，主要集中在2 501～2 800 m和2 901～3 100 m 5個(gè)梯度內(nèi)，說(shuō)明建設(shè)用地在海拔梯度上的空間孤立相比空間集聚較弱；空間異常LH單元所在海拔區(qū)域土地面積占總土地面積的21.36%，為1 701～1 800 m和2 201～2 400 m 3個(gè)梯度，表明這些梯度的建設(shè)用地面積較小，其周邊梯度則建設(shè)用地面積較大。

3）園地局部垂直相關(guān)性分析。尋甸縣園地不顯著相關(guān)單元類型有0個(gè)，說(shuō)明尋甸縣園地在垂直方向上的空間自相關(guān)性較強(qiáng)。其中HH單元所在海拔區(qū)域土地面積占總土地面積的57.87%，主要集中在1 801～2 200 m 4個(gè)梯度內(nèi)，為較低的海拔梯度，園地在這些海拔梯度內(nèi)面積較大，周圍海拔梯度內(nèi)園地面積也較大，說(shuō)明大部分園地在海拔梯度上呈空間集聚；LL單元有8個(gè)，這8個(gè)單元所在海拔區(qū)域土地面積占總土地面積的20.76%，說(shuō)明園地在海拔梯度上的空間孤立相比空間集聚較弱；空間異常LH單元所在海拔區(qū)域土地面積占總土地面積的21.36%，為1 701～1 800 m和2 201～2 400 m 3個(gè)梯度，主要集中于HH海拔梯度周邊，這些梯度的園地面積較小，其周邊梯度如1 801～1 900 m則園地面積較大。

4）林地局部垂直相關(guān)性分析。尋甸縣林地不顯著相關(guān)單元類型有8個(gè)，這8個(gè)海拔區(qū)域土地面積占總土地面積的31.61%，說(shuō)明尋甸縣林地在海拔垂直方向上的空間自相關(guān)性較強(qiáng)。其中HH單元所在海拔區(qū)域土地面積占總土地面積的52.78%，主要集中在1 901～2 000 m和2 101～2 400 m 4個(gè)梯度內(nèi)，說(shuō)明大部分林地在海拔梯度上呈現(xiàn)空間集聚；LL單元所在海拔區(qū)域土地面積占總土地面積的9.17%，主要集中在中等海拔梯度內(nèi)，這些梯度內(nèi)林地面積較小，周邊梯度林地面積也較少，說(shuō)明林地在海拔梯度上空間孤立相比空間聚集較弱；空間異常LH單元所在海拔區(qū)域土地面積占總土地面積的6.44%，為" " "1 801～1 900 m海拔梯度，位于HH海拔梯度周邊，這些梯度的林地面積較小，其周邊梯度如1 901～2 000 m則林地面積較大。

5）草地局部垂直相關(guān)性分析。尋甸縣草地不顯著相關(guān)單元類型有10個(gè)，這10個(gè)海拔區(qū)域土地面積占總土地面積的40.78%，說(shuō)明尋甸縣草地在垂直方向上的空間自相關(guān)性較強(qiáng)。其中HH單元所在海拔區(qū)域土地面積占總土地面積的52.78%，主要集中在1 901～2 000 m和2 101～2 400 m 4個(gè)梯度內(nèi)，為中低海拔梯度，草地在這些海拔梯度內(nèi)面積較大，周圍海拔梯度內(nèi)草地面積也較大，說(shuō)明大部分草地在海拔梯度上呈現(xiàn)空間集聚；空間異常LH單元所在海拔區(qū)域土地面積占總土地面積的6.44%，為1 801～" " "1 900 m海拔梯度，位于HH海拔梯度周邊，這些梯度的草地面積較小，其周邊梯度如1 901～2 000 m則草地面積較大。

6）水域及水利設(shè)施用地局部垂直相關(guān)性分析。尋甸縣水域及水利設(shè)施用地不顯著相關(guān)單元類型有4個(gè)，這4個(gè)海拔區(qū)域土地面積占總土地面積的52.78%，說(shuō)明尋甸縣水域及水利設(shè)施用地在小部分海拔梯度上呈垂直空間自相關(guān)。其中HH單元所在海拔區(qū)域土地面積占總土地面積的26.04%，主要集中在1 801～1 900 m和2 001～2 100 m共計(jì)2個(gè)梯度內(nèi)，為中低海拔梯度，水域及水利設(shè)施用地在這些海拔梯度內(nèi)面積較大，周圍海拔梯度內(nèi)水域及水利設(shè)施用地面積也較大，說(shuō)明少數(shù)水域及水利設(shè)施用地在海拔梯度上呈現(xiàn)空間集聚；LL單元所在海拔區(qū)域土地面積占總土地面積的20.76%，主要集中在最低和中高海拔8個(gè)梯度內(nèi)，這些梯度內(nèi)水域及水利設(shè)施用地面積較小，周邊梯度水域及水利設(shè)施用地面積也較小，說(shuō)明水域及水利設(shè)施用地在海拔梯度上的空間孤立相比空間集聚較弱；空間異常LH單元所在海拔區(qū)域土地面積占總土地面積的0.43%，為" 1 701～1 800 m海拔梯度，位于HH海拔梯度周邊，這些梯度的水域及水利設(shè)施用地面積較小，其周邊梯度如1 801～1 900 m則水域及水利設(shè)施用地面積較大。

7）工礦倉(cāng)儲(chǔ)用地局部垂直相關(guān)性分析。尋甸縣工礦倉(cāng)儲(chǔ)用地不顯著相關(guān)單元類型有0個(gè)，說(shuō)明尋甸縣工礦倉(cāng)儲(chǔ)用地在海拔垂直方向上的空間自相關(guān)性較強(qiáng)。其中HH單元所在海拔區(qū)域土地面積占總土地面積的72.37%，主要集中在1 901～2 400 m 5個(gè)梯度內(nèi)，為中等海拔梯度，工礦倉(cāng)儲(chǔ)用地在這些海拔梯度內(nèi)面積較大，周圍海拔梯度內(nèi)工礦倉(cāng)儲(chǔ)用地面積也較大，說(shuō)明大部分工礦倉(cāng)儲(chǔ)用地在海拔梯度上呈現(xiàn)空間集聚；LL單元所在海拔區(qū)域土地面積占總土地面積的20.76%，主要集中在最低和中高海拔梯度內(nèi)，周邊梯度工礦倉(cāng)儲(chǔ)用地面積也較小，說(shuō)明工礦倉(cāng)儲(chǔ)用地在海拔梯度上的空間孤立相比空間集聚較弱；空間異常LH單元所在海拔區(qū)域土地面積占總土地面積的6.87%，為1 701～1 900 m 2個(gè)梯度，位于HH海拔梯度周邊，這些梯度的工礦倉(cāng)儲(chǔ)用地面積較小，其周邊梯度如1 901～2 000 m則工礦倉(cāng)儲(chǔ)用地面積較大。

8）其他用地局部垂直自相關(guān)性分析。尋甸縣其他用地不顯著相關(guān)單元類型有2個(gè)，這2個(gè)海拔區(qū)域土地面積占總土地面積的6.69%，說(shuō)明尋甸縣其他用地在海拔垂直方向上的空間自相關(guān)性較強(qiáng)。其中HH單元所在海拔區(qū)域土地面積占總土地面積的72.37%，主要集中在1 901～2 400 m共計(jì)5個(gè)梯度內(nèi)，為中等海拔梯度，其他用地在這些海拔梯度內(nèi)面積較大，周圍海拔梯度內(nèi)其他用地面積也較大，說(shuō)明大部分其他用地在海拔梯度上呈現(xiàn)空間集聚；LL單元所在海拔區(qū)域土地面積占總土地面積的14.07%，主要集中在最低和中高海拔梯度內(nèi)，這些梯度內(nèi)其他用地面積較小，周邊梯度其他用地面積也較小，說(shuō)明其他用地在海拔梯度上的空間孤立相比空間聚集較弱；空間異常LH單元只有2個(gè)，這2個(gè)海拔區(qū)域土地面積占總土地面積的6.87%，為1 701～1 900 m 2個(gè)梯度，位于HH海拔梯度周邊，這些梯度的其他用地面積較小，其周邊梯度如1 901～2 000 m則其他用地面積較大。

4 小結(jié)與討論

4.1 小結(jié)

選取尋甸縣作為典型縣域進(jìn)行研究，根據(jù)尋甸縣2015年土地利用現(xiàn)狀，把尋甸縣土地利用分為耕地、建設(shè)用地、園地、林地、草地、水域及水利設(shè)施用地、工礦倉(cāng)儲(chǔ)用地和其他用地8種類型，運(yùn)用ArcGIS 10.4軟件得到尋甸縣10個(gè)鄉(xiāng)鎮(zhèn)土地利用數(shù)據(jù)，借助GeoDa軟件對(duì)8種土地利用類型在空間上的分布以及集聚和孤立情況進(jìn)行了分析，結(jié)果如下。

1）在垂直方向上，8種土地利用類型在15個(gè)海拔梯度內(nèi)分布不均勻，基本都是隨著海拔的上升，面積呈先增加后逐漸減少的趨勢(shì)，在15個(gè)海拔梯度內(nèi)的土地利用結(jié)構(gòu)存在差異，在較低和較高海拔梯度內(nèi)土地利用結(jié)構(gòu)比較簡(jiǎn)單。

2）尋甸縣8種土地利用類型在垂直海拔上均呈空間相關(guān)性，其中其他用地的空間自相關(guān)性最大，其次分別為園地、林地、工礦倉(cāng)儲(chǔ)用地、耕地、建設(shè)用地和水域及水利設(shè)施用地，草地空間自相關(guān)性最小。

3）尋甸縣耕地局部垂直相關(guān)HH單元占78.81%，建設(shè)用地局部垂直相關(guān)HH單元占43.61%，園地局部垂直相關(guān)HH單元占57.87%，林地局部垂直相關(guān)HH單元占52.78%，草地局部垂直相關(guān)HH單元占52.78%，水域及水利設(shè)施用地局部垂直相關(guān)HH單元占26.04%，工礦倉(cāng)儲(chǔ)用地局部垂直相關(guān)HH單元占72.37%，其他用地局部垂直相關(guān)HH單元占72.37%。

4.2 討論

在《云南省鄉(xiāng)村振興戰(zhàn)略規(guī)劃（2018—2022年）》分類推進(jìn)鄉(xiāng)村發(fā)展中，提出按集聚提升類村莊、城郊融合類村莊、特色保護(hù)類村莊、搬遷撤并類村莊和其他類型村莊5種類型分類推進(jìn)鄉(xiāng)村振興，對(duì)實(shí)現(xiàn)國(guó)家鄉(xiāng)村振興戰(zhàn)略的階段性目標(biāo)具有重要意義，有利于各類村莊精準(zhǔn)實(shí)施發(fā)展策略，有序發(fā)展，為實(shí)現(xiàn)鄉(xiāng)村振興戰(zhàn)略找準(zhǔn)有利基點(diǎn)，對(duì)有序推進(jìn)鄉(xiāng)村振興、推動(dòng)農(nóng)業(yè)全面升級(jí)、農(nóng)村全面進(jìn)步、農(nóng)民全面發(fā)展具有重要指導(dǎo)意義。通過(guò)分析尋甸縣8種土地利用類型海拔梯度上的垂直空間相關(guān)性，可以結(jié)合尋甸縣自然、社會(huì)、經(jīng)濟(jì)因素對(duì)尋甸縣的村落進(jìn)行分類，進(jìn)行鄉(xiāng)村空間布局優(yōu)化。下一步研究將尋甸縣水平尺度與本次研究結(jié)論相結(jié)合，再進(jìn)一步對(duì)尋甸縣土地利用格局進(jìn)行多時(shí)段和多因子驅(qū)動(dòng)力分析，以便更好地揭示尋甸縣土地利用類型的動(dòng)態(tài)演化與高程的關(guān)系，為尋甸縣土地利用類型科學(xué)合理布局提供規(guī)劃參考，同時(shí)為尋甸縣環(huán)境保護(hù)與經(jīng)濟(jì)發(fā)展相協(xié)調(diào)提供科學(xué)依據(jù)，從而助力尋甸縣發(fā)展與鄉(xiāng)村振興。

參考文獻(xiàn)：

［1］ 幸瑞燊，周啟剛，李 輝，等.基于地形梯度的三峽庫(kù)區(qū)萬(wàn)州區(qū)土地利用時(shí)空變化分析［J］.水土保持研究，2019，26（2）：297-304.

［2］ 石運(yùn)杰，蔡 葵，吳文春，等.基于地形特征分異的土地利用分布研究——以珠穆朗瑪峰國(guó)家級(jí)自然保護(hù)區(qū)為例［J］.水土保持究，2021，28（2）：394-400.

［3］ 錢麗波.促進(jìn)云南烏蒙山區(qū)脫貧攻堅(jiān)與鄉(xiāng)村振興有效銜接［J］.創(chuàng)造，2021，29（5）：43-49.

［4］ 李曉娜，王 超，萬(wàn)秀云，等.延懷盆地不同土地利用類型土壤風(fēng)蝕物特征［J］.水土保持研究，2022，29（1）：14-20，27.

［5］ 施政樂(lè)，張建軍，申明爽，等.晉西黃土區(qū)不同土地利用類型對(duì)土壤水分的影響［J］.水土保持學(xué)報(bào)，2021，35（6）：190-197.

［6］ 宋靜雯，張學(xué)霞，姜東旸，等.不同尺度土地利用方式對(duì)地表水環(huán)境質(zhì)量的影響及驅(qū)動(dòng)機(jī)制［J］.環(huán)境科學(xué)，2022，43（6）：3016-3026.

［7］ 劉 然，徐國(guó)策，李占斌，等.陜西大理河流域枯水期土地利用結(jié)構(gòu)對(duì)河流水質(zhì)的影響［J］.地球科學(xué)與環(huán)境學(xué)報(bào)，2022，44（2）：327-337.

［8］ 鄧宣凱.武漢市土地利用碳排放的影響因素研究——基于擴(kuò)展的Kaya等式和LMDI分解方法［J］.農(nóng)業(yè)與技術(shù)，2021，41（20）：104-109.

［9］ 劉春暉.基于土地利用變化的相對(duì)貧困地區(qū)現(xiàn)狀研究——以陜西省山陽(yáng)縣為例［J］.山西農(nóng)經(jīng)，2021（20）：50-51.

［10］ 王 思，張路路，林偉彪，等.基于MODIS-歸一化植被指數(shù)的廣東省植被覆蓋與土地利用變化研究［J］.生態(tài)學(xué)報(bào)，2022，42（6）：2149-2163..

［11］ PANAYOTOU T. Empirical tests and policy analysis of environmental degradation at different stages of economic development［M］.Geneva，Switzerland：International labour office，1993.

［12］ 張緒清.生態(tài)旅游：烏蒙山區(qū)生態(tài)修復(fù)與脫貧發(fā)展新思路［J］.安徽農(nóng)業(yè)科學(xué)，2010，38（13）：6981-6983.

［13］ 楊盡利，李旭東，邱 添.貴州烏蒙山區(qū)耕地與人口及經(jīng)濟(jì)空間變化探討［J］.貴州師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2012，30（4）：30-34.

［14］ LINDERT P H. Shifting ground： The changing agricultural soils of China and Indonesia［M］.Cambridge， Massachusetts， USA：The MIT press，2000.

［15］ BURGHARDT W. Soil in urban and industrial environments［J］.Journal of plant nutrition and soil science，1994，157（3）：205-214.

［16］ 張?zhí)熘?，張鳳榮，謝 臻，等.精準(zhǔn)扶貧背景下云南少數(shù)民族山區(qū)農(nóng)村居民點(diǎn)空間格局演變［J］.農(nóng)業(yè)工程學(xué)報(bào)，2019，35（9）：246-254.

［17］ 董振華，張繼權(quán)，阿 榮，等.基于GIS的內(nèi)蒙古道路密度空間分異特征研究［J］.東北師大學(xué)報(bào)（自然科學(xué)版），2017，49（1）：154-161.