吳雪梅 梁長(zhǎng)江 張大斌 喻麗華 張富貴

(貴州大學(xué)機(jī)械工程學(xué)院， 貴陽(yáng) 550025)

0 引言

覆膜技術(shù)能夠有效優(yōu)化和改善栽培條件，其蓄水、保溫、保墑、保藥、除草等作用顯著，利于作物早熟，提高農(nóng)作物產(chǎn)量[1-2]。然而，由于地膜回收技術(shù)不成熟、農(nóng)戶缺乏環(huán)保意識(shí)等，造成大量地膜殘留在農(nóng)田中，形成殘膜，破壞了土壤的連續(xù)性，降低土壤的含水率、透氣性和肥力水平，導(dǎo)致農(nóng)作物減產(chǎn)[3-5]。掌握農(nóng)田殘膜時(shí)空分布圖、精確獲取殘膜時(shí)空分布信息及時(shí)空變化特征，可為相關(guān)部門開展農(nóng)田環(huán)境健康評(píng)估、地膜回收管理、農(nóng)作物監(jiān)測(cè)、殘膜污染防控等提供管理決策依據(jù)。傳統(tǒng)方法采用實(shí)地調(diào)查獲取地膜覆蓋和殘膜分布信息，該方法主觀性強(qiáng)，耗時(shí)、費(fèi)力。衛(wèi)星遙感技術(shù)可克服現(xiàn)場(chǎng)調(diào)研的不足，具有覆蓋面積廣、檢測(cè)成本低等優(yōu)勢(shì)。因此，基于遙感的農(nóng)田地膜覆蓋監(jiān)測(cè)和地膜識(shí)別技術(shù)成為研究熱點(diǎn)。

李佳雨等[6]基于Landsat衛(wèi)星和資源三號(hào)衛(wèi)星，獲取了甘肅中部地區(qū)6月的覆膜農(nóng)田圖像，分析了圖像的光譜特征、紋理特征、指數(shù)特征和幾何特征，利用隨機(jī)森林分類法訓(xùn)練地膜樣本并進(jìn)行分類，實(shí)現(xiàn)了地膜的快速識(shí)別，精度達(dá)90%。KOC-SAN[7]利用Word View-2衛(wèi)星采集安塔利亞研究區(qū)影像，利用隨機(jī)森林法、最大似然法和支持向量機(jī)等分類方法識(shí)別塑料薄膜，識(shí)別精度均高于90%。HASI等[8]基于Landsat-8遙感數(shù)據(jù)，提取了4月的地膜影像數(shù)據(jù)，結(jié)合光譜與紋理的組合特征，對(duì)農(nóng)田中的地膜分別采用支持向量機(jī)(SVM)、最大似然分類器(MLC)和最小距離分類器(MDC)方法進(jìn)行了分類識(shí)別，發(fā)現(xiàn)SVM分類法精度最高，達(dá)94%。杭丹維[9]利用MODIS遙感影像，提取了4月和5月的覆膜影像數(shù)據(jù)，結(jié)合改進(jìn)的亞像元空間引力模型(SSSAM)，完成了覆膜農(nóng)田分布信息的高精度識(shí)別。高夢(mèng)婕等[10]基于GF-2遙感衛(wèi)星，獲取了赤峰市研究區(qū)塑料大棚覆蓋薄膜的階段遙感影像，提取了圖像的最優(yōu)分割參數(shù)，對(duì)比了決策樹、隨機(jī)森林和支持向量機(jī)3種分類方法對(duì)大棚目標(biāo)的分類效果，發(fā)現(xiàn)隨機(jī)森林法能提取連片塑料大棚，分類精度達(dá)89.65%。沙先麗[11]利用Landsat TM獲取了不同年份4月和5月上旬新疆棉區(qū)的遙感影像，分析影像的遙感特征、紋理特征和指數(shù)特征，并建立決策樹模型，實(shí)現(xiàn)了覆膜棉田的分類識(shí)別，總體精度達(dá)90%。TARANTINA等[12]采用面向?qū)ο蠓椒ɑ诤娇沼跋駲z測(cè)了塑料大棚葡萄園。

