宋卓沁張軍龍*李建軍

1）中國地震局地震預測研究所，北京 100036

2）中北大學理學院，太原 030051

基于模糊聚類和空間自相關(guān)性的河流階地自動提取方法
——以岷江上游漳臘盆地為例

宋卓沁1）張軍龍1）*李建軍2）

1）中國地震局地震預測研究所，北京 100036

2）中北大學理學院，太原 030051

傳統(tǒng)研究方法以點、線代表區(qū)域地貌面，數(shù)據(jù)缺乏連續(xù)性和系統(tǒng)性，限制了對階地演化的進一步研究。高分辨率DEM處理方法的發(fā)展使廣泛精確獲取地貌面的定量信息、重建階地模型、進一步研究階地演化成為可能。以岷江上游漳臘盆地為研究區(qū)域，采用SPOT5高分辨率數(shù)據(jù)影像構(gòu)建DEM，將統(tǒng)一標度后的階地橫剖面作為基本研究單元。通過在不同維度使用模糊C均值聚類，得到研究區(qū)域內(nèi)階地總體形態(tài)特征。在提取各級階地面的過程中，將階地剖面間的空間相關(guān)性與階地總體形態(tài)特征納入分析，最終建立了一種基于區(qū)域總體特征和區(qū)域內(nèi)部自相關(guān)的階地提取方法。對不同級階地面提取的精度評價表明，其結(jié)果與野外實測數(shù)據(jù)有較高的一致性。文中以實驗結(jié)果為依據(jù)，進一步討論了漳臘盆地內(nèi)階地的形態(tài)特征、河流下切速率和區(qū)域隆升之間的關(guān)系。因此，將模糊聚類分析和空間自相關(guān)性結(jié)合的區(qū)域地貌面提取方法是簡捷易行的，適用于階地等層狀地貌面發(fā)育區(qū)域構(gòu)造-地貌等相關(guān)問題的研究。

河流階地剖面 模糊C均值聚類分析 數(shù)字高程模型 空間相關(guān)性 漳臘盆地

0 引言

河流階地是構(gòu)造地貌研究的重要內(nèi)容之一（程紹平等，2008）。因其形態(tài)、結(jié)構(gòu)特征可以反映構(gòu)造作用的強度而受到廣泛的關(guān)注，尤其是在構(gòu)造作用強烈的青藏高原周緣深切割區(qū)域（Kirby et al.，2000）。

然而，傳統(tǒng)的測量方法往往通過使用單條剖面或者一系列高程相位點來描述1個區(qū)域階地的特點。該方法獲得階地空間參數(shù)的效率較低，限制了階地研究的進一步發(fā)展。隨著GPS技術(shù)的發(fā)展，DGPS方法的應用提高了空間參數(shù)獲取的效率，該方法是一種精確測量的有效手段，在較小范圍內(nèi)，運用DGPS方法可以獲得精確的層狀地質(zhì)地貌面數(shù)據(jù)，通過數(shù)據(jù)分析，可以得到新構(gòu)造時期不同時間尺度內(nèi)構(gòu)造演化過程（張軍龍等，2008），促進了河流階地的研究。即便如此，也僅能通過多條階地剖面對比，得到區(qū)域內(nèi)階地的相對關(guān)系。由于受到野外地形地貌復雜性的限制，該方法記錄河谷盆地的地形地貌可能丟失部分有效信息，此時將無法滿足快速精細刻畫區(qū)域構(gòu)造活動的要求。

目前，高分辨率遙感技術(shù)的發(fā)展和DEM數(shù)據(jù)精度的不斷提高，促進了快速獲得河流地貌空間參數(shù)技術(shù)的發(fā)展。Straumann等（2008）提出了區(qū)域生長法，該方法以河網(wǎng)中的點作為種子生長的起始點，通過基于對象的自頂向下的算法對相鄰的像元經(jīng)過相似性判定，根據(jù)是否滿足坡度閾值進行區(qū)域的劃分，對河谷區(qū)域、河道、臺地面、山脊等進行劃分，初步實現(xiàn)在地貌面上對河谷地貌的識別。Demoulin等（2007）將河拔高度重新投影到高程-坡度圖上，得到高度集中的平坦階地面的高程值，重建了連續(xù)階地縱剖面，該方法將階地的形態(tài)特征抽象成高程和坡度這2個可以量化的值，并將聚類分析的數(shù)學思想引入到階地劃分中；宮會玲等（2008）結(jié)合Demoulin等（2007）提出的方法對安寧河斷裂紫馬跨地區(qū)的高程進行分類，并將分類得到的河流階地結(jié)果表示在區(qū)域面上；Bi等（2012）使用K-mean聚類分析，對分維模型得到的具有2維參量的單位區(qū)域進行聚類，同時考慮了地形起伏頻率和起伏度對地貌類型和地層分布的影響，并將結(jié)果直觀地表現(xiàn)在區(qū)域面上。這些研究促進了包括階地在內(nèi)的地貌面的提取和分析。

