王健，孫林

(山東科技大學(xué)測繪科學(xué)與工程學(xué)院，山東 青島 266590)

0 引言

地表反照率(Albedo)是指太陽短波波段(0.3～3μm)在半球空間的所有地表入射太陽輻射的反射通量與入射的太陽輻射通量之比。地表反照率在地球能量平衡(EBR)之中是一個重要參數(shù)，并且也是中長期天氣預(yù)報與全球變化研究中的重要參數(shù)，對地表反照率的研究就顯得尤為重要。[1－5]

當(dāng)前的地表反照率產(chǎn)品存在諸多不足，如受反演算法的影響，需要較長周期的觀測數(shù)據(jù)來進(jìn)行反演，時間分辨率低，受云的影響數(shù)據(jù)缺失較多等[6]。為了彌補(bǔ)這方面的不足，GLASS地表反照率產(chǎn)品應(yīng)運(yùn)而生，作為一個全新的科學(xué)數(shù)據(jù)產(chǎn)品，該產(chǎn)品通過算法改進(jìn)，使產(chǎn)品時間分辨率僅為1天，極大的提高了產(chǎn)品的使用范圍。

一種科學(xué)產(chǎn)品的應(yīng)用與推廣最主要的影響因素就為其產(chǎn)品精度，所以對GLASS數(shù)據(jù)的精度驗證就顯得尤為重要[7]。為了驗證GLASS全球地表反照率產(chǎn)品精度，文章使用全世界廣泛分布的以通量塔為基礎(chǔ)的全球通量觀測網(wǎng)絡(luò)(FLUXNET)站點地表實測數(shù)據(jù)，對多個區(qū)域不同地表覆蓋類型的GLASS地表反照率產(chǎn)品開展驗證工作。

1 GLASS反照率產(chǎn)品

由NASA發(fā)布的MCD43產(chǎn)品作為當(dāng)前應(yīng)用最廣泛的地表反照率產(chǎn)品，與其它遙感產(chǎn)品一樣，該地表反照率產(chǎn)品也會受到云像元的影響而出現(xiàn)地表數(shù)據(jù)的缺失現(xiàn)象，目前多數(shù)地表參數(shù)產(chǎn)品通過多天(比如8天、16天)數(shù)據(jù)合成從一定程度上減少缺失數(shù)據(jù)的比例，但是并不能徹底消除缺失。

GLASS地表反照率產(chǎn)品基本思想是通過融合不同數(shù)據(jù)源、不同反演算法產(chǎn)生的初級產(chǎn)品，開發(fā)與現(xiàn)有地表反照率遙感反演算法具有互補(bǔ)性的新反演算法，得到高精度、時空連續(xù)一致的全球反照率合成產(chǎn)品。GLASS地表反照率產(chǎn)品反演算法主要采用了基于MODIS地表反射率的反照率反演算法[8－9]AngularBin1(AB1)算法和基于MODIS大氣層頂反射率的反照率反演 AngularBin2(AB2)算法[10]。AB1算法是在太陽/觀測角度空間網(wǎng)格化的基礎(chǔ)上，建立MODIS波段地表方向反射率與寬波段反照率之間的線性回歸關(guān)系，訓(xùn)練數(shù)據(jù)來自POLDERBRDF數(shù)據(jù)集;AB2算法則進(jìn)一步直接建立MODIS波段大氣層頂方向反射率與寬波段反照率之間的線性回歸關(guān)系。

AB1算法的整體思路是建立MODIS地表方向反射率與地表寬波段反照率之間分格網(wǎng)(angular bin)的線性回歸關(guān)系[7]。可以實現(xiàn)利用Terra/Aqua平臺上的MODIS傳感器每天經(jīng)過大氣校正后的方向反射率數(shù)據(jù)直接反演地表寬波段反照率，從而大大縮短MODIS地表反照率產(chǎn)品生成所需的時間窗口，實現(xiàn)日地表反照率產(chǎn)品的生成。AB1算法在太陽/觀測角度空間網(wǎng)格化的基礎(chǔ)上，對每一個網(wǎng)格建立MODIS波段地表方向反射率與寬波段反照率之間的線性回歸關(guān)系。具有算法簡單、對輸入數(shù)據(jù)要求低的優(yōu)點，同時充分考慮了地表的BRDF和波譜特性。