近年來(lái)，無(wú)人機(jī)作為遙感平臺(tái)，在作物生長(zhǎng)監(jiān)測(cè)方面顯示出巨大的潛力[13-16]，對(duì)比衛(wèi)星遙感影像，無(wú)人機(jī)遙感影像分辨率高、質(zhì)量好、成本低。但基于無(wú)人機(jī)的農(nóng)田地膜識(shí)別研究較少。朱秀芳等[17]以無(wú)人機(jī)為遙感平臺(tái)，獲取了5月初期云南省覆膜農(nóng)田圖像，利用隨機(jī)森林優(yōu)選圖像紋理特征，并進(jìn)行最大似然法分類，結(jié)合形態(tài)學(xué)算法與面積閾值分割處理提取覆膜農(nóng)田面積，識(shí)別精度達(dá)94.84%。孫鈺等[18]利用無(wú)人機(jī)采集了赤峰市大棚和覆膜農(nóng)田圖像，結(jié)合神經(jīng)網(wǎng)絡(luò)和多尺度融合等深度學(xué)習(xí)方法，快速識(shí)別大棚和覆膜農(nóng)田，識(shí)別率高達(dá)97%。

以上基于衛(wèi)星和無(wú)人機(jī)遙感的農(nóng)田地膜識(shí)別，識(shí)別目標(biāo)有地膜和塑料大棚；特征提取包括地膜光譜特征、指數(shù)特征、幾何特征和紋理特征；地膜影像獲取時(shí)間大多為4月或5月，即農(nóng)田地膜覆膜初期或作物生長(zhǎng)初期的遙感影像。此時(shí)，地膜表面受裸土的干擾小，地膜具有覆蓋面積大、連續(xù)成條狀、特征明顯等特點(diǎn)，利用隨機(jī)森林法、最大似然法、閾值分割和支持向量機(jī)等分類方法可獲得較高的識(shí)別精度。對(duì)農(nóng)田地膜覆蓋后期，即作物收獲期后的農(nóng)田地膜，因地膜未能及時(shí)撿拾而在農(nóng)田間殘留，形成殘膜，殘膜受裸土和風(fēng)雨侵蝕的影響，不再連接成條狀，殘膜與裸土的紋理和形狀等特征差異不明顯，利用以上的特征和方法進(jìn)行地膜識(shí)別分類，精度較低。為此，本文基于無(wú)人機(jī)遙感可見光圖像，針對(duì)秋收后期的農(nóng)田殘膜和土壤數(shù)據(jù)集特點(diǎn)，提出一種基于顏色特征的殘膜識(shí)別方法，旨在克服光照強(qiáng)度對(duì)目標(biāo)識(shí)別的影響，獲得更為準(zhǔn)確的殘膜覆蓋面積及其地理分布信息。

1 材料和方法

1.1 無(wú)人機(jī)圖像獲取

試驗(yàn)數(shù)據(jù)為2018年10—12月進(jìn)行航拍的無(wú)人機(jī)影像，航拍地點(diǎn)為貴州省安順市平壩區(qū)，北緯26°22′，東經(jīng)106°31′，海拔1 217 m，田間土壤以黃壤為主，種植作物為煙葉。所使用的無(wú)人機(jī)為大疆精靈4(Phantom4)四旋翼無(wú)人機(jī)，該機(jī)飛行荷載1.38 kg，最大上升速度為21.6 km/h，飛行速度為72 km/h，最大下降速度為14.4 km/h，最大任務(wù)載荷為6 kg，續(xù)航時(shí)間為28 min，控制距離5 000 m。搭載的相機(jī)為內(nèi)置DJI Lightbridge高清視頻傳輸系統(tǒng)，傳感器型號(hào)為1/2.3英寸 CMOS，有效像素為1 240萬(wàn)，影像尺寸4 000像素×3 000像素。航拍高度為18～55 m，旁向重疊率30%，航向重疊率70%。航拍時(shí)天氣晴朗、微風(fēng)。為研究農(nóng)田覆膜整個(gè)生命周期的地膜識(shí)別方法，針對(duì)煙葉不同生長(zhǎng)期，采集煙葉6葉期、煙葉采摘后、煙桿拔除后的農(nóng)田地膜或殘膜圖像，不同飛行高度各采集20幅圖像，共1 500幅。采集到的數(shù)據(jù)源為JPG格式影像數(shù)據(jù)。無(wú)人機(jī)采集影像前，在飛行區(qū)域內(nèi)布置黑白板，用于相機(jī)標(biāo)定，降低系統(tǒng)誤差。無(wú)人機(jī)與田間殘膜圖像如圖1所示。

