環境衛星CCD傳感器場地輻射定標與交叉定標的比較*

余曉磊，巫兆聰

(武漢大學遙感信息工程學院，武漢 430079)

環境衛星(HJ－1A/B)是我國新近發射的減災防災衛星星座，每顆衛星并排搭載兩臺CCD傳感器，每臺傳感器有三個可見光譜段和一個近紅外譜段(具體參數參見表1)。為深入其定量化研究，需要對傳感器定期定標，以提供準確的定標系數。輻射定標就是要建立傳感器的每個探測元件輸出的數字量化信號(DN)與它所對應視場中輸出輻亮度值之間的定量關系［1－2］。本文利用2009年8月份敦煌衛星輻射校正場同步測量的地物光譜反射率數據和大氣光學特征參數，計算HJ－1A衛星CCD2傳感器4個波段的定標系數。針對單一地物定標法，以及環境衛星缺乏傳感器冷空間計數值的問題，使用同期過境的LANDSAT－5衛星TM影像進行交叉定標。結果表明，場地絕對輻射定標和交叉定標的平均差異在10%以內，具有良好的一致性。同時環境衛星CCD傳感器自發射以來，各波段有不同程度的退化，其中藍光通道的退化最為嚴重。

1 定標原理與算法

1.1 反射率法定標原理

基于反射率法的場地絕對輻射定標是在衛星傳感器過境時同步測量地面目標反射率因子和大氣光學特征參數(如大氣氣溶膠光學厚度、大氣柱狀水汽含量等)，然后利用輻射傳輸模型(其一般輸入數據有分子散射率和氣溶膠光學厚度、氣溶膠的復折射指數和Mie散射指數等)計算出傳感器入瞳處的反射率值［1，3－4］。該方法相比于輻亮度法、輻照度法，精度較高，測量的參數相對較少，但是必須保證所需參量都是同步觀測時獲取［2］，其原理如下:

其中ai，bi分別表示該通道的增益和偏置。進一步傳感器單通道入瞳處的大氣頂層表觀反射率可以表示為:

式中d表示衛星獲取影像時真實的日地距離和日地平均距離之比;θs表示衛星獲取影像時的太陽天頂角;E0，i表示傳感器單通道大氣層外太陽輻照度，又稱為波段太陽常數。根據大氣輻射傳輸理論，傳感器在大氣上層測量的目標反射率可以表示為:

其中Tg表示大氣透射率，ρA+R表示大氣分子自身含氣溶膠層在內的反射率，Tdown，Tup表示有太陽到地表以及地表到傳感器的大氣透射率，S表示大氣半球反照率表示朗伯體地物目標的反射率因子經過第i通道的傳感器光譜響應函數積分后得到的地表反射率。如前所述通過大氣輻射傳輸模型(常見的例如6S、LOWTRAN、MODTRAN)可以通過地表反射率、大氣特征參數以及觀測幾何條件計算得到傳感器入瞳處的大氣頂層的表觀反射率［5－6］。

1.2 交叉定標原理

根據公式1、2、3可知，為了計算傳感器的定標系數，即單通道的偏置和增益ai、bi至少需要兩組計數值及對應的地物反射率。對于我國先后發射的資源衛星和環境減災衛星的可見光、近紅外通道，其CCD線陣掃描成像式傳感器不具備觀測冷空間的能力，也沒有星上定標裝置，所以在進行單一地物法場地絕對輻射定標時，通常認為傳感器單通道的偏置值bi已知，且為發射前實驗室檢校得到的結果。但事實上，由于傳感器所處的外太空環境復雜以及長期工作后設備老化，其偏置往往會發生比較大的變化，美國NOVAA系列衛星AVHRR傳感器的定標結果表明，其各通道的偏置隨時間存在較大的變化［6－7］。為此本研究使用 LANDSAT－5 衛星 TM 影像對環境衛星CCD影像進行交叉定標。

交叉定標是使用較高精度的衛星通道定標較低精度的衛星通道的一種方法，已經廣泛的用于AVHRR，MODIS等傳感器的輻射定標［8］，其基本原理如下:

其中k表示兩種不同傳感器之間對應通道間的光譜匹配系數［9］:

式中:ρs(λ)為典型地物光譜反射率，f(λ)與f'(λ)分別為兩類不同傳感器分波段的光譜響應函數，λ1～λ2和λ'1～λ'2分別為兩傳感器分波段光譜的起止波長。

