中國5類典型區(qū)域常用DEM高程精度評價

蔣廣鑫,謝元禮,3,高志遠,周 鵬

(1.西北大學(xué) 城市與環(huán)境學(xué)院,西安 710127;2.陜西省地表系統(tǒng)與環(huán)境承載力重點實驗室,西安 710127;3.西北大學(xué) 地表系統(tǒng)與災(zāi)害研究院,西安 710127)

數(shù)字高程模型(Digital Elevation Model,DEM)作為地表高程信息的數(shù)據(jù)源,對于水土流失以及地形因子提取等地學(xué)研究具有關(guān)鍵性作用。作為地學(xué)研究中廣泛應(yīng)用的基礎(chǔ)數(shù)據(jù)之一,其精度的優(yōu)劣直接關(guān)系到科研結(jié)果的可信度及科學(xué)性。現(xiàn)階段,針對DEM精度評價的方法主要有2種:(1) 利用更高精度的地面測量控制點對DEM數(shù)據(jù)進行對比分析,如GPS實測點數(shù)據(jù)[1-3]和高程控制點數(shù)據(jù)[4-5];(2) 利用更高精度的DEM數(shù)據(jù)與之進行對比,獲得精度評價,例如最常用的方法,通過與大比例尺地形圖生成的DEM對數(shù)據(jù)進行分析評價[6-8]。然而當(dāng)研究范圍是大區(qū)域時,將會面臨以下問題:第一,無法大范圍地獲取地面測量數(shù)據(jù)。第二,由于研究區(qū)域的范圍較廣以及數(shù)據(jù)的保密性,無法獲取更高的DEM數(shù)據(jù)來評價已有數(shù)據(jù)?；诖?目前已有不少學(xué)者以ICESAT/GLAS高程數(shù)據(jù)作為地面參考點數(shù)據(jù)對大尺度區(qū)域DEM進行精度分析。胡加佩等[9]利用GLAH14數(shù)據(jù)對全國范圍的SRTM V4和ASTER GDEM數(shù)據(jù)進行了對比,系統(tǒng)分析了兩種DEM的誤差統(tǒng)計特征和空間分布特征,結(jié)果表明,兩種DEM都與地形有著較強的相關(guān)性。武文嬌等[10]以山西省為研究區(qū),基于ICESat/GLAS GLAH14測高數(shù)據(jù)對SRTMGL1和ASTER GDEM V2數(shù)據(jù)的垂直精度進行了對比,結(jié)合地形因子以及利用地形剖面分析的方法分析了數(shù)據(jù)差異。結(jié)果表明了SRTM1 DEM垂直精度高于ASTER GDEM V2,并且兩種數(shù)據(jù)均受坡度影響嚴(yán)重,土地利用類型和地貌類型影響較小。

目前的DEM精度評價研究多是針對某一類DEM精度的討論,或兩類DEM之間的對比,研究目的更多的是突出說明某種DEM數(shù)據(jù)的整體優(yōu)越性。常見評價的DEM種類主要集中在ASTER GDEM V2,SRTMGL1和SRTM V4等3類高程數(shù)據(jù),且研究范圍上,以全國范圍或者局部小范圍區(qū)域為主,缺乏對典型區(qū)域精度的研究?；诖?本文的研究目的是通過探究常用DEM在不同區(qū)域的精度表現(xiàn),以便為之后的地學(xué)研究數(shù)據(jù)選取提供一定的指導(dǎo)作用,因此本文選取這3種DEM在5類地貌典型區(qū)進行對比研究。

本文以ICESAT/GLAS GLAH14點數(shù)據(jù)為參照,以地貌為特征,選取中國5類典型區(qū)域作為研究區(qū)域,對SRTMGL1,SRTM V4,ASTER V2數(shù)據(jù)的質(zhì)量進行定量化分析,以便于了解3種高程產(chǎn)品在不同研究區(qū)的誤差分布特征,為數(shù)據(jù)選用及誤差修正提供依據(jù)。

1 試驗區(qū)與數(shù)據(jù)

1.1 研究區(qū)概況