圖1 無(wú)人機(jī)與田間殘膜圖像Fig.1 UAV and image of plastic film residual

1.2 研究方法

1.2.1殘膜識(shí)別流程

本文旨在提取田間殘膜，即當(dāng)季收獲期后的地膜。晴朗或多云天氣，隨著太陽(yáng)入射角度的變化，田間會(huì)出現(xiàn)陰影區(qū)，陰影區(qū)和陽(yáng)光直射區(qū)的殘膜和土壤顏色特征有明顯差異，為識(shí)別出目標(biāo)，需分析出該差異的特點(diǎn)，再分割目標(biāo)。為此，本文設(shè)計(jì)殘膜識(shí)別流程，如圖2所示。主要步驟如下：①對(duì)UAV航拍的圖像進(jìn)行畸變校正。②進(jìn)行3×3中值濾波和直方圖均衡化處理。③提取HSV空間的S分量和RGB空間的B分量作為分割算子。④基于顏色分割算子，開展不同分割算法的分割效果對(duì)比，優(yōu)選最佳分割方法和分割算子。⑤提出殘膜識(shí)別方法，對(duì)分割后圖像進(jìn)行形態(tài)學(xué)處理，提取殘膜計(jì)算面積和空間分布信息。⑥計(jì)算識(shí)別率，進(jìn)行精度驗(yàn)證，輸出識(shí)別結(jié)果。

1.2.2鏡頭畸變校正

利用AgisoftPhotoScan軟件對(duì)采集的圖像進(jìn)行拼接處理和幾何校正，基于幾何參考板GPS 控制點(diǎn)進(jìn)行圖像幾何校正，去除無(wú)人機(jī)姿態(tài)變化、大氣折射等影響[19]。

1.2.3殘膜識(shí)別分割方法

(1)顏色分量

RGB顏色模型有紅(Red)、綠(Green)、藍(lán)(Blue)3個(gè)通道，取值范圍為0～255，各個(gè)分量的數(shù)值越小，亮度越低；數(shù)值越大，亮度越高。HSV(hue, saturation, value)是RGB 模型的一種非線性變換，接近人類的色彩感知。其中H表示色相(hue)，是由物體反射或透射出來(lái)的顏色，在色相環(huán)上用0°～360°來(lái)表示；S表示飽和度(saturation)，是指顏色的強(qiáng)度，用色相中彩色成分所占的比例來(lái)表示，取值范圍為0～1，值越大顏色越飽和；V表示明度(value)，反映顏色的相對(duì)明暗程度，取值范圍為0～255，取0時(shí)為黑色，取255時(shí)為白色。

(2)顏色特征提取

煙桿拔除后，農(nóng)田的主要景物為土壤和殘膜。殘膜被壟體遮擋時(shí)，無(wú)陽(yáng)光直射，稱為陰影區(qū)，有陽(yáng)光直射，稱為直射區(qū)，如圖3所示。相同目標(biāo)在直射區(qū)和陰影區(qū)的顏色灰度有差異。為探求目標(biāo)和背景的顏色特征，將圖像中的景物分為4類：直射區(qū)殘膜、直射區(qū)土壤、陰影區(qū)殘膜和陰影區(qū)土壤。

