天繪一號衛(wèi)星相機(jī)的輻射交叉定標(biāo)研究

竇乃英，叢 琳，黃令勇

（1.中國天繪衛(wèi)星中心，北京 102102）

天繪一號（TH-1）衛(wèi)星開展長期在軌業(yè)務(wù)運(yùn)行，已為國民經(jīng)濟(jì)建設(shè)提供了準(zhǔn)確的空間信息源。由于衛(wèi)星在軌運(yùn)行期間外部環(huán)境將發(fā)生變化，所以必須對其進(jìn)行絕對輻射定標(biāo)。我國主要通過衛(wèi)星輻射定標(biāo)場對遙感衛(wèi)星開展場地定標(biāo)實(shí)驗，但運(yùn)用場地法定標(biāo)需要大量的地面實(shí)際測量數(shù)據(jù)，財力、物力、人力消耗較大，不能滿足TH-1衛(wèi)星業(yè)務(wù)化定標(biāo)的需要[1]。

目前，國內(nèi)外積極開展交叉定標(biāo)實(shí)驗，進(jìn)行多傳感器、多時相衛(wèi)星數(shù)據(jù)的交叉定標(biāo)，以滿足衛(wèi)星業(yè)務(wù)化定標(biāo)需求。國內(nèi)外專家[2-4]均利用交叉定標(biāo)方法對遙感衛(wèi)星進(jìn)行絕對輻射定標(biāo)。Landsat8是國際上業(yè)務(wù)穩(wěn)定運(yùn)行、數(shù)據(jù)廣泛應(yīng)用的高精度遙感衛(wèi)星，絕對輻射定標(biāo)精度為5%。利用Landsat8的觀測數(shù)據(jù)對TH-1衛(wèi)星相機(jī)進(jìn)行交叉定標(biāo)，建立穩(wěn)定可靠的TH-1衛(wèi)星與Landsat8之間的交叉定標(biāo)模型，能得到可靠的定標(biāo)系數(shù)，以滿足TH-1衛(wèi)星業(yè)務(wù)化定標(biāo)需求。

1 交叉輻射定標(biāo)原理

交叉定標(biāo)法是以定標(biāo)精度高的傳感器為參考標(biāo)準(zhǔn)，選擇相同或相近觀測的衛(wèi)星影像，建立參考傳感器入瞳輻亮度與目標(biāo)傳感器之間的關(guān)系[5]?；贚andsat8的TH-1衛(wèi)星交叉定標(biāo)方法主要分為3個步驟：①數(shù)據(jù)選取，選擇時間、空間和觀測幾何接近的數(shù)據(jù)；②分析Landsat 8傳感器和TH-1傳感器的光譜差異，確定TH-1傳感器的入瞳輻亮度值；③定標(biāo)計算，將經(jīng)過挑選和修正的樣本數(shù)據(jù)建立關(guān)系，獲取定標(biāo)系數(shù)。利用Landsat8標(biāo)定TH-1衛(wèi)星的流程如圖1所示。

圖1 TH-1衛(wèi)星和Landsat8交叉定標(biāo)流程圖

利用Landsat8的OLI 相機(jī)對TH-1衛(wèi)星交叉定標(biāo)，TH-1衛(wèi)星的相機(jī)輻射特性采用線性設(shè)計，即具有式（1）的關(guān)系，則可通過計算TH-1衛(wèi)星相機(jī)波段i的表觀輻亮度和對應(yīng)波段的灰度值DN，得到相機(jī)波段i的絕對輻射定標(biāo)系數(shù)G和b，從而完成相機(jī)的絕對輻射定標(biāo)。

式中，G為定標(biāo)增益；b為定標(biāo)截距，單位為W.m-2.s-1.um-1；DNi為TH-1衛(wèi)星相機(jī)波段i的灰度值。

由于Landsat8的OLI相機(jī)與TH-1衛(wèi)星相機(jī)的光譜響應(yīng)函數(shù)存在差異，因此需建立兩個相機(jī)光譜輻亮度之間的關(guān)系?？紤]到兩個相機(jī)對大氣條件、地物響應(yīng)的差異和不同觀測角度等因素，得到兩個相機(jī)入瞳輻亮度的關(guān)系。TH-1衛(wèi)星多光譜相機(jī)表觀輻亮度的絕對輻射定標(biāo)公式為：

