何靈莉 胡秀清 王玲 陳林 徐娜 張鵬

基于自動(dòng)化方法的光學(xué)載荷定標(biāo)跟蹤

何靈莉1,3胡秀清2,3*王玲2,3陳林2,3徐娜2,3張鵬2,3

（1 中國氣象科學(xué)研究院，中國氣象局，北京 100081）（2 中國遙感衛(wèi)星輻射測量和定標(biāo)重點(diǎn)開放實(shí)驗(yàn)室，中國氣象局，北京 100081）（3 國家衛(wèi)星氣象中心，中國氣象局，北京 100081）

為了實(shí)現(xiàn)對光學(xué)傳感器在軌性能的連續(xù)監(jiān)測，文章使用敦煌場地的自動(dòng)化觀測數(shù)據(jù)對“風(fēng)云三號(hào)”衛(wèi)星可見光紅外掃描輻射計(jì)（Visible and Infrared Radiometer，VIRR）近三年的數(shù)據(jù)進(jìn)行了連續(xù)定標(biāo)跟蹤，推算VIRR各通道的定標(biāo)斜率，基于敦煌輻射校正場自動(dòng)化觀測的定標(biāo)相較于傳統(tǒng)現(xiàn)場試驗(yàn)定標(biāo)方法具有明顯提高定標(biāo)頻次優(yōu)勢。采用搭載于Aqua衛(wèi)星的高精度的中分辨率成像光譜儀（Aqua Moderate Resolution Imaging Spectroradiometer，Aqua MODIS）觀測對自動(dòng)觀測定標(biāo)方法進(jìn)行精度驗(yàn)證，表明該方法在可見光近紅外（Visible and Near Infrared，VNIR）具有3%的定標(biāo)精度，在短波紅外（Short Wave Infrared，SWIR）具有5%的定標(biāo)精度。將該方法的結(jié)果與業(yè)務(wù)多場地方法的結(jié)果進(jìn)行比較，在3年尺度上均值相對偏差整體低于2%。對VIRR的跟蹤結(jié)果表明：基于場地自動(dòng)化觀測能夠?qū)崿F(xiàn)長期定標(biāo)監(jiān)測，可將該方法推廣應(yīng)用于其他衛(wèi)星光學(xué)載荷以真正實(shí)現(xiàn)多載荷的場地自動(dòng)化定標(biāo)業(yè)務(wù)化。

可見光紅外掃描輻射計(jì) 太陽反射波段 自動(dòng)化輻射定標(biāo) 反射率基法 “風(fēng)云三號(hào)”衛(wèi)星 航天遙感

0 引言

對地遙感觀測衛(wèi)星的傳感器發(fā)射進(jìn)入預(yù)定軌道后，隨著時(shí)間的推移，復(fù)雜惡劣的太空環(huán)境和元器件老化會(huì)導(dǎo)致傳感器探測性能衰減。在沒有星上定標(biāo)設(shè)備的情況下，最常用的替代定標(biāo)方法就是采用地面穩(wěn)定目標(biāo)場對傳感器進(jìn)行絕對輻射定標(biāo)。搭載于我國“風(fēng)云三號(hào)”系列極軌氣象衛(wèi)星（FY-3）的可見光紅外掃描輻射計(jì)（VIRR）沒有配備可見光星上定標(biāo)設(shè)備，它的太陽反射波段定標(biāo)主要采用敦煌輻射定標(biāo)場的場地替代定標(biāo)方法[1]。該方法是早期“風(fēng)云”衛(wèi)星太陽波段的基線定標(biāo)方法，一般每年夏季在敦煌進(jìn)行場地定標(biāo)試驗(yàn)[2]，定標(biāo)精度約5%~6%，同樣的方法計(jì)算的中分辨率成像光譜儀（Moderate Resolution Imaging Spectroradiometer，MODIS）表觀反射率，偏差約3%左右[3]。但是這種一年一次的定標(biāo)試驗(yàn)滿足不了可見光紅外掃描輻射計(jì)（Visible and Infrared Radiometer，VIRR）輻射性能的在軌漂移連續(xù)性監(jiān)測和評估。

