耦合京津冀氣溶膠模式的HJ-1衛(wèi)星CCD數(shù)據(jù)大氣校正

謝東海, 程天海, 吳 俁, 余 潔, 郭 紅, 鐘若飛

1.首都師范大學(xué)城市環(huán)境過程與數(shù)字模擬國家重點實驗室培育基地，北京 100048 2.中國科學(xué)院遙感與數(shù)字地球研究所，北京 100101 3.首都師范大學(xué)三維信息獲取與應(yīng)用教育部重點實驗室，北京 100048

耦合京津冀氣溶膠模式的HJ-1衛(wèi)星CCD數(shù)據(jù)大氣校正

謝東海1,3, 程天海2*, 吳 俁2, 余 潔1,3, 郭 紅2, 鐘若飛1,3

1.首都師范大學(xué)城市環(huán)境過程與數(shù)字模擬國家重點實驗室培育基地，北京 100048 2.中國科學(xué)院遙感與數(shù)字地球研究所，北京 100101 3.首都師范大學(xué)三維信息獲取與應(yīng)用教育部重點實驗室，北京 100048

針對目前HJ-1 CCD大氣校正沒有考慮中國地區(qū)氣溶膠模式的問題，提出一種耦合中國地區(qū)局部氣溶膠模式的大氣校正方法。以京津冀地區(qū)作為研究區(qū)域，該方法對地基北京城區(qū)和香河站點反演的氣溶膠模式參數(shù)進行聚類，得到京津冀地區(qū)具有代表性的四類氣溶膠模式，并根據(jù)四類氣溶膠模式來建立查找表進行氣溶膠光學(xué)厚度的反演。HJ-1 CCD數(shù)據(jù)沒有短波紅外波段(2.12 μm)，無法采用MODIS的氣溶膠算法中獲得地表反射率的方法來計算藍紅波段的反射率，本文在氣溶膠光學(xué)厚度的反演中采用HJ-1衛(wèi)星的藍色(0.43～0.52 μm)和紅色(0.63～0.69 μm)波段的反射率比值作為誤差方程的依據(jù)，不需要輸入地表目標(biāo)的反射率?；诜囱莺蟮墓鈱W(xué)厚度對HJ-1 CCD數(shù)據(jù)進行大氣校正，并與ASD光譜輻射計測量數(shù)據(jù)以及MODIS地表產(chǎn)品數(shù)據(jù)(MOD09)進行對比。結(jié)果表明，該方法得到的大氣校正結(jié)果與ASD測量結(jié)果接近，并與MOD09有較強的相關(guān)性，紅色波段的平均相關(guān)系數(shù)達到了0.8以上，受氣溶膠影響最大的藍色波段平均的相關(guān)系數(shù)也達到了0.75左右。

大氣校正； 氣溶膠； 京津冀地區(qū)； 環(huán)境一號衛(wèi)星

引 言

衛(wèi)星圖像的亮度值，不僅包括來自地表反射太陽的貢獻，還包括大氣中分子和氣溶膠的貢獻。大氣還會對地表反射的電磁波進行散射和吸收，導(dǎo)致地表反射光線在進入傳感器之前被衰減。根據(jù)模型的不同，大氣校正的方法可以分為三種[1]，基于圖像信息的方法，基于衛(wèi)星過境時同步觀測的地面數(shù)據(jù)的方法，基于輻射傳輸模型的方法?；趫D像信息的方法利用圖像自身的信息來進行大氣校正，不需要知道大氣參數(shù)，是最簡單的大氣校正方法，但缺點是校正的精度較低。代表性的方法是黑暗像元法[2-3]?；谛l(wèi)星過境時同步觀測的地面數(shù)據(jù)的方法利用野外實測的波譜信息，建立經(jīng)驗線性方程完成大氣校正，包括比值法，不變目標(biāo)法，直方圖匹配法和經(jīng)驗線性校正[4]?；谳椛鋫鬏斈Ｐ偷拇髿庑Ｕ椒ㄊ紫冗M行氣溶膠參數(shù)的反演，然后根據(jù)嚴(yán)格的物理模型來解算地表反射率，其校正的精度較高，能獲得地表的真實反射率(如MODIS的地表反射率產(chǎn)品MOD09)。

