李陽(yáng)陽(yáng)，張 軍，劉陳立，楊旭超，李 杰

(云南大學(xué)資源環(huán)境與地球科學(xué)學(xué)院，云南 昆明 650504)

老撾北部5省橡膠林提取及時(shí)空擴(kuò)張研究

李陽(yáng)陽(yáng)，張 軍*，劉陳立，楊旭超，李 杰

(云南大學(xué)資源環(huán)境與地球科學(xué)學(xué)院，云南 昆明 650504)

目的分析老撾北部“金三角”地區(qū)橡膠林提取方法及橡膠林時(shí)空擴(kuò)張?zhí)卣?，為合理?guī)劃橡膠種植、區(qū)域生態(tài)環(huán)境保護(hù)等提供方法和依據(jù)。方法基于MODIS影像和橡膠林物候特征，選取不同地類(lèi)特征參數(shù)，實(shí)現(xiàn)橡膠林空間分布快速可重復(fù)的決策樹(shù)提取新方法，并基于單一土地利用率模型研究橡膠林?jǐn)U張機(jī)制。結(jié)果(1)橡膠林提取時(shí)間窗口為1月至3月，基于不同地類(lèi)特征參數(shù)生成的橡膠成林、幼林決策樹(shù)提取方法總精度為80.67%，高分一號(hào)影像通過(guò)面向?qū)ο蠓椒?yàn)證MODIS影像提取結(jié)果精度可靠。(2)2010—2015年，老撾北部5省橡膠林由12.6萬(wàn)hm2擴(kuò)張到16.96萬(wàn)hm2，面積變化率最大的是烏多姆賽省，豐沙里省擴(kuò)張速度最慢。結(jié)論MODIS影像提取橡膠林方法和結(jié)果簡(jiǎn)單有效。橡膠林分布在海拔和坡度上呈現(xiàn)不斷上升的趨勢(shì)，海拔1 500 m以上零星分布，坡度上分布上限為35°，而坡向擴(kuò)張無(wú)明顯規(guī)律。從LUCC視角分析得到橡膠林時(shí)空擴(kuò)張的主要影響因素是社會(huì)政策和經(jīng)濟(jì)因素。

多源遙感；老撾北部；橡膠林；時(shí)空擴(kuò)張

Abstract：[Objective] The aim of this research is to analyze the extraction method of rubber forest in the “Golden Triangle” region of northern Laos, to monitor the rubber forest and its spatial expansion dynamically, and to provide methods and lay a basis for rational planning of rubber planting and regional ecological environment protection. [Method]A new method for rapid and repeatable extraction of rubber forest spatial distribution was used based on MODIS remote sensing data, the phenological characteristics and the geometric parameters of rubber forest, and the expansion mechanism of rubber forest was analyzed based on single land use rate model. [Result] (1) The rubber tree extraction ranged from January to March. The precision of extraction method based on phenological characteristics and geometric parameters was 80.67%. GF-1 images were classified by the object-oriented method to verify the MODIS data which proved to be reliable. (2) The area of rubber forest in five provinces of northern Laos was expanded from 12.6×104hm2to 16.96 ×104hm2between 2010 and 2015. The area of rubber forest changed the most in Oudomxay province, and the least in Phongsaly province. [Conclusion] The method to extract rubber forest by MODIS data is simple, valid and effective. The distribution of rubber forests increases with the elevation and slope and scatters over the elevation of 1500m, the upper limits of slope is 35°. While, no obvious rule about distribution on the aspect is found. The main factors influencing the spatial and temporal expansion of rubber forest are social policy and economics from the perspective of LUCC.

Keywords: Multi-source remote sensing; Laos; rubber forest; spatial and temporal expansion

天然橡膠是國(guó)防和經(jīng)濟(jì)發(fā)展不可或缺的戰(zhàn)略物資[1]，在“金三角”貧窮國(guó)家發(fā)展替代罌粟的橡膠種植業(yè)[2-4]，也是我國(guó)促進(jìn)邊疆繁榮穩(wěn)定的國(guó)家戰(zhàn)略。地處廣義“金三角”老撾北部5省(南塔省、豐沙里省、烏多姆塞省、波喬省和瑯勃拉邦省)充足的水熱條件為該區(qū)橡膠種植提供了良好的環(huán)境。隨著橡膠林不斷擴(kuò)張，原生植被的破壞與不合理的橡膠種植將損害中老邊境地區(qū)生物多樣性與生態(tài)安全，因此，掌握橡膠林分布和時(shí)空擴(kuò)張監(jiān)測(cè)機(jī)制對(duì)橡膠林種植規(guī)劃、橡膠產(chǎn)品的宏觀調(diào)控和進(jìn)出口貿(mào)易等都具有重要作用。

