肖艷芳，朱 琳，趙文吉

（首都師范大學(xué) 三維信息獲取與應(yīng)用教育部重點實驗室，北京100048）

HJ－1A星輻射定標及近海水體大氣校正的精度驗證＊

肖艷芳，朱 琳，趙文吉

（首都師范大學(xué) 三維信息獲取與應(yīng)用教育部重點實驗室，北京100048）

以應(yīng)用HJ－1A星數(shù)據(jù)進行近海水體定量反演為目的，對HJ－1A星CCD的星上定標系數(shù)的精度進行驗證，并基于Landsat－5TM數(shù)據(jù)對其進行交叉定標。在此基礎(chǔ)上反演近海水體的反射率，將反演結(jié)果與實測水體反射率進行比較發(fā)現(xiàn)：對HJ－1A星CCD數(shù)據(jù)進行交叉定標能在一定程度上能提高近海水體反射率反演的精度。結(jié)果表明，HJ－1A星CCD影像的星上定標系數(shù)不夠精確，在利用該數(shù)據(jù)進行近海水體的定量反演研究時，需要利用交叉定標對星上定標參數(shù)進行修正。

交叉定標；星上定標；反射率反演；近海水體；HJ－1A；Landsat－5TM

HJ－1A／B衛(wèi)星是HJ－1星座的2顆光學(xué)小衛(wèi)星，2008－09－06在太原衛(wèi)星發(fā)射中心成功發(fā)射。2顆衛(wèi)星上都搭載2臺寬覆蓋多光譜可見光CCD相機。同時，A星上搭載一臺超廣譜成像儀，B星上搭載一臺紅外相機。CCD相機共設(shè)置4個波段，地面像元分辨率為30m，單臺相機幅寬360km，48h即可對全球覆蓋一次。作為我國首顆環(huán)境與災(zāi)害監(jiān)測國際合作衛(wèi)星，HJ－1星座在環(huán)境與災(zāi)害事件的應(yīng)急觀測與快速評估、資源調(diào)查與管理、林業(yè)與農(nóng)業(yè)生產(chǎn)等領(lǐng)域具有很好的應(yīng)用前景，對我國及周邊國家的遙感應(yīng)用起到極大的作用。

衛(wèi)星發(fā)射后，一些研究者對衛(wèi)星的應(yīng)用潛力進行了分析［1－2］，目前環(huán)境一號小衛(wèi)星的數(shù)據(jù)已經(jīng)應(yīng)用在生態(tài)環(huán)境監(jiān)測［3］、水質(zhì)監(jiān)測［4－5］、大氣監(jiān)測［6］和農(nóng)作物病蟲害監(jiān)測［7］方面，對于數(shù)據(jù)星上定標質(zhì)量的評價主要集中熱紅外數(shù)據(jù)［8－9］和高光譜數(shù)據(jù)［10］上，關(guān)于CCD相機星上定標精度的驗證鮮有研究。本研究利用Landsat－5的TM數(shù)據(jù)對HJ－1A星CCD相機的星上定標結(jié)果進行驗證，并對其進行基于TM的交叉定標計算；利用基于MODIS數(shù)據(jù)的大氣校正方法，反演近海水體的反射率?；趯崪y的水體反射率，分析交叉定標對HJ－1A星的CCD影像反演近海水體反射率的影響。

1 輻射定標

輻射定標主要是指將星上觀測到的遙感計數(shù)值轉(zhuǎn)換為物理量絕對值的過程［11］。在大氣層外，太陽的輻照度可以認為是一個常數(shù)，因此可以選擇太陽作為基準光源，通過太陽定標系統(tǒng)對星載成像光譜儀器進行輻射定標，即星上定標，如TM、MODIS等數(shù)據(jù)。在這類數(shù)據(jù)的頭文件中都會提供相應(yīng)的定標系數(shù)，HJ－1A的CCD也提供了定標系數(shù)，但定標系數(shù)的精確程度是否能夠滿足近海水體定量反演的要求還有待商討。對沒有提供定標系數(shù)或定標系數(shù)不精確的衛(wèi)星影像，需要用定標精度較高的衛(wèi)星對其進行交叉定標［12］。Landsat－5衛(wèi)星經(jīng)過精確的標定和嚴格的在軌測試，運行期間為了保持數(shù)據(jù)的一致性和數(shù)據(jù)精度而不斷更新定標算法。同時，HJ－1ACCD的波段設(shè)置與Landsat－5TM前4個波段的設(shè)置相似（表1）。

