基于物候特征參數(shù)的山西煤礦區(qū)典型復(fù)墾植被分類

張彥彬， 安楠， 劉佩艷， 賈坤， 姚云軍

(1.山西省自動化研究所,太原 030012； 2.堪薩斯州立大學(xué)農(nóng)學(xué)系環(huán)境與農(nóng)業(yè)空間分析實驗室,

堪薩斯州 66506,美國； 3.北京師范大學(xué)地理學(xué)與遙感科學(xué)學(xué)院,北京 100875)

基于物候特征參數(shù)的山西煤礦區(qū)典型復(fù)墾植被分類

張彥彬1， 安楠2， 劉佩艷1， 賈坤3， 姚云軍3

(1.山西省自動化研究所,太原 030012； 2.堪薩斯州立大學(xué)農(nóng)學(xué)系環(huán)境與農(nóng)業(yè)空間分析實驗室,

堪薩斯州 66506,美國； 3.北京師范大學(xué)地理學(xué)與遙感科學(xué)學(xué)院,北京 100875)

基于2001―2013年獲取的MOD13Q1 NDVI數(shù)據(jù)，采用低通平滑Savitzky-Golay(S-G)濾波方法、插值法及切比雪夫多項式(Chebyshev Polynomial)擬合對NDVI時序數(shù)據(jù)進(jìn)行重構(gòu)； 通過提取植被生長季開始日期、生長季長度、生長季結(jié)束日期、生長季NDVI最大值及NDVI最大值出現(xiàn)日期等關(guān)鍵物候特征參數(shù)，對研究區(qū)典型復(fù)墾植被類型進(jìn)行分類。結(jié)果表明： 研究區(qū)不同植被的物候特征具有顯著差異，從生長季開始日期及NDVI最大值出現(xiàn)日期來看，農(nóng)作物較有規(guī)律； 而林地的生長季NDVI累積總值則明顯區(qū)別于農(nóng)作物及草地； 農(nóng)作物、草地和林地基于植被物候特征參數(shù)分類取得了較好結(jié)果，總體分類精度達(dá)到89.67%，優(yōu)于采用多時相非監(jiān)督分類的結(jié)果； 該研究為山西省煤炭礦區(qū)生態(tài)環(huán)境恢復(fù)評價提供了一定的數(shù)據(jù)基礎(chǔ)。

遙感數(shù)據(jù)； MODIS； NDVI； 物候參數(shù)； 植被分類

0 引言

及時對煤礦復(fù)墾狀況進(jìn)行動態(tài)監(jiān)測及評估是煤礦復(fù)墾的重要需求。遙感數(shù)據(jù)分類是獲取區(qū)域土地覆蓋信息的重要手段，為土地覆蓋相關(guān)研究和應(yīng)用提供了不可或缺的重要數(shù)據(jù)基礎(chǔ)[1-2]。近年來，由于不同植被類型的生長特征存在明顯的差異，國內(nèi)外學(xué)者廣泛采用時序特征進(jìn)行植被覆蓋分類，提高了分類精度[3]，如Pringle等[4-7]利用中分辨率成像光譜儀歸一化植被指數(shù)(MODIS NDVI)數(shù)據(jù)對土地覆蓋分類進(jìn)行了研究，并且針對不同類型植被特征信息的提取，采用了不同的研究方法，如傅里葉變換法[8]、指數(shù)模型法[9]和閾值法[10]等。但上述方法均缺乏對植被物候特征及生理特性的分析。物候包括植物的發(fā)芽、開花、結(jié)果和落葉等[11]。通過時序植被指數(shù)的變化規(guī)律分析植被物候特征參數(shù)[12]，提取對應(yīng)的關(guān)鍵時間節(jié)點及特征值，如生長季開始日期、生長季結(jié)束日期、生長季長度、生長季NDVI最大值、生長季NDVI最大值對應(yīng)日期及生長季NDVI累計總值等物候參數(shù)[13]。國內(nèi)外利用時序遙感數(shù)據(jù)進(jìn)行大尺度的物候特征研究取得了很大進(jìn)展，如Beck等[14]利用函數(shù)擬合方法研究植被物候特征變化； 宮攀[15]利用MODIS數(shù)據(jù)關(guān)鍵物候特征參數(shù)對東北地區(qū)植被覆蓋進(jìn)行分類研究； Murakami等[16]和趙延征等[17]利用SPOT_VGT數(shù)據(jù)提取物候信息，對大面積農(nóng)作物空間分布進(jìn)行研究。這些研究表明，同一區(qū)域內(nèi)不同植被的物候特征具有明顯的差異，因此利用物候特征參數(shù)對土地覆蓋分類是可行的。但是，將物候特征參數(shù)應(yīng)用于煤礦復(fù)墾植被分類方面的研究相對較少。本文通過對物候參數(shù)的提取，充分利用不同植被類型的物候特征實現(xiàn)山西省6大煤田區(qū)高精度的典型復(fù)墾植被分類，為進(jìn)一步開展我國礦區(qū)生態(tài)環(huán)境評價提供支撐數(shù)據(jù)。