基于提高定標(biāo)頻次和降低定標(biāo)成本的考慮，亞利桑那大學(xué)遙感組（Remote Sensing Group，RSG）在2004年提出了自動(dòng)化的輻射定標(biāo)試驗(yàn)場（Radiometric Calibration Test Site，RadCaTS）概念，在沒有人員參與的情況下對定標(biāo)場的地表反射率和大氣參數(shù)進(jìn)行自動(dòng)化觀測，從而實(shí)現(xiàn)對空載和星載傳感器進(jìn)行輻射定標(biāo)。RSG在內(nèi)華達(dá)州的鐵路谷試驗(yàn)場對RadCaTS進(jìn)行了驗(yàn)證，將計(jì)算得到的大氣層頂（Top of Atmosphere，TOA）光譜輻亮度與MODIS觀測進(jìn)行比較偏差小于3.6%[4]。其后，RSG致力于不斷精細(xì)化和完善RadCaTS定標(biāo)的各個(gè)環(huán)節(jié)[5-8]。使用傳統(tǒng)跑場方法和RadCaTS同時(shí)對業(yè)務(wù)陸地成像儀（Operational Land Imager，OLI）、“地球觀測1號(hào)”Hyperion和先進(jìn)陸地成像儀（Advanced Land Imager，ALI）進(jìn)行絕對輻射定標(biāo)，輻射穩(wěn)定性和精度評價(jià)，RadCaTS方法整體優(yōu)于人工測量方法[9-10]。RSG和噴氣推進(jìn)實(shí)驗(yàn)室的研究均證實(shí)了自動(dòng)化輻射定標(biāo)方法對于星載傳感器太陽反射波段定標(biāo)的可行性，這為后來自動(dòng)化輻射定標(biāo)技術(shù)的進(jìn)一步應(yīng)用和布局提供了經(jīng)驗(yàn)參考。經(jīng)過十幾年的發(fā)展，世界上基于裝備自動(dòng)化觀測儀器的輻射定標(biāo)場網(wǎng)（Radiometric Calibration Network，RadCalNet）初步形成，該定標(biāo)場網(wǎng)是對LANDNET定標(biāo)場網(wǎng)概念的繼承和延伸，它在數(shù)據(jù)品質(zhì)和場地管理上更加的系統(tǒng)和流程化，并為衛(wèi)星運(yùn)營方提供可溯源至國際單位SI的TOA光譜反射率，中國包頭定標(biāo)場是其成員之一，敦煌場正在加入過程中。

我國場地自動(dòng)化觀測定標(biāo)始于2014年，李新等對國外自動(dòng)化定標(biāo)方法和技術(shù)的現(xiàn)狀進(jìn)行分析，提出了我國自動(dòng)化定標(biāo)設(shè)備的設(shè)計(jì)方案和關(guān)鍵技術(shù)[11]；隨后，尹亞鵬等對場地自動(dòng)化試驗(yàn)場地輻射計(jì)（Automatic Test Radiometer，ATR）進(jìn)行了設(shè)計(jì)，實(shí)驗(yàn)室定標(biāo)和基于敦煌場地的試驗(yàn)，并初步驗(yàn)證了輻射計(jì)的設(shè)計(jì)合理性[12]；李新等基于場地試驗(yàn)結(jié)果驗(yàn)證了ATR的觀測性能可靠，隨著對場地ATR的不斷驗(yàn)證和改進(jìn)，其性能不斷提升，成為了敦煌自動(dòng)化輻射定標(biāo)體系的核心設(shè)備[13]；邱剛剛、呂佳彥等均以MODIS的觀測為輻射基準(zhǔn)，對基于敦煌場的自動(dòng)化定標(biāo)技術(shù)進(jìn)行了精度驗(yàn)證[14-15]。結(jié)果表明該項(xiàng)技術(shù)能夠在保證傳統(tǒng)定標(biāo)方法精度的基礎(chǔ)上有效增加定標(biāo)頻次和降低定標(biāo)成本。由于單一ATR觀測對整個(gè)敦煌場地表反射率的表征有限，安徽光學(xué)精密儀器研究所于2018年4月又在敦煌場地布設(shè)了兩臺(tái)新的ATR，張孟等基于2018年5月~11月的新ATR自動(dòng)化觀測數(shù)據(jù)實(shí)現(xiàn)了對可見光紅外成像儀（Visible Infrared Imaging Radiometer Suite，VIIRS）的高頻次在軌輻射定標(biāo)試驗(yàn)[16]，但目前為止尚沒有采用敦煌場長期自動(dòng)化觀測對中國遙感衛(wèi)星開展長時(shí)間連續(xù)定標(biāo)和評估。本研究主要是基于敦煌場2016年9月以來近三年的數(shù)據(jù)開展FY-3C的VIRR開展長期自動(dòng)定標(biāo)試驗(yàn)。

多場地定標(biāo)方法基于具有不同亮度等級(jí)的地表穩(wěn)定目標(biāo)和輻射傳輸模式對衛(wèi)星傳感器進(jìn)行定標(biāo)，孫凌等使用該方法對FY-3A中分辨率光譜成像儀（Medium Resolution Spectral Imager，MERSI）進(jìn)行了定標(biāo)跟蹤[17]；王玲等使用該方法對FY-3A、FY-3B、FY-3C VIRR進(jìn)行了再定標(biāo)[18]，均表明該方法具有高頻和高精度優(yōu)勢。

本文基于2016年9月~2019年6月的敦煌場自動(dòng)化觀測數(shù)據(jù)，使用自動(dòng)化輻射定標(biāo)方法對FY-3C VIRR進(jìn)行了長時(shí)間序列輻射定標(biāo)，并將所得的定標(biāo)結(jié)果與多場地定標(biāo)方法的結(jié)果進(jìn)行比較。

1 場地與儀器介紹

1.1 敦煌場地及觀測儀器