橡膠林提取早期多以中低分辨率遙感影像通過(guò)監(jiān)督分類(lèi)進(jìn)行，分類(lèi)精度可達(dá)80%左右[5-7]。近年來(lái)，以中高分辨率遙感影像為主要數(shù)據(jù)源的面向?qū)ο蠓椒ㄖ饾u成熟，Ekadinata等[8]、劉少軍等[9]、劉曉娜等等[10]用面向?qū)ο蠓椒ㄌ崛〉南鹉z林總精度在85%左右。決策樹(shù)方法需要對(duì)比多時(shí)相數(shù)據(jù)分析橡膠林提取時(shí)間窗口，基于光譜特征和其他特征提高分類(lèi)精度[11]，目前尚處于探索階段。綜合而言，目前難以形成高效、穩(wěn)定、可靠的橡膠林提取算法，快速了解大尺度橡膠林分布面積及時(shí)空特征仍是瓶頸。面向?qū)ο蠓椒ㄟm用高分辨率影像，在保證提取精度的前提下較難適用于大尺度范圍，但可作為其他方法分類(lèi)的驗(yàn)證；而基于物候特征、像元和光譜特征、紋理特征等的優(yōu)化決策樹(shù)分類(lèi)方法受影像數(shù)據(jù)源限制較小，分類(lèi)精度較高[12]，所以，在大尺度范圍內(nèi)具有更好的可操作性和普適性。研究基于MODIS影像，根據(jù)老撾北部橡膠林物候特征構(gòu)建橡膠林的空間分布遙感提取算法，進(jìn)行橡膠林時(shí)空擴(kuò)張監(jiān)測(cè)研究。掌握老撾北部五省橡膠種植狀況，為老撾北部國(guó)土資源開(kāi)發(fā)、生態(tài)環(huán)境監(jiān)測(cè)、橡膠產(chǎn)業(yè)可持續(xù)發(fā)展及社會(huì)經(jīng)濟(jì)發(fā)展提供參考。

1 研究區(qū)概況

圖1 研究區(qū)位示意圖Fig.1 Location of study area

2 材料與方法

2.1數(shù)據(jù)源及預(yù)處理

遙感數(shù)據(jù)包括8天合成的2015全年和2010年1至3月的MODIS地表反射率產(chǎn)品MOD09A1、MOD09Q1，用于決策樹(shù)橡膠林提取方法的構(gòu)建；2015年3月2日的高分一號(hào)寬幅(GF-1WFV)、2016年1月16日高分辨率(GF-1PMS)影像用來(lái)驗(yàn)證MODIS影像的提取效果；谷歌高清影像輔助決策地類(lèi)樣本的選擇。此外，還包括來(lái)源于地理空間數(shù)據(jù)云的DEM和2次野外調(diào)查數(shù)據(jù)。

研究區(qū)涉及到MODIS影像圖幅分別為h27v06，h27v07，采用MRT軟件基于命令行批處理實(shí)現(xiàn)影像拼接裁剪和投影轉(zhuǎn)換。GF-1WFV、GF-1PMS影像景序列號(hào)分別為961759，2001118，利用野外調(diào)查控制點(diǎn)進(jìn)行正射校正，并在大氣校正后融合GF-1PMS影像生成2m多光譜遙感影像。

2.2研究方法

2.2.1橡膠林提取方法 老撾北部熱帶季風(fēng)氣候使得區(qū)域內(nèi)形成豐富的植被覆蓋地類(lèi)。本研究在國(guó)際生物圈地圈計(jì)劃(IGBP)分類(lèi)體系基礎(chǔ)上，結(jié)合研究區(qū)實(shí)際情況，將植被覆蓋地類(lèi)分為橡膠林地、有林地、灌木林地、旱地、水田，其中橡膠林地根據(jù)不同林齡的光譜差異，分為橡膠幼林(<6a)和橡膠成林(≥6a)。研究借助野外調(diào)查數(shù)據(jù)和谷歌高清影像建立5個(gè)省各地類(lèi)的遙感解譯標(biāo)志，共300個(gè)樣本。不同植被覆蓋地類(lèi)特征參數(shù)的選擇是構(gòu)建決策樹(shù)的關(guān)鍵，研究區(qū)橡膠林獨(dú)特的落葉現(xiàn)象易于在特定時(shí)間區(qū)別其他地類(lèi)，因此，研究基于物候特征、歸一化指數(shù)、主成分分析、紋理信息、原始光譜特征等采用決策樹(shù)分類(lèi)方法，建立老撾北部5省橡膠幼林、橡膠成林決策樹(shù)提取模型。不同植被覆蓋地類(lèi)具體光譜差異及橡膠林提取將在下述研究中闡述。