圖3 直射區(qū)與陰影區(qū)示意圖與實(shí)物圖Fig.3 Direct light area and shadow area schematic and color image

直方圖分析法統(tǒng)計(jì)圖像中景物灰度出現(xiàn)的范圍和頻率，卻忽略了不同灰度所對(duì)應(yīng)景物的信息，針對(duì)目標(biāo)和背景灰度受光照影響大的特點(diǎn)，截取圖像中直射區(qū)殘膜、直射區(qū)土壤、陰影區(qū)殘膜和陰影區(qū)土壤4個(gè)目標(biāo)，分別統(tǒng)計(jì)4個(gè)目標(biāo)的顏色分量灰度，并繪制箱線圖，目的是通過(guò)對(duì)比分析，獲取目標(biāo)的顏色特點(diǎn)，從而確定用于圖像分割的顏色分量。

(3)分割方法

圖像分割依據(jù)圖像中目標(biāo)和背景在形狀、顏色、灰度、紋理和空間等特征上的差異，分割目標(biāo)和背景，分割特征具有相同的特點(diǎn)：在目標(biāo)區(qū)域內(nèi)具有較大的一致性和相似性，在背景區(qū)域內(nèi)具有極大的差異性[20-21]。為獲取較高的識(shí)別精度，對(duì)比分析了手動(dòng)閾值法、迭代閾值法、最大類間方差法、最大熵值法、K-means均值聚類和脈沖耦合神經(jīng)網(wǎng)絡(luò)法。

圖4 殘膜與土壤的灰度分布箱線圖Fig.4 Gray-scale distribution boxes of residual film and soil

(4)精度驗(yàn)證

殘膜實(shí)際面積采用實(shí)測(cè)法獲取，將1 m2的黑色方框置于煙田，截取黑框內(nèi)表面裸露的殘膜于實(shí)驗(yàn)室測(cè)量其表面積[2]，測(cè)量面積即為殘膜實(shí)際面積。

采用本文方法識(shí)別出該黑色方框內(nèi)的殘膜，利用面積占優(yōu)法計(jì)算殘膜面積。統(tǒng)計(jì)5個(gè)黑色方框的實(shí)際面積和識(shí)別面積，利用識(shí)別率評(píng)價(jià)識(shí)別精度，識(shí)別率為殘膜識(shí)別面積與殘膜實(shí)際面積之比，5個(gè)方框的平均識(shí)別率越大，則識(shí)別精度越高。

2 結(jié)果與討論

2.1 顏色分量提取

殘膜圖像灰度統(tǒng)計(jì)箱線圖如圖4所示，直射區(qū)殘膜的R、G、B分量灰度最高，陰影區(qū)土壤灰度最低，陰影區(qū)殘膜灰度處于直射區(qū)土壤和陰影區(qū)土壤之間，且陰影區(qū)土壤與陰影區(qū)殘膜灰度重疊較少，區(qū)分度較高。R分量圖中，直射區(qū)土壤、陰影區(qū)殘膜和陰影區(qū)土壤的灰度區(qū)間重疊較少，但與直射區(qū)殘膜重疊較多；G、B分量的直射區(qū)殘膜灰度高于其他3種景物，G分量的陰影區(qū)殘膜和直射區(qū)土壤灰度區(qū)間重合較多，B分量圖中4種景物的灰度區(qū)間相互重疊較多。就HSV顏色模型而言，H、V分量圖中直射區(qū)殘膜、直射區(qū)土壤、陰影區(qū)殘膜和陰影區(qū)土壤的灰度依次減小，其中，H分量圖陰影區(qū)殘膜和陰影區(qū)土壤灰度區(qū)間重疊較多，V分量圖直射區(qū)殘膜和直射區(qū)土壤灰度區(qū)間有部分重疊，S分量圖直射區(qū)殘膜與直射區(qū)土壤以及陰影區(qū)殘膜與陰影區(qū)土壤的灰度區(qū)間重疊較少。

