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

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

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

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

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