羅丹 王濤 張庭瑜

摘要 以山西省垣曲縣為研究區(qū)，基于研究區(qū)土地利用現(xiàn)狀和DEM數(shù)據(jù)，利用ArcGIS 10.3對(duì)研究區(qū)高程、坡度、坡向等地形因子進(jìn)行提取，將地形特征與土地利用現(xiàn)狀疊置分析，得到該地區(qū)不同高程、坡度、坡向上的土地類型分布情況。結(jié)果表明，土地利用的空間分布與地形因子具有顯著相關(guān)性。在低海拔、小坡度的無坡向和陽坡地區(qū)，土地利用類型多受到人類活動(dòng)的影響，耕地、園地、居民點(diǎn)及工礦用地多分布于這一地區(qū);高海拔、大坡度的地形對(duì)林地、草地的變化影響較大;水域多分布于低海拔、地勢(shì)平緩地區(qū)，受坡向影響不大。

關(guān)鍵詞 土地利用分布;地形因子;ArcGIS;DEM分布指數(shù);垣曲縣

中圖分類號(hào) F 301.2? 文獻(xiàn)標(biāo)識(shí)碼 A

文章編號(hào) 0517-6611（2021）20-0076-04

doi：10.3969/j.issn.0517-6611.2021.20.021

開放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：

Research on the Relationship between Land Use Distribution and Terrain Factors Based on DEM—A Case Study of Yuanqu County，Shanxi Province

LUO Dan WANG Tao ZHANG Ting-yu 3，4，5 （1.Land Surveying，Planning and Design Institute of Shaanxi Provincial Land Engineering Construction Group Co.，Ltd.，Xian，Shaanxi 710075;2.Institute of Land Engineering and Technology，Shaanxi Provincial Land Engineering Construction Group Co.，Ltd.，Xian，Shaanxi 710075;3.Shaanxi Provincial Land Engineering Construction Group Co.，Ltd.，Xian，Shaanxi 710075;4.Key Laboratory of Degraded and Unused Land Consolidation Engineering，Ministry of Natural Resources，Xian，Shaanxi 710075;5.Shaanxi Provincial Land Consolidation Engineering Technology Research Center，Xian，Shaanxi 710075）

Abstract Taking Yuanqu County of Shanxi Province as the research area，based on the current land use status and DEM data of the research area，ArcGIS 10.3 is used to extract terrain factors such as elevation，slope，and aspect of the research area，and the topographic features and current land use status are superimposed and analyzed to obtain different elevations，slopes，and slope upwards in the area.The results showed that the spatial distribution of land use had a significant correlation with topographic factors.In low-altitude，low-slope，non-slope and sunny slope areas，land use types were mostly affected by human activities.Cultivated land，garden plots，residential areas，and industrial and mining land were mostly distributed in this area.High-altitude，large-slope terrain had a greater impact on the changes of woodland and grassland;waters were mostly distributed in low-altitude，gentle terrain，and were not affected by the slope aspect.

Key words Land use distribution;Topographic factors;ArcGIS;DEM distribution index;Yuanqu County

基金項(xiàng)目 長(zhǎng)安大學(xué)中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金項(xiàng)目（300102351502）;陜西省土地工程建設(shè)集團(tuán)有限責(zé)任公司科研項(xiàng)目（DJNY2021-10）。

作者簡(jiǎn)介 羅丹（1991—），女，山西運(yùn)城人，工程師，碩士，從事遙感與GIS應(yīng)用研究。*通信作者，高級(jí)工程師，從事土地測(cè)繪應(yīng)用研究。

收稿日期 2021-05-21;修回日期 2021-06-27

土地是人類賴以生存和發(fā)展的最基本的自然資源[1-2]。土地利用是自然環(huán)境和人類活動(dòng)相互作用的一個(gè)綜合復(fù)雜體[3-4]，是人類主要生產(chǎn)活動(dòng)方式的反映。地形因子作為自然環(huán)境中的重要因素，是土地利用格局空間分布的重要影響因子，深刻影響著人類利用土地的方式和程序。通過對(duì)土地利用類型與地形因子關(guān)系的研究，能夠更加清楚地了解地形因子對(duì)土地利用類型分布的影響規(guī)律[5]，從而更好地掌握土地利用的現(xiàn)狀和未來趨勢(shì)[6]。

