利用雷達(dá)攝影測(cè)量方法提取DEM及其精度評(píng)價(jià)

1、利用雷達(dá)攝影測(cè)量方法提取DEM 及其精度評(píng)價(jià)鄧?yán)? 陳悅竹, 陳云浩, 胡德勇, 李京(北京師范大學(xué)資源學(xué)院, 資源技術(shù)與工程研究所, 北京100875; 中華測(cè)繪技術(shù)服務(wù)公司, 北京,100088摘要:為了解雷達(dá)立體攝影測(cè)量中各項(xiàng)因素對(duì)最終DEM 精度的影響, 采用馬來(lái)西亞熱帶雨林地區(qū)具有不同波束模式和入射角的6對(duì)Radarsat -1影像, 應(yīng)用兩種不同的SAR 成像模型, 即距離/多普勒模型和等效共線方程模型, 對(duì)使用雷達(dá)攝影測(cè)量方法提取DEM 進(jìn)行了試驗(yàn)。首先分析了兩種成像模型不同的物理基礎(chǔ), 然后通過(guò)比較由它們得到的DEM 的精度, 發(fā)現(xiàn)利用距離/多普勒模型提取DEM 的精度優(yōu)于等效

2、共線方程模型; 然后, 利用研究區(qū)的數(shù)字地形圖等參考資料, 分析了不同軌道、不同模式、不同分辨率、不同交角的立體像對(duì)組合以及地形因素對(duì)DEM 精度的影響, 要得到高精度的DEM , 必須綜合考慮立體像對(duì)的選取和研究區(qū)的地形、地貌等特征。關(guān)鍵詞:雷達(dá)攝影測(cè)量;Radarsat ; 合成孔徑雷達(dá)圖像;DEM ; 精度評(píng)價(jià)中圖分類(lèi)號(hào):P237. 9文獻(xiàn)標(biāo)識(shí)碼:A 文章編號(hào):1000-3177(2006 85-0037-041引言DEM (Digital Elevation Models 是地表形態(tài)的數(shù)字形式, 它由規(guī)則水平間隔處地面點(diǎn)的抽樣高程矩陣組成, 具有廣泛的應(yīng)用潛力1。隨著遙感技術(shù)的發(fā)展,

3、像提取DEM 已經(jīng)得到廣泛應(yīng)用。如利用(Spaceborne Thermal Emission , 65m 2, 利用Radarsat SAR 數(shù)生成精度約2060m 的DEM 3, 利用SPO T 數(shù)據(jù)提取的DEM 精度可達(dá)10m 以?xún)?nèi)4, 而利用InSAR 生成的DEM 精度可達(dá)米級(jí)5。SAR 是一種工作在微波波段的主動(dòng)式傳感器, 能夠不受日照和天氣的限制, 全天候、全天時(shí)地對(duì)地觀測(cè), 同時(shí)由于衛(wèi)星影像成像范圍大, 可重復(fù)對(duì)地觀測(cè), 因此SAR 衛(wèi)星影像的應(yīng)用越來(lái)越受重視。目前利用SAR 傳感器進(jìn)行DEM 提取主要分為三種路線69:雷達(dá)攝影測(cè)量技術(shù)(StereoSAR 和Radarclin

4、ometry , 干涉測(cè)量技術(shù)(InSAR 和極化測(cè)量技術(shù)(Polarimetry 。盡管InSAR 在精度方面比StereoSAR 具有絕對(duì)優(yōu)勢(shì)10, 但是綜合科技水平、經(jīng)濟(jì)實(shí)力、應(yīng)用效益等條件,InSAR 技術(shù)對(duì)許多國(guó)家來(lái)說(shuō)仍難以實(shí)現(xiàn), 而且StereoSAR 與InSAR 兩種技術(shù)之間存在相互補(bǔ)充和融合的領(lǐng)域, 比如, 采用InSAR 獲取DEM 時(shí), 由于大氣相對(duì)濕度變化導(dǎo)致干涉相位扭曲, 得到的DEM , 通過(guò)將InSAR 和StereoSAR , InSAR 數(shù)據(jù)中-1衛(wèi)星的SAR 傳感、異軌獲取立體影像的特點(diǎn), 從SAR 影像變得簡(jiǎn)單。本文使用熱6對(duì)Radarsat 影像構(gòu)成的立

5、體像對(duì), 利用兩種基于不同原理的SAR 圖像成像模型, 進(jìn)行立體SAR 影像提取DEM 方法的研究。通過(guò)試驗(yàn), 將兩種成像模型提取的DEM 的精度進(jìn)行了比較, 同時(shí)分析了不同的立體像對(duì)組合對(duì)DEM 精度的影響。2利用雷達(dá)攝影測(cè)量方法提取DEM立體像對(duì)是由不同攝站獲取的同一地區(qū)的具有一定重疊度的兩張像片, 雷達(dá)影像存在高差引起的投影差, 同一高程的地物在兩張雷達(dá)影像上投影差不同, 當(dāng)用兩像片進(jìn)行立體觀察時(shí)將引起視差, 產(chǎn)生立體感。雷達(dá)立體成像的方式有兩種, 即在目標(biāo)同側(cè)成像構(gòu)成立體像對(duì)和在目標(biāo)異側(cè)成像構(gòu)成立體像對(duì), 如圖1所示。異側(cè)獲取的立體像對(duì), 圖像幾何變形差異較大, 而且影像色調(diào)差異也大,

6、 給立體模型的建立和影像匹配帶來(lái)困難。而同側(cè)獲取的立體像對(duì), 兩景影像上的變形基本相同, 透視收縮、疊掩和陰影在兩景影像出現(xiàn)的情形大致相同, 引起的色調(diào)差異也較小, 因此, 利用雷達(dá)立體圖1立體像對(duì)的構(gòu)成方式11收稿日期:2005-12-07修訂日期:2006-03-25基金項(xiàng)目:高等學(xué)校博士點(diǎn)專(zhuān)項(xiàng)科研基金(20030027014 , 三維信息獲取與應(yīng)用教育部重點(diǎn)實(shí)驗(yàn)室開(kāi)放基金聯(lián)合資助。作者簡(jiǎn)介:鄧?yán)?1976 , 博士生, 主要從事遙感圖像處理研究. E 2m ail :DengL ires. cn732006. 4應(yīng)用技術(shù) 遙感信息影像進(jìn)行立體測(cè)量多采用同側(cè)方式得到的像對(duì)。在提取DEM 的

7、過(guò)程中, 雷達(dá)圖像成像模型的選擇對(duì)最終結(jié)果有很大的影響。本研究中, 分別采用了基于距離/多普勒SAR 成像模型和基于等效共線方程的模型1213。2. 1距離/多普勒SAR 成像模型距離/多普勒成像模型是根據(jù)雷達(dá)圖像像點(diǎn)的距離條件和多普勒條件表達(dá)像點(diǎn)、物點(diǎn)和雷達(dá)成像參數(shù)之間的關(guān)系。對(duì)于地距圖像, 其距離條件公式為:(X -X s 2+(Y -Y s 2+(Z -Z s 2=R 2(1式中, (X , Y , Z 為地面點(diǎn)坐標(biāo), (Xs , Ys , Zs 為地面點(diǎn)對(duì)應(yīng)的雷達(dá)天線的物方空間坐標(biāo), R 是要量測(cè)的點(diǎn)到傳感器的斜距。對(duì)于方位向, 采用零多普勒條件, 即:當(dāng)衛(wèi)星飛行速度矢量與天線至地面點(diǎn)

8、距離矢量間的角度為90°時(shí), 多普勒頻率為零, 如公式(2 所示:Vsx (X -X s +Vsy (Y -Y s +V S Z (Z -Z s =0(2 Vsx Vsy Vsz 為衛(wèi)星飛行速度矢量的三個(gè)分量。根據(jù)距離/多普勒構(gòu)像方程, 左片和右片可以各列2個(gè)方程:(X -X s 1 2+(Y -Y s 1 2+(Z -Z s 1 2-R 21=0V sx 1(X -X s 1 +V sy 1(Y -Y s 1 +V sz 1(Z -Z s 1 =0(X -X s 2 2+(Y -Y s 2 2+(Z -Z s 2 2-R 22=0V sx 2(X -X s 2 +V sy 2(Y

9、-Y s 2 +V sz 2(Z -Z s 2 = 0(3 X s 1Y s 1Z s 1V sx 1V sy 1V sz 1是左片的衛(wèi)星位置矢量和速度矢量; X s 2Y s 2Z s 2V sx 2V sy 2V sz 2是右片的衛(wèi)星位置矢量和速度矢量; 求解方程組(3 , 即可以求得地面點(diǎn)的坐標(biāo)(X , Y , Z 。2. 2等效共線方程模型SAR 圖像是一種距離投影圖像, 它與可見(jiàn)光圖像的中心投影有很大的差別。K onecny 于1988年首先綜合分析了兩種投影之間的不同, 并提出了利用等效共線方程來(lái)描述SAR 成像模型。對(duì)于雷達(dá)成像, 其等效共線方程為:x gr =0=-f x( (

10、 ( c 1j (X i -X -X sj +c 2j (Y i -Y -Y sj c -Z sj y gr =0=-f y ( c 1j (X i -X -X sj (Y -sj c 3j (Z -Z sj (4X i 、Y i 、Z i 是地面i 點(diǎn)坐標(biāo), X sj 、Y sj 、Z sj 在j 時(shí)刻的位置, f x 、f y 是等效焦距, x i 點(diǎn)坐標(biāo), a 1j 、a 2j 、c 3j 、組成的方向余弦, SAR 圖像, 需要求解6(X sj , Y sj , Z sj , , , 。3實(shí)驗(yàn)結(jié)果本文選擇馬來(lái)西半島中部的Raub 地區(qū)進(jìn)行DEM 提取試驗(yàn)并進(jìn)行精度驗(yàn)證。該研究區(qū)的地理位

11、置在N 3°16-3°47,E 101°37-102°14之間, 覆蓋近3000km 2的區(qū)域, 最高峰2048m 。全境處于熱帶氣候區(qū), 植被覆蓋情況復(fù)雜, 既有大片茂密的原始熱帶雨林, 又有大量的人工種植的油棕和橡膠樹(shù)等植物。該研究區(qū)地勢(shì)起伏顯著, 有大片的平原區(qū), 同時(shí)擁有一系列不同形狀、不同大小、不同高低的山脈集合而成的山地地貌。本次研究使用了研究區(qū)的6對(duì)Radarsat 數(shù)據(jù), 其詳細(xì)情況如表1所示。精度驗(yàn)證的數(shù)據(jù)有:該地區(qū)1:50,000地形圖和由等高線提取的圖像分辨率為30m 的DEM 。不同的SAR 成像模型對(duì)地面控制點(diǎn)(GCP 的數(shù)目要

12、求有所不同。而GCP 的采集, 一般來(lái)說(shuō)費(fèi)時(shí)費(fèi)力, 而且在SAR圖像上定位GCP 有一定的難度, 因此, 對(duì)GCP 的依賴(lài)越少越好。當(dāng)使用距離/多普勒模型時(shí), 由于Radarsat 的軌道數(shù)據(jù)相對(duì)比較穩(wěn)定, 甚至可以不用地面控制點(diǎn), 完全依靠軌道數(shù)據(jù)進(jìn)行SAR 影像成像模型的建立。但為了提高結(jié)果的精度, 在本研究中使用Trimble 4000SE 進(jìn)行野外差分GPS 測(cè)量, 得到高精度的地面控制點(diǎn)(GCP , 利用這些控制點(diǎn)對(duì)衛(wèi)星軌道參數(shù)進(jìn)行糾正。對(duì)于等效共線方程模型而言, 為有效, 一般要求至少有6個(gè)均勻分布的GCP , 但通過(guò)試驗(yàn)我們發(fā)現(xiàn), 一般取多于20個(gè)均勻分布的GCP 才能獲得比較好

13、的效果。表2列出了每個(gè)立體像對(duì)組合使用的GCP 數(shù)目。表1研究區(qū)內(nèi)的Radarsat 數(shù)據(jù)NO. 波束軌道入射角(. 分辨率(m (距離×方位 獲取時(shí)間1S6升軌454920×272005-06-16W2升軌31-3925×272005-06-062S230-3720272005-02-08S7升軌39-4219×272005-02 -103S3降軌30-3725×272005-02-10S5降軌36-4222×272005-02-174S5升軌36-4222×272005-02-18S7升軌39-4219×272

14、005-02-115S1降軌20-2724×272005-02-27S3降軌30-3725×272005-02-106F1降軌37-418. 3×8. 42005-03-13F4降軌43-467. 3×8. 42005-02-24表2不同的立體像對(duì)使用的地面控制點(diǎn)(GCP 數(shù)目S1/S3S3/S5S5/S7S2/S7W2/S6F1/F4距離/多普勒566764等效共線方程263331261714試驗(yàn)提取的DEM 像素分辨率均為30m 。插頁(yè)2圖1顯示了由W2/S6像對(duì)、S2/S7像對(duì)組合得到的DEM 和相應(yīng)的地貌圖。4精度評(píng)價(jià)與分析在實(shí)際應(yīng)用中, 常用的

15、DEM 精度評(píng)價(jià)方法有檢查點(diǎn)法、83遙感信息應(yīng)用技術(shù)2006. 4剖面法、等高線回放法等14。本文主要采用檢查點(diǎn)法和剖面法。按照地理坐標(biāo), 在6個(gè)DEM 中, 隨機(jī)選取大小相同的6塊進(jìn)行驗(yàn)證, 大小為26km ×43km 。結(jié)合DEM 和地形圖取得該區(qū)域內(nèi)270個(gè)點(diǎn)的高程值作為檢查點(diǎn), 使用總均方根差(RMSE 對(duì)DEM 進(jìn)行精度檢驗(yàn)。其定義為:RMS E =n(y i -y i 2n (4 式中, y i 和y i 分別為測(cè)定值和預(yù)測(cè)值, n 為樣本數(shù)。表3列出了使用檢查點(diǎn)對(duì)DEM 進(jìn)行精度檢驗(yàn)的統(tǒng)計(jì)結(jié)果。表3DEM 精度檢驗(yàn)結(jié)果(單位:m 像對(duì)交角RMSE 均值標(biāo)準(zhǔn)差等效共線方程

16、W2/S610. 041. 251. 2541. 23S2/S719. 236. 9714. 3134. 08S1/S310. 855. 85-16. 1353. 47S3/S55. 178. 3714. 5677. 00S5/S77. 749. 0717. 6645. 78F1/F45. 846. 20-5. 9345. 82距離/多普勒W2/S610. 024. 250. 8124. 23S2/S719. 224. 589. 1122. 83S1/S310. 830. 426. 3829. 74S3/S55. 149. 0914. 4446. 92S5/S77. 739. 490. 503

17、9. 49F1/F45. 825. 41-0. 0325. 41地形圖N/A 8. 206. 98觀察表3, 對(duì)于距離/, 和S7組合得到的DEM , 其分別為24. 25m 和24. 58m 。,S2/S7得到的DEM 的精度最高, 其RMSE 為36. 97m 。圖2顯示了使用不同SAR 成像模型和不同的像對(duì)組合得到的DEM 的RMSE 。圖2不同像對(duì)使用不同成像模型得到的DEM 的精度比較為了進(jìn)一步分析不同的地形對(duì)DEM 的精度的影響, 我們利用剖面法在使用不同SAR 成像模型提取的DEM 上任意做4條剖面, 并計(jì)算其RMSE 。結(jié)果見(jiàn)插頁(yè)2圖2, 限于篇幅, 只列出了S2/S7組合的剖

18、面圖。綜合分析表3、圖2及插頁(yè)2圖2, 我們發(fā)現(xiàn):(1 對(duì)于各種立體像對(duì)組合, 由距離/多普勒模型提取的DEM 的精度均高于由等效共線方程模型提取的DEM 的精度;(2 由插頁(yè)2圖2所顯示的剖面曲線可以看出, 在地形比較平緩的地區(qū), 距離/多普勒模型提取的DEM 的精度高于等效共線方程模型(如插頁(yè)2圖2(a (b 所示 ; 而在地形陡峭的地區(qū), 由等效共線方程模型提取的DEM 的精度要高于距離/多普勒模型(如插頁(yè)2圖2(c (d 所示 ; (3 在各種地形上, 由R -D 模型提取的DEM 的剖面曲線更加接近于由等高線生成的DEM 的剖面曲線。通過(guò)以上分析我們可以發(fā)現(xiàn), 由距離/多普勒模型提取

19、的DEM 的精度要優(yōu)于由等效共線方程提取的DEM 的精度。下面我們將進(jìn)一步分析地形與不同立體像對(duì)之間相互作用的關(guān)系并由此判斷不同的立體像對(duì)提取的DEM 的優(yōu)劣。插頁(yè)2圖3顯示了6個(gè)DEM 上任意4條剖面曲線, 圖3顯示了地形因素對(duì)提取的DEM 精度的影響。為了分析的方便起見(jiàn), 這里使用的DEM 都是由距離/多普勒模型提取的 。DEM 精度的影響23, 插頁(yè)2圖3, 可以發(fā)現(xiàn):(, 得到的DEM 的精度隨之下降, ;(2 地形因素對(duì)DEM 的精度有極大的影響, 由圖5可以看出:隨著地形由平緩變得陡峭,DEM 的精度基本是線性下降的; 同時(shí), 地形因素對(duì)于交角大的立體像對(duì)的影響更大, 比如從圖5中

20、看到,S2/S7在平坦地區(qū), 其精度高于大部分其他立體像對(duì), 但是隨著地形的變陡, 其精度下降的速度遠(yuǎn)遠(yuǎn)高于其他立體像對(duì)組合;(3 在地形相對(duì)平緩地區(qū), 交角大的立體像對(duì)組合提取的DEM 精度比較高, 而在地形陡峭地區(qū), 交角小的組合則相對(duì)較優(yōu)。比如插頁(yè)2圖3(a 顯示, 在高山地區(qū), S2/S7的精度較低, 而其他交角相對(duì)較小的組合(如W2/S6, S1/S3 的DEM 精度高于S2/S7; 而在相對(duì)平緩地區(qū), 如插頁(yè)2圖3(b 、(d 所示,S2/S7提取的DEM 有不錯(cuò)的精度;(4 影像分辨率的提高, 并不能大幅度的提高最終DEM 的精度。比如,F1/F4組合提取的DEM 的精度與S2/

21、S7組合的DEM 具有類(lèi)似的精度, 而由表1可以看到, 精細(xì)模式的空間分辨率是標(biāo)準(zhǔn)模式的3倍; (5 從插頁(yè)2圖3還可以發(fā)現(xiàn), 所有組合提取的DEM 均比由地形圖得到的DEM 高。這是由于Radarsat -1采用的是C 波段的SAR , 而C 波段無(wú)法穿透濃密的植被, 因此, 由該方法得到的DEM 包含植被的高度。通過(guò)以上分析我們可以發(fā)現(xiàn), 由雷達(dá)攝影測(cè)量的方法提取DEM 的精度受到多種因素的影響, 要得到高精度的DEM , 必須綜合考慮立體像對(duì)的選取和研究區(qū)的地形、地貌等特征。932006. 4應(yīng)用技術(shù)遙感信息5結(jié)論與討論隨著遙感技術(shù)的發(fā)展, 越來(lái)越多載有SAR 傳感器的衛(wèi)星進(jìn)入軌道, 從

22、而使得利用衛(wèi)星獲取立體SAR 影像變得簡(jiǎn)單, 數(shù)據(jù)源非常豐富, 為解決那些常年多云霧、多陰雨和植被覆蓋茂密, 光學(xué)和光電成像困難地區(qū)的測(cè)繪難題提供了一種新的途徑。加拿大的Radarsat -1衛(wèi)星的SAR 傳感器具有波束模式入射角可調(diào)、異軌獲取立體影像的特點(diǎn), 本文通過(guò)使用同一地區(qū)的具有不同交角的6對(duì)Radarsat 影像, 使用2種不同的SAR成像模型, 對(duì)利用雷達(dá)攝影測(cè)量方法提取DEM 的方法進(jìn)行了研究, 并將結(jié)果與地形圖等資料進(jìn)行了比較, 結(jié)果表明, 由該方法獲取的DEM 的垂直精度在25m 左右, 可以滿(mǎn)足一般應(yīng)用的需要。同時(shí)發(fā)現(xiàn), 距離/多普勒立體模型基于SAR 成像的物理原理, 其

23、構(gòu)成簡(jiǎn)單, 對(duì)地面控制點(diǎn)要求不高, 并且得到的DEM 的精度高于等效共線方程模型; 地形因素和立體像對(duì)的交角共同影響著最終DEM 的精度, 大的交角適合于地勢(shì)平緩地區(qū), 而對(duì)于陡峭地區(qū), 則應(yīng)選擇具有較小交角的立體像對(duì)組合; 影像分辨率的提高并不能大幅度提高DEM 的精度, 而且, 隨著影像分辨率的提高, 影像中斑點(diǎn)噪聲的影響也隨之增強(qiáng), 增加了處理的難度; 并且, 由于傳感器自身的限制, 生成的DEM 含有系統(tǒng)誤差, 即包含植被的高度。致謝:感謝馬來(lái)西亞國(guó)家遙感中心(Malaysia Centre for Remote Sensing :MACRES 為本次試驗(yàn)提供Radarsat 衛(wèi)星影像

24、以及參考數(shù)據(jù)。參考文獻(xiàn)1李志林, 朱慶. 數(shù)字高程模型M .武漢:武漢大學(xué)出版社, 2001.2惠鳳鳴, 田慶久, 李英成. Aster 數(shù)據(jù)的DEM 生產(chǎn)及精度評(píng)價(jià). 遙感信息, 2004(1 :1418.3Thiery Toutin. Evaluating Radarsat stereoscopic pairs for DEM generation J/OL .Proc. IEEE IGRASS 98, 2:10341036.4黃玉琪. SPO T 影像的DEM 自動(dòng)生成J.測(cè)繪通報(bào), 1998(9 :1316.5張永紅, 張繼賢, 林宗堅(jiān). 由星載INSAR 生成DEM 的理論誤差分析J

25、2 :6Radargrammetry. http :/155. 207. 39. 4/orfeas/radargrammetry. html.7Sophie Paquerault , Henri Maitre , Jean 2Marie. Radarclinometry 1.IG ARSS 96, 1996(1 :503505.8Thierry Toutin , Laurence Gray. State 2of 2the 2SAR dataJ.ISPRS. 2000, 55(1 , 1933.9Schler , Jong 2Sen Lee , G ianfranco topography us

26、ing polarimetric SAR imagesJ.IEEE Trans. On G eoscience And remote sensing , 1996, (126610Anselm Schmieder , Automatic derivation of generalized contour lines for topographic maps using high 2resolution airborne inter 2ferometric radar data C.International Archives of Photogrammetry and Remote Sensi

27、ng , 2000, XXXII , Part B1.11劉明軍. 利用Radarsat 立體影像提取DEM 攝影測(cè)量方法研究D.武漢:武漢大學(xué). 2004.12Leberl , F. W. Radargrammetric image processingM .Artech House , Norwood , 1990.13G. K onecny , W. Schuhr. Reliability of radar image data C.International Archives of Photogrammetry and Remote Sensing , 1988, . XXXII ,P

28、art B2.14唐新明, 林宗堅(jiān), 吳嵐. 基于等高線和高程點(diǎn)建立DEM 的精度評(píng)價(jià)方法探討J.遙感信息, 1999(3 :710DEM G eneration and Accuracy Assessment Using R adargrammetryDEN G Lei , CHEN Yue 2zhu , CHEN Yun 2hao , HU De 2yong , L I Jing (College of Resources Science and Technology , Beijing Norm al U niversity , Beijing 100875;Technology Service of S urveying &M apping Co. of China , Beijing 100088Abstract :To understand the relationship between DEM accuracy and some major factors in radargrammetry , six pairs of Radarsat 21image with different beam