在前人研究的基礎(chǔ)上，本文基于岷江上游漳臘盆地高精DEM數(shù)據(jù)，對復雜地貌環(huán)境中的細節(jié)信息進行識別。首先，采用條件模擬方法對數(shù)據(jù)進行預處理，減少人工因素的影響；再以地形橫剖面為基本研究單元，進行模糊C均值聚類分析，該方法根據(jù)階地的高度和點的位置特征將橫剖面上的點自然劃分為不同的類別。并將地貌特征的空間相關(guān)性納入分析過程中，充分考慮了相鄰位置地貌的差異性和相似性。建立了一種河流階地自動提取方法。根據(jù)實驗結(jié)果，并結(jié)合測年數(shù)據(jù)，得到了區(qū)域形態(tài)特征和河流下切速率，并嘗試解釋區(qū)域隆升和河流下切之間的關(guān)系。

1 研究區(qū)構(gòu)造背景

研究區(qū)位于青藏高原東緣岷山隆起帶內(nèi)。岷江斷裂主斷裂為岷山隆起的西邊界斷裂，主體沿岷江展布，全長約110km，走向為N15°～20°E，傾向NW，傾角不定，為逆沖兼左旋斷層（楊景春等，1979；趙小麟等，1994；錢洪等，1995），水平滑動速率為1～2mm／a（趙小麟等，1994；唐文清等，2004），垂直滑動速率為0.37～0.53mm／a（周榮軍等，2000；張軍龍等，2013）。第四紀以來，岷江斷裂北段（貢嘎嶺—虹橋關(guān)段）構(gòu)造活動較為頻繁，活動性較強。岷江斷裂北段又以紅橋關(guān)隆起和貢嘎嶺隆起為界，可分為3亞段。各亞段的力學特征與地貌單元形態(tài)吻合（唐榮昌等，1991；趙小麟等，1994；周榮軍等，2000），其中，南亞段長約10km，走向N10°～15°E，漳臘盆地的展布主要受其控制（圖1）。

作為岷江流域上游的獨立新近紀侵蝕盆地之一，漳臘盆地沉積主體呈楔形，西厚東薄，礫徑自西向東變細（楊農(nóng)等，2003）。盆地發(fā)育3～4級不對稱階地，階地基座上覆沉積物分層明顯。沉積物上部主要為黃灰色細礫石夾砂質(zhì)黏土層，下部主要為河漫灘相礫石夾粉細砂。盆地內(nèi)岷江河道長約10km，河流縱剖面高度自北向南降低了約100m，坡比降近1%，具有較強的侵蝕力。

2 階地面提取方法

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

實驗選擇spot5異軌立體影像對為數(shù)據(jù)源。Spot5全色影像（分辨率為5m）由搭載在Spot5衛(wèi)星上的高分辨率立體成像儀（HRS）獲取，通過對相鄰軌道拍攝的立體相對進行處理，可以獲得2.5m的高分辨率DEM。通過對成像時間為2008年1月和2010年1月的spot5異軌影像對（影像中心：103.63°E，32.85°N）進行模型解算、地面校驗點驗證，生成DEM精度為5～10m（圖2）。由于低階地面容易受到建筑物等因素影響，為了降低噪聲，對于研究區(qū)域中較大范圍居民地，采用條件模擬法（conditional stimulation）（Strebelle et al.，2002）預先對DEM數(shù)據(jù)進行填補處理：利用居民地周圍未受明顯人為活動改動的高程數(shù)據(jù)作為條件，估算居民地在未受人為活動影響情況下的高程值，以此消除人為活動對自然地貌的影響。

2.2 河流階地面提取原理

假設(shè)實驗區(qū)域內(nèi)，不考慮人為因素，階地發(fā)育連續(xù)，各級階地內(nèi)部高度不發(fā)生突變，即使是地貌面起伏較大的高階地，高度變化也應是連續(xù)的。

2.2.1 模糊C均值聚類算法

實驗中以階地橫剖面不同位置的高程值（相關(guān)方法參考本文2.3.1節(jié)和2.3.2節(jié)）為研究數(shù)據(jù)，而高程值相對連續(xù)故無法得到明確的類別劃分。針對這一特點，使用模糊C均值聚類算法（FCM）（Pal et al.，1995）提取區(qū)域階地總體特征。

聚類分析是統(tǒng)計分析中一類經(jīng)典算法，目的是在缺少先驗知識的情況下，按照不同數(shù)據(jù)對象自身的特征自然地將其劃分到不同的類別中。本文中將高程值劃分為相對集中的幾類，各類的中心可以作為階地高程的代表值。模糊聚類能夠描述樣本類屬的中介性，即一個樣本可以屬于多個聚類并計算屬于各類別的概率，從而避免對高程值的硬劃分。FCM算法是一個迭代過程，通過以下步驟確定n個數(shù)據(jù)點的聚類中心ci和隸屬度矩陣U。

