基于均值變點分析方法的起伏度研究

魏蘭 解博 唐夢鴿 羅明良

摘 要：基于甘孜地區的SRTM3—DEM數據，運用ArcGIS提取遞增窗口3×3，5×5，7×7，9×9，…，75×75對應地勢起伏度，并通過均值變點分析法確定最佳統計單元，輸出地勢起伏度分級圖。結果表明，甘孜地區地形起伏度最佳分析窗口為29×29，對應的最佳統計面積為6.8km2。地勢起伏度分級圖中以中大起伏山地的分布最廣，與甘孜地區多山地貌特征相符，可能受地質構造運動內力作用與冰川、流水侵蝕等外力的共同作用。

關鍵詞：地勢起伏度;ArcGIS;均值變點;甘孜

中圖分類號 P208 文獻標識碼 A文章編號 1007-7731（2019）17-0106-04

Study on Relief Amplitude Based on Mean Change Point Analysis Method

——Taking the Ganzi Area as an Example

Wei Lan1，2 et al.

（1Land and Resources School，China West Normal University，Nanchong 637009，China;2Institute of Landform Surface processes and Environmental Changes，China west normal University，Nanchong 637009，China）

Abstract：Based on SRTM3-DEM data in Ganzi area，ArcGIS is used to extract the incremental windows 3×3，5×5，7×7，9×9，...，75×75 corresponding to the relief amplitude，and the best is determined by mean point analysis.The statistical unit outputs a gradient map of the topography.The results show that the optimal analysis window for relief amplitude in Ganzi area is 29×29，and the best statistical area is 6.8km2.The relief amplitude map is the most widely distributed in the undulating hills，which is consistent with the mountainous features of the Ganzi area.It may be affected by the internal forces of geological tectonic movements and external forces such as glaciers and water erosion.

Key words：Relief amplitude;ArcGIS;Mean change point;Ganzi

地表是人類生產生活的重要場所，地表形態的高低，對人類生產生活、人口分布、規劃建設等的影響尤為重要[1-3]。通常地表起伏走向，我們用地勢高低來形容，地勢起伏度是指研究區域的高程差，地勢起伏度可以反映研究區地表起伏狀況、地貌特征，同時也可以為地形地貌劃分提供定量依據[4-6]。涂漢明等通過大量實驗，證明中國地勢起伏度的最佳統計單元面積為21km2，對于不同規模的地勢起伏度研究時，分別適用2，6，16，20，22km2等5種最佳分析區域[7]。張錦明等通過不同區域、不同地貌的64幅樣區研究了地貌形態與地勢起伏研究的相關關系[8]。唐飛等基于1∶25萬新疆DEM數據的研究表明，適合當地地勢研究的最佳分析單元面積為4km2[9];在相同1∶25萬比例尺尺度下，不同學者對各自所研究區域進行了研究，張磊研究得出9.61km2的范圍適合京津冀地區;王雷研究認為昆明地區最佳分析區域為16km2;郎玲玲等研究認為福建省最佳統計單元為4.41km2[10-13]。不同地區、不同地貌所需采用的分析單元不同，研究區最佳統計單元的確定是對該地區地形起伏度研究正確開展的關鍵，均值變點法能快速確定研究數據拐點所在，即地勢起伏度中對應的最佳統計單元。目前，該方法已廣被廣泛使用，在地勢起伏度的研究中具有獨特的優勢[3，14]。

甘孜藏族自治州，簡稱甘孜州，地理位置特殊，地貌類型復雜，高山峽谷分布，地勢起伏大，與我國其他地區地形地貌差異巨大。而目前多數學者研究區域多分布于我國一、二級階梯，地勢相對比較平緩，而甘孜所處地區地形起伏較大，地貌類型復雜，研究相對較少。為此，本研究以甘孜州區域SRTM3—DEM為數據源，運用統計法、均值變點法等對川西的甘孜地區地貌的地勢起伏度進行了研究，輸出甘孜地勢分級圖，對川西甘孜地區的地形進行了定量分析。

1 研究區域與數據

