李朋澤，王新生,2，郭光毅，徐 靜，許 亮，謝 凱

（1. 湖北大學(xué) 資源環(huán)境學(xué)院，湖北 武漢430062；2. 農(nóng)業(yè)部 遙感應(yīng)用中心 武漢分中心，湖北 武漢 430062）

基于MODIS的江漢平原比輻射率季節(jié)變化規(guī)律研究

李朋澤1，王新生1,2，郭光毅1，徐 靜1，許 亮1，謝 凱1

（1. 湖北大學(xué) 資源環(huán)境學(xué)院，湖北 武漢430062；2. 農(nóng)業(yè)部 遙感應(yīng)用中心 武漢分中心，湖北 武漢 430062）

基于Terra-MODIS L3級(jí)產(chǎn)品（MOD11C3），研究江漢平原區(qū)域比輻射率在一年四季的變化特征。結(jié)果表明，江漢平原區(qū)域的比輻射率值在1 a內(nèi)整體處于一個(gè)穩(wěn)定的區(qū)間,相比其他三季，冬季高比輻射率區(qū)間所占面積有明顯增加，同時(shí)存在多種類別的比輻射率值；春季和冬季的平均比輻射率值最高，夏季最低，秋季比輻射率值居夏冬兩季之中。對(duì)不同土地覆被類型的比輻射率值的研究表明，水田和水域的比輻射率值最高，建設(shè)用地比輻射率值最低但波動(dòng)較大，林地和草地的比輻射率值最為穩(wěn)定，耕地的比輻射率值明顯高于自然植被；各種地類的比輻射率年變化趨勢(shì)基本一致，均在夏季出現(xiàn)波谷，秋季呈上升趨勢(shì)，冬季有小幅的波動(dòng)，春季最穩(wěn)定。

江漢平原；MODIS；地表比輻射率；季節(jié)變化 ；土地覆被類型

1 MODIS數(shù)據(jù)及其預(yù)處理

中分辨率成像光譜儀MODIS 是Terra和Aqua衛(wèi)星上搭載的主要傳感器之一，其標(biāo)準(zhǔn)數(shù)據(jù)產(chǎn)品根據(jù)內(nèi)容的不同分為0級(jí)、1級(jí)數(shù)據(jù)產(chǎn)品，在1 B級(jí)數(shù)據(jù)產(chǎn)品之后，劃分2～4級(jí)數(shù)據(jù)產(chǎn)品，包括陸地標(biāo)準(zhǔn)數(shù)據(jù)產(chǎn)品、大氣標(biāo)準(zhǔn)數(shù)據(jù)產(chǎn)品和海洋標(biāo)準(zhǔn)數(shù)據(jù)產(chǎn)品等3種主要標(biāo)準(zhǔn)數(shù)據(jù)產(chǎn)品類型，總計(jì)分解為44種標(biāo)準(zhǔn)數(shù)據(jù)產(chǎn)品類型。本文使用的MOD11C3數(shù)據(jù)為陸地3級(jí)標(biāo)準(zhǔn)數(shù)據(jù)產(chǎn)品，內(nèi)容為地表溫度和輻射率的全球0．05°×0．05°格網(wǎng)月值數(shù)據(jù)集[1-7]，選取的時(shí)間范圍為2011-03～2014-02。

運(yùn)用MRT（MODIS reprojection tool）軟件將數(shù)據(jù)重新投影為Albers等積投影下的GeoTIFF格式，輸出所需的31、32兩波段值，通過Erdas軟件對(duì)所得數(shù)據(jù)進(jìn)行轉(zhuǎn)換處理，得到兩波段的比輻射率值：

ε31= DN × 0.002 + 0.49 (1)

再采用梁順林提出的線性擬合方法得出寬波段的比輻射率[8]，線性擬合公式為：

ε = 0.314ε31+ 0.411ε32+ 0.261 (2)式中，ε31和ε32分別為MODIS31波段和MODIS32波段的比輻射率；ε為寬波段的比輻射率，通過裁剪得到江漢平原區(qū)域比輻射率月產(chǎn)品時(shí)序圖像。

2 研究區(qū)概況

江漢平原位于湖北省中南部，由長江與漢江沖積而成，是中國海拔最低的平原之一。介于北緯29°26'～31°10'，東經(jīng)111°45'～114°16'。在行政區(qū)劃上含荊州市（原沙市市、江陵縣）、石首、漢川、監(jiān)利、潛江、云夢(mèng)、松滋、應(yīng)城、天門、當(dāng)陽、洪湖、仙桃、公安、枝江，總面積為28 235．12 km2，是長江中游重要的產(chǎn)糧區(qū)和農(nóng)副產(chǎn)品基地。該區(qū)氣候溫暖濕潤，無霜期長達(dá)243～275 d；活動(dòng)積溫5 100℃～5 400℃；年降水量為980～1 300 mm，4～9月降水量占全年雨量2/3以上，雨熱同期[9]，植被呈現(xiàn)明顯的季節(jié)變化。