敦煌場（東經(jīng)94.27°，北緯40.18°）是中國輻射校正場地之一，坐落在中國甘肅省敦煌市中心西部35km處。整個(gè)場地約30km×30km，位于穩(wěn)定的沖積扇中，地表沒有植被，主要由混合的砂礫和黏土等構(gòu)成。在場地中心的10km×10km區(qū)域內(nèi)（東經(jīng)94.2°~東經(jīng)94.3°，北緯40.05°~北緯40.15°），表征場地空間均勻性和時(shí)間穩(wěn)定性的變異系數(shù)（標(biāo)準(zhǔn)差/平均值）分別小于2%和3%；可見光近紅外波段的地表反射率大約在15%~30%；場地地表具有非朗伯特性；場地大氣干潔氣溶膠含量較低，除沙塵季節(jié)外，550nm的平均光學(xué)厚度小于0.2。敦煌場的上述各項(xiàng)特性，使其成為較為理想的替代定標(biāo)場地。2008年，敦煌場被地球觀測衛(wèi)星委員會(huì)（CEOS）的定標(biāo)驗(yàn)證工作組（WGCV）選為了儀器參考場地[3]。

2016年9月，在敦煌定標(biāo)場的試驗(yàn)觀測基地樣地里布設(shè)了1臺(tái)自動(dòng)化觀測的通道式輻射計(jì) ATR00，其波段覆蓋了可見光近紅外波段的8個(gè)通道（405nm、450nm、555nm、610nm、700nm、808nm、980nm、1 540nm），光學(xué)頭部位于1.8m高的支架上，采用天底觀測方式獲取正下方直徑為30cm圓形區(qū)域內(nèi)的上行輻射亮度，視場角為10°，每三分鐘采集一次數(shù)據(jù)。2018年4月在敦煌場中心布設(shè)了另外一臺(tái)儀器ATR01。

1.2 VIRR簡介

可見光紅外掃描輻射計(jì)VIRR，具有10個(gè)波段，其中7個(gè)太陽反射波段覆蓋了從0.4~1.65μm的波段范圍，3個(gè)熱發(fā)射通道覆蓋了3.5~12.5μm的波段范圍（如表1所示）。VIRR進(jìn)行對地觀測時(shí)使用45°的掃描鏡，采用跨軌掃描方式，對地掃描角為±55.1°，星下點(diǎn)空間分辨率為1 100m，幅寬在跨軌方向?yàn)? 800km，較大的幅寬使其能夠在一天內(nèi)對全球進(jìn)行覆蓋。

2 場地均勻性驗(yàn)證

ATR對直徑為30cm的圓形區(qū)域進(jìn)行采樣，局部區(qū)域地表反射率對整個(gè)場地反射率的可表征性需要進(jìn)行驗(yàn)證。可以通過使用MODIS衛(wèi)星觀測對場地空間均勻性進(jìn)行評估，并對比相距較遠(yuǎn)的ATR00和ATR01采集的地表反射率，雙管齊下地評價(jià)局部觀測對場地整體的可表征性。

對場地空間均勻性進(jìn)行評價(jià)必須考慮尺度效應(yīng)，在不同的尺度下評價(jià)均勻性會(huì)有不同的結(jié)果。為適應(yīng)VIRR 1.1km的空間分辨率評價(jià)需求，使用1km空間分辨率的MODIS地表反射率產(chǎn)品MOD09對敦煌場的空間均勻性進(jìn)行評價(jià)，評價(jià)使用空間變異系數(shù)（Coefficient of Variation，CV）作為指標(biāo)，計(jì)算時(shí)的滑動(dòng)窗口為10km×10km。在可見光近紅外波段，除了通道3（469nm）的空間變異系數(shù)在10%以上，通道4（555nm）在4%左右，其余通道的空間變異系數(shù)均小于3%（如圖1）。在進(jìn)行輻射定標(biāo)場的空間均勻性評價(jià)中，通常以空間變異系數(shù)小于3%認(rèn)定場地的均勻性[19-20]。通道3表現(xiàn)出的不均勻性可能是大氣強(qiáng)散射作用造成的。

表1 VIRR通道特性

Tab.1 The VIRR channel characteristics

注：3、4、5通道的靈敏度單位含義為300K的環(huán)境下變化溫度。

圖1 0.858μm和1.640μm空間均勻性評價(jià)結(jié)果

選取2018年4月~2019年4月滿足定標(biāo)要求的ATR00和ATR01觀測數(shù)據(jù)，計(jì)算ATR通道地表反射率并進(jìn)行比較（圖2），通道1~3和通道8的百分偏差整體低于5%，年平均百分偏差在2%~3.8%之間；通道4~7整體高于5%，年均值在6%~8%。百分偏差與季節(jié)變化相關(guān)，敦煌場在8~10月份云量較少，大氣較為穩(wěn)定，最為適合定標(biāo)，因此偏差相對較小，定標(biāo)結(jié)果較為穩(wěn)定，波動(dòng)性較小。

