詹 蕾,湯國安,楊 昕

(1.南京師范大學(xué)虛擬地理環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室,江蘇南京 210046;2.西安文理學(xué)院資源環(huán)境與旅游系,陜西西安 710065）

SRTM DEM高程精度評價

詹 蕾1,2,湯國安1＊,楊 昕1

(1.南京師范大學(xué)虛擬地理環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室,江蘇南京 210046;2.西安文理學(xué)院資源環(huán)境與旅游系,陜西西安 710065）

為了全面認(rèn)識SRTM DEM數(shù)據(jù)精度特征并完善SRTM DEM數(shù)據(jù)精度評定方法,該文以我國1∶5萬比例尺DEM為參考數(shù)據(jù),以具有多種地貌類型的陜西省為實(shí)驗(yàn)樣區(qū),利用高程中誤差模型及空間插值方法對 SRTM DEM進(jìn)行高程精度分析。結(jié)果表明:陜西省的SRTM DEM高程中誤差在3.5～60.7 m,呈現(xiàn)出較為顯著的空間分異特征;并且高程中誤差與實(shí)驗(yàn)樣區(qū)平均坡度有較強(qiáng)的指數(shù)相關(guān)性,擬合的指數(shù)函數(shù)具有較高的模擬精度。

SRTM DEM;高程精度;中誤差;內(nèi)插

0 引言

2000年2月11-22日,美國奮進(jìn)號航天飛機(jī)的雷達(dá)地形測繪(Shuttle Radar Topography M ission, SRTM）獲得了地球表面60°N至56°S、覆蓋陸地表面80%以上的三維雷達(dá)數(shù)據(jù)[1],其中3″×3″(相當(dāng)于90 m柵格分辨率）的SRTM 3已對全球免費(fèi)發(fā)布。這一數(shù)據(jù)將推動空間地理科學(xué)的發(fā)展,并為數(shù)字地形分析提供數(shù)據(jù)支持。國外關(guān)于SRTM數(shù)據(jù)的質(zhì)量研究較多,集中在高程采樣誤差評價方面,部分學(xué)者對SRTM提取的地形參數(shù)描述的精度也進(jìn)行了研究。如 Koch等根據(jù)空間相似性變換原理得到SRTM與其他高精度高程數(shù)據(jù)的匹配算法,可根據(jù)某一參考DTM評估SRTM的精度[2];Sun等采用SLA-02和SRTM數(shù)據(jù)進(jìn)行交叉驗(yàn)證,得到SRTM在表面無遮蔽區(qū)的精度高于其標(biāo)稱精度,在森林區(qū)精度低于標(biāo)稱精度的結(jié)論[3];Brow n等通過計(jì)算GPS點(diǎn)與SRTM相應(yīng)點(diǎn)的高程絕對誤差和相對誤差評價 SRTM的數(shù)據(jù)精度[4];M iliaresis等利用誤差圖將 SRTM與1∶25萬DEM比較,計(jì)算其誤差[5,6];Ludw ig等利用高程剖面和誤差場分析SRTM DEM與參考DEM間的高程差異[7]。

SRTM數(shù)據(jù)有現(xiàn)實(shí)性強(qiáng)、免費(fèi)獲取等優(yōu)點(diǎn),但其特殊的數(shù)據(jù)獲取方式、特定的空間分辨率都給利用SRTM DEM進(jìn)行數(shù)字地形分析帶來諸多不便,有必要深入研究SRTM DEM的數(shù)據(jù)質(zhì)量。本文以我國標(biāo)準(zhǔn)化1∶5萬DEM為對比數(shù)據(jù),以具有多種地貌類型的陜西省為實(shí)驗(yàn)樣區(qū),利用高程中誤差模型及空間插值方法對SRTM DEM進(jìn)行了高程精度分析。

1 實(shí)驗(yàn)樣區(qū)及數(shù)據(jù)

選取陜西省82個地貌類型區(qū)作為實(shí)驗(yàn)樣區(qū)(圖1）。陜西省位于我國中部偏北的黃河中游(東經(jīng)105° 29′～111°15′,北緯31°42′～39°35′）,東西寬200～500 km,面積為20.56萬km2。全省地域狹長,境內(nèi)山塬起伏、河川縱橫、地形復(fù)雜多樣,地勢南北高、中間低,西部高、東部低,地形地貌具有很強(qiáng)的代表性。

圖1 實(shí)驗(yàn)樣區(qū)分布概況Fig.1 General introduction of test areas

實(shí)驗(yàn)采用國家基礎(chǔ)地理數(shù)據(jù)庫1∶5萬DEM數(shù)據(jù)和SRTM DEM數(shù)據(jù)。前者采用1980國家大地坐標(biāo)系,高斯-克呂格投影(Gauss-Kruger）,柵格分辨率為25 m;后者采用WGS1984基準(zhǔn)面,其數(shù)據(jù)下載及基本情況參見美國NASA、CGIAR及USGS官方網(wǎng)站。由于兩者采用的基準(zhǔn)面不同,地理位置會產(chǎn)生系統(tǒng)性偏差。為能將兩種DEM疊合比較,須將其轉(zhuǎn)換為同一投影坐標(biāo)系,并消除SRTM DEM數(shù)據(jù)參考橢球體基準(zhǔn)面不同帶來的系統(tǒng)偏差。我國高程基準(zhǔn)是1985國家高程基準(zhǔn),而SRTM DEM采用EGM 96高程基準(zhǔn),高程基準(zhǔn)不同導(dǎo)致高程測量值會產(chǎn)生垂直偏差。郭海榮等計(jì)算了我國1985國家高程基準(zhǔn)與EGM 96之間的系統(tǒng)差(垂直偏差約35.7 cm）,且系統(tǒng)差自東向西增大[8]。根據(jù)這一規(guī)律,對SRTM DEM數(shù)據(jù)進(jìn)行垂直誤差校正。完成投影轉(zhuǎn)換、水平校正和垂直校正后,依據(jù)1∶5萬DEM數(shù)據(jù)的區(qū)域范圍,對82個實(shí)驗(yàn)樣區(qū)SRTM DEM進(jìn)行裁切,以供進(jìn)一步研究。

2 SRTM DEM高程中誤差空間內(nèi)插方法

2.1 SRTM DEM高程中誤差

目前常用的DEM高程精度評價標(biāo)準(zhǔn)為DEM高程中誤差模型,它是用子樣方差表示母體方差。DEM中誤差的計(jì)算方法是,假設(shè)檢驗(yàn)點(diǎn)的高程為Zk(k=1,2,…,n）,在建立的DEM上對應(yīng)這些點(diǎn)的高程為zk,則DEM的中誤差為:

我國1∶5萬DEM采用28個分布在圖幅內(nèi)和圖幅邊緣的檢驗(yàn)點(diǎn),按上述方法對DEM質(zhì)量進(jìn)行大體的精度評定;實(shí)際操作中,28個檢驗(yàn)點(diǎn)過于稀少,具有很大的偶然性,因此有必要增加檢驗(yàn)點(diǎn)數(shù)目。本實(shí)驗(yàn)基于A rcGIS軟件,將1∶5萬DEM原始等高線轉(zhuǎn)為節(jié)點(diǎn)數(shù)據(jù),作為計(jì)算中誤差的檢驗(yàn)點(diǎn),再轉(zhuǎn)換為與SRTM DEM相對應(yīng)的90 m分辨率的柵格數(shù)據(jù);然后編寫AML程序,對SRTM DEM和得到的柵格圖層進(jìn)行逐柵格運(yùn)算,利用中誤差公式計(jì)算SRTM DEM與此柵格數(shù)據(jù)的高程中誤差。

2.2 SRTM DEM高程中誤差空間內(nèi)插

反距離加權(quán)插值法(Inverse Distance Weighted,IDW）基于相似相近原理,即兩個物體距離越近,其性質(zhì)越相似。該方法內(nèi)插結(jié)果較穩(wěn)定,也接近實(shí)際情況。本文采用反距離加權(quán)插值法對SRTM DEM中誤差進(jìn)行空間內(nèi)插,即以待插點(diǎn)為中心,在待插點(diǎn) X、Y處內(nèi)插變量Z值時,在局部鄰域中計(jì)算各數(shù)據(jù)點(diǎn)的加權(quán)平均值。公式如下:

式中:zi是離散點(diǎn)高程值,d2i是距離的平方倒數(shù)權(quán),一般將冪設(shè)為2。

3 SRTM DEM中誤差模擬結(jié)果與精度分析

分別計(jì)算82個樣區(qū)的高程中誤差,同時,利用A rcGIS中的地統(tǒng)計(jì)分析模塊(Geostatistical Analyst）進(jìn)行反距離加權(quán)插值,得到陜西省SRTM DEM的中誤差空間分異圖(圖2）,發(fā)現(xiàn)陜西省 SRTM DEM高程中誤差由南向北呈現(xiàn)從秦巴山地的較大值到關(guān)中盆地的較小值,到黃土丘陵區(qū)中誤差值逐漸增大,到風(fēng)沙河谷區(qū)中誤差逐漸減小。

圖2 SRTM DEM高程中誤差空間分異Fig.2 Spacial distribution of elevation RMSEs of the SRTM DEM

陜西省南部的秦巴山區(qū)SRTM DEM的高程誤差較大,這是因?yàn)樯降貐^(qū)地形結(jié)構(gòu)復(fù)雜,地形起伏度、切割深度均較大,SRTM DEM 90 m的分辨率無法細(xì)致地刻畫地形起伏變化。另外,由于雷達(dá)干涉測量本身的缺陷,在山體陡峭的華山、太白山、石泉樣區(qū),山體背坡坡度較大時可能會產(chǎn)生雷達(dá)陰影從而造成數(shù)據(jù)缺失,僅通過內(nèi)插等方法填補(bǔ)這些空洞點(diǎn),可能會出現(xiàn)更大的誤差。秦嶺—巴山中部的漢中盆地、安康低山丘陵區(qū)平均海拔較低,且地勢變化較平緩,因此高程中誤差值較周圍的中高山區(qū)小。處于關(guān)中腹地的渭河平原地形平緩,中誤差較小;但該地形區(qū)面積相對狹小,實(shí)驗(yàn)樣區(qū)受到限制,受南部秦嶺山區(qū)和渭北“北山”的影響,該區(qū)的中誤差值比實(shí)際值偏大。在黃土塬區(qū)和黃土丘陵溝壑區(qū),地形起伏度、切割深度均逐漸增大,高程中誤差也逐漸增大。在風(fēng)沙河谷區(qū)有連片的低緩沙丘,樣區(qū)內(nèi)高差較小,地形起伏平緩,中誤差值最小。經(jīng)過多次試驗(yàn),發(fā)現(xiàn)中誤差與實(shí)驗(yàn)樣區(qū)平均坡度有較強(qiáng)的指數(shù)相關(guān)性(圖3）,擬合公式如下:

其中slpm是樣區(qū)平均坡度。

圖3 中誤差—坡度關(guān)系示意Fig.3 The relationship of RMSEsand average slope of samples

為了驗(yàn)證模擬結(jié)果的可靠性,將未參與模擬的6個樣區(qū)SRTM DEM作為檢驗(yàn)樣區(qū),其中誤差計(jì)算值及模擬值見表1。中誤差的模擬值與計(jì)算值的最大絕對誤差為5.98 m,最小誤差僅為1.29 m,最大相對誤差為0.29,最小僅為0.05。可見,式(2）的模擬結(jié)果精度較高。根據(jù)以上分析,可以認(rèn)為坡度越大、起伏越大的區(qū)域,其高程中誤差也越大,而坡度越小、越平緩的區(qū)域,其高程中誤差越小。

表1 SRTM DEM中誤差模擬結(jié)果檢驗(yàn)Table 1 Testing of the RMSE simulation result of SRTM DEM

4 結(jié)語

本文采用中誤差模型,通過內(nèi)插得到陜西省SRTM DEM高程中誤差空間分異圖。結(jié)果表明:陜西省SRTM DEM的高程中誤差在3.5～60.7 m,其在山地區(qū)最大,在黃土丘陵區(qū)較大,在關(guān)中盆地和風(fēng)沙河谷區(qū)較小;SRTM DEM高程中誤差與實(shí)驗(yàn)樣區(qū)平均坡度有較強(qiáng)的相關(guān)性,可以用指數(shù)函數(shù)模擬,且模擬精度較高。本研究以我國1∶5萬DEM為假定真值,代替了實(shí)測的點(diǎn)高程值。盡管我國1∶5萬DEM覆蓋了各種地貌類型,且精度較高,但與實(shí)測的點(diǎn)高程值還有一定誤差。在今后的研究中,應(yīng)進(jìn)行野外實(shí)地采樣,以期得到更精確的SRTM DEM。

[1] 汪凌.美國航天飛機(jī)雷達(dá)地形測繪使命簡介[J].測繪通報(bào), 2000(12）:38-40.

[2] KOCH A,HEIPKEC.Quality assessment of digital surfacemodels derived f rom the Shuttle Radar Topography M ission(SRTM） [A].Proc.of the IEEE 2001 International Geoscience and Remote Sensing Symposium[C].Sydney,Australia:2001.

[3] SUN G,RANSON K J,KHARU K V I,et al.Validation of surface height from Shuttle Radar Topography M ission using shuttle laser altimeter[J].Remote Sensing of Environment,2003 (88）:401-411.

[4] BROWN C G,SARABANDIJ K,PIERCE L E.Validation of the Shuttle Radar Topography Mission height data[J].IEEE Transactions on Geoscience and Remote Sensing,2005,43(8）: 1707-1715.

[5] M ILIARESIS G C,PARASCHOU C V.Vertical accuracy of the SRTM DTED level 1 of Crete[J].International Journalof Applied Earth Observation and Geoinformation,2005(7）:49-59.

[6] MA L,L I Y.Evaluation of SRTM DEM over China[A].IEEE Geoscience and Remote Sensing Symposium[C].Denver:2006.

[7] LUDW IGR,SCHNEIDER P.Validation of digital elevation models from SRTM X-SAR for applications in hydrologic modeling [J].ISPRS Journal of Photogrammetry and Remote Sensing, 2006(60）:339-358.

[8] 郭海榮,焦文海,楊元喜.1985國家高程基準(zhǔn)與全球似大地水準(zhǔn)面之間的系統(tǒng)差及其分布規(guī)律[J].測繪學(xué)報(bào),2004,33(2）:100-104.

Evaluation of SRTM DEM s′Elevation Accuracy:A Case Study in Shaanxi Province

ZHAN Lei1,2,TANG Guo-an1,YANG Xin1
(1.Key Laboratory of Virtual Geographic Environment,M inistry of Education,Nanjing Normal University,Nanjing 210046;
2.Department of Resources Environment&Tourism,Xi′an University of Arts&Science,Xi′an 710065,China）

Treating the state 1∶50 000 DEM as reference contrastive data,investigating in Shaanxi Province w hich hasmany different geomo rphological types,the SRTM DEM′s elevation accuracy are researched in this paper,using elevation RMSEmodel and method of spatial interpolation.Above all,the elevation RMSE of each test area is calculated through w riting p rogram s in A rcGISaccording to RM SEmodel.Then,interpolating the calculated RMSEs using IDW method,the spatial differentiation of SRTM DEM′s RM SEof w hole Shaanxi Province is gained.The SRTM DEM′s elevation RMSEs are from 3.5 m to 60.7 m and have themuch mo re notable special distribution characteristic in Shaanxi Province:the SRTM DEM′s elevation RMSEs are biggest in the mountain areas,bigger in loess hill and tableland areas and smaller in p lain areas.There is the exponential correlation between the elevation RM SEs and the average slope of samples.The fitting exponential function has better accuracy.This study has certain reference significance to know ing the characteristic of SRTM DEM accuracy tho roughly and imp roving the evaluation method of SRTM DEM accuracy.

SRTM DEM;elevation accuracy;RMSE;spatial interpolation

P208

A

1672-0504(2010）01-0034-03

2009-07-22;

2009-09-24

國家自然科學(xué)基金項(xiàng)目(40671148）

詹蕾(1983-）,女,碩士研究生,主要從事地理信息系統(tǒng)方面的研究。＊通訊作者E-mail:tangguoan@njnu.edu.cn