錢志奇，朱曉雯，房瑩瑩，李文林，陸殿梅

(江蘇省地質(zhì)勘查技術(shù)院，江蘇 南京 210049)

0 引言

資源一號02C衛(wèi)星(ZY－1 02C星)發(fā)射于2011年12月22日，廣泛應(yīng)用于土地資源、礦產(chǎn)資源、地質(zhì)環(huán)境調(diào)查以及國土資源、地質(zhì)災(zāi)害應(yīng)急監(jiān)測等(李德仁等，2006;王海俠等，2011)。同部分高分辨率商業(yè)遙感衛(wèi)星(如IKONOS、QuickBird等)一樣，其傳感器信息暫時未向普通用戶公開(李德仁等，2006;楊國榮，2007;鄭琳等，2007;王海俠等，2011)。由于缺乏衛(wèi)星軌道星歷參數(shù)和傳感器參數(shù)，嚴格的共線方程模型便不適用，通常應(yīng)用有理函數(shù)模型(Rational Function Model，RFM)來解決從二維圖像反演三維信息的問題。近年來，RFM模型應(yīng)用于遙感技術(shù)日漸成熟，取得了較高的反演精度(張永生等，2004;胡堃，2009;李慶鵬等，2011;邵俊等，2010;袁修孝等，2012)。ZY－1 02C星HR影像數(shù)據(jù)提供了與QuickBird數(shù)據(jù)格式相似的RPB文件與經(jīng)過輻射校正后的柵格數(shù)據(jù)，可以使用RFM模型進行系統(tǒng)幾何校正或正射校正(胡鳳偉等，2012;胡鋒，2013;楊博等，2013;周偉等，2014)。常規(guī)的影像處理軟件(如ENVI、Erdas等)具有RPC正射校正功能，但數(shù)據(jù)文檔選擇繁瑣，影像讀寫耗時，且未提供自動化的批處理功能，對于大范圍海量數(shù)據(jù)的處理，不但耗時、耗力而且容易出錯。本次研究旨在利用IDL語言，為需要處理ZY－1 02C星高分影像的用戶提供一種快速有效的數(shù)據(jù)預(yù)處理途徑，推進國產(chǎn)資源衛(wèi)星數(shù)據(jù)的廣泛應(yīng)用。

1 RFM模型幾何校正

1.1 ZY－1 02C星HR數(shù)據(jù)

ZY－1 02C衛(wèi)星是我國高分辨率遙感數(shù)據(jù)衛(wèi)星之一，裝有全色多光譜相機和全色高分辨率相機HR，主要任務(wù)是獲取全色和多光譜圖像數(shù)據(jù)，數(shù)據(jù)可以用于1∶2.5萬和1∶5萬比例尺的礦產(chǎn)資源、土地資源、地質(zhì)環(huán)境調(diào)查業(yè)務(wù)。衛(wèi)星上搭載有2臺空間分辨率2.36 m的HR相機，拼接幅寬54 km。搭載有空間分辨率5 m的全色及10 m的多光譜相機，幅寬為60 km(胡鳳偉等，2012)。中國資源衛(wèi)星應(yīng)用中心提供的數(shù)據(jù)，包括0級的原始數(shù)據(jù)以及標準產(chǎn)品數(shù)據(jù)。標準產(chǎn)品包括1級的輻射校正產(chǎn)品和2級的系統(tǒng)幾何校正產(chǎn)品。0級數(shù)據(jù)經(jīng)輻射校正、數(shù)據(jù)質(zhì)量增強與簡單的數(shù)據(jù)篩選后生成1級輻射校正產(chǎn)品，期間將排除云量覆蓋較大，數(shù)據(jù)質(zhì)量不理想等無效數(shù)據(jù)。1級輻射校正產(chǎn)品無地圖投影信息，但提供RPC模型參數(shù)。2級系統(tǒng)幾何校正產(chǎn)品為1級輻射校正產(chǎn)品經(jīng)系統(tǒng)幾何校正處理后生成的數(shù)據(jù)產(chǎn)品，具有地圖投影信息。本次研究主要研究1級輻射校正產(chǎn)品，其文件夾內(nèi)關(guān)鍵文件見表1。

表1 主要使用文件Table 1 Main used files

由表1可知，HR數(shù)據(jù)包括2個HR傳感器的HR1與HR2影像，可經(jīng)幾何校正處理后鑲嵌而得。HR1與HR2影像的快速幾何校正與鑲嵌便是此次研究的重點內(nèi)容。

1.2 RFM模型