本文依據(jù)研究區(qū)內(nèi)地貌類型構(gòu)成的不同,選擇中國的5類典型地貌區(qū):華北平原、黃土高原、青藏高原、云貴高原和塔里木盆地。其中華北平原區(qū)域地貌主要由平原區(qū)和低山丘陵區(qū)組成;黃土高原其間溝壑縱橫,由山地區(qū)、黃土丘陵區(qū)、黃土塬區(qū)、黃土臺塬區(qū)、河谷平原區(qū)組成,是地貌研究的熱點區(qū)域;青藏高原主要由高原山區(qū)和高原平原組成,平均海拔超過了5 000 m;塔里木盆地位于中國新疆南部,是中國面積最大的內(nèi)陸盆地,中間盆地,邊緣礫石沙漠;云貴高原作為中國四大高原之一,區(qū)域內(nèi)大部分為山區(qū)地形,是世界上喀斯特地貌最發(fā)育的典型地區(qū)之一。這5個區(qū)域都屬于地貌研究熱點區(qū)域,具有較高的研究價值和地貌代表性。

1.2 研究數(shù)據(jù)

本文研究數(shù)據(jù)分為:

(1) 分析數(shù)據(jù),包括SRTMGL1,SRTM V4和ASTER GDEM。其中,SRTMGL1即SRTM 1 Arc-Second Global數(shù)據(jù),對應(yīng)的平面分辨率精度為30 m,數(shù)據(jù)下載網(wǎng)址為http:∥earthexplorer.usgs.gov/。SRTM V4采用的是90 m分辨率SRTM V4.1版本數(shù)據(jù),來源于(http:∥srtm.csi.cgiar.org/);ASTER GDEM選擇版本為ASTERV2,來源于(http:∥www.gscloud.cn/),平面分辨率為30 m,是由對地觀測衛(wèi)星Terra的近紅外波段垂直向下和后視兩個傳感器獲取的立體像對生成的DEM數(shù)據(jù)。

(2) 參考數(shù)據(jù),包括ICESAT/GLAS數(shù)據(jù)和中國地貌類型數(shù)據(jù)。ICESAT/GLAS數(shù)據(jù)為2003—2009年觀測并經(jīng)處理后的全球(南北緯86°之間)地球表面的高程散點數(shù)據(jù),具有覆蓋范圍廣、精度高的特點。本文采用的是GLAH14數(shù)據(jù),來源于美國國家冰雪數(shù)據(jù)中心(http:∥nsidc.org/data)。中國地貌類型數(shù)據(jù)包括2005年中國土地利用現(xiàn)狀遙感監(jiān)測數(shù)據(jù)和中國100萬地貌類型空間分布數(shù)據(jù)。2005年中國土地利用現(xiàn)狀遙感監(jiān)測數(shù)據(jù)是以各期Landsat TM/ETM遙感影像為主要數(shù)據(jù)源,通過人工目視解譯生成,空間分辨率為1 km,數(shù)據(jù)來源為中國科學(xué)院資源環(huán)境數(shù)據(jù)云平臺(http:∥www.resdc.cn/)。中國100萬地貌類型空間分布數(shù)據(jù)來源于《中華人民共和國地貌圖集(1∶100萬)》,數(shù)據(jù)來源為中國科學(xué)院資源環(huán)境數(shù)據(jù)云平臺(http:∥www.resdc.cn/)。

1.3 數(shù)據(jù)預(yù)處理

(1) 格式統(tǒng)一。SRTMGL1,SRTM V4,ASTER GDEM V2均為TIF格式的柵格數(shù)據(jù),GLA14為二進制點云數(shù)據(jù),利用HDFVier工具,從GLA14原始二進制文件中讀取測高數(shù)據(jù),并將其轉(zhuǎn)換成shape文件[11]。

(2) 基準(zhǔn)統(tǒng)一。在高程基準(zhǔn)方面,可以利用式(1) 將GLAH14的高程統(tǒng)一到WGS84參考橢球高程下。

HWGS84=h-N-offset

(1)

式中:h是激光光斑中心相對于T/P橢球的高程;N是EGM96水準(zhǔn)面與T/P橢球面的差值;offset是兩個橢球體高程差異;offset一般取值為0.7 m;h，N使用HDFVier工具直接提取。

