常直楊，孫偉紅

(1.南京旅游職業(yè)學(xué)院，江蘇 南京 211100； 2.江蘇省測(cè)繪工程院，江蘇 南京 210013)

白龍江流域地勢(shì)起伏度與微觀地形因子之間的關(guān)系探討

常直楊1，孫偉紅2

(1.南京旅游職業(yè)學(xué)院，江蘇 南京 211100； 2.江蘇省測(cè)繪工程院，江蘇 南京 210013)

地勢(shì)起伏度；地形因子；最佳分析窗口；白龍江流域

基于SRTM-DEM數(shù)據(jù)，以青藏高原東緣白龍江流域33個(gè)子流域?yàn)槔肁rcGIS的鄰域分析及均值變點(diǎn)分析方法確定了白龍江流域地勢(shì)起伏度提取的最佳分析窗口面積，并分析了子流域地勢(shì)起伏度與坡度、坡度變率、高差等微觀地形因子之間的關(guān)系，結(jié)果表明：白龍江流域地勢(shì)起伏度的最佳分析窗口面積為2.340 9 km2；白龍江不同子流域的地勢(shì)起伏度可依據(jù)流域高差、平均坡度、平均坡度變率等微觀地形因子建立相應(yīng)的計(jì)算模型。

不同的地形因子可以從不同側(cè)面反映地形特征。根據(jù)其所描述的空間區(qū)域范圍，可以將常用的地形因子劃分為微觀地形因子與宏觀地形因子兩種基本類型[1]。微觀地形因子反映地形局部特征，宏觀地形因子反映一定區(qū)域上地形的總體特征。在研究宏觀地形因子中需要首先界定其最佳分析面積。關(guān)于宏觀地形因子的尺度效應(yīng)研究較多的是地勢(shì)起伏度。國(guó)際地理學(xué)會(huì)地貌調(diào)查制圖委員會(huì)采用16 km2作為歐洲地貌圖的統(tǒng)計(jì)面積；涂漢明等[2]指出中國(guó)地勢(shì)起伏度的最佳分析面積為21 km2；郎玲玲等[3]計(jì)算得出4.41 km2與0.40 km2分別為基于1∶25萬(wàn)和1∶10萬(wàn)兩種DEM數(shù)據(jù)的福建地區(qū)地勢(shì)起伏度最佳統(tǒng)計(jì)窗口面積；程維明等[4]指出我國(guó)存在0.4、4、12、18、21 km2五種不同規(guī)模的地勢(shì)起伏度最佳單元，并分別對(duì)應(yīng)著不同的比例尺；于慧等[5]利用比例尺相當(dāng)于1∶10萬(wàn)的ASTER GDEM確定三峽庫(kù)區(qū)地勢(shì)起伏度的最佳統(tǒng)計(jì)單元為0.15 km2；張偉等[6]的研究結(jié)果顯示，全國(guó)1∶25萬(wàn)和1∶10萬(wàn)DEM數(shù)據(jù)的最佳統(tǒng)計(jì)窗口面積分別為4.72和3.20 km2；張錦明等[7]在中國(guó)區(qū)域內(nèi)隨機(jī)選取78個(gè)矩形試驗(yàn)區(qū)域的DEM數(shù)據(jù)作為試驗(yàn)對(duì)象，建立了地形起伏度最佳分析區(qū)域預(yù)測(cè)模型。隨著流域在地貌學(xué)中的應(yīng)用研究越來(lái)越深入，以流域?yàn)閱挝婚_(kāi)展地貌學(xué)研究已成為一種趨勢(shì)[8-9]。而以流域?yàn)檠芯繉?duì)象時(shí)，不同的研究區(qū)域提取宏觀地形因子地勢(shì)起伏度的最佳分析窗口面積是否一致，地勢(shì)起伏度與微觀地形因子之間有何種關(guān)系尚缺乏研究。本研究以青藏高原東緣白龍江流域32個(gè)子流域?yàn)槔?，嘗試對(duì)上述問(wèn)題進(jìn)行探討。

1 研究區(qū)概況