有理多項式函數(shù)模型(Rational Polynomial Coefficients，RPC)本質(zhì)上就是有理函數(shù)模型(Rational Function Model，RFM)，是傳感器幾何模型的一種抽象表達方式，適用于多種航空航天傳感器，是近似糾正模型更精確的形式(張永生等，2004;張過等，2007，2010;胡堃，2009;祝小勇等，2009;李慶鵬等，2011;宋妍等，2013)。采用有理多項式模型擬合傳感器內(nèi)外方位元素，可以將地面點空間D(X，Y，Z)作為自變量與對應(yīng)的像點坐標d(Line，Sample)采用比值多項式的方式關(guān)聯(lián)起來。模型的函數(shù)多項式系數(shù)即為RPC參數(shù)，該參數(shù)隨傳感器瞬時姿態(tài)變化。在沒有RPC參數(shù)的情況下，影像一般不會用于測繪(李德仁等，2006;蔡耀君等，2013)等高精度應(yīng)用。影像比值多項式如式(1)。

通常為增強模型參數(shù)的求解穩(wěn)定性，需將地面坐標和影像坐標正則化至［－1，l］范圍。其中(P，L，H)為正則化的地面坐標，(X，Y)為正則化的影像坐標，下標S代表列，下標L代表行，正則化表達式如式(2)所示：

式(2)中，Offset與Scale為地面坐標與影像坐標的正則化參數(shù)。式(1)每個多項式有20個多項式系數(shù)，與式(3)所示類似：

式(3)中，ai(i=0，1，2，…，19)即為各多項式系數(shù)，為保證比值多項式式(1)恒有意義，分母與常數(shù)項通常為1。式(1)中有4個多項式共包含80個多項式參數(shù)。利用該80個多項式參數(shù)與10個正則化參數(shù)便能進行物方空間坐標與像方空間坐標的轉(zhuǎn)化，實現(xiàn)地面坐標與影像坐標之間的轉(zhuǎn)換。

1.3 RFM模型應(yīng)用

本次研究主要以ZY－1 02C星的HR數(shù)據(jù)的1級輻射校正產(chǎn)品的幾何校正為主。采用RFM模型首先需要RPC參數(shù)，RPC標準化參數(shù)存放于同一文件夾內(nèi)的影像同名文件“*.rpb”文件中，該文件類型與QuickBird衛(wèi)星影像附帶文件類型相同，使用ENVI等軟件處理時，可以采用與QuickBird處理相同的方法進行數(shù)據(jù)顯示與正射校正等。首先以規(guī)則數(shù)據(jù)結(jié)構(gòu)讀取rpb文件中相關(guān)RPC參數(shù)信息(表2)，通過RFM模型反演影像實際地面坐標(P，L，H)，篩選地面點作為影像地面控制點進行影像幾何糾正，后以HR1為基準影像，對HR2進行基于HR1像元控制點的幾何配準，再將HR1與HR2進行鑲嵌等預(yù)處理，最后通過頭文件等信息初步判斷合成影像的合理性。詳細流程如圖1所示。

表2 RPC參數(shù)信息表Table 2 List of the RPC parameters

圖1 基于RPC的ZY－1 02C星HR數(shù)據(jù)預(yù)處理流程圖Fig.1 Flowchart showing pre-processing of the ZY-1 02C data based on RPC

2 RFM幾何校正模型實現(xiàn)

2.1 RFM模型的IDL程序設(shè)計

IDL(Interactive Data Language)語言是面向?qū)ο?、面向矩陣操作的交互式?shù)據(jù)語言(韓培友，2006)。利用其可視化控件IDLDrawWidget與非可視化控件COM_IDL_Connect可以實現(xiàn)與大部分主流開發(fā)平臺的交互與集成。該語言針對遙感影像處理封裝了大量函數(shù)供用戶調(diào)用，減少用戶從底層進行開發(fā)的麻煩(錢志奇等，2011;尹新沆等，2014)。考慮ZY－1 02C星HR數(shù)據(jù)可使用ENVI完成幾何校正等預(yù)處理，而ENVI正是以IDL語言開發(fā)的遙感應(yīng)用軟件，所以為實現(xiàn)影像大數(shù)據(jù)量批處理功能，采用IDL語言進行功能函數(shù)開發(fā)與整合。

根據(jù)流程編寫了4個自定義函數(shù)，分別是：READ_RPBFILE，用于讀取rpb文件內(nèi)RPC參數(shù);ZY_02C_AUTOREGIST_HR，用于通過RPC參數(shù)實現(xiàn)控制點反演與影像幾何校正;ZY_02C_MOSAIC，用于影像鑲嵌;READ_XML，用于讀取XML頭文件信息。由于篇幅有限，在此僅對調(diào)用的部分關(guān)鍵函數(shù)進行說明。