對于不同類型的地物而言，其光譜匹配因子一般不同，選取兩類或兩類上的典型地物，代入式(4)中，便可建立方程組求解傳感器的定標參數［10］。根據表1、2所示環境衛星CCD傳感器與LANDSAT－5衛星TM傳感器的光譜參數對比可知兩類傳感器的1～4波段設置類似，且二者的空間分辨率接近(30 m)，量化等級相同(8 bit)。而TM傳感器經過嚴格的發射前定標，長期在軌監測和定期場地絕對定標，有很高的輻射定標精度，所以用TM影像對環境衛星CCD影像進行交叉定標是最佳的選擇。

表1 HJ衛星CCD傳感器參數

2 定標數據的采集與處理

敦煌衛星輻射定標場是我國衛星輻射定標的主要場地，自1995年建立以來承擔了我國多種類型衛星的場地輻射定標任務。它位于東經94.4°，北緯40.08°的戈壁上，整個場區平坦均勻，基本無植被覆蓋，具有穩定的地表光學特性，年變化很小，海拔高度1 300 m左右，大氣干潔，受云霧等影響較?。?1］。

2.1 地表反射率測量

本次定標所需的地面實測數據，為中國資源衛星應用中心2009年8月26日于甘肅敦煌衛星定標場實地測量所得，戈壁校正場地表反射率的測量采用550 m×550 m的測試區，每隔50 m測量一個點，在衛星過境一個半小時內利用ASD FR型便攜式光譜儀共計進行144個點的測量，每個測量點對地物目標測量5次，同時對參考白板測量2次。

按照以下公式處理得到定標場地物目標的反射率:

式中ρs(λ)表示參考白板的反射率，通常為實驗室檢校的結果，V(λ)表示地物目標測量值，Vs(λ)表示參考白板測量值，ρ(λ)表示地物目標的反射率。計算測量點的平均反射率如圖1所示。

圖1 2009年8月26日敦煌定標場反射率

2.2 大氣參數測量

在衛星過境時刻，使用CE318太陽光度計同步采集了場地上空的氣溶膠光學厚度數據。太陽光度計采用蘭立(LangLey)法定標［12］。氣溶膠光學厚度的測量按照大氣渾濁度觀測標準進行，其結果如圖2所示。氣象數據來自敦煌地區氣象站提供的測量結果，包括環境溫濕度、大氣壓、露點溫度等。

圖2 2009年8月26日大氣氣溶膠特征

由于本文使用的6S輻射傳輸模型需要550 nm處的氣溶膠光學厚度作為模型的輸入參數，而CE318的標準濾光鏡沒有對應的波段，所以測量的結果需要按照Angstrom公式進行換算:

式中τ為氣溶膠光學厚度，λ為波長，β為渾濁度系數，α為Angstrom波長指數，利用測量的任意兩個波長的氣溶膠光學厚度，即可計算出α與β，繼而求出550 nm處的大氣氣溶膠光學厚度。

2.3 定標影像數據處理

場地絕對定標采用的影像數據為2009年8月26日過境的HJ－1A衛星CCD2影像，以保證與地面實測嚴格同步。交叉定標所采用的影像分別為2009年6月21日過境的HJ－1A衛星CCD2影像和6月22日過境的LANDSAT－5衛星TM影像。其中8月26日的影像利用1∶50 000的數字地形圖進行幾何校正，6月21日的環境衛星影像和6月22日的LANDSAT－5衛星TM影像，進行了相對幾何校正，幾何校正的誤差控制在1個像素范圍內，以保證定位精度。

本文選取敦煌定標場內戈壁灘和敦煌地區的農作物這兩類具有明顯差別的地物進行交叉定標。地物的光譜反射率來源于中國典型地物波譜數據集中對敦煌地區長期測量的平均值。計算兩類傳感器對應通道之間的光譜匹配系數k所需的光譜響應函數和典型地物光譜反射率如圖3所示。計算所得的匹配系數如表2所列。

圖3 HJ－1A－CCD2 和 Landsat－TM5 光譜響應函數

表2 HJ－1A－CCD2與Landsat－TM5匹配通道及匹配因子

3 輻射定標系數計算與結果分析

場地絕對輻射定標和交叉輻射定標的工作流程分別如圖4和圖5所示。本次絕對輻射定標的結果和08年中國資源衛星應用中心公布的發射前定標系數，以及相對輻射定標的結果如表3所示。

圖4 絕對輻射定標工作流程

圖5 交叉輻射定標工作流程

表3 HJ－1A－CCD2傳感器定標結果對比

3.1 場地定標與發射前定標系數的比較