黃土高原地形復(fù)雜多樣，地形因子與土地利用空間分布有著密切的關(guān)系。筆者以山西省垣曲縣土地利用數(shù)據(jù)為基礎(chǔ)數(shù)據(jù)，應(yīng)用DEM和ArcGIS軟件，通過空間疊置分析，準(zhǔn)確且快速得到土地利用空間格局與地形因子的關(guān)系[7];探討土地利用空間格局的分布特征，研究土地利用的空間分布規(guī)律及其與地形因子之間的相互關(guān)系，優(yōu)化不同地形上土地利用類型的合理布局，以期為區(qū)域土地可持續(xù)利用、土地管理決策以及生態(tài)環(huán)境建設(shè)提供理論支持[8-9]。

1 資料與方法

1.1 研究區(qū)概況

垣曲縣位于山西省運(yùn)城市東端，地處111°31′～112°10′E、35°00′～35°39′N。東北與陽城、沁水兩縣毗連，北與翼城、絳縣接壤，西與聞喜交界，西南連接夏縣，東鄰河南省的濟(jì)源市，南與河南省三門峽市的澠池縣、洛陽市的新安縣隔河相望。背靠中條山，西倚中條支脈，東鄰太行山、王屋山，南臨黃河，呈現(xiàn)為三面環(huán)山、一面臨水的丘嶺小盆地特有的多樣化小氣候，屬于暖溫帶半濕潤(rùn)大陸性季風(fēng)氣候，四季分明，年均氣溫13.5 ℃，年均降水量631 mm，無霜期230 d。垣曲縣地勢(shì)西北高、東南低，中部為垣曲山間盆地，呈北寬南窄的葫蘆形展布。境內(nèi)較大河流有毫清河、允西河、板澗河、西陽河、五福澗河等，均由北而南注入黃河。河流兩岸多為狹窄的河谷平原，灌溉方便，人煙稠密，為主要農(nóng)作區(qū)。

1.2 數(shù)據(jù)來源及處理 利用ArcGIS 10.3對(duì)垣曲縣2018年1∶5萬土地利用現(xiàn)狀圖進(jìn)行數(shù)字化，以《全國土地分類》為標(biāo)準(zhǔn)[10]，結(jié)合實(shí)際情況及研究需要，將土地利用現(xiàn)狀分為7類：耕地、園地、林地、草地、居民點(diǎn)及工礦用地、水域、其他用地。DEM數(shù)據(jù)來源于地理空間數(shù)據(jù)云，對(duì)數(shù)據(jù)進(jìn)行拼接裁剪，用Spatial Analyst模塊的Slope、Aspect，生成坡度、坡向圖;Reclassify工具對(duì)高程、坡度、坡向進(jìn)行重分類;將地類圖層與地形因子分級(jí)成果進(jìn)行疊置分析，得到地形因子與土地利用的分布情況（圖1）。

1.3 地形分布指數(shù)

為了消除面積差異量綱的影響[11]，引入地形分布指數(shù)來描述地類在該地形因子級(jí)別上是否呈優(yōu)勢(shì)分布。計(jì)算公式如下：

P=（Sie/Si）×（S/Se）

式中，P為分布指數(shù);Sie為在e地形因子該等級(jí)下的i地類面積;Si為i地類的面積;Se為整個(gè)區(qū)域e地形因子該等級(jí)下的總面積，S為整個(gè)區(qū)域的面積。土地利用類型的地形分布指數(shù)對(duì)于其在地形上的空間分布具有指示意義[12]：P> 表示該地類在該級(jí)地形上的分布高于其他地類在該區(qū)域的平均分布水平，屬于優(yōu)勢(shì)分布，且P值越大，該地類在該地形等級(jí)上優(yōu)勢(shì)度越高;P< 則該地類在該地形因子級(jí)別上的分布屬劣勢(shì)分布，說明該地類不適宜分布在該地形上[13]。

2 結(jié)果與分析