AB2算法的整體思路是建立MODIS大氣層頂方向反射率與地表寬波段反照率之間分格網(wǎng)(angular bin)的線性回歸關(guān)系。該算法可以實現(xiàn)利用Terra/Aqua平臺上的MODIS傳感器每天獲取的大氣層頂方向反射率數(shù)據(jù)直接反演地表寬波段反照率，從而大大縮短MODIS地表反照率產(chǎn)品生成所需的時間窗口，實現(xiàn)日地表反照率產(chǎn)品的生成。

數(shù)據(jù)組在生產(chǎn)過程中收集了1985—1999年間的AVHRR和2000—2010年間的MODIS共15年的全球地表數(shù)據(jù)，其中基于 AVHRR數(shù)據(jù)產(chǎn)品的GLASS反照率產(chǎn)品分辨率為5km，基于MODIS數(shù)據(jù)產(chǎn)品為1km分辨率。

2 通量網(wǎng)(FLUXNET)

圖1 FLUXNET觀測站點全球分布圖

圖2 觀測塔站、區(qū)域網(wǎng)絡(luò)、全球網(wǎng)絡(luò)關(guān)系圖

FLUXNET是一個以全世界廣泛分布的通量塔為基礎(chǔ)的全球通量觀測網(wǎng)絡(luò)。通量塔網(wǎng)絡(luò)遍布世界各大主要國家，在北美、歐洲、亞洲、非洲都有它的子網(wǎng)絡(luò)存在(如美洲通量網(wǎng) Ameriflux，亞洲通量網(wǎng)Asiaflux等)。到2009年為止，共有超過500個可長期觀測數(shù)據(jù)的通量塔站點加入到 FLUXNET中，覆蓋的經(jīng)緯度范圍從北緯70°至南緯30°，用戶可以從中獲得這些站點的植被、土壤、水文、氣象等信息。其中圖1所示為FLUXNET站點各子網(wǎng)絡(luò)全球分布圖，從圖2中可以看到FLUXNET站點分三級網(wǎng)絡(luò)，自上而下依次為全球網(wǎng)絡(luò)FLUXNET，中間為區(qū)域網(wǎng)絡(luò)，如美洲通量網(wǎng)，中國通量網(wǎng)等。最底層為各觀測站點的通量塔。

FLUXNET提供的數(shù)據(jù)中有相當(dāng)一部分站點進(jìn)行了太陽短波輻射入射與出射的觀測(提供每半小時的觀測數(shù)據(jù))，可以應(yīng)用于地表反照率數(shù)據(jù)的驗證。

文章采用的地面驗證數(shù)據(jù)來自全球通量網(wǎng)輻射觀測站點。輻射觀測站點每10s觀測一次輻射通量數(shù)據(jù)，觀測數(shù)據(jù)為短波波段(280～3000nm)的上行和下行輻射通量。文章采用地方時11:30—12:30間的短波波段輻射通量數(shù)據(jù)計算反照率，把該時間段內(nèi)的上行和下行輻射通量取平均后做比值，得到的值即作為該地當(dāng)天的短波波段反照率真實值。

3 精度驗證

本研究共選擇4個FLUXNET站點的數(shù)據(jù)進(jìn)行驗證。這些站點分布在北美洲和亞洲，地表覆蓋類型包括草地和針葉以及混交林地。各站點的基本信息(如表1所示)，其中山東禹城站點和北京大興站點為農(nóng)田區(qū)域，美國華盛頓州的Wind River站點植被類型為常綠針葉林，威斯康辛州的Lost Creek植被類型為混交林。圖3a為禹城輻射觀測站點，圖3(b)為地表覆蓋物分布圖，圖4(a)為Wind River輻射觀測站點，圖4(b)為地表覆蓋物分布。

表1 各觀測站點基本信息

對選定的4個FLUXNET站點中的反照率驗證產(chǎn)品，按照影像獲取時間以及其經(jīng)緯度坐標(biāo)與站點實測數(shù)據(jù)獲取時間及所在位置匹配對應(yīng)并進(jìn)行對比驗證。提取像元的反照率值與地面觀測數(shù)據(jù)進(jìn)行長時間序列的對比分析。使用2005年全年觀測數(shù)據(jù)，剔除掉其中數(shù)據(jù)缺失，受云像元干擾以及有明顯錯誤的數(shù)據(jù)，平均每個站點的有效數(shù)據(jù)可達(dá)330天。各個站點分別比較驗證結(jié)果如圖5所示。

