劉志鵬 王 義
(中鐵二院工程集團有限責任公司,成都 610031)
資源三號衛(wèi)星是中國第一顆自主發(fā)射的民用高分辨率立體測繪衛(wèi)星,該衛(wèi)星攝影范圍較寬,單景影像覆蓋面積達50km×50km。資源三號衛(wèi)星歷經(jīng)兩代,其中,其前后視立體影像分辨率由3.5m提升到2.5m,中視空間分辨率為2.1m。通過資源三號衛(wèi)星影像的逐步推廣應用,基本實現(xiàn)了國產(chǎn)遙感衛(wèi)星從“有”到“好用”,從示范應用到業(yè)務化運行的根本轉(zhuǎn)變[1]。然而,將資源三號衛(wèi)星影像應用于復雜艱險山區(qū)鐵路制圖的研究相對較少。
另一方面,我國西南地區(qū)地形復雜,氣象條件較差,常年云霧彌漫,航空攝影極其困難。國外高分辨率衛(wèi)星影像雖然能滿足大、中比例地形圖制作,但存檔數(shù)據(jù)較少,覆蓋長大鐵路干線項目的范圍有限。為了解決該區(qū)域攝影資料缺乏等問題,利用多景資源三號02號衛(wèi)星影像,研究在無像控條件下的數(shù)字正射影像圖的快速制作方法,并驗證布設少量控制點制作1∶10000地形圖精度的可靠性。
有理函數(shù)模型簡稱RFM(Rational Function Model),是通用傳感器模型中的一種,它的實質(zhì)是通過復雜的數(shù)學模型來描述嚴格成像過程,建立物方空間坐標與像方像點坐標的直接對應關(guān)系。RFM可以有效針對衛(wèi)星軌道參數(shù)、傳感器屬性等涉密參數(shù)進行封裝,用戶可在不知道相關(guān)信息的前提下,直接根據(jù)RFM模型及參數(shù)解算出像點的空間坐標,從而進行各種應用。RFM 一般采用 RPC參數(shù)(Rational Polynomial Coefficients,有理多項式系數(shù))來表示。
式中,a0…d19為參數(shù)文件中提供的有理多項式參數(shù);xn,yn為圖像上對應點的規(guī)格化坐標[2]。
RPC文件中詳細說明了影像中心點地理坐標,平面坐標用經(jīng)緯度表示,高程為大地高,影像的行列數(shù)以像素為單位,公式中羅列了物方與像方坐標轉(zhuǎn)換的80個有理函數(shù)系數(shù)。通過定向或者空三加密即可求出任意像點的大地坐標[3]。
資源三號衛(wèi)星影像RPC參數(shù)文件格式與常用的Ikonos、Worldview 、Geo-eye、Irs-P5、Quickbird等衛(wèi)星影像一致[4],預處理及空三加密方法與這幾種衛(wèi)星影像相同。具體處理流程見圖1。
圖1 資源三號衛(wèi)星影像處理流程
資源三號衛(wèi)星影像數(shù)據(jù)以16位存儲,影像數(shù)據(jù)量較大,灰度較暗。必須利用PhotoShop、Eardas等軟件進行灰度拉伸,以提高亮度和對比度,增強人眼辨識度。并利用Rescale模塊對影像進行降位處理,從16位降位到8 位,降位后幾何精度不變,但影像數(shù)據(jù)量減少1/2,以方便快速打開影像,實施漫游測圖[5]。
資源三號衛(wèi)星影像為同軌成像模式,為恢復立體模型,一般需對影像及RPC參數(shù)作90°旋轉(zhuǎn)。可利用自主開發(fā)軟件實施影像及RPC參數(shù)旋轉(zhuǎn),當RPC參數(shù)與影像完全匹配時,可進行立體測圖。
衛(wèi)星影像RPC參數(shù)是對衛(wèi)星影像的粗略定位,若定位精度不高,必須利用空三加密、模型連接、像控點控制,區(qū)域網(wǎng)平差等方法提高定位精度,從而達到對RPC參數(shù)重新糾正的目的[6]。
針對四川某鐵路項目,購買了3景資源三號衛(wèi)星影像,開展衛(wèi)星影像區(qū)域網(wǎng)空三加密及制圖研究工作,首先利用ERDAS的IMAGINE Photogrammetry軟件對資源三號衛(wèi)星影像作空三處理,衛(wèi)星模型選擇ZY-3 RPC,設置測圖坐標系、中央子午線及投影方式。
然后,加載衛(wèi)星影像,選擇與影像匹配的RPC地理定位數(shù)據(jù),通過Compute Pyramids and Statistics模塊批量建立金字塔影像,并設置匹配策略,匹配點個數(shù)設置為整幅影像點的1/10,立體像對自動匹配連接點,匹配完成后,剔除大于1個像素的粗差點,并補充適當人工點,調(diào)整上下視差,使其符合限差要求。調(diào)整后,像點殘差精度報告見表1。
表1 像點坐標精度報告 像素
由表1可知,在全自動匹配情況下,影像匹配中誤差為0.4個像素,最大誤差為0.7個像素,匹配效果較好。證明前、后兩視角影像姿態(tài)都較穩(wěn)定,可以用于后續(xù)測圖工作。利用Geometric Calibration下的 RPC Generation模塊糾正RPC參數(shù),輸出準確的RPC參數(shù),從而完成自由網(wǎng)平差[7]。
基于RPC參數(shù)進行自由網(wǎng)平差解算,平面坐標為地理坐標系,格式為經(jīng)緯度,測圖時需對該坐標進行投影轉(zhuǎn)換[8];其高程為大地高,需通過EGM2008高程模型對測區(qū)作高程異常改正,將高程納入到正常高系統(tǒng),從而提高無控測圖精度[9]。
通過自由網(wǎng)平差后,資源三號衛(wèi)星影像可以制作成數(shù)字正射影像刺點及調(diào)繪片,數(shù)字正射影像附有坐標,以外業(yè)導航尋路,從而提高作業(yè)效率。另外,用于刺點及調(diào)繪的數(shù)字正射影像精度要求不高,可以下載SRTM 1-ARC Second 30m格網(wǎng)間距的DEM數(shù)據(jù)作高程糾正,平面基于有理函數(shù)模型糾正,即利用RPC參數(shù)文件進行平面定位[10]。
資源三號衛(wèi)星影像全色中視片空間分辨率為2.1m,分辨率優(yōu)于前、后視,正射影像制作優(yōu)先選擇該視角影像進行。為提高運行速度,實現(xiàn)多線程多任務,可運用Batch工具實現(xiàn)多景影像批處理(可同時運行4個任務)。
基于糾正好的全色正射影像對多光譜影像進行自動配準,自動配準過程中,自動匹配同名點,計算誤差并剔除粗差點,對多光譜影像進行正射重采樣,使用的模塊為IMAGINE Photogrammetry下的AutoSync工具。
勻光勻色有兩個目的:①使整個測區(qū)內(nèi)色彩基本一致;②使每張影像內(nèi)部色調(diào)一致。具體操作時,應選擇色調(diào)具有代表性的清晰影像作為基準影像,對測區(qū)所有影像進行勻色處理,從而達到整個測區(qū)影像色彩一致,并利用軟件作輻射校正和去除影像中的輕霧,對自動效果不滿意的區(qū)域,可進行局部人工調(diào)整。
多張正射影像制作完成后,通過mosaic模塊完成影像鑲嵌,設置自動選擇拼接線,以及拼接線的羽化寬度,保證拼接處影像灰度平滑過渡自然。鑲嵌完成后,對正射影像進行分幅裁剪,輸出GeoTiff影像,再利用具有坐標格式的影像開展數(shù)字刺點和調(diào)繪工作[11]。
利用前期完成的自由網(wǎng)平差工程,3景影像中平均每景影像各量測了5個平高控制點(四角分別布設了1個控制點,像對中間布設1個控制點),旁向區(qū)域的4個控制點共用,同時增加了旁向連接,以保證像對之間的接邊精度。最后對整個工程進行了平差解算,完成衛(wèi)星影像區(qū)域網(wǎng)空三加密工作[12],像控點坐標殘差見表2。
表2 像控點坐標殘差 m
續(xù)表2
由表2可知,像控點殘差平面和高程誤差都優(yōu)于2m,滿足TB10050—2010《鐵路工程攝影測量規(guī)范》關(guān)于1∶10000空三加密精度要求[13]。
通過空三解算結(jié)果校正RPC參數(shù)。在全數(shù)字攝影測量工作站中,通過加載處理后的衛(wèi)星影像、經(jīng)空三平差校正完成的RPC參數(shù)建立核線影像,恢復立體模型,并設置中央子午線及投影面大地高,最后將模型坐標系歸化到控制點測圖坐標系中。定向完成后,將外業(yè)檢查點導入矢量測圖中進行精度統(tǒng)計[14]。一共檢查了20個外業(yè)檢查點,其統(tǒng)計報告見表3。
表3 模型量測外業(yè)檢查點誤差統(tǒng)計 m
由表3可知,利用資源三號衛(wèi)星影像進行測圖,其平面誤差優(yōu)于2.5m,高程誤差優(yōu)于2m,平面和高程誤差基本滿足《鐵路工程攝影測量規(guī)范》中高山地1∶10000地形圖精度要求[15]。
資源三號衛(wèi)星存檔數(shù)據(jù)較多,且影像姿態(tài)較穩(wěn)定,但分辨率一般。利用自由網(wǎng)平差及SRTM130m格網(wǎng)間距的DEM數(shù)據(jù),可快速完成數(shù)字正射影像圖的制作,其帶有坐標的數(shù)字影像,可與谷歌影像聯(lián)動,使導航定位更加準確,從而提高了刺點和調(diào)繪的作業(yè)效率。大面積衛(wèi)星數(shù)字正射影像有利于地勘專業(yè)進行室內(nèi)地質(zhì)判釋,同時為線路專業(yè)大范圍選線設計提供了有力的數(shù)據(jù)保障。
通過在每景影像布設5個控制點(四角各布設1個,中間再布設1個點)的技術(shù)方案,完成了1∶10000地形圖制作,并利用外業(yè)檢查點進行精度檢查,研究表明,其結(jié)果基本滿足《鐵路工程攝影測量規(guī)范》中高山地1∶10000地形圖精度要求。