劉 斌，史 云，吳文斌，段玉林，趙立成

(中國(guó)農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所/農(nóng)業(yè)農(nóng)村部信息技術(shù)重點(diǎn)實(shí)驗(yàn)室，北京 100081)

0 引言

無(wú)人機(jī)遙感具備高分辨率、操作簡(jiǎn)單、獲取數(shù)據(jù)快及低成本等特點(diǎn)，可以快速針對(duì)某一區(qū)域進(jìn)行影像采集工作，獲取更精確的作物分布信息，成為衛(wèi)星遙感和航空遙感的重要補(bǔ)充，對(duì)作物監(jiān)測(cè)技術(shù)的發(fā)展和應(yīng)用具有重大意義[1-2]。

準(zhǔn)確的作物分類仍然是一項(xiàng)具有挑戰(zhàn)性的任務(wù)，有效進(jìn)行多種特征變量結(jié)合并使用是提高農(nóng)作物遙感分類精度的關(guān)鍵。目前已有學(xué)者基于無(wú)人機(jī)影像針對(duì)農(nóng)作物分類信息提取作了大量的研究，并提出相關(guān)的技術(shù)和方法。利用無(wú)人機(jī)影像進(jìn)行農(nóng)作物分類方法主要分為光譜特征及紋理特征兩類?；诠庾V特征的農(nóng)作物分類方法主要是通過(guò)計(jì)算各類別光譜統(tǒng)計(jì)量，構(gòu)造植被差異指數(shù)，并在這些植被差異指數(shù)的基礎(chǔ)上通過(guò)統(tǒng)計(jì)學(xué)習(xí)實(shí)現(xiàn)農(nóng)作物遙感分類。其中，汪小欽等[3]提出了一種新的可見(jiàn)光波段差異植被指數(shù)，實(shí)現(xiàn)了健康綠色植被信息的提取。李鑫[4]提出了基于紅綠波段的增強(qiáng)型紅綠差值指數(shù)，對(duì)藍(lán)藻信息識(shí)別可信度高。針對(duì)多光譜數(shù)據(jù)，劉偉等[5]結(jié)合NDVI、NDWI、及均值特征確定了最佳特征組合，改善了裸地與沙石的區(qū)分。Pena等[6]使用六波段多光譜影像，實(shí)現(xiàn)了在玉米生長(zhǎng)早期的雜草制圖。Doi[7]對(duì)多光譜影像進(jìn)行彩色合成，增加了相似像元的可區(qū)分性。由于單純依賴多光譜特征的農(nóng)作物遙感分類存在同物異譜及同譜異物的現(xiàn)象，此類研究多適用于分類對(duì)象單一，分類特征針對(duì)特定對(duì)象的情形。

另外，通過(guò)影像濾波提取紋理特征，可以實(shí)現(xiàn)在遙感影像中對(duì)特定作物的識(shí)別或分類。如李蘊(yùn)雅等[8]建立相應(yīng)的紋理規(guī)則集，對(duì)無(wú)人機(jī)影像進(jìn)行樹(shù)冠的提取。郭鵬[9]使用可見(jiàn)光影像，選取亮度、飽和度和紅色二階矩作為最優(yōu)分類特征對(duì)農(nóng)田作物進(jìn)行分類，分類精度明顯高于顏色指數(shù)方法。李宗南等[10]選取最適宜區(qū)分正常和倒伏玉米的特征為灰度共生矩陣的紅、綠、藍(lán)色均值特征。通過(guò)比較基于色彩特征或紋理特征的倒伏玉米面積提取結(jié)果，基于紋理特征分類結(jié)果更準(zhǔn)確，提取誤差顯著低于基于色彩特征提取方法。但是，農(nóng)作物紋理特征隨著農(nóng)作物、遙感影像空間分辨率及波段的變化而變化，因此在復(fù)雜種植結(jié)構(gòu)區(qū)域影像中存在紋理特征確定困難等問(wèn)題。