2008年9月6日，環(huán)境一號衛(wèi)星(代號HJ-1)一期工程“2+1”星座中的兩顆光學(xué)小衛(wèi)星(HJ-1A和HJ-1B)在太原衛(wèi)星發(fā)射中心“一箭雙星”成功發(fā)射。在經(jīng)歷了在軌測試后，2009年3月30正式交付衛(wèi)星用戶使用。環(huán)境一號衛(wèi)星是為環(huán)境監(jiān)測定做的衛(wèi)星，配備了寬覆蓋CCD相機、紅外相機、高光譜成像儀、合成孔徑雷達(SAR)等多種類型傳感器，星下點的分辨率為30 m，可以實現(xiàn)48 h重復(fù)觀測。其中CCD相機的有效視場角為31°，具有三個可見光波段，即藍波段(0.43～0.52 μm)，綠波段(0.52～0.60 μm)，紅波段(0.63～0.69 μm)和一個近紅波段(0.76～0.90 μm)。CCD相機的相對定標(biāo)精度指標(biāo)設(shè)計為5%，絕對定標(biāo)精度指標(biāo)設(shè)計為10%。

國內(nèi)學(xué)者針對HJ-1 CCD數(shù)據(jù)的大氣校正提出了多種模型和方法。王愛春[5]針對HJ-1 CCD的波段特點，對暗目標(biāo)法進行了改進，提出了利用比值植被指數(shù)(RVI)、土壤調(diào)整植被指數(shù)(SAVI)和歸一化水體指數(shù)(NDWI)的綜合分析法確定暗像元自動提取的方法。其查找表采用MODTRAN來構(gòu)建，考慮了三種氣溶膠模式：大陸模式、城市工業(yè)模式和沙塵模式。李莘莘基于北京與珠三角地區(qū)的地物光譜試驗，獲得暗像元的植被指數(shù)與紅藍波段反射率比值，構(gòu)建基于輻射傳輸模型的大氣校正算法。孫長奎[6]提出一種基于MODIS反射率庫的大氣校正方法，還從氣溶膠光學(xué)厚度的反演精度、大氣水汽含量的變化、輻射定標(biāo)精度、海拔高度等方面對大氣校正的不確定性進行了分析。其查找表采用6S[7]來構(gòu)建，氣溶膠模式只考慮了大陸型。方莉[8]提出一種基于純像元提取的氣溶膠光學(xué)厚度反演方法并應(yīng)用于大氣校正中，該方法的查找表采用6S來構(gòu)建。王中挺[9]采用深藍算法來反演氣溶膠光學(xué)厚度，然后對HJ-1號數(shù)據(jù)進行大氣校正，并利用核驅(qū)動模型實現(xiàn)了BRDF的校正。在誤差分析中, 研究了氣溶膠類型對大氣校正的影響，結(jié)果顯示近紅外波段的誤差較大，達到5%左右。其查找表也是用6S來建立。

目前的HJ-1 CCD數(shù)據(jù)大氣校正的方法都是基于輻射傳輸模型來實現(xiàn)的，氣溶膠光學(xué)厚度基于HJ-1衛(wèi)星CCD數(shù)據(jù)反演來獲得，反演過程都需要查找表的支持。但查找表的建立時只考慮了國外輻射傳輸軟件(MODTRAN，6S)中內(nèi)嵌的幾種通用的氣溶膠模式，沒有對中國區(qū)域的氣溶膠模式進行定量分析。由于氣溶膠反演是大氣校正的核心，其參數(shù)的精度直接影響了大氣校正的精度。而目前基于查找表的氣溶膠反演算法[10-11, 13-14, 19]中，氣溶膠模式是影響氣溶膠光學(xué)厚度反演精度的兩個重要因素之一(另外一個因素是地表反射率)。針對目前HJ-1數(shù)據(jù)大氣校正只考慮國外通用氣溶膠模式的問題，本文提出耦合中國地區(qū)局部氣溶膠模式的大氣校正方法，并以京津冀地區(qū)為例，對地基反演的氣溶膠模式參數(shù)進行聚類，得到待反演區(qū)域的典型氣溶膠模式并基于RT3[12]構(gòu)建查找表。以藍紅波段比值作為誤差依據(jù)，對HJ-1 CCD數(shù)據(jù)進行氣溶膠光學(xué)厚度的反演。然后將反演結(jié)果輸入大氣校正模塊，得到地表反射率值。最后利用ASD和MODIS地表反射率產(chǎn)品來對校正結(jié)果進行了驗證。

