趙 蓮,張錦水1,,胡潭高,陳聯(lián)裙,李 樂

(1.農(nóng)業(yè)部資源遙感與數(shù)字農(nóng)業(yè)重點開放實驗室,北京 100081;2.北京師范大學(xué)資源學(xué)院,北京 100875)

變端元混合像元分解冬小麥種植面積測量方法

趙 蓮2,張錦水1,2,胡潭高2,陳聯(lián)裙2,李 樂2

(1.農(nóng)業(yè)部資源遙感與數(shù)字農(nóng)業(yè)重點開放實驗室,北京 100081;2.北京師范大學(xué)資源學(xué)院,北京 100875)

針對線性混合像元分解 (Linear Spectral Unmixing,LSU)在端元 (Endmember)個數(shù)不變情況下常會出現(xiàn)端元分解過剩現(xiàn)象導(dǎo)致分解結(jié)果精度不高的問題,以地物分布的聚集性特征為基礎(chǔ),提出了基于格網(wǎng)的變端元線性混合像元分解 (Dynamic Endmember LSU,DELSU)方法。以冬小麥為研究目標(biāo),采用 Landsat TM圖像為實驗數(shù)據(jù)、高分QuickBird圖像目視解譯冬小麥結(jié)果為真值精度評價數(shù)據(jù),利用本文提出的 DELSU方法進行冬小麥提取。實驗結(jié)果表明:該方法比最大似然方法、LSU方法更能準(zhǔn)確地獲取冬小麥面積,在一定程度上吸收了傳統(tǒng)分類方法的優(yōu)點,提高了目標(biāo)地物的測量精度;同時作為一種改進的LSU方法也適用于其他土地利用 /覆蓋類型的測量。

DELSU;LSU;格網(wǎng);冬小麥;遙感

0 引言

利用遙感圖像進行土地利用 /覆蓋監(jiān)測中,混合像元現(xiàn)象不可避免,尤其是中分辨率遙感圖像中混合像元問題十分嚴(yán)重[1-3]。通過一定的方法計算出混合像元的典型地物組成比例,可解決混合像元問題,提高定性和定量遙感測量精度,在農(nóng)作物種植面積監(jiān)測中有著重要的應(yīng)用價值[4-5]?；旌舷裨纸?Spectral Unmixing,SU)是提取植被覆蓋度的主要方法[6],其測量結(jié)果為各種地物的豐度,其模型可歸結(jié)為 5類[7]:線性模型、概率模型、幾何光學(xué)模型、隨機幾何模型和模糊分析模型。

線性混合像元分解 (Linear Spectral Unmixing,LSU)模型因其簡單實用而被廣泛應(yīng)用[8-9],特別是在圖像波段數(shù)目較少、光譜分辨率不高的情況下[10]。在像元分解中,用傳統(tǒng)的LSU模型從圖像上選取所有端元進行分解,但實際上大多數(shù)圖像區(qū)域或混合像元只是由特定幾種端元組成。許多學(xué)者采用可變的端元模式[11-13]以提高測量精度,如Roberts等[14]提出了多端元混合像元分析 (Multiple End member Spectral Mixture Analysis,M ESMA)方法 ,將各種地物的光譜值組成一系列的候選模型 (Candidate Model),用每個候選模型分別對圖像進行分解,通過比較均方根誤差 (RM S)來確定何種模型入選,有效地提高了分解精度;但在全區(qū)范圍內(nèi)選擇候選模型,采用窮舉方法確定端元,計算量較大,影響其應(yīng)用效果。叢浩等[10]提出了基于混合像元的光譜響應(yīng)特征和地物分布的集聚性來實現(xiàn)端元可變的方法,可以保證光譜矢量最接近的端元入選并參與運算,并考慮到地物分布的聚集性,計算以像元為中心的 3像元 ×3像元模板內(nèi)的像元的端元加權(quán)和,通過設(shè)定閾值確定入選端元;但由于遙感圖像光譜的不確定性,3像元 ×3像元模板并不一定適合。

