季建萬, 沙晉明, 金 彪,2, 包忠聰,3

1(福建師范大學(xué) 地理科學(xué)學(xué)院,福州 350007)

2(福建師范大學(xué) 軟件學(xué)院,福州 350108)

3(福州市勘測(cè)院,福州 350108)

1 引言

土地利用是土地可持續(xù)發(fā)展和規(guī)劃的至關(guān)重要的組成部分,以前的土地利用圖大多基于野外調(diào)查與室內(nèi)繪制,這不僅耗費(fèi)大量的時(shí)間和金錢,而且效率低下.因此,一般大范圍的土地利用圖年代比較久遠(yuǎn),遠(yuǎn)不能滿足現(xiàn)代化建設(shè)的需求[1].

遙感技術(shù)因其低成本、高效率的特征,已經(jīng)成為土地利用信息提取的重要途徑. 許多學(xué)者都已經(jīng)運(yùn)用不同的基于像元的分類方法提取得到土地利用信息[2].比如孫琳等基于HJ-1B遙感影像,運(yùn)用最大似然法和支持向量機(jī)方提取太湖流域土地利用信息[3],馬友平等基于變差函數(shù)紋理和BP人工神經(jīng)網(wǎng)絡(luò)方法對(duì)QuickBird影像進(jìn)行分類研究[4],陳文倩等運(yùn)用結(jié)合影像紋理信息的支持向量機(jī)法對(duì)新疆艾比湖濕地保護(hù)區(qū)進(jìn)行土地利用覆蓋分類研究[5]等. 但是對(duì)于高分辨率的遙感影像,運(yùn)用常規(guī)的基于像元的分類方法效果卻不是非常理想,而面向?qū)ο蠓诸惙椒ú粌H考慮影像的光譜特征,而且考慮影像的紋理特征與幾何特征,因此分類結(jié)果具有良好的整體性.

本研究旨在利用面向?qū)ο蠓椒ê退姆N傳統(tǒng)的基于像元的分類方法分別提取研究區(qū)土地利用信息,并以目視解譯結(jié)果為參考從分配分歧和數(shù)量分歧兩方面對(duì)各種分類方法進(jìn)行精度評(píng)價(jià),為今后的土地利用信息提取提供一定的建議.

2 研究區(qū)概況

研究區(qū)位于福州市平潭縣中北部,地理位置為119°45′4″E～119°46′18″E,25°33′54″N～25°35′40″N,總面積為6.56 km2(圖1). 研究區(qū)基本覆蓋了平潭縣典型地物類型,且為了減少數(shù)據(jù)量,提高分析效率,選擇該研究區(qū)作為研究對(duì)象進(jìn)行分析.

圖1 研究區(qū)位置(左)和標(biāo)準(zhǔn)假彩色影像(右)

3 數(shù)據(jù)與方法

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

本研究采用WorldView-2衛(wèi)星的2 m多光譜影像為數(shù)據(jù)源(影像獲取時(shí)間為2011年4月11日,03：16：33,投影坐標(biāo)系為 WGS_1984 UTM_Zone_50N).WorldView-2影像有4個(gè)多光譜波段,波段范圍分別為： 藍(lán)光波段 450～510 nm,綠光波段 510～580 nm,紅光波段為630～590 nm和近紅外波段770～895 nm. 通過ENVI軟件對(duì)WorldView-2影像進(jìn)行幾何校正、輻射定標(biāo)、FLAASH大氣校正和影像裁剪等預(yù)處理操作,并利用預(yù)處理后的遙感影像為數(shù)據(jù)源展開研究.

3.2 研究區(qū)土地利用分類系統(tǒng)

依據(jù)GB/T21010-2007《土地利用現(xiàn)狀分類》系統(tǒng)[6],并結(jié)合研究區(qū)實(shí)際情況,將研究區(qū)分為7類,如表1所示.

表1 研究區(qū)土地利用分類系統(tǒng)

3.3 研究區(qū)土地利用分類方法