1 研究區(qū)概況與數(shù)據(jù)源

1.1 研究區(qū)概況

山西省地處黃河中游東岸、黃土高原東部，地理位置在E110°14′～114°33′，N34°34′～40°43′之間，總面積約15.6萬km2，占全國總面積的1.63%。按照地質(zhì)構(gòu)造、地理位置和規(guī)模大小，全省煤礦區(qū)劃分為大同、寧武、河?xùn)|、太原西山、霍西和沁水6大煤田(圖1)。研究區(qū)屬溫帶大陸性氣候，冬寒夏暖，4季分明，氣溫和降水的空間分布差異較大，植被類型多樣且分布差異大。

圖1 山西省6大煤田分布圖

由于煤炭資源大量開采破壞了煤礦區(qū)的生態(tài)環(huán)境，近年來山西省各級政府積極倡導(dǎo)煤礦區(qū)復(fù)墾，以治理和改善其脆弱的生態(tài)環(huán)境。研究該區(qū)域典型復(fù)墾植被類型對評價生態(tài)環(huán)境改善狀況以及政府部門制定相應(yīng)政策都具有重要意義。

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

本文使用的遙感數(shù)據(jù)來源于美國地球資源觀測系統(tǒng)(EROS)數(shù)據(jù)中心提供的時間跨度為2001―2013年、16 d合成、空間分辨率為250 m的MOD13Q1 NDVI影像。

對多時相MODIS NDVI時序數(shù)據(jù)進(jìn)行平滑處理是提取植被物候特征參數(shù)的基礎(chǔ)，可最大限度地減弱噪聲干擾。目前基于濾波函數(shù)的NDVI時序曲線平滑法主要有Savitzky-Golay(S-G)濾波、均值迭代濾波、傅里葉變換和時間序列諧波分析等。邊金虎等[18]曾利用S-G濾波法對MODIS NDVI時序數(shù)據(jù)進(jìn)行了重構(gòu)，并與傅立葉變換和均值迭代濾波相比較，結(jié)果表明采用S-G濾波法的擬合效果較好。S-G濾波算法保留了數(shù)據(jù)最大值和最小值的特性，并能很好地反映時序數(shù)據(jù)的局部突變信息及變化趨勢，因此本文采用S-G濾波法對時間序列NDVI進(jìn)行平滑處理。S-G濾波的基本公式為

(1)

式中：Y為NDVI原始值；Y*為NDVI擬合值；j為原始NDVI數(shù)組的系數(shù)；Ci為第i個NDVI值濾波時的系數(shù)；m為濾波窗口的寬度；N為滑動窗口所包括的數(shù)據(jù)點，N=2m+1。

2 研究方法

2.1 物候特征參數(shù)提取

為了更準(zhǔn)確地提取物候參數(shù)，本文首先對S-G濾波算法平滑后的曲線進(jìn)行空間插值，然后利用切比雪夫多項式(Chebyshev Polynomial)對曲線進(jìn)行擬合。利用切比雪夫插值法對插值節(jié)點尋優(yōu)，進(jìn)行函數(shù)最佳逼近，可有效地提高參數(shù)提取精度。

設(shè)x0,x1,…,xn,x為區(qū)間[a,b]上(n+1)個互不相同的點，f(x)∈Cn+1[a,b]，則對任何x∈[a,b]存在εx∈[x0,x1,…,xn,x]，使得拉格朗日插值余項Rn(x)=f(x)-Ln(x)滿足

(2)

其中

[x0,x1,…,xn,x]=[min{x0,x1,…,xn,x}, max{x0,x1,…,xn,x}]∈[a,b],

(3)

(4)

要使拉格朗日插值多項式Ln(x)盡量逼近f(x)，就要使余項Rn(x)盡量小。在Rn(x)中，f(x)是固定的，而εx是未知數(shù)； 所以要減小Rn(x)，只有恰當(dāng)選擇節(jié)點集，使得在插值區(qū)間內(nèi)余項的最大值為極小值。為了應(yīng)用切比雪夫多項式，首先應(yīng)將插值區(qū)間[a,b]通過簡單變換歸一化到區(qū)間[-1,1]，做變換