(3) GLA14點云篩選。由于ICESAT數(shù)據(jù)存在粗差,因此,基于ICESAT數(shù)據(jù)對SRTM,ASTER GDEM等DEM精度進行評估時,需要對ICESAT點云數(shù)據(jù)進行篩選。通過對已有的研究成果分析并結(jié)合本文研究區(qū)的地形特征,ICESAT數(shù)據(jù)粗差剔除原則[9]為:ICESAT數(shù)據(jù)中,云頂端的回波信號可導(dǎo)致GLA14數(shù)據(jù)含有超過±100 m左右的粗差,考慮到SRTM數(shù)據(jù)精度優(yōu)于ASTER數(shù)據(jù),故通過計算對應(yīng)位置上SRTM DEM和ICESat/GLAS差值,剔除差值在±100 m以外的ICESat點[12-13]。

(4) 點數(shù)據(jù)地類屬性提取。在GLA14點云篩選結(jié)果的基礎(chǔ)上,采用雙線性內(nèi)插的方法獲取對應(yīng)位置的屬性值。

2 數(shù)據(jù)精度分析方法

研究中精度指標(biāo)選取了高程誤差d,平均誤差Mean,標(biāo)準(zhǔn)偏差SD和中誤差RSME作為精度評價指標(biāo)。

d=hDEM-hICESat/GLAH14

(2)

(3)

(4)

(5)

式中:利用高程誤差d表示每個實際點位DEM與ICESat/GLAS GLAH14數(shù)據(jù)之間的差值;利用平均誤差Mean表示DEM數(shù)據(jù)集相較于ICESat/GLAS GLAH14數(shù)據(jù)集的數(shù)據(jù)整體精度。標(biāo)準(zhǔn)差SD則代表了高程誤差間的離散程度。中誤差RSME則代表了測量值與真值的偏離程度,也是評價數(shù)據(jù)精度的最直接標(biāo)準(zhǔn)。

利用以上幾種精度評價指標(biāo),首先對3種DEM在不同典型區(qū)域的總體精度進行統(tǒng)計,然后,結(jié)合坡度、坡向、土地覆蓋、地貌因素等進行3種DEM精度的分析。坡度因子分析中,將所有典型區(qū)域坡度進行統(tǒng)一的離散化分級,級差劃分為0°～1°,1°～5°,5°～10°,10°～15°,15°～20°,20°～25°,25°～30°,30°～35°,35°～40°,以及>40°,分別計算每個坡度分級里不同DEM精度表現(xiàn)。坡向因子方面,將坡向以45°為一個區(qū)間分為8個坡向區(qū)域,計算每個坡向區(qū)間內(nèi)各個DEM的精度表現(xiàn)。土地覆蓋和地貌因素方面,計算各個DEM在不同地貌因素和土地覆蓋中的SD和RSME,對3種DEM進行精度比對分析。

3 結(jié)果及分析

3.1 不同區(qū)域的DEM精度總體評價

3種DEM在5類典型區(qū)域的精度評價結(jié)果見表1,結(jié)論如下:(1) 總體上,SRTMGL1的數(shù)據(jù)精度明顯優(yōu)于SRTM V4和ASTER V2,SRTMGL1的平均誤差為-2.083 m,系統(tǒng)誤差為(2.083±8.201)m,數(shù)據(jù)精度為8.462 m,離散程度較另外2種DEM也偏小;ASTER V2總體精度則略微優(yōu)于SRTM V4,ASTER V2平均精度為2.540 m,系統(tǒng)誤差為(2.540±16.312)m,數(shù)據(jù)精度為16.509 m,SRTM V4平均精度為2.208 m,系統(tǒng)誤差為(2.208±17.579)m,數(shù)據(jù)精度為17.717 m;(2) 結(jié)合試驗區(qū),SRTMGL1在5類試驗區(qū)中精度表現(xiàn)均為最優(yōu),數(shù)據(jù)精度保持在10 m左右,SRTM V4和ASTER V2則在不同的試驗區(qū)精度表現(xiàn)各有優(yōu)劣。華北平原和塔里木盆地兩個試驗區(qū)中,SRTM V4精度優(yōu)于ASTER V2;而在黃土高原,青藏高原和云貴高原3個試驗區(qū)中,SRTM V4精度迅速下降,精度明顯低于ASTER V2,分析其原因,極大程度是因為ASTER GDEM通過對立體影像對解譯獲得,其精度易受云或者地標(biāo)建筑等影響,而SRTM V4受地形起伏或是植被覆蓋度等影響,因此,在平坦區(qū)域,SRTM V4精度優(yōu)于ASTER V2,當(dāng)?shù)匦纹鸱^大,如山地地區(qū),ASTER V2精度優(yōu)于SRTM V4。(3) 3種DEM隨著試驗區(qū)的不同,精度表現(xiàn)均有著不同程度的變化,華北平原和塔里木盆地,3種DEM均具有較好的精度,三者分異性不大,當(dāng)試驗區(qū)為高原地區(qū),3種DEM則產(chǎn)生了較大差異,充分說明,3種DEM與地形有著密切的關(guān)系[14]。