3 江漢平原比輻射率季節(jié)變化規(guī)律

通過對(duì)數(shù)據(jù)的處理分析，得到了2011-03～2014-02每月的比輻射率值。為了對(duì)江漢平原區(qū)域進(jìn)行季相分析，選取每年的3～5月、6～8月、9～11月、12月～次年2月作為春、夏、秋、冬四季，所以時(shí)間上是從2011-03開始到2014-02結(jié)束。

將每月的比輻射率值做成折線圖（見圖1），可發(fā)現(xiàn)江漢平原區(qū)域的比輻射率值在1 a內(nèi)呈現(xiàn)“兩頭高中間低”的趨勢(shì)：春冬兩季明顯比夏秋兩季更高；春季的比輻射率最穩(wěn)定，變化不大，夏季起伏最大。

圖1 江漢平原地表比輻射率月份變化

在對(duì)江漢平原區(qū)域進(jìn)行季相分析之前，需要將各季節(jié)對(duì)應(yīng)的3個(gè)月份比輻射率值做平均值，得到對(duì)應(yīng)年份的四季值。為了能夠更好地表現(xiàn)出比輻射率的變化，結(jié)合月比輻射率值的實(shí)際情況，將其分為5個(gè)區(qū)間：低比輻射率區(qū)（0．970 0～0．971 5）；中低比輻射率區(qū)（0．971 6～0．972 5）；中比輻射率區(qū)（0．972 6～0．073 5）；中高比輻射率區(qū)（0．973 6～0．974 5）；高比輻射率區(qū)（0．974 6～0．976 1）。

通過對(duì)比分析3年各個(gè)季節(jié)的比輻射率平均值可以發(fā)現(xiàn)，江漢平原區(qū)域的比輻射率值基本處于一個(gè)穩(wěn)定的區(qū)間（0．970 0～0．980 0之間），各季節(jié)的分布范圍也比較固定，主要原因可能與其地勢(shì)平坦且一年四季農(nóng)作物覆蓋廣泛、土壤水分充足有一定的關(guān)系?？v向?qū)Ρ韧攴莶煌竟?jié)的比輻射率，由圖2可以看出，0．973 6～0．974 5的中高值區(qū)和0．974 6～0．976 1的高值區(qū)所占的面積比例最大，二者之和可以達(dá)到90%以上；冬季相比其他三季，不同區(qū)間比輻射率的分布更加豐富一些，不過仍以中高區(qū)和高區(qū)為主，高比輻射率區(qū)間所占面積有明顯增大。

圖2 江漢平原（2011～2013年）四季比輻射率

橫向?qū)Ρ炔煌攴菹嗤竟?jié)的各個(gè)比輻射率區(qū)間的空間分布可以發(fā)現(xiàn)，雖然存在一定的差異，但在整體上具有一定的吻合度。例如在冬季，比輻射率的高值區(qū)域與低值區(qū)域很明顯在一定程度上保持一致。為了能進(jìn)一步探討江漢平原區(qū)域比輻射率在季節(jié)上的變化規(guī)律，對(duì)各個(gè)季節(jié)的比輻射率圖計(jì)算出均值，制作均值在不同季節(jié)變化的折線圖（如圖3）。

由圖3可以看出，江漢平原區(qū)域比輻射率四季的均值整體維持在一個(gè)穩(wěn)定的范圍（0．972 0～0．976 0之間），3條曲線的變化呈現(xiàn)出一定的規(guī)律：無論是哪一年，曲線的波峰都是在春季和冬季，波谷都出現(xiàn)在夏季，秋季比輻射率值則都處在了夏冬兩季的中間區(qū)域。夏季是我國地表溫度達(dá)到最高值的時(shí)間段，江漢平原位于長江中下游區(qū)域，更是我國的高溫區(qū)，但其比輻射率值卻在四季中降到谷底，相反春冬兩季地表溫度低得多，其比輻射率值卻上升到一定的高度，氣溫的變化或許對(duì)地物比輻射率的變化產(chǎn)生了一定的影響。

圖3 江漢平原（2011～2013年）四季平均比輻射率變化曲線