(5)

所以插值節(jié)點應(yīng)取為

(6)

式中：k為切比雪夫節(jié)點序號，k=0,1,2,…,n-1；a和b分別為插值區(qū)間的下限和上限。

所以需討論區(qū)間[-1,1]上函數(shù)的切比雪夫點Xk，即

(7)

則有

(8)

(9)

故切比雪夫插值法可使余項的最大值極小化，得到較佳逼近的多項式。2種不同類型植被NDVI曲線經(jīng)過S-G濾波及切比雪夫插值擬合的結(jié)果見圖2。

圖2 NDVI時間序列曲線重構(gòu)

本文選取生長季開始日期、生長季結(jié)束日期、生長季長度、生長季峰值、NDVI最大值對應(yīng)日期和生長季NDVI累計值等參數(shù)為分類物候特征參量(圖3)。

A： 生長季開始日期；B： 生長季結(jié)束日期；C： 生長季長度；D： 生長季峰值；E： NDVI最大值對應(yīng)日期；F： 生長季NDVI累計值(灰色部分)

圖3 植被物候特征定義

Fig.3 Definitions of vegetation phonological parameters

生長季開始日期定義為擬合曲線上升速率急劇增加的點A所對應(yīng)的日期。對擬合曲線計算一階導(dǎo)數(shù)，一階導(dǎo)數(shù)曲線出現(xiàn)的第一個峰值表示原擬合曲線上升率最大的點，這個點所對應(yīng)的日期即為生長季開始日期A′?？紤]到MODIS 16 d合成數(shù)據(jù)不利于研究植被生長季開始日期，本文利用空間插值法模擬植被真實發(fā)芽日期。同理，可得到生長季結(jié)束日期B′。由生長季開始日期和生長季結(jié)束日期可得到生長季長度C。通過監(jiān)測擬合曲線波峰，比較波峰處的NDVI值，提取NDVI最大值D，D點所對應(yīng)的日期即為NDVI最大值出現(xiàn)日期E。F為植被在1 a中生長期內(nèi)的生長量，即通過計算得到的植被生長期內(nèi)的NDVI累計總值。從圖3可以看出，擬合曲線與坐標(biāo)橫軸之間的面積(灰色部分)為NDVI累計總值。根據(jù)以上物候期關(guān)鍵時間節(jié)點及特征值，通過梯形公式(10)積分計算切比雪夫多項式曲線下方的面積得出生長季累計總值F。在利用梯形公式進(jìn)行積分計算時，假設(shè)積分區(qū)間[a,b]平均分成N份，即

(10)

式中：a和b分別為積分區(qū)間的下限和上限；N為區(qū)間被分割數(shù)；K為分點序號，K=1,2,…,n-1。

2.2 典型復(fù)墾植被分類

本文利用MODIS的1 a 23個時相的NDVI時序數(shù)據(jù)來描述NDVI年內(nèi)動態(tài)變化特征，該動態(tài)變化曲線可反映植被的整個生長過程(即從返青/發(fā)芽到枯萎/收割的周期性規(guī)律)。利用不同植被類型具有的不同生長規(guī)律特征，可以實現(xiàn)典型植被的分類。典型復(fù)墾植被類型包括耕地、林地和草地，因此本文分類的類型主要包括城鎮(zhèn)及非植被覆蓋區(qū)、耕地、林地和草地。

2.2.1 城鎮(zhèn)及非植被覆蓋信息提取

由于城鎮(zhèn)中建筑物和植被交錯分布，因此250 m空間分辨率的影像數(shù)據(jù)大多數(shù)屬于混合像元。建筑物的NDVI值較低且隨季節(jié)變化小，即城鎮(zhèn)NDVI值趨于平緩，沒有明顯峰值，且NDVI<0.4，可根據(jù)此特征提取城鎮(zhèn)及非植被覆蓋信息。

2.2.2 耕地信息提取

因為山西省農(nóng)作物種植模式只有1年1熟制和1年2熟制。如果1 a有2個波峰，則可以判定植被類型為耕地且為1年2季作物。本文采用二次差分法提取NDVI時序曲線極大值，進(jìn)而計算波峰個數(shù)，判斷是否為1年2熟制耕地，二次差分法表達(dá)式為

(11)

式中：S1為前后2時相NDVI的差值序列，若S1>0，則把S2賦值為1； 若S1<0，則把S2賦值為-1,得到序列S2；S3為S2中前后2元素差值序列，當(dāng)S3為-2時，表示為時序數(shù)據(jù)的極大值(即該時相對應(yīng)的NDVI為波峰)；S3=2時，表示為時序數(shù)據(jù)的極小值(即NDVI在該位置出現(xiàn)波谷)； 可依據(jù)S3提取像元波峰出現(xiàn)的時間及NDVI值。