2.2.2橡膠林時(shí)空擴(kuò)張分析方法 研究在2010年、2015年橡膠林空間分布疊加基礎(chǔ)上，監(jiān)測(cè)老撾北部5省橡膠林面積時(shí)空擴(kuò)張過(guò)程、擴(kuò)張的地形因子變化特點(diǎn)和研究區(qū)橡膠林?jǐn)U張的驅(qū)動(dòng)因素。橡膠林面積擴(kuò)張分析采用單一土地利用類(lèi)型相對(duì)變化率模型[13]，該模型可定量表達(dá)一類(lèi)土地利用類(lèi)型數(shù)量變化在不同區(qū)域間的差異性。表達(dá)式如下：

(1)

式中：Ka和Kb分別表示某一區(qū)域內(nèi)部某種土地利用類(lèi)型在研究初期和末期的總面積，Ca和Cb分別表示整個(gè)研究區(qū)該土地利用類(lèi)型在研究初期和末期的總面積。研究結(jié)合地形因素的3個(gè)主要因子(海拔、坡度、坡向)構(gòu)建2010年到2015年橡膠林時(shí)空動(dòng)態(tài)變化分析模型，基于該模型掌握5年間老撾北部5省地區(qū)橡膠林時(shí)空變化特征，并從土地利用/土地覆被變化(LUCC)角度探索橡膠林時(shí)空變化驅(qū)動(dòng)機(jī)制。

3 結(jié)果與分析

3.1NDVI時(shí)間序列曲線(xiàn)

歸一化植被指數(shù)(NDVI)可間接反映植被在不同季節(jié)、生育期的不同生理特征，多時(shí)相NDVI時(shí)間序列曲線(xiàn)能夠表現(xiàn)植被在不同生長(zhǎng)時(shí)期的差異性，以及同一時(shí)期不同植被在曲線(xiàn)上的區(qū)別，即植被的物候特征。研究通過(guò)IDL語(yǔ)言批處理生成6個(gè)地類(lèi)樣本的MODIS-NDVI時(shí)間序列，并基于時(shí)間序列諧波分析法(Hants)對(duì)時(shí)間序列曲線(xiàn)平滑去噪(圖2)。

圖2 不同植被覆蓋地類(lèi)NDVI時(shí)間序列Hants濾波曲線(xiàn)Fig.2 Time series NDVI curves of different vegetation coverage land types by Hants

6類(lèi)植被NDVI總體值在8月進(jìn)入低谷，結(jié)合研究區(qū)8月份降水值可知：老撾北部8月份以來(lái)頻頻遭受洪澇災(zāi)害，使得大氣引起的NDVI值總體偏低，目前各種濾波方法對(duì)NDVI時(shí)間序列的單個(gè)大氣噪點(diǎn)平滑效果較好[14]，而對(duì)多期數(shù)據(jù)出現(xiàn)的大氣噪點(diǎn)普遍效果不理想[15]。依據(jù)該區(qū)水田生長(zhǎng)的時(shí)間特點(diǎn)，Hants濾波曲線(xiàn)除8月份外都客觀反映了區(qū)域植被生長(zhǎng)隨時(shí)間變化的特點(diǎn)。濾波后的NDVI時(shí)間序列結(jié)果顯示：水田和旱地 NDVI 時(shí)間序列曲線(xiàn)除去8月份的低谷噪點(diǎn)后為單峰單谷，水田NDVI曲線(xiàn)也表現(xiàn)為5、6月份上升，10月份后迅速下降，與實(shí)地調(diào)查中北部5省水田5月份種植，10月份收割的實(shí)際狀況一致。有林地和橡膠成林的曲線(xiàn)最為接近，有林地NDVI值從11月至次年3月高于橡膠成林且為區(qū)域內(nèi)NDVI值最大的植被地類(lèi)，而在4月下旬至9月下旬的雨季間，其一直低于橡膠成林，表明水分條件是影響橡膠林生長(zhǎng)的重要因子。橡膠幼林和灌木林地的NDVI曲線(xiàn)趨勢(shì)非常相近，且有部分重疊，但在1月上旬至4月上旬，二者的NDVI值差異明顯。綜合橡膠林的落葉物候特征，12月至1月下旬橡膠林NDVI下降至最低后在2月上旬開(kāi)始增長(zhǎng)至4月上旬，表明此間橡膠林正經(jīng)歷著落葉及生長(zhǎng)期，其N(xiāo)DVI值與其他植被區(qū)分明顯。橡膠林物候特征和NDVI曲線(xiàn)地類(lèi)差異表明：1月到3月是橡膠林提取的時(shí)間窗口?；诖耍狙芯拷?jīng)過(guò)對(duì)比選擇8d合成的MODIS-NDVI時(shí)間序列的第5、6、7、8、9、10波段數(shù)據(jù)，即DOY33(DOY表示年積日，從當(dāng)前1月1日開(kāi)始計(jì)天數(shù)，即DOY33表示第33天。)、DOY41、DOY49、DOY57、DOY65、DOY73六期影像平均值，并通過(guò)六期影像相對(duì)應(yīng)的MOD09A1、MOD09Q1第5、6、7、8、9、10波段的數(shù)據(jù)計(jì)算其他指數(shù)平均值，分類(lèi)中同時(shí)采用其他季節(jié)的影像輔助識(shí)別以提高分類(lèi)精度。