（1）利用0～1間的隨機數(shù)初始化隸屬度矩陣U，每個數(shù)據(jù)點屬于各個模糊組的隸屬度之和等于1：

圖1 岷江流域漳臘盆地位置（a）及地質(zhì)圖（b）（改自Kirby et al.，2000）Fig.1 The location of Zhangla Basin（a）and regional tectonic geologicalmap（b）（after Kirby et al.，2000）.

圖2 利用spot5異軌立體像對生成岷江上游DEM（斷層位置根據(jù)張軍龍等，2013）Fig.2 The DEM of upper Minjiang River created by SPOT5 across-track stereo images.

（2）計算c個聚類中心ci：

（3）計算價值函數(shù)

Uij∈（0，1）之間，ci為模糊組i的聚類中心，dij=‖ci-xj‖為第i個聚類中心與第j個數(shù)據(jù)點間的歐幾里得距離，m屬于［1，∞］是一個加權(quán)指數(shù)。如果它小于某個確定的閾值或者相對上一次價值函數(shù)的改變量小于某個閾值，則算法停止，使上式達到最小值的必要條件為λj（j=1，…，n）是式（1）的n個約束式的拉格朗日乘子。

否則，d使用

計算新的U矩陣，返回步驟2。

由于算法的性能依賴于初始聚類中心及其初始類別數(shù)，而初始聚類中心的確定卻是隨機的，所以不能確保FCM收斂于一個全局最優(yōu)解。但是由于聚類中心是多次迭代的結(jié)果，雖然有少數(shù)結(jié)果僅收斂于局部最優(yōu)解，但大部分迭代的結(jié)果均收斂于同一聚類中心。實驗中通過多次運行FCM，仍能確定全局最優(yōu)解，對劃分階地的級別影響不大。

上述聚類之后得到了漳臘盆地各級階地特征高程值和階地面分布比例（特征階地剖面），該結(jié)果為全局經(jīng)驗最優(yōu)值，是一個統(tǒng)計量，它考慮了區(qū)域內(nèi)階地面特征的統(tǒng)計結(jié)果，反映了整個區(qū)域的階地總體規(guī)律。

2.2.2 加權(quán)平方差和

但是特征階地剖面并非任意一條階地橫剖面的實際情況，無法通過特征剖面線和局部剖面線之間的距離來丈量相關(guān)性。因此，本文采用類似于監(jiān)督分類的思路，通過選擇與特征剖面較為接近的樣本剖面，并引入專家知識，微調(diào)得到符合實際情況的樣本剖面線。

而樣本剖面線與特征階地剖面之間的相似程度可以通過各點的加權(quán)平方差和來度量，加權(quán)平方差和的公式為

式（6）中：i為階地級數(shù)，n為階地級數(shù)的總和，m為每一級標準階地包含的點數(shù)量，λ為經(jīng)驗權(quán)值。由于高階地形成時間較早，侵蝕時間更長，地貌面起伏較大，所以，其高度的變化將會主導加權(quán)平方差和的值。由于不同級階地面提取的準確性對結(jié)果準確性的貢獻是一致的，為了平衡高階地和低階地與標準剖面的偏差，考慮到階地的拔河高度與年齡呈正相關(guān)，標準剖面的偏差值隨階地高度增加而變大，所以將各級階地的高度占總拔河高度的比例的倒數(shù)作為經(jīng)驗權(quán)值λ，以平衡高階地和低階地之間偏差量的不同。

2.2.3 空間自相關(guān)的理論基礎(chǔ)

地理事物在空間分布上存在著一定規(guī)律，地表所有事物和現(xiàn)象在空間上都是關(guān)聯(lián)的（Tobler，1970）。距離越近，關(guān)聯(lián)程度就越強，距離越遠，關(guān)聯(lián)程度就越弱。階地發(fā)育作為一種空間現(xiàn)象，在平面內(nèi)的變化具有連續(xù)性，即在不考慮人為改造地貌因素的前提下，每一條局部階地剖面都應和相鄰位置局部階地剖面具有高度的相關(guān)性，而隨著局部剖面之間距離的增加，二者之間的相關(guān)性將明顯下降。上文得到的幾條樣本剖面，既具有了區(qū)域內(nèi)特征階地剖面的共性特點又具有樣本所在位置的局部特征，將作為判斷這種地理相關(guān)性的起始位置。

將相鄰未人工劃分的局部剖面與已劃分的剖面進行比較，根據(jù)二者的相關(guān)性及通過坡度變率、高度、起伏度等局部地貌特征（Hengl et al.，2008）劃定其余局部剖面。而設(shè)置多條樣本剖面線避免了直接由單一特征剖面線確定大量局部剖面線帶來的誤差。

2.3 實驗步驟

以岷江上游漳臘盆地為研究區(qū)域，利用該地區(qū)高分辨率DEM（圖2）提取橫剖面并對其進行處理。

2.3.1 提取階地橫向剖面

