基于多遙感器的GF4交叉輻射定標

田傳陽 謝勇 宦海 邵雯

摘? 要： 在定量遙感的分析與應用中，輻射定標是基礎。文中以敦煌輻射校正場為實驗場地開展GF1，MODIS，LandSat8與GF4可見光及近紅外的交叉定標，中波紅外波段采用基于Atmospheric Infrared Sounder（AIRS）的中波紅外交叉定標。通過對4個遙感器的可見光及近紅外波段的匹配以及GF4與AIRS中波紅外波段的匹配，選取遙感器相互匹配的波段。對衛星拍攝的同時對相同區域的影像開展交叉定標，結果表明，可見光及近紅外交叉定標獲得的定標系數誤差均在5%以內，中波紅外交叉輻射定標的定標系數誤差在8%以內。通過GF4可見光及近紅外、中波紅外的定標系數誤差，表明GF4的定標精度相對較高，滿足遙感定量的需求。

關鍵詞： 交叉定標; 輻射定標; 遙感圖像; 遙感器; 波段匹配; 定標精度

中圖分類號： TN215?34; TP79? ? ? ? ? ? ? ? ? ?文獻標識碼： A? ? ? ? ? ? ? ? ? ? ? ?文章編號： 1004?373X（2020）10?0001?05

GF4 cross radiometric calibration based on multiple remote sensors

TIAN Chuanyang1， XIE Yong2，3， HUAN Hai1， SHAO Wen1

（1. School of Electronic & Information Engineering， Nanjing University of Information Science and Technology， Nanjing 210044， China;

2. School of Geography & Science， Nanjing University of Information Science and Technology， Nanjing 210044， China;

3. Nanjing Research Center， National Engineering Laboratory for Remote Sensing Satellite Applications， Nanjing 210044， China）

Abstract： The radiometric calibration is the basis for the analysis and application of quantitative remote sensing. The GF4 cross calibration in near infrared and visible light wavebands are carried out in the Dunhuang radiometric calibration field， in which GF1， moderate?resolution imaging spectroradiometer （MODIS） and LandSat8 are taken as the reference satellites， and the mid?wave infrared cross calibration in the mid?wave infrared waveband is based on atmospheric infrared sounder （AIRS）. The matching wavebands of the remote sensors are selected by matching the visible light and near?infrared wavebands of the four remote sensors， and the mid?wave infrared wavebands of GF4 and AIRS. The cross?calibration of the image taken by the satellites in the same time and area is carried out. The results show that the calibration coefficient errors obtained by cross?calibration of visible light and near?infrared wavebands are all within 5%， and the cross radiometric calibration of mid?wave infrared wavebands is within 8%. The calibration coefficient errors of GF4 at visible light， near infrared and mid?wave infrared wavebands indicate that the calibration accuracy of GF4 is relatively high， which can meet the requirements of remote sensing quantification.

Keywords： cross calibration; radiometric calibration; remote sensing image; remote sensor; wavebands matching; calibration precision

“高分四號”（GF4）于2015年成功發射，經過4次變軌到達預定軌道，于2016年6月開始投入使用[1]。GF4是“國家高分辨率對地觀測系統重大專項”中唯一的一顆高軌遙感衛星，軌道高度為36 000 km，具備全色、多光譜和中波紅外譜段成像的能力，搭載一臺可見光50 m分辨率、中波紅外400 m分辨率、幅寬400 km以上的凝視相機。該相機具有3種成像方式：凝視成像、區域成像、巡查成像[2]，與此同時，GF4具有長期駐留固定區域的能力。利用GF4這一特性和快速指向調整的優勢，可快速對我國及周邊地區的觀測任務進行響應，實現高時間分辨率和較高空間分辨率的結合。目前GF4主要應用于減災、林業、地震、氣象、火災等方面，并可為海洋、國土、水利等方面提供遙感數據[3]。

截至目前，GF4在軌運行已有4年之久。衛星上搭載的光學器件隨著在軌運行時間的推移以及太空環境的變化，光學元件逐漸老化[4?5]，譜段的輻射特性發生改變，發射前實驗室定標的結果已不能反映此時GF4光學器件的輻射特性。而且星上搭載定標黑體的光路結構差異、性能穩定性下降以及性能變化監測手段的匱乏等造成了星上定標結果精度受限和時間有效性的縮短，這些都導致輻射性能發生變化[6]。通過交叉定標可以隨時對傳感器進行性能監測，進而根據監測結果進行調整與改進，從而減少誤差對定標的影響。