對比場地絕對輻射定標結果和發射前傳感器的定標系數可知，自發射以來，HJ－1A衛星CCD2傳感器的各波段都存在不同程度的退化，其中以第一波段(藍光通道)的退化最為嚴重，達6.7%。而第四波段(近紅外波段)的變化相對較小，為1.8%，第二(綠光)和第三(紅光)波段的絕對定標系數變化居中，分別為3.4%和4.8%。表明衛星發射升空和在軌運行以來，傳感器的性能相對實驗室檢定而言，發生了較大的變化，這主要是由于傳感器隨運載火箭發射升空時的復雜狀態和傳感器長期在軌工作造成的性能退化。

3.2 交叉輻射定標的驗證

為驗證交叉輻射定標的結果，通過選取2009年6月21日HJ－1A衛星CCD2影像上的四類典型地物目標(水體、城區、植被、干河床)，分別用絕對定標系數和交叉定標系數按照式1求出地物在對應波段的輻亮度，并計算其相對偏差，結果如表4所示。比較發現，交叉定標和場地絕對定標結果，對于城區和植被，其符合性較好，平均偏差在5%以內;對于水體和干河床，其最大偏差分別達到了6.62%和8.38%。這主要是由于，對城區而言，其光譜反射率曲線的特征在傳感器的四個通道上都變化都比較平緩，而植被的光譜反射率曲線在藍光波段較低而在近紅外和綠光波段上相對較高;對于干凈的內陸水體，其光譜反射率曲線在藍光波段上的瑞麗散射占總信號的80%［13］;而干河床在所有波段上的反射率都相對較大。所以才會帶來其絕對定標和交叉定標的部分差異。但就整體而言，交叉定標和絕對定標的結果偏差在10%以內。這表明，交叉輻射定標可以作為替代定標的手段，進行傳感器的長期在軌狀態監測。

表4 HJ－1A－CCD2傳感器交叉定標結果驗證

4 結束語

通過對環境衛星CCD傳感器的絕對輻射定標和交叉定標的比較，分析傳感器性能的變化，驗證了交叉定標與絕對定標的一致性。對于監測在軌運行的星載傳感器而言，還需要做進一步的真實性檢驗和長期的定標實驗，以了解其運行狀況，促進國產資源衛星影像的定量化應用。

［1］鞏慧，田國良，余濤，等.CBERS02B衛星CCD相機在軌輻射定標與真實性檢驗［J］.遙感學報，2010，14(1):7－12.

［2］童慶禧，張兵，鄭蘭芬，等.高光譜遙感［M］.北京:高等教育出版社出版社，2006，32－35，68－74.

［3］傅俏燕，閔祥軍，李杏朝，等.敦煌場地CBERS－02CCD傳感器絕對輻射定標研究［J］.遙感學報，2010，10(4):433－439.

［4］李小燕，顧行發，閔祥軍，等.利用MODIS對CBERS－02衛星CCD相機進行輻射交叉定標［J］.中國科學 E輯，2005，35(增刊):41－58.

［5］Kastner C J，Slater P N.In-Flight Radiometric Calibration of Advance Remote Systems［J］.SPIE，1982，356:158－165.

［6］顧明澧.星載遙感器在飛行時的絕對輻射定標方法［J］.航天返回與遙感，2000，21(1):16－23.

［7］Black S E，Helder D L，Schiller S J.Irradiance-Based Cross-Calibration of Landsat-5 and Landsat-7 Thematic Mapper Sensors［J］.IJRS，2003，24(2):287－304.

［8］陳正超，劉翔，李俊生，等.北京一號衛星多光譜遙感器交叉定標［J］.宇航學報，2008，29(2):638－645.

［9］唐軍武，顧行發，牛生麗，等.基于水體目標的對CBERS－02衛星CCD相機與MODIS的交叉輻射定標［J］.中國科學E輯，2005，35(增刊):59－69.

［10］彭光雄，何宇華，李京，等.中巴資源02星CCD圖像交叉定標與大氣校正研究［J］.紅外與毫米波學報，2007，26(1):29－34.

［11］中國遙感衛星輻射校正場專題論證組.中國遙感衛星輻射校正場技術方案［R］.北京:中國資源衛星應用中心，1994，7－8.

［12］中華人名共和國氣象行業標準:QXT69－2007，大氣渾濁度觀測－太陽光度計方法［S］.

［13］孫凌，張杰，郭茂華.針對HY－1A CCD數據處理的瑞利查找表［J］.遙感學報，2006，10(3):306－311.