1.1 研究區概況 甘孜藏族自治州（簡稱甘孜州），位于四川西南部，地處東經98～102°、北緯28～34°。地勢整體由西向東南傾斜，地貌類型大致可分，西北部高原區、折多山以西山原地區、大雪山以東的高山峽谷區3個部分。高原地區以丘狀高原為主，金沙江、大渡河等河流支流上游遷回其間，切割較淺，河谷寬平，相對高差較小，頂面仍較寬坦，都比較完好地保存有準平原時期的“平坦”地面。高山峽谷區，受冰川流水主要作用，山高谷深組合形成高差極大的高山峽谷區，其中不少山峰海拔高于6000m，如雀兒山、貢嘎嶺等。

1.2 研究數據 甘孜DEM數據源于CGIAR（國際農業研究磋商組織）網站（http：//srtm.csi.cgiar.org/）下載，其水平基準為WGS84_1984_Albers，按照所需區域進行拼接投影裁剪為所需區域，基準面選擇D_WGS_1984。

2 研究方法

2.1 地形起伏度 地形起伏度是研究區域最高和最低點的高差，研究地勢起伏度可以為地形地貌定量描述提供依據[15]。地勢起伏度計算表達式為[15]：

[Δh=hxy，max-hx，y，min]? （1）

[S=k=1nΔhx，y] （2）

[M=S/N] （3）

式中，hx，y為領域內像元的高程值，hx，y，max 和hxi，y，min 分別代表目標柵格最大和最小高程值（x，y=1，2，3，…，n）;Δh表示高程值差;N為柵格的個數;S為柵格的起伏度和;M為平均起伏度值。

2.2 均值變點法 均值變點法是一種數理統計方法，研究數據擬合曲線應為非線性，并且只有1個變點。在地勢起伏度研究中，該方法可以快速確定地勢平均起伏度的增長速度拐點，即最佳統計單元[16，17]。在該系統模型統計中，當某個時刻突然變化，則該時刻為所求的變點[15]。計算過程大致如下[18]：

將樣本數據i=2，…，N，分為X1，X2，…，Xi-1和Xi，Xi+1，…，Xt。

[Si=t=1i-1（Xt-Xi1）2+t=iN（Xt-Xi2）2] （4）

[S=t=iN（Xt-X）2] （5）

式中：[Xi1]和[Xi2]代表每段樣本的算術平均值，Si為樣本數據分段后統計量，S表示原始樣本所得統計量，t表示目標柵格總數量。當S與Si的差值達到最大時，則該時刻對應為研究區地勢起伏度變化拐點，為所求最佳統計單元。

3 結果與分析

本文基于ArcMap的toolbo模塊，利用Focal Statistics和Neighborhood工具模型，計算甘孜地區的地勢起伏度。以目標柵格為中心，選取矩形分析窗口，以3×3、5×5，7×7，…為半徑，提取得到其各窗口地勢起伏度平均值，將該值與各單元窗口對應面積作擬合[19，20]。

表1為甘孜地區不同網格地勢起伏度值，表格羅列從3×3…75×75遞增窗口及對應地勢起伏度值，平均地勢起伏度最低為99.70m，最高達1146.20m。對應面積分別為7.29×104m2和4556.25×104m2，通過將平均地勢起伏度與格網面積進行擬合（圖1），可以看出擬合曲線總體呈增長趨勢，R2值為0.98，說明該曲線擬合可信度較高，從曲線變化情況可以看出，曲線開始呈急劇增長，說明隨窗口大小增大，平均地勢起伏度增長速度明顯，當單元面積值大約介于（500～2000）×104m2時，擬合曲線增加幅度稍微變緩;當單元面積值大于2000×104m2時，曲線相對變化更加平緩，平均地勢起伏度增長速度明顯變慢。

以上曲線拐點受主觀影響，可能存在出入，故由表1數據運用均值變點分析法進行科學曲線擬合，尋找拐點位[10]。根據表1的平均起伏度數據，構建序列并和實際統計值S與樣本統計值Si差值進行曲線擬合（表2、圖2）。

根據均值變點法分析，擬合曲線只存在一個拐點，拐點為點數序列為14時，S和Si的差值達到最大，對應表1分析窗口第14個為29×29，故確定甘孜地區地勢起伏度最佳統計面積對應為6.8km2。