由于二次差分法對離散點的峰值特別敏感，會將每一個微小的“峰”提取出來，因此需要利用一定的約束條件對干擾波峰進(jìn)行取舍。范錦龍等[19]采用閾值法將裸地覆被形成的微小的“峰”舍去，同時控制時序長度來屏蔽生長季之外的“峰”和“冬前峰”。根據(jù)已有研究結(jié)果，判斷植被波峰的NDVI值不小于0.4且生長季長度應(yīng)達(dá)到90 d。如果1 a有2個波峰，則可判定植被類型為耕地且為1年2季農(nóng)作物種植模式； 如果1 a有1個波峰，則要根據(jù)植被的物候期特點做進(jìn)一步的判斷。

研究區(qū)1年1季春季作物的返青/出芽日期最早出現(xiàn)在2月下旬―3月上旬(如冬小麥)，NDVI最大值出現(xiàn)在4月下旬―5月上旬； 秋季作物NDVI最大值則出現(xiàn)在10月下旬―11月份上旬。根據(jù)研究區(qū)農(nóng)作物物候特征參數(shù)，可將耕地與其他植被類型明顯區(qū)分。

2.2.3 林地信息提取

根據(jù)林地的物候期特點，在其整個生長季中，NDVI曲線波峰平緩且持續(xù)時間最長、振幅最大，即生長季NDVI累計值最大。通過對2004―2013年10 a數(shù)據(jù)進(jìn)行反復(fù)試驗證明，生長季NDVI累計值在7 000處能與其他植被類型明顯區(qū)分。

2.2.4 草地信息提取

參照NDVI最大值分布特征，草地在1 a之內(nèi)0.4≤NDVI最大值≤0.6，并結(jié)合生長季NDVI累計值，草地的NDVI累計值應(yīng)該大于城鎮(zhèn)及非植被覆蓋地的NDVI累計值而小于農(nóng)作物的NDVI累計值。據(jù)此，可將草地與其他植被類型區(qū)分出來。

2.2.5 多時相NDVI數(shù)據(jù)非監(jiān)督分類

為了比較本文方法的分類效果，選取2013年獲取的MODIS13Q1 NDVI時間序列數(shù)據(jù)進(jìn)行非監(jiān)督分類，并進(jìn)行精度評價。對2013年23個時相數(shù)據(jù)進(jìn)行集成，利用ISODATA算法進(jìn)行分類，輸入?yún)?shù)分別設(shè)置為： 初始分類為50類，最大迭代次數(shù)20次，循環(huán)收斂閾值為0.998，生成研究區(qū)的土地覆蓋分類圖。

2.3 分類精度評價

精度評價是遙感數(shù)據(jù)分類的重要步驟，可定量地對分類精度進(jìn)行分析。本文采用Congalton[20]提出的誤差矩陣(亦稱混淆矩陣)方法衡量分類精度。混淆矩陣的主要評價指標(biāo)包括總體精度、用戶精度及Kappa系數(shù)等。其中總體精度是所有分類正確的樣本數(shù)占樣本總數(shù)的百分比，反映分類結(jié)果總的正確程度； Kappa系數(shù)是反映整個誤差矩陣的精度系數(shù)，其計算公式為

(12)

式中：K為Kappa 系數(shù)；r為誤差矩陣的列數(shù)；Xii為誤差矩陣第i行i列(主對角線)上的值；Xi+與X+i分別為分類誤差矩陣第i行像元數(shù)的和與第i列像元數(shù)的和；N為檢驗樣本總數(shù)。

3 結(jié)果與分析

3.1 研究區(qū)物候分布特征

為了獲取研究區(qū)植被主要物候期空間分布特征，本文主要計算了2001―2013年植被物候期特征均值，提取了研究區(qū)植被關(guān)鍵物候特征參數(shù)(圖4)。

通過觀察可以看出，不同植被表現(xiàn)出不同的物候特征。對于生長季開始日期，本文研究植被可能出現(xiàn)發(fā)芽生長的時間周期，即從第45—300 d。

從圖4(a)可以看出，植被生長季開始日期較早的區(qū)域集中分布在霍西煤田南部的冬小麥種植區(qū)； 寧武煤田中部、河?xùn)|煤田北部主要是1年1熟模式種植區(qū)，生長季開始日期較晚； 沁水煤田西南部、河?xùn)|煤田東南部及霍西煤田西南部主要為山地，土地覆蓋類型多為林地，生長季開始日期也較晚。