綜上，敦煌定標(biāo)場空間均勻，且ATR00與ATR01觀測的地表反射率相對百分偏差在通道1~3和8小于5%，表明ATR00觀測的地表反射率能夠以小于5%的精度表征以ATR00為中心10km×10km區(qū)域的反射率。雖然通道4~7具有6%~8%的偏差，但可以減小這幾個(gè)通道在后續(xù)光譜平移中的權(quán)重以消除影響。

圖2 ATR00相對于ATR01地表反射率的長序列（2018年4月~2019年4月）百分偏差

3 自動(dòng)化輻射定標(biāo)方法

場地自動(dòng)化觀測的輻射定標(biāo)方法與傳統(tǒng)的場地地表試驗(yàn)沒有本質(zhì)差異，仍舊采用反射率基法，包括地表反射率計(jì)算，大氣參數(shù)獲取，輻射傳輸計(jì)算等關(guān)鍵環(huán)節(jié)，自動(dòng)化輻射定標(biāo)技術(shù)流程如圖3所示。其中MYD08_D3為MODIS的大氣參數(shù)日產(chǎn)品，TOA表示大氣層頂，BRDF（Bidirectional Reflection Distribution Function）為雙向反射分布函數(shù)，6S為Second Simulation of Satellite Signal in the Solar Spectrum的簡寫。

圖3 自動(dòng)化輻射定標(biāo)處理流程圖

3.1 數(shù)據(jù)品質(zhì)控制及篩選

數(shù)據(jù)品質(zhì)控制主要包括以下三個(gè)方面：判斷衛(wèi)星觀測是否被云、積雪等影響；輔助大氣參數(shù)是否滿足適用條件；自動(dòng)化儀器觀測數(shù)據(jù)是否異常。要求VIRR衛(wèi)星數(shù)據(jù)觀測天頂角小于30°，定位誤差小于1km，觀測點(diǎn)10像元×10像元窗口的變異系數(shù)小于3%，較大的變異系數(shù)表明可能有云、沙塵暴、積雪等影響空間均勻性。較高的氣溶膠含量會(huì)增加輻射傳輸計(jì)算的不確定性，若氣溶膠含量高于0.2便不適于定標(biāo)。ATR是自動(dòng)化輻射定標(biāo)方法的核心儀器，其數(shù)據(jù)品質(zhì)對定標(biāo)精度有較大影響，其數(shù)據(jù)品質(zhì)控制主要包括剔除負(fù)值、跳變等，如圖4所示（a）是ATR的原始計(jì)數(shù)值，可以發(fā)現(xiàn)具有非常明顯的數(shù)據(jù)跳變，（b）是進(jìn)行品質(zhì)控制后的數(shù)據(jù)，線型連續(xù)且平滑。基于數(shù)據(jù)品質(zhì)控制后的ATR觀測數(shù)據(jù)進(jìn)行可定標(biāo)性判定，其思路是利用衛(wèi)星過境前后1小時(shí)的ATR觀測數(shù)據(jù)波動(dòng)性判斷大氣狀況穩(wěn)定性。大氣質(zhì)量是指光穿過大氣的路徑長度，大氣不穩(wěn)定會(huì)破壞ATR觀測計(jì)數(shù)值的對數(shù)值與大氣質(zhì)量之間的強(qiáng)線性關(guān)系，可以通過二者線性關(guān)系的強(qiáng)弱程度判斷當(dāng)日是否適合定標(biāo)。

圖4 數(shù)據(jù)品質(zhì)控制過程

3.2 計(jì)算ATR通道等效地表反射率

對衛(wèi)星過境前后10min的ATR觀測進(jìn)行輻射校正獲得地表上行輻射亮度，利用6S大氣輻射傳輸模式計(jì)算總下行輻照度，模式運(yùn)算所需要的氣溶膠、水汽和臭氧含量等大氣參數(shù)來自于Aqua MODIS的全球日平均產(chǎn)品MYD08_D3，地表反射率的計(jì)算可表達(dá)為式（1）

3.3 獲取高光譜地表反射率

由ATR通道地表反射率推導(dǎo)衛(wèi)星觀測方向高光譜地表反射率，必須首先明確一個(gè)基本的前提：敦煌場地的地表特性具有時(shí)間穩(wěn)定性，即隨著時(shí)間變化，地表反射率光譜只產(chǎn)生幅度的變化，而形狀不發(fā)生變化，這是基于通道儀器觀測推算高光譜反射率的一個(gè)假定條件[21]?；谠撉疤?，可以通過平移先驗(yàn)地表反射率光譜擬合ATR推算的通道地表反射率，獲得平移光譜。由于敦煌場地不是朗伯表面，反射率具有方向特性，需要利用BRDF將天底觀測的平移光譜校正到衛(wèi)星觀測方向獲得衛(wèi)星觀測方向光譜（見圖5）。這里使用的BRDF模型是中國氣象局國家衛(wèi)星氣象中心基于2008年夏季的場地實(shí)測數(shù)據(jù)和Ross-Li驅(qū)動(dòng)核構(gòu)建的[3]。