比輻射率的變化也與環(huán)境因素的變化有關(guān)系。雖然冬季樹木草地枯萎，農(nóng)作物也收割完畢，本應(yīng)隨著裸地面積的增加比輻射率下降，但由于江漢平原是我國重要的農(nóng)作物種植基地，農(nóng)田在冬季一般都會(huì)進(jìn)行灌水或者保持很高水分，使農(nóng)田比輻射率表現(xiàn)出水的特征，而水具有很高的比輻射率，導(dǎo)致冬季比輻射率整體較高。到了夏季，隨著農(nóng)作物和植被的生長，下墊面產(chǎn)生了改變，江漢平原以農(nóng)作物為主，從而比輻射率下降。

4 江漢平原不同土地覆被類型比輻射率變化分析

地表比輻射率除具有明顯的波譜特征外[10]，其主要取決于植被覆被、土壤濕度、土壤紋理、礦物質(zhì)組分以及冰雪[11]。不同的土地覆被類型是這些因素的綜合反映結(jié)果。地表比輻射率隨著土地覆被類型變化而改變，并與下墊面植被和土壤比例密切相關(guān)[12,13]。因此，為了更好地探討江漢平原區(qū)域地表比輻射率的變化機(jī)理，將江漢平原土地類型分為城鎮(zhèn)用地、耕地、林地、水域、草地和未利用地6大類，其中，耕地又分成水田和旱地2類，研究江漢平原區(qū)域不同的土地覆被類型對(duì)其比輻射率變化產(chǎn)生的影響。

由圖4可以看出，耕地占江漢平原區(qū)域大部分面積，經(jīng)統(tǒng)計(jì)達(dá)到了70%。所以將耕地分為水田和旱地2 類進(jìn)行分析，能更加準(zhǔn)確地對(duì)其變化規(guī)律進(jìn)行判斷。這里，水田指有水源保證和灌溉設(shè)施，在一般年景能正常灌溉，用以種植水稻、蓮藕等水生農(nóng)作物的耕地，包括實(shí)行水稻和旱地作物輪種的耕地；旱地指無灌溉水源及設(shè)施，靠天然降水生長作物的耕地，或有水源和澆灌設(shè)施，在一般年景下能正常灌溉的旱作物耕地，以及以種菜為主的耕地，正常輪作的休閑地和輪歇地。

圖4 江漢平原土地分類圖

水田和水域的比輻射率值最高；建設(shè)用地比輻射率值最低，且其波動(dòng)最大；林地和草地的比輻射率值最為穩(wěn)定；耕地的比輻射率值明顯高于自然植被；各種地類的比輻射率月變化趨勢(shì)一致，均在夏季出現(xiàn)波谷，秋季呈上升趨勢(shì)，冬季有小幅的波動(dòng)，春季最穩(wěn)定。由于不同的土地利用類型的比輻射率表現(xiàn)出的波動(dòng)性除了與自身組成物質(zhì)的性質(zhì)有關(guān)，還受人類活動(dòng)影響，人類活動(dòng)程度越大，其波動(dòng)的程度越明顯[14]，如圖5所示。

通過耕地發(fā)現(xiàn)，在4月份有個(gè)升高的小趨勢(shì)，這是由于農(nóng)作物返青的原因。5月份降低，是因?yàn)槎咀魑锏氖崭睢＝瓭h平原大范圍種植水稻，水稻種植時(shí)，稻田被水面覆蓋，這些農(nóng)用地會(huì)表現(xiàn)為水的特征。農(nóng)用地的比輻射率在7、8月份有個(gè)急劇降低的過程，可能與夏季雨水多、云的影響，或夏季水汽含量高，蒸騰作用強(qiáng)等原因有關(guān)，需要進(jìn)一步研究。

圖5 江漢平原不同土地類型比輻射率變化

根據(jù)覃志豪等人研究成果，水體、植被、裸土和建筑的比輻射率分別可以用0．995、0．986、0．972 15和0．97來作為估計(jì)[15]。以此作為參考，可以發(fā)現(xiàn)MODIS Day/Light算法對(duì)水域的計(jì)算偏低，水體的比輻射率值基本小于農(nóng)用地。整體在夏季值偏低，表現(xiàn)出較小的季節(jié)性?；谶@種算法得到的比輻射率對(duì)土地類型比較敏感，各地類類型比輻射率值界限明顯，但同一地類的比輻射率差異小。農(nóng)作物在不同生長期變化大，特別是水稻種植，稻田被水面覆蓋，表現(xiàn)為水的特征，與林地和草地相比，變化更大。因此MODIS Day/ Light 算法結(jié)果用來計(jì)算農(nóng)用地比輻射率值以及溫度會(huì)帶來一定的誤差。