本文針對 LSU的不足,提出變端元線性混合像元分解 (Dynamic Endmember LSU,DELSU)方法。以冬小麥為研究目標(biāo),根據(jù)地物分布具有一定空間聚集性的特征,采用基于格網(wǎng)的 DELSU,以格網(wǎng)為單元利用局部端元類型和數(shù)目進行變端元方式的混合像元分解,消除混合像元分解端元數(shù)目“過?！痹斐傻挠绊?分析格網(wǎng)大小對測量精度的影響,以實現(xiàn)在最優(yōu)格網(wǎng)尺寸下進行 DELSU,提高測量精度。

1 變端元線性混合像元分解模型

1.1 LSU模型

LSU模型是混合像元分解常用方法,該模型假設(shè)像元的輻射值 (或反射率)可以表示成像元內(nèi)各端元的光譜特征及其所占面積百分比的線性函數(shù)[15-16],模型定義為

式中,Rb為被處理像元在第 b波段的反射率;N為端元數(shù)目;fi為端元 i的面積組成百分比;Rib為端元 i在第 b波段的反射率;eb為該像元在第 b波段的分解殘差。其中

式中,S是 m行 ×n列的矩陣,列是 n個端元的光譜矢量;f是 n行 ×1列的矩陣,列是 n個端元在混合像元 R中的百分含量?？紤]到 fi非負(fù)約束的復(fù)雜性,本文采用簡化方法將 fi中的負(fù)值歸零并對其余 fi簡單地線性壓縮使其重新歸一化。

1.2 DELSU模型

在LSU中常遇見的情況是:整幅圖像中存在若干類型的端元,但作為某一個混合像元,其端元的組成類型是特定的。LSU采用圖像上所有端元進行混合像元分解,即 S不變,這容易造成端元分解過剩,導(dǎo)致分解精度降低。在圖 1中,分別標(biāo)示出冬小麥、水體、建筑用地和裸地 4種地物類型?；旌舷裨啻嬖谟诓煌匚锔采w類型的邊界處 (如圖 1中 A、B、C和 D處)。LSU在分解圖 1的像元時,采用全區(qū)端元進行分解,如像元 A僅由建筑用地和冬小麥組成,而 LSU卻要采用 4種端元分解,于是就產(chǎn)生了端元過剩的問題。

圖 1 變端元混合像元分解方法示意圖Fig.1 Sketch map of DELSU method

為此,本文提出 DELSU模型,采用格網(wǎng)的形式將圖像分成若干個子區(qū),在每個子區(qū)內(nèi)確定端元的類型和數(shù)目 (根據(jù)參考圖像的分類結(jié)果確定)。分解原則為:若格網(wǎng)內(nèi)的分類結(jié)果中僅含一種地物類型則不進行線性分解,直接賦值為 1,其他類型賦 0;若格網(wǎng)內(nèi)分類結(jié)果類型為 n種,則確定端元數(shù)目為n,用其含有的端元類型組成端元矩陣 SN,并用 SN對混合像元進行線性分解 (N為圖像上總的地物類型個數(shù))。于是式 (2)變?yōu)?/p>

這樣將能有效地去除不相關(guān)端元,如圖 1中像元 A所處的 3×3格網(wǎng)中只含有冬小麥和建筑用地2種地物類型,因而像元 A用冬小麥和建筑用地 2種端元組成的矩陣 S2進行線性分解;同理,像元 B直接賦值為冬小麥;像元 C用冬小麥、裸地和建筑用地 3種端元組成的矩陣 S3進行分解;像元 D則用冬小麥、裸地、建筑用地和水體 4種端元組成的矩陣 S4進行分解。

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