表1 SRTMGL1,SRTM V4和ASTER V2精度指標(biāo) m

3.2 坡度對DEM精度影響

3種DEM中誤差在不同研究區(qū)域隨坡度分級的變化特征如圖1所示,3種DEM在所有區(qū)域都保持著隨著坡度的增加,中誤差逐漸遞增的趨勢。在所有的研究區(qū)域中,SRTMGL1的數(shù)據(jù)精度都要優(yōu)于SRTM V4,在地勢較為平坦的地區(qū)(坡度<1°),SRTMGL1和SRTM V4中誤差較為接近,其中,華北平原兩種DEM中誤差之差為0.02 m,塔里木盆地兩種DEM中誤差之差為0.05 m,黃土高原、青藏高原、云貴高原兩類DEM中誤差之差依次為0.25 m,0.4 m,1.5 m,精度差異細微,2類DEM可以互為補充,在此坡度分級下,由于成像方式的不同,ASTER V2與SRTM DEM存在著3～13 m的中誤差之差,精度較差。隨著研究區(qū)內(nèi)坡度增大,不同研究區(qū)內(nèi)3種DEM中誤差變化規(guī)律逐漸不同,SRTMGL1和ASTER V2在各研究區(qū)內(nèi)都保持著相對平穩(wěn)的中誤差之差,華北平原、黃土高原、青藏高原、塔里木盆地、云貴高原2類數(shù)據(jù)之差依次為(0.6±4.0)m,(2.8±0.6),(11.4±1.4)m,(7.6±3.8)m,(6.5±3.3)m。分析DEM在所有5個試驗區(qū)中精度表現(xiàn),SRTMGL1一般保持著最高精度,但是在華北平原地區(qū),當(dāng)坡度>30°時,SRTMGL1精度迅速下降,ASTER V2精度優(yōu)于SRTMGL1,為3種數(shù)據(jù)精度最高。而SRTM V4和ASTER V2隨著坡度分級,各有優(yōu)劣:在坡度較小(<10°)時,SRTM V4精度優(yōu)于ASTER V2,當(dāng)坡度>10°,ASTER V2精度開始優(yōu)于SRTM V4,且SRTM V4精度隨著坡度的增加,分異程度明顯大于ASTER V2,遞增速率明顯,精度較差,ASTER V2則與SRTMGL1保持著接近的遞增速率,相對SRTM V4而言具有較高的精度[14]。

圖1 各研究區(qū)3種DEM中誤差和坡度的關(guān)系

3.3 坡向?qū)EM精度影響

將坡向以45°為一個區(qū)間,劃分為8個坡向級,北方向(337.5°～22.5°)、東北方向(22.5°～67.5°)、東方向(67.5°～112.5°)、東南方向(112.5°～157.5°)、南方向(157.5°～202.5°)、西南方向(202.5°～247.5°)、西方向(247.5°～292.5°)、西北方向(292.5°～337.5°),參考Zhang等[15]在研究地表覆蓋對SRTM DEM精度的影響時控制變量的思想,為了研究3種DEM在不同區(qū)域的精度表現(xiàn),將3種DEM在5類典型區(qū)域的各坡向中誤差歸一化到區(qū)間內(nèi),繪制在雷達圖上,以研究3種DEM在不同區(qū)域的精度表現(xiàn),結(jié)果如圖2所示。其中,ASTER V2在5個區(qū)域的西南坡向上均出現(xiàn)了較大的中誤差,SRTM V4的中誤差主要出現(xiàn)在正西,正東,正北,正南方向。SRTMGL1則呈現(xiàn)不規(guī)律性。