階地橫剖面角度的選取與河道流向相關(guān)，考慮階地剖面方向近垂直于河道方向，但自然河流線形狀復雜，大量河曲導致剖面線平行河流甚至反向。為防止該情況發(fā)生，在剖面提取之前對岷江河道進行裁彎取直，得到較為平直的河流線。階地橫剖面長度選取則以地質(zhì)圖中第四紀沉積物的分布范圍為剖面線的終點。最終，分別在東西兩岸自動提取垂直于河流線、等間距分布的245條局部階地橫剖面，每條剖面均由間隔10m的一系列高程點組成（圖3）。

圖3 垂直裁彎取直后的岷江河道，等間距提取階地橫剖面；對剖面進行重采樣并投影在長度-河拔圖上Fig.3 After cutting off the Minjiang River，fluvial terraces are perpendicularly obtained with cross sections evenly spaced.Each cross section is resampled and projected onto the length-height of river sketch.

2.3.2 橫剖面數(shù)據(jù)處理

在上述（2.2）平穩(wěn)變換的假設(shè)下進一步分析：以河拔高度0為起始位置，將245條階地剖面的長度重新采樣（張維等，2011），以點到河的垂直距離與剖面長度的比值代替固定距離作為兩點之間的間隔，形成以長度比例為橫軸、河拔高度為縱軸的階地剖面線（圖3）；對組成剖面線的高程點進行篩選，按照前人經(jīng)驗（Demoulin et al.，2007），保留角度＜13°的點，組成平緩的階地面序列；將這些點投影在同一長度-拔河圖上進行進一步分析（圖3）。

2.3.3 獲得特征階地剖面

為排除離散點群，對圖3中每一縱列點的高程值取置信度為80%時的平均值作為該位置的階地高度的指示量。為了區(qū)分各級階地，對所有位置的高度均值統(tǒng)計值進行聚類，在無先驗知識情況下，將高度均值分成6類，再根據(jù)先驗知識分別對這6類結(jié)果進一步合并或劃分。高階地內(nèi)河流及其分支下切造成的沖溝和風化侵蝕往往對聚類的數(shù)值計算產(chǎn)生較大的影響，因此需要對較高幾類的聚類結(jié)果進行合并；而低階地在類間差距上小于高階地，因而往往將低階地聚為1類，所以在第1個聚類內(nèi)部再次進行劃分。漳臘盆地內(nèi)發(fā)育3級低階地，故將聚類中心數(shù)設(shè)為3。經(jīng)過上述實驗步驟，在左右兩岸均形成高度和位置上集中的4個類別。聚類中心為該級階地面的高度，每個類別的起始值對應階地陡坎的位置，每一類中的最大值和最小值則是相應階地高度的范圍（圖4）。

圖4 a漳臘盆地岷江西岸2次聚類后結(jié)果圖；b漳臘盆地岷江東岸2次聚類后結(jié)果；c FCM聚類算法提取河流剖面（黑色）與野外實測剖面（灰色）的對比Fig.4 Results of twice clustering on the west（a）and east（b）bank of M injiang River，and the comparison（c）between the profile generated by FCM clustering（solid line）and field survey（dotted line）.

表1 東岸第160，20，80條樣本剖面的加權(quán)平方差Table 1 WSSD of sections No.20，80 and 160 on the east side of Minjiang River

2.3.4 選取樣本剖面

根據(jù)特征剖面中各級階地的垂直高度和水平分布，衡量每條局部剖面與特征階地剖面的相似性（式（6），表1），在兩岸各選擇若干條與特征剖面差異較小的局部剖面作為待選樣本剖面。由于剖面分布的連續(xù)性，直接選擇差異最小的局部橫剖面極有可能是連續(xù)分布的（第159與第160、第161條階地），不利于誤差的控制。故再在這些候選剖面中選擇3條離散分布的局部剖面，經(jīng)過人工糾正階地剖面產(chǎn)生的錯分、漏分，確定每一級階地的邊界，形成3條具有特征剖面特征和該位置特點的樣本階地剖面（圖5）。對岷江左岸也做相應處理。

圖5 糾正后得到的樣本剖面Fig.5 The sample sections after correcting errors.

假設(shè)圖6中的a，b，c，d為4條連續(xù)排列的階地橫剖面，所有點在水平方向上均勻分布，劃分步驟如下：

（a）第1條為樣本剖面，根據(jù)特征剖面和專家知識劃分各級階地。統(tǒng)計每一級階地的最高點和最低點，確定該階地在垂直方向的范圍，而起點和終點則確定了階地在水平方向上的范圍（灰色區(qū)域）。

（b）第2條為與樣本剖面相鄰的局部剖面，虛線框為樣本剖面a的各級階地范圍，如果在b剖面的點同樣屬于這個區(qū)域，則認為這些點的歸屬和上一剖面相同。對于不在這一區(qū)域的點，若點的數(shù)目≤2，則計算坡度變率，如果坡度變率＜1，則屬于高一級的階地，反之則屬于較低一級。

（c）第3條局部剖面，使用b中所用方法判斷每個點的階地歸屬，通過確定每一級階地的最高點和最低點，以確定階地的范圍。