研究區(qū)位于北京市通州區(qū)東北角與河北省的交界處。該區(qū)地形平坦,主要地物類型為冬小麥、裸地、建筑用地和水體。該區(qū)農(nóng)業(yè)發(fā)達,冬小麥一般在4月份進入抽穗期,其光譜特征明顯地從土壤背景中突顯出來,有利于進行冬小麥種植面積遙感估測。

本研究從 123/032的 Landsat TM圖像 (獲取時間為 2006年 4月 7號,覆蓋范圍為 10 km ×10 km,空間分辨率為 30m,含有 6個波段)上選取質(zhì)量較好、無云的子區(qū)作為研究區(qū) (圖 2(左))。

利用與研究區(qū)范圍相同、時相相近的QuickBird多光譜圖像 (獲取時間為 2006年 5月 2號,空間分辨率為 2.4 m)驗證本文方法的實驗結(jié)果 (圖 2(右))。由圖 2可見,4種地物特征明顯,能夠清晰獲得冬小麥的分布信息。

3 實驗

3.1 流程設(shè)計

本文實驗流程主要包括圖像處理、端元選取、監(jiān)督分類、DELSU和結(jié)果分析 (圖 3)。

圖 3 實驗技術(shù)流程Fig.3 Flow char t of the test

3.2 數(shù)據(jù)預(yù)處理

首先,將遙感器獲得的輻射亮度 DN值轉(zhuǎn)換為反射率值[17,18],得到反射率值的 TM圖像,并對其進行大氣輻射校正和幾何糾正。將處理后的 TM圖像與QuickB ird圖像進行幾何配準(zhǔn),誤差控制在一個像元之內(nèi)。對QuickBird圖像進行人機交互解譯,并將解譯結(jié)果數(shù)字化,得到研究區(qū)冬小麥分布的矢量地塊數(shù)據(jù);將矢量地塊數(shù)據(jù)轉(zhuǎn)化成與 TM數(shù)據(jù)相同分辨率的豐度圖,用作真值對最大似然分類、LSU和DELSU方法進行精度評價。

3.3 端元選擇

端元選擇是本研究的關(guān)鍵所在,包括端元值和端元類型的確定。端元值通過在整個圖像中選擇純凈端元并計算其各波段的反射率得到;端元類型則通過監(jiān)督分類、根據(jù)格網(wǎng)內(nèi)統(tǒng)計分類結(jié)果類型確定。

3.3.1 端元值的確定

最小噪聲分離變換 (Minimum Noise Fraction,MNF)技術(shù)是對高光譜遙感圖像進行降維的有效方法,目的是將噪音與有用信息分離[19]。利用 MNF變換對 TM圖像進行去噪處理;對圖像經(jīng)MNF變換后生成的二維散點圖 (圖 4(a))進行分析,確定最終要選擇的端元。選取散點圖犄角處的點,在 TM影像上分別為冬小麥、建筑用地、亮裸地和水體。經(jīng)過對實驗數(shù)據(jù)所覆蓋地物類型 (圖 2(右))的光譜特征分析,原始圖像的亮裸地覆蓋范圍較小,在圖像上手動選擇較暗處的裸地作為端元之一,并命名為裸地;最終選擇冬小麥、建筑用地、裸地和水體 4種典型地物作為端元類型[20],獲得 TM的 6個波段圖像上每種端元組分的反射率 (圖 4(b))。

圖 4 端元光譜反射率Fig.4 Spectral reflectance of the end members

由圖 4(b)可以看出,冬小麥與其他 3種端元的反射率分離性較好,可供較高精度地提取冬小麥的種植面積;建筑用地和裸地在 5、7波段的反射率值有交叉,這會在一定程度上造成 2類地物的提取誤差,但均與冬小麥有較大的分離性,對冬小麥提取豐度的影響不大,故上述端元選擇較合理。

3.3.2 端元類型的確定

端元類型的確定是通過在格網(wǎng)內(nèi)統(tǒng)計最大似然分類結(jié)果得到的,監(jiān)督分類結(jié)果對端元選擇的正確性有很大影響,因此必須將監(jiān)督分類結(jié)果與真值數(shù)據(jù)進行比較,對分類精度進行評估。