B、S分量中殘膜和土壤的灰度區(qū)間重合較少，區(qū)分度較高，利于分割，故本文以B、S分量作為分割算子。

2.2 殘膜圖像分割

基于B、S分量分別進(jìn)行手動(dòng)閾值分割、迭代閾值分割、最大類間方差法分割、最大熵值法分割、K-means均值聚類分割和脈沖耦合神經(jīng)網(wǎng)絡(luò)分割，通過(guò)對(duì)比研究獲取最佳殘膜分割方法。

2.2.1基于B顏色分量的殘膜圖像分割

對(duì)無(wú)人機(jī)采集的田間殘膜圖像進(jìn)行鏡頭畸變校正、3×3模板中值濾波和灰度直方圖均衡化處理，然后提取B分量。

(1)手動(dòng)閾值分割

基于B分量進(jìn)行手動(dòng)閾值分割處理，如圖5所示，可分割出直射區(qū)殘膜，但不能分割出陰影區(qū)殘膜，閾值為150時(shí)，存在過(guò)分割現(xiàn)象，把直射區(qū)土壤識(shí)別為殘膜，如圖5中的圓圈所示，閾值為200時(shí)，無(wú)過(guò)分割，直射區(qū)殘膜分割效果較好。

圖5 基于B分量的手動(dòng)閾值分割結(jié)果Fig.5 Manual threshold segmentation results based on B component

(2)迭代閾值分割

基于B分量進(jìn)行不同迭代條件的迭代閾值分割處理，結(jié)果如圖6所示。與圖5a相比，當(dāng)?shù)撝礣<0.1、T<1、T<10和T<15時(shí)，均能分割出直射區(qū)殘膜，不能識(shí)別出陰影區(qū)殘膜。

圖6 基于B分量的迭代閾值分割結(jié)果Fig.6 Iterative threshold segmentation results based on B component

(3)其他方法分割

最大類間方差法、最大熵值法、K-means均值聚類法和脈沖耦合神經(jīng)網(wǎng)絡(luò)法分割效果如圖7所示。對(duì)比圖5a發(fā)現(xiàn)，方框標(biāo)記處的陰影區(qū)殘膜均出現(xiàn)了欠分割現(xiàn)象，4種方法都不能識(shí)別出陰影區(qū)目標(biāo)，只能識(shí)別出直射區(qū)目標(biāo)，其中，脈沖耦合神經(jīng)網(wǎng)絡(luò)分割法識(shí)別效果相對(duì)較好，識(shí)別率較高。

綜合對(duì)比6種分割法，迭代閾值法和脈沖耦合神經(jīng)網(wǎng)絡(luò)法的識(shí)別效果較好。但因陰影區(qū)殘膜灰度介于直射區(qū)和陰影區(qū)土壤灰度間，重疊區(qū)域大，故基于B分量不能分割陰影區(qū)目標(biāo)。

圖7 基于B分量的4種分割方法分割效果對(duì)比Fig.7 Segmentation results of maximum inter-class variance, maximum entropy method, K-means clustering method and pulse-based neural network method based on B component

2.2.2基于S分量的殘膜圖像分割

基于S分量的手動(dòng)閾值分割方法結(jié)果如圖8所示。當(dāng)手動(dòng)閾值為125、120、110時(shí)，直射區(qū)和部分陰影區(qū)殘膜可識(shí)別，圖8方框內(nèi)目標(biāo)均被分割出，但陰影區(qū)目標(biāo)有欠分割現(xiàn)象，隨著分割閾值的降低，欠分割程度加強(qiáng)，如圖8中的橢圓標(biāo)記區(qū)的目標(biāo)均未被識(shí)別。基于S分量的手動(dòng)閾值法分割效果欠佳。