2.1 土地利用現(xiàn)狀分析 由于垣曲縣地處山地，土地利用現(xiàn)狀結(jié)構(gòu)中林地所占面積最大，占整個(gè)垣曲縣土地總面積的57.66%。耕地在土地利用現(xiàn)狀結(jié)構(gòu)中僅次于林地，占總面積的19.15%，其中旱地面積占耕地面積的90%以上。草地面積次于耕地，其中99%以上為其他草地，人工牧草地所占比例較小。園地所占面積較小，僅占土地總面積的0.45%。剩下的依次為居民點(diǎn)及工礦用地、水域、其他用地（表1）。

2.2 高程、坡向和坡度分級(jí)

將研究區(qū)地形因子進(jìn)行分級(jí)，結(jié)果見圖2。垣曲縣海拔最高點(diǎn)為2 306 m，最低點(diǎn)為194 m，平均海拔為1 250.5 m。將研究區(qū)分為7個(gè)高程帶：400 m以下為第1高程帶，400～600 m為第2高程帶，>600～800 m為第3高程帶，>800～1 000 m為第4高程帶，>1 000～1 300 m為第5高程帶，>1 300～1 700 m為第6高程帶，1 700 m以上為第7高程帶。第2高程帶土地面積最大，占總面積的27.27%，第3高程帶次之，占總面積的21.06%，第1高程帶占總面積的15.41%，第3高程帶之后土地面積依次遞減，第7高程帶所占比例最小?？梢娧芯繀^(qū)地形起伏較大，在各高程帶均有土地分布。

垣曲縣坡度為0～63.33°，將坡度劃分為6級(jí)：0～2°為平地，>2°～6°為緩平地，>6°～15°為緩地，>15°～25°為斜坡，>25°～35°為緩陡坡，>35°為陡坡[14]。緩地所占面積最大，占總面積的36.57%，斜坡次之，占總面積的27.06%，平地和陡坡所占比例相當(dāng)，均為4%左右。研究區(qū)地形較為復(fù)雜，以丘陵地帶為主，坡度主要集中在6°～25°。為了便于計(jì)算，將坡向合并為5個(gè)：平地即坡向界線為- 為無坡向;南、西南為陽坡;北、東北為陰坡;西、東南為半陽坡;西北、東為半陰坡[15]。研究區(qū)無坡向較少，僅占總面積的2.59%，陰坡、半陰坡、半陽坡、陽坡依次遞增，所占比例分別為19.76%、23.02%、26.77%、27.86%。

2.3 地形因子影響下的土地利用空間分布特征

從圖3和表2可以明顯看出，第6、第7高程帶主要分布在研究區(qū)西北部，而第1、第2高程帶主要分布在東南部，表明研究區(qū)地勢(shì)西北高、東南低。各種土地類型隨著海拔的變化，具有非常明顯的分布層次。耕地在第1、第2高程帶所占比例占整個(gè)耕地面積的81.05%，在海拔超過600 m的地區(qū)分布指數(shù)均小于 且呈依次遞減趨勢(shì)，不適宜耕種。水域在第1高程帶上分布指數(shù)為2.70，優(yōu)勢(shì)性顯著。居民點(diǎn)及工礦用地在海拔200～800 m的分布指數(shù)大于 在海拔高于800 m地區(qū)，分布指數(shù)隨海拔的升高而降低;海拔較低的地區(qū)具有更有利于人類生活以及從事生產(chǎn)活動(dòng)的地理?xiàng)l件和水土資源條件，熱量和水分等對(duì)于作物的生長(zhǎng)十分有利，適宜耕種;而隨著海拔的升高，氣候條件、地形地貌不利于人類生活和從事生產(chǎn)活動(dòng)。園地主要分布在第1、第2高程帶，占園地總面積的86.17%，海拔高于1 300 m的地區(qū)分布指數(shù)為0，高海拔不利于果樹生長(zhǎng)發(fā)育。林地的面積分布呈現(xiàn)峰值現(xiàn)象，在第4高程帶面積最大，在海拔高于800 m的地區(qū)地形分布指數(shù)大于 在第7高程帶除林地和其他用地外，其他地類分布指數(shù)均小于 且分布指數(shù)高達(dá)5.28，可見林地在第7高程帶上具有較大優(yōu)勢(shì)。其他用地在第7高程帶分布指數(shù)高達(dá)20.2 具有絕對(duì)優(yōu)勢(shì)。這是由于景區(qū)大部分位于高山上，歷山國家級(jí)森林公園就坐落于垣曲縣西北部，整個(gè)景區(qū)受人類經(jīng)濟(jì)活動(dòng)干擾較少，自然生態(tài)環(huán)境保持相對(duì)完好。草地在第2高程帶分布最多，在第4、第5高程帶上分布指數(shù)大于 呈優(yōu)勢(shì)分布，在第4高程帶以上地區(qū)面積比例逐漸下降。

