黎　武,王汝蘭,慕　凱,余林艷,賀春明,徐　珍

(西華師范大學 國土資源學院，四川 南充 637009)

?

基于DEM的湖北省蘄春縣地形起伏度分析

黎武,王汝蘭,慕凱,余林艷,賀春明,徐珍

(西華師范大學 國土資源學院，四川 南充 637009)

基于大別山區(qū)的湖北省蘄春縣ASTER GDEM(空間分辨率30m×30m)數(shù)據(jù)，運用ArcGIS 10.2的窗口分析法依次計算3×3、5×5、7×7、…、65×65窗口下的地形起伏度，最后用均值變點分析法確定最佳分析窗口，利用最佳分析窗口制作蘄春縣地形起伏度分級圖。研究結(jié)果表明：蘄春縣最佳分析窗口為21×21，最佳統(tǒng)計面積為106.20×104m2，地形起伏度可分為5級，其中丘陵分布最廣占總面積的33.75%，平原占16.61%，臺地占32.96%，小起伏山地占16.67%，中起伏山地分布最少占0.02%；地形的起伏是造成水土流失的重要因素，蘄春縣地形起伏度較小，地勢低平，地形條件對水土流失影響較輕；地形起伏導致土壤侵蝕退化，土壤肥力下降，制約農(nóng)業(yè)的發(fā)展。針對起伏度分級的不同地貌類型區(qū)域，因地制宜的進行農(nóng)業(yè)生產(chǎn)活動。臺地和丘陵地區(qū)坡改梯，防治坡耕地水土流失，中起伏和小起伏山地列入保護用地，適當發(fā)展經(jīng)濟林，防止?jié)撛诘乃亮魇А?/p>

蘄春縣；DEM；地形起伏度；均值變點分析法；水土流失

地形起伏度 (Relief Amplitude)是指某一確定面積內(nèi)最高點和最低點之高差[1]。它反映區(qū)域內(nèi)地面的起伏特征，是定量描述地貌形態(tài)、劃分地貌類型的重要指標[2]。地形起伏度是重要的地形因子，它不僅反映了地表切割剝蝕程度，也是導致水土流失最直接的因素[2]。近年來，隨著3S技術(shù)的發(fā)展，用數(shù)字高程模型(DEM)研究區(qū)域地形起伏度逐漸興起[3-5]，地形起伏度已經(jīng)應用于人口分布(人口分布與經(jīng)濟發(fā)展)[5-8]土地利用時空變化[9]，水土流失定量評價[10]，地質(zhì)災害評價(滑坡)[11-12]，生態(tài)環(huán)境(人居環(huán)境)[13-14]等方面。

目前對地形起伏度的研究大多集中在較大范圍[4，15-17]。對局部小范圍的研究相對較少，如陳珂[18]以華鎣市DEM、丁賢法[19]以云南省富寧縣DEM為數(shù)據(jù)源，分析地形起伏度。局部小范圍特別是縣域范圍下的研究也有重要的意義。本文以ASTER GDEM數(shù)據(jù)為基礎(chǔ)數(shù)據(jù)源，湖北省蘄春縣為研究區(qū)，采用GIS空間分析(Spatial Analyst)中的窗口分析法提取蘄春縣地形起伏度，最后利用均值變點分析法確定最佳分析窗口，依據(jù)數(shù)字地貌制圖規(guī)范，制作地形起伏度分級圖，據(jù)分級圖對該縣地形起伏度進行初步分析。通過分析蘄春縣地形起伏度的空間分布規(guī)律，不僅可以認識該地區(qū)宏觀地形特征與水土流失情況，也能夠為土地資源的合理利用和水土流失的防范治理提供重要的理論依據(jù)。

1　研究區(qū)概況

蘄春縣在115°12′～115°55′E、29°59′～30°41′N之間。位于湖北省東南部，大別山南麓，東臨安徽省，長江中游北岸，隸屬黃岡市(如圖1所示)。地勢北高南低，地貌復雜。北部為山區(qū)，中部為丘陵地帶，南部為平原。屬于亞熱帶季風氣候，四季分明，氣候溫和，光照充足，雨量充沛，四季分明。歷年平均無霜期249.1天，日照時數(shù)2 025.8小時，年平均氣溫16.8 ℃，年降水量1 341.7 mm[20]。蘄春版圖面積2 398平方千米，下轄15個鄉(xiāng)鎮(zhèn)辦、2個省級園區(qū)(開發(fā)區(qū))和1個國家級濕地公園，人口約101.1萬(2014年統(tǒng)計)。

