線陣垂軌環(huán)掃式光學遙感衛(wèi)星影像幾何糾正

薛 武，王 鵬，鐘靈毓

（航天工程大學，北京 101416）

1 引言

地面空間分辨率和成像幅寬是現(xiàn)代光學遙感衛(wèi)星設計中十分重要的指標，對于線陣推掃式衛(wèi)星而言，二者通常互相矛盾，在實際應用中需要進行折衷考慮。為了在不降低空間分辨率的前提下提高影像的幅寬，主流的線陣推掃式光學遙感衛(wèi)星通常采用多片CCD 拼接、多相機多角度視場拼接、衛(wèi)星姿態(tài)快速機動下的寬覆蓋影像獲取和多星組網(wǎng)等解決方案［1-8］。以上方式雖然在一定程度上提高了衛(wèi)星的成像幅寬，但存在成本高、連續(xù)大范圍成像能力不足、幅寬提升有限等問題，沒有從根本上解決影像的地面空間分辨率與成像幅寬之間的矛盾。

隨著衛(wèi)星平臺技術(shù)、姿態(tài)控制技術(shù)以及光學傳感器技術(shù)的發(fā)展與進步，線陣環(huán)掃成像的設計思想被提出和實踐。通過將多條線陣CCD 沿衛(wèi)星飛行方向安置，相機鏡頭繞飛行方向360°旋轉(zhuǎn)，旋轉(zhuǎn)過程中完成對地面的掃描成像，在保證地面分辨率的前提下可以大幅拓展影像的幅寬，這就是線陣環(huán)掃式衛(wèi)星的成像基本原理［8-9］。以軌道高度500 km 的光學遙感衛(wèi)星為例，當掃描區(qū)域張角達到130°時（即左右側(cè)擺角度65°時），其成像幅寬可以接近3 000 km，相對線陣推掃式光學遙感衛(wèi)星提升了近10 倍，大大增強了廣域態(tài)勢感知能力，具有重要的應用價值。

由于采用環(huán)掃成像的方式，同一景影像內(nèi)不同像元成像的角度差異較大，導致其對應的地面分辨差異較大，最大可達5 倍。特別是靠近影像邊緣處，由于大角度側(cè)視成像、地形起伏等導致影像產(chǎn)生嚴重的幾何畸變，而且?guī)缀位冊谟跋裆戏植疾灰?guī)律，現(xiàn)有的糾正模型難以描述其分布特點，如不進行針對性處理則無法進行影像判讀解譯等工作。

現(xiàn)有的光學遙感衛(wèi)星影像幾何糾正方法主要針對線陣推掃式或面陣影像，此類影像的幅寬相對較小，同一景影像內(nèi)成像視角差異不大，幾何畸變較小，影像幾何糾正理論相對比較成熟。特別是有理函數(shù)模型的廣泛應用，其簡潔的表達形式為影像幾何糾正提供了便利條件［10］。在傳感器成像方式?jīng)]有發(fā)生顯著變化的近10 年中，光學衛(wèi)星影像幾何糾正的理論沒有大的變化和發(fā)展。隨著線陣擺掃成像理念的提出和工程實踐，影像幾何糾正理論與方法也要與時俱進，以滿足影像高精度處理與應用的需求。線陣環(huán)掃影像的成像模型與現(xiàn)有推掃成像衛(wèi)星區(qū)別較大，加上幅寬達到幾千公里，地形起伏跨度大，影像幾何變形大且變形分布不規(guī)律，直接采用現(xiàn)有的理論方法無法得到高精度的糾正結(jié)果，必須研究建立相應的幾何糾正模型和方法。

本文在深入分析線陣環(huán)掃衛(wèi)星成像原理和畸變特點的基礎上，構(gòu)建了其嚴格成像模型，依據(jù)嚴格成像模型，利用正射影像和公開DEM 數(shù)據(jù)進行了影像模擬生成。在深入分析影像成像特點的基礎上，提出了一套適用于線陣環(huán)掃衛(wèi)星影像的幾何糾正方法，并通過實驗驗證了其有效性。

2 線陣環(huán)掃式光學遙感衛(wèi)星成像原理、嚴格成像模型構(gòu)建及影像模擬

2.1 成像原理