表1 HJ－1A星CCD和Landsat－5TM對應(yīng)波段帶寬（μm）Table 1 Bandwidth of CCD and Landsat－5TM（μm）

1.1 星上定標精度驗證

選取2009－06－25北京延慶盆地的HJ－1A星CCD數(shù)據(jù)和同日該地區(qū)的Landsat－5TM數(shù)據(jù)進行定標比對及交叉定標研究。兩景影像資料過境時間基本一致，且天空晴朗無云，完全可以用于進行對比分析。

利用式（1）將HJ－1A星CCD圖像的計數(shù)值轉(zhuǎn)換為輻亮度值：

式中，L為某個波段的光譜輻亮度；DN為圖像計數(shù)值；g和L0為輻射定標參數(shù)，可在圖像數(shù)據(jù)的頭文件中查找，見表2。

表2 HJ－1A星CCD波段1～4的星上定標系數(shù)Table 2 Radiance calibration coefficients of HJ－1Aband 1～4

本次研究使用的Landsat－5TM數(shù)據(jù)為1級產(chǎn)品，1級產(chǎn)品的DN值是由輻射亮度經(jīng)線性變換得到的，輻亮度的計算公式如下

式中，DN代表1級產(chǎn)品的像元值；Grescale和Brescale分別為變換的增益和偏置，見表3。

表3 Landsat－5TM波段1～4的輻射定標系數(shù)Table 3 Radiance calibration coefficients of Landsat－5TM band 1～4

以Landsat－5TM為基準，對兩景影像進行精確的幾何糾正，精度控制在0.5個像元以內(nèi)，確保2種數(shù)據(jù)的目標基本一致。然后在兩景影像對應(yīng)區(qū)域選取建筑物、裸地、水體和植被四類地物，每類地物分別采集10個樣點，樣點各波段的平均值作為此類地物在相應(yīng)波段的輻亮度值。由表4可知，對HJ－1A的CCD影像進行星上定標得到的輻亮度與Landsat－5TM的輻亮度存在較大差異，這說明HJ－1A的CCD影像所提供的星上定標系數(shù)達不到進行一系列水體定量反演的精度要求。

表4 HJ－1A星CCD星上定標與Landsat－5TM輻射定標結(jié)果比較Table 4 Comparison between the satellite Radiance calibration of CCD mounted on Satellite HJ－1A and the Radiance calibration of Landsat－5TM

1.2 交叉定標及精度分析

分別在兩景影像對應(yīng)位置選擇高、低計數(shù)值均勻的兩塊區(qū)域（每塊區(qū)域多于50個像元，均方根差小于3）［13－14］，應(yīng)用兩點法，建立Landsat－5TM 和 HJ－1A的CCD各對應(yīng)波段DN值之間的一元線性回歸方程如下

利用式（1），（3），（4），（5），（6）可計算得到 HJ－1ACCD數(shù)據(jù)1～4波段的交叉輻射定標系數(shù)，見表5。

表5 HJ－1A星CCD2數(shù)據(jù)1～4波段交叉定標系數(shù)Table 5 Coefficients of cross radiance calibration for HJ－1ACCD 1～4bands

2 基于MODIS影像的近岸水體反射率反演

2.1 模型計算

衛(wèi)星接收的總輻射亮度可表示為［15］

式中，Lr為大氣分子Rayleigh散射；La為大氣氣溶膠散射；Lw為離水輻亮度；Lf為海面白帽信號；Lg為海面對太陽光的直射反射；Lb為水底的反射；t為大氣漫射透過率；T為大氣直射透過率；λ為波段。

在不考慮來自水底反射的影響，同時避開太陽耀斑和白帽的條件下，衛(wèi)星接收到的總輻亮度可簡化為

二類水體大氣校正的目的就是要從衛(wèi)星接收的總輻亮度中去除大氣散射Lr，氣溶膠散射La，得到離水輻亮度Lw。