此外，不同農(nóng)作物間的冠層結(jié)構(gòu)差異非常明顯，作物理性狀也可用于分類，如高度、形態(tài)、葉傾角等。但傳統(tǒng)遙感無(wú)法獲得高分辨率冠層結(jié)構(gòu)數(shù)據(jù)，隨著無(wú)人機(jī)遙感及攝影測(cè)量技術(shù)的快速發(fā)展，從而能夠獲取更多類型的數(shù)據(jù)，為農(nóng)作物遙感分類提供了新的發(fā)展空間，因此如何從無(wú)人機(jī)影像中提取作物物理性狀成為新的研究熱點(diǎn)[11]。其中，Kim等[12]通過(guò)手動(dòng)采集地面控制點(diǎn)生成TIN(Triangulated Irregular Network)模型進(jìn)而生成DEM(Digital Elevation Model)數(shù)據(jù)，計(jì)算了DSM和DEM的差值nDSM(normalized DSM)，選擇RGB波段、nDSM和改進(jìn)NDVI為特征，實(shí)現(xiàn)了地表覆蓋類型分類。楊琦等[13]采集了蔗糖全生育期的高清數(shù)碼影像，建立各生育期作物表面模型CSMs(crop surface models)，并提取株高，該方法表明CSMs提取的株高擁有較高的精度。Bendig等[14]證實(shí)了從CSMs提取的株高具有良好的精度，并建立了大麥株高與生物量的估算模型。研究表明，利用無(wú)人機(jī)獲取的數(shù)字表面模型可以提取地表信息，如作物株高信息等。在農(nóng)田區(qū)域下，提取分類特征nDSM需要DEM，控制點(diǎn)布控及量測(cè)較難，難以生成高精度DEM。同時(shí)，目前利用DSM數(shù)據(jù)提取分類特征的研究較少。

基于此，文章的核心思路是圍繞如何利用作物在光譜和空間維度上的聯(lián)合特征尤其是作物高程特征以提高農(nóng)作物分類精度。在該研究中，通過(guò)無(wú)人機(jī)搭載可見(jiàn)光相機(jī)獲得高分辨率圖像和數(shù)字表面模型，進(jìn)而利用兩個(gè)不同時(shí)期作物的DSM數(shù)據(jù)，生成差異數(shù)字表面模型(DDSM)，突出作物生長(zhǎng)差異特征，并將該特征引入農(nóng)作物分類中。同時(shí)，分析并提取農(nóng)作物在光譜及空間維度上的聯(lián)合特征，通過(guò)比較分析變異系數(shù)和類間差異系數(shù)，構(gòu)建了基于色彩、紋理、高度特征組合的SVM分類模型用于農(nóng)作物精細(xì)分類。

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

研究區(qū)位于黑龍江省農(nóng)業(yè)科學(xué)院民主研究基地，基地面積約為7.5km2，地形有較大起伏，測(cè)區(qū)內(nèi)包含水稻、玉米、大豆及馬鈴薯等多類農(nóng)作物。測(cè)區(qū)面積適中，作物類型多樣，種植結(jié)構(gòu)復(fù)雜度高，是檢驗(yàn)作物分類算法的理想實(shí)驗(yàn)區(qū)域。該研究中所涉及的實(shí)驗(yàn)數(shù)據(jù)分兩期在該區(qū)域進(jìn)行采集，分別是2017年8月及2018年6月。圖像獲取設(shè)備為SONY單反相機(jī)，選用固定翼無(wú)人機(jī)作為相機(jī)的搭載平臺(tái)。無(wú)人機(jī)獲取的影像為RGB圖像，且在無(wú)人機(jī)航拍過(guò)程中，旁向、航向上圖像的重疊率均達(dá)到80%，滿足生成RGB正射影像的要求。采用Smart3D軟件對(duì)獲取的圖像進(jìn)行處理，輸出實(shí)驗(yàn)對(duì)象區(qū)域的RGB正射影像和數(shù)字地表模型。實(shí)驗(yàn)基地如圖1。

圖1 研究區(qū)域

圖2 原始RGB影像及DSM

2 基于無(wú)人機(jī)影像的農(nóng)作物分類方法

2.1 RGB正射影像及數(shù)字地表模型獲取

該研究中研究對(duì)象區(qū)域的RGB正射影像及數(shù)字地表模型見(jiàn)圖2。其像大小為6 000×6 000像素，地面分辨率為0.1m。其中，在2018年6月采集實(shí)驗(yàn)數(shù)據(jù)時(shí)，試驗(yàn)區(qū)域作物尚處在作物生長(zhǎng)初期，苗高較低，可視為實(shí)驗(yàn)區(qū)域的背景數(shù)據(jù)。

2.2 無(wú)人機(jī)影像農(nóng)作物分類方法

圖3 技術(shù)路線