針對分類結(jié)果,采用隨機評價方法,用分層隨機抽樣 (Stratified Random)方法在研究區(qū)布置 256個隨機點[21]。以 QuickBird數(shù)據(jù)為參考,確定每個樣本點的真實地類 (分類誤差矩陣見表 1)。

表 1 最大似然分類誤差矩陣Tab.1 Error matrix of the maximum likelihood classification

由表 1可以看出,研究目標(biāo)冬小麥的分類精度達到 91%(70/78),比較理想;建筑用地和裸地的反射率較接近而導(dǎo)致這 2類地物的錯分率偏大,但是總分類精度達到了 88%,精度較高?？傊?用最大似然分類的結(jié)果作為窗口內(nèi)端元類型選擇的依據(jù)是可行的。

3.4 最優(yōu)格網(wǎng)尺寸分析

DELSU模式是在每個格網(wǎng)內(nèi)確定端元的類型和數(shù)目,顯然,格網(wǎng)大小對保證端元類型和數(shù)目的正確性有重要的作用。

由于通過地面實際測量得到混合像元的比例非常困難[22],本文采用 QuickB ird圖像分類結(jié)果作為真值與 TM圖像分類結(jié)果進行比較。采用均方根誤差 (RM SE)作為精度評價指標(biāo),其定義為[23]

式中,x(i)為 TM圖像單個像元中冬小麥所占比例;y(i)為與 TM圖像對應(yīng)分辨率的QuickB ird圖像單個像元中冬小麥所占比例;n為 TM圖像中像元的總數(shù)。RM SE值越小,表示與真值越接近,說明混合像元分解精度越高。

為有效地進行精度評價,將評價區(qū)域劃分為 4種情況:全區(qū)、冬小麥區(qū)、冬小麥純凈區(qū)和冬小麥混合區(qū)。其中,全區(qū)為整個研究區(qū);其余 3種情況是以QuickBird影像目視解譯得到的冬小麥豐度圖為參考標(biāo)準(zhǔn),定義QuickBird提取的冬小麥豐度圖中不等于0的地區(qū)為冬小麥區(qū),大于 0小于 100%的地區(qū)為冬小麥混合區(qū),等于 100%的地區(qū)為冬小麥純凈區(qū)。

鑒于本文研究目標(biāo)是冬小麥,由于地物在空間分布上表現(xiàn)出一定的聚集性以及遙感本身的不確定性,分別選取 3像元 ×3像元、5像元 ×5像元、7像元 ×7像元、9像元 ×9像元和 11像元 ×11像元大小的格網(wǎng)劃分研究區(qū);當(dāng)格網(wǎng)尺寸大于 11像元 ×11像元時,多數(shù)格網(wǎng)內(nèi)端元數(shù)目為 4種,格網(wǎng)過大設(shè)置就失去了意義。采用以上幾種格網(wǎng)在全區(qū)、冬小麥區(qū)、冬小麥純凈區(qū)和冬小麥混合區(qū)進行冬小麥面積提取,其 RM SE如表 2所示。

表 2 不同尺寸格網(wǎng)的 RMSEFig.2 RM SE computed in serial grid scales

圖 5(a)～(d)分別為 QuickBird圖像及其 3像元 ×3像元、5像元 ×5像元、7像元 ×7像元格網(wǎng)分類結(jié)果的局部放大圖。

圖 5 冬小麥空間分布局部放大圖Fig.5 Spatial distribution of the winter wheat