圖4(b)為植被在1 a之中NDVI最大值的分布情況。沁水煤田西南部、河?xùn)|煤田東南部及霍西煤田西南部區(qū)域林地覆蓋密集，NDVI值最大； 霍西煤田南部、沁水煤田東南部耕地分布密集，1 a中NDVI最大值次之； 在大同煤田北部、河?xùn)|煤田中部、霍西煤田西北部、寧武煤田南部植被覆蓋稀疏，主要土地覆蓋類型為草地，NDVI最大值較?。?而寧武煤田北部分布有平朔煤炭礦區(qū)(露天煤礦)及零星分布在各煤田的居民點，在1 a中的NDVI最大值最小。

(a) 生長季開始日期分布 (b) NDVI最大值分布

(c) NDVI最大值出現(xiàn)日期分布 (d) 生長季NDVI累計總值分布

圖4 植被物候參數(shù)分布

Fig.4 Distribution of vegetation phonological parameters

根據(jù)MOD13Q1 16 d合成數(shù)據(jù)的特點，本文將1 a中植被覆蓋數(shù)據(jù)按每16 d為1個周期進(jìn)行分割，對這23個周期的NDVI最大值出現(xiàn)日期進(jìn)行分析。根據(jù)山西省植被生長規(guī)律及農(nóng)作物種植結(jié)構(gòu)，本文僅研究第1—18周期NDVI最大值出現(xiàn)日期分布情況。從圖4(c)看出，NDVI最大值出現(xiàn)在第7―8周期(4月上旬或下旬)的區(qū)域主要集中分布在霍西煤田南部，該區(qū)域的土地覆蓋類型為1年2熟的春季作物。NDVI最大值出現(xiàn)在第16周期(8月下旬)的區(qū)域主要分布在霍西煤田南部，土地覆蓋類型以1年2熟的夏季作物為主。NDVI最大值出現(xiàn)在第13―14周期(7月上旬或下旬)的區(qū)域密集集中在沁水煤田西南部、河?xùn)|煤田東南部及霍西煤田西南部，主要土地覆蓋類型為林地。

從圖4(d)可以看出，NDVI累計總值大的區(qū)域(即紅色區(qū)域)主要集中分布在霍西煤田的東部、河?xùn)|煤田的南部及呂梁山的南端以及西山煤田的中部； 霍西煤田南部、沁水煤田東南部的NDVI累計總值次之； 而河?xùn)|煤田北部、寧武煤田北部及大同煤田的NDVI累計總值較小。

3.2 典型復(fù)墾植被分類結(jié)果與驗證

采用植被物候特征參數(shù)的分類結(jié)果和多時相非監(jiān)督分類結(jié)果如圖5所示。

從基于植被物候特征參數(shù)的土地覆蓋分類結(jié)果(圖5(a))來看，土地覆蓋類型的主要分布與實際情況相符。林地主要分布在霍西煤田的東部(即太行山西南邊緣)、河?xùn)|煤田的南部、呂梁山的南端及西山煤田的中部(即中條山中部)； 耕地分布在霍西煤田南部(即晉南盆地)、沁水煤田的東南部分布密集； 草地主要分布于河?xùn)|煤田的山壑、丘陵及河流兩岸。從基于多時相非監(jiān)督分類的分類結(jié)果(圖5(b))的空間分布來看，煤礦區(qū)部分區(qū)域出現(xiàn)成片土地覆蓋混合為1種類別，主要分布在河?xùn)|煤田北部，與實際情況出入較大； 沁水煤田西南部、河?xùn)|煤田東南部及霍西煤田西南部區(qū)域林地覆蓋密集，耕地主要分布在霍西煤田南部，與實際情況基本吻合，但分類效果不如基于物候特征參數(shù)的分類結(jié)果。

(a) 基于植被物候特征參數(shù)(b) 基于多時相非監(jiān)督分類

圖5 土地覆蓋分類結(jié)果

Fig.5 Results of land cover classification

為了定量分析2種分類方法的分類精度，從上述4種土地覆蓋類型中分別隨機(jī)選取200個樣本點對分類結(jié)果進(jìn)行驗證，分類混淆矩陣如表1和表2所示。本文方法總體分類精度為89.67%，Kappa系數(shù)為0.82； 而多時相非監(jiān)督分類總體精度為78%，Kappa系數(shù)為0.74。

表1 物候特征參數(shù)分類混淆矩陣

表2 多時相非監(jiān)督分類混淆矩陣