通過對比，我們可以看出在 Wind River站點(圖5(c))精度最高，均方根誤差約為0.0565，禹城站點(圖5(a))的精度其次，均方根誤差約為0.0776，Lost Creek 站點(圖 5(d))均方根誤差為0.0809，大興站點(圖5(b))的精度最低，均方根誤差為0.164。雖然同屬農(nóng)田區(qū)域大興站點與禹城相比卻有明顯的區(qū)別，在Google Earth上目視解譯可以明顯看出，大興地區(qū)的地表相對禹城來說，覆蓋類型并不均一。

圖3 山東禹城輻射觀測站點

圖4 Wind River輻射觀測站點

4 結(jié)論

文章首先利用全球通量網(wǎng)輻射站點的觀測數(shù)據(jù)，處理得到地面觀測反照率，作為地表反照率的真值。之后利用GLASS 2005全年的時間序列的反照率產(chǎn)品與對應(yīng)地點的實測產(chǎn)品進(jìn)行精度驗證。通過大量數(shù)據(jù)的對比表明，直接使用地表實測數(shù)據(jù)對GLASS反照率產(chǎn)品精度驗證，兩者的一致性較差。對于兩者數(shù)據(jù)的尺度方面，GLASS數(shù)據(jù)可以看做一個面像元，而實測數(shù)據(jù)僅為一個點像元，地面站點實測值與低空間分辨率反照率數(shù)據(jù)直接對比驗證，對于地表均一、大面積的區(qū)域驗證效果較好，對于混合像元嚴(yán)重的區(qū)域，二者的直接驗證和比較是不可行的。從Google Earth上對4各站點的目視解譯中我們看出，大興站點的地面覆蓋類型較復(fù)雜，當(dāng)使用這個點數(shù)據(jù)對面數(shù)據(jù)的驗證時，就會出現(xiàn)較大的誤差，另外三個站點的覆蓋較為均一，實測數(shù)據(jù)與反演數(shù)據(jù)較為接近。GLASS反照率產(chǎn)品的算法中使用了濾波融合，可能會濾掉一些符合情況但變化較大的點，相對地面觀測數(shù)據(jù)的一些峰值點GLASS表現(xiàn)的都較為平滑。

圖5 GLASS數(shù)據(jù)與實測數(shù)據(jù)對比圖

隨著遙感反演精度的進(jìn)一步提高和業(yè)務(wù)化的逐漸推進(jìn)，在不遠(yuǎn)的將來完全可以將衛(wèi)星遙感反照率產(chǎn)品直接引入模式應(yīng)用[10]，提高陸面過程模式和氣候模式的性能。

[1]Dickinson R..Land processes in climate models[J].Remote Sensing of Environment，1995(51):27 －38.

[2]Liang Shunlin，Wang Kaicun，Zhang Xiaotong，et al.Review of estimation of land surface radiation and energy budgets from ground measurements，remote sensing and modal simulation[J].IEEE JSTARS，2010，3(3):225 －240.

[3]鄭國強(qiáng)，魯敏，張濤，等，地表比輻射率求算對濟(jì)南市地表溫度反演結(jié)果的影響[J].山東建筑大學(xué)學(xué)報，2010，25(5):438－442.

[4]程潔，孫林，朱壽東，等，寬波段發(fā)射率對長波凈輻射估算影響分析[J]山東建筑大學(xué)學(xué)報，2012，27(4):355－359.

[5]馬春鋒，王維真，吳月茹，等.基于MODIS數(shù)據(jù)的黑河流域土壤熱慣量反演研究[J].遙感技術(shù)與應(yīng)用，2012，27(2):197－207.

[6]Liang Shunlin.A direct algorithm for estimating land surface broadband albedos from MODIS imagery[J].IEEE Transactions on Geosciences and Remote Sensing，2003(41):136－145.

[7]Russell M.，Nunez M.，Chladil M.，et al.Conversion of nadir narrowband reflectance in red and near－infrared channels to hemispherical surface albedo[J].Remote Sensing of Environment，1997(61):16－23.

[8]Schaaf C.，Gao Feng，Strahler A.，et al.First operational BRDF，albedo nadir reflectance products from MODIS[J].Remote Sensing of Environment，2002(83):135 －148.

[9]Lucht，W.，Schaaf，C.，Strahler，A..An algorithm for the retrieval of albedo from space using semiempirical BRDF models[J].Geoscience and Remote Sensing，IEEE Transactions on，2002(38):977－998.

[10]王介民，高峰.關(guān)于地表反照率遙感反演的幾個問題[J].遙感技術(shù)與應(yīng)用，2004，19(5):295 －300.