2.1.1 大氣散射Lr

根據(jù) Gordon單次散射，可推導(dǎo)出 Rayleigh散射［16－17］：

式中，cos（α±）＝±cosθcosθ0＋sinθsinθ0cos（φ－φ0），θ，φ分別為對應(yīng)像元到遙感器之間連線所構(gòu)成的向量的天頂角和方位角；θ0，φ0分別為對應(yīng)像元與太陽之間連線所構(gòu)成的向量的天頂角和方位角；α為兩個向量之間的夾角。在本研究中，以影像頭文件中衛(wèi)星的天頂角和方位角近似代替θ和φ，以太陽的天頂角和方位角近似代替θ0和φ0；ω0為瑞利單次散射反照率，在避開水汽和氧氣吸收的波段可近似為1［18］。

Pr（α）為瑞利散射相函數(shù)，公式如下

式中，τr為大氣瑞利光學(xué)厚度：τr（λ）＝（P／P0）τr0（λ）；P0為標準大氣壓；P為實際大氣壓［19］。

在本研究中，所使用的瑞利光學(xué)厚度應(yīng)為“波段等效光學(xué)厚度”［18］：

式中，Si（λ）為HJ－1A星CCD數(shù)據(jù)的波段響應(yīng)函數(shù)；λ1，λ2分別為波譜響應(yīng)函數(shù)中對應(yīng)波段的最大值和最小值；為經(jīng)過臭氧吸收修訂的大氣層外太陽輻照度［20］。

式中，＜F0＞為平均日地距離處的大氣層外太陽輻照度；D為從1月1日算起的相應(yīng)日期的公歷天數(shù)；τoz（λ）為臭氧層光學(xué)厚度，可通過NASA網(wǎng)站的實時大氣臭氧含量計算得到［18］。

ρ為氣水表面對天空光的反射率，由于清潔水體對近紅外波段幾乎全部吸收，這些波段的水體信號就是天空光的反射信號，因此可以通過式（14）確定氣水表面的反射率［21］：

本研究于衛(wèi)星過境當天在曹妃甸附近海域進行了水體光譜的測量，選取其中最為清潔的水體計算得到ρ為0.029；Lsw為光譜儀測量得到的水體輻亮度；Lsky為光譜儀測量得到的天空光輻亮度。

2.1.2 氣溶膠散射La

選取研究區(qū)內(nèi)遠離海岸線的清潔水體像元，其近紅外波段的離水輻亮度近似為0，即Lw（λ4）～0，忽略太陽直射反射和白帽的影響，第4波段的氣溶膠散射為

其他波段氣溶膠散射的計算根據(jù)Wang和Gordon修改的Angstrom定律［22］：

式中，c為常數(shù)，可利用同一塊清潔水體的HJ－1A星CCD圖像近紅外波段和對應(yīng)MODIS的近紅外波段數(shù)據(jù)得到［23－24］。

2.1.3 大氣漫射透過率t

式中，τr（λ）為瑞利散射光學(xué)厚度；τoz（λ）為臭氧光學(xué)厚度；ta（θ，λ）為氣溶膠散射透過率，由于其與氣溶膠光學(xué)厚度關(guān)系不大，通常都將其作為一個常數(shù)，取值為1［22］。

遙感反射率Rrs可通過式（18）計算：

Lwn（λ）為歸一化離水輻亮度：

2.2 數(shù)據(jù)處理

本研究以位于渤海灣西北的曹妃甸附近海域為研究區(qū)，區(qū)域示意圖及實測點位置如圖1所示。本研究的第1部分對2009－06－25的HJ－1A星CCD影像分別進行了基于星上定標系數(shù)的輻射定標和基于Landsat－5TM的交叉輻射定標，對上述2種輻射定標的結(jié)果，分別利用第2部分的基于MODIS的近海水體大氣校正模型反演水體反射率。

圖1 研究區(qū)及采樣點示意圖Fig.1 A sketch map of study area and sampling locations

水體光譜的測量與分析參照唐軍武的“水面以上測量法”［21］。對于未嚴格定標的光譜儀，遙感反射率可通過式（20）計算：