首先RPB文件結(jié)構(gòu)固定，數(shù)據(jù)量小，采用簡單文本逐行讀取即可，返回RPC數(shù)據(jù)結(jié)構(gòu)數(shù)組供使用，主要使用函數(shù)如表3。

表3 READ_RPBFILE函數(shù)主要代碼表Table 3 List of main codes for READ_RPBFILE function

在獲取RPC參數(shù)后，通過該參數(shù)反演地面控制點并進行投影變換與幾何校正，根據(jù)數(shù)據(jù)特點設(shè)計了ZY_02C_AUTOREGIST_HR功能函數(shù)。函數(shù)主要分為2個部分：一部分為RPC地面控制點反演，主要使用如表4所示函數(shù);另一部分為投影變換與幾何校正，主要使用了ENVI功能函數(shù)ENVI_DOIT下的ENVI_REGISTER_DOIT函數(shù)，該函數(shù)使用形式如下，具體參數(shù)說明可以查詢ENVI幫助文檔(下同)：

影像校正后，需要將HR1與HR2進行鑲嵌處理。為避免影像誤差過大導(dǎo)致鑲嵌效果不佳，故以HR2投影為基準(因為 HR2為左影像)，首先將HR1進行基于像元的控制點幾何校正，主要功能函數(shù)為調(diào)用 ENVI_DOIT下的 ENVI_AUTO_TIE_POINT函數(shù)，再采用ENVI_REGISTER_DOIT進行校正，函數(shù)使用形式如下：

表4 ZY_02C_AUTOREGIST_HR函數(shù)主要代碼表Table 4 List of main codes for ZY_02C_AUTOREGIST_HR function

將校正后的HR2與HR1進行鑲嵌處理，主要功能函數(shù)為調(diào)用ENVI_DOIT下的MOSAIC_DOIT函數(shù)，函數(shù)使用形式如下：

查詢XML頭文件信息，通常采用IDL下新建'IDLffXMLDOMDocument'對象的方法，通過XML數(shù)據(jù)結(jié)構(gòu)下對應(yīng)名稱查詢信息，使用形式如下：

此處獲取的頭文件信息是指HR1與HR2各自對應(yīng)頭文件，研究中主要用于獲取其相應(yīng)4個角點的經(jīng)緯度信息，用于驗證以上數(shù)據(jù)處理結(jié)果的正確性。

2.2 HR影像處理實例

IDL可以通過DIALOG_PICKFILE功能函數(shù)快速獲取指定文件夾內(nèi)特定數(shù)據(jù)類型數(shù)據(jù)，實驗僅以其中部分影像結(jié)果來判斷程序可靠性。圖2所示為HR2與HR1幾何校正鑲嵌過程示意圖。其中鑲嵌邊緣局部放大圖上可見，鑲嵌邊緣無明顯偏差痕跡，顏色過渡平滑，符合常規(guī)影像使用要求;由圖可知，經(jīng)過投影變換后，HR2與HR1鑲嵌結(jié)果并非為四邊形，在重疊區(qū)南北兩端存有偏差，若采用簡單像元線性匹配將會造成較大誤差。在不考慮羽化邊緣時，1景HR數(shù)據(jù)程序若采用ENVI軟件人工操作處理時需要幾十分鐘，而采用此程序運行僅需數(shù)十秒即完成鑲嵌結(jié)果，結(jié)果精度與人工操作無異。

3 結(jié)論

本次研究主要針對ZY－1 02C衛(wèi)星的1級輻射校正數(shù)據(jù)進行，所得結(jié)果相當(dāng)于其2級系統(tǒng)幾何校正數(shù)據(jù)，并未進行精幾何校正，數(shù)據(jù)精度取決于RPC參數(shù)精度。此程序能夠很好地實現(xiàn)HR數(shù)據(jù)的批處理功能，對于只有1級輻射校正數(shù)據(jù)的用戶可以參考進行數(shù)據(jù)幾何校正。同時，由于該衛(wèi)星多光譜數(shù)據(jù)具有相同數(shù)據(jù)形式，可修改后用作多光譜數(shù)據(jù)處理，并集成數(shù)據(jù)融合等預(yù)處理功能。本次研究亦有以下部分需要改進。

圖2 HR2與HR1幾何校正鑲嵌過程示意圖Fig.2 Schematic mosaic process of HR1 and HR2 geometric correction

(1)并非所有的1級輻射校正數(shù)據(jù)中均有DEM數(shù)據(jù)，所以未將DEM數(shù)據(jù)納入考慮范圍，故正射校正效果無法判斷，僅實現(xiàn)了系統(tǒng)幾何校正功能。

(2)采用反演控制點的方法能夠同時將控制點存儲，但缺乏對控制點的篩選，對于某些特殊地區(qū)，也許控制點會產(chǎn)生異常問題，還需改進。

