災害風險地圖中黃土高原地貌自動綜合研究

潘東華，王靜愛，賈慧聰

(1.民政部國家減災中心，北京 100124； 2.民政部衛星減災應用中心, 北京 100124； 3.北京師范大學地理學與遙感科學學院, 北京 100875; 4.北京師范大學地表過程與資源生態國家重點實驗室，北京 100875；5.中國科學院遙感與數字地球研究所, 北京 100094； 6.中國科學院數字地球重點實驗室，北京 100094)

?

災害風險地圖中黃土高原地貌自動綜合研究

潘東華1,2,3，王靜愛3,4，賈慧聰5,6

(1.民政部國家減災中心，北京 100124； 2.民政部衛星減災應用中心, 北京 100124； 3.北京師范大學地理學與遙感科學學院, 北京 100875; 4.北京師范大學地表過程與資源生態國家重點實驗室，北京 100875；5.中國科學院遙感與數字地球研究所, 北京 100094； 6.中國科學院數字地球重點實驗室，北京 100094)

探討了結合圖層約束理論和小波分析的黃土高原地區DEM孕災環境的自動綜合技術，構建了基于不同尺度的水系、坡度等約束圖層信息，并通過小波方法重構了不同尺度、多圖層約束下的黃土高原地貌,為不同區域尺度下孕災環境的最優表達探索了一種新的數據處理模型和表示方法。

制圖綜合；圖層約束；黃土高原；小波分析；自然災害風險地圖

一、引 言

黃土高原是世界上水土流失最嚴重的地區，是我國水土保持的重點區域，也是防治自然災害和次生災害的重點區域；同時，本地區氣候資源變幅較寬，年際變化差異大，氣候災害發生頻率高[1]。因此，開展本地區的自然災害風險防范具有重要的現實意義。自然災害風險地圖作為一項非工程性措施，是區域災害風險防范與決策管理的重要參考依據[2]。由于不同區域尺度的地理規律差異明顯，不同區域層次應用上的要求也不相同，因此不同尺度的地貌綜合成為黃土高原研究的一種必然需求。本文在當前災害風險制圖綜合主要技術方法的基礎上[3-6]，結合小波分析方法，開展黃土高原地貌在災害風險地圖中的自動綜合研究，對全方位認識該地區的空間結構特征具有重要指導意義，同時也為特殊地貌孕災環境提供了新的數據挖掘途徑。因此，多尺度的地貌綜合技術可為區域災害風險管理提供決策支持。

二、面狀自然災害風險地圖自動綜合的理論與方法

1.理論框架

本研究基于災害系統理論，并結合地理學和地圖學知識，探討了面狀孕災環境(黃土高原地貌)在自然災害風險地圖中的自動綜合問題，在揭示孕災環境對災害風險表達明辨性、邏輯一致性和結構整體性的前提下[7]，考慮與黃土孕災環境相關的自然要素，提取了該地區多尺度水系、坡度等地貌自動綜合的主要約束圖層(LC)，并結合小波分析，通過二維小波變換，作了不同尺度下的分解與重構，從而優化了多尺度、多圖層約束下的黃土高原地區DEM的表達，為黃土高原地區地貌災害與孕災環境的尺度約束提供了一定的參考，對地表形態的多尺度表達和自動綜合及其綜合程度評估方法作了初步的研究。

2.小波變換在地貌自動綜合中的優勢

小波理論是當前時空—頻率分析工具，由于其具有自適應性和“數學顯微鏡”性質而成為許多學科研究的工具。黃土高原DEM自動綜合的小波變換(WT)模型著重回答了在地圖綜合過程中“怎么做”的問題，即當生成小比例尺地圖或低分辨率數據庫時為不同層次決策提供適量信息而作出的方法選擇。WT模型具有時頻域的良好局部化特性，隨著信號不同頻率成分在時/空域取樣的疏密自動調節，可以在任意尺度觀察函數(信號、圖像等)的任意細節并加以分析。本研究將黃土地貌看成一個信號場，它由基本骨架(低頻信息)和局部細節(高頻信息)構成[8]，同時結合約束圖層[9-10]，對因小波分析被舍去的在災害系統中具有重要作用的關鍵細節加以補充。本研究選取了多級水系和坡度作為約束圖層，因此，基于圖層約束的WT模型具有較好的空間互補性。