表1 農(nóng)田地物樣本數(shù)

類別訓(xùn)練樣本(個(gè))測(cè)試樣本(個(gè))池塘300100水稻200100玉米300200大豆400150樹(shù)300100建筑250100道路400150馬鈴薯150100亞麻400250草地300100裸地200100合計(jì)3 2001 450

由于無(wú)人機(jī)遙感可見(jiàn)光影像僅含有紅、綠、藍(lán)3通道的灰度信息，難以實(shí)現(xiàn)對(duì)農(nóng)作物類別的區(qū)分，因此該研究首先通過(guò)紋理濾波獲得研究區(qū)域農(nóng)作物紋理信息； 然后利用兩個(gè)不同時(shí)期作物的DSM數(shù)據(jù)，生成差異數(shù)字表面模型(DDSM，Difference of DSM)，分析并提取農(nóng)作物在光譜及空間維度上的聯(lián)合特征，通過(guò)比較分析變異系數(shù)和類間差異系數(shù)，構(gòu)建了彩色特征、紋理特征、高度特征組合； 分別基于的組合訓(xùn)練SVM分類器，進(jìn)行農(nóng)作物分類； 以地面真值結(jié)合目視解譯所得到的樣點(diǎn)農(nóng)作物類別的實(shí)測(cè)值，以檢驗(yàn)分類結(jié)果的精度。技術(shù)路線如圖3所示。

2.3 分類樣本選擇

根據(jù)目視分析和地面實(shí)際調(diào)查，確定研究區(qū)域內(nèi)地物分為11類，包括池塘、水稻、玉米、大豆、樹(shù)、建筑、道路、馬鈴薯、亞麻、草地和裸地。根據(jù)地物的面積比例大小選取樣本點(diǎn)數(shù)，如表1所示。

2.4 特征分析及選取

該文通過(guò)提取無(wú)人機(jī)影像中的光譜、紋理、高度等特征，分析并篩選適于農(nóng)作物分類的特征。濾波處理可以獲得農(nóng)作物的紋理特征，有助于解決同譜異物和同物異譜的問(wèn)題，進(jìn)而提高農(nóng)作物的分類精度。該研究中濾波處理在ENVI5.3環(huán)境下進(jìn)行，濾波窗口大小為27×27。通過(guò)濾波處理，可得到相應(yīng)波段的均值、方差、協(xié)同性、對(duì)比度、相異性、信息熵、二階矩和相關(guān)性等特征信息。

進(jìn)一步地，利用兩個(gè)不同時(shí)期作物的DSM數(shù)據(jù)，生成差異數(shù)字表面模型(DDSM)，突出作物生長(zhǎng)差異特征，并將該特征引入農(nóng)作物分類中。在該研究中，在2018年6月采集實(shí)驗(yàn)數(shù)據(jù)時(shí)，在東北地區(qū)該時(shí)期為常見(jiàn)農(nóng)作物生長(zhǎng)初期，苗高較低，可視為裸地。2017年8月份是常見(jiàn)作物生長(zhǎng)中后期，因此，將兩個(gè)時(shí)期的數(shù)字表面模型進(jìn)行差值計(jì)算，可獲得差值數(shù)字表面模型DDSM，得到不同作物生長(zhǎng)高度信息。

同時(shí)，考慮到不同作物株型不同，成片種植的作物在空間上呈現(xiàn)的局部鄰域高度信息不盡相同，因此再將DDSM影像進(jìn)行濾波處理，獲得鄰域波動(dòng)信息，并進(jìn)行特征評(píng)價(jià)挑選出適合農(nóng)作物分類的特征。

DDSM=DSM1-DSM2

(1)

式(1)中，DSM1為2018年6月15號(hào)DSM數(shù)據(jù)；DSM2為2017年8月3號(hào)DSM數(shù)據(jù)。根據(jù)雙時(shí)相DSM計(jì)算的DDSM，其處理結(jié)果如圖3所示。圖片中顏色越深代表高度差異越大。DDSM隨機(jī)誤差出現(xiàn)在樹(shù)木、房屋邊緣，以及水體表面漏洞。

圖4 差異數(shù)字表面模型(DDSM)

統(tǒng)計(jì)RGB影像及DDSM影像共32項(xiàng)濾波影像中各類樣本的均值，方差(D)，變異系數(shù)(V)和類間差異系數(shù)Dw。計(jì)算公式為：