當前國際上在軌輻射定標主要有場地定標[7]、星上定標[8]和交叉定標[9]等3種方法，官方通常以場地定標法為基準，獲取場地輻射定標系數。但由于場地定標的影響因素較多，費時費力，國外通常1年只進行兩次場地定標實驗，國內為一次。場地定標的次數有限，無法動態地監測傳感器的輻射性能變化[10]。交叉定標是使用高輻射精度的衛星來校準輻射精度較低的衛星，是目前檢查與驗證傳感器測量精度和穩定性的有效方法。現階段GF4每年只開展一次場地定標，不能動態地監測傳感器的輻射特性。而采用交叉輻射定標可以及時了解衛星傳感器的性能變化情況，有助于提高定標精度。本文選取我國的GF1以及國外的MODIS（Moderate?resolution Imaging Spectroradiometer），LandSat8，AIRS[11]作為參考衛星，開展與GF4的交叉定標。通過多顆高精度衛星對GF4進行交叉輻射定標，綜合對比多顆衛星的定標結果，全面了解GF4運行的輻射特性。

1? 基本原理

1.1? 交叉輻射定標方法

交叉輻射定標是一種定標方法，該方法是以精度較高的衛星作為參考衛星，以同時或近同時的相同影像區域為基準，進行傳感器間的光譜通道匹配，建立參考傳感器與待定標傳感器之間的線性關系，通過線性定標系數，實現交叉輻射定標。

交叉輻射定標流程如圖1所示。

交叉定標方法只需要通過影像數據即可完成定標，不需要實地觀測數據，提高了定標頻率，減少了人力物力的損耗，是目前較方便的定標方法。而且該方法還可用于對過去影像的定標，了解傳感器定標隨時間的變化。GF1，GF4，MODIS，LandSat8衛星各波段對比如表1所示。

在進行交叉定標時，MODIS和LandSat8的表觀輻射亮度公式為：

[RiMLGF1=Wigain（DNiMLGF1-Bioffset）] （1）

式中：[RiMLGF1]是MODIS，LandSat8，GF1衛星影像的第i波段的表觀輻射亮度;[Wigain]是第i波段的增益系數;[Bioffset]是第i波段的偏移量;[DNiMLGF1]是第i波段的數字計數值（MODIS的可見近紅外波段是1，2，3，4，LandSat8是波段2，3，4，5）。再根據GF4的DN（Digital Number）建立一元線性擬合，獲取定標系數，方程為：

[RjMLGF1=gainj·DNjGF4-offsetj] （2）

式中：[gainj]是第j波段的增益系數;offsetj是第j波段的偏移量;[DNjGF4]是GF4第j波段的數字計數值，j=1，2，3，4。

在進行中波紅外定標時，由于GF4中波紅外波段帶寬較寬，與GF4相匹配的AIRS有234個波段，其中有效波段有222個。本文根據GF4的中心波長為3.80 ?m，選取AIRS在3.80 ?m附近的3個波段，根據3個波段與GF4光譜響應函數的交點分別計算對應的權重，以及將3個波段擬合成1個波段與GF4進行交叉輻射定標。

權重計算公式為：

[bx=yxy1+y2+y3，? x=1，2，3] （3）

式中：x表示波段;[bx]表示第x波段的權重;y表示對應的3個光譜響應函數的交點值。

擬合輻亮度計算公式為：

[Rfitting=x=13bx·rx] （4）

式中：[Rfitting]表示三波段擬合計算的輻亮度值;[bx]表示第x波段的權重;rx表示第x波段的輻亮度值。

1.2? 研究區域

本文選取敦煌輻射校正場（93°42′56″E～95°3′35″E，39°34′8″N～40°38′28″N）作為可見光及近紅外定標場地[4]。作為我國國家級輻射校正場，敦煌地處甘肅，明顯的特點是氣候干燥、降雨量少、蒸發量大、晝夜溫差大。敦煌區域地表均勻，生物多樣性較少，輻射特性較為穩定不易變化，適合進行絕對輻射定標。中波紅外定標場選取的是河北的草場和福建的沙灘等地，都具有區域輻射特性均勻穩定的特性，適合開展定標研究。