三、黃土高原地貌孕災環境的自動綜合應用

本文利用coif2小波進行多尺度分解與重構，派生了5個尺度的DEM(為了便于對比，縮放為同樣大小)，并將相應尺度的水系和坡度周邊的DEM作為細節補充，彌補因小波分析被過濾掉的重要信息，為客觀反映災害發生及形成機制提供了更準確的尺度依據。

1.多尺度DEM的分解與重構

本文選取黃土高原1 km×1 km的DEM作為數據源，利用具有相對較好的對稱性和消失距的coif2小波進行多尺度分解與重構，派生5個尺度的DEM，大小為原尺寸的2-k倍，其比例尺相應逐級降低1/2，并分別從水平、垂直、對角線3個方向作高頻信息過濾，特征信息隨著尺度的降低逐漸減少，細小地貌特征被逐步濾掉，保留主要特征。隨著分解尺度的增加及DEM分辨率的降低，黃土高原地貌的整體骨架特征得以顯現。這表明小波分析技術具有較好的多尺度表達能力。

2.水系與坡度的細節補充

為了準確反映孕災環境對區域災害形成機制的影響，本文對小波變換后的DEM進行了再處理，用水系和坡度兩個圖層對該地區孕災環境的地貌形態和特征加以約束，從而實現了地貌孕災環境基本結構特征的保留。

為了反映不同區域尺度水系與坡度對地貌孕災環境空間結構特征的影響，本文將水系按照干支流分為1～5個級別，在小波變換L1、L3、L53個尺度下選取水系級別分別為R_rank>5，3，2級水系周圍1 km×1 km、2 km×2 km、4 km×4 km的DEM作為細節補充(如圖1—圖3所示)；同樣，將坡度也作了類似處理，將坡度S>2，4，6周圍的DEM作為細節補充(如圖4—圖6所示)，從而彌補在小波分析中可能被當作高頻信息而被刪除的溝壑地貌的細節部分，更客觀地反映了黃土高原特殊孕災環境與災害成因機制之間的復雜關系。水系和坡度的細節圖層對黃土高原特殊地貌結構的保持具有較好的空間約束力。

3.基于LC-WT多尺度DEM的自動綜合表達

將小波分析的初步結果(本文選取了L1、L3、L53個尺度)與相應尺度下的水系、坡度的地貌作為細節補充，即水系1(R_rank>5)、水系2(R_rank>3)、水系3(R_rank>2級)，以及坡度1(S>2)、坡度2(S>4)、坡度3(S>6)。由于小波技術剔除地貌高頻信息(細節地貌)是一個機械的數理過程，沒有考慮區域地貌的整體結構，其中細節部分可能被當作高頻信息剔除，因此，通過不同尺度水系和坡度圖層約束下的地貌細節補充，可以最大限度地保持黃土高原特殊地貌的空間結構，準確反映孕災環境與該地區地貌災害發生及成因上的客觀聯系，從而為不同區域尺度下的防災減災管理提供較合理的地理信息概括。

圖1 河流分布區DEM(R_rank>5)

圖2 河流分布區DEM(R_rank>3)

圖3 河流分布區DEM(R_rank>2)

圖4 坡度分布(S>2)

圖5 坡度分布(S>4)

圖6 坡度分布(S>6)

四、結果分析

建立滿足一定條件的抽象模型是實現空間數據多尺度表達的核心。本文構建的小波綜合方法可以優化表示地表形態的關鍵特征。基于該模型可以實現多個尺度的地貌自動制圖綜合，綜合程度可得到相應的定量化評估。通過LC-WT自動綜合模型對黃土高原地區DEM 3個尺度自動綜合結果的高程直方圖統計可以看出(如圖7—圖10所示)，整體地貌結構特征得到保持，即該地區的整體地貌骨架(低頻信息)總體趨勢沒有發生改變，基本地貌單元特征得到較明顯表達。

圖7 原始DEM直方圖統計

圖8 L1尺度綜合下的直方圖統計

圖9 L3尺度綜合下的直方圖統計

圖10 L5尺度綜合下的直方圖統計