多條線陣CCD 拼接為一條長線陣CCD，沿飛行方向安置，在衛(wèi)星飛行過程中，探測器光軸繞衛(wèi)星的前進方向360°旋轉(zhuǎn)環(huán)掃，在掃描方向上視場不斷疊加，并沿著飛行方向推進向前平移，以達到大視場寬幅成像的目的，其成像原理如圖1 所示。這種成像方式綜合利用了衛(wèi)星的機動能力和成像能力，思路新穎，在保證滿足一定空間分辨率要求的前提下擴大了衛(wèi)星成像幅寬。

圖1 線陣環(huán)掃式成像原理示意圖Fig.1 Schematic diagram of linear array circular scanning imaging principle

在掃描過程中，傳感器在特定角度范圍內(nèi)成像，該過程類似于單程擺掃的線陣擺掃成像，但不存在掃描反射鏡的相對運動，攝影中心S與像點p、物點P之間滿足攝影測量共線條件關系，每一條掃描行均滿足的線中心投影關系［8-9］。因此，成像過程可描述為光線的坐標系轉(zhuǎn)動，可以參考線陣推掃式傳感器的嚴格成像模型，構(gòu)建垂軌環(huán)掃式傳感器的坐標轉(zhuǎn)換關系，建立相應嚴格成像模型。

2.2 嚴格成像模型構(gòu)建

線陣垂軌環(huán)掃的特殊成像方式，導致現(xiàn)有的光學衛(wèi)星影像成像模型不再適用，必須建立專屬的成像模型為目標定位、影像幾何糾正提供依據(jù)。嚴格成型模型的構(gòu)建主要是依據(jù)在成像時投影中心、像點、物方點三點共線的幾何關系，構(gòu)建矢量共線的方程，即在同一坐標系下，投影中心至像點和物方點的矢量滿足共線關系［11-13］。具體如下：

在地心地固坐標系（WGS84）下，從GPS 相位中心到物方點的矢量為：

在像空間坐標系下，從投影中心到像方點的矢量為：

構(gòu)建嚴格成像模型，需要進行一系列的坐標轉(zhuǎn)換，將兩個矢量統(tǒng)一到相同坐標系下。首先，將物方矢量從地心地固坐標系轉(zhuǎn)換到天球坐標系：

再將物方矢量旋轉(zhuǎn)到衛(wèi)星本體坐標系：

最后，將其轉(zhuǎn)換到像空間坐標系：

根據(jù)物方矢量和像方矢量共線條件，可得如下等式：

式（7）即為線陣環(huán)掃式衛(wèi)星影像的嚴格成像模型，亦可將像方矢量通過一系列的坐標轉(zhuǎn)換到地心地固坐標系，二者是等價的，此處不再贅述。

需要說明的是：線陣垂軌環(huán)掃衛(wèi)星影像的嚴格成像模型雖然在形式上與線陣推掃式衛(wèi)星比較相似，但部分符號的含義有顯著的差異，特別是像空間坐標；另外，由于相機始終處于擺掃狀態(tài)，Rcbaomdy是隨時間一直在變化的變量，這些在應用中都要注意。

2.3 影像模擬生成

盡吾所知，根據(jù)公開報道，目前國內(nèi)外尚無類似成像方式的光學遙感衛(wèi)星影像可供開展幾何糾正參考使用。為獲取實驗數(shù)據(jù)，開展了影像模擬生成研究，開發(fā)了線陣環(huán)掃式衛(wèi)星影像模擬軟件，軟件主界面如圖2 所示。該模擬軟件能夠根據(jù)線陣環(huán)掃衛(wèi)星的成像模式和指標參數(shù)，考慮影響影像幾何質(zhì)量和輻射質(zhì)量的主要因素，基于衛(wèi)星成像過程全鏈路物理模型，生成反映衛(wèi)星在軌成像特性、真實性好的模擬影像，為衛(wèi)星入軌前開展數(shù)據(jù)處理研究提供真實性好的數(shù)據(jù)源。軟件的主要功能包括模擬參數(shù)設置（參數(shù)包括衛(wèi)星軌道根數(shù)、成像時長、大氣條件、相機幾何和輻射參數(shù)、姿態(tài)模型、姿態(tài)抖動等）、模擬影像生成與瀏覽展示、成像過程動畫展示等。由于影像模擬過程進行了光學成像的全鏈路仿真，計算復雜，加上影像數(shù)據(jù)量大，單景影像數(shù)據(jù)量超過300 GB，模擬軟件采用了CPU-GPU 協(xié)同并行加速的計算架構(gòu)，充分挖掘硬件的計算能力，大幅提高了影像模擬的效率。