D=S2

(2)

(3)

(4)

式(2)～(4)中：D為方差，S為標(biāo)準(zhǔn)差，Dw為類間變異系數(shù)，M為均值，M1為第一類均值，M2為第二類均值。

3 農(nóng)作物分類結(jié)果

3.1 農(nóng)作物分類特征

分析各類樣本的均值，方差(D)，變異系數(shù)(V)和類間差異系數(shù)Dw。亞麻與其他作物可見(jiàn)光紋理特征比較結(jié)果如表2，圖5所示。

亞麻與其他作物DDSM濾波特征比較結(jié)果如表3，圖6所示。

如表2、表3所示，各項(xiàng)特征的變異系數(shù)和相對(duì)差異有很大差異，在紋理濾波中，亞麻最大變異系數(shù)為0.76，最小變異系數(shù)為0.10，亞麻與大豆類間差異系數(shù)最大為930.12%，

表2 亞麻與其他作物可見(jiàn)光紋理特征比較 %

與建筑類間差異系數(shù)最小為0.23%； 在DDSM濾波中，亞麻的特征變異系數(shù)最大為2.34，最小為0.01，與樹(shù)木類間差異系數(shù)最大為209.20%，最小為-0.89%。部分濾波特征變異系數(shù)較大或類間差異系數(shù)較小，并非所有濾波特征都能增強(qiáng)農(nóng)作物之間的差異。

圖5 亞麻與其他作物可見(jiàn)光紋理特征比較

表3 亞麻與其他作物DDSM濾波特征比較 %

圖6 亞麻與其他作物DDSM濾波特征比較

區(qū)分不同作物的特征存在交叉，適用于分類的特征需要具有較小的變異系數(shù)和較大的類間差異系數(shù)。變異系數(shù)反應(yīng)的是離散程度，變異系數(shù)越小，越能根據(jù)該特征區(qū)分該類地物。類間差異系數(shù)反應(yīng)兩種地物在該特征上的差異，差異越大越能區(qū)分兩種地物。如圖5所示，以亞麻為例，RGB影像的二階矩紋理濾波變異系數(shù)較小，與其他作物類間差異系數(shù)較大。其中相比紅波段、藍(lán)波段二階矩濾波特征，亞麻的綠波段二階矩與其他作物類間差異系數(shù)較大并且變異系數(shù)較小，能夠較好區(qū)分11類不同類型地物，可作為有效特征。

如圖6所示，作物在同一時(shí)期高度不同，DDSM可以作為高度特征參加分類。因作物的株型不同，作物表層波動(dòng)情況表現(xiàn)為DDSM方差濾波特征和DDSM對(duì)比度濾波特征。

因此，該文分為3類組合特征：

①直接使用紅波段(R)、綠波段(G)、藍(lán)波段(B)特征進(jìn)行分類；

②將作物的可見(jiàn)光光譜信息和紋理特征進(jìn)行組合，即使用紅波段(R)、綠波段(G)、藍(lán)波段(B)和綠波段(G)二階矩濾波為組合特征進(jìn)行SVM分類；

③將作物的可見(jiàn)光光譜、紋理及高度特征進(jìn)行組合，即使用紅波段(R)、綠波段(G)、藍(lán)波段(B)、綠波段(G)二階矩濾波和DDSM方差、DDSM對(duì)比度濾波特征進(jìn)行分類。

3.2 實(shí)驗(yàn)結(jié)果

基于以上特征組合，對(duì)2017年8月無(wú)人機(jī)影像進(jìn)行分類，并采用總分類精度和Kappa系數(shù)兩項(xiàng)指標(biāo)評(píng)價(jià)分類結(jié)果。上述3種組合的分類結(jié)果及統(tǒng)計(jì)數(shù)據(jù)見(jiàn)圖7、表4。

圖7 基于特征組合的分類

表4 不同特征組合的分類精度及Kappa

特征精度KappaRGB76.000.73RGB、G二階矩82.890.81RGB、DDSM86.550.85RGB、DDSM、G二階矩、DDSM方差、DDSM對(duì)比度91.900.90