（d）第4條局部剖面，如果虛線區(qū)域外有多個點，則計算最后一個點的坡度變率，判斷階地歸屬。由于實驗中，較低階地更容易受到干擾，故在判斷歸屬時，如屬于高一級的階地，則再判斷前一個點的階地歸屬，以此類推；而如點屬于低一級階地，則其之前的點均屬于低一級階地。

2.3.5 劃分各局部剖面

從3條樣本剖面起算，根據(jù)相鄰位置階地剖面的相似性及各局部剖面之間的差異性：對于特征相似的部分，直接套用樣本剖面的分割位置對局部剖面進行劃分；而對于形狀差異的部分，則根據(jù)其樣本剖面的分割位置、波動程度確定搜索半徑大小，在這一范圍內(nèi)，通過坡度變率、高度、起伏度等地貌參數(shù)判斷不同級階地之間的分隔點，并進一步統(tǒng)計該位置的剖面特征，用來輔助劃分相鄰階地剖面（圖6）。

分別記錄每條橫剖面線內(nèi)各級階地的始末兩點，在沿河流方向依次連接后得到區(qū)域內(nèi)各級階地面的分布（圖7）。

3 實驗結(jié)果精度評價

通過對漳臘盆地245條階地橫剖面進行分析處理，形成漳臘盆地的特征剖面（圖4c），黑色實線為聚類分析結(jié)果，灰色虛線為野外實測的階地剖面，兩者具有高度的相似性，且每一級實測階地高程都在相應統(tǒng)計得到的高程范圍內(nèi)。因此可以得出結(jié)論，該統(tǒng)計結(jié)果能夠代表岷江流域漳臘盆地內(nèi)階地的總體特征。

通過野外實測，確定了漳臘盆地內(nèi)部分區(qū)域的階地分布（103°39′4.18″E，32°49′5.96″N），本文通過實測結(jié)果，對計算機進行提取結(jié)果驗證和精度評價（圖8）。

驗證結(jié)果中，不同級階地區(qū)分的誤差主要集中于Ⅰ級階地和右岸Ⅲ級階地，除實驗設(shè)計中的局限性外，影響實驗結(jié)果的因素主要有：Ⅰ級階地自身較低且人工建筑較多，以致階地高度、坡度都發(fā)生了較大的、非自然狀況下的變化；而右岸Ⅲ級階地分布較窄，很多區(qū)域僅有5～10m，且由于數(shù)據(jù)源的分辨率和精度的原因，極易與相鄰階地陡坎混淆，所以階地劃分結(jié)果不理想。但總體來說，各級階地的空間分布均與野外實測結(jié)果較為吻合。

為了進一步對計算機自動提取的精度進行評價，根據(jù)混淆矩陣法（Stehman et al.，1997）在研究區(qū)域內(nèi)隨機取300個校驗點，判斷在校驗點上階地分級的一致性。得到校驗點一致性的混淆矩陣（表2），計算得出，各級階地的總體分類精度為82.67%，kappa系數(shù)為0.765，說明實驗結(jié)果和真實值具有高度的一致性（kappa系數(shù)為0.61～0.80）（Landis et al.，1977）。

圖6 逐條階地局部剖面劃分示意Fig.6 The process of segmentation on every single profile.

4 存在的問題與討論

實驗中將第四紀沉積物邊界作為研究區(qū)域范圍，在其他研究區(qū)域內(nèi)可根據(jù)研究目的、方法的不同，設(shè)置不同的邊界。

圖7 漳臘盆地各級地貌面提取結(jié)果Fig.7 The result of geomorphic surfaces extraction in the Zhangla Basin.1～3為階地面，4為Q2p高階地

圖8 小范圍內(nèi)區(qū)分各級階地的階地面提取結(jié)果與實際結(jié)果的對比Fig.8 Contrast of terrace surfaces extraction and the real situation in a small range.圖中綠色區(qū)域為實際結(jié)果和計算機提取結(jié)果一致的部分，紅色區(qū)域為兩者不相符的部分

表2 樣本點（300個）一致性的混淆矩陣Table 2 The confusion matrix to illustrate the consistency of the 300 sample points

模糊C均值聚類算法作為一個廣泛應用的統(tǒng)計分析方法，使得聚類結(jié)果中類內(nèi)（同級階地）差異最小而類間（不同級階地）差異最大，實現(xiàn)了統(tǒng)計意義上的階地劃分。統(tǒng)計值概括了樣本中最主要的特點，能夠代表相應區(qū)域內(nèi)階地的特征。但是自然現(xiàn)象是一種復雜的現(xiàn)象，統(tǒng)計方法僅從數(shù)學角度對其進行了抽象，卻無法考慮空間特征，如階地分布的連續(xù)性。直接使用統(tǒng)計分析結(jié)果，極有可能與真實情況存在偏差。因此，本文在提取階地面的過程中，有限地考慮了階地分布的空間連續(xù)性。

