王忠武 劉順喜 戴建旺 尤淑撐 孟超
(中國土地勘測規劃院,北京 100035)
“高分二號”(GF-2)衛星于2014年8月19日成功發射,標志著我國遙感衛星進入亞米級高分辨率時代。GF-2衛星搭載了2臺全色多光譜傳感器,每臺傳感器可提供幅寬不低于23km、空間分辨率優于1m的全色影像和優于4m的多光譜影像,經數據處理獲得的高分辨率多光譜影像可作為土地利用現狀調查、土地利用動態監測等高精度土地基礎數據采集工作的基礎數據源之一。由于多光譜與全色影像高精度配準是保證影像數據應用效率和效果的前提[1],因此,多光譜與全色影像配準策略對大規模、業務化影像預處理效率影響顯著[2-5]。
影像配準方法主要有自動配準方法和糾正后配準方法,自動配準方法采用區域灰度相關或空間特征匹配同名點自動提取,并采用薄板樣條插值技術,實現影像間自動配準;糾正后配準方法主要對經過幾何糾正后的影像進行配準[6-8]。在規模化應用中,需根據影像配準精度和處理時間確定配準方法,本文研究了 GF-2衛星全色與多光譜影像配準方法,為形成業務化應用中的全色和多光譜影像配準策略提供技術依據。
依據土地資源遙感調查監測技術規范,采用多景覆蓋平原區、山區的 GF-2衛星完整多光譜與全色影像,開展原始多光譜與全色影像自動配準和糾正后配準方法的配準精度對比試驗與分析,總結形成規模化數據應用中較好的配準策略,具體技術流程如圖1所示。

圖1 試驗技術流程Fig.1 Workflow of experiments
根據文獻[9-11]的要求,多光譜與全色影像配準要求配準結果目視無重影,配準控制點殘差需滿足平地、丘陵地不超過1倍基準影像分辨率,山地、高山地不超過1.5倍基準影像分辨率。
影像配準方法包括自動配準方法和糾正后配準方法。
(1)自動配準方法
影像自動配準通常包括匹配特征提取、匹配關系計算這2個步驟[6]。匹配特征提取主要包括基于區域的方法、基于特征的方法、結合區域和特征的方法等三類。基于區域的方法利用一定區域中對應像素灰度的統計特征來識別同名特征點,例如差平方和、相關函數、歸一化互相關系數、互信息、梯度、相位相關等信息。基于特征的方法利用影像間的點、線、面等同名特征描述參數來配準特征。結合區域和特征的方法將以上2種方法的優點結合起來,綜合利用灰度統計信息和空間特征進行匹配特征提取。匹配關系計算主要利用提取的匹配特征,基于薄板樣條等假設計算待配準影像間的函數關系,實現影像配準。通常來說,自動配準方法對于相同或相近幾何變形的待配準影像具有較好的配準效果。
ERDAS IMAGINE軟件提供了AutoSysnc模塊[12],該模塊主要利用基于區域的方法,通過點匹配算法自動產生大量同名特征點,進而基于同名特征點構造待配準影像間的數學關系函數模型,實現影像自動配準。該模塊可基本實現全自動影像配準,具有較高的配準精度,是目前普遍采用的自動配準工具之一。
(2)糾正后配準方法
對于具有不同幾何變形,或影像自動配準方法難以達到較高配準精度的待配準影像,則需要首先對待配準影像進行幾何糾正,基本消除各待配準影像內部的幾何變形,或將待配準影像的不同幾何變形糾正為相近幾何變形,然后進行配準。采用該方法的影像制圖步驟依次為:全色影像糾正、多光譜影像糾正、全色多光譜影像配準、全色多光譜影像融合、融合影像糾正、鑲嵌勻光裁切,其中融合影像糾正為根據需要選擇是否再次進行影像糾正。
GF-2衛星1A級全色和多光譜影像分別自帶RPC(rational polynomial coefficient)參數,適合于采用 RPC模型進行幾何糾正。RPC模型的實質是利用有理多項式模擬影像成像瞬間的幾何關系[13-15],將地面點大地坐標X、Y、Z與其對應像點坐標x、y進行關聯。RPC模型可表示為