圖2 各研究區(qū)3種DEM中誤差和坡向的關(guān)系

在黃土高原研究區(qū),3種DEM均呈現(xiàn)對稱分布。而在地勢平坦的華北平原地區(qū),SRTM V4和SRTMGL1精度隨坡向分異性大致相似,在地勢復(fù)雜的云貴高原研究區(qū),3種DEM精度分異性大體一致,分析其原因,可能是因為云貴高原的高密度植被覆蓋導(dǎo)致DEM的精度均有所下降。在地勢復(fù)雜而地表覆蓋不高的青藏高原研究區(qū),3種DEM分別呈現(xiàn)出截然不同的分異性,SRTM V4,ASTER V2和SRTMGL1分別在東南,西南,西北方向精度最低,并明顯規(guī)律。隨坡向的變化,3種DEM表現(xiàn)出高程測量偏離值分異特征,原因可能與衛(wèi)星傳感器在上升軌道和下降軌道的航向以及SRTM傳感器雷達與地表的入射角度有關(guān)。

將每一種DEM在不同研究區(qū)的隨坡向的精度表現(xiàn)繪制雷達圖,結(jié)果如圖3所示3種DEM在華北平原、塔里木盆地、黃土高原和青藏高原4個研究區(qū)中,隨著坡向的變化,每一種DEM在4個研究區(qū)的精度變化程度大體相當(dāng),這表明在同一地區(qū),坡向?qū)EM數(shù)據(jù)精度的影響不大。云貴高原則在東方向,東南方向以及東北方向出現(xiàn)了中誤差變小的情況,而SRTMGL1和SRTM V4兩種數(shù)據(jù)在5個研究區(qū)內(nèi)誤差分布相似,推測可能與SRTM傳感器與地表入射角度有關(guān)。

圖3 各DEM不同研究區(qū)中誤差和坡向的關(guān)系

3.4 土地利用類型對DEM精度影響

基于土地類型對研究區(qū)域進行分類,將試驗區(qū)分為耕地,林地,草地,水域,人造表面,未利用土地等5種類型。分別計算每個試驗區(qū)各土地類型范圍內(nèi)的3種DEM與ICESat/GLAS控制點的中誤差,進行統(tǒng)計,繪制柱狀圖(圖4)以及將每種DEM的5個試驗區(qū)隨土地類型的變化進行縱向比較,制成折線圖(圖5)。

圖4 各研究區(qū)3種DEM中誤差和土地利用類型的關(guān)系

圖5 各DEM不同研究區(qū)中誤差和土地利用類型的關(guān)系

試驗結(jié)果表明,同一試驗區(qū)內(nèi),隨著土地覆蓋類型的變化,SRTMGL1精度均優(yōu)于SRTM V4和ASTER V2。地勢平坦地區(qū)(華北平原和塔里木盆地),SRTM V4和SRTMGL1兩種DEM中誤差之差較為接近,除林地類型以外,中誤差之差最大值為2.2 m,最小值僅為0.01 m,充分說明平坦地區(qū),不同的土地覆蓋類型,SRTM V4與SRTMGL1數(shù)據(jù)精度相當(dāng),ASTER V2精度最低。而在高原地區(qū),SRTM V4數(shù)據(jù)集精度略微低于ASTER V2,且在林地類型區(qū)域,SRTM V4數(shù)據(jù)偏差均高于ASTER V2,說明在高密度植被覆蓋區(qū)域,ASTER V2的精度優(yōu)于SRTM V4[16]。2種DEM雖然生成原理不同,但均受植被覆蓋影響嚴(yán)重,因此林地的精度較低。

在對3種DEM的橫向?qū)Ρ戎锌梢园l(fā)現(xiàn),SRTM V4數(shù)據(jù)集隨土地利用類型的精度變化,由高到低依次為:人造表面,耕地,未利用土地,水域,草地,林地。其中SRTMGL1與SRTM V4保持統(tǒng)一的精度變化趨勢,但是SRTMGL1數(shù)據(jù)集整體精度高于SRTM V4。ASTER V2的精度隨著土地利用覆蓋類型的變化,中誤差隨著人造表面,未利用土地,水域,草地,耕地,林地而依次增大[17]。