分析表 2和圖 5可以看出:①冬小麥純凈區(qū)的RM SE值隨格網(wǎng)尺寸的變大逐漸升高。這是由于格網(wǎng)變大后,格網(wǎng)內(nèi)部的端元類型變多,純凈的像元被分解,使精度降低 (圖 5(a)的 A處明顯為純凈的冬小麥區(qū),而 5(d)的 A處卻被分解);②冬小麥混合區(qū)的 RM SE值與格網(wǎng)尺寸沒有明顯的比例關(guān)系,這是由于網(wǎng)格過大或過小均會對混合區(qū)分類精度造成影響;其中 3像元 ×3像元格網(wǎng)的 RM SE值最高,原因是格網(wǎng)選擇過小,最大似然分類結(jié)果在格網(wǎng)內(nèi)部以含有一種地物的情況居多,導(dǎo)致混合像元不被分解,而被當(dāng)作純凈像元處理 (如圖 5(a)的 B處明顯為道路,而圖 5(b)的 B處卻被分類為冬小麥)。由此可知,格網(wǎng)尺寸過大會導(dǎo)致端元數(shù)目選擇“過?！?而尺寸過小又會出現(xiàn)混合像元不被分解的情況,因此選擇最優(yōu)格網(wǎng)尺寸十分必要。

綜上所述,以 5像元 ×5像元格網(wǎng)劃分所選研究區(qū)域,全區(qū)和冬小麥區(qū)的 RM SE值均為最小,冬小麥純凈區(qū)和混合區(qū)的誤差都相對合理 (如圖 5(c))。3像元 ×3像元格網(wǎng)雖然在冬小麥純凈區(qū)RM SE值比較小,但全區(qū)和冬小麥混合區(qū) RM SE值均偏大。格網(wǎng)大于 5像元 ×5像元以后,由于端元過剩,RMSE值均有增大的趨勢;但冬小麥混合區(qū)精度略微變高,這是由于研究區(qū)中的冬小麥混合區(qū)多集中于地塊邊界像元比例相對較低的道路處,監(jiān)督分類時此部分往往被誤分為冬小麥,導(dǎo)致在小尺寸格網(wǎng)內(nèi)選擇端元時出現(xiàn)誤選現(xiàn)象。因此,本文確定最優(yōu)格網(wǎng)尺寸為 5像元 ×5像元。

4 結(jié)果分析

4.1 精度分析

以 5像元 ×5像元為格網(wǎng),采用以下 3種方法對研究區(qū) TM影像進行分類:①最大似然方法;②LSU;③DELSU。分別在全區(qū)、冬小麥區(qū),冬小麥純凈區(qū)和冬小麥混合區(qū)進行冬小麥面積提取并進行不同分類方法的精度比較 (表 3)。

表 3 不同分類方法的 RMSE值Tab.3 RM SE computed in serial classification methods

從表 3可以看出:①DELSU在全區(qū)規(guī)模上的測量精度高于最大似然分類和 LSU;②DELSU在冬小麥區(qū)的測量精度高于最大似然分類和 LSU;③DELSU在冬小麥純凈區(qū)的測量精度明顯高于 LSU;④DELSU在冬小麥混合區(qū)的測量精度優(yōu)于最大似然分類方法;⑤DELSU在冬小麥混合區(qū)的測量精度要略低于LSU,主要原因是監(jiān)督分類的誤分現(xiàn)象而導(dǎo)致在格網(wǎng)內(nèi)選擇端元類型和數(shù)目時出現(xiàn)誤選的現(xiàn)象。但是在端元選擇正確地方 (如圖 5(c))的冬小麥面積提取精度是要比 LSU精度高的。故用此方法進行冬小麥面積測量會取得比較好的結(jié)果。

4.2 空間分布特征分析

下面從空間分布特征方面,評價 LSU和 DELSU分類結(jié)果 (圖 6、7)。

圖 6 LSU分類結(jié)果Fig.6 The result of LSU classification

圖 7 DELSU分類結(jié)果Fig.7 The result of DELSU classification