TH-1衛(wèi)星多光譜相機(jī)表觀反射率的絕對輻射定標(biāo)公式為：

表觀反射率與表觀輻亮度的關(guān)系為：

即

由式（2）、式（5）可得：

Landsat8表觀反射率的絕對輻射定標(biāo)關(guān)系為：

式中，為光譜匹配因子，本文利用6S輻射傳輸模型計算得到；= ρ cosθ 、= ρ c osθ 分 TH i TH Li L別為TH-1衛(wèi)星和Landsat8多光譜相機(jī)過頂時通道i的表觀反射率；θTH、θL分別為TH-1衛(wèi)星和Landsat8相機(jī)過頂時的太陽天頂角；ECi為第i通道的日地平均距離處太陽輻照度；d2為太陽輻照度的日地距離修正因子；CLi和PLi為Landsat8的表觀反射率定標(biāo)系數(shù)，可在Landsat8數(shù)據(jù)頭文件中獲?。籇NLi、DNTHi分別為Landsat8和TH-1衛(wèi)星對應(yīng)波段處的灰度值；GTHi、bTHi分別為TH-1衛(wèi)星絕對輻射定標(biāo)的增益和偏置。

2 定標(biāo)實(shí)驗與結(jié)果

2.1 定標(biāo)靶區(qū)的選擇

在提取影像中定標(biāo)點(diǎn)數(shù)字計數(shù)值的過程中，像元定位會存在一定偏差，因此交叉輻射定標(biāo)的定標(biāo)靶區(qū)一般選用大部分均勻場地，地勢平坦，地物變化和BRDF的影響較小，具有較高的穩(wěn)定性和均勻性[6]。本文選用的定標(biāo)靶區(qū)為云南均勻場，場點(diǎn)屬于高原山地氣候兼亞熱帶季風(fēng)氣候，屬海洋型氣溶膠模式，地物的光譜特性如圖2所示。云南均勻場不僅能滿足交叉定標(biāo)對場地光譜特性等方面的要求，而且每年國內(nèi)相關(guān)部門在該地均有實(shí)地同步數(shù)據(jù)測量，因此可以有大量歷史數(shù)據(jù)進(jìn)行參考和使用。

圖2 交叉定標(biāo)靶區(qū)光譜反射率

2.2 匹配數(shù)據(jù)的獲取

TH-1衛(wèi)星和Landsat8均為太陽同步軌道，綜合考慮地面場地條件和衛(wèi)星觀測幾何等因素，以TH-1衛(wèi)星的實(shí)際過場時間為標(biāo)準(zhǔn)，在USGS官網(wǎng)上下載最相鄰時間的Landsat8過場影像；并下載同時間的MODIS氣溶膠和反射率產(chǎn)品。目前，MODIS數(shù)據(jù)產(chǎn)品已經(jīng)過多次對比和驗證，確保了其數(shù)據(jù)的準(zhǔn)確性和適用性，因此本文以它提供的氣溶膠和反射率數(shù)據(jù)作為已知可信數(shù)據(jù)使用[7]。TH-1衛(wèi)星與Landsat8的過境時間與觀測幾何如表1所示。

表1 TH-1衛(wèi)星與Landsat8的過境時間與觀測幾何

2.3 光譜匹配因子的計算

光譜匹配因子計算是交叉定標(biāo)的關(guān)鍵環(huán)節(jié)。TH-1衛(wèi)星與Landsat8相匹配的4個多光譜波段的光譜范圍如表2所示。

表2 Landsat8與TH-1衛(wèi)星相應(yīng)的光譜范圍

為了比較不同相機(jī)光譜響應(yīng)對表觀輻亮度的影響，利用6S輻射傳輸模型計算對應(yīng)各通道的表觀反射率，其中定標(biāo)點(diǎn)的海拔確定為2.476 km，地表反射率和550 nm處的氣溶膠光學(xué)厚度采用定標(biāo)場的實(shí)時MODIS數(shù)據(jù)。

