□ 夏立福 顧劍華 王雷 孫鑫

（1.黑龍江第三測繪工程院，黑龍江哈爾濱150025；2.黑龍江第一測繪工程院，黑龍江哈爾濱150025）

1.引言

為保持國家基礎1∶50000地形圖的現(xiàn)勢性，保障國民經(jīng)濟部門對基礎地形圖的需求，國家測繪地理信息局決定在“十二五”期間對“十一五”期間生產(chǎn)的1∶50000地形數(shù)據(jù)庫進行動態(tài)更新。為使動態(tài)更新工作快捷、有序進行，作者根據(jù)自身生產(chǎn)實際，選取天津、黑龍江區(qū)域的四景影像通過不同的方法進行DOM產(chǎn)品的生產(chǎn)試驗，尋找資源三號衛(wèi)星影像糾正的最佳方法，以確保達到精度要求，為要素變化發(fā)現(xiàn)提供依據(jù)。

2.DOM糾正技術流程與DEM選取

DOM糾正實驗選取四景影像，充分考慮地形類別及已有基礎地理信息數(shù)據(jù)，其分布情況及時相如下：天津Z389512725P20121031地形類別為平地覆蓋在東麗區(qū)、津南區(qū)及濱海新區(qū)；黑龍江877109_20120404地形類別為平地丘陵結(jié)合區(qū)域，876110_20120612為山地，875112_20120419為丘陵，主要覆蓋區(qū)域在哈爾濱市周邊五常市、呼蘭縣、阿城、賓縣等。影像均無云雪覆蓋，信息豐富滿足影像制作要求，且時相分別為2012年10月31日、2012年4月4日、2012年6月12日、2012年 4月19日。根據(jù)已有“十一五”期間生產(chǎn)的1∶50000地形數(shù)據(jù)庫3D產(chǎn)品為依托，分別利用 RPC+DEM、RPC+DEM+一次多項式、RPC+DEM+二次多項式模型進行糾正精度的對比。經(jīng)分析研究利用RPC+DEM+一次多項式的糾正方法精度較好，滿足1∶50000地形數(shù)據(jù)庫動態(tài)更新的需求，其DOM糾正具體技術流程如下：

DEM的格網(wǎng)間距理論上應與地形坡度成反比。在平原地區(qū)，地面坡度平緩，DEM數(shù)據(jù)可以采用較大的格網(wǎng)間距；而在地形破碎的丘陵地區(qū)和山區(qū)，地形陡峭，DEM格網(wǎng)間距需要采用較小的格網(wǎng)間距，才能保證地形地貌的特征與實地相符。本次實驗DEM格網(wǎng)間距為25米，滿足不同地形類別影像的糾正要求。為提高影像質(zhì)量精度，在地形復雜、高差大的山區(qū)，需考慮衛(wèi)星數(shù)據(jù)的成像的側(cè)視角。

3.DOM成果的精度評定和分析

RPC+DEM+一次多項式的糾正方法即滿足1∶50000地形數(shù)據(jù)庫動態(tài)更新的需求，可從其DOM成果的精度檢測中體現(xiàn)。DOM成果的精度評定是依據(jù)1∶50000DLG在DOM影像上選取檢查點位置并讀取坐標值，計算其較差的中誤差。檢查點分布均勻且在影像上能準確定位。本次實驗分別選取研究區(qū)的每景25個及以上特征地物點作為檢查點，所有檢查點均勻分布在整個研究區(qū)，是影像和DLG上都容易辨認準確定位的點。其精度統(tǒng)計包括同精度檢測（內(nèi)業(yè)精度檢測）及高精度檢測（外業(yè)實測檢測）。

同精度檢測

高精度檢測

其中，RMS為檢查點較差中誤差，n為檢查點個數(shù)，xi，yi表示 DOM影像上檢查點的坐標，Xi，Yi為 GPS外業(yè)檢查點或圖形檢查點的坐標。

本實驗內(nèi)業(yè)檢測點及外業(yè)檢測點分布如下（圓形為內(nèi)業(yè)檢查點；三角形為外業(yè)檢查點）：

圖1 895127 內(nèi)、外業(yè)檢查點分布圖

圖2 877109內(nèi)、外業(yè)檢查點分布圖

圖3 876110內(nèi)、外業(yè)檢查點分布圖

圖4 875112內(nèi)、外業(yè)檢查點分布圖

3.1 內(nèi)業(yè)相對精度分析

內(nèi)業(yè)檢測精度是1∶50000數(shù)字正射影像地物點相對于附近1∶50000地形數(shù)據(jù)的點位中誤差。

精度統(tǒng)計公式如下：

其中，RMS為檢查點較差中誤差，n為檢查點個數(shù)，xi，yi表示DOM影像上檢查點的坐標，Xi，Yi為GPS外業(yè)檢查點或圖形檢查點的坐標。

按照同精度檢測公式分別統(tǒng)計四景影像的內(nèi)符合精度：

表1 內(nèi)符合精度統(tǒng)計表

3.2 外業(yè)實測精度分析

外業(yè)實測精度是：1∶50000數(shù)字正射影像地物點相對于附近野外控制點的點位中誤差不得大于表1規(guī)定：

表2 2.5米分辨率數(shù)字正射影像平面中誤差

精度統(tǒng)計公式如下：

其中，RMS為檢查點較差中誤差，n為檢查點個數(shù)，xi，yi表示DOM影像上檢查點的坐標，Xi，Yi為GPS外業(yè)檢查點或圖形檢查點的坐標。

按照高精度檢測公式分別統(tǒng)計四景影像的外符合精度：

表3 外符合精度統(tǒng)計表

3.3 接邊精度

接邊精度的檢查：取相鄰兩景正射影像重疊區(qū)域處同名點作為檢查點，分別量取兩同名點的距離，或者是讀取同名點的坐標，算出兩點間的距離，檢查同名點的距離或坐標差是否符合限差。本次試驗所選取的四景影像在與相鄰影像間經(jīng)過接邊精度的檢查，坐標差均在限差范圍之內(nèi)，故符合接邊精度的要求。

4.結(jié)論

此次實驗選取的四景資源三號衛(wèi)星影像，地形類別包括平地、丘陵、山地?？刂泣c均勻分布選擇在有線狀地物交叉部位，盡可能選正十字或丁字交叉口。經(jīng)過內(nèi)業(yè)相對精度分析和外業(yè)實測像控點精度分析（詳見衛(wèi)星影像高精度、同精度檢測統(tǒng)計表），四景影像均滿足精度要求，故此次利用ERDAS的RPC+DEM+一次多項式的糾正方法適用且滿足1∶5萬地形數(shù)據(jù)庫動態(tài)更新工程正射影像糾正的精度要求。

【1】馬治，熊康軍，陶源，王鐵軍.基于資源三號衛(wèi)星影像的1∶50000數(shù)據(jù)生產(chǎn)試驗[J].測繪信息與工程.2012，35（11）：214-216.

【2】景耀全，陳笑峰，諶佳.DEM精度對高分辨率衛(wèi)星影像糾正的影響[J].測繪.2011，34（2）：70-73.

【3】張過，李德仁，袁修孝，張春玲.衛(wèi)星遙感影像的區(qū)域網(wǎng)平差成圖精度[J].測繪科學技術學報.2006，23（4）：239-241.

【4】楊金香，程學豐.基于ERDAS IMAGINE操作平臺的遙感圖像處理[J].測繪與空間地理信息.2006，29（2）：107-110.

【5】馮秀麗，王珂，樓立明.基于ERDAS的遙感影像正射圖的制作[J].遙感技術與應用.2003，18（3）：176-179.