假設(shè)漳臘盆地內(nèi)部階地高度總體穩(wěn)定，分布連續(xù)。基于空間連續(xù)性的階地面提取，從方法上，僅能概括地形連續(xù)變化的這一地學要素，尚需將其他地學要素及空間特征納入分析。但是從總體思路上，該方法將現(xiàn)實中自然要素的地貌特點納入到階地的劃分，改變了聚類算法僅從數(shù)學角度對地貌進行抽象的方法，具有一定的應用意義。

理論上，剖面截取的間隔越小，地理上的關(guān)聯(lián)程度就越強，實驗結(jié)果也將更加連續(xù)。但若進一步減小橫剖面的采樣間隔，必然導致計算量的大量增加。受限于計算效率以及說明實驗效果的目的，本文采樣了40m這一適中的采樣尺度。

漳臘盆地各級階地的平均統(tǒng)計河拔高度分別為：岷江西岸，T1（5.14m）、T2（29.06m）、T3（70.42m）、T4（204m），T1、T2較窄，一般僅占階地總寬度的20%；岷江東岸，T1（4.84m）、T2（24.28m）、T3（61.67m）、T4（194m），其中T1、T3分布較窄。各級階地均連續(xù)不對稱分布。根據(jù)統(tǒng)計結(jié)果，岷江西岸各級階地均明顯高于東岸，兩者間的差值總體上隨階地的增高而增大。形成這樣的形態(tài)特征主要是由于受到岷江斷裂由西向東的逆沖推覆作用（張岳橋等，2012），形成整體西高東低的形態(tài)特征，是岷江河谷演化的外在表現(xiàn)。

岷江在漳臘盆地段的下切速率可以通過河拔高度和階地年齡估算，其中漳臘盆地內(nèi)各級階地面的拔河高度通過兩岸特征剖面的統(tǒng)計高度估算，根據(jù)前人在該區(qū)域中測年統(tǒng)計結(jié)果（張軍龍，2013），選擇了階地二元結(jié)構(gòu)上層河漫灘堆積中接近階地面但沒有后期改造的位置的年齡數(shù)據(jù)作為階地年齡（圖9），計算得出漳臘Ⅰ級階地河流下切速率為0.64～1.78m／ka；Ⅱ級階地為1.13～2.61m／ka；Ⅲ級為2.19～2.79m／ka；Ⅳ級為0.78～1.27m／ka，晚更新世以來，岷山地區(qū)持續(xù)隆升。上述階地下切速率是青藏高原東緣構(gòu)造和氣候共同作用的結(jié)果，兩種作用的比例仍需進一步研究。

漳臘盆地東岸階地發(fā)生明顯傾斜，主要是由于岷山隆起東邊界的虎牙斷裂晚新生代以來垂向活動量相對西邊界岷江斷裂較大（周榮軍等，2000），產(chǎn)生了整體向西掀斜的結(jié)果。

圖9 河流階地測年數(shù)據(jù)所在拔河高度（［1］趙小麟等，1994；［2］Kirby et al.，2000；［3］周榮軍等，2000；［4］張軍龍，2013）Fig.9 Dating age data of river terraces at different water levels（after［1］ZHAO Xiao-lin et al.，1994；［2］Kirby et al.，2000；［3］ZHOU Rong-jun et al.，2000；［4］ZHANG Jun-long，2013）.

5 結(jié)論

本文以岷江上游漳臘地區(qū)階地為例，提出了基于模糊聚類和空間自相關(guān)性的河流階地提取方法，并結(jié)合年齡數(shù)據(jù)討論了河流下切速率和區(qū)域隆升之間的關(guān)系，得到了以下結(jié)論：

（1）通過對高分辨率DEM進行分析處理，以階地剖面線為研究單元，將地貌面的空間相關(guān)性納入數(shù)學統(tǒng)計分析，嘗試通過計算機自動提取自然階地面，并運用該方法對岷江流域漳臘盆地內(nèi)階地進行提取，得到了各級階地面的垂直和水平分布，后續(xù)的精度評價證明該方法得到的階地面分布與野外實測結(jié)果具有較高的一致性?；谀：垲惡涂臻g相關(guān)性的河流階地自動提取方法可以應用于其他階地發(fā)育的區(qū)域，為快速、高效、定量地描述和研究階地及其他面狀地貌的形態(tài)學特征提供了新的思路。

（2）岷江流域漳臘盆地內(nèi)主要發(fā)育4級階地，連續(xù)分布于河流兩側(cè)，西岸各級階地高程高于東岸，分布范圍小于東岸，兩岸階地不對稱發(fā)育，計算機提取的地貌形態(tài)特征與野外地質(zhì)調(diào)查結(jié)果相似。岷江西岸展布了西傾逆沖-推覆斷層（張岳橋等，2012），受青藏高原向東運動影響，第四紀以來西部物質(zhì)沿斷層向東仰沖，造成上盤階地抬升。