圖1　研究區(qū)地理位置

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

2.1數(shù)據(jù)來源與處理

本文采用ASTER GDEM數(shù)據(jù)，該數(shù)據(jù)來自于地理空間數(shù)據(jù)云平臺網(wǎng)站(http：//www.gscloud.cn/)，該數(shù)據(jù)是在2009年6月由日本經(jīng)濟產(chǎn)業(yè)省(METI)和美國國家航空航天局(NASA)共同發(fā)布?？臻g分辨率為30 m×30 m，數(shù)據(jù)覆蓋范圍為83°S到83°N之間的所有陸地區(qū)域，達到了地球陸地表面覆蓋率的99%，其平面基準為WGS84，置信度為95%[21]。在ArcGIS10.2空間分析中，對獲取的蘄春縣DEM數(shù)據(jù)進行拼接處理，然后在裁剪工具(Clip)中，用從國家基礎(chǔ)地理信息中心下載1∶400萬的蘄春縣矢量邊界圖去裁剪拼接好的DEM數(shù)據(jù)，得到蘄春縣縣域范圍內(nèi)的DEM數(shù)據(jù)。

2.2研究方法

2.2.1窗口分析法提取地形起伏度

在ArcGIS 10.2的空間分析(Spatial Analyst)中，使用窗口遞增分析法提取蘄春縣地形的起伏度。窗口類型為矩形，窗口依次為3×3、5×5、7×7…，移動步距為2，終止窗口為65×65。用柵格計算器(Raster Calculator)計算每個分析窗口內(nèi)最大高程與最小高程之差，每個矩形窗口內(nèi)的高差作為起伏度，窗口下平均值即是平均地形起伏度，由此計算出每個分析窗口的起伏度。每個窗口內(nèi)的最大高程圖層與最小高程圖層的差值所得出的新圖層為該窗口內(nèi)的地形起伏度圖層。計算公式為[16]：

(1)

式中：ΔHi表示以第i個矩形窗口內(nèi)的地形起伏度，Hmax、Hmin分別表示該窗口內(nèi)的最大、最小高程值。

2.2.2均值變點分析法

基于DEM提取地形起伏度，關(guān)鍵在于確定最佳分析窗口[1]。本文使用數(shù)理統(tǒng)計上的均值變點分析法確定最佳分析窗口，該方法的原理是變點的存在會使原始樣本的統(tǒng)計量S與樣本分段后的統(tǒng)計量Si之間的差距增大。令i=2,3…,n，對每個i將樣本分為兩段：X1, X2,…,Xi-1和Xi,Xi+1,…,Xn，計算步驟大致如下[16，22]。

(a)計算每段樣本的算術(shù)平均值Xi1、Xi2及統(tǒng)計量Si：

(2)

(b)計算統(tǒng)計量S：

(3)

(c)計算期望值E(S-Si)：

(4)

3　結(jié)果與分析

3.1不同窗口與平均起伏度的對應關(guān)系

利用GIS的窗口分析法計算蘄春縣地形起伏度，依次計算3×3、5×5、7×7、…、65×65窗口下的地形起伏度，取每個分析窗口的地形起伏度的平均值。窗口大小、窗口面積、平均起伏度的關(guān)系(如表1所示)。

表1　分析窗口面積與平均起伏度對應關(guān)系圖

對每個分析窗口下的平均地形起伏度與各分析窗口面積進行對數(shù)擬合(y=AlnX+B)，得出的擬合曲線(圖2)擬合系數(shù)為0.95，擬合效果較好。由圖2可以看出：隨著窗口面積的增加，地形起伏度的值也會增加，最終會趨向穩(wěn)定，由快到慢的變化過程中會存在一個最佳統(tǒng)計單元。圖2中窗口面積在0～50×104m2平均地形起伏度迅速增大，之后增速減緩趨于平穩(wěn)。擬合曲線由陡變緩的點對應的窗口面積就是所要尋找的最佳統(tǒng)計單元。

圖2　分析窗口面積與平均起伏度擬合曲線

用窗口遞增法得到的每個分析窗口的地形起伏度隨分析窗口的變化曲線，可以從曲線的變化規(guī)律中確定最佳統(tǒng)計單元。先前的研究多采用目視判斷，但是這樣往往會受到主觀因素的影響，影響結(jié)果的準確性。為了克服主觀因素的影響，本文采用均值變點分析法確定最佳統(tǒng)計單元。

3.2均值變點分析法確定最佳統(tǒng)計單元

用均值變點法對窗口分析法得到的數(shù)據(jù)結(jié)果進行處理。處理步驟如下：

(a)首先用表1中的數(shù)據(jù)結(jié)果計算各分析窗口下單位面積的平均地形起伏度(即平均地勢起伏度與窗口面積的比值)，記為Ti(i=1,3,5,7,9…65)。

(b)然后對T以ln為底取對數(shù)構(gòu)建樣本序列X(X={Xi,i=3,5,7…,63})。用均值變點法公式計算樣本序列X的統(tǒng)計量S與Si值，S值為28.95，并計算樣本序列的Si與S-Si，其Si、S-Si值見表2。

表2　均值變點分析統(tǒng)計結(jié)果

對S-Si值進行擬合，擬合曲線見圖3。從圖3可以看出，在第10個點時S-Si值達到最大，即為所求的變點(由陡變緩的點)。該點所對應的分析窗口大小為21×21，矩形窗口面積106.20×104m2，即蘄春縣最佳統(tǒng)計面積為106.20×104m2。

圖3　S與Si差值的變化曲線

3.3地形起伏度分析

以分析窗口為21×21，面積為106.20×104m2作為蘄春縣起伏度的最佳分析窗口，用最佳分析窗口提取蘄春縣地形起伏度為0～569 m(如圖4所示)。

根據(jù)數(shù)字地貌制圖規(guī)范，地形起伏度可劃分為7級，即：平原(<30 m)、臺地(30～70 m)、丘陵(70～200 m)、小起伏山地(200～500 m)、中起伏山地(500～1000 m)、大起伏山地(1000～2 500 m)和極大起伏山地(>2 500 m)[16,23]。利用最佳窗口(21×21)提取的蘄春縣地形起伏度0～569 m(圖4)，依據(jù)此分類標準，可將蘄春縣地貌劃分為平原、臺地、丘陵、小起伏山地、中起伏山地5種類型。在ArcGIS 10.2的空間分析(Spatial Analyst)中，用重分類(Reclassify)工具對蘄春縣地形起伏度進行分級，制作地形起伏度分級圖見圖5。

圖4　用最佳窗口(21×21)提取的蘄春縣地形起伏圖

圖5　蘄春縣地勢起伏度分級圖

從空間分布上來看：蘄春縣地處大別山區(qū)，整個地勢東北高、西南低，由東北向西南漸次傾斜，地形起伏度也呈現(xiàn)東北向西南傾斜的空間分布規(guī)律。據(jù)圖5可以看出，該縣地形起伏度空間分布整體趨勢為北部高于南部，東部高于西部。起伏度<70 m的主要分布西部，包括漕河鎮(zhèn)、赤東鎮(zhèn)、葛州鎮(zhèn)、八里湖農(nóng)場、橫車鎮(zhèn)、彭思鎮(zhèn)、管窯鎮(zhèn)、劉河鎮(zhèn)北部、青石鎮(zhèn)北部。這些區(qū)域海拔較低，起伏度較小，地形主要以平原為主，是蘄春縣人口分布最多的區(qū)域。70 m<起伏度<1000 m主要分布在與英山縣、武穴市、黃梅縣和安徽省的交界處，包括譚林鎮(zhèn)、大同鎮(zhèn)、向橋鄉(xiāng)、青石鎮(zhèn)南部、劉河鎮(zhèn)南部、株林鎮(zhèn)、漕河鎮(zhèn)東部。這些區(qū)域地處大別山區(qū)，地勢較高，起伏度較大。大起伏度的山區(qū)地表起伏，地形破碎，溝壑縱橫，土壤結(jié)構(gòu)松散，抗沖性弱，蘄春縣屬于亞熱帶季風氣候區(qū)，降水豐沛，加之人類活動的干擾，易對當?shù)剞r(nóng)業(yè)活動、工業(yè)生產(chǎn)、交通以及地質(zhì)災害防治、生態(tài)環(huán)境治理產(chǎn)生較大影響。由于蘄春縣最佳分析窗口為21×21，起伏度為0～569 m，依據(jù)數(shù)字地貌制圖規(guī)范的分類標準，所以該區(qū)域沒有出現(xiàn)大起伏山地和極大起伏山地。