采用2014年10月25日、11月10日、12月4日接收的10景GF-2衛星1A級全色與多光譜影像作為試驗數據,見表1。影像側視角在4.950 22°~11.000 90°之間,進行影像自動配準和糾正后配準試驗。試驗數據中,第一個傳感器影像2景,第二個傳感器影像8景;4景位于內蒙古自治區赤峰市元寶山區,試驗區地形主要為山區;4景位于江蘇省無錫市市轄區市,試驗區地形主要為平原區;2景位于浙江省寧波市奉化市,試驗區地形主要為山區。

表1 影像配準試驗數據Tab.1 Test data for image registration
采用自動配準方法,以原始全色影像為參考,采用ERDAS軟件的AutoSysnc模塊對多光譜影像進行自動配準,在原始全色影像和配準后多光譜影像上均勻選擇10個同名點作為檢查點,評價影像自動配準精度,結果如表2所示,可以看出,山區6景影像的配準中誤差分別為0.94、0.93、0.88、0.83、0.90、0.88個像元,最大誤差分別為1.50、1.49、1.97、1.27、1.30、1.53個像元;平原區4景影像的配準中誤差分別為0.99、0.88、0.96、0.84個像元,最大誤差分別為1.41、2.45、1.44、2.93個像元。總體來說,配準中誤差為山區0.83~0.94個像元,平原區0.84~0.99個像元,滿足影像配準技術指標要求。

表2 原始影像自動配準精度Tab.2 Automatic registration precision of original images像元
采用無控制點的RPC模型分別對全色和多光譜影像進行幾何糾正,在糾正后的全色影像和多光譜影像上均勻選取10個同名地物點作為檢查點,以全色影像為參考值,量測多光譜影像的偏移量,評價糾正后影像配準精度,結果如表3所示,可以看出,山區 6景影像的配準中誤差分別為2.69、1.62、1.92、2.45、1.86、2.97個像元,最大誤差分別為 4.90、2.83、3.70、4.29、2.86、4.02個像元;平原區 4景影像的配準中誤差分別為1.56、2.26、1.97、3.24個像元,最大誤差分別為2.05、4.66、3.06、5.46個像元。總體來說,配準中誤差為山區1.62~2.97個像元、平原區1.56~3.24個像元,不滿足影像配準技術指標要求。