從精度評價結(jié)果可以看出，參照植被物候特征參數(shù)取得的分類結(jié)果中，除林地的生產(chǎn)者精度與非監(jiān)督分類結(jié)果持平外，其他3種植被類型的精度均較非監(jiān)督分類結(jié)果有明顯提高。物候特征參數(shù)分類中，與其他植被類型相比，農(nóng)作物的返青/發(fā)芽日期、NDVI最大值出現(xiàn)日期等物候期關(guān)鍵參數(shù)具有較為明顯的差異，所以利用物候特征參數(shù)的分類方法對耕地的辨識精度很高。在整個生長季，林地在生長季長度及NDVI波峰持續(xù)時間上與其他植被明顯不同，并且由這些關(guān)鍵節(jié)點所確定的生長季NDVI累計總值遠(yuǎn)遠(yuǎn)大于其他植被類型。草地由于在物候特征參數(shù)上與林地、農(nóng)作物特征值分布比較相近，且農(nóng)作物的物候特征參數(shù)值方差較大，部分農(nóng)作物和林地被錯分為草地，導(dǎo)致草地的分類精度略低； 而城鎮(zhèn)及其他非植被覆蓋地類的物候參數(shù)特征并不明顯，且城鎮(zhèn)內(nèi)部建筑物和植被交錯分布，對于250 m空間分辨率的MODIS NDVI影像數(shù)據(jù)來說，均作為混合像元存在，這也成為影響分類精度的一個重要因素。多時相非監(jiān)督分類結(jié)果雖然比單一時相分類精度有了較大提高，但由于“同物異譜”或“異物同譜”及波段冗余等問題的出現(xiàn)，分類結(jié)果精度仍不盡人意，最明顯的是將耕地錯分為無植被覆蓋或者草地，將灌木林地錯分為草地。分類精度評價結(jié)果表明，基于物候特征參數(shù)的分類結(jié)果在耕地和草地類別上明顯優(yōu)于非監(jiān)督分類結(jié)果，說明基于物候特征參數(shù)的分類方法能夠?qū)崿F(xiàn)山西煤礦區(qū)典型復(fù)墾植被類型較好識別效果。

煤礦區(qū)典型復(fù)墾植被類型分類結(jié)果為煤礦區(qū)復(fù)墾的生態(tài)恢復(fù)評估提供了重要的數(shù)據(jù)基礎(chǔ)。本文研究結(jié)果結(jié)合煤炭礦區(qū)植被的動態(tài)監(jiān)測，對評估煤礦區(qū)復(fù)墾狀況(包括煤礦區(qū)有無進(jìn)行復(fù)墾、復(fù)墾處于哪個階段及復(fù)墾植被種類等)具有重要意義。

4 結(jié)論

本文以2001―2013年獲取的MOD13Q1 NDVI為主要數(shù)據(jù)源，利用S-G濾波算法對時間序列NDVI曲線進(jìn)行平滑，并通過插值及切比雪夫多項式進(jìn)行曲線擬合； 通過提取研究區(qū)物候特征參數(shù)實現(xiàn)對山西省6大煤田復(fù)墾典型植被類型的分類。主要結(jié)論如下：

1)利用切比雪夫多項式對S-G濾波平滑后的曲線進(jìn)行插值及擬合，達(dá)到了更為精確提取植被物候特征參數(shù)的效果。

2)依據(jù)物候特征參數(shù)對研究區(qū)主要復(fù)墾典型植被類型的總體分類精度達(dá)到了89.67%，明顯優(yōu)于利用多時相NDVI數(shù)據(jù)非監(jiān)督分類。

本文尚未考慮混合像元的影響及非遙感數(shù)據(jù)(如氣象數(shù)據(jù))對植被物候參數(shù)的影響。在今后的研究中應(yīng)進(jìn)行混合像元分解研究，并將遙感數(shù)據(jù)與非遙感數(shù)據(jù)(如氣溫、降水等因素)相結(jié)合，進(jìn)行更為精細(xì)的分類。

[1] 賈坤,李強(qiáng)子,田亦陳,等.遙感影像分類方法研究進(jìn)展[J].光譜學(xué)與光譜分析,2011,31(10):2618-2623. Jia K,Li Q Z,Tian Y C,et al.A review of classification methods of remote sensing imagery[J].Spectroscopy and Spectral Analysis,2011,31(10):2618-2623.

[2] 賈坤,李強(qiáng)子.農(nóng)作物遙感分類特征變量選擇研究現(xiàn)狀與展望[J].資源科學(xué),2013,35(12):2507-2516. Jia K,Li Q Z.Review of features selection in crop classification using remote sensing data[J].Resources Science,2013,35(12):2507-2516.