圖2 影像模擬軟件主界面Fig.2 Main interface of image simulation software

基于影像模擬軟件生成了豐富的模擬影像數(shù)據(jù)，主要步驟包括衛(wèi)星軌道參數(shù)設置、相機成像參數(shù)設置、大氣條件設置、輸入DOM 和DEM、模擬影像生成、影像預處理等。模擬影像充分反映了衛(wèi)星成像特點，具有嚴密的幾何關系，成像區(qū)域選擇在印度北部地區(qū)，影像幅寬約2 000 km，地面場景涵蓋多種地貌地物類型，其基本參數(shù)如表1 所示。印度北部地區(qū)2 景模擬影像如圖3（a）和圖3（b）所示，影像細節(jié)如圖3（c）所示。需要說明的是，由于生成模擬影像所采用的DOM 在個別小區(qū)域有缺失，所以模擬影像上對應區(qū)域也存在缺失問題，但并不影響本文的研究。

圖3 模擬影像及細節(jié)展示Fig.3 Simulated images and detail display

表1 模擬影像主要參數(shù)Tab.1 Main parameters of simulated images

3 影像幾何糾正關鍵技術(shù)

3.1 基于嚴格成像模型的幾何粗糾正

在構(gòu)建的嚴格成像模型的基礎上，采用間接法對模擬影像進行數(shù)字微分糾正［14］，主要流程如圖4 所示，基本步驟如下：

圖4 幾何粗糾正流程圖Fig.4 Flow chart of geometric rough correction

（1）計算影像地面覆蓋范圍：利用相機成像參數(shù)、影像軌道姿態(tài)參數(shù)和DEM 數(shù)據(jù)計算影像4個角點坐標，得到影像的地面覆蓋范圍；

（2）根據(jù)輸出分辨率要求，逐像素計算輸出影像對應的地面坐標；

（3）利用DEM 內(nèi)插高程得到當前像元對應地面點的高程值；

（4）利用嚴格成像模型計算當前點在原始影像上的像素坐標；

（5）采用雙線性內(nèi)插得到當前點的灰度值；

（6）逐像素計算糾正后影像的灰度值，得到糾正后的影像。

為提高幾何粗糾正的效率，實際開發(fā)中采用一系列優(yōu)化策略，包括影像灰度賦值算法優(yōu)化、CPU-GPU 協(xié)同并行加速等。同時，考慮到多種場景下的需求，提供了按照指定地理坐標范圍糾正和輸出影像產(chǎn)品的功能，可以快速輸出感興趣區(qū)域的局部影像。

幾何粗糾正僅需要相機幾何參數(shù)、衛(wèi)星軌道姿態(tài)、DEM 等數(shù)據(jù)的支持，具有原理簡單、處理速度較快、無需人工干預等優(yōu)點，能夠滿足國土資源普查、災害評估、地理國情監(jiān)測等對精度要求不高的場景。

3.2 基于影像匹配的幾何精糾正

由于幾何粗糾正的精度有限，難以滿足國土測繪、目標地理定位、變化檢測等高精度應用的需求，需要在粗糾正的基礎上進行影像精糾正，提供高精度影像產(chǎn)品。影像精糾正的基本思路是在粗糾正的基礎上，通過與參考正射影像進行匹配，得到大量同名點坐標作為約束條件對影像進行精糾正，具體步驟如下：

3.2.1 金字塔逐級引導的影像快速匹配

線陣垂軌環(huán)掃影像幅寬達數(shù)千公里，影像數(shù)據(jù)量很大，印度地區(qū)1 景原始影像大小約312 GB，幾何粗糾正后約583 GB，加上影像側(cè)擺角度最大達60°，成像角度和正射影像差別較大，而且影像的地面分辨率不均勻，星下點與邊緣處差別最大可達5 倍，這些因素給影像匹配算法帶來了很大的挑戰(zhàn)。本文設計了一種金字塔逐級引導的影像快速匹配算法，有效解決了線陣垂軌環(huán)掃影像匹配中的幾何變形大、輻射差異大、分辨率差異大、視角差異大以及數(shù)據(jù)量大等系列問題［15］，取得了較好的效果，算法的具體流程如下：

（1）分別建立粗糾正結(jié)果影像和對應區(qū)域參考正射影像的金字塔；