3.3.1 研究區(qū)土地利用分類流程

研究表明,利用影像紋理特征輔助多光譜特征對(duì)影像進(jìn)行分類能夠有效提高影像分類精度[7]. 本文通過提取WorldView-2遙感影像近紅外波段、紅波段、綠波段和藍(lán)波段紋理信息,將其與經(jīng)過大氣校正后的波段進(jìn)行疊加,并利用ENVI和eCognition軟件完成研究區(qū)土地利用信息的提取,具體流程如圖2所示.

3.3.2 影像紋理信息的提取

紋理特征提取的主要方法有統(tǒng)計(jì)方法、模型方法、信號(hào)處理方法和結(jié)構(gòu)方法,其中統(tǒng)計(jì)方法是基于像元及其領(lǐng)域的灰度屬性,研究紋理區(qū)域中的統(tǒng)計(jì)特性[8]. 而灰度共生矩陣法(GLCM)正是通過研究灰度的空間相關(guān)特性從而實(shí)現(xiàn)影像紋理的描述與提取[9]. 本文運(yùn)用灰度共生矩陣法提取WorldView-2影像的紋理信息,具體包括信息熵、角二階矩、對(duì)比度、相關(guān)性、方差和相異性,其中紋理窗口以及移動(dòng)步長(zhǎng)的選擇對(duì)紋理信息的提取結(jié)果及后續(xù)的影像分類有著非常重要的影響,較大的窗口保證了同類地物紋理統(tǒng)計(jì)的代表性,適合光譜異質(zhì)性強(qiáng)的類別; 較小的窗口保證了窗口內(nèi)像素的同質(zhì)性,適合光譜均勻性強(qiáng)的地物. 本文在ENVI軟件下,先固定移動(dòng)步長(zhǎng)為(1,1),以2為間隔,得到最佳的紋理窗口為5×5,同理,固定最佳紋理窗口,以2為步長(zhǎng),得到最佳的移動(dòng)步長(zhǎng)為(1,1),圖3是以WorldView-2影像藍(lán)光波段為例,得到的6種紋理信息提取結(jié)果.

圖2 WorldView-2遙感影像土地利用信息提取流程圖

圖3 藍(lán)光波段紋理特征提取結(jié)果

3.3.3 面向?qū)ο蠓诸惙椒?/h4>

利用eCognition軟件對(duì)研究區(qū)進(jìn)行面向?qū)ο蠓诸?主要步驟有確定分割尺度、建立分類規(guī)則、分類后處理等. 尺度分割是面向?qū)ο蠓诸惍?dāng)中重要的一步,分割尺度的選擇直接影響到分類結(jié)果,分割尺度過小,影像被分割的過細(xì),分割尺度過大,影像被分割的過疏. 本研究采用 Lucian Drǎgu?提出的 ESP (Estimation of Scale Parameter)尺度分割工具來確定最佳分割尺度[10],該工具是以局部方差的標(biāo)準(zhǔn)偏差均值來代表某塊影像檢查窗口中對(duì)象的異質(zhì)性,并通過迭代過程得到最優(yōu)分割尺度[11],實(shí)驗(yàn)結(jié)果表明研究區(qū)最佳分割尺度為19. 分割完成后,先固定形狀指數(shù)(shape),以0.1為步長(zhǎng),得到最佳的緊致度指數(shù)(compactness)為0.5,光滑度指數(shù)(roundness)為0.5; 同理,固定最佳的緊致度指數(shù),以0.1為步長(zhǎng),得到最佳的形狀指數(shù)(shape)為0.2,色彩指數(shù)(color)為0.8. 分類規(guī)則的建立是面向?qū)ο蠓诸愡^程當(dāng)中又一極其重要的一步. 高分辨率影像由于地物的光譜、形狀、紋理等差異性增大,因此必須綜合多種特征和知識(shí)來進(jìn)行判別和提取. 本文所參與土地利用信息提取的光譜特征包括歸一化土壤指數(shù)(NDSI)、歸一化植被指數(shù)(NDVI)、歸一化水體指數(shù)(NDWI)、陰影指數(shù)(SI)、各波段貢獻(xiàn)率、各波段標(biāo)準(zhǔn)差、各波段平均值、各波段偏度與亮度(Brightness);紋理特征包括角二階矩、信息熵、相關(guān)性、平均值及對(duì)比度; 形狀特征包括面積、長(zhǎng)寬比、形狀指數(shù)、密度、寬度、邊界長(zhǎng)度.