平地、緩平地主要集中在研究區(qū)中部和南部，緩陡坡、陡坡主要分布在研究區(qū)東北部和西南部。平地占總面積比例最小，但在這一地區(qū)耕地所占面積最大，園地分布指數(shù)最高。耕地和園地集中分布在平地、緩平地、緩地上，分別達(dá)到總面積的81.2%和90.51%，分布指數(shù)隨坡度的上升而下降，主要原因是坡度影響土壤肥力，且坡度越高，投入耕種成本越高，加大了人類從事農(nóng)業(yè)生產(chǎn)活動(dòng)的難度。水域在平地分布指數(shù)為4.25，優(yōu)勢(shì)性顯著，隨著坡度升高分布指數(shù)逐漸降低，由于受到水體自身特征的影響，水域多分布在平地上。居民點(diǎn)及工礦用地分布指數(shù)在平地對(duì)比其他地區(qū)達(dá)到最高，分布面積隨坡度上升而減少，因?yàn)槿祟悘氖律詈蜕a(chǎn)活動(dòng)多集中在地勢(shì)平坦的地區(qū)。林地和草地集中分布于坡度大于15°的地區(qū)，達(dá)到總面積的90%以上，分布指數(shù)呈峰值現(xiàn)象，在斜坡地區(qū)達(dá)到峰值，在緩陡坡地區(qū)分布指數(shù)大于 屬于優(yōu)勢(shì)分布。林地和草地在坡度高的地區(qū)優(yōu)勢(shì)性明顯，對(duì)高坡度的環(huán)境具有較強(qiáng)的適應(yīng)性。在陡坡地區(qū)，只有其他用地分布指數(shù)大于 說明在這一坡度界限內(nèi)具有絕對(duì)優(yōu)勢(shì)。

坡向?qū)τ谏降厣鷳B(tài)有著較大的作用，山地的方位對(duì)日照時(shí)數(shù)和太陽輻射強(qiáng)度有影響，對(duì)耕地和園地的分布有很大影響。無坡向地區(qū)土地面積最小，占總面積的2.59%，但在這一地區(qū)耕地、園地分布指數(shù)相比其他坡向?yàn)樽畲笾担页送庵挥性诎腙柶潞完柶路植贾笖?shù)大于 說明耕地、園地在無坡向地區(qū)和陽坡地區(qū)優(yōu)勢(shì)性顯著，這是因?yàn)檫@些地區(qū)可為作物生長(zhǎng)提供足夠的光照，適宜耕種。林地、草地分布受坡向影響不明顯，除在無坡向地區(qū)分布較少外，在其他各坡向分布較為均勻，在半陰坡和陰坡分布指數(shù)大于 說明在這些地區(qū)屬優(yōu)勢(shì)分布。居民點(diǎn)及工礦用地在無坡向和陽坡地區(qū)分布指數(shù)大于 有75%以上分布在這些地區(qū)，這是由于人類從事生產(chǎn)和生活活動(dòng)多集中在平地，且有一定的采光要求，所以這些地區(qū)相對(duì)其他區(qū)域具有優(yōu)勢(shì)性。坡向?qū)λ蚍植紵o明顯影響。

3 結(jié)論