3.2橡膠林提取決策樹(shù)模型構(gòu)建

研究在橡膠林提取的最佳時(shí)間窗口基礎(chǔ)上采用所有植被類(lèi)NDVI值大于0.35的特點(diǎn)，以0.35作為閾值提取了2015年植被區(qū)。將對(duì)裸地信息敏感NDBI作為高、低密度植被覆蓋區(qū)的區(qū)分標(biāo)準(zhǔn)，圖3表明：有林地和橡膠成林的NDBI指數(shù)小于0，灌木林地、橡膠幼林、旱地和水田NDBI指數(shù)均大于0，即NDBI值大于0即為低密度覆蓋植被，小于0為高密度覆蓋植被。

圖3 不同地類(lèi)歸一化建筑指數(shù)(NPBI)Fig.3 NDBI of different land types

圖4 不同地類(lèi)在第一主成分上的特征Fig.4 PCA1 of different land types

主成分分析(PCA)是多光譜數(shù)據(jù)的壓縮和去噪工具[16]。6類(lèi)植被第一主成分分量(PCA1)值(圖4)表明：橡膠成林、有林地PCA1值大于0，其余地類(lèi)均小于0，且有林地在PCA1上的值遠(yuǎn)大于0.04，因此，PCA1既可以用來(lái)輔助分離高密度和低密度覆蓋植被，同時(shí)也是區(qū)分橡膠成林與有林地的有效方法。

紋理特征是影像信息提取的重要方法之一。4種主要的紋理特征中，角二階矩反映了圖像灰度分布的均勻程度和紋理的粗細(xì)度，熵是與角二階矩相關(guān)的一個(gè)測(cè)度值，對(duì)比度對(duì)于檢測(cè)圖像邊緣特征和圖像反差具有良好的效果，同質(zhì)性可度量圖像紋理的局部變化大小[17]。影像4種紋理特征(圖5)對(duì)比發(fā)現(xiàn)，有林地熵值和對(duì)比度值皆最??；不同低密度植被在角二階矩和同質(zhì)性上的特征相近，表明其在圖像灰度分布上均勻度一致；對(duì)比度和熵上的變化特征明顯，尤其是對(duì)比度值大于1時(shí)，水田分離效果最好。

圖5 不同地類(lèi)的紋理特征Fig.5 Texture features of different land types

圖6 不同地類(lèi)的原始光譜Fig.6 Original spectral characteristics of different land types

植被原始光譜(圖6)中第1、3、4、7波段上各地類(lèi)平均光譜值較為接近，而第2波段(B2)上灌木林平均反射率最高且大于0.25，旱地在第5波段(B5)上的平均反射率最大且大于0.3，因此，第2波段和第5波段可作為唯一區(qū)分旱地和灌木林的光譜特征?；谏鲜鲅芯恐兄饘犹崛〉南鹉z成林和橡膠幼林，研究最終基于橡膠林落葉期的關(guān)鍵物候特征得到老撾北部橡膠成林、橡膠幼林決策樹(shù)提取方法(圖7)。

圖7 基于MODIS影像的決策樹(shù)分類(lèi)模型Fig.7 Decision tree classification modelbased on MOIDS

3.3精度驗(yàn)證

3.3.1混淆矩陣驗(yàn)證 采用2次野外調(diào)查得到的實(shí)地橡膠林點(diǎn)分布數(shù)據(jù)，結(jié)合最新Google earth高清影像獲取不同地類(lèi)樣本(ROI)。根據(jù)ROI生成的混淆矩陣評(píng)價(jià)老撾北部5省MODIS影像橡膠林分類(lèi)總精度達(dá)到80.67%，Kappa系數(shù)為0.8，橡膠成林和橡膠幼林分類(lèi)總精度分別達(dá)86%和78%(表1)。

表1 分類(lèi)精度評(píng)價(jià)