將WorldView-2影像的各種特征信息導(dǎo)入到Salford Predictive Modeler(SPM)軟件中,采用CART決策樹方法[12],經(jīng)過多次的試驗(yàn),得到最終的研究區(qū)WorldView-2影像的分類規(guī)則,如圖4所示. 分類完成后,仍需要對(duì)分類結(jié)果進(jìn)行后處理操作,具體包括類別的合并,細(xì)小對(duì)象的修正以及其他后處理操作,最終的分類結(jié)果如圖5(a)所示.

圖4 WorldView-2遙感影像面向?qū)ο蠓诸愐?guī)則

3.3.4 基于像元的分類方法

基于像元的分類方法是指以像元為分類基本單元并基于一定的計(jì)算規(guī)則從而得到分類結(jié)果. 本文以神經(jīng)網(wǎng)絡(luò)分類法、支持向量機(jī)法、最大似然法和馬氏距離法作為代表,在ENVI軟件中,利用疊加紋理波段后影像作為輸入影像完成影像的分類工作. 各種方法的分類結(jié)果如圖5(b)-圖5(e)所示.

3.3.5 目視解譯法

基于團(tuán)隊(duì)成員實(shí)際野外考察結(jié)果,并利用人工目視解譯法,經(jīng)過多次修改,得到研究區(qū)的目視解譯結(jié)果,如圖5(f)所示.

3.4 分類后精度評(píng)價(jià)

本文利用400個(gè)樣本點(diǎn)進(jìn)行精度檢驗(yàn),由于不同類別面積不同,采用分層比例采樣法得到各類別參與檢驗(yàn)的樣本數(shù)量. 具體處理方法為依據(jù)目視解譯圖,計(jì)算研究區(qū)各類別面積和所占比例,并依據(jù)比例得到各類別樣本數(shù)量,結(jié)果如表2所示.

表2 各類別隨機(jī)驗(yàn)證樣點(diǎn)數(shù)量統(tǒng)計(jì)

利用ArcGIS 10.3軟件為各樣本點(diǎn)賦予相對(duì)應(yīng)的屬性值,依據(jù)文獻(xiàn)[13]提出的計(jì)算方法,得到不同分類方法的樣本矩陣,并依據(jù)式(1)將其輸出為總體矩陣,因?yàn)閷?duì)于一幅高分辨率遙感影像而言,包含數(shù)量巨大的柵格點(diǎn),而樣本矩陣在實(shí)際操作中很難考慮所有樣本點(diǎn),因此通過式(1)將其轉(zhuǎn)換到比例尺度上進(jìn)行分析,得到更具代表性的整體評(píng)價(jià)矩陣.

其中pij代表整個(gè)研究區(qū)的評(píng)價(jià)比例,i是指對(duì)比圖中的第i類,j是指參考圖中的第j類,J是指類別總數(shù),Ni是指對(duì)比圖中第i類的樣本數(shù),nij是指樣本矩陣中分類為第i行第j列的樣本數(shù)量[13].