(3)影像鑲嵌時需考慮羽化鑲嵌邊緣問題，否則對于某些區(qū)域在鑲嵌邊緣將產(chǎn)生色差。

(4)研究發(fā)現(xiàn)某些1級輻射校正數(shù)據(jù)邊緣有1條數(shù)個像元寬度的固定數(shù)值黑邊，數(shù)值會使鑲嵌結(jié)果形成1條明顯的鑲嵌黑邊。黑邊問題產(chǎn)生原因尚待研究，面對該問題實驗依賴人工裁剪，后續(xù)工作需要設(shè)計代碼忽略該黑邊值，使程序完善。

蔡耀君，華璀，盧遠，等.2013.基于缺失RPC文件的衛(wèi)星影像正射校正研究［J］.測繪與空間地理信息，36(12)：68－70.

韓培友.2006.IDL可視化分析與應(yīng)用［M］.陜西西安：西北工業(yè)大學(xué)出版社.

胡堃.2009.基于RFM通用成像模型的區(qū)域網(wǎng)平差研究［J］.現(xiàn)代商貿(mào)工業(yè)，21(5)：305－306.

胡鳳偉，胡龍華，李濤.2012.資源一號02C衛(wèi)星PMS數(shù)據(jù)處理方法研究［J］.華北科技學(xué)院學(xué)報，9(4)：40－44.

胡鋒.2013.衛(wèi)片RPC參數(shù)解算淺談［J］.江西測繪，(1)：14－39.

李德仁，張過，江萬壽，等.2006.缺少控制點的SPOT－5 HRS影像RPC模型區(qū)域網(wǎng)平差［J］.武漢大學(xué)學(xué)報：信息科學(xué)版，31(5)：377－381.

李慶鵬，王志剛，陳琦.2011.基于嚴格仿射變換模型的遙感影像RPC參數(shù)求解［J］.測繪信息與工程，36(3)：1－4.

李慶鵬，王志剛，雷玉飛.2011.CBERS－02B衛(wèi)星HR相機遙感影像區(qū)域網(wǎng)平差技術(shù)［J］.航天返回與遙感，32(4)：30－36.

錢志奇，張友靜，鄧世贊，等.2011.基于組件技術(shù)的地表蒸散發(fā)遙感監(jiān)測系統(tǒng)設(shè)計與應(yīng)用［J］.地球信息科學(xué)學(xué)報，13(5)：679－686.

邵俊，郭建寧.2010.CBERS－02B衛(wèi)星HR相機遙感影像RPC校正方法［J］.航天返回與遙感，31(2)：29－37.

宋妍，范高晶，左佳.2013.多源高空間分辨率遙感圖像平差方法［J］.國土資源遙感，25(2)：21－26.

王海俠，高飛，胡小華.2011.基于RPC模型的QuickBird影像幾何精校正的研究與探討［J］.城市勘測，(3)：54－58.

楊國榮.2007.IKONOS立體像對提取DEM及正射糾正的實驗研究［J］.現(xiàn)代測繪，30(5)：26－29.

袁修孝，汪韜陽.2012.CBERS－02B衛(wèi)星遙感影像的區(qū)域網(wǎng)平差［J］.遙感學(xué)報，16(2)：310－324.

楊博，王密.2013.資源一號02C衛(wèi)星全色相機在軌幾何定標方法［J］.遙感學(xué)報，17(5)：1175－1190.

尹新沆，張友靜，錢志奇，等.2014.HJ－1B星地表蒸散發(fā)遙感估算系統(tǒng)的設(shè)計與應(yīng)用［J］.地理空間信息，12(1)：89－92.

張永生，劉軍.2004.高分辨率遙感衛(wèi)星立體影像RPC模型定位的算法及其優(yōu)化［J］.測繪工程，13(1)：1－4.

張過，李德仁.2007.衛(wèi)星遙感影像RPC參數(shù)求解算法研究［J］.中國圖象圖形學(xué)報，12(12)：2080 －2088.

鄭琳，陳鷹，林怡.2007.SPOT影像的RPC模型糾正［J］.測繪與空間地理信息，30(2)：16－19.

祝小勇，張過，秦緒文.2009.國產(chǎn)光學(xué)衛(wèi)星影像RPC制作［J］.國土資源遙感，21(2)：32 －40.

張過，費文波，李貞，等.2010.用RPC替代星載SAR嚴密成像幾何模型的試驗與分析［J］.測繪學(xué)報，39(3)：264－270.

周偉，黃煒，王彥佐，等.2014.資源一號02C衛(wèi)星數(shù)據(jù)管理與服務(wù)系統(tǒng)研建［J］.國土資源遙感，26(1)：179－185.