[3] Jia K,Liang S L,Wei X Q,et al.Land cover classification of landsat data with phenological features extracted from time series MODIS NDVI data[J].Remote Sensing,2014,6(11):11518-11532.

[4] 郭芬芬,范建容,邊金虎,等.基于MODIS NDVI時間序列數(shù)據(jù)的藏北草地類型識別[J].遙感技術(shù)與應(yīng)用,2011,26(6):821-826. Guo F F,Fan J R,Bian J H,et al.Grassland types identification based on time-series MODIS NDVI data in northern Tibet[J].Remote Sensing Technology and Application,2011,26(6):821-826.

[5] 宮攀,陳仲新,唐華俊,等.基于MODIS溫度/植被指數(shù)的東北地區(qū)土地覆蓋分類[J].農(nóng)業(yè)工程學(xué)報,2006,22(9):94-99. Gong P,Chen Z X,Tang H J,et al.Land cover classification based on MODIS temperature-vegetation index time-series data in northeastern China[J].Transactions of the CSAE,2006,22(9):94-99.

[6] Pringle M J,Denham R J,Devadas R.Identification of cropping activity in central and southern Queensland,Australia,with the aid of MODIS MOD13Q1 imagery[J].International Journal of Applied Earth Observation and Geoinformation,2012,19:276-285.

[7] Mkhabela M S,Bullock P,Raj S,et al.Crop yield forecasting on the Canadian Prairies using MODIS NDVI data[J].Agricultural and Forest Meteorology,2011,151(3):385-393.

[8] 那曉東,張樹清,李曉峰,等.MODIS NDVI時間序列在三江平原濕地植被信息提取中的應(yīng)用[J].濕地科學(xué),2007,5(3):227-236. Na X D,Zhang S Q,Li X F,et al.Application of MODIS NDVI time series to extracting wetland vegetation information in the Sanjiang Plain[J].Wetland Science,2007,5(3):227-236.

[9] 潘耀忠,李樂,張錦水,等.基于典型物候特征的MODIS-EVI時間序列數(shù)據(jù)農(nóng)作物種植面積提取方法——小區(qū)域冬小麥實驗研究[J].遙感學(xué)報,2011,15(3):578-594. Pan Y Z,Li L,Zhang J S,et al.Crop area estimation based on MODIS-EVI time series according to distinct characteristics of key phenology phases:A case study of winter wheat area estimation in small-scale area[J].Journal of Remote Sensing,2011,15(3):578-594.

[10]郝鵬宇,牛錚,王力,等.基于歷史時序植被指數(shù)庫的多源數(shù)據(jù)作物面積自動提取方法[J].農(nóng)業(yè)工程學(xué)報,2012,28(23):123-131. Hao P Y,Niu Z,Wang L,et al.Multi-source automatic crop pattern mapping based on historical vegetation index profiles[J].Transactions of the Chinese Society of Agricultural Engineering,2012,28(23):123-131.

[11]康峻,侯學(xué)會,牛錚,等.基于擬合物候參數(shù)的植被遙感決策樹分類[J].農(nóng)業(yè)工程學(xué)報,2014,30(9):148-156. Kang J,Hou X H,Niu Z,et al.Decision tree classification based on fitted phenology parameters from remotely sensed vegetation data[J].Transactions of the Chinese Society of Agricultural Engineering,2014,30(9):148-156.

[12]夏傳福,李靜,柳欽火.植被物候遙感監(jiān)測研究進(jìn)展[J].遙感學(xué)報,2013,17(1):1-16. Xia C F,Li J,Liu Q H.Review of advances in vegetation phenology monitoring by remote sensing[J].Journal of Remote Sensing,2013,17(1):1-16.

[13]李正國,楊鵬,周清波,等.基于時序植被指數(shù)的華北地區(qū)作物物候期/種植制度的時空格局特征[J].生態(tài)學(xué)報,2009,29(11):6216-6226. Li Z G,Yang P,Zhou Q B,et al.Research on spatiotemporal pattern of crop phenological characteristics and cropping system in north China based on NDVI time series data[J].Acta Ecologica Sinica,2009,29(11):6216-6226.

[14]Beck P S A,Atzberger C,H?gda K A,et al.Improved monitoring of vegetation dynamics at very high latitudes:A new method using MODIS NDVI[J].Remote Sensing of Environment,2006,100(3):321-334.

[15]宮攀.基于MODIS數(shù)據(jù)關(guān)鍵物候特征參數(shù)的東北地區(qū)植被覆蓋分類[J].資源科學(xué),2010,32(6):1154-1160. Gong P.Vegetation classification based on phenology indices derived from MODIS data in northeastern China[J].Resources Science,2010,32(6):1154-1160.