基于以上地勢起伏度最佳單元的確定，將地勢起伏度按數字地貌制圖標準分別以30、70、200、500、1000、2500m為界限劃分為7種類型，為平原、臺地、丘陵、小起伏山地、中起伏山地、大起伏山地和極大起伏山地[21，22]。運用ARCGIS軟件按照所求最佳29×29單元求得地勢起伏圖，并按以上標準劃分輸出甘孜地區地勢分級圖（圖3）。

4 結論與討論

（1）通過ARCGIS、均值變點法等分析得出，甘孜地區最佳分析窗口為29×29，最佳統計面積為6.8km2，將起伏度按標準分級，可以看出，中大起伏山地分布最廣，與當地山高谷深地埋特征符合。

（2）研究區地處四川西南部，處于青藏高原與四川盆地過渡，且大部分分布于橫斷山區，中大起伏山地的分布主要可能受地質構造作用與冰川流水的外力侵蝕影響。

（3）甘孜地區地勢起伏研究，可以為西南其他地區地勢起伏度研究提供參考，進一步探討地域性地勢起伏度對實際生活的影響。

參考文獻

[1]朱紅春，陳楠，劉海英，等.自1∶10000比例尺DEM提取地形起伏度——以陜北黃土高原的實驗為例[J].測繪科學，2005，30（04）：86-88，7.

[2]MESSERLI B，WEINGARTNER V R.A New Typology for Mountains and Other Relief Classes：An Application to Global Continental Water Resources and Population Distribution[J].Mountain Research and Development，2009，21：34-45

[3]張靜靜，李艷紅，朱連奇.豫西山區縣域地形起伏度與人口、經濟活動分布的關系[J].地域研究與開發，2019，38（02）：55-60.

[4]程維明，周成虎.多尺度數字地貌等級分類方法[J].地理科學進展，2014，33（01）：23-33.

[5]龔熊波，胡紅兵.黃岡市地形起伏度及其與人口、經濟的相關性研究[J].綠色科技，2018（02）：151-154.

[6]后立勝，蔡運龍.土地利用/覆被變化研究的實質分析與進展評述[J].地理科學進展，2004，23（06）：96-104.

[7]涂漢明，劉振東.中國地勢起伏度最佳統計單元的求證[J].湖北大學學報（自然科學版），1990，12（03）：266-271.

[8]張錦明，游雄.地形起伏度最佳分析區域研究[J].測繪科學技術學報，2011，28（05）：369-373.

[9]鄭颯颯，楊佑洪，劉志，等.四川省冰雹分布與地形因子關系分析[J].氣象科技，2018，46（06）：1280-1286.

[10]唐飛，陳曦，程維明，等.基于DEM的準噶爾盆地及其西北山區地勢起伏度研究[J].干旱區地理，2006，29（03）：388-392.

[11]王雷，朱杰勇，周燕.基于1∶25萬DEM昆明地區地貌形態特征分析[J].昆明理工大學學報（理工版），2007，32（01）：6-9，14.

[12]張磊.基于地形起伏度的地貌形態劃分研究[D].石家莊：河北師范大學，2009.

[13]郎玲玲，程維明，朱啟疆，等.多尺度DEM提取地勢起伏度的對比分析——以福建低山丘陵區為例[J].地球信息科學，2007，9（06）：1-6，135-136.

[14]王巖，劉少峰.基于DEM的青海貴德地區地形起伏度的研究[J].地質通報，2008，27（12）：2117-2121.

[15]楊艷林，邵長生.長江中游地形起伏度分析研究[J].人民長江，2018，49（02）.

[16]王玲，同小娟.基于變點分析的地形起伏度研究[J].地理與地理信息科學，2007，23（06）：65-67.

[17]陳學兄，常慶瑞，畢如田，等.地形起伏度最佳統計單元算法的比較研究[J].水土保持研究，2018，25（01）：52-56.

[18]項靜恬.非線性系統中數據處理的統計方法[J].北京：科學出版社，1997.

[19]張軍，李曉東，陳春艷，等.新疆地勢起伏度的分析研究[J].蘭州大學學報（自然科學版），2008（S1）：10-13，19.

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

[21]王康，何俊仕，于德浩，等.采用ArcGIS平臺的地勢起伏度自動提取技術研究[J].沈陽理工大學學報，2013，32（02）：63-67.

[22]程維明.中國1∶100萬數字地貌圖研究;proceedings of the 中國地理學會百年慶典，F，2009[C]. （責編：張宏民）