表3 糾正后影像配準精度Tab.3 Registration precision of corrected images像元
通過對10景覆蓋山區和平原區的GF-2衛星完整景1A級多光譜與全色影像進行的配準試驗,可以發現,相對于無控制點RPC模型糾正后配準,原始影像自動配準的精度能滿足相關技術規范要求,同時采用原始影像自動配準方法的影像預處理步驟為“全色多光譜影像配準、全色多光譜影像融合、融合影像糾正”,相對于“全色影像糾正、多光譜影像糾正、全色多光譜配準、全色多光譜影像融合”的處理步驟,能有效縮短影像處理時間,是規模化數據應用中影像配準的較好方法。需要指出的是,雖然本文采用多景真實影像進行了試驗,但總體上試驗數據數量和類型有限,測試結果的客觀性會受到一定影響,相關結論還需要在后續大量數據應用中進一步驗證與完善。
References)
[1]柴淵, 李萬東. 土地利用動態監測技術與方法[M]. 北京: 地質出版社, 2011.CHAI Yuan, LI Wandong. Land Use Dynamic Monitoring Technology and Method[M]. Beijing: Geological Press, 2011. (in Chinese)
[2]王洪海, 徐建艷, 龍小祥. 一種衛星多光譜圖像亞像元波段配準精度自動評價方法[J]. 航天返回與遙感, 2010, 31(3):55-63.WANG Honghai, XU Jianyan, LONG Xiaoxiang. An Automated Evaluation Method of Sub-pixel Band-to-band Registration Accuracy for Satellite Multispectral Image[J]. Spacecraft Recovery amp; Remote Sensing, 2010, 31(3): 55-63. (in Chinese)
[3]朱近, 司美玲, 夏德深, 等. 用局部灰度極值方法進行多光譜圖像配準算法研究[J]. 航天返回與遙感, 2010, 31(5): 32-37.ZHU Jin, SI Meiling, XIA Deshen, etal. Research on Multi-spectral Image Registration with Local Gray Extremism Method[J]. Spacecraft Recovery amp; Remote Sensing, 2010, 31(5): 32-37. (in Chinese)
[4]張國瑞. 多光譜多鏡頭CCD相機的光學配準[J]. 航天返回與遙感, 2001, 22(1): 51-54.ZHANG Guorui. Optical Registration for Multi-spectral and Multi-lens CCD Camera[J]. Spacecraft Recovery amp; Remote Sensing, 2001, 22(1): 51-54. (in Chinese)
[5]郭悅, 楊樺. 多光譜CCD相機的光學配準[J]. 航天返回與遙感, 2003, 24(1): 29-32.GUO Yue, YANG Hua. Optical Registration of the Multi-spectral CCD Camera[J]. Spacecraft Recovery amp; Remote Sensing,2003, 24(1): 29-32.(in Chinese)
[6]王瑞瑞, 馬建文, 陳雪. 多源遙感影像自動配準技術的研究進展[J]. 遙感信息, 2011(3): 121-127.WANG Ruirui, MA Jianwen, CHEN Xue. The Research State of the Automatic Registration Techniques of Multi-sensor Remote Sensing Imagery[J]. Remote Sensing Information, 2011(3): 121-127.(in Chinese)
[7]王瑞瑞, 馬建文. 多尺度、多角度異源遙感影像配準方法研究[J]. 計算機工程與應用, 2012, 48(30): 5-9.WANG Ruirui, MA Jianwen. Registration Method for Heterogeneous Remote Sensing Images with Multi-scale and Multi-angle[J]. Computer Engineering and Applications, 2012, 48(30): 5-9.(in Chinese)
[8]田上成, 張喬, 劉保成, 等. 基于傳感器幾何特性和圖像特征的影像配準方法[J]. 航天返回與遙感, 2013, 34(5): 85-89.TIAN Shangcheng, ZHANG Qiao, LIU Baocheng, etal. An Improved Matching Method for Image Based on Sensor Geometric Features and Image Features[J]. Spacecraft Recovery amp; Remote Sensing, 2013, 34(5): 85-89. (in Chinese)
[9]中華人民共和國國土資源部. TD/T 1014-2007. 第二次全國土地利用調查技術規程[S]. 北京:中華人民共和國國土資源部, 2007.MLR of PRC. TD/T 1014-2007. Technical Regulation for the Second National Land-use Investigation[S]. Beijing: MLR of China. 2007 (in Chinese)
[10]贠小蘇. 第二次全國土地調查培訓教材[M]. 北京: 中國農業出版社, 2007.YUN Xiaosu. The Second National Land Survey Training Materials[M]. Beijing: Chinese Agriculture Press, 2007. (in Chinese)
[11]高延利, 朱留華, 王廣華, 等. TD-1010土地利用動態遙感監測規程[S]. 北京: 中華人民共和國國土資源部, 1999.GAO Yanli, ZHU Liuhua, WANG Guanghua, etal. TD-1010 Regulation of Land-use Dynamic Remote Sensing Monitoring[S]. Beijing: MLR of China, 1999. (in Chinese)
[12]萬里紅, 施艷麗, 楊武年. 基于AutoSync的快速影像配準研究[J]. 測繪, 2010, 33(3): 103-105.WAN Lihong, SHI Yanli, YANG Wunian. Research of Quick Image Registration Based on AutoSync[J]. Surveying and Mapping, 2010, 33(3): 103-105. (in Chinese)
[13]張過. 缺少控制點的高分辨率衛星遙感影像幾何糾正[D]. 武漢: 武漢大學, 2005.ZHANG Guo. Rectification for High Resolution Remote Sensing Image under Lack of Ground Control Points[D]. Wuhan:Wuhan University, 2005. (in Chinese)
[14]王忠武, 郝偉, 劉順喜, 等. 萊索托地區北京一號全色影像幾何校正方法研究[J]. 中國土地科學, 2012, 26(9): 78-82.WANG Zhongwu, HAO Wei, LIU Shunxi. Research on Geometric Correction Method for the BJ-1 Panchromatic Image Covering Kingdom of Lesotho[J]. China Land Sciences, 2012, 26(9): 78-82. (in Chinese)
[15]張過, 李揚, 祝小勇, 等. 有理函數模型在光學衛星影像幾何糾正中的應用[J]. 航天返回與遙感, 2010, 31(4): 51-57.ZHANG Guo, LI Yang, ZHU Xiaoyong, etal. Application of RFM in Geometric Rectification of Optical Satellite Image[J].Spacecraft Recovery amp; Remote Sensing, 2010, 31(4): 51-57. (in Chinese)