1.3? 參考衛星的選擇

本文選取國外的MODIS和LandSat8衛星作為可見光及近紅外波段的參考衛星。在Terra和Aqua衛星上都有搭載MODIS相機，有36個波段，光譜范圍較寬[12]，單波段帶寬較窄，且具有星上定標的能力，誤差在2%左右;LandSat8衛星的定標精度為5%左右，精度較高，選取作為參考遙感器。AIRS是國際上第一個真正意義上的星載高光譜紅外大氣探測儀器。其光譜分辨率優于1 200 m，輻射測量精度優于0.2 K，使星載大氣探測儀探測能力滿足當前精度水平，適合作為參考衛星[13]。GF1衛星的遙感數據通過多種參考目標法實現寬動態定標，定標精度相對較高，精度為5%，也被選作為參考衛星[14]，進行國內衛星之間的交叉定標。

1.4? 數據的選取

在開展交叉輻射定標之前，根據可見光及近紅外波段選擇GF1，GF4，MODIS，LandSat8四個衛星的有效影像對，根據中波紅外波段選擇GF4與AIRS匹配的影像對。其中數據選取有以下幾個條件：相同區域的時間差要控制在1 h以內，時間差過長將會導致區域氣候發生變化以及光譜特性改變，對定標產生影響;要求區域上方無云，云會對光的傳輸產生影響而且云在短時間內可能會發生劇烈變化;盡量選取影像的中心區域[15]。

1.5? 光譜響應函數匹配

本文針對GF4的可見光及近紅外通道與GF1，MODIS，LandSat8進行光譜匹配，光譜匹配的關鍵環節是選擇與GF4衛星相機的4個多光譜波段相匹配的GF1，MODIS和LandSat8衛星波段。

從圖3四顆衛星的光譜響應曲線可以看出，GF1和GF4前3個波段和MODIS，LandSat8的波段匹配比較好，第4波段匹配相對有所偏移，但是GF1和GF4的波段匹配的很好。同時也可以看出，MODIS與LandSat8的波段相對匹配較好[16]，帶寬相對國產衛星的較窄。

而GF4中波紅外通道與AIRS進行波段匹配時，由圖3可知在GF4中波紅外通道范圍內，AIRS有234個波段（12個波段損壞）。由于匹配通道較多，本文選取GF4中心波長處的3個波段。

1.6? 數據的處理

由于GF4/PMS的分辨率為50 m，MODIS的分辨率為250 m，LandSat8的空間分辨率為30 m，GF1的分辨率為16 m，因此需要統一4個傳感器的空間分辨率。通過3次內插方法把GF4，GF1，MODIS，LandSat8空間分辨率統一降為250 m。在進行中波紅外定標時將GF4空間分辨率降為1 000 m，在此基礎上與AIRS進行定標。本文采用SIFT算法進行GF4與MODIS，LandSat8，GF1以及AIRS的數據匹配，通過該算法選取極值點（具有極強的局部穩定性的點）作為影像對的同名點。

2? 實驗結果

依據式（1）、式（2）計算出GF1，MODIS，LandSat8，AIRS的輻亮度值以及傳感器相應的增益和偏移量。圖4～圖6是GF4衛星和MODIS，LandSat8，GF1的交叉輻射定標擬合圖。圖中，x軸表示GF4的DN值，y軸表示MODIS、LandSat8、GF1的表觀輻射亮度值，擬合系數及誤差如表2所示。

定標 增益 誤差 /% 增益 誤差 /% 增益 誤差 /% Band2 0.174 9 0.170 5 2.50 0.177 5 1.49 0.180 5 3.20 Band3 0.190 3 0.182 3 4.20 0.192 2 0.99 0.187 2 1.63 Band4 0.153 2 0.148 6 3.00 0.157 4 2.74 0.160 2 4.57 Band5 0.107 3 0.103 4 3.58 0.110 3 2.80 0.105 7 1.52 ]

圖7是GF4與AIRS在河北等地的交叉輻射定標結果。在圖7中，x軸表示GF4 DN值，y軸表示三波段的擬合輻亮度值。從表2中可以看出，通過GF4與MODIS，LandSat 8，GF1交叉定標結果，可以得到以下結論：