（2）在金字塔頂層影像上均勻劃分格網(wǎng)，在格網(wǎng)內(nèi)提取Harris 特征點，每個格網(wǎng)內(nèi)保留特征最明顯的特征點，作為候選匹配點；

（3）基于影像已有的地理坐標信息，合理設置搜索空間，采用相關系數(shù)的匹配方法，進行影像匹配，得到初始種子點；

（4）將種子點坐標傳遞至影像金字塔中分辨率更高的一層，利用種子點構(gòu)建不規(guī)則三角網(wǎng)，然后在每個三角網(wǎng)內(nèi)部提取Harris 特征點，每個格網(wǎng)內(nèi)保留特征最明顯的特征點，作為候選匹配點；

（5）基于不規(guī)則三角網(wǎng)的約束關系，采用相關系數(shù)的匹配方法進行影像匹配，得到加密后的種子點；

（6）重復步驟（4）～（5），直至原始分辨率影像，得到粗糾正后影像和參考正射影像上大量的同名點坐標；

（7）采用虛警值最小化的誤匹配剔除技術(shù)剔除誤匹配點。

在特征提取與匹配的過程中，也采用了多線程加速等策略對算法進行了優(yōu)化，實現(xiàn)了超大影像的快速匹配。

3.2.2 基于樣條函數(shù)的幾何精糾正

精糾正的基本原理是利用影像匹配得到的大量同名點作為控制點對粗糾正后影像上殘留的幾何畸變進行擬合，求出變換參數(shù)，再利用變換參數(shù)改正影像上所有像點殘留的畸變值。

傳統(tǒng)光學衛(wèi)星影像通常采用一定階數(shù)的多項式來進行幾何糾正模型參數(shù)的擬合，能夠較好補償影像的系統(tǒng)性幾何畸變。但由于成像視角差異較大、影像幅寬很大、地形起伏等原因，粗糾正后的線陣環(huán)掃式衛(wèi)星影像中殘留的誤差不僅有系統(tǒng)性誤差，還有部分非系統(tǒng)性誤差，幾何變形不規(guī)律，采用多項式擬合存在精度有限和振蕩現(xiàn)象，總體效果不夠理想［16］。

由于樣條函數(shù)構(gòu)造簡單，使用方便，擬合準確，克服了高次多項式的振蕩現(xiàn)象，具有較好的數(shù)值穩(wěn)定性和收斂性，本文采用三次樣條函數(shù)。其實質(zhì)是由多段三次多項式拼合成的曲線，在拼接處，不僅函數(shù)本身是連續(xù)的，其一階和二階導數(shù)也是連續(xù)的。三次樣條函數(shù)的具體構(gòu)建方法及參數(shù)解算可參考文獻［17］，此處不再贅述。

利用樣條函數(shù)得到粗糾正后影像和參考影像之間的變換參數(shù)后，采用間接法幾何糾正的方法對粗糾正結(jié)果影像進行重采樣，得到精糾正結(jié)果影像。

3.3 幾何糾正精度評價

為了評價幾何糾正的精度，將糾正后的影像和基準影像進行特征匹配，得到同名點作為檢查點，統(tǒng)計檢查點坐標差值中誤差，即得到幾何糾正的精度。用(xi，yi)表示糾正后影像上檢查點坐標，(xi0，yi0)表示參考正射影像上同名點坐標，n表示檢查點的數(shù)量，幾何糾正中誤差σ的具體計算過程如下：

4 實驗與結(jié)果

4.1 幾何糾正實驗

利用印度北部地區(qū)的兩景模擬影像開展了線陣環(huán)掃式衛(wèi)星影像幾何糾正實驗。實驗的主要過程包括：首先，利用SRTM 開源DEM 數(shù)據(jù)和影像的軌道姿態(tài)數(shù)據(jù)進行影像粗糾正；然后利用3.2 節(jié)的方法進行影像精糾正；最后，利用3.3節(jié)的方法分別對粗糾正和精糾正的結(jié)果進行精度評價。粗糾正影像與正射影像的匹配結(jié)果如圖5 所示（綠色點表示同名點），圖中部分影像同名點如圖6 所示（黃色“×”表示同名點位置），匹配時間統(tǒng)計如表2 所示。影像幾何精糾正結(jié)果疊加在參考正射影像上的效果如圖7 所示。幾何糾正前后影像對比如圖8 所示。其中，圖8（a）是距離星下點約400 km 處某工廠影像；幾何糾正前，其形狀發(fā)生了較為明顯的變化，通過幾何糾正后恢復了其正確的形狀，如圖8（b）所示，疊加在參考正射影像上卷簾顯示如圖8（c）所示。圖8（d）是距離星下點約800 km 處某機場影像。幾何糾正前，其幾何變形十分嚴重，若不進行幾何糾正很難正確判明影像中目標的具體類型；經(jīng)過幾何糾正后恢復了其正確的形狀，如圖8（e）所示，疊加在參考正射影像上卷簾顯示如圖8（f）所示。印度地區(qū)2 景影像幾何粗糾正和精糾正的精度統(tǒng)計結(jié)果如表3～6 所示。