致謝 中國科學院遙感與數(shù)字地球研究所研究生李龍翔在論文撰寫中提供了幫助，審稿專家對本文提出了寶貴意見，在此一并致謝。

程紹平，楊桂枝.2008.國外新構(gòu)造研究進展評述［J］.地震地質(zhì)，30（1）：31—43.

CHENG Shao-ping，YANG Gui-zhi.2008.Current advance of overseas research on neotectonics：A review and comments［J］.Seismology and Geology，30（1）：31—43（in Chinese）.

宮會玲，冉勇康，陳立春.2008.基于DEM的階地分析方法：以安寧河斷裂紫馬跨地區(qū)為例［J］.地震地質(zhì)，30（1）：339—348.

GONG Hui-ling，RAN Yong-kang，CHEN Li-chun.2008.Themethod of terrace analysis based on DEM：A case study in Zimakua of Anninghe Fault［J］.Seismology and Geology，30（1）：339—348（in Chinese）.

錢洪，馬聲浩，龔宇.1995.關(guān)于岷江斷裂若干問題的討論［J］.中國地震，2（2）：140—146.

QIAN Hong，MA Sheng-hao，GONG Yu.1995.Discussions on the Minjiang Fault［J］.Earthquake Research in China，2（2）：140—146（in Chinese）.

唐榮昌，文德華，黃祖智，等.1991.松潘-龍門山地區(qū)主要活動斷裂帶第四紀活動特征［J］.中國地震，7（3）：64—71.

TANG Rong-chang，WEN De-hua，HUANG Zu-zhi，et al.1991.The Quaternary activity characteristics of several major active faults in the Songpan-Longmenshan region［J］.Earthquake Research in China，7（3）：64—71（in Chinese）.

唐文清，劉宇平，陳智梁，等.2004.龍門山斷裂構(gòu)造帶GPS研究［J］.大地測量與地球動力學，24（3）：57—59.

TAND Wei-qing，LIU YU-ping，CHEN Zhi-liang，et al.2004.GPS study on Longmenshan fault zone［J］.Journal of Geodesy and Geodynam ics，24（3）：57—59.

楊景春，鄧天崗，王元海，等.1979.岷江上游地區(qū)第四紀構(gòu)造應力狀態(tài)及其與地震的關(guān)系［J］.地震地質(zhì)，1（3）：68—75.

YANG Jing-chun，DENG Tian-gang，WANG Yuan-hai，et al.1979.The Quaternary tectonic stress over the up-stream area of Minjiang River in Sichuan and its relations to earthquakes［J］.Seismology and Geology，1（3）：68—75（in Chinese）.

楊農(nóng)，張岳橋，孟輝，等.2003.川西高原岷江上游河流階地初步研究［J］.地質(zhì)力學學報，9（4）：363—370.

YANG Nong，ZHANG Yue-qiao，MENG Hui，et al.2003.Study of the Minjiang River terraces in the western Sichuan plateau［J］.Journal of Geomechanics，9（4）：363—370（in Chinese）.

趙小麟，鄧起東，陳社發(fā).1994.岷山隆起的構(gòu)造地貌學探究［J］.地震地質(zhì)，16（4）：429—439.

ZHAO Xiao-lin，DENG Qi-dong，CHEN She-fa.1994.Tectonic geomorphology of the M inshan uplift in western Sichuan，southwestern China［J］.Seismology and Geology，16（4）：429—439（in Chinese）.

張軍龍，任金衛(wèi)，陳長云，等.2013.網(wǎng)格DEM剖面方法獲取岷江上游層狀地貌參數(shù)［J］.武漢大學學報（信息科學版），38：102—108.

ZHANG Jun-long，REN Jin-wei，CHEN Chang-yun，et al.2013.Acquisition of parameters of stratified reliefs by grid DEM profilesmethod on the upper Minjiang［J］.Journal of Wuhan University（Natural Science Edition），38：102—108（in Chinese）.

張軍龍，田勤儉，張小龍，等.2008.DGPS方法在新構(gòu)造研究中的應用探討［J］.地學前緣，15（4）：290—297.

ZHANG Jun-long，TIAN Qin-jian，ZHANG Xiao-long，et al.2008.The applications of DGPSmethod in neo-tectonic research［J］.Earth Science Frontiers，15（4）：290—297（in Chinese）.

張維.2011.基于DEM的陜北黃土高原流域剖面譜研究［D］：［學位論文］.南京：南京師范大學.

ZHANG Wei.2011.Research of catchment profile spectrum on northern Shaanxi loess plateau，China［D］. Dissertation.Nanjing Normal University，Nanjing.

張岳橋，李海龍，吳滿璐，等.2012.岷江斷裂帶晚新生代逆沖推覆構(gòu)造：來自鉆孔的證據(jù)［J］.地質(zhì)評論，58（2）：215—223.

ZHANG Yue-qiao，LI Hai-long，WU Man-lu，et al.2012.Late Cenozoic thrust and nappe structure along the Minjiang：Evidence from a drill hole［J］.Geological Review，58（2）：215—223（in Chinese）.