1 大氣校正的原理和方法

衛(wèi)星獲得的表觀反射率來源比較復(fù)雜，既有來自地表和大氣的直接反射，也有大氣多次散射，以及目標(biāo)鄰域的間接反射等。Tanré等[13]對上述各種貢獻進行抽象和綜合，得到表觀反射率的公式

(1)

(2)

(3)

(4)

2 京津冀地區(qū)氣溶膠模式分析

京津冀地區(qū)位于中國華北，包括北京、天津和河北省，是我國繼珠江三角洲和長江三角洲之后又一經(jīng)濟快速發(fā)展的地區(qū)。該地區(qū)春季易受來自北方的沙塵影響，夏季雨水較多，秋冬干燥。北京，天津是我國直轄市，機動車保有量達到數(shù)百萬量，河北省包含了我國重要的一些工業(yè)城市(如石家莊，唐山，保定等)。京津冀地區(qū)一年四季都會有大量來自人為活動的氣溶膠產(chǎn)生，大氣污染較為嚴(yán)重，氣溶膠的組成復(fù)雜而且四季變化明顯。Lee等[16]通過對東亞地區(qū)20個AERONET站點觀測數(shù)據(jù)的分析，聚類得到6種典型的氣溶膠模式。采用類似的方法，針對研究區(qū)域的地理范圍，選取AERONET在北京城區(qū)(116.381 37°，39.976 89°)和河北香河縣(116.961 50°，39.753 60°)設(shè)立站點的反演數(shù)據(jù)來對氣溶膠模式進行統(tǒng)計分析。

圖1中的填充區(qū)域為京津冀地區(qū)，兩個三角形是本文研究使用的AERONET地基觀測站點，北京站點位于北京城區(qū)，香河站點位于北京和天津之間。從圖1可以看出，河北省包圍了北京和天津市，這兩個站點位于河北省的地理中心區(qū)域，因此可以用來反映京津冀地區(qū)的氣溶膠分布特性。統(tǒng)計的數(shù)據(jù)為2013年1月到12月的氣溶膠反演數(shù)據(jù)，反演的級別為Level 1.5，北京站點能進行反演的天數(shù)為110 d，香河站點能進行反演的天數(shù)為78 d。

圖1 研究區(qū)域和地基站點

參考Omar[17]提出的聚類方法，對北京和香河兩個站點的氣溶膠模式進行自動聚類。參與聚類的參數(shù)包括譜分布參數(shù)(6個參數(shù))，復(fù)折射指數(shù)(4個波段共8個參數(shù))，單次散射反照率(SSA)(4個波段共4個參數(shù))，不對稱因子(4個波段共4個參數(shù))，參與聚類的參數(shù)數(shù)量為22個。本文自動聚類采用迭代自組織數(shù)據(jù)分析算法(iterative self-organizing data analysis technique algorithm，ISODATA)。ISODATA是廣泛使用的一種非監(jiān)督分類方法，通過設(shè)置的參數(shù)控制分裂與合并來自動進行聚類。算法中的關(guān)鍵輸入?yún)?shù)包括：初始的聚類數(shù)量(IC)，期待得到的聚類數(shù)量(EC)，最大迭代次數(shù)(MN)，每一類最少的樣本數(shù)(MS)，分裂閾值(ST)，合并閾值(MT)。根據(jù)輸入聚類算法的氣溶膠模式參數(shù)特點，本文ISODATA的輸入?yún)?shù)設(shè)置見表1。