表2 影像匹配耗時統(tǒng)計Tab.2 Image matching time-consuming statistics

表3 模擬影像1 粗糾正結(jié)果精度統(tǒng)計Tab.3 Accuracy statistics of coarse correction results of simulation image 1

圖5 影像匹配結(jié)果Fig.5 Result of image matching

圖6 影像同名點Fig.6 Corresponding points on images

圖7 模擬影像精糾正結(jié)果Fig.7 Accurate correction results of simulated images

圖8 幾何精糾正前后對比Fig.8 Comparison before and after geometric accurate correction

表4 模擬影像1 精糾正結(jié)果精度統(tǒng)計Tab.4 Accuracy statistics of accurate correction results of simulation image 1

表5 模擬影像2 粗糾正結(jié)果精度統(tǒng)計Tab.5 Accuracy statistics of coarse correction results of simulation image 2

表6 模擬影像2 精正結(jié)果精度統(tǒng)計Tab.6 Accuracy statistics of accurate correction results of simulation image 2

4.2 實驗分析

從影像匹配結(jié)果來看，本文設計的金字塔逐級引導的影像快速匹配算法，有效解決了線陣垂軌環(huán)掃影像匹配中的幾何變形大、輻射差異大、分辨率差異大、視角差異大以及數(shù)據(jù)量大等系列問題，得到了數(shù)量充足、分布均勻、可靠性高的同名點，在荒漠、建筑物、植被等多種地表覆蓋場景中均取得了較好的效果。600 GB 左右的粗糾正結(jié)果影像與正射參考影像匹配1.5 萬個同名點耗時僅約30 min，計算效率顯著高于市面商業(yè)軟件，滿足實際生產(chǎn)作業(yè)需求。

從影像幾何糾正效果圖來看，幾何糾正能夠較好的解決影像幾何畸變嚴重的問題，糾正后的影像呈現(xiàn)中間窄、兩側(cè)寬的特點，影像形狀與衛(wèi)星成像設計以及可視化仿真結(jié)果高度一致。影像中由于成像角度變化、地形起伏產(chǎn)生的幾何畸變得到了較好的糾正，與參考正射影像卷簾檢查時，具有很好的吻合性。特別是在影像成像邊緣處，原始影像由于幾何畸變嚴重幾乎無法判明目標的類型，經(jīng)過幾何糾正后，較好地消除了幾何變形，能夠滿足判讀分析的應用要求。分析幾何糾正精度統(tǒng)計結(jié)果，幾何粗糾正精度優(yōu)于5 Pix?el，精糾正精度達到1 Pixel，能夠滿足實際應用的需求，也驗證了本文所提出的線陣環(huán)掃式光學遙感衛(wèi)星影像幾何糾正方法的有效性和正確性。

5 結(jié)論

線陣環(huán)掃式衛(wèi)星實現(xiàn)了超大幅寬和高分辨率對地觀測的兼顧，但影像幾何畸變較為嚴重。為了解決影像成像方式復雜、影像幾何畸變大的問題，本文提出了適用于線陣環(huán)掃衛(wèi)星影像的幾何糾正方法，并利用模擬影像進行了試驗驗證。該方法分為粗糾正和精糾正兩個步驟，實際應用中可以根據(jù)精度需求進行選擇使用。實驗結(jié)果表明：粗糾正精度優(yōu)于5 Pixel，精糾正精度達到1 Pixel，能夠滿足衛(wèi)星應用的需求，說明本文提出方法的正確性和有效性。下一步將研究如何在保證幾何糾正精度的同時，提高幾何糾正的處理效率，提高影像應用的時效性。