數(shù)量分歧(quantity disagreement)是指由于類別比例未達(dá)到最佳匹配而導(dǎo)致參考圖與對(duì)比圖之間的差異;分配分歧(allocation disagreement)是指由于類別空間分布未達(dá)到最佳匹配而導(dǎo)致參考圖與對(duì)比圖之間的差異[12]. 國(guó)外已有學(xué)者應(yīng)用數(shù)量分歧和分配分歧指標(biāo)對(duì)影像分類方法結(jié)果進(jìn)行評(píng)價(jià)[14-16],而國(guó)內(nèi)學(xué)者一般仍以總體分類精度的高低作為評(píng)價(jià)影像分類結(jié)果的標(biāo)準(zhǔn),但僅以總體分類精度作為評(píng)價(jià)標(biāo)準(zhǔn)并不能全面地了解分類結(jié)果的內(nèi)部情況,因此本文結(jié)合總體分類精度和數(shù)量分歧與分配分歧指標(biāo)對(duì)各種分類方法進(jìn)行精度評(píng)價(jià),式(2)至式(5)為計(jì)算公式.

其中g(shù)是指任意類別,qg是任意g類別的數(shù)量分歧,Q是所有分類類別總的數(shù)量分歧;ag是任意g類別的分配分歧,A是所有分類類別總的分配分歧[13].

4 結(jié)果與分析

4.1 不同分類方法精度比較

在總體矩陣的基礎(chǔ)上,計(jì)算得到各分類方法的制圖精度、用戶精度與總體精度,如表3所示.

表3 不同土地利用分類方法精度統(tǒng)計(jì)(單位： %)

結(jié)合圖5與表3可以得到,5種分類方法平均總體精度為75.00%,其中最大似然法分類效果最差,總體精度僅為62.00%,將其與目視解譯圖相比,得到田坎與耕地混淆現(xiàn)象嚴(yán)重,可能是因?yàn)樘锟惨话惴植荚诟刂車?且易生長(zhǎng)雜草,使兩者的地物光譜呈現(xiàn)較為相似的特征; 也有可能是因?yàn)椴煌母胤N植的作物以及收割時(shí)間一般不同步,造成同種地物類型呈現(xiàn)不同的光譜特征. 此外,灌木林地與農(nóng)村道路混淆嚴(yán)重,這可能是因?yàn)樵诘缆穬膳陨L(zhǎng)有植被,使其在影像上容易混淆. 面向?qū)ο蠓诸惙ǚ诸惪傮w精度最高,為84.25%,其分類結(jié)果圖中,田坎與耕地也比較容易相混淆,表明高分辨率的影像因其像元分辨率的提高,增加了同類之間的光譜異質(zhì)性,減少了異類之間的光譜同質(zhì)性,造成了“同物異譜”、“同譜異物”地物的區(qū)分難度[17]. 對(duì)比面向?qū)ο蠓ㄅc傳統(tǒng)的四種基于像元的分類方法,可以得到基于像元的分類結(jié)果存在明顯的“椒鹽效應(yīng)”,而面向?qū)ο蠓诸惤Y(jié)果卻有效避免了“椒鹽效應(yīng)”,這是因?yàn)槊嫦驅(qū)ο蠓诸惙椒ㄊ且詫?duì)象為分類基本單元,綜合利用光譜、紋理和形狀特征對(duì)影像進(jìn)行分割和分類,從而使其分類結(jié)果更具有完整性.

4.2 不同分類方法整體數(shù)量分歧和分配分歧結(jié)果比較

運(yùn)用數(shù)量分歧和分配分歧兩種指標(biāo)對(duì)不同分類方法整體進(jìn)行評(píng)價(jià),評(píng)價(jià)結(jié)果如圖6所示.

圖6 不同分類方法整體精度評(píng)價(jià)結(jié)果