3.3.2 高分一號(hào)影像驗(yàn)證 進(jìn)而通過(guò)GF-1 WFV影像采用面向?qū)ο蠓椒ㄌ崛±蠐氡辈磕纤∠鹉z林，在影像分割尺度80的基礎(chǔ)上，以NDVI、RVI、DEM和紋理特征建立隸屬度函數(shù)進(jìn)行模糊分類(lèi)，提取老撾北部南塔省橡膠林，提取結(jié)果采用混淆矩陣驗(yàn)證得到的總精度達(dá)88.76%，kappa系數(shù)為0.84，橡膠成林的分類(lèi)精度91.65%，橡膠幼林為81.28%。選取若干像元用GF-1 WFV影像對(duì)MODIS影像提取的2015年南塔省橡膠林分布進(jìn)行對(duì)比驗(yàn)證，結(jié)果表明：基于MODIS數(shù)據(jù)的橡膠林空間分布提取精度為GF-1 WFV影像的91%，同時(shí),利用GF-1 PMS影像通過(guò)面向?qū)ο蠛腿藱C(jī)交互的解譯方法提取南塔市周邊橡膠林，作為GF-1 WFV影像驗(yàn)證結(jié)果的補(bǔ)充，混淆矩陣計(jì)算的總體精度為96.76%，kappa系數(shù)為0.93。研究在GF-1 PMS高精度的分類(lèi)提取基礎(chǔ)上，分別從像元和總面積2個(gè)角度驗(yàn)證了MODIS影像提取橡膠林的精度：(1)當(dāng)GF-1 PMS影像上確定為橡膠林分布的區(qū)域時(shí)，在對(duì)應(yīng)的MODIS影像上隨機(jī)抽取若干像元驗(yàn)證基于MODIS影像的橡膠林空間分布提取精度為82.49%。(2)GF-1 PMS影像覆蓋范圍內(nèi)橡膠林分布面積為1.37萬(wàn)hm2，對(duì)應(yīng)范圍內(nèi)的MODIS提取橡膠林面積為 1.63萬(wàn)hm2，GF-1 PMS影像提取的橡膠林面積為MODIS影像的83.95%。像元和面積2個(gè)角度的驗(yàn)證結(jié)果均與混淆矩陣和基于GF-1 WFV影像驗(yàn)證結(jié)果近似，因此，對(duì)于大面積橡膠林提取，MODIS影像可以在保證時(shí)效性的同時(shí)保證提取精度。基于驗(yàn)證的橡膠林提取算法，研究分別提取了2010年和2015年的橡膠林空間分布(圖8)。

圖8 基于MODIS影像提取的老撾北部5省橡膠林分布示意圖Fig.8 Distribution of rubber plantations in the five provinces of northern Laos based on MODIS

3.4橡膠林時(shí)空擴(kuò)張分析

3.4.1面積擴(kuò)張分析 研究區(qū)橡膠林主要分布于老撾北部南塔省、豐沙里省西部和瑯勃拉邦省北部，同時(shí)在波喬省和烏多姆賽省西部有較小范圍的零星散落分布。2015年橡膠成林面積和橡膠幼林面積分別為14.17、2.79萬(wàn)hm2。2010年橡膠林成林、橡膠幼林面積分別為1.47、11.13萬(wàn)hm2。5年間橡膠林面積增加了4.36萬(wàn)hm2(圖9)，其中，橡膠成林面積從2010年的1.47萬(wàn)hm2增加到2015年的14.17萬(wàn)hm2，而橡膠幼林面積大幅減少了8.34hm2，這表明老撾北部5省2010年的橡膠幼林在2015年大多轉(zhuǎn)化為橡膠成林，而這5年間橡膠幼林面積新增較少表明橡膠種植熱度減退。

圖9 2010—2015年新增種植橡膠林分布示意圖 Fig.9 Newly distributed rubber forest between 2010 and 2015

根據(jù)單一土地利用類(lèi)型相對(duì)變化率模型得到老撾北部5省2010到2015年期間5省橡膠林土地利用的相對(duì)變化率(圖10)。5年間橡膠林土地利用在數(shù)量方面的區(qū)域差異性非常明顯，豐沙里省擴(kuò)張速度最慢，相對(duì)變化率只有58.88%；烏多姆賽省橡膠種植面積擴(kuò)張速度最快，相對(duì)變化率達(dá)到328.85%，其次是波喬省，這說(shuō)明烏多姆賽和波喬成為橡膠種植發(fā)展的新區(qū)域；南塔省和瑯勃拉邦省橡膠林相對(duì)變化率分別為70.16%和74.30%，表明2省橡膠林持續(xù)緩慢增長(zhǎng)。

圖10 2010—2015年橡膠林地利用相對(duì)變化率Fig.10 Rubber forest land use relative change rate between 2010 and 2015

3.4.2地形擴(kuò)張分析 地形分析中原始坡向圖只代表坡面的朝向，值域?yàn)閇0,360][18]，而由于坡向不同所造成的太陽(yáng)輻射能量不同是影響橡膠林分布的主要原因，因此，研究對(duì)坡向轉(zhuǎn)換后生成Trasp圖[19]，Trasp計(jì)算公式如下：