4　結(jié)論與討論

(1)基于均值變點分析法提取蘄春縣地形起伏度，最佳統(tǒng)計面積為106.20×104m2，地形起伏度可分為5級，地形起伏度30 m～70 m的區(qū)域占總面積的32.96%，起伏度70 m～200 m占33.75%。地形起伏度>200 m占全縣面積的16.69%，起伏度大的區(qū)域較少，水土流失情況較輕。地形起伏度<70 m的集中在西部，小起伏與中起伏山地集中在東部，東部地形起伏度大于西部，北部大于南部。

(2)采用GIS的窗口分析方法提取蘄春縣的地形起伏度，均值變點分析方法確定最佳統(tǒng)計單元，這種方法具有簡單、準確、直觀、快速的優(yōu)點。用最佳統(tǒng)計單元提取的蘄春縣(0～569 m)的起伏度圖，ArcGIS軟件的重分類工具制作的起伏度分級圖能清晰的反映蘄春縣地勢起伏的空間變化和整體變化。

(3)小區(qū)域較高分辨率的DEM研究地形起伏度有待深入，本文在縣域尺度下，選取空間分辨率為(30 m×30 m)的ASTER GDEM數(shù)據(jù)來提取蘄春縣的地勢起伏度。不僅能較好的認識蘄春縣的地貌形態(tài)，也能為縣域尺度下地形起伏度的進一步研究提供借鑒意義。