白龍江是嘉陵江上游最大的一級(jí)支流，發(fā)源于甘川交界的岷迭山系西段郎木寺以西的郭爾莽梁北麓，流經(jīng)甘川高原、岷迭山系和西秦嶺，主流向由西北至東南，穿越龍門(mén)山后在四川昭華匯入嘉陵江。干流全長(zhǎng)576 km，流域面積31 808 km2，天然落差達(dá)2 783 m。地勢(shì)西北高東南低，水系不對(duì)稱發(fā)育，一級(jí)支流有白水江、達(dá)拉溝、臘子溝、拱壩河等，河谷兩岸多分布夷平面及河流階地[10]。流域處于青藏高原與秦嶺造山帶的過(guò)渡地帶，晚新生代以來(lái)，伴隨著青藏高原的不斷隆升與擴(kuò)展，構(gòu)造活動(dòng)強(qiáng)烈。利用ArcMap的水文分析模塊提取了白龍江流域的33個(gè)子流域(均為白龍江的一級(jí)支流)，其中干流右側(cè)14個(gè)、左側(cè)19個(gè)(圖1)。

圖1 白龍江流域子流域劃分

2 材料與方法

2.1 數(shù)據(jù)來(lái)源

本研究所使用的SRTM 90 m DEM數(shù)據(jù)下載自中國(guó)科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心國(guó)際科學(xué)數(shù)據(jù)服務(wù)平臺(tái)(http://datamirror.csdb.cn)，再經(jīng)過(guò)投影轉(zhuǎn)換、拼接、裁切、高程異常值重新填補(bǔ)而成。其平面基準(zhǔn)為WGS 84，高程基準(zhǔn)為EGM 96，垂直精度為±20 m，水平精度為±16 m。

2.2 地勢(shì)起伏度提取原則

地勢(shì)起伏度是反映地勢(shì)起伏的宏觀地貌因子，被定義為某一確定面積內(nèi)最高點(diǎn)和最低點(diǎn)之高差，用公式可表示為

R=Hmax-Hmin

(1)

式中：R為分析區(qū)域內(nèi)的地勢(shì)起伏度值；Hmax、Hmin分別為分析區(qū)域內(nèi)的最高高程值和最低高程值。

涂漢明等[2]在分析中國(guó)1∶100萬(wàn)最佳地勢(shì)起伏度計(jì)算面積時(shí)，指出地勢(shì)起伏度最佳分析區(qū)域面積的提取需要滿足山體完整性和區(qū)域普適性兩個(gè)原則。山體完整性是指分析區(qū)域大小合適地反映了山體的完整性，區(qū)域普適性是指求取出滿足最大山體完整性的分析區(qū)域。

2.3 鄰域分析法

目前基于地勢(shì)起伏度的計(jì)算絕大部分是利用ArcMap中的鄰域分析法來(lái)實(shí)現(xiàn)的。其計(jì)算過(guò)程是以目標(biāo)柵格為中心，開(kāi)辟3×3、5×5或更大分析半徑的矩形窗口，求取分析窗口內(nèi)的高差，作為目標(biāo)柵格的起伏度。以DEM上的每個(gè)柵格作為目標(biāo)柵格，用這個(gè)分析窗口對(duì)研究區(qū)逐柵格求取高差，就得到了地勢(shì)起伏度的柵格數(shù)字矩陣。同時(shí)，在計(jì)算某一分析區(qū)域尺度的地勢(shì)起伏度時(shí)需考慮邊界效應(yīng)，刪除不能構(gòu)成完整分析區(qū)域的柵格。選取了大于白龍江流域范圍的DEM進(jìn)行不同分析窗口起伏度的提取，再利用白龍江流域邊界進(jìn)行裁切。

2.4 最佳分析窗口的計(jì)算

在確定所提取地勢(shì)起伏度的最佳分析窗口面積時(shí)，常用的方法有最大高差法、模糊數(shù)學(xué)法[2]和均值變點(diǎn)法[11-13]等。本研究采用的是均值變點(diǎn)法，該方法對(duì)恰有一個(gè)變點(diǎn)的檢驗(yàn)最為有效[14]。

3 結(jié)果與討論

3.1 白龍江流域最佳地勢(shì)起伏度的確定

為了確定白龍江33個(gè)子流域地勢(shì)起伏度的最佳分析窗口大小，利用鄰域分析法提取了SRTM 90 m DEM尺度下不同分析窗口的最大地勢(shì)起伏度值，如表1所示。

表1 白龍江流域SRTM 90 m DEM提取的地勢(shì)起伏度值

地勢(shì)起伏度與分析區(qū)域邊長(zhǎng)之間關(guān)系的散點(diǎn)圖見(jiàn)圖2，采用對(duì)數(shù)函數(shù)進(jìn)行擬合，公式為

Y=665.35lnx-3 127.2 (R2=0.973 2)

(2)

式中：Y為最大地勢(shì)起伏度值，m；x為分析區(qū)域邊長(zhǎng)，m。

圖2 白龍江流域SRTM 90 m DEM數(shù)據(jù)計(jì)算得到的散點(diǎn)圖

運(yùn)用均值變點(diǎn)分析法[14]計(jì)算得出，SRTM 90 m DEM尺度下白龍江流域地勢(shì)起伏度的最佳分析窗口為17×17(表1)，對(duì)應(yīng)的分析窗口面積為2.340 9 km2。圖3為計(jì)算所得的白龍江流域地勢(shì)起伏度空間分布圖，可以看出白龍江流域地勢(shì)起伏度較大值分布在干流及部分支流的深切河谷兩岸，這些地區(qū)斷裂帶分布較多，屬于滑坡發(fā)育的密集區(qū)域，水土流失和地質(zhì)災(zāi)害嚴(yán)重，反映出新構(gòu)造運(yùn)動(dòng)強(qiáng)烈的特征。

圖3 白龍江流域地勢(shì)起伏度

3.2 地勢(shì)起伏度與微觀地形因子之間的計(jì)算模型

在提取最佳地勢(shì)起伏度時(shí)，需要逐步增大分析窗口，而微觀地形因子的計(jì)算主要是基于3×3分析窗口，并不需要增大分析窗口，那么地勢(shì)起伏度是否與微觀地形因子之間存在某種關(guān)系，從而使地勢(shì)起伏度的計(jì)算得以簡(jiǎn)化，這是本小節(jié)要探討的問(wèn)題。前人計(jì)算表明，地勢(shì)起伏度和分析區(qū)域邊長(zhǎng)之間存在對(duì)數(shù)函數(shù)關(guān)系[2-3,6-7,11-12]，即Y=alnx±b，其中a、b為擬合參數(shù)。

為了嘗試建立地勢(shì)起伏度與微觀地形因子之間的關(guān)系，采用同樣的方法擬合了白龍江流域33個(gè)子流域的最大地勢(shì)起伏度與分析區(qū)域邊長(zhǎng)之間的對(duì)數(shù)函數(shù)關(guān)系，并計(jì)算了各子流域的微觀地形因子，包括流域的最小高程、最大高程、平均高程、流域高差、平均坡度、平均坡度變率。各微觀地形因子的統(tǒng)計(jì)及擬合函數(shù)系數(shù)見(jiàn)表2(以編號(hào)為1～5的流域?yàn)槔?，擬合系數(shù)均比較高，擬合效果好。

表2 白龍江流域典型子流域微觀地形因子統(tǒng)計(jì)及擬合系數(shù)

運(yùn)用相關(guān)分析法分析參數(shù)a、b和微觀地形因子之間的相關(guān)關(guān)系(表3)發(fā)現(xiàn)，a、b和最大高程、平均高程、流域高差、平均坡度存在一定的相關(guān)性。但是由于a、b值與這些微觀地形因子的可決系數(shù)較低，單獨(dú)使用一種地形因子擬合未知參數(shù)達(dá)不到理想效果，因此運(yùn)用逐步回歸分析法[7]建立了a、b和最小高程、最大高程、平均高程、流域高差、平均坡度、平均坡度變率的多元線性回歸函數(shù)。結(jié)果發(fā)現(xiàn)，流域高差、平均坡度、平均坡度變率可以參與對(duì)a、b值的擬合，即當(dāng)置信水平為0.05時(shí)，它們對(duì)a、b值具有顯著性影響。建立的回歸方程為

a=540.466+0.077H+28.273S-297.894SOS

(R2=0.711)

(3)

b=3 545.169+0.31H+185.41S-2 079.465SOS

(R2= 0.626)

(4)

式中：S為流域平均坡度，(°)；SOS為流域平均坡度變率；H為流域高差,m。

表3 擬合參數(shù)a、b和微觀地形因子相關(guān)系數(shù)

3.3 計(jì)算模型的驗(yàn)證

為了驗(yàn)證建立的基于微觀地形因子的最佳分析區(qū)域計(jì)算模型是否正確，提取白龍江流域相關(guān)微觀地形因子，統(tǒng)計(jì)白龍江流域的流域高差、平均坡度、平均坡度變率分別為4 533 m、25.9°、3.12。將各因子值帶入式(3)、(4)，計(jì)算可得a=683.26，b=3 234.58。式(2)中，擬合所得參數(shù)a=665.35(在0.05的置信水平時(shí)置信區(qū)間為[642.35，688.35])，b=3 228.05(在0.05的置信水平時(shí)置信區(qū)間為[2 952，3 302])。兩種方法計(jì)算所得a、b值的對(duì)比驗(yàn)證結(jié)果表明，所建立的白龍江流域最佳分析窗口計(jì)算模型是正確的，白龍江流域的地勢(shì)起伏度計(jì)算模型為

Y= (540.466+0.077H+28.273S-297.894SOS)×

lnx+3 545.169+0.31H+185.41S-

2 079.465SOS

(5)

4 結(jié) 論

(1)在水平分辨率SRTM90mDEM尺度下，應(yīng)用均值變點(diǎn)分析法得出白龍江流域地勢(shì)起伏度的最佳分析窗口面積為2.340 9km2，為該流域地貌類型的空間劃分提供了依據(jù)。

(2)白龍江不同子流域的地勢(shì)起伏度與微觀地形因子之間可以建立Y=(540.466+0.077H+28.273S-297.894SOS)lnx+3 545.169+0.31H+185.41S-2 079.465SOS的計(jì)算模型。

[1] 湯國(guó)安,楊昕. ArcGIS 地理信息系統(tǒng)空間分析實(shí)驗(yàn)教程[M].北京:科學(xué)出版社,2006:1-32.

[2] 涂漢明,劉振東.中國(guó)地勢(shì)起伏度最佳統(tǒng)計(jì)單元的求證[J].湖北大學(xué)學(xué)報(bào):自然科學(xué)版,1990,12(3):266-271.

[3] 郎玲玲,程維明,朱啟疆,等.多尺度DEM提取地勢(shì)起伏度的對(duì)比分析——以福建低山丘陵區(qū)為例[J].地球信息科學(xué),2007,9(6):1-6.

[4] 程維明,周成虎,蔡慧霞,等.中國(guó)陸地地貌基本形態(tài)類型定量提取與分析[J].地球信息科學(xué)學(xué)報(bào),2009,11(6):725-736.

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

[6] 張偉,李愛(ài)農(nóng).基于DEM的中國(guó)地勢(shì)起伏度適宜計(jì)算尺度研究[J].地理與地理信息科學(xué),2012,28(4):8-12.

[7] 張錦明,游雄.地形起伏度最佳分析區(qū)域預(yù)測(cè)模型[J].遙感學(xué)報(bào),2013,17(4):728-741.

[8] 朱紅春,陳楠,劉海英,等.自1∶10 000比例尺DEM提取地形起伏度——以陜北黃土高原的實(shí)驗(yàn)為例[J].測(cè)繪科學(xué),2005,30(4):86-88.

[9] 祝士杰,湯國(guó)安,李發(fā)源,等.基于DEM的黃土高原面積高程積分研究[J].地理學(xué)報(bào),2013,68(7):921-932.

[10] 郭進(jìn)京,韓文峰,梁收運(yùn),等.西秦嶺岷縣-宕昌地區(qū)洮河和岷江階地特征對(duì)比研究——兼論中國(guó)南北構(gòu)造帶在西秦嶺的地貌響應(yīng)[J].地質(zhì)調(diào)查與研究,2006,29(4):271-278.

[11] 王玲,呂新.基于DEM的新疆地勢(shì)起伏度分析[J].測(cè)繪科學(xué),2009,34(1):113-116.

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

[13] 張學(xué)儒,官冬杰,牟鳳云,等.基于ASTER GDEM數(shù)據(jù)的青藏高原東部山區(qū)地勢(shì)起伏度分析[J].地理與地理信息科學(xué),2012,28(3):11-14.

[14] 項(xiàng)靜恬,史久恩,李笑吟.非線性系統(tǒng)中數(shù)據(jù)處理的統(tǒng)計(jì)方法[M].北京:科學(xué)出版社,2000:3-6.

(責(zé)任編輯 李楊楊)

江蘇省高校自然科學(xué)研究面上項(xiàng)目(16KJB170014)

S157

A

1000-0941(2017)03-0038-04

常直楊(1987—)，男，河南濟(jì)源市人，講師，博士，主要從事環(huán)境演變及GIS在地貌中的應(yīng)用研究。

2015-10-15