表1 氣溶膠模式聚類分析的輸入?yún)?shù)

聚類時，地基數(shù)據(jù)每天的氣溶膠反演結(jié)果用一個22維的矢量來描述，不同矢量之間的距離公式采用最常用的歐式距離。

圖2為ISODATA算法自動聚類后得到的四種氣溶膠模式體積譜分布效果圖，表2中列舉了每個類別包含的樣本的數(shù)量。

圖2 聚類后得到的京津冀地區(qū)四種氣溶膠模式

表2 自動聚類后每類包含的樣本數(shù)量

從表2可以看出：第一類氣溶膠包括的樣本最多，即這一類氣溶膠在該地區(qū)出現(xiàn)的概率最大，其粗粒子的比例略高于細粒子。第二類氣溶膠的粗細粒子比例相當(dāng)，但是濃度都較大，意味著這類氣溶膠對應(yīng)的地基反演結(jié)果的光學(xué)厚度較大，主要來源是京津冀地區(qū)的污染天氣(比如霾)。第三類氣溶膠包括的樣本數(shù)量的比例最低，其中粗粒子的濃度要明顯高于細粒子，意味著該類氣溶膠中粗粒子占優(yōu)，主要來源于該地區(qū)的沙塵天氣。第四類氣溶膠的形狀與第一類比較相似，但粗粒子的比例要高于第一類，而且該類的粗粒子半徑是四個類別中最大的。本文將利用這四種氣溶膠類型的參數(shù)來建立查找表，并對京津冀地區(qū)的HJ-1號衛(wèi)星數(shù)據(jù)進行大氣校正。

3 大氣校正流程

在氣溶膠模式確定后，可以利用輻射傳輸軟件來生成查找表。已有的HJ-1大氣校正算法一般基于6S[7]，但是考慮到查找表的生成效率，本文采用RT3來生成查找表。RT3內(nèi)部沒有集成氣溶膠模式，而是直接將散射相函數(shù)的勒讓德系數(shù)作為輸入。散射相函數(shù)根據(jù)粒子的形狀可以利用MIE代碼或者非球形粒子相函數(shù)解算軟件(如T-Matrix[18])來解算。RT3的優(yōu)點是能夠靈活的集成不同類型的氣溶膠模式，且計算速度快。查找表生成后，可以根據(jù)圖3所示的流程來實現(xiàn)大氣校正。

圖3 大氣校正流程

從圖3可以看出，基于輻射傳輸?shù)墓鈱W(xué)衛(wèi)星影像大氣校正的前提是獲取氣溶膠的模式參數(shù)以及光學(xué)參數(shù)(氣溶膠光學(xué)厚度)，這些參數(shù)可以利用氣溶膠反演來獲取。由于HJ-1衛(wèi)星缺少2.12 μm波段，不能根據(jù)經(jīng)驗公式來獲得藍紅波段的地表反射率[14]，所以本文基于反射率比值[19]來建立誤差方程進行反演。對MODIS數(shù)據(jù)進行大氣校正并進行大量統(tǒng)計表明，地表藍紅波段反射率的比值為0.49，而且?guī)缀醪皇苌⑸浣嵌鹊挠绊慬14-15, 20]，因此誤差方程可以寫為

(5)

該方法的優(yōu)點是不需要知道紅藍波段反射率的真實值，只需要知道它們之間比值即可。本文采用NDVI來判斷暗目標(biāo)，由于NDVI會受到大氣的影響，因此我們選擇晴朗天氣的數(shù)據(jù)來計算NDVI?？紤]到暗目標(biāo)在一定時間段內(nèi)不會變化(如8 d或者一個月)，本文利用HJ-1歷史數(shù)據(jù)來計算某段時間內(nèi)晴朗天氣的NDVI平均值。在對衛(wèi)星影像進行逐點的大氣校正時，需要知道每個地表目標(biāo)上空的氣溶膠光學(xué)厚度。本文的氣溶膠反演算法只能反演暗目標(biāo)上空的氣溶膠光學(xué)厚度，考慮到氣溶膠在空間分布上具有一定的連續(xù)性，采用不規(guī)則點的插值方法來進行插值得到連續(xù)分布的氣溶膠光學(xué)厚度。

