饒 喆 張靜遠 馮 煒

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

0 引 言

隨著海洋資源勘探、開發(fā)和軍事應(yīng)用需求不斷提高，水下地形匹配制導(dǎo)技術(shù)已成為水下航行體實現(xiàn)精確導(dǎo)航與制導(dǎo)的重要手段［1－2］.作為水下地形匹配制導(dǎo)技術(shù)的重要組成部分，水下地形的研究日益深入［3］.數(shù)字高程模型（digital elevation model，DEM）是一種地形數(shù)字化表達方式，能夠客觀描述和有效提取地形信息，是數(shù)字地形研究的有力工具.分辨率作為DEM的基本屬性，不僅表征了DEM對地形的描述精度，而且還體現(xiàn)了DEM包含的信息量，是確定地形參數(shù)和應(yīng)用尺度的重要指標［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］在實際應(yīng)用中由于坡度信息提取存在不確定性，導(dǎo)致難以確定適宜分辨率；文獻［9－11］沒有確定特定地形下的DEM適宜分辨率；因此，根據(jù)文獻［12］，本文提出了基于局部方差粗糙度的方法以確定水下地形的DEM適宜分辨率.

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

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

根據(jù)地形信息的尺度效應(yīng)，適宜分辨率確定的基本方法是通過定量刻畫尺度效應(yīng)曲線來找出定量指標最大值對應(yīng)的分辨率.因此，確定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］指出，定量指標反映地表的起伏程度，指標值越大表明起伏越劇烈，其最大值對應(yīng)的分辨率即為適宜分辨率.本文結(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ù).在實際應(yīng)用中分析窗口的尺寸通常選定為網(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 實例應(yīng)用

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具有較低的適宜分辨率，且其地形局部方差粗糙度計算值較大，這是由于地形起伏度變化程度更劇烈導(dǎo)致.

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

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

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

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

4 結(jié)束語

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

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

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

［2］劉興堂，戴革林.精確制導(dǎo)武器與精確制導(dǎo)控制技術(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.