由圖6可以得到,在基于像元的分類方法中,神經(jīng)網(wǎng)絡(luò)分類方法具有最低的數(shù)量分歧,約11.25%,其次為支持向量機(jī)(15.00%)、馬氏距離法(21.75%)和最大似然法(31.00%). 與其相比,面向?qū)ο蠓诸惙椒ň哂凶畹偷臄?shù)量分歧,僅4.25%. 從分配分歧結(jié)果來看,依次為面向?qū)ο蠓?11.50%)>馬氏距離法(8.75%)>神經(jīng)網(wǎng)絡(luò)法(8.50%)>最大似然法(7.00%)>支持向量機(jī)法(5.75%). 總體來說,數(shù)量分歧是這五種分類方法精度差異的主要原因,因?yàn)閿?shù)量分歧值越高,表明分類錯(cuò)誤的像元數(shù)量越多. 以神經(jīng)網(wǎng)絡(luò)法與面向?qū)ο蠓ǚ诸惤Y(jié)果為例,兩者的分配分歧僅有3%的差距,但是數(shù)量分歧卻相差7%,即神經(jīng)網(wǎng)絡(luò)分類結(jié)果中分類錯(cuò)誤的像元數(shù)量比面向?qū)ο蠓诸惤Y(jié)果中分類錯(cuò)誤的像元數(shù)量要多7%.

4.3 不同分類方法各類別數(shù)量分歧和分配分歧結(jié)果比較

運(yùn)用數(shù)量分歧和分配分歧對(duì)不同分類方法各地物類別進(jìn)行分析,評(píng)價(jià)結(jié)果如圖7所示.

圖7 不同分類方法各地物類別精度評(píng)價(jià)結(jié)果

分析單類別可以得到,不同的分類方法當(dāng)中,耕地的數(shù)量分歧大小依次為最大似然法(28.75%)>馬氏距離法(21.50%)>支持向量機(jī)法(14.75%)>神經(jīng)網(wǎng)絡(luò)法(11.00%)>面向?qū)ο蠓?3.00%); 田坎的數(shù)量分歧比例大小依次為最大似然法(13.50%)>馬氏距離法(6.75%)>支持向量機(jī)法(6.00%)>神經(jīng)網(wǎng)絡(luò)法(4.50%)>面向?qū)ο蠓?3.00%). 對(duì)比可以得到田坎與耕地的數(shù)量分歧大小排序相同,表明這兩種地物類型的數(shù)量分歧對(duì)不同分類方法最終精度影像較大,分析可能是因?yàn)椴煌馗髌诘牟煌蛊湓诟叻直媛视跋裆铣尸F(xiàn)較大的光譜差異; 田坎則主要是因?yàn)槠淇臻g分布的特殊性,大部分位于耕地與耕地之間,且田坎一般寬度較窄,易與周圍地物相混合而形成混合像元,使兩者的分類結(jié)果較差. 從分配分歧結(jié)果分析,耕地的分配分歧大小依次為面向?qū)ο蠓?10.50%)>神經(jīng)網(wǎng)絡(luò)法(5.00%)>支持向量機(jī)法(1.50%)>最大似然法(0.50%)>馬氏距離法(0.00%),灌木林地的分配分歧比例大小依次為馬氏距離法(9.50%)>最大似然法(7.50%)>面向?qū)ο蠓?4.50%)>支持向量機(jī)法(4.00%)>神經(jīng)網(wǎng)絡(luò)法(3.50%),表明面向?qū)ο蠓傮w分類精度主要受耕地與灌木林地分配分歧的影響.

5 結(jié)論與討論

本文利用WorldView-2高分辨率影像數(shù)據(jù)不同的分類方法對(duì)平潭縣中北部地區(qū)進(jìn)行土地利用信息提取,并從數(shù)量分歧和分配分歧兩方面對(duì)分類結(jié)果進(jìn)行精度評(píng)價(jià),結(jié)論如下：

(1)不同分類方法平均總體精度75.00%,其中總體精度最高的是面向?qū)ο蠓诸惙椒?為84.25%,最低的是最大似然法,僅為62.00%.

(2)面向?qū)ο蠓诸惙椒ň哂凶畹偷臄?shù)量分歧,為4.25%,其次依次為神經(jīng)網(wǎng)絡(luò)法＜支持向量機(jī)法＜馬氏距離法＜最大似然法. 在分配分歧方面,支持向量機(jī)方法其值最低,為5.75%,其次依次為最大似然法＜神經(jīng)網(wǎng)絡(luò)法＜馬氏距離法＜面向?qū)ο蠓?