[16]Murakami T,Ogawa S,Ishitsuka N,et al.Crop discrimination with multitemporal SPOT/HRV data in the Saga Plains,Japan[J].International Journal of Remote Sensing,2001,22(7):1335-1348.

[17]楊延征,趙鵬翔,郝紅科,等.基于SPOT-VGT NDVI的陜北植被覆蓋時空變化[J].應(yīng)用生態(tài)學(xué)報,2012,23(7):1897-1903. Yang Y Z,Zhao P X,Hao H K,et al.Spatiotemporal variation of vegetation in northern Shaanxi of northwest China based on SPOT-VGT NDVI[J].Chinese Journal of Applied Ecology,2012,23(7):1897-1903.

[18]邊金虎,李愛農(nóng),宋孟強(qiáng),等.MODIS植被指數(shù)時間序列Savitzky-Golay濾波算法重構(gòu)[J].遙感學(xué)報,2010,14(4):725-741. Bian J H,Li A N,Song M Q,et al.Reconstruction of NDVI time-series datasets of MODIS based on Savitzky-Golay filter[J].Journal of Remote Sensing,2010,14(4):725-741.

[19]范錦龍,吳炳方.復(fù)種指數(shù)遙感監(jiān)測方法[J].遙感學(xué)報,2004,8(6):628-636. Fan J L,Wu B F.A methodology for retrieving cropping index from NDVI profile[J].Journal of Remote Sensing,2004,8(6):628-636.

[20]Congalton R G.A review of assessing the accuracy of classifications of remotely sensed data[J].Remote Sensing of Environment,1991,37(1):35-46.

(責(zé)任編輯： 邢宇)

Typical reclamation vegetation classification based on phenological feature parameters for coalfields in Shanxi Province

ZHANG Yanbin1， AN Nan2， LIU Peiyan1， JIA Kun3， YAO Yunjun3

(1.ShanxiAutomationResearchInstitute，Taiyuan030012，China； 2.Ecology&AgricultureSpatialAnalysisLaboratory，DepartmentofAgronomy，KansasStateUniversity，Kansas66506,USA； 3.SchoolofGeography，BeijingNormalUniversity，Beijing100875，China)

In this paper, the authors reconstructed MOD13Q1 time-series NDVI data from 2001 to 2013 using Savitzky-Golay filter and Chebyshev Polynomial methods for classifying vegetation types in the six coalfields in Shanxi Province. The key phenological parameters were extracted from the reconstructed NDVI data, such as the beginning dates of the growing season, length of the growing season, the ending dates of the growing season, the maximum NDVI value and the responding dates. The results show that different vegetation types of the six major coalfields in Shanxi have different phenological features. Cropland has distinguishable differences from grass and forest. Similarly, forest is distinguished from grass and cropland by integration of total growth. It is shown that the classification of vegetation types can achieve better results by extracting and analyzing the phonological parameters compared with multi-temporal unsupervised classification. The overall classification accuracy reaches 89.67%. This study provides a robust method for assessing long-term ecological conditions and monitoring vegetation coverage changes of the six major coalfields in Shanxi Province.

remote sensing data； MODIS； NDVI； phenological parameters； vegetation classification

10.6046/gtzyyg.2017.01.26

張彥彬,安楠,劉佩艷,等.基于物候特征參數(shù)的山西煤礦區(qū)典型復(fù)墾植被分類[J].國土資源遙感,2017,29(1):170-177.(Zhang Y B,An N,Liu P Y,et al.Typical reclamation vegetation classification based on phenological feature parameters for coalfields in Shanxi Province[J].Remote Sensing for Land and Resources,2017,29(1):170-177.)

2015-09-29；

2015-12-18

國際合作項目“利用衛(wèi)星遙感技術(shù)對煤礦復(fù)墾生態(tài)環(huán)境的動態(tài)監(jiān)測及分析”(編號： 2013DFA91870)和中國科學(xué)院數(shù)字地球重點實驗室開放基金項目“低空間分辨率遙感數(shù)據(jù)時相特征改善高分辨率數(shù)據(jù)農(nóng)作物分類精度研究” (編號： 2014LDE011)共同資助。

張彥彬(1967-),男,高級工程師,主要從事遙感應(yīng)用方面的研究。Email： zyb9633@163.com。

賈坤(1983-),男,博士,副教授,主要從事定量遙感和土地覆蓋分類方面的研究。 Email： jiakun@bnu.edu.cn。

TP 79; S 127

A

1001-070X(2017)01-0170-08