2） 比較GF4?MODIS，GF4?GF1和GF4?LandSat8交叉定標獲得的定標系數及誤差可知，以LandSat8為參考衛星的可見光及近紅外的定標結果最好，且都在3%以內。

3） 通過GF1衛星對GF4衛星進行定標，結果與官方定標系數相比誤差分別是3.2%，1.63%，4.57%，1.52%。并由圖6可以看出，GF4和GF1的光譜匹配效果比其他兩顆衛星要好，且都是國產衛星沒有跨時區的問題，適合開展我國自主同類衛星遙感器間的交叉定標研究。

4） 從圖7可知GF4中波紅外波段的定標誤差分別為6.6%，1.3%，1.8%，7.8%，表明GF4中波紅外波段的精度是比較高的。

3? 結? 論

通過對GF4可見光及近紅外、中波紅外波段的交叉輻射定標研究，結果表明GF4的輻射定標精度是比較高的。無論是可見光及近紅外波段還是中波紅外波段，定標精度都滿足定量遙感的需求。同時，通過GF1對GF4的定標結果來看，GF1的定標精度在5%以內，而且GF1與GF4波段匹配較好，且在相同的時區，適合開展我國自主同類衛星遙感器間的交叉定標研究。目前為止，針對GF4與AIRS的中波紅外定標研究還有許多不足，本文只是對AIRS的222個匹配波段中的3個波段，做了一個擬合定標，并未對所有匹配波段進行細致的研究，接下來將會對AIRS的222個有效波段進行總體的分析與研究。

參考文獻

[1] 練敏隆，石志城，王躍，等.“高分四號”衛星凝視相機設計與驗證[J].航天返回與遙感，2016（4）：32?39.

[2] 馬文坡，練敏隆.“高分四號”衛星凝視相機的技術特點[J]，航天返回與遙感，2016（4）：26?31.

[3] 劉鳳晶，李果，于登云，等.高分四號衛星及應用概況[J].衛星應用，2018（12）：12?18.

[4] 丁闖，謝勇，宦海，等.基于GF?1/WFV與MODIS和LandSat8相機交叉定標[J].現代電子技術，2019，42（8）：37?41.

[5] 胡曄，尤建潔，戴奇燕，等.圖像邊緣信號能量與MTF的相關研究[J].航天返回與遙感，2005（4）：20?24.

[6] 韓杰，謝勇.GF?1衛星WFV影像間勻色方法[J].測繪學報，2016（12）：1423?1433.

[7] 周珂，劉李，余濤，等.光譜匹配因子對GF?1/WFV時間序列交叉定標的影響分析[J].光譜學與光譜分析，2017，37（12）：3809?3813.

[8] XIONG X X， SUN J Q， BARNES W， et al. Multiyear on?orbit calibration and performance of terra MODIS reflective solar bands [J]. IEEE transactions on geoscience and remote sensing， 2007， 45（4）： 879?889.

[9] 韓啟金，潘志強，張學文.基于均勻目標的GF?1WFV與SJ?9A PMS相機交叉定標[J].光學學報，2014（6）：292?299.

[10] 鐘雯嬌，謝勇，宦海，等.GF?4可見光及近紅外譜段的輻射定標[J].現代電子技術，2019，42（8）：113?117.

[11] 劉李，傅俏燕，史婷婷，等.利用AIRS高光譜數據開展HJ?1B熱紅外通道的多點交叉定標及驗證[J].中國科學：技術科學，2015（1）：103?110.

[12] 張玉環，毛慧琴，王中挺，等.基于光譜與幾何匹配的GOCI與MODIS交叉輻射定標[J].大氣與環境光學學報，2016（6）：412?422.

[13] 劉李，顧行發，余濤，等.HJ?1B衛星熱紅外通道在軌場地定標與驗證[J].紅外與激光工程，2012（5）：1119?1125.

[14] 杜麗麗，易維寧，王昱，等.“天繪一號”衛星多傳感器協同輻射定標方法[J].光學學報，2019（4）：42?49.

[15] 殷麗梅，喬兵，劉俊池，等.地基紅外輻射測量系統聯合輻射定標法[J].光學學報，2018（4）：209?214.

[16] 呂佳彥，何明元，陳林，等.基于敦煌輻射校正場的自動化輻射定標方法[J].光學學報，2017（8）：25?32.