(2)

式(2)中：Trasp值域?yàn)閇0,1]，0表示接收到的太陽(yáng)輻射能最少，1表示接收的太陽(yáng)輻射能最大，aspect為原始坡向圖。

表2 2010—2015年橡膠林在不同海拔梯度上的空間分布比例

表3 2010—2015年橡膠林在不同坡度區(qū)間上的空間分布比例

表4 2010—2015年橡膠林在不同Trsap上的空間分布比例

3.4.3 擴(kuò)張驅(qū)動(dòng)因素分析 2015年橡膠成林源于2010年橡膠幼林，2015年橡膠幼林面積大幅減少表明2010年后老撾北部5省橡膠林定植面積很少。橡膠林在面積和地形上的擴(kuò)張?zhí)卣黠@示其生長(zhǎng)的立地條件不斷變化。橡膠林?jǐn)U張變化是LUCC的一部分，自然環(huán)境和人類(lèi)活動(dòng)的干擾[21]是LUCC變化的驅(qū)動(dòng)因素，其中，人類(lèi)活動(dòng)是LUCC變化的主要驅(qū)動(dòng)因素，核心是土地利用政策和管理主體的變化[22]，總結(jié)老撾北部5省橡膠林面積變化的核心包括社會(huì)政策和社會(huì)經(jīng)濟(jì)兩方面。

(1)為從源頭上根除老撾北部地區(qū)罌粟對(duì)中國(guó)及國(guó)際社會(huì)的威脅，中國(guó)政府積極在老撾北部地區(qū)發(fā)展“罌粟替代種植”，在經(jīng)濟(jì)利益和工業(yè)發(fā)展的驅(qū)動(dòng)下，橡膠林從傳統(tǒng)種植區(qū)不斷向高緯度地區(qū)擴(kuò)展，因此，2005年之后，老撾北部橡膠種植勢(shì)如破竹；但瑯勃拉邦省、南塔省分別于2008年2月和11月宣布暫停審批橡膠種植項(xiàng)目，老撾政府首先認(rèn)為橡膠規(guī)模種植所帶來(lái)的經(jīng)濟(jì)和社會(huì)效益并不明顯，需要認(rèn)真評(píng)價(jià)橡膠種植的經(jīng)濟(jì)、社會(huì)效益；其次由于2008年2月份橡膠價(jià)格降幅超過(guò)30%，政府擔(dān)心橡膠開(kāi)割期依然低迷的橡膠市場(chǎng)，鼓勵(lì)當(dāng)?shù)胤N植見(jiàn)效快的農(nóng)作物；再者，老撾國(guó)家土地署表明，老撾勞動(dòng)力嚴(yán)重不足，前期的橡膠種植對(duì)經(jīng)濟(jì)、社會(huì)、環(huán)境影響沒(méi)有做充分的可研，由此老撾北部南塔省和瑯勃拉邦省2015年橡膠幼林面積非常少。此外，2012年老撾政府暫停全國(guó)橡膠、桉樹(shù)種植、礦產(chǎn)開(kāi)采的土地許可，導(dǎo)致2012年之后老撾北部5省橡膠幼林少有增加，2015年橡膠幼林面積相比2010年大幅減少。

(2)根據(jù)天然橡膠生產(chǎn)國(guó)組織發(fā)布的數(shù)據(jù)得到2009年至2016年平均天然標(biāo)準(zhǔn)橡膠價(jià)格(圖11)。從2009年開(kāi)始，天然標(biāo)準(zhǔn)橡膠價(jià)格經(jīng)過(guò)2年多的持續(xù)上升到2011年3月突破歷史最高紀(jì)錄，達(dá)到每噸32 000元·t-1，漲幅高達(dá)158%，此時(shí)部分膠農(nóng)受價(jià)格利益驅(qū)動(dòng)有少量橡膠定植，而后天然橡膠價(jià)格一路大幅跌落，老撾北部橡膠種植不再快速蔓延。因此，基于社會(huì)政策和社會(huì)經(jīng)濟(jì)因素，老撾北部5省2015年橡膠幼林面積僅有2.79萬(wàn)hm2，占2015年橡膠林總面積的16.46%。

圖11 2009至2016年間標(biāo)準(zhǔn)膠價(jià)格Fig.11 Standard rubber price between 2009 and 2016

4 討論