式中，Ssw（λ），Ssky（λ），Sp（λ）分別為光譜儀測得的水體、天空光以及白板的信號值；rs為標準板的反射率；ρ為水氣界面反射率，取值0.029，實測水體樣點的反射率曲線見圖2。

由于光譜響應(yīng)函數(shù)的不同，對于相同的入射光，不同傳感器在同一波長位置處接收的輻亮度會有較大差異。因此，為了便于比較，需要計算實測水體的等效波段反射率（圖3）：

式中，Si（λ）為HJ－1A星CCD的波段響應(yīng)函數(shù)；λ1，λ2分別為響應(yīng)函數(shù)對應(yīng)波段的最大值和最小值。

2.3 大氣校正結(jié)果分析

為了驗證交叉輻射定標對HJ－1A星CCD影像近海水體反射率反演的影響，下面結(jié)合實測水體反射率對大氣校正的結(jié)果進行比較和驗證。圖4和圖5分別與圖3相比顯示，交叉定標結(jié)果經(jīng)大氣校正后得到的水體反射率的波段1和波段4的反射率明顯降低，比星上定標結(jié)果經(jīng)大氣校正得到的水體反射率更接近于實測的水體反射率；波段2和波段3的差異不明顯。從水體反射率曲線上看，交叉定標得到的反射率曲線更接近實測的反射率曲線。表6給出了每個樣點的星上定標反射率和交叉定標反射率與實測反射率的相對誤差，其中，對HJ－1A星CCD影像進行交叉定標后，反演得到的波段1的反射率比進行星上定標得到的結(jié)果相對誤差降低了10%，波段4的反射率相對誤差則降低了120%左右，波段2和波段3沒有較大變化。可見，HJ－1A星CCD影像的星上定標參數(shù)對于反演近海水體反射率是相對不穩(wěn)定的，在進行近海水體一系列定量反演時，對其進行交叉輻射定標是必要的。

表6 HJ－1A星CCD星上定標反射率、交叉定標反射率與實測反射率的結(jié)果比較Table 6 Comparison among the reflectivities resulting from satellite radiance calibration of CCD，cross radiance calibration and in－situ measurements

3 結(jié) 語

本研究基于Landsat－5TM數(shù)據(jù)對HJ－1A星CCD的定標參數(shù)進行驗證，認為在利用HJ－1A星CCD數(shù)據(jù)進行高精度水體反演時直接使用星上定標系數(shù)會存在較大誤差，需要對星上定標參數(shù)進行修正，基于Landsat－5TM數(shù)據(jù)對HJ－1A星CCD數(shù)據(jù)進行交叉定標在一定程度上提高了水體反射率反演的精度。由于輻射定標是進行定量遙感需要解決的首要問題，因此本研究為今后利用HJ－1A星CCD數(shù)據(jù)進行近海水體懸浮泥沙等的定量反演提供了基礎(chǔ)。

（References）：

［1］LI C R，TANG L L，HU J，et al.Potential applications of HJ－1optical satellites［J］.Science ＆Technology Review，2008，26（13）：56－59.李傳榮，唐伶俐，胡堅，等.HJ－1光學(xué)衛(wèi)星應(yīng)用潛力［J］.科技導(dǎo)報，2008，26（13）：56－59.

［2］WANG Q，WU C Q，LI Q.Environment satellite 1and its application in environmental monitoring［J］.Journal of Remote Sensing，2010，1：104－118.王橋，吳傳慶，厲青.環(huán)境一號衛(wèi)星及其在環(huán)境監(jiān)測中的應(yīng)用［J］.遙感學(xué)報，2010，1：104－118.

［3］LIU X M，WANG Q，SUN Z P，et al.Assessment of ecosystem health of nature reserve based on HJ－1remote sensing imagery［J］.China Environmental Science，2011，31（5）：863－870.劉曉曼，王橋，孫中平，等.基于環(huán)境一號衛(wèi)星的自然保護區(qū)生態(tài)系統(tǒng)健康評價［J］.中國環(huán)境科學(xué)，2011，31（5）：863－870.