5 結(jié) 語

江漢平原區(qū)域四季的平均比輻射率值具有一定的規(guī)律性：春季和冬季的平均比輻射率值最高，夏季最低，秋季比輻射率值在夏冬兩季的中間區(qū)域。MODIS Day/ Light算法得出的比輻射率對(duì)土地類型比較敏感，但表現(xiàn)出較小的季節(jié)性。對(duì)于大面積農(nóng)用地來說，農(nóng)作物在不同生長期變化大，MODIS Day/Light 算法得到的MODIS 6 km產(chǎn)品可以大面積獲取地表比輻射率，大尺度監(jiān)測地球溫度信息，其影像和產(chǎn)品意義重大，但MODISC3產(chǎn)品的分辨率低，只能滿足全球范圍的溫度估計(jì)與分析，對(duì)于區(qū)域范圍的比輻射率的反映不夠靈敏。因此，建立不同地表下墊面比輻射率遙感反演算法具有十分重要的意義。

[1] 沈斌, 閻廣建, 馮雪，等． 熱紅外比輻射率光譜野外測量方法的對(duì)比研究[J]． 北京師范大學(xué)學(xué)報(bào)：自然科學(xué)版, 2007, 43(3): 264-368[2] 梁順林． 定量遙感[M]． 北京：科學(xué)出版社, 2009

[3] 徐靜, 王新生, 高守杰，等． 地物比輻射率數(shù)據(jù)分析[J]． 遙感技術(shù)與應(yīng)用, 2013, 28(5):815-823

[4] 劉三超, 柳欽火, 高懋芳，等． 波譜響應(yīng)函數(shù)和波寬對(duì)地表溫度反演的影響[J]． 遙感信息, 2007(5): 3-6

[5] 張仁華． 對(duì)于定量熱紅外遙感的一些思考[J]． 國土資源遙感, 1999 (1): 1-6

[6] 王新生, 徐靜, 柳菲, 等． 近十年我國地表比輻射率的時(shí)空變化[J]． 地理學(xué)報(bào), 2012, 67(1): 93-100

[7] 侯鵬, 曹廣真, 蔣衛(wèi)國，等． 城市復(fù)雜地表TM溫度反演及其與MODIS產(chǎn)品的比較[J]． 自然災(zāi)害學(xué)報(bào), 2009, 18(5):113-119

[8] Wan Z M． New Refinements and Validation of the MODIS Land-surface Temperature/Emissivity Products[J]． Remote Sensing of Environment, 2008(112)：59-74

[9] 李仁東, 程學(xué)軍, 隋曉麗． 江漢平原土地利用的時(shí)空變化及其驅(qū)動(dòng)因素分析[J]． 地理研究, 2003, 22(4): 423-432

[10] Nerry F, Labed J, Stoll MP． Spectral Properties of Land Surfaces in the Thermal Infrared ． 1． Laboratory Measurements of Absolute Spectral Emissivity Signatures[J]． Journal of Geophysical Research Solid Earth and Planets, 1990, 95(B5):7 027-7 044

[11] Van de Griend A, Owe M , Groen M, et al． Measurement and Spatial Variation of Thermal Infrared Surface Emissivity in a Savanna Environment[J]． Water Resources Research, 1991, 27(3):371-379

[12] Pinker RT, Sun D L, Hung M P, et al． Evaluation of Satellite Estimates of Land Surface Temperature from Goes over the Anited States[J]． Journal of Applied Meteorology and Climatology, 2009, 48(1):167-180

[13] 翟俊, 劉紀(jì)元． 2000～2011年中國地表比輻射率時(shí)空格局及影響因素分析[J]． 資源科學(xué), 2013, 35(10): 2 094-2 103

[14] 侯光雷, 張洪巖, 王野喬,等． 基于MODIS數(shù)據(jù)的吉林省中部地表溫度反演及空間分布研究[J]． 地理科學(xué), 2010, 30(3):421-427

[15] 覃志豪, 李文娟, 徐斌,等． 陸地衛(wèi)星TM6波段范圍內(nèi)地表比輻射率的估計(jì)[J]． 國土資源遙感, 2004(3):27-32

P237.9

B

1672-4623（2015）05-0069-03

10.3969/j.issn.1672-4623.2015.05.023

李朋澤，碩士，主要從事農(nóng)業(yè)遙感和比輻射率遙感反演方面的研究。

2014-08-28。

項(xiàng)目來源：全球變化研究國家重大科學(xué)研究計(jì)劃資助項(xiàng)目（2010CB950902）。