4 結(jié)果與討論

4.1 京津冀地區(qū)大氣校正

本文選擇2010年8月16日的HJ-1衛(wèi)星數(shù)據(jù)來進行大氣校正。為了研究京津冀地區(qū)的大氣校正效果，首先從HJ-1 CCD圖像中將京津冀地區(qū)裁剪出來。由于HJ-1衛(wèi)星的覆蓋范圍有限，一景圖像不能覆蓋京津冀地區(qū)，所以選擇連續(xù)拍攝的三景圖像來進行實驗。氣溶膠的反演以及大氣校正分別在每景圖像中完成，然后將三景圖像拼接起來得到整體的反演和大氣校正圖。圖4是將三景影像拼接后得到的氣溶膠光學(xué)厚度，以及基于氣溶膠光學(xué)厚度進行逐點大氣校正后得到的地表反射率圖。從圖4可以看出，原始的圖像由于氣體分子和氣溶膠的影響，對比度較低，而且整體偏藍。利用本文的大氣校正方法能夠較好的恢復(fù)地表的反射率信息，植被在校正后恢復(fù)為正常的顏色。

4.2 與地面測量比較

2013年7月28日，我們在北京奧林匹克森林公園的南門用ASD光譜輻射計對草地和水泥地目標(biāo)進行了測量。為了配合HJ-1衛(wèi)星過境的時間，測量時間選擇在上午10：30左右。與該天時間最接近，而且覆蓋了北京地區(qū)的HJ-1號衛(wèi)星數(shù)據(jù)是在2013年7月24日拍攝得到的。采用該天的數(shù)據(jù)來進行大氣校正(北京地區(qū)進行大氣校正后的效果見圖5)，然后選取相同位置的數(shù)據(jù)來與ASD光譜輻射計測量結(jié)果進行對比。

圖4 拼接后的大氣校正效果圖

圖5 北京地區(qū)大氣校正效果圖(2013.7.24)

圖6中的藍色曲線是ASD光譜輻射計測量的結(jié)果，圓點是HJ-1 CCD大氣校正后的反射率。橫坐標(biāo)為波長，ASD測量的范圍為350～2 500 nm，考慮到HJ-1 CCD相機只有四個波段，圖中截取了350～1 000 nm的范圍來顯示。從圖6可以看出，HJ-1號衛(wèi)星的四個波段大氣校正后的反射率數(shù)據(jù)與ASD觀測的結(jié)果接近。

4.3 與MODIS地表反射率產(chǎn)品進行對比

選擇2010年8月16日，9月30日和10月29日的HJ-1衛(wèi)星大氣校正后的反射率數(shù)據(jù)來與MODIS地表反射率產(chǎn)品進行對比。MODIS根據(jù)8天數(shù)據(jù)來進行大氣校正后得到了地表反射率產(chǎn)品(MOD09)，其中500 m的產(chǎn)品中包括了前7個波段。京津冀地區(qū)的MOD09(500 m)的效果見圖7(a)。

為了方便進行對比，將HJ-1的大氣校正數(shù)據(jù)的分辨率從30 m降低到500 m。對比前首先進行云檢測，方法是如果藍，綠，紅三個波段反射率都大于0.2，就識別為云。圖7為2010年8月16日的環(huán)境衛(wèi)星大氣校正后的反射率圖像[圖7(b)]與對應(yīng)的MOD09產(chǎn)品[圖7(a)]的對比圖，其余兩天的數(shù)據(jù)也采用相同的方法來進行處理。

圖6 HJ-1衛(wèi)星影像大氣校正結(jié)果與ASD測量結(jié)果對比