精度在遙感信息提取中是重要的描述參數(shù)。研究采用250 m分辨率結(jié)合500 m分辨率的MODIS遙感影像提取橡膠林，精度達(dá)到80.67%，在同類(lèi)影像中相對(duì)較高；而老撾北部統(tǒng)計(jì)數(shù)據(jù)顯示，2010年和2015年橡膠林面積分別為10.61、14.68萬(wàn)hm2，較研究提取的2010年和2015年橡膠林面積12.6、16.96萬(wàn)hm2相比偏小，可能原因有兩方面：一是MODIS影像屬低分辨率遙感影像，一個(gè)像元的大小為250 m×250 m或500 m×500 m，因此，提取過(guò)程中容易將橡膠林林間路、林間地和較小的山體溝壑部分都誤分成橡膠林；二是統(tǒng)計(jì)數(shù)據(jù)存在欠缺狀況，政府部門(mén)的統(tǒng)計(jì)數(shù)據(jù)大多為大規(guī)模橡膠園，農(nóng)民零散種植的面積較小加之上報(bào)數(shù)據(jù)存在偏差，導(dǎo)致往往存在漏統(tǒng)現(xiàn)象，因此研究中提取的橡膠林面積偏大。

研究從面積和地形兩方面提取了老撾北部5省橡膠林空間分布及時(shí)空變化特征。下一步的研究將考慮更多影響橡膠林空間分布的自然因子，如地貌因素中的土壤因子，氣候因素中的溫度、降水、濕度因子等，從而分析橡膠林空間分布的多自然因子變化特征，進(jìn)而定量分析自然因子對(duì)橡膠林空間分布的影響，得到橡膠林空間分布驅(qū)動(dòng)模型。

隨著高分辨率遙感影像的易獲取性，橡膠林空間分布提取在決策樹(shù)分類(lèi)模型基礎(chǔ)上，可融入專(zhuān)家知識(shí)，非遙感決策函數(shù)等建立更高精度區(qū)域性強(qiáng)普適性的遙感橡膠林提取模型，提高橡膠林分類(lèi)精度，以滿(mǎn)足橡膠林分布區(qū)土地資源管理和環(huán)境保護(hù)的需要。

5 結(jié)論

研究采用遙感數(shù)據(jù)結(jié)合橡膠林物候特征實(shí)現(xiàn)了橡膠林空間分布快速可重復(fù)提取的新方法，并基于社會(huì)和經(jīng)濟(jì)數(shù)據(jù)建立了橡膠林時(shí)空擴(kuò)張模型。在此基礎(chǔ)上，探索了橡膠林時(shí)空擴(kuò)張?zhí)卣骷捌潋?qū)動(dòng)機(jī)制。

(1)1月至3月是橡膠林提取的時(shí)間窗口，結(jié)合NDVI、NDBI、PCA、影像的紋理特征和原始光譜值等構(gòu)建的橡膠林提取決策樹(shù)算法分類(lèi)總精度較高，達(dá)到80.67%，Kappa系數(shù)為0.8。

(2)基于高分一號(hào)GF-1WFV和GF-1PMS影像對(duì)MODIS影像提取的2015年南塔省橡膠林分布從像元和面積2個(gè)角度進(jìn)行驗(yàn)證，驗(yàn)證結(jié)果在80.67%左右浮動(dòng)，不同驗(yàn)證方法充分表明對(duì)大面積橡膠林提取MODIS數(shù)據(jù)可以在保證時(shí)效性的同時(shí)兼顧精度。

(3)2010年到2015年橡膠林時(shí)空擴(kuò)張?zhí)攸c(diǎn)表明，老撾北部5省橡膠林空間分布在海拔和坡度上總體呈現(xiàn)不斷升高的趨勢(shì)，海拔1 500 m以上只有零星分布，坡度分布上限為35°，坡向上分布無(wú)明顯規(guī)律。從土地利用變化角度分析得到，老撾北部5省的社會(huì)政策和經(jīng)濟(jì)因素是影響橡膠林時(shí)空擴(kuò)張的主要驅(qū)動(dòng)因素。

[1] 柯佑鵬, 過(guò)建春, 張賽麗. 世界天然像膠經(jīng)濟(jì)研究[M]. 北京: 經(jīng)濟(jì)科學(xué)出版社, 2010.

[2] 姚元園. 東南亞天然橡膠產(chǎn)業(yè)研究[D]. 廈門(mén): 廈門(mén)大學(xué), 2014.

[3] 劉洪江, 蘭恒星, 張 軍，等. 老撾北部罌粟替代種植高分辨率遙感調(diào)查評(píng)價(jià)與分析[J]. 資源科學(xué),2010,32(7):1425-1432.

[4] 董加相,巖坎糯. 西雙版納對(duì)老撾罌粟替代種植面臨的機(jī)遇與挑戰(zhàn)[J]. 時(shí)代金融,2007(1):126-127.

[5] Li H M, Aide T M, Ma Y X,etal. Demand for rubber is causing the loss of high diversity rain forest in SW China[J]. Biodiversity and Conservation, 2007, 16(6):1731-1745.