交叉定標(biāo)包括單點(diǎn)法交叉定標(biāo)、雙點(diǎn)法交叉定標(biāo)和多點(diǎn)法交叉定標(biāo)3種方法，若傳感器各通道的偏移量為零，則只需一個點(diǎn)即可計算傳感器某通道的輻亮度定標(biāo)系數(shù)增益；若傳感器各通道的偏移量不為零，那么至少需要兩個點(diǎn)才能計算輻亮度定標(biāo)系數(shù)增益。由于TH-1衛(wèi)星交叉定標(biāo)系統(tǒng)的研究尚處于起步階段，自動化程度極低，為了使結(jié)果更具科學(xué)性和普適性，兼顧工作量能夠負(fù)擔(dān)的情況下，本文在云南均勻場選取了10對廣泛離散分布的同名點(diǎn)，光譜匹配因子計算結(jié)果如表3所示。

通過光譜匹配因子的計算，即可實(shí)現(xiàn)利用歸一化的表觀反射率代替表觀輻亮度來進(jìn)行交叉定標(biāo)，從而消除了TH-1衛(wèi)星和Landsat8兩個傳感器不同成像時間對太陽天頂角余弦值的影響以及不同的光譜響應(yīng)帶來的表觀輻亮度差異[2]。

表3 TH-1衛(wèi)星與Landsat8相應(yīng)通道光譜匹配因子值

2.4 定標(biāo)系數(shù)的計算

根據(jù)影像成像日期計算得到日地平均距離修正因子為0.966 8；日地平均距離處TH-1衛(wèi)星CCD各通道等效太陽輻照度ECi的計算結(jié)果分別為1 913.75、1 840.08、1 671.59和1 412.09，具體計算過程在此不贅述。本文在定標(biāo)區(qū)域采集10個樣本點(diǎn)數(shù)據(jù)，根據(jù)計算公式，最終得到TH-1衛(wèi)星的等效入瞳輻亮度以及對應(yīng)的DN值如表4所示。

表4 10個采樣點(diǎn)各通道的入瞳輻亮度以及對應(yīng)的DN值

根據(jù)入瞳輻亮度以及與對應(yīng)傳感器計數(shù)值之間的線性關(guān)系，利用線性回歸分析，得到絕對定標(biāo)系數(shù)結(jié)果如表5所示。

表5 各波段定標(biāo)系數(shù)結(jié)果

2015年9月13日國內(nèi)相關(guān)部門基于固定輻射靶標(biāo)對TH-1衛(wèi)星的多光譜載荷進(jìn)行了在軌絕對輻射定標(biāo)，定標(biāo)結(jié)果如表6所示。

表6 TH-1衛(wèi)星場地定標(biāo)結(jié)果

由表5、6可知，交叉定標(biāo)結(jié)果和野外實(shí)驗檢測得到的在軌絕對輻射定標(biāo)系數(shù)K和B都相差不大，增益的相對誤差分別為5.5%、7.1%、7%和4.2%。筆者認(rèn)為出現(xiàn)誤差可能的原因包括：①采樣點(diǎn)數(shù)非最優(yōu)；②Landsat8標(biāo)準(zhǔn)傳感器自身定標(biāo)系數(shù)誤差；③像元匹配誤差；④光譜匹配因子計算誤差，其中替代模型使得參數(shù)與真實(shí)值之間存在偏差；⑤MODIS數(shù)據(jù)產(chǎn)品誤差；⑥傳感器之間存在差異，如獲取時間、拍攝角度等。因此，本文在進(jìn)行影像對選擇時盡量保證成像時間和觀測幾何一致，但由于不同系列衛(wèi)星其相應(yīng)條件不可能完全一致，因此其相關(guān)性上必然存在一定誤差。

3 結(jié) 語

基于定標(biāo)靶標(biāo)場的均勻場地交叉定標(biāo)，綜合考慮了光譜響應(yīng)差異、觀測時間等因素，獲得的交叉定標(biāo)系數(shù)精度可靠；因此可認(rèn)為以定標(biāo)精度較高的標(biāo)準(zhǔn)傳感器為標(biāo)準(zhǔn)，利用線性回歸方法獲取的相機(jī)交叉定標(biāo)精度與TH-1衛(wèi)星在軌場地定標(biāo)精度有很好的一致性，達(dá)到了業(yè)務(wù)化運(yùn)行的定標(biāo)精度要求。今后的研究中，將在不確定度的分析、定標(biāo)點(diǎn)數(shù)最優(yōu)選擇上下功夫，以期進(jìn)一步提高交叉定標(biāo)的普適性。