圖5 ATR多通道地表反射率擬合推導(dǎo)衛(wèi)星觀測方向地表反射率光譜

3.4 計(jì)算定標(biāo)斜率

計(jì)算定標(biāo)斜率需要將衛(wèi)星觀測方向的地表光譜卷積為VIRR通道等效地表反射率，并輸入6S輻射傳輸模式，計(jì)算獲得VIRR通道等效TOA反射率TOA。定標(biāo)斜率的計(jì)算如式（2）

求取定標(biāo)系數(shù)涉及衛(wèi)星和地面觀測的空間匹配。一般來講，定標(biāo)精度同具體的定標(biāo)點(diǎn)選取相關(guān)[22]，本文中采用ATR00觀測地表反射率正演的VIRR通道等效TOA反射率與以ATR01為中心10km×10km區(qū)域的VIRR衛(wèi)星觀測進(jìn)行匹配。中心像元的經(jīng)緯度定位誤差要小于1km，窗口的變異系數(shù)需小于3%。

4 精度驗(yàn)證

本研究使用自動(dòng)化輻射定標(biāo)方法在近三年的時(shí)間跨度上對VIRR進(jìn)行定標(biāo)跟蹤，該方法的不確定性可能來自于ATR儀器性能衰減、輻射傳輸模式計(jì)算、BRDF模型在敦煌場地的季節(jié)變化[23]等。為了驗(yàn)證本研究方法的精度和不確定性，使用搭載于Aqua衛(wèi)星的MODIS的1級(jí)產(chǎn)品對方法整體精度進(jìn)行獨(dú)立驗(yàn)證?；谠摲椒ㄓ?jì)算Aqua MODIS的表觀反射率相對于衛(wèi)星觀測值的百分偏差。從表2可以發(fā)現(xiàn)在VNIR具有3%的定標(biāo)精度，在SWIR具有5%的定標(biāo)精度，標(biāo)準(zhǔn)差范圍從1.95%到3.67%。

表2 MODIS驗(yàn)證結(jié)果

Tab.2 The MODIS validation results %

5 結(jié)果與討論

5.1 時(shí)間序列定標(biāo)斜率

利用2016年9月~2019年6月的敦煌場自動(dòng)化觀測數(shù)據(jù)對FY-3C VIRR太陽反射波段進(jìn)行定標(biāo)，通道10受水氣影響較大，此處不納入討論。自動(dòng)化輻射定標(biāo)的結(jié)果表現(xiàn)出不平穩(wěn)性是因?yàn)槎?biāo)鏈路中不確定性來源較多，如：自動(dòng)化儀器觀測和大氣參數(shù)的準(zhǔn)確性、地表滯留水等，長時(shí)序變化突顯了這種不確定性。雖然衛(wèi)星傳感器的輻射性能會(huì)隨時(shí)間產(chǎn)生一定程度的衰減，但在短時(shí)間內(nèi)卻可以假設(shè)它保持不變。基于上述假設(shè)以可定標(biāo)天為中心對前后15天之內(nèi)的結(jié)果進(jìn)行平均，平滑后的定標(biāo)斜率時(shí)序變化如圖6所示。在近3年的時(shí)間尺度上，F(xiàn)Y-3C VIRR在敦煌場共有753次觀測，最終可定標(biāo)100天。敦煌自動(dòng)化輻射定標(biāo)方法的頻次主要受天氣狀況制約，冬季受降雪影響，春季受高發(fā)沙塵暴影響，夏秋季節(jié)云少大氣穩(wěn)定，因而具有較高的定標(biāo)頻次。各通道定標(biāo)斜率隨著時(shí)間變化趨勢一致，存在季節(jié)性波動(dòng)，但變化趨勢整體平穩(wěn)，顯示出FY-3C VIRR輻射響應(yīng)性能的穩(wěn)定，表明自動(dòng)化輻射定標(biāo)方法具有對VIRR傳感器進(jìn)行長時(shí)間序列定標(biāo)跟蹤的能力。

圖6 FY-3C VIRR的定標(biāo)斜率時(shí)序變化圖

對于敦煌場，季節(jié)性的降雪、沙塵等會(huì)對地表均勻性產(chǎn)生影響，其BRDF也存在季節(jié)性變化[23]，季節(jié)性沙塵還會(huì)造成氣溶膠含量的季節(jié)性變化，當(dāng)自動(dòng)化輻射定標(biāo)算法對這些因素控制不夠全面時(shí)，結(jié)果就呈現(xiàn)出一定的季節(jié)性。通道1和2相較于其他通道波動(dòng)性更明顯，因?yàn)檫@兩個(gè)通道對云等更加敏感。本研究選取ATR00測量點(diǎn)的觀測與以ATR01為中心的10km×10km區(qū)域進(jìn)行匹配，因?yàn)锳TR00與場地有一定的距離，匹配時(shí)對空間均勻性有一定要求，理論上選取ATR01測量點(diǎn)更為合理，但對于跟蹤衛(wèi)星傳感器的輻射性能來講，使用具有更長時(shí)間序列可用數(shù)據(jù)的ATR00更有意義，如能夠使用10km×10km區(qū)域內(nèi)的多臺(tái)ATR觀測平均，想必能夠進(jìn)一步提高定標(biāo)精度與定標(biāo)穩(wěn)定性。前文在1km的空間分辨率下對10km×10km的區(qū)域進(jìn)行空間均勻性驗(yàn)證，結(jié)果表明敦煌場足夠均勻，能夠以局部觀測表征場地整體的平均反射率。由于尺度效應(yīng)的存在，對于具有更高空間分辨率的衛(wèi)星傳感器，不均勻性會(huì)增加，因此需要重新驗(yàn)證場地均勻性。

