饒 喆 張靜遠 馮 煒

（海軍工程大學兵器工程系 武漢 430033）

0 引 言

隨著海洋資源勘探、開發(fā)和軍事應用需求不斷提高，水下地形匹配制導技術(shù)已成為水下航行體實現(xiàn)精確導航與制導的重要手段［1－2］.作為水下地形匹配制導技術(shù)的重要組成部分，水下地形的研究日益深入［3］.數(shù)字高程模型（digital elevation model，DEM）是一種地形數(shù)字化表達方式，能夠客觀描述和有效提取地形信息，是數(shù)字地形研究的有力工具.分辨率作為DEM的基本屬性，不僅表征了DEM對地形的描述精度，而且還體現(xiàn)了DEM包含的信息量，是確定地形參數(shù)和應用尺度的重要指標［4－6］.因此，選取最大限度反映原始數(shù)據(jù)所包含地形信息量的最低分辨率，即適宜分辨率，具有重要意義.

自上個世紀90年代起，國內(nèi)外學者對DEM適宜分辨率的確定問題展開了一系列研究.Hutchinson［7］提出了一種基于坡度中誤差的DEM適宜分辨率的確定方法；Florinsky等［8］探討了DEM適宜分辨率的求解方法；David等［9］對DEM數(shù)據(jù)分辨率對地形模型的影響進行了分析；Walker和Willgoose分析了DEM分辨率大小對地形指數(shù)的影響［10］；劉學軍等［11］研究了 DEM 網(wǎng)格的分辨率、方向和高程數(shù)據(jù)的準確度對地形表達的影響.但文獻［7－8］在實際應用中由于坡度信息提取存在不確定性，導致難以確定適宜分辨率；文獻［9－11］沒有確定特定地形下的DEM適宜分辨率；因此，根據(jù)文獻［12］，本文提出了基于局部方差粗糙度的方法以確定水下地形的DEM適宜分辨率.

1 局部方差粗糙度確定適宜分辨率的理論方法

1.1 確定適宜分辨率的基本思路

根據(jù)地形信息的尺度效應，適宜分辨率確定的基本方法是通過定量刻畫尺度效應曲線來找出定量指標最大值對應的分辨率.因此，確定DEM適宜分辨率的基本思路是將地形因子作為地形信息，通過計算地形因子的局部方差粗糙度與網(wǎng)格尺寸的關(guān)系曲線來確定適宜分辨率.目前，確定DEM適宜分辨率的一般步驟為：

1）生成多分辨率DEM數(shù)據(jù)，即通過原始DEM數(shù)據(jù)生成多分辨率的DEM數(shù)據(jù).

2）計算地形因子的局部方差粗糙度與網(wǎng)格尺寸關(guān)系曲線，即利用第一步生成的多分辨率DEM數(shù)據(jù)計算出多分辨率的地形因子值，再以某定量指標作為地形信息的度量，得出定量指標與分辨率的關(guān)系曲線.

3）確定DEM適宜分辨率，即根據(jù)定量指標的地學意義，從定量指標與分辨率的關(guān)系曲線中找出適宜分辨率.

1.2 基于局部方差粗糙度的適宜分辨率確定方法

文獻［12］指出，定量指標反映地表的起伏程度，指標值越大表明起伏越劇烈，其最大值對應的分辨率即為適宜分辨率.本文結(jié)合地形粗糙度的概念，提出了局部方差粗糙度，并利用局部方差粗糙度作為定量指標來確定適宜分辨率.

對一塊地形區(qū)域，進行數(shù)字化后生成規(guī)則網(wǎng)格 DEM，大小為 M×N，z（i，j）為格網(wǎng)（i，j）的高程值，其中i＝1，2，…，M；j＝1，2，…，N.局部方差粗糙度的具體計算方法如下.

1）計算地形因子局部方差

式中：V（i，j）為網(wǎng)格（i，j）的地形因子局部方差；xk為分析窗口中第k個網(wǎng)格的地形因子值為分析窗口中所有網(wǎng)格的地形因子均值；n0為分析窗口中的網(wǎng)格數(shù).在實際應用中分析窗口的尺寸通常選定為網(wǎng)格大小的奇數(shù)倍，因此，文中選定分析窗口尺寸為3乘3的網(wǎng)格，此時n0＝9.

2）計算地形因子局部方差粗糙度 根據(jù)分析窗口尺寸，原地形區(qū)域被分析窗口劃分以后，通過式（1）計算可得到m×n個地形因子局部方差粗糙度rV，則

式中：D（x）為地形因子基準方差；p，q為不同分辨率下DEM的格網(wǎng)行列數(shù)；ˉx為所有網(wǎng)格的地形因子均值.

2 實例應用

2.1 實驗樣區(qū)基本信息

本實驗采用武漢木蘭湖實測水下地形深度數(shù)據(jù)，選取區(qū)域A和B作為實驗樣區(qū)，其三維地形圖見圖1.根據(jù)該數(shù)據(jù)制作的DEM數(shù)據(jù)的初始分辨率為5m，然后通過局部平均法生成10m，15 m，20m，25m，30m等多分辨率DEM數(shù)據(jù).

圖1 實驗樣區(qū)三維地形圖

2.2 設(shè)計思路

根據(jù)該數(shù)據(jù)制作的多分辨率DEM數(shù)據(jù)，利用上述方法分析并確定DEM適宜分辨率.設(shè)計思路見圖2.

圖2 設(shè)計思路圖

2.3 參數(shù)計算與曲線繪制