(3)耕地的數(shù)量分歧大小依次為最大似然法(28.75%)>馬氏距離法(21.50%)>支持向量機(jī)法(14.75%)>神經(jīng)網(wǎng)絡(luò)法(11.00%)>面向?qū)ο蠓?3.00%).表明在影像中,耕地像元被錯(cuò)誤分為其他土地利用類別的數(shù)量比例最大,這是因?yàn)樵谟跋癞?dāng)中,耕地呈現(xiàn)出不同的耕作狀態(tài),造成其同物異譜現(xiàn)象嚴(yán)重. 在分配分歧方面,耕地的分配分歧比例大小依次為面向?qū)ο蠓?10.50%)>神經(jīng)網(wǎng)絡(luò)法(5.00%)>支持向量機(jī)法(1.50%)>最大似然法(0.50%)>馬氏距離法(0.00%). 表明耕地同樣是影響不同分類方法分配分歧大小的重要類別,與數(shù)量分歧不同,分配分歧對(duì)面向?qū)ο蠓ǖ挠绊懨黠@高于其他方法,究其原因可能是面向?qū)ο蠓ㄊ且苑指顚?duì)象為分類基礎(chǔ)單元,其忽略了部分光譜異質(zhì)性信息.

結(jié)合上文,面向?qū)ο蠓椒ㄆ淇傮w精度最高,且其分類過程中綜合考慮了影像的紋理特征、光譜特征和形狀特征,因此對(duì)于WorldView-2影像而言,能夠比傳統(tǒng)的基于像元的分類方法更充分地利用影像包含的信息.同時(shí)需要注意的是,與四種傳統(tǒng)的基于像元的分類方法相比,雖然面向?qū)ο蠓ㄆ浞诸惤Y(jié)果往往具有完整性且避免了“椒鹽效應(yīng)”,但是也易存在丟失一些光譜信息. 因?yàn)樵谀承┑匚镱愋蛢?nèi)部可能存在一些單獨(dú)像元,其在現(xiàn)實(shí)中與周圍地物類型并不相同,而在面向?qū)ο蠓诸愡^程中,也較容易忽視這些光譜信息,這在未來的工作中應(yīng)該進(jìn)行更深一步的研究.

對(duì)于傳統(tǒng)的基于像元的分類方法,適合于分析多光譜影像和中低尺度分辨率的遙感影像,尤其在分析大范圍區(qū)域面積變化等方面具有一定的優(yōu)勢(shì),且大多數(shù)遙感軟件已經(jīng)有相關(guān)的集成操作,在提取中低分辨率遙感影像中效率較高且效果也較好,但是對(duì)于高分辨率遙感影像,豐富的光譜信息使得明顯存在“異物同譜”、“同譜異物”的現(xiàn)象,因此在未來的遙感影像信息提取過程中,需要更加充分地挖掘高分辨率遙感影像中的光譜信息及紋理等其他信息.

此外,本文僅選取了一處研究區(qū)進(jìn)行分析,其結(jié)果可能存在一定的誤差,在后續(xù)的工作中,應(yīng)該對(duì)研究區(qū)的代表性進(jìn)行更加深入的分析. 此外,還應(yīng)該選取檢驗(yàn)區(qū)域使結(jié)果更具有科學(xué)性. 在精度分析方面,應(yīng)考慮結(jié)合多種評(píng)價(jià)指標(biāo)對(duì)分類結(jié)果進(jìn)行更加全面的判斷分析.此外,目視解譯圖作為參考圖,其結(jié)果準(zhǔn)確與否與最終的評(píng)價(jià)結(jié)果密切相關(guān),后續(xù)應(yīng)該繼續(xù)提高目視解譯結(jié)果的準(zhǔn)確性,使其更加接近真實(shí)地物類別.

1Estoque RC,Estoque RS,Murayama YJ. Prioritizing areas for rehabilitation by monitoring change in barangay-based vegetation cover. ISPRS International Journal of Geo-Information,2012,1(3)： 46-68. [doi： 10.3390/ijgi1010046]