從圖 6、7中可以清晰地看出,DELSU的分類結(jié)果圖像上不含該類端元的區(qū)域值為 0,但是 LSU的分類結(jié)果幾乎覆蓋整個圖像,明顯受到端元過剩的影響。圖 6(b)表現(xiàn)尤為明顯,在 QuickBird圖像 (圖2(右))上可看出,水體的覆蓋范圍很小,但是在LSU的分類結(jié)果圖像上卻存在大范圍不為 0的區(qū)域,這顯然與實際情況不符。另外,對于組成分圖像上近似認(rèn)為是純凈像元的區(qū)域 (比如純凈的冬小麥和水體部分)在圖 7上的亮度明顯高于圖 6,而邊緣帶混合像元被分解,這正好與圖像特征相符合。

圖 8為冬小麥區(qū)局部放大圖。從圖 8中的A區(qū)可以看出,最大似然分類方法擴大了冬小麥的光譜信息,道路被誤分為冬小麥,擴大了冬小麥面積的提取結(jié)果;而DELSU和 LSU進行混合像元分解可以體現(xiàn)道路的光譜特征。B區(qū)為長勢好且大面積分布的冬小麥生長區(qū),從 QuickB ird圖像看出,該區(qū)冬小麥豐度為 100%;LSU對所有像元進行無區(qū)別分解,在該區(qū)得到冬小麥豐度在 40%～80%,顯然與實際情況不符;DELSU在局部確定端元,因其內(nèi)部只有一種端元,該區(qū)提取結(jié)果為 100%,與實際相符,此方法在該區(qū)域吸收了大片種植特征,從而在一定程度上抵制了光譜不穩(wěn)定性對混合像元分解造成的影響,分解結(jié)果與最大似然分類方法相同。C區(qū)為冬小麥零星分布區(qū)域,直接計算 C區(qū)冬小麥真值平均豐度為 0.10,DELSU對此區(qū)的提取豐度是 0.12,LSU為 0.06,顯然 DELSU精度高于 LSU,這是因為此區(qū)的最大似然分類結(jié)果中只有冬小麥和裸地 2種地物,DELSU用這 2種端元進行線性分解,有效地去除了不相干端元的影響,提高了冬小麥提取精度。

圖 8 不同方法提取的冬小麥區(qū)分布對比圖 (局部放大)Fig.8 Spatial distribution of the winter wheat extraced by different methods

由空間分布特征分析可以看出,DELSU在混合區(qū)的冬小麥提取精度大于最大似然分類,在純凈區(qū)的提取精度優(yōu)于 LSU,在端元提取準(zhǔn)確的情況下混合區(qū)提取精度也優(yōu)于LSU分類結(jié)果。因此,DELSU兼顧了傳統(tǒng)分類方法的優(yōu)點,既能夠使冬小麥純凈區(qū)在一定程度上避免了光譜不穩(wěn)定的影響,又能夠在冬小麥混合區(qū)利用光譜的混合特性進行分解,達到了提高冬小麥面積測量精度的目的。

5 結(jié)論

(1)本文提出的 DELSU,利用端元在空間分布上具有一定的形狀和聚集性的特點,動態(tài)調(diào)整端元矩陣,有效地去除了不相干端元,提高了冬小麥種植面積測量的精度。該方法不但適合冬小麥,對其他地物類型的提取也有一定的借鑒意義。

(2)從全區(qū)覆蓋來看,DELSU優(yōu)化了冬小麥面積測量結(jié)果。與冬小麥真值豐度圖相比,DELSU的RM SE=0.20,優(yōu)于最大似然分類 (RM SE=0.25)和LSU(RM SE=0.27)。

(3)DELSU綜合了傳統(tǒng)分類方法的優(yōu)點,在不同的區(qū)域采用不同的分類方法,提高了測量精度。DELUS在冬小麥混合區(qū)的提取精度 (RM SE=0.36)優(yōu)于最大似然分類 (RM SE=0.50);在冬小麥純凈區(qū)的提取精度 (RM SE=0.25)優(yōu)于 LSU(RM SE=0.40);在端元選擇準(zhǔn)確的情況下對冬小麥混合區(qū)的提取精度也優(yōu)于LSU;DELSU對冬小麥區(qū)的提取精度(RM SE=0.31)優(yōu)于最大似然分類 (RM SE=0.38)和 LSU(RM SE=0.37)。

