對(duì)于MOD12Q1 土地覆蓋分類數(shù)據(jù)，本文選擇國(guó)際地圈生物圈計(jì)劃( IGBP) 全球植被分類方案，數(shù)據(jù)處理主要包括投影轉(zhuǎn)換和土地利用類型合并。文中參照中國(guó)土地資源分類系統(tǒng)對(duì) IGBP 分類體系進(jìn)行類別合并( 表 4) ，最終分為水體、林地、草地、耕地、城市用地及未利用地 6 類，以此提取林地信息。MOD09A1 地表反射率數(shù)據(jù)的處理主要為去云處理，文中對(duì)每年獲取的 4 景 MOD09A1 數(shù)據(jù)利用最小值合成法合成一幅無(wú)云影像，計(jì)算歸一化建筑指數(shù)( NDBI) ［31］，突出影像中的城市用地信息; 城市信息提取，文中考慮到城市用地面積相對(duì)較小，具有整體集中、局部分散的特點(diǎn)，分布散落的城市在遙感影像中受周?chē)匚锏挠绊戄^大，在分類過(guò)程中容易出現(xiàn)錯(cuò)分現(xiàn)象，導(dǎo)致分類精度較低。因此，為保證城市用地的提取精度，本文采用遙感目視解譯與參考研究區(qū) NDBI、高分辨率遙感影像等輔助數(shù)據(jù)相結(jié)合的手段提取城市用地信息。NDBI 的計(jì)算公式如下所示:

式中，ＲSWIＲ和 ＲNIＲ表示短波紅外和近紅外波段的光譜反射率，分別對(duì)應(yīng) MOD09A1 數(shù)據(jù)的第 6 和第 2 波段。

土地分類結(jié)果更新，文中將遙感目視解譯提取的城市用地替代 MOD12Q1 中的城市用地，同時(shí)更新分類結(jié)果，得到山東省 3 個(gè)典型城市 2001—2010 年土地利用分類圖( 圖 3) ，由圖 3 可以看出，3 個(gè)城市主要以耕地為主，所占面積達(dá)到 80% 以上，分布較為穩(wěn)定; 水體、草地、未利用地所占面積較小，年際變化不明顯，而林地與城市用地的變化相對(duì)較為明顯; 信息統(tǒng)計(jì)，文中分別統(tǒng)計(jì)青島、濟(jì)南、淄博 10 年間林地和城市用地的面積，進(jìn)行百分比統(tǒng)計(jì)，表 5 為統(tǒng)計(jì)結(jié)果。

聲明：參考《生態(tài)學(xué)報(bào)》第30卷16期，如涉版權(quán)請(qǐng)聯(lián)系刪除

相關(guān)鏈接：土地覆蓋，植被，林地