利用迭代閾值法、最大熵值法、最大類間方差法、脈沖耦合神經(jīng)網(wǎng)絡(luò)法和K-means均值聚類法分割圖像，效果如圖9所示，5種方法均能識(shí)別出圖中方框標(biāo)記的陰影區(qū)和直射區(qū)殘膜，但最大熵值法和K-means均值聚類法欠分割較嚴(yán)重，如圖9b、9e中橢圓所示，識(shí)別出的陰影區(qū)殘膜面積比實(shí)際面積小。基于脈沖耦合神經(jīng)網(wǎng)絡(luò)法的識(shí)別效果最佳，最大類間方差法和迭代閾值法相當(dāng)，K-means均值聚類法最差，5種分割方法均優(yōu)于手動(dòng)閾值法。

綜上，基于S分量的脈沖耦合神經(jīng)網(wǎng)絡(luò)法的分割效果最好，可識(shí)別出陽(yáng)光直射區(qū)和陰影區(qū)目標(biāo)，準(zhǔn)確度高。該方法可消除光照強(qiáng)度對(duì)目標(biāo)識(shí)別的影響，解決收獲后田間殘膜識(shí)別的問(wèn)題。

圖8 基于S分量的手動(dòng)閾值分割效果對(duì)比Fig.8 Manual threshold segmentation effects comparison based on S component

圖9 基于S分量的不同方法分割效果對(duì)比Fig.9 Segmentation effects comparison of different methods based on S component

2.3 殘膜識(shí)別率驗(yàn)證

對(duì)不同時(shí)期農(nóng)田殘膜進(jìn)行識(shí)別，統(tǒng)計(jì)殘膜識(shí)別面積、實(shí)際面積和識(shí)別率，結(jié)果如表1所示。6葉期的識(shí)別率最高，為96.99%，該時(shí)期煙田背景簡(jiǎn)單，地膜覆蓋時(shí)間短，裸土對(duì)地膜的影響小，目標(biāo)和背景的區(qū)分度高，便于分割。煙葉收獲后未拔桿時(shí)的殘膜識(shí)別率為69.47%，該時(shí)期煙桿及頂部葉芽遮擋殘膜，造成識(shí)別率低。煙桿拔除后殘膜識(shí)別率為93.55%，低于6葉期地膜識(shí)別率，因地膜長(zhǎng)時(shí)間覆蓋后，顏色變黃，與土壤灰度交叉，導(dǎo)致識(shí)別率低。冬季空閑期的識(shí)別率為88.95%，此時(shí)雜草枯黃，殘膜上方混雜土壤，顏色差異不顯著。地膜覆蓋周期的平均識(shí)別率為87.49%。

表1 不同時(shí)期田間殘膜識(shí)別率Tab.1 Identification accuracy of different growth stages

3 結(jié)論

(1)基于無(wú)人機(jī)遙感影像，以RGB、HSV顏色空間模型的B、S分量作為分割算子，分割田間殘膜目標(biāo)，S分量可消除光照強(qiáng)度對(duì)目標(biāo)識(shí)別的影響。

(2)基于S、B分割算子，利用迭代閾值法、最大熵值法、脈沖耦合神經(jīng)網(wǎng)絡(luò)法、手動(dòng)閾值法、最大類間方差法、K-means均值聚類法分別對(duì)田間殘膜圖像進(jìn)行分割。基于B分量的6種分割法均能分割出直射區(qū)殘膜，但不能分割陰影區(qū)殘膜；基于S分量的6種分割方法均能分割出直射區(qū)殘膜和陰影區(qū)殘膜，基于脈沖耦合神經(jīng)網(wǎng)絡(luò)方法的分割效果最佳。

(3)以S分量作為分割算子，利用基于脈沖耦合神經(jīng)網(wǎng)絡(luò)法識(shí)別不同時(shí)期的田間殘膜，識(shí)別率有較大差異，6葉期識(shí)別率最高，煙葉收獲后的識(shí)別率最低，地膜覆蓋周期的平均識(shí)別率為87.49%。