5.2 與多場地定標(biāo)方法的比較

多場地定標(biāo)方法基于多個(gè)地球穩(wěn)定目標(biāo)對傳感器進(jìn)行定標(biāo)，該方法同樣是對反射率基法的應(yīng)用，在假設(shè)傳感器定標(biāo)系數(shù)在短時(shí)間內(nèi)（如30天）不變的前提下，模擬全球多個(gè)場地的TOA反射率，并與衛(wèi)星實(shí)際觀測TOA反射率進(jìn)行線性回歸，得到目標(biāo)傳感器的定標(biāo)系數(shù)。它與本文方法都具有較高的定標(biāo)頻次。

圖7是FY-3C VIRR自動(dòng)化輻射定標(biāo)方法與多場地定標(biāo)方法的時(shí)間序列定標(biāo)斜率對比結(jié)果。表3是自動(dòng)化方法相對多場地方法的百分偏差統(tǒng)計(jì)結(jié)果。多場地定標(biāo)方法使用不同定標(biāo)場地30天內(nèi)的TOA反射率來計(jì)算定標(biāo)系數(shù)，而自動(dòng)化輻射定標(biāo)方法使用敦煌單場地單天的觀測數(shù)據(jù)，所以自動(dòng)化方法的斜率沒有多場地定標(biāo)斜率平滑，按天比較兩種方法的結(jié)果具有較大差異，但在整個(gè)時(shí)間尺度上的平均值卻具有較好的一致性。除通道8（0.505μm）和9（0.555μm）在5%左右，其余通道平均相對百分偏差在2%。

表3 FY-3C VIRR多場地與自動(dòng)化定標(biāo)方法比較的統(tǒng)計(jì)結(jié)果

Tab.3 Statistical results of the comparison between the multi-site and automated radiometric calibration methods of FY-3C VIRR %

6 結(jié)束語

本文采用基于反射率基法的場地自動(dòng)化觀測輻射定標(biāo)方法對FY-3C VIRR進(jìn)行了2016年9月~2019年6月的連續(xù)定標(biāo)跟蹤，結(jié)果顯示平均每7天能實(shí)現(xiàn)一次有效定標(biāo)，有效定標(biāo)主要集中在夏秋季節(jié)，部分通道的定標(biāo)系數(shù)具有季節(jié)波動(dòng)性。與傳統(tǒng)場地定標(biāo)方法相比該方法極大地提高了定標(biāo)頻次，節(jié)約了人力財(cái)力物力。使用Aqua MODIS觀測對本文所使用的方法進(jìn)行了精度驗(yàn)證，大部分通道精度在3%以內(nèi)。將FY-3C VIRR的自動(dòng)化定標(biāo)結(jié)果與多場地定標(biāo)結(jié)果相比較，相對百分偏差除8、9通道外，其余均在2%以下。以上結(jié)果表明：敦煌場自動(dòng)化輻射定標(biāo)方法具有在可靠精度下對衛(wèi)星光學(xué)載荷探測性能進(jìn)行長時(shí)間序列定標(biāo)和跟蹤的能力，可以推廣應(yīng)用于其它同類型衛(wèi)星載荷。跟蹤結(jié)果顯示出一定的季節(jié)性，原因可能來自諸多方面，如：場地光學(xué)特性的季節(jié)變化、數(shù)據(jù)品質(zhì)控制不夠嚴(yán)格等等，具體還需做進(jìn)一步分析。

致 謝

感謝中國科學(xué)院安徽光學(xué)精密機(jī)械研究所張艷娜、韋瑋等提供的數(shù)據(jù)支持和幫助。

[1] LI Y, RONG Z, ZHENG Z. et al. Post Launch Site Calibration of Visible and Near-infrared Channels of FY-3A Visible and Infrared Radiometers[J]. Optics and Precision Engineering, 2009, 17(12): 2966-2973.

[2] HU X, LIU J, SUN L, et al. Characterization of CRCS Dunhuang Test Site and Vicarious Calibration Utilization for Fengyun (FY) Series Sensors[J]. Canadian Journal of Remote Sensing, 2010, 36(5): 566-582.

[3] SUN L, GUO M, XU N, et al. On-orbit Response Variation Analysis of FY-3 MERSI Reflective Solar Bands Based on Dunhuang Site Calibration[J]. Spectroscopy and Spectral Analysis, 2012, 32(7): 1869-1877.