(4)地形是最重要的自然地理要素，地形起伏度可以定量描述地形地貌特征。地表起伏度的大小反映了地表的平整與破碎程度，分析地形起伏度空間分布對我們認識區(qū)域人口分布、經(jīng)濟活動、工農(nóng)業(yè)活動、自然災害的規(guī)律有著重要的意義。地形的起伏不僅是導致水土流失最直接的原因，而且也是我們認識區(qū)域地貌類型特征的重要評價因子。因此，地形起伏度因子的研究對區(qū)域水土流失宏觀評價和地貌類型的認識有著重要意義。

(5)由于認識水平有限和數(shù)據(jù)的限制，本文主要選取區(qū)域宏觀地形因子中的地形起伏度來分析蘄春縣的地貌類型。在今后的研究中不僅要選取更多的區(qū)域宏觀地形因子，如地表粗糙度、地表切割度等因子，而且也應當分析地形起伏度對地面相關(guān)事物的影響，以及其在社會經(jīng)濟發(fā)展中的意義。

[1]涂漢明, 劉振東. 中國地勢起伏度最佳統(tǒng)計單元的求證[J].湖北大學學報(自然科學版),1990,12(03):266-271.

[2]王玲,呂新. 基于DEM的新疆地勢起伏度分析[J].測繪科學,2009,34(01):113-116.

[3]王巖,劉少峰. 基于DEM的青海貴德地區(qū)地形起伏度的研究[J].地質(zhì)通報,2008,27(12):2117-2121.

[4]陳學兄,常慶瑞,郭碧云,等. 基于SRTM DEM數(shù)據(jù)的中國地形起伏度分析研究[J].應用基礎(chǔ)與工程科學學報,2013,21(04):670-678.

[5]封志明, 唐焰, 楊艷昭, 等. 中國地形起伏度及其與人口分布的相關(guān)性[J].地理學報,2007,62(10):1073-1082.

[6]楊雪婷, 蘇維詞. 云南省縣域地形起伏度與人口分布的耦合關(guān)系研究[J].重慶師范大學學報(自然科學版),2016,33(01):123-129.

[7]周自翔, 李晶, 任志遠. 基于GIS的關(guān)中-天水經(jīng)濟區(qū)地形起伏度與人口分布研究[J].地理科學,2012,32(08):951-957.

[8]于慧, 鄧偉, 劉邵權(quán). 地勢起伏度對三峽庫區(qū)人口及經(jīng)濟發(fā)展水平的影響[J].長江流域資源與環(huán)境,2013,22(06):686-690.

[9]胡業(yè)翠, 鄭新奇, 吳家章. 土地利用時空變化與地形起伏度的關(guān)系——以浙江省溫州都市圈為例[J].國土資源科技管理,2011,28(04): 1-6.

[10]劉新華, 楊勤科, 湯國安. 中國地形起伏度的提取及在水土流失定量評價中的應用[J]. 水土保持通報,2001,21(01):57-62.