2韓濤. 遙感監(jiān)測(cè)土地覆蓋變化的方法及研究進(jìn)展. 干旱氣象,2004,22(2)： 76-81.

3孫琳,宋愛紅,高文秀. 基于最大似然法和SVM法的太湖流域HJ-1B影像分類. 測(cè)繪信息與工程,2012,37(3)：30-33.

4馬友平. 基于變差函數(shù)紋理和BP人工神經(jīng)網(wǎng)絡(luò)的QuickBird影像分類研究. 遙感技術(shù)與應(yīng)用,2010,25(4)： 540-546.[doi： 10.11873/j.issn.1004-0323.2010.4.540]

5陳文倩,丁建麗,王嬌,等. 基于高分一號(hào)影像的土地覆被分類方法初探. 干旱區(qū)地理,2016,39(1)： 182-189.

6中華人民共和國(guó)國(guó)家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局,中國(guó)國(guó)家標(biāo)準(zhǔn)化管理委員會(huì). GB/T 21010-2007 土地利用現(xiàn)狀分類.北京： 中國(guó)標(biāo)準(zhǔn)出版社,2007.

7胡玉福,鄧良基,匡先輝,等. 基于紋理特征的高分辨率遙感圖像土地利用分類研究. 地理與地理信息科學(xué),2011,27(5)： 42-45.

8高程程,惠曉威. 基于灰度共生矩陣的紋理特征提取. 計(jì)算機(jī)系統(tǒng)應(yīng)用,2010,19(6)： 195-198.

9賈永紅. 計(jì)算機(jī)圖像處理與分析. 武漢： 武漢大學(xué)出版社,2001： 231-235.

10Drǎgu? L,Tiede D,Levick SR. ESP： A tool to estimate scale parameter for multiresolution image segmentationof remotely sensed data. International Journal of Geographical Information Science,2010,24(6)： 859-871. [doi： 10.1080/13658810903174803]

11金文韜,許捍衛(wèi),李正源,等. 面向?qū)ο蟮母叻直媛视跋裢恋乩眯畔⑻崛? 地理空間信息,2016,14(9)： 104-106.

12張超,李智曉,李鵬山,等. 基于高分辨率遙感影像分類的城鎮(zhèn)土地利用規(guī)劃監(jiān)測(cè). 農(nóng)業(yè)機(jī)械學(xué)報(bào),2015,46(11)：323-329. [doi： 10.6041/j.issn.1000-1298.2015.11.044]

13Pontius Jr RG,Millones M. Death to Kappa： Birth of quantity disagreement and allocation disagreement for accuracy assessment. International Journal of Remote Sensing,2011,32(15)： 4407-4429. [doi： 10.1080/01431161.2011.552923]

14Nourqolipour R,Shariff ARBM,Balasundram ESK,et al. A GIS-based model to analyze the spatial and temporal development of oil palm land use in Kuala Langat district,Malaysia. Environmental Earth Sciences,2015,73(4)：1687-1700. [doi： 10.1007/s12665-014-3521-1]

15Alonso-Benito A,Arroyo LA,Arbelo M,et al. Pixel and object-based classification approaches for mapping forest fuel types in Tenerife Island from ASTER data. International Journal of Wildland Fire,2013,22(3)： 306-317. [doi： 10.1071/WF11068]

16Amiri F,Rahdari V,Najafabadi SM,et al. Multi-temporal landsat images based on eco-environmental change analysis in and around Chah Nimeh reservoir,Balochestan (Iran).Environmental Earth Sciences,2014,72(3)： 801-809. [doi：10.1007/s12665-013-3004-9]

17張俊,于慶國(guó),侯家槐. 面向?qū)ο蟮母叻直媛视跋穹诸惻c信息提取. 遙感技術(shù)與應(yīng)用,2010,25(1)： 112-117. [doi：10.11873/j.issn.1004-0323.2010.1.112]