在有限的二維表示平面內，只有準確表達地貌特征，才能客觀地反映地表形態。因此，必須對地貌實施結構化的綜合。對于黃土地貌孕災環境而言，就是要綜合考慮黃土高原的特殊地貌空間結構。通過L1、L3、L53個尺度的綜合結果，以及黃土高原地區地貌的形成過程和自然特征的差異，可以將該地區的地貌孕災環境大致分為以下幾個主要地貌單元：①隴中盆地；②渭河地塹平原；③隴東、陜北高原；④山西高原；⑤鄂爾多斯高原地區(如圖11—圖14所示)。

圖11 黃土高原地區原始DEM

圖12 自動綜合結果(L1)

圖13 自動綜合結果(L3)

圖14 自動綜合結果(L5)

總體而言，黃土高原地區呈現“環狀”空間結構特征(如圖15所示)。該地區自然環境的特點是：四周環山(秦嶺—太行山—陰山—賀蘭山—日月山)，構成外環，是本地區的綠色天然屏障；山前斷裂帶所形成的一系列平原盆地(渭河平原—運城盆地—臨汾盆地—太原盆地—大同盆地—內蒙古河套平原—寧夏河套平原)構成內環；環帶中心為地臺，由于自然侵蝕與人為加速侵蝕疊加，區內自然災害類型多樣，且發生于臺地邊緣地區。通過不同尺度的自動綜合，區域孕災環境的空間結構特征得到較明顯的反映，災害發生規律區間差異較明顯，為不同尺度的區域風險管理與減災規劃提供了科學依據。

圖15 黃土高原地區地貌的環狀空間結構

五、結束語

本文基于地理學、地圖學和災害學理論基礎，初步構建了圖層約束—小波分析(LC-WT)的面狀自然災害風險制圖中孕災環境的自動綜合模型，并以我國黃土高原地區的DEM為例，進行了多尺度的地貌自動綜合表達研究，對結果的有效性作了分析，即是否符合該地區基本地貌結構的地域分異規律，并就多尺度轉換中黃土高原地區地貌孕災環境的區域分異及空間結構特征作了概括和總結，為多尺度的區域災害風險管理與決策提供理論與方法支持。同時，LC-WT自動模型具有較好的空間互補性，對科學認識區域災害系統提供了新的途徑。

[1] 中國科學院黃土高原綜合科學考察隊.黃土高原地區自然環境及其演變[M].北京:科學出版社,1991.

[2] 王靜愛，史培軍，王瑛，等.基于災害系統論的《中國自然災害系統地圖集》編制[J].自然災害學報,2003,12(4):1-8.

[3] TOPFER F.Generalization[M].Beijing: Surveying and Mapping Press,1982.

[4] 武芳,錢海忠,鄧紅艷.面向地圖自動綜合的空間信息智能處理[M].北京:科學出版社, 2008.

[5] 潘東華,王靜愛,賈慧聰,等.圖層約束下的點狀自然災害風險地圖綜合[J].中國圖象圖形學報，2013, 18(4):429-435.

[6] 吳凡,祝國瑞.基于小波分析的地貌多尺度表達與自動綜合[J].武漢大學學報：信息科學版,2001,26(2):170-176.

[7] 潘東華,王靜愛,賈慧聰,等.自然災害風險地圖中的制圖綜合研究——以點狀承災體為例[J].武漢大學學報：信息科學版,2011,36(1):51-55.

[8] 李霖,吳凡.空間數據多尺度表達模型及其可視化[M].北京:科學出版社,2005.

[9] 曹亞妮,吳芳華,江南,等.電子專題地圖多尺度表達的相關問題及策略[J].測繪通報,2014(3):58-62.

[10] 季民,靳奉祥,李云嶺,等.海洋漁業專題屬性數據多尺度綜合與表達[J].測繪通報,2004(6):28-31.

On Automated Generalization for DEM of Loess Plateau in Risk Mapping of Natural Hazard

PAN Donghua，WANG Jingai,JIA Huicong

潘東華，王靜愛，賈慧聰.災害風險地圖中黃土高原地貌自動綜合研究[J].測繪通報，2015(7)：58-61.

10.13474/j.cnki.11-2246.2015.0211

2014-04-22

潘東華(1981—)，男，博士，副研究員，主要從事GIS及RS在自然災害風險分析與制圖中的應用研究。E-mail:eliteast@mail.bnu.edu.cn

P208

：B

：0494-0911(2015)07-0058-04