[11]郭芳芳, 楊農(nóng), 孟暉, 等. 地形起伏度和坡度分析在區(qū)域滑坡災害評價中的應用[J]. 中國地質(zhì),2008,35(01):131-143.

[12]郭芳芳, 楊農(nóng), 張岳橋, 等. 基于GIS的滑坡地質(zhì)災害地貌因素分析[J]. 地質(zhì)力學學報,2008,14(01):87-96.

[13]王慶安, 張翔, 毛竹. 汶川地震生態(tài)環(huán)境破壞嚴重區(qū)地形高程和起伏度分布特點[J]. 科技導報,2009,53(16):53-59.

[14]查瑞生, 陳夢琳, 趙曉雪. 基于地形起伏度的南川區(qū)人居環(huán)境地形適宜性評價[J]. 西南大學學報(自然科學版),2014,36(10):150-156.

[15]張偉, 李愛農(nóng). 基于DEM的中國地形起伏度適宜計算尺度研究[J]. 地理與地理信息科學,2012,28(04):8-12.

[16]韓海輝, 高婷, 易歡, 等. 基于變點分析法提取地勢起伏度——以青藏高原為例[J]. 地理科學,2012,32(01):101-104.

[17]王玲, 同小娟. 基于變點分析的地形起伏度研究[J]. 地理與地理信息科學,2007,23(06):65-67.

[18]陳珂, 黃小羽, 羅明良. 基于ASTER GDEM的縣域尺度地形起伏度分析——以華鎣市為例[J].遙感信息,2014,29(05):69-100.

[19]丁賢法. 基于SRTM DEM與變點分析法的云南省富寧縣地貌形態(tài)研究[J].測繪與空間地理信息,2014,37(11):98-100.

[20]曾春陽, 唐代生, 唐嘉鍇. 3S技術(shù)在湖北蘄春縣森林資源信息提取中的應用研究[J].湖北林業(yè)科技,2008,149(01):26-38.

[21]康曉偉, 馮鐘葵. ASTER GDEM數(shù)據(jù)介紹與程序讀取[J].遙感信息,2011,69(06):69-72.

[22]項靜恬，史久恩. 非線性系統(tǒng)中數(shù)據(jù)處理的統(tǒng)計方法[M].北京:科學出版社,1997.

[23]中國科學院地貌圖集編輯委員會. 中華人民共和國地貌圖集(1∶100萬)[M]. 北京: 科學出版社, 2009.

責任編輯王菊平

A DEM-based analysis of topographic relief in Qichun County of Hubei Province

LI Wu, WANG Ru-lan, MU Kai, YU Lin-yan, HE Chun-ming, XU Zhen

(School of Land and Resources, China West Normal University, Nanchong 637009，Sichuan，China)

Based on the ASTER GDEM (spatial resolution 30 m*30 m) data of the Dabie Mountain Area in Qichun County, Hubei Province, the window analysis method of ArcGIS 10.2 was used to calculate the 3×5, 3×5, 7×7, … 65×65 window of the terrain respectively, and then the method of mean change point analysis to determine the best analysis window which was used to make the Qichun County topographic relief map. The results showed that the best analysis window was 21*21 Qichun, best statistical area 106.20×104m2. Relief degree could be divided into 5 levels, including hills, the most widespread accounting for 33.75% of the total area; plain, accounting for 16.61%; platform, accounting for 32.96%; small rolling mountain, accounting for 16.67%; undulating mountain, the least distribution accounting for 0.02%. The ups and downs of terrain is one of the major factors causing soil and water loss. Qichun’s relief degree is quite low and thus its terrain condition has a minor impact on soil and water loss. Relief leads to degradation of soil erosion and soil fertility decline which restrict the development of agriculture. Agricultural production activities should be carried out differently according to the landscapes with different reliefs. Platforms and hills should be converted into terraces to prevent farmland soil and water loss. The ups and downs and small undulating hills should be protected. Economic forest should be developed properly to prevent potential soil erosion.

Qichun County; digital elevation model (DEM); relief; change point analysis; soil erosion

P942；TP79

A

1003-8078(2016)03-0088-07

2015-09-01

10.3969/j.issn.1003-8078.2016.03.22

黎武，男，湖北黃岡人，碩士研究生，主要研究方向為土壤侵蝕與水土保持。