周榮軍，蒲曉虹，何玉林，等.2000.四川岷江斷裂帶北段的新活動、岷山斷塊的隆起及其與地震活動的關(guān)系［J］.地震地質(zhì)，22（3）：285—294.

ZHOU Rong-jun，PU Xiao-hong，HE Yu-lin，et al.2000.Recent activity of M injiang Fault zone，uplift of Minshan block and their relationship with seismicity of Sichuan［J］.Seismology and Geology，22（3）：285—294（in Chinese）.

Bi Lisi，He Honglin，Wei Zhanyu，et al.2012.Fractal properties of landforms in the Ordos block and surrounding areas，China［J］.Geomorphology，175-176：152—162.

Demoulin A，Bovy B，Rixhon G，et al.2007.An automated method to extract fluvial terraces from digital elevation model：The Vesdre valley，a case study in eastern Belgium［J］.Geomorphology，91：51—64.

Hengl T，Reuter H I.2008.Geomorphometry：Concepts，Software，Applications（Vol 33：Developments in Soil Science）［M］.Elsevier Science，772.

Kirby E，Whipple K X，Burchfiel B C，et al.2000.Neotectonics of the Min Shan，China：Implications formechanisms driving Quaternary deformation along the eastern margin of the Tibetan plateau［J］.GSA Bulletin，112（3）：375—393.

Landis J R，Koch G G.1977.An application of hierarchical kappa-type statistics in the assessment of majority agreement amongmultiple observers［J］.Biometrics，33（2）：363—374.

Li Yong，Densmore A L，Zhou R J，et al.2006.Profiles of digital elevation models（DEM）crossing the eastern margin of the Tibetan plateau and their constraints on dissection depths and incision rates of late Cenozoic rivers［J］. Quaternary Sciences，79（1）：28—37.

Pal N R，Bezdek JC.1995.On cluster validity for the fuzzy C-means model［J］.IEEE Trans Fuzzy Syst，3（3）：370—379.

Stehman S V.1997.Selecting and interpreting measures of thematic classification accuracy［J］.Remote SensingEnviron，62：77—89.

Straumann R K，Purves R S.2008.Delineation of valleys and valley floors［J］.Geographic Information Science，5266：320—336.

Strebelle S.2002.Conditional simulation of complex geological structures using multip le-point statistics［J］. Mathematical Geology，34（1）：1—21.

Tobler W R.1970.A computer movie simulating urban growth in the Detroit region［J］.Economic Geography，46：234—240.

A CLUSTERING APPROACH FOR INCORPORATING SPATIAL DEPENDENCE INTO THE AUTOMATIC FLUVIAL TERRACE EXTRACTION FROM DIG ITAL ELEVATION MODEL：A CASE STUDY FROM THE ZHANGLA BASIN ALONG THE UPSTREAM OF THE M INJIANG RIVER

SONG Zhuo-qin1）ZHANG Jun-long1）LI Jian-jun2）
1）Institute of Earthquake Science，China Earthquake Administration，Beijing 100036，China
2）School of Science，North University of China，Taiyuan 030051，China

The terraces along the upper reaches of the Minjiang River record rapid uplift around the Tibetan plateau since the Quaternary.However，the common method to obtain the terraces elevation data always relies on single point or line to represent the whole landform.The available results and further analyses are usually not continuous and systematic.A ll of these，therefore，restrict the development of further study.A large number of qualitative information and the digital terracemodels are relied on the improvement of DEM processing technology.Our present paper applies fuzzy C-means algorithm to a bunch of the Minjiang River terrace cross sections in different dimensions to get the integrated distribution characteristics of the terraces in the Zhangla Basin.SPOT5 image is adopted when generating the DEM.Then we consider both the spatial correlation and distribution characteristics of terrace cross sections into analysis and we consequently find a successfulway to extract different levels of terraces based on both whole indicator and internal correlation.The precision evaluation suggests that result of terrace extraction is highly consistent with the field survey data.Based on these，we discuss relationship between characteristics of terraces，the incision rates and the regional uplift patterns.We suggest that a clustering approach for incorporating spatial dependence into the automatic fluvial terrace extraction can be used to study and understand the regional tectonic and geomorphic features in the layered landscapes.

fluvial terrace，fuzzy C-means analysis，DEM，spatial correlation，Zhangla Basin

P315.2

A

0253-4967（2014）04-1029-14

宋卓沁，女，1988年生，現(xiàn)為中國地震局地震預測研究所在讀碩士研究生，研究方向為構(gòu)造地貌，電話：010-88015246，E-mail：176578721＠qq.com。

10.3969／j.issn.0253-4967.2014.04.008

2013-08-30收稿，2014-04-01改回。

國家自然科學基金（41372215）、中國地震局地震行業(yè)科研專項（201408023）和山西省基金資助項目（2012012003）共同資助。

*通訊作者：張軍龍，副研究員，E-mail：zhjulo_2002＠163.com。