［4］XIA R，LI Y M，WU C Q，et al.Spatial distribution and variation of concentration of suspended solids in Taihu Lake based on HJ－1satellite data［J］.Scientia Geographica Sinica，2011，31（2）：197－203.夏叡，李云梅，吳傳慶，等.基于 HJ－1號衛(wèi)星數(shù)據(jù)的太湖懸浮物濃度空間分布和變異研究［J］.地理科學(xué)，2011，31（2）：197－203.

［5］KUANG D，HAN X Z，LIU X，et al.Quantitative estimation of Taihu chlorophyll－a concentration using HJ－1Aand 1BCCD imagery［J］.China Environmental Science，2010，30（9）：1268－1273.曠達，韓秀珍，劉翔，等.基于環(huán)境一號衛(wèi)星的太湖葉綠素a濃度提?。跩］.中國環(huán)境科學(xué)，2010，30（9）：1268－1273.

［6］WANG Z T，WANG Z F，LI Q，et al.Monitoring of PM＿（10）from HJ－1CCD data［J］.China Environmental Science，2011，31（2）：202－206.王中挺，王子峰，厲青，等.環(huán)境一號衛(wèi)星在監(jiān)測大氣PM＿（10）中的應(yīng)用［J］.中國環(huán)境科學(xué)，2011，31（2）：202－206.

［7］FENG L，WU W，CHEN X L，et al.Diseases and insect pests area monitoring for winter wheat based on HJ－1CCD imagery［J］.Transactions of the Chinese Society of Agricultural Engineering，2010，26（7）：213－219.馮煉，吳煒，陳曉玲，等.基于HJ衛(wèi)星CCD數(shù)據(jù)的冬小麥病蟲害面積監(jiān)測［J］.農(nóng)業(yè)工程學(xué)報，2010，26（7）：213－219.

［8］LI J G，GU X F，LI X Y，et al.Validation of HJ－1Bthermal band on－board calibration and its sensitivity analysis［J］.Remote Sensing Information，2011，1：3－8.李家國，顧行發(fā)，李小英，等.HJ－1B熱紅外通道星上定標精度檢驗與敏感性分析［J］.遙感信息，2011，1：3－8.

［9］HAN Q J，MIN X J，F(xiàn)U Q Y.Research of onboard calibration and cross－calibration of thermal infrared channel of HJ－1B［J］.Spacecraft Recoverty ＆ Remote Sensing，2009，30（4）：42－48.韓啟金，閔祥軍，傅俏燕.“環(huán)境一號”熱紅外通道星上定標與交叉定標研究［J］.航天返回與遙感，2009，30（4）：42－48.

［10］GUO X J，WANG Y C，WANG A H，et al.A strip noise removal method of hyperspectral data HJ－1A［J］.Remote Sensing Information 2011，1：54－58.郭興杰，王陽春，汪愛華，等.HJ－1A高光譜數(shù)據(jù)的條帶噪聲去除方法研究［J］.遙感信息，2011，1：54－58

［11］ADEL A M S，REINLDO A V.A statistical model for indoor multipath propagation［J］.IEEE Journal on selected areas in communications，1987，5（2）：128－137.

［12］TURIN G L，CLAPP F D，JOHNSON T L，et al.A statistical model of urban multipath propagation［J］.IEEE Trans.Veh.Techno.，1972，（21）：1－9.

［13］YANG Z D，GU S Y，QIU H，et al.CBERS－1’s CCD image quality evaluating and cross calibrating study［J］.Journal of Remote Sensing，2004，8（2）：113－120.楊忠東，谷松巖，邱紅，等.中巴地球資源衛(wèi)星一號衛(wèi)星CCD圖像質(zhì)量評價和交叉定標研究［J］.遙感學(xué)報，2004，8（2）：113－120.

［14］PENG G X，HE Y H，LI J，et al.Study on CBERS’s CCD image cross calibration and atmospheric correction［J］.J.Infrared Millim Waves，2007，26（1）：22－25.彭光雄，何宇華，李京，等.中巴資源地球02星CCD圖像交叉定標與大氣校正研究［J］.紅外與毫米波學(xué)報，2007，26（1）：22－25.

［15］GORDON H R，WANG M.Retrieval of water－leaving radiance and aerosol optical thickness over the oceans with SeaWiFS：apreliminary algorithm［J］.Applied Optics，1994，33（3）：443－42.