[4] CZAPLAMYERS J S, THOME K J, LEISSO N P. Radiometric Calibration of Earth-observing Sensors Using an Automated Test Site at Railroad Valley, Nevada[J]. Canadian Journal of Remote Sensing, 2010, 36(5): 474-487.

[5] CZAPLAMYERS J S, THOME K J, BIGGAR S F. Design, Calibration, and Characterization of a Field Radiometer Using Light-emitting Diodes as Detectors[J]. Applied Optics, 2009, 47(36): 6753-6762.

[6] LEISSO N, CZAPLAMYERS J S. Comparison of Diffuse Sky Irradiance Calculation Methods and Effect on Surface Reflectance Retrieval from an Automated Radiometric Calibration Test Site[J]. Proceedings of SPIE, 2011, 8153(2): 815310-815310-11.

[7] ANDERSON N, CZAPLAMYERS J S, LEISSO N, et al. Design and Calibration of Field Deployable Ground-viewing Radiometers[J]. Applied Optics, 2013, 52(2): 231-240.

[8] ANDERSON N, CZAPLAMYERS J S. Ground Viewing Radiometer Characterization, Implementation and Calibration Applications—A Summary After Two Years of Field Deployment[J]. Proceedings of SPIE, 2013, 8866: 88660N-88660N-10.

[9] CZAPLAMYERS J S, ANDERSON N, THOME K J, et al. The Absolute Radiometric Calibration of the Landsat 8 Operational Land Imager Using the Reflectance-based Approach and the Radiometric Calibration Test Site (RadCaTS)[J]. Proceedings of SPIE, 2014, 9218: 921819-921819-9.