由研究區(qū)土地利用現(xiàn)狀數(shù)據(jù)可知，該地區(qū)土地利用類型主要為林地、耕地、草地，分別達(dá)到總面積的57.66%、19.15%、14.33%，占整個(gè)地區(qū)土地面積的90%以上，其他土地利用類型所占比例較小。土地利用類型的空間分布明顯受到地形因素的影響和制約，土地利用和地形因子具有顯著的相關(guān)性。耕地和園地主要分布在高程小于1 100 m且坡度小于15°的平地、緩平地和緩地上，在無坡向和陽坡地區(qū)分布較多;林地主要分布在海拔500 m以上且坡度大于15°的斜坡和陡坡上，受坡向影響不大;草地主要分布在海拔800 m以下且坡度大于15°的斜坡和陡坡上，陽坡上分布大于陰坡;居民點(diǎn)及工礦用地多分布于海拔低于800 m且坡度小于15°的平地和緩地上，在無坡向和陽坡上分布較多;水域主要分布在海拔低于500 m且坡度小于6°的緩地和緩平地上，受坡向影響不大。綜上所述，在低海拔、地勢(shì)平緩、陽光充足且水源豐富的地區(qū);適宜人類進(jìn)行各種生產(chǎn)和生活活動(dòng)，耕地和園地多分布于這一類地區(qū);而在高海拔、地勢(shì)陡峭的陰面地區(qū)，惡劣的外部環(huán)境限制人類活動(dòng)，故林地和草地多分布于這些地區(qū)。

參考文獻(xiàn)

[1] 劉紀(jì)遠(yuǎn).中國資源環(huán)境遙感宏觀調(diào)查與動(dòng)態(tài)研究[M].北京：中國科學(xué)技術(shù)出版社，1996.

[2] 李碩，沈占鋒，劉克儉，等.大清河流域土地利用變化的地形梯度效應(yīng)分析[J].農(nóng)業(yè)工程學(xué)報(bào)，202 37（5）：275-284.

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

[4] 李玲，何方，吳楠，等.淮河上游山丘區(qū)不同地形因子下的土地利用變化研究[J].安徽農(nóng)業(yè)科學(xué)，2009，37（34）：16973-16975.

[5] 李益敏，王東馳，袁靜，等.基于地形梯度的高山峽谷區(qū)土地利用時(shí)空格局及功能分區(qū)[J].水土保持通報(bào)，2020，40（6）：303-31 331.

[6] 郝雅珺，常慶瑞，劉海飛.基于GIS的縣域土地利用與地形因子關(guān)系研究：以陜西省黃陵縣為例[J].水土保持通報(bào)，201 32（6）：134-137，141.

[7] 楊洋，畢如田.黃土高原典型地貌類型的土地利用分形特征[J].地理與地理信息科學(xué)，201 27（1）：101-104.

[8] 呂靜，常慶瑞.黃土高原土石山區(qū)土地利用空間分布研究[J].干旱區(qū)資源與環(huán)境，201 26（8）：93-96.

[9] 趙益禎，曹建農(nóng)，張曉棟，等.寶塔區(qū)土地利用的地形梯度效應(yīng)與空間格局[J].干旱區(qū)地理，2020，43（5）：1307-1315.

[10] 孫丕苓，許月卿，王數(shù).環(huán)京津貧困帶土地利用變化的地形梯度效應(yīng)分析[J].農(nóng)業(yè)工程學(xué)報(bào)，2014，30（14）：277-288.

[11] 羅丹，常慶瑞.渭北黃土臺(tái)塬區(qū)土地利用空間分布的地形特征研究[J].西部大開發(fā)（土地開發(fā)工程研究），2018，3（1）：1-8.

[12] 代冬梅，吳國平，司惠超.基于DEM的瀏陽市土地利用空間格局分析[J].安徽農(nóng)業(yè)科學(xué)，2009，37（3）：1231-1234.

[13] 馬廷剛，秦占飛，常慶瑞，等.縣域土地利用的地形特征分析：以陜西省戶縣為例[J].干旱地區(qū)農(nóng)業(yè)研究，201 29（6）：191-194，221.

[14] 羅丹，王濤，常慶瑞.縣域農(nóng)村居民點(diǎn)適宜性評(píng)價(jià)：以陜西省隴縣為例[J].中國農(nóng)學(xué)通報(bào)，2019，35（14）：157-164.

[15] 林超，李昌文.陰陽坡在山地地理研究中的意義[J].地理學(xué)報(bào)，1985，40（1）：20-28.