3.5 地貌類型對DEM精度影響

基于地貌類型對SRTMGL1,SRTM V4,ASTER V2進行分類,根據(jù)中國100萬地貌類型空間分布數(shù)據(jù)集分類標(biāo)準(zhǔn)將地貌類型分為平原、臺地、丘陵、小起伏山脈、中起伏山脈、大起伏山脈、極大起伏山脈,研究在5個試驗區(qū)內(nèi),地貌類型的變化對于3種DEM的精度影響(圖6和圖7)。

圖6 各研究區(qū)3種DEM中誤差和地貌類型的關(guān)系

圖7 各DEM不同研究區(qū)中誤差和地貌的關(guān)系

結(jié)果表明,SRTMGL1在5個試驗區(qū)均為精度最高的數(shù)據(jù)集,且ASTER V2在地貌類型為平原的地區(qū)精度是3種DEM中最低的。5個試驗區(qū)中,華北平原、青藏高原、塔里木盆地、云貴高原4個試驗區(qū)中SRTM V4在平原、臺地、丘陵的精度均高于ASTER V2。當(dāng)?shù)孛差愋妥優(yōu)樯矫},且隨著海拔的升高,ASTER V2精度高于SRTM V4。

總體而言,3種DEM均隨著海拔的升高,中誤差逐漸遞增,其中,SRTM V4隨著海拔的升高,中誤差遞增速率較快,在高原試驗區(qū)的山脈地區(qū)中誤差升高迅速,精度下降速度較快。SRTMGL1和ASTER V2中誤差遞增速率則呈現(xiàn)緩慢的遞增趨勢。3種DEM在5個試驗區(qū)的臺地地貌中精度明顯下降,中誤差明顯高于平原和丘陵地貌,在華北平原試驗區(qū)較為明顯[18]。

4 結(jié) 論

(1) 對比5個研究區(qū)共979 490個ICESat/GLAS GLAH14數(shù)據(jù)點,SRTM V4,ASTER V2和SRTMGL1數(shù)據(jù)的總體平均誤差分別為2.2 m,2.5 m和2.1 m,誤差標(biāo)準(zhǔn)差為17.5 m,16.3 m和8.2 m,中誤差分別為17.7 m,16.5 m,8.4 m?？傮w而言,SRTMGL1數(shù)據(jù)集的垂直精度要高于ASTER V2,而ASTER V2的總體高程精度則要略微優(yōu)于SRTM V4。

(2) SRTMGL1,ASTER V2和SRTM V4均受到坡度、土地利用類型和地貌類型影響較大,而坡向?qū)τ贒EM的影響較小,且誤差分布較為均勻。坡度對3種DEM的影響總體趨勢為:坡度越高,精度越差。SRTMGL1在各個坡度分級區(qū)間數(shù)據(jù)精度均為最高,ASTER V2和SRTM V4在坡度較小(<15°)時,ASTER V2精度低于SRTM V4,當(dāng)坡度>15°,ASTER V2精度開始優(yōu)于SRTM V4。3種DEM在耕地和人造表面覆蓋區(qū)域的精度最高,在林地覆蓋區(qū)域,DEM精度最低[19]。

(3) 不同地貌分區(qū)對3種DEM的精度具有較大的影響。SRTM V4隨著海拔和地形起伏度的增大,數(shù)據(jù)中誤差遞增明顯,SRTMGL1則隨著地形起伏的變化,差值變化幅度很小。3種DEM在臺地類型區(qū)域數(shù)據(jù)精度均低于平原和丘陵。ASTER V2隨著地貌變化(在山脈類型區(qū)域),精度變化幅度較SRTM V4更加平穩(wěn),數(shù)據(jù)精度更高。

本文對SRTMGL1,SRTM V4和ASTER GDEM V2的精度進行了詳細的對比分析,可以為這3種DEM數(shù)據(jù)在不同類型區(qū)域的應(yīng)用提供一定的指導(dǎo)作用。數(shù)據(jù)選取的總體原則為:SRTMGL1最優(yōu)選,其次,平原為主地區(qū)選擇SRTM V4,高原以及植被覆蓋度高的地形復(fù)雜區(qū)域盡量選擇ASTER GDEM V2.