圖7 2010.8.16京津冀地區(qū)MODIS與HJ-1衛(wèi)星反射率對比圖

圖8為本文大氣校正后的HJ-1可見光數(shù)據(jù)與對應(yīng)的MOD09數(shù)據(jù)的相關(guān)性分析結(jié)果。根據(jù)計算的相關(guān)系數(shù)可以看出，紅波段的大氣校正后的反射率與MOD09的相關(guān)性整體比較高，三天的數(shù)據(jù)都達到了0.8以上。藍色波段大氣校正后的反射率與MOD09的相關(guān)性整體最低，三天平均后為0.75左右。綠色波段的相關(guān)性在2010年8月16日和9月30日都位于藍紅波段之間，但是10月29日的相關(guān)性最高，達到了0.893。紅波段大氣校正效果較好的原因是該波段的氣溶膠光學(xué)厚度較低，大氣影響較小。而藍色波段的氣溶膠光學(xué)厚度最高，因此相比紅色波段，藍色波段光學(xué)厚度的反演精度對大氣校正的影響較大。MODIS數(shù)據(jù)的質(zhì)量較高，其產(chǎn)品在定量遙感領(lǐng)域也被廣泛使用。本文提出的大氣校正算法與MOD09數(shù)據(jù)能夠達到較高的相關(guān)性，證明了本文算法的有效性。

圖8 HJ-1與MOD09京津冀地區(qū)地表反射率對比散點圖, r為相關(guān)系數(shù)

5 結(jié) 論

在對京津冀地區(qū)的氣溶膠模式進行分析的基礎(chǔ)上，提出了一種耦合局部氣溶膠模式的HJ-1衛(wèi)星大氣校正方法，并利用地面ASD測量數(shù)據(jù)和MOD09數(shù)據(jù)對校正結(jié)果進行了初步的驗證。在基于查找表的氣溶膠光學(xué)厚度反演算法中，氣溶膠模式會影響到氣溶膠光學(xué)厚度反演的精度。通過對地基反演數(shù)據(jù)的分析而獲得的氣溶膠模式參數(shù)更接近實際的氣溶膠分布，驗證分析表明，將局部氣溶膠模式與輻射傳輸進行耦合后對HJ-1衛(wèi)星數(shù)據(jù)進行大氣校正，能得到較好的校正效果。后續(xù)的研究中，需要對算法繼續(xù)進行改進。具體包括：采用更長時間序列的地基數(shù)據(jù)來對氣溶膠模式進行分析，以得到更準(zhǔn)確的氣溶膠模式參數(shù)； 利用大量的數(shù)據(jù)來進行光學(xué)厚度的反演，然后與地基數(shù)據(jù)進行對比，對氣溶膠光學(xué)厚度的精度進行分析； 利用模擬數(shù)據(jù)來定量分析氣溶膠模式對大氣校正的精度影響。

[1] Hadjimitsis D G, Clayton C R I，Hope V S.International Journal of Remote Sensing, 2004, 25(18): 3651.

[2] Liang S, Fang H，Chen M.IEEE Transactions on Geoscience and Remote Sensing, 2001, 39(11): 2490.

[3] Hadjimitsis D G, et al.Natural Hazards and Earth System Science, 2010，10(1): 89.

[4] ZHAO Ying-shi, et al(趙英時，等).Principles and Methods of Remote Sensing Application and Analysis(遙感應(yīng)用分析原理與方法).Beijing：Science Press(北京：科學(xué)出版社)，2013.502.

[5] WANG Ai-chun, et al(王愛春，等).Science China: Information Sciences(中國科學(xué): 信息科學(xué)), 2011，(S1): 76.

[6] SUN Chang-kui, et al(孫長奎, 等).Journal of Remote Sensing(遙感學(xué)報), 2012, (4): 826.

[7] Vermote E F, et al.IEEE Transactions on Geoscience and Remote Sensing, 1997，35(3): 675.