[6] 陳匯林, 陳小敏, 陳珍麗, 等. 基于MODIS遙感數(shù)據(jù)提取海南橡膠信息初步研究[J]. 熱帶作物學(xué)報(bào), 2010, 31(7): 1181-1185.

[7] 張京紅, 陶忠良, 陳少軍, 等. 基于TM影像的海南島橡膠種植面積信息提取[J].熱帶作物學(xué)報(bào), 2010, 31(4): 661-665.

[8] Ekadinata A, Widayati A, Vincent G. Rubber agroforest identification using object-based classification in Bungo District, Jambi, Indonesia. Proceedings of the 25th Asian Conference on Remote Sensing[C]. Chiang Mai, Thailand: ACRS, 2004: 551-556.

[9] 劉少軍, 張京紅, 何政偉, 等. 基于面向?qū)ο蟮南鹉z分布面積估算研究[J]. 廣東農(nóng)業(yè)科學(xué),2010, 37(1): 168-170.

[10] 劉曉娜, 封志明, 姜魯光, 等. 西雙版納橡膠林地的遙感識(shí)別與數(shù)字制圖[J]. 資源科學(xué), 2012, 34(9): 1769-1780.

[11] 李亞飛, 劉高煥, 黃 翀. 基于HJ-1 CCD數(shù)據(jù)的西雙版納地區(qū)橡膠林分布特征[J]. 中國(guó)科學(xué): 信息科學(xué), 2011, 41(增刊): 166-176.

[12] 任 沖,鞠洪波,張懷清,等. 天水市近30年林地動(dòng)態(tài)變化遙感監(jiān)測(cè)研究[J]. 林業(yè)科學(xué)研究,2017,30(1):25-33.

[13] 王 越,宋 戈,王盼盼. 松嫩高平原土地利用類(lèi)型變化特征及時(shí)空格局研究[J]. 干旱區(qū)資源與環(huán)境,2017,31(4):61-66.

[14] 王 丹,姜小光,唐伶俐,等. 利用時(shí)間序列傅立葉分析重構(gòu)無(wú)云NDVI圖像[J]. 國(guó)土資源遙感,2005,17(2):29-32.

[15] 侯光雷,張洪巖,王野喬,等. 基于時(shí)間序列諧波分析的東北地區(qū)耕地資源提取[J]. 自然資源學(xué)報(bào),2010,25(9):1607-1617.

[16] Zhao G, Maclean A L. A Comparison of Canonical Discriminant Analysis and Principal Component Analysis for Spectral Transformation[J]. Photogrmmetric Engineering &Remote Sensing, 2000,66(7): 841-847.

[17] 李智峰,朱谷昌,董泰鋒. 基于灰度共生矩陣的圖像紋理特征地物分類(lèi)應(yīng)用[J]. 地質(zhì)與勘探,2011,47(3):456-461.

[18] Leng W F, He H S, Bu R C. The spatial distribution of constructive species of Northeast forest under the climate changing[J]. Acta Ecol Sin, 2006, 26(12): 4258-4265.

[19] Xu Y J, Chen Y N, Li W H,etal. Distribution pattern of plant species diversity in the mountainous region of I li River Valley, Xinjiang[J]. Environmental Monitoring and Assessment, 2010,177(4), 681-694.

[20] GB/T 16453.2-2008.中華人民共和國(guó)國(guó)家標(biāo)準(zhǔn).水土保持綜合治理[S].

[21] Sun L, Wei J, Duan D H,etal. Impact of Land-Use and Land-Cover Change on urban air quality in representative cities of China[J]. Journal of Atmospheric and Solar-Terrestrial Physics,2016,142:43-54.

[22] 方廣玲,吳 斌,張宇清,等. 鹽池縣土地利用/覆被變化的人文因素分析[J]. 中國(guó)水土保持科學(xué),2011,9(1):74-81.

(責(zé)任編輯：張 玲)

ResearchonExtractionandSpatial-TemporalExpansionofRubberForestinFiveProvincesofNorthernLaosBasedonMulti-sourceRemoteSensing

LIYang-yang，ZHANGJun，LIUChen-li，YANGXu-chao，LIJie

(Department of Environment and Resource, Yunnan University, Kunming 650504, Yunnan, China)

S794.1

A

1001-1498(2017)05-0709-09

10.13275/j.cnki.lykxyj.2017.05.002

2017-04-15

云南省應(yīng)用基礎(chǔ)研究計(jì)劃項(xiàng)目“中高分辨率遙感數(shù)據(jù)支持下的橡膠林地信息定量識(shí)別方法研究(2013FZ002)”

李陽(yáng)陽(yáng)(1990—)，女. 碩士，主要研究方向?yàn)檫b感與地理信息系統(tǒng)應(yīng)用. E-mail：978098813@qq.com

* 通訊作者：張 軍(1972—)，男，博士，副研究員. E-mail：50981534@qq.com