[10] CZAPLAMYERS J S, ONG L, THOME K J, et al. Validation of EO-1 Hyperion and Advanced Land Imager Using the Radiometric Calibration Test Site at Railroad Valley, Nevada[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2016, 9(2): 816-826.

[11] 李新, 鄭小兵, 尹亞鵬. 場地自動(dòng)化定標(biāo)技術(shù)進(jìn)展[J]. 大氣與環(huán)境光學(xué)學(xué)報(bào), 2014, 9(1): 17-21. LI Xin, ZHENG Xiaobing, YIN Yapeng. Progress in Site Automation Calibration Technology[J]. Journal of Atmospheric and Environmental Optics, 2014, 9(1): 17-21. (in Chinese)

[12] 尹亞鵬, 李新, 鄭小兵, 等. 場地自動(dòng)化觀測輻射計(jì)的設(shè)計(jì)與實(shí)現(xiàn)[J]. 大氣與環(huán)境光學(xué)學(xué)報(bào), 2016, 11(1): 44-50. YIN Yapeng, LI Xin, ZHENG Xiaobing, et al. Design and Implementation of Site Automation Observation Radiometer[J]. Journal of Atmospheric and Environmental Optics, 2016, 11(1): 44-50. (in Chinese)

[13] LI X, YIN Y P, LIU E C, et al. Automated Test-site Radiometer for Vicarious Calibration[C]//SPIE Asia-Pacific Remote Sensing, 2014, Beijing, China. SPIE: 2014: 92640V.

[14] 邱剛剛, 李新, 韋瑋, 等. 基于場地自動(dòng)化觀測技術(shù)的遙感器在軌輻射定標(biāo)試驗(yàn)與分析[J]. 光學(xué)學(xué)報(bào), 2016(7): 1-9. QIU Ganggang, LI Xin, WEI Wei, et al. In-orbit Radiation Calibration Test and Analysis of Remote Sensing Based on Site Automation Observation Technology[J]. Acta Optica Sinica, 2016(7): 1-9. (in Chinese)

[15] 呂佳彥, 何明元, 陳林, 等. 基于敦煌輻射校正場的自動(dòng)化輻射定標(biāo)方法[J]. 光學(xué)學(xué)報(bào), 2017(8): 25-32. LYU Jiayan, HE Mingyuan, CHEN Lin, et al. Automated Radiometric Calibration Method Based on Dunhuang Radiometric Test Site[J]. Acta Optica Sinica, 2017(8): 25-32. (in Chinese)

[16] 張孟, 韋瑋, 張艷娜, 等. 基于場地自動(dòng)化觀測技術(shù)的SNPP VIIRS高頻次在軌輻射定標(biāo)[J]. 光子學(xué)報(bào), 2019, 48(4): 149-160. ZHANG Meng, WEI Wei, ZHANG Yanna, et al. High-frequency On-orbit Radiometric Calibration of SNPP VIIRS Based on In-site Automatied Observation Technology[J]. Acta Photonica Sinica, 2019, 48(4): 149-160. (in Chinese)

[17] 孫凌, 胡秀清, 郭茂華, 等. 風(fēng)云三號(hào)A星中分辨率光譜成像儀反射太陽波段的多場地定標(biāo)跟蹤[J]. 氣象科技進(jìn)展, 2013(4): 86-98. SUN Ling, HU Xiuqing, GUO Maohua, et al. Multisite Calibration Tracking for FY-3A MERSI Solar Bands[J]. Meteorological Science and Technology Progress, 2013(4): 86-98. (in Chinese)

[18] WANG L, HU X, CHEN L, et al. Consistent Calibration of VIRR Reflective Solar Channels Onboard FY-3A, FY-3B, and FY-3C Using a Multisite Calibration Method[J]. Remote Sensing, 2018, 10(9): 1336.

[19] COSNEFROY H, BRIOTTET X, LEROY M. Selection and Characterization of Saharan and Arabian Desert Sites for the Calibration of Optical Satellite Sensors[J]. Remote Sensing of Environment, 1996, 58(1): 101-114.

[20] ODONGO V O, HAMM N A S, MILTON E J. Spatio-temporal Assessment of Tuz G?lü, Turkey as a Potential Radiometric Vicarious Calibration Site[J]. Remote Sensing, 2014, 6(3): 2494-2513.

[21] 顧行發(fā). 航天光學(xué)遙感器輻射定標(biāo)原理與方法[M]. 北京: 科學(xué)出版社, 2013.GU Xingfa. Principle and Method of Radiation Calibration for Space Optical Remote Sensors[M]. Beijing: Science Press, 2013. (in Chinese)

[22] 龍亮, 王世濤, 周峰. 航天遙感器輻射定標(biāo)精度計(jì)算方法研究[J]. 航天返回與遙感, 2012, 33(6): 51-55. LONG Liang, WANG Shitao, ZHOU Feng. Research on Calculation Method of Radiation Calibration Accuracy of Space Remote Sensor[J]. Spacecraft Recovery & Remote Sensing, 2012, 33(6): 51-55. (in Chinese)

[23] 余譚其, 韋瑋, 張艷娜, 等. 敦煌輻射校正場春季BRDF特性分析[J]. 光子學(xué)報(bào), 2018, 47(6): 197-206. YU Tanqi, WEI Wei, ZHANG Yanna, et al. Analysis of Spring BRDF Characteristics of Dunhuang Radiometric Calibration Site[J]. Acta Photonica Sinica, 2018, 47(6): 197-206. (in Chinese)

Calibration Tracking of the Optical Payload Based on Automatic Methods

HE Lingli1,3HU Xiuqing2,3*WANG Ling2,3CHEN Lin2,3XU Na2,3ZHANG Peng2,3

（1 Chinese Academy of Meteorological Science, China Meteorological Administration, Beijing 100081, China）（2 Key Laboratory of Radiometric Calibration and Validation for Environmental Satellite, China Meteorological Administration, Beijing 100081, China）（3 National Satellite Meteorological Center, China Meteorological Administration, Beijing 100081, China）

In order to achieve continuous performance monitoring for the optical sensor in orbit, this paper uses the automatic observation data of the Dunhuang site to continuously calibrate and track the data of the Fengyun-3 satellite Visible and Infrared Radiometer (VIRR) for three years and calculate the calibration slope for each channel of VIRR. The calibration based on the automatic observation has the advantage in significantly increasing the calibration frequency over the traditional calibration methods. The accuracy of the automatic calibration method is verified by the observation data of American high-precision Aqua Moderate Resolution Imaging Spectroradiometer (MODIS), which shows that the method has a 3% calibration accuracy in the visible and near infrared (VNIR) and a 5% calibration accuracy in short wave infrared (SWIR). Comparing to the operational multi-site method, the relative deviation of mean value on the three-year scale is less than 2%. The tracking results of VIRR show that long-term calibration monitoring can be achieved based on automatic site observations. This method can be extended to other space optical payloads to truly realize the operation of automatic site calibration with multiple payloads.

Visible and Infrared Radiometer; reflective solar bands; automatic radiometric calibration; reflectance-based method; FY-3 satellite; space remote sensing

TP722.4

A

1009-8518(2020)06-0103-11

10.3969/j.issn.1009-8518.2020.06.010

2019-12-03

國家重點(diǎn)研發(fā)計(jì)劃（2018YFB0504900，2018YFB0504901）；國家自然科學(xué)基金（41871249，41471302）

何靈莉, 胡秀清, 王玲, 等. 基于自動(dòng)化方法的光學(xué)載荷定標(biāo)跟蹤[J]. 航天返回與遙感, 2020, 41(6): 103-113.

HE Lingli, HU Xiuqing, WANG Ling, et al. Calibration Tracking of the Optical Payload Based on Automatic Methods[J]. Spacecraft Recovery & Remote Sensing, 2020, 41(6): 103-113. (in Chinese)

何靈莉，女，1994年生，2020年獲中國氣象科學(xué)研究院大氣物理與大氣環(huán)境專業(yè)碩士學(xué)位，現(xiàn)在中國氣象科學(xué)研究院攻讀氣象博士學(xué)位。研究方向?yàn)楹Ｑ筝椛鋫鬏斃碚?。E-mail：hllsat@163.com。

（編輯：龐冰）