首先，基于可見(jiàn)光RGB波段的支持向量機(jī)分類的分類精度為76.00%，kappa系數(shù)為0.73，其中馬鈴薯未能分出，裸地、水稻錯(cuò)分較少，而馬鈴薯、玉米、樹(shù)木、草地錯(cuò)分較多，其中馬鈴薯錯(cuò)分為玉米和樹(shù)木。其次，在添加綠波段二階矩濾波特征后精度分類精度為82.89%，Kappa為0.81，總體精度提升，玉米、樹(shù)木、水體分類中椒鹽現(xiàn)象減少，能區(qū)分出馬鈴薯，但仍存在馬鈴薯和大豆存在錯(cuò)分情況； 再次，由RGB添加DDSM形成RGB-DDSM的 4維特征分類精度為86.55%，Kappa為0.85，其中玉米和樹(shù)木的區(qū)分效果明顯提升，地塊邊界更明顯，各類地塊內(nèi)椒鹽現(xiàn)象減少，并減少了馬鈴薯的錯(cuò)分； 最后，在RGB-DDSM波段基礎(chǔ)上添加綠波段的二階矩濾波特征和DDSM方差、對(duì)比度濾波特征后，分類精度為91.90%，Kappa系數(shù)為0.90，玉米、亞麻、大豆、水稻地塊邊界區(qū)分明顯，玉米、大豆內(nèi)噪點(diǎn)現(xiàn)象減少，草地及樹(shù)木因陰影導(dǎo)致的錯(cuò)分減少。盡管在添加綠色波段二階矩濾波、DDSM及DDSM方差和對(duì)比度后總體分類精度最高，但是分類時(shí)把部分樹(shù)的邊緣錯(cuò)分為了建筑。

4 結(jié)論與討論

在該研究中，采用無(wú)人機(jī)搭載數(shù)碼相機(jī)以獲取研究區(qū)域的RGB圖像，與側(cè)重于從高分辨率RGB圖像中提取紋理特征的分類方法所不同，該研究中利用兩個(gè)不同時(shí)期作物的DSM數(shù)據(jù)，生成差異數(shù)字表面模型DDSM，突出作物生長(zhǎng)差異特征，并創(chuàng)新性地將該特征引入農(nóng)作物分類中。進(jìn)一步地，該研究對(duì)RGB及DDSM數(shù)據(jù)進(jìn)行濾波，詳細(xì)分析并提取了農(nóng)作物在光譜及空間維度(特別是作物高度)上的聯(lián)合特征，采用SVM分類模型用于農(nóng)作物分類。實(shí)驗(yàn)結(jié)果表明，相比于僅采用光譜特征，在采用光譜、紋理、高度特征等多維特征組合時(shí)，農(nóng)作物的分類總體精度從76.00%提高到91.90%，Kappa系數(shù)由0.73提高到0.90。特別是在加入高度差這一空間特征時(shí)，整體分類精度得到顯著提高，證實(shí)了作物高程特征對(duì)提高農(nóng)作物分類精度的有效性。

通過(guò)實(shí)驗(yàn)結(jié)果發(fā)現(xiàn)，作物高程特征在農(nóng)作物分類中的主要貢獻(xiàn)如下。

(1)除了提高整體分類準(zhǔn)確性外，各類分類結(jié)果中椒鹽現(xiàn)象減少；

(2)當(dāng)引入高度特征時(shí)，在RGB影像中難以分離的類，例如樹(shù)木，草地和玉米分類精度提升；

(3)借助高度特征，RGB圖像中的由陰影引起的錯(cuò)分減少。

盡管采用該文所提出的方法取得了較高的分類精度，但仍存在以下幾個(gè)問(wèn)題。例如，在該文的實(shí)驗(yàn)結(jié)果中，在樹(shù)的邊緣誤分類現(xiàn)象較明顯。通過(guò)分析發(fā)現(xiàn)，在圖像處理過(guò)程中生成地物數(shù)字表面模型時(shí)，通常在物體的邊緣即高差發(fā)生顯著變化的區(qū)域，所生成的高度值的誤差顯著增加，因而導(dǎo)致了誤分類現(xiàn)象。因此，如何提高地物數(shù)字表面模型數(shù)據(jù)的質(zhì)量是今后需要進(jìn)行進(jìn)一步的研究的課題。同時(shí)，在該研究中僅使用了可見(jiàn)光相機(jī)，所獲取的作物特征有限，在今后的研究中將進(jìn)一步嘗試與其他類型的傳感器相結(jié)合，如低成本多光譜傳感器，增加光譜維度上的特征，以進(jìn)一步提高農(nóng)作物的分類精度。