[8] FANG Li, et al(方 莉，等).Journal of Remote Sensing(遙感學(xué)報), 2013, (1): 151.

[9] WANG Zhong-ting, et al(王中挺，等).Spectroscopy and Spectral Analysis(光譜學(xué)與光譜分析), 2014, 34(3): 729.

[10] Kaufman Y J, et al.Journal of Geophysical Research, 1997，102(D14): 17051.

[11] Levy R C, et al.Atmospheric Measurement Techniques Discussions, 2013，6(1): 159.

[12] Levy R C, Remer L A，Dubovik O.Journal of Geophysical Research: Atmospheres, 2007，112(D13): .

[13] Remer L A, et al.Journal of the Atmospheric Sciences, 2005，62(4): 947.

[14] Kaufman Y J, et al.IEEE Transactions on Geoscience and Remote Sensing, 1997，35(5): 1286.

[15] Evans K F，Stephens G L.Journal of Quantitative Spectroscopy and Radiative Transfer, 1991, 46(5): 413.

[16] Tanre D, Herman M，Deschamps P Y.Applied Optics, 1981，20(20): 3676.

[17] Lee K H，Kim Y J.Atmospheric Measurement Techniques, 2010, 3(6): 1771.

[18] Omar A H.Journal of Geophysical Research, 2005，110(D10).

[19] Mishchenko M I，Travis L D.Journal of Quantitative Spectroscopy & Radiative Transfer, 1998，60(3): 309.

[20] SHENG Li, et al(盛 莉，等).Journal of Infrared and Millimeter Waves(紅外與毫米波學(xué)報), 2013，(3): 277.

*Corresponding author

(Received Dec.30, 2014; accepted Mar.6, 2015)

Atmospheric Correction for HJ-1 CCD Data Coupling with Aerosol Models of Beijing-Tianjin-Hebei Region

XIE Dong-hai1, 3, CHEN Tian-hai2*, WU Yu2, YU Jie1, 3, GUO Hong2, ZHONG Ruo-fei1, 3

1.State Key Laboratory Incubation Base of Urban Environmental Processes and Digital Simulation, Capital Normal University, Beijing 100048, 2.State Key Laboratory of Remote Sensing Science, Jointly Sponsored by the Institute of Remote Sensing Applications of Chinese Academy of Sciences and Beijing Normal University, Beijing 100101, 3.Beijing Key Laboratory of Resource Environment and Geographic Information System, Capital Normal University, Beijing 100048, China

Current atmospheric correction of HJ-1 CCD does not take the aerosol models of Chinese specific regions into consideration.This paper proposes a method of atmospheric correction coupling local aerosol models.Choosing Beijing-Tianjin-Hebei region as the study area, aerosol model parameters of the method on the foundation of AERONET inversion are analyzed and four types of aerosol model parameters of Beijing-Tianjin- Hebei regions are clustered to build the lookup table for the inversion of aerosol optical thickness.Lack of SWIR (2.12 μm) in the bands of HJ-1 CCD , the blue(0.43～0.52 μm) and red(0.63～0.69 μm) cannot be calculated like the MODIS.So the ratio of blue and red is used as the basis of error equation for aerosol optical depth retrieval, with no need of ground target reflectance.HJ-1 CCD data after atmospheric correction are compared with MODIS surface product data (MOD09) and ASD measurements.The results show that the reflectance obtained by the atmospheric correction is close to the ASD measurement results, and there is a strong correlation with MOD09, the average correlation coefficient in the red band reached more than 0.8, the one of blue band affected by molecule mostly is up to 0.75.

Atmospheric correction; Aerosol; Beijing-Tianjin-Hebei region; HJ-1 satellite

2014-12-30，

2015-03-06

國家自然科學(xué)基金項目(41401386)資助

謝東海，1978年生，首都師范大學(xué)資源環(huán)境與旅游學(xué)院講師 e-mail: xdhbj@126.com *通訊聯(lián)系人 e-mail: cthy@irsa.ac.cn

P407.6

A

10.3964/j.issn.1000-0593(2016)05-1284-07