(4)本文方法僅是利用變端元進行混合像元分解的一種嘗試,仍存在一些問題有待進一步深入探討:①本文只解決了端元類型和數(shù)目的選擇問題,但是全區(qū)共用一套平均端元光譜值進行線性分解,還不能有效地解決格網(wǎng)內(nèi)部的“同物異譜”問題;②最大似然分類的結(jié)果作為窗口內(nèi)端元選擇的依據(jù),其測量精度受監(jiān)督分類結(jié)果的影響較大;③以格網(wǎng)為單元確定端元類型和數(shù)目的方法會導(dǎo)致各網(wǎng)格結(jié)果的不連續(xù)問題,今后的研究中可考慮從研究像元為中心的周圍像元中獲得端元類型和數(shù)目的方法來解決此類問題。

[1] 胡健波,陳 瑋,李小玉,等.基于線性混合像元分解的沈陽市三環(huán)內(nèi)城市植被覆蓋變化[J].應(yīng)用生態(tài)學(xué)報,2009,20(5):1140-1146.

[2] 許文波,張國平,范錦龍,等.利用MOD IS遙感數(shù)據(jù)檢測冬小麥種植面積[J].農(nóng)業(yè)工程學(xué)報,2007,23(12):144-149.

[3] 李 霞,王 飛,徐德斌,等.基于混合像元分解提取大豆種植面積的應(yīng)用探討[J].農(nóng)業(yè)工程學(xué)報,2008,24(1):213-217.

[4] 鄒金秋,陳佑啟,Uchida S,等.利用 Terra/MOD IS數(shù)據(jù)提取冬小麥面積及精度分析[J].農(nóng)業(yè)工程學(xué)報,2007,23(11):195-200.

[5] 武永利,王云峰,張建新,等.應(yīng)用線性混合模型遙感監(jiān)測遙感監(jiān)測冬小麥種植面積[J].農(nóng)業(yè)工程學(xué)報,2009,25(2):136-140.

[6] 田 靜,閻 雨,陳圣波.植被覆蓋率的遙感研究進展[J].國土資源遙感,2004(1):1-5.

[7] Ichoku C,Kamieli A.A Review of Modeling Technique for Subpixel L land Cover Estimation[J].Remote Sensing Reviews,1996,13(3):161-186.

[8] Chabrillat S,Pinet PC,Ceuleneer G,et al.Ronda Peridotite Massif:Methodology for Its Geological Mapping and Lithological Discrimination from Airborne Hyper-spectral Data[J].International Journal of Remote Sensing,2000,21:2363-2388.

[9] 李 惠,陳健飛,余 明.線性光譜混合模型的 ASTER影像植被應(yīng)用分析[J].地球信息科學(xué),2005,7(1):103-106.

[10]叢 浩,張良培,李平湘.一種端元可變的混合像元分解方法[J].中國圖象圖形學(xué)報,2006,11(8):1092-1096.

[11]Asner G P,Heidebrecht KB.Spectral Unmixing of Vegetation,Soil and Dry Carbon Cover in A rid Regions:Comparing Multispectral and Hyper-spectral Observations[J].International Journal of Remote Sensing,2002,23:3939-3958.

[12]Song C H.Spectral Mixture Analysis for Subpixel Vegetation Fractions in the Urban Environment:How To Incorporate Endmember Variability?[J].Remote Sensing of Environment,2005,95:248-263.

[13]Bates on C A,Asner G P,Wessman C A.Endmember Bundles:A New Approach to Incorporating Endmember Variability into Spectral Mixture Analysis[J].IEEE Transactions on Geoscience and Remote Sensing,2000,38:1083-1094.