選取高程、坡度、平面曲率、剖面曲率四種地形因子來確定DEM適宜分辨率.在坡度坡向計算方面，Burrough［13］提出的三階反距離平方權(quán)差分算法能夠減小高程誤差對坡度的影響，故本文采用該方法計算坡度值.由于 Wood［14］提出的限制二次曲面擬合法具有地形曲面擬合精確、能夠平滑DEM數(shù)據(jù)誤差等優(yōu)點，因此，本文采用該方法計算平面曲率、剖面曲率.

最后，在不同分辨率下，利用式（2）計算出以上4種地形因子的局部方差粗糙度，并繪制關(guān)系曲線見圖3.

圖3 “局部方差粗糙度－分辨率”曲線圖

3 實驗結(jié)果分析及驗證

由圖3可知，對于實驗樣區(qū)A，4種地形因子的局部方差粗糙度在分辨率為15m時達到最大，而實驗樣區(qū)B在分辨率為20m時達到最大，因此，實驗樣區(qū)A和B的適宜分辨率分別為15m和20m.相比于實驗樣區(qū)A，實驗樣區(qū)B具有較低的適宜分辨率，且其地形局部方差粗糙度計算值較大，這是由于地形起伏度變化程度更劇烈導致.

為驗證本文方法的正確性，運用基于信息含量的分析原理［15］，對上述地形數(shù)據(jù)做同樣分析，通過找到DEM信息含量最大時對應的最低分辨率來確定適宜分辨率.

根據(jù)文獻［16］，利用信息熵的分析方法，分別計算高程信息熵，坡度信息熵，平面曲率信息熵和剖面曲率信息熵，得到4種地形因子信息熵隨分辨率變化的曲線見圖4.

圖4 “信息熵－分辨率”曲線圖

根據(jù)圖4得到的適宜分辨率與圖3一致，從而證明文中的分析方法合理、有效.

4 結(jié)束語

本文提出局部方差粗糙度的概念，并將其作為定量指標，繪制局部方差粗糙度與分辨率的關(guān)系曲線，找出其中最大值對應的分辨率作為適宜分辨率.

從計算結(jié)果可知，隨著分辨率的提高，地形局部方差粗糙度逐漸增大，當分辨率達到某一定值時，地形局部方差粗糙度達到最大值，此后地形局部方差粗糙度隨分辨率的提高逐漸減小，從而得到了本實驗數(shù)據(jù)下的DEM適宜分辨率.對于深海海底地形，利用現(xiàn)階段測繪技術(shù)難以獲取高精度、高密度的地形水深數(shù)據(jù)，同時深海海底地形較為復雜，因此該方法是否適用于深海海底地形適宜分辨率的確定，有待進一步的實驗和分析論證.

［1］NYGREN I.Terrain navigation for underwater vehicles［D］.Sweden：Royal Institute of Technology，2005.

［2］劉興堂，戴革林.精確制導武器與精確制導控制技術(shù)［M］.西安：西北工業(yè)出版社，2009.

［3］MEDUNA D K.Terrain relative navigation for sensor－limited systems with application to underwater vehicles［D］.America Stanford：Stanford University，2011.

［4］湯國安，劉學軍，閭國年.數(shù)字高程模型及地學分析的原理與方法［M］.北京：科學出版社，2005.

［5］徐 靜，顧留碗，張耀民，等.規(guī)則網(wǎng)格DEM地形綜合方法分析［J］.測繪工程，2014，23（9）：26－31.

［6］徐道柱，張衛(wèi)柱，田 浩，等.基于等高線約束三角網(wǎng)的格網(wǎng)DEM精度評估方法［J］.測繪工程，2013，22，（2）：75－78.

［7］HUTCHINSON M F.A locally adaptive approach to the interpo－lation of digital elevation models［C］∥NCGIA ed.Proceeding of the T－h(huán)ird International Conference Integrating GIS and Environmental Models.Santa Fe，New Mexico，1996（1）：21－25.

［8］FOLRINSKY I V，KURYAKOVA G A.Determination of grid size for digital terrain modeling in landscape investigations exemplified by soil moisture distribution at a mirco－scale［J］.Internat－ional Journal of Geographical Information Science，2000，14（8）：815－832.

［9］DAVID M W，CUAIS V P.Effects of distal elevation model map scale and data resolution on a topography－based watershed mode［J］.Water Resources Research，1994，30（11）：3041－3052.

［10］WALKER J P，WILLGOOSE G R.On the effect of digital elevation model accuracy on hydrology and geomorphology［J］. Water Resources Research，1999，35（7）：2259－2268.

［11］劉學軍，龔健雅，周啟鳴，等.DEM結(jié)構(gòu)特征對坡度坡向的影響分析［J］.地理與地理信息科學，2004，20（6）：1－6.

［12］秦承志，呼雪梅.柵格數(shù)字地形分析中的尺度問題研究方法［J］.地理研究，2014，33（2）：270－283.

［13］BURROUGH P A，McDONELL R A.Principles of geographical information systems［M］.New York：Oxford University Press，1998.

［14］WOOD J.The geomorphological characterization of digital elevation model［D］.Leicester，UK：University of Leicester，1996.

［15］楊勤科，張彩霞，李領(lǐng)濤，等.基于信息含量分析法確定DEM分辨率的方法研究［J］.長江科學院院報，2006，23（5）：21－24.

［16］朱 偉，王東華，周曉光.基于信息熵的DEM適宜分辨率確定方法研究［J］.遙感信息，2008（5）：79－82.