本發(fā)明涉及精準農(nóng)業(yè)領域和無人機農(nóng)業(yè)應用領域,尤其是涉及一種基于航拍圖像的煙草植株識別與計數(shù)方法。
背景技術:
:煙草是我國主要的經(jīng)濟作物和重要的稅收來源,我國的煙草總產(chǎn)量占全世界的41.5%。為了對煙草種植進行宏觀調控和精準管理,需要對煙草種植產(chǎn)量進行精準的估算。傳統(tǒng)的方法主要有兩種:一是通過測量煙草種植的面積來估算煙草產(chǎn)量;二是通過專業(yè)的煙草技術人員人工清點煙草的株數(shù),但是這些傳統(tǒng)方法的效率低、準確度低、數(shù)據(jù)可靠性差。然而,精準農(nóng)業(yè)是當前全球農(nóng)業(yè)發(fā)展的大勢所趨,其中對農(nóng)作物的產(chǎn)量的精準估計是實現(xiàn)精準農(nóng)業(yè)管理的一個重要方面。很顯然,煙草植株產(chǎn)量估計的傳統(tǒng)方法已經(jīng)無法達到精準農(nóng)業(yè)的要求。此外,無人機技術的迅速崛起促進了我國農(nóng)業(yè)現(xiàn)代化的進程,其中最具有代表性的農(nóng)業(yè)無人機是深圳市大疆創(chuàng)新科技有限公司2015年推出的一款智能農(nóng)業(yè)噴灑防治無人機MG-1,它的作業(yè)效率高,能夠節(jié)省大量的人力物力,有利于環(huán)境保護。無人機在精準農(nóng)業(yè)領域有著廣泛的應用,無人機可以搭載各種各樣的傳感器采集環(huán)境中的數(shù)據(jù),其中搭載攝像頭進行航拍是常用的一個功能。本發(fā)明對航拍的煙草圖像進行分析,識別煙草植株并且計數(shù),可以提高統(tǒng)計效率和精確性。航拍圖片中一般包含的場景非常復雜,圖片中會有各種各樣的植物,需要從復雜的場景中識別出煙草植株。還有煙草植株與植株之間相互交叉在一起,需要把交叉在一起的煙草植株分成單株煙草。如何提取煙草植株特征,有效區(qū)分植株之間,進行分區(qū)域統(tǒng)計,是實現(xiàn)航拍圖像的煙草植株識別與計數(shù)難點和重點。技術實現(xiàn)要素:本發(fā)明提供了一種基于航拍圖像的煙草植株識別與計數(shù)方法,使用圖像處理算法來識別與計數(shù)航拍圖像中的煙草植株,從而可以進一步實現(xiàn)對煙草產(chǎn)量的精準估計。為了達到上述發(fā)明的目的,采用以下技術方案:為解決上述現(xiàn)有技術問題,本發(fā)明提供的一種基于航拍圖像的煙草植株識別與計數(shù)方法,包括以下步驟:S1:使用無人機拍攝煙草種植區(qū)域,得到煙草植株航拍圖像;S2:對航拍圖像進行預處理,分割出煙草植株的候選區(qū)域;S3:提取煙草植株候選區(qū)域的顏色特征和紋理特征用于分類器分類;S4:根據(jù)提取的煙草植株候選區(qū)域特征,使用分類器對煙草植株候選區(qū)域進行分類;S5:統(tǒng)計分類結果,在原航拍圖像中標記出檢測到的煙草植株。進一步的,所述S1中,無人機拍攝時,拍攝角度在向下,拍攝角度在與垂直方向相差30度的范圍之內。高度一般為15到25米,無人機的高度和拍攝角度必須確保拍攝到的煙草植株圖像清晰。進一步的,所述S2是利用煙草植株的中心區(qū)域與周圍區(qū)域的色差分割出煙草植株的中心區(qū)域,再對分割出的煙草植株的中心區(qū)域的圖像進行距離變換,再使用過綠操作去除非綠部分,最終獲得煙草植株候選區(qū)域。進一步的,所述的利用煙草植株中心區(qū)域與周圍區(qū)域的色差分割出煙草植株的中心區(qū)域時,將航拍圖像轉換到Lab空間,選取Lab空間中的a通道或者b通道,煙草植株的中心區(qū)域為所選取通道的圖像中的局部極值,分割出局部極值區(qū)域,分割出來的局部極值區(qū)域就是煙草植株的中心區(qū)域。在對該航拍圖像進行預處理時,可以選取Lab空間的a通道或者b通道,還可以選擇其它能體現(xiàn)色差的顏色空間中的通道,如YCBCR顏色空間中的通道。進一步的,在分割局部極值區(qū)域前需要進行降噪處理,采用的降噪方法為對已選通道的圖像進行腐蝕,然后進行形態(tài)學重建。降噪方式還可以選用其他方法,如高斯濾波器,中值濾波器。進一步的,所述對分割出煙草植株的中心區(qū)域的圖像進行距離變換,是將非煙草植株中心區(qū)域像的像素分配給離它最近的煙草中心區(qū)域。進一步的,所述S3提取的煙草植株候選區(qū)域特征包括顏色特征和紋理特征,顏色特征包括綠色通道的直方圖以及RGB空間、HSV空間每個通道的平均值和方差;紋理特征包括綠色通道的一階導數(shù)的和二階導數(shù)的平均值和方差。進一步的,所述S4中,選用的SVM分類器對候選區(qū)域進行分類,也可以選用其他類型的分類器。與傳統(tǒng)方法相比,本發(fā)明通過使用無人機拍攝煙草種植區(qū)域,利用圖像識別方法來識別與計數(shù)煙草植株,效率更高、準確度更高、數(shù)據(jù)更可靠、操作簡單;本發(fā)明將無人機用于現(xiàn)代精準農(nóng)業(yè),可以實現(xiàn)對煙草產(chǎn)量的精準估計,能夠用于煙草種植的精準管理。附圖說明圖1為本發(fā)明基于航拍圖像的煙草植株識別與計數(shù)方法的步驟流程圖;圖2為本發(fā)明利用無人機拍攝獲取的煙草植株圖像;圖3為本發(fā)明對所獲取的煙草植株圖像進行預處理的流程圖;圖4為本發(fā)明基于圖2的煙草植株圖像的局部放大圖;圖5為本發(fā)明所獲取的煙草植株圖像轉換到Lab空間中選取b通道的圖像;圖6為本發(fā)明對所獲取的煙草植株圖像進行預處理得到的煙草植株中心區(qū)域圖像;圖7為本發(fā)明基于圖4進行距離變換的示意圖;圖8為本發(fā)明所獲得的煙草植株候選區(qū)域的示意圖。圖9為本發(fā)明基于圖2的煙草植株圖像的識別與計數(shù)結果圖的局部放大圖。具體實施方式本發(fā)明有多種不同形式的實施例,下面結合附圖和實施例對本發(fā)明作進一步的描述。實施例1如圖1所示,本發(fā)明提供一種基于航拍圖像的煙草植株識別與計數(shù)方法,該方法包括以下步驟:S1,使用無人機拍攝煙草種植區(qū)域,得到煙草植株航拍圖像;S2:對航拍圖像進行預處理,分割出煙草植株的候選區(qū)域;S3:提取煙草植株候選區(qū)域的顏色特征和紋理特征用于分類器分類;S4:根據(jù)提取的煙草植株候選區(qū)域特征,使用分類器對煙草植株候選區(qū)域進行分類;S5:統(tǒng)計分類結果,在原航拍圖像中標記出檢測到的煙草植株。如圖2所示,首先使用大疆無人機精靈4在煙草種植田的上空進行航拍獲取煙草植株圖像,總共拍攝了12幅圖像。為確保所獲得的該煙草植株圖像中的煙草植株的清晰度,拍攝前對無人機的高度和角度進行調整,高度為15到25米,拍攝角度垂直向下。如圖3所示,按照圖3的步驟對航拍圖像進行預處理,利用煙草植株的中心區(qū)域與周圍區(qū)域的色差分割出煙草植株的中心區(qū)域。如圖4所示,由于煙草植株的中心區(qū)域32都是嫩葉,所以煙草植株的中心區(qū)域與周圍區(qū)域31存在色差,在該實施例中,將該航拍圖像從RGB空間轉換到Lab空間,選取Lab空間的b通道。如圖5所示,可以觀察到該煙草植株的中心區(qū)域為局部最亮區(qū)域。對所選取的通道的圖像使用大小為5的結構元進行腐蝕,然后對進行過腐蝕的圖像進行形態(tài)學重建,去除噪聲。對進行形態(tài)學重建后的圖像分割局部最大值,如圖6所示,分割出的每一個連通區(qū)域表示一棵煙草植株的中心區(qū)域。如圖7所示,對分割出的該煙草植株中心區(qū)域的圖像進行距離變換,將圖像中的非煙草中心區(qū)域的像素分配給離它最近的煙草中心區(qū)域的植株形成一個完整的煙草植株區(qū)域。如圖8所示,使用過綠操作去除該完整煙草植株區(qū)域中的非綠色部分,最終得到該煙草植株的候選區(qū)域。提取上述的煙草植株候選區(qū)域的特征。提取的特征主要包括顏色特征和紋理特征,該顏色特征包括綠色通道的直方圖以及RGB空間、HSV空間每個通道的平均值和方差;該紋理特征包括綠色通道的一階導數(shù)的和二階導數(shù)的平均值和方差。根據(jù)提取的煙草植株候選區(qū)域特征,使用分類器對該煙草植株候選區(qū)域進行分類,并統(tǒng)計分類結果,在該航拍圖像中標記出檢測到的煙草植株。在該實施例中,將無人機拍攝的12幅圖像分成訓練集和測試集,訓練集2幅,測試集10幅。首先,隨機的對訓練集中的該煙草植株候選區(qū)域進行人工標記樣本,煙草植株區(qū)域標記為正樣本,非煙草植株標記為負樣本,正負樣本的數(shù)量各1500個,使用這些樣本的特征訓練出SVM分類模型。然后,使用SVM分類模型對該煙草植株候選區(qū)域進行分類,提取每一顆煙草植株候選區(qū)域的特征,根據(jù)每一顆煙草植株候選區(qū)域的特征使用訓練好的該SVM分類模型進行分類。如圖9所示,黑色小方塊表示檢測到的煙草植株。請參閱表1,統(tǒng)計分類的結果為正的區(qū)域數(shù),則為煙草植株的數(shù)量。表1圖像編號真實數(shù)量(株)檢測到的數(shù)量(株)訓練集0117541809訓練集0217671813測試集0124572413測試集0230992718測試集0328662648測試集0422972348測試集0516771796測試集0625642429測試集0720541995測試集0812081263測試集0932823460測試集1034303195當前第1頁1 2 3