[14]Roberts D A,Gardner M,Church R,et al.Mapping Chaparral in the Santa Monica Mountains Using Multiple Endmember Spectral Mixture Models.Remote Sensing of Environment[J].1998,65:267-279.

[15]Joseph M,Piwowar.Temporal Mixture Analysis of Arcticsea Ice Imagery:A New Approach for Monitoring Environmental Change.Remote Sensing of Environment[J].1998,63:195-207.

[16]Zhou X T,Folving S.Application of Spectral Mixture Modeling to the Regional Assessment and Degradation:A Case Study from Basilicata Italy[J].Land Degradation&Rehabilitation,1994,5:215-222.

[17]劉建貴,吳長山,張 兵,等.PH I成像光譜圖像反射率轉(zhuǎn)換[J].遙感學(xué)報,1999,3(4):290-294.

[18]韋玉春,黃家柱.Landsat5圖像的增益、偏置取值及其對行星反射率計算分析[J].地球信息科學(xué),2006(1):110-113.

[19]楊可明,陳云浩,郭達志,等.基于端元提取的高光譜影像特定目標(biāo)識別[J].金屬礦山,2006(6):48-52.

[20]Ridd M K.ExploringV-I-S(Vegetation-Impervious Surface-Soil)Model for Urban Ecosystem Analysis through Remote Sensing:Comparative Anatomy for Cities[J].International Journal of Remote sensing,1995,16:2165-2185.

[21]劉詠梅,楊勤科,湯國安.陜北黃土丘陵地區(qū)坡耕地遙感分類方法研究[J].水土保持通報,2004,24(4):51-54.

[22]吳 波,張良培,李平湘.基于支撐向量機概率輸出的高光譜影像混合像元分解[J].武漢大學(xué)學(xué)報 (信息科學(xué)版),2006(1):51-54.

[23]胡潭高,張錦水,賈 斌,等.不同分辨率遙感圖像混合像元線性分解方法研究[J].地理與地理信息科學(xué),2008,24(3):20-24.

(責(zé)任編輯:劉心季)

The Application of the Dynamic Endmember Linear Spectral Unmixing Model to Winter Wheat Area Estimation

ZHAO Lian2,ZHANG Jin-shui1,2,HU Tan-gao2,CHEN Lian-qun2,L ILe2
(1.State Key Laboratory of Earth Surface Processes and Resource Eco logy,Beijing 100081,China;2.College of Resources Science and Technology,Beijing Normal University,Beijing 100875,Ch ina)

Linear spectral unmixing(LSU)is the most common method for solving mixed pixel problem;nevertheless,if the number of the pixels’end member is regarded as unchangeable,the traditional pixel unmixing algorithm cannotattain a good result.In the light of the characteristic that pixelsusually have the same composition as their neighboring pixels,the authors proposed a grid-based dynamic end member linear spectral unmixing(DELSU)model.The land cover classification experiment was conducted with the adoption of the Landsat TM im age as the experimental data.The abundance map of winter wheat derived from the visual interpretation of the Quick Bird image was used for accuracy evaluation.The experimental results show that the use of the DELSU model could extract the area of winter wheat at a precision higher than that of the traditional maximum likelihood method and the LSU model.This model absorbs the traditional classification advantages and imp roves the measurement accuracy of the target features.A s an imp roved method of LSU,DELSU is also applicable to the measurement of other land use/cover types.

DELSU;LSU;Grid;Winter wheat;Remote sensing

趙 蓮 (1987-),女,地圖制圖學(xué)與地理信息工程專業(yè)在讀碩士研究生,主要從事統(tǒng)計遙感方面的研究。

TP 751.1

A

1001-070X(2011)01-0066-07

2010-04-15;

2010-06-17

農(nóng)業(yè)部資源遙感與數(shù)字農(nóng)業(yè)重點實驗室開放基金項目 (編號:RDA 0807)和國家高技術(shù)研究發(fā)展計劃資助項目 (編號:2006AA 120103、2006AA 120101)共同資助。