［16］GORDON H R.Radiative transfer in the atmosphere for correction of ocean color remote sensing［M］∥BARALE V，SCHLITTENHARDT P M.Ocean Color：Theory and Applications in a Decade of CZCS Experience.London：Kluwer Academic Publishers，1993：33－77.

［17］GORDON H R，CASTANO D J.Coastal zone color scanner atmospheric correction algorithm：Multiple scattering effects［J］.Applied Optics，1987，26（11）：2111－2122.

［18］TANG J W，MA C F，NIU S L，et al.The study on quantitative retrieval of water composition based on CBERS－02CCD camera data［J］.Science in China Ser.E：Information Sciences（Suppl.1），2005，35：156－170.唐軍武，馬超飛，牛生麗，等.CBERS－02衛(wèi)星CCD相機資料定量化反演水體成分初探.中國科學(xué)E輯：信息科學(xué)，2005，35（增刊1）：156－170.

［19］GORDON H R，CLARK D K.Clear water radiances for atmospheric correction of coastal zone color scanner image［J］.Applied Optics，1981，20（24）：4175－4180.

［20］NECKEL H，LABS D.The solar radiation between 3，300and 12，500［J］.Solar Physics，1984，90：205－258.

［21］TANG J W，TIAN G L，WANG X Y，et al.The methods of water spectra measurement and analysisⅠ：above－water method［J］.Journal of Remote Sensing，2004，8（1）：37－44.唐軍武，田國良，汪小勇，等.水體光譜測量與分析Ⅰ：水面以上測量法［J］.遙感學(xué)報，2004，8（1）：37－44.

［22］WANG M H，GORDON H R.A simple moderate accurate atmospheric correction algorithm for SeaWiFS.Remote Sensing of Environment，1994，50：231－239.

［23］LI J S，ZHANG B，CHEN Z C，et al.The study on atmospheric correction of CBERS image by using MODIS［J］.Science in China：Ser.E Technological Sciences，2006，36（Suppl.1）：141－150.李俊生，張兵，陳正超，等.MODIS數(shù)據(jù)輔助中巴資源衛(wèi)星圖像大氣校正研究.中國科學(xué)：E輯，2006，36（增刊1）：141－150.

［24］GONG S Q.The study on water constituents retrieval of Taihu Lake based on a semi－empirical model［D］.Nanjing：School of Geography Science in Nanjing Normal University，2007.龔紹琦.基于半分析模型的太湖水體組分濃度反演研究［D］.南京：南京師范大學(xué)地理科學(xué)學(xué)院，2007.

Precision Verifications of Radiance Calibration and Coastal Water Atmospheric Correction for Satellite HJ－1A

XIAO Yan－fang，ZHU Lin，ZHAO Wen－ji
（Key Lab of Three Dimension Information Acquisition and Application，Capital Normal University，Beijing 100048，China）

For the purpose of quantitative inversion in coastal waters by using the CCD data of Satellite HJ－1A，the precision of radiance calibration coefficients of the CCD mounted on Satellite HJ－1Awas validated and a cross radiance calibration was performed based on the data from Landsat－5TM.And then，an inversion was made of the reflectivity of coastal waters on the basis of all these data.By comparing the results from the inversion and those from the measurements，it has been found that the cross radiance calibration for the CCT data of Satellite HJ－1Acan，to a certain extent，improve the precision of the inversion of the reflectivity of coastal waters.These results indicate that the onboard radiance calibration coefficients of the CCD images on Satellite HJ－1Aare not precise so enough that the cross radiance calibration should be also applied when the CCD data from Satellite HJ－1Aare used to carry out the quantitative inversion in the coastal waters.

cross radiance calibration；onboard radiance calibration；reflectivity inversion；coastal waters；HJ－1A；Landsat－5TM

September 6，2010

TP79

A

1671－6647（2012）01－0054－09

2010－09－06

國家科技支撐計劃重點項目——環(huán)北京區(qū)域遙感綜合監(jiān)測信息服務(wù)系統(tǒng)（2007BAH15B02）

肖艷芳（1985－），女，山東濰坊人，博士研究生，主要從事生態(tài)遙感方面研究.E－mail：xiaoyanfang2006＠126.com

（張 騫 編輯）