寧 川,趙慶展,韓 峰
(1.石河子大學 信息科學與技術學院,新疆 石河子 832000;2.新疆生產建設兵團空間信息工程技術研究中心,新疆 石河子 832000)
近年來,我國人口持續增長,農業耕地面積卻日益減少,加劇了農產品需求與農業生產力之間的矛盾。因此,在現有的耕地上如何獲得最大的產出,成為一個重要的研究課題。為了緩解這個矛盾,人們在原有的生產方式基礎上引入各種新型技術,探索出更加高效的農業生產模式。精準農業是新型農業生產模式的一種,它從農業可持續發展的基本點出發,努力提高現有土地資源的利用率,使農業生產資料得到更為合理的配置,從而優化農業生產環境[1]。
精準農業實施的前提是實時準確地獲得作物生長信息和環境信息,經過分析診斷后精確調控各種環境因素,促進農業物資合理利用并實現高產[2]。在現代農業體系中,作物生長的環境信息可以通過物聯網和無線傳感網絡獲取,采集的難度不大。崔世鋼等在Android系統中設計了一個植物生長柜的監控程序,可以對多種環境因子進行監測和控制[3]。施連敏等利用無線傳感網絡構建了蔬菜生長環境的監測系統,還具有環境信息顯示和存儲功能,可以作為智能決策和控制的技術支撐[4]。
與環境條件相比,作物的生長狀態信息種類繁多,且相互之間關聯復雜,對采集方式的要求更高。在各種生長信息中,作物長勢相關參數是研究的重點。作物的長勢包括外部性狀和內部性狀,外部性狀是植株的大小、形狀和顏色等,可以直接觀察的外觀特征,是作物健康狀況的直觀反映和產量預測的基礎;內部性狀包括各種生理代謝狀況和生理生化指標,是外部性狀形成的原因。劉業亮和賈彪分別基于ZigBee無線傳感器網絡和“互聯網+”技術建立蔬菜和玉米的遠程監控平臺,可以采集作物的生長態勢,為診斷決策提供依據[5-6]。目前已有無線傳感器和多維圖像等多種先進技術應用到了作物長勢監控系統中,覆蓋了玉米和小麥等作物[7-9]。
大部分的作物長勢監控系統是依托設施農業建設的,安裝在特定的位置,檢測固定區域內的作物生長。這些系統的穩定性好,檢測的性狀種類較多,但受到覆蓋范圍的限制,無法進行大規模的作物長勢調查。無人機的全稱是無人駕駛航空飛行器,在20世紀初誕生,發展至今已經形成了具有多種型號和功能的家族。無人機的飛行速度快,機動靈活性強,非常適合作為作物信息采集的平臺。基于無人機的作物信息采集方法能夠適應農田的復雜環境,在作業效率和成本上也有著明顯的優勢,是農田作物信息采集的重要手段。王妮利用無人機和多光譜傳感器提取不同種植條件下的小麥冠層反射率,構建了小麥的長勢監測和產量預測模型[10]。周祥利用無人機搭載的多種相機拍攝水稻不同時期的影像,通過光譜和紋理等特征分析植被指數與葉面積指數和產量之間的關系,構建了水稻產量的預測模型[11]。裴浩杰等利用無人機拍攝冬小麥的高光譜影像,分析了各種內部和外部性狀以評價作物長勢,篩選到多個能夠準確反映綜合生長指標的光譜參數[12]。隨著無人機技術的不斷提升,其在農業領域的應用也將日益廣泛。
無人機所搭載的傳感器有多種類型,包括數碼相機、紅外成像儀、光譜相機和激光雷達等。這些設備針對的作物性狀不同,采集的精度也存在差異,有各自的適用范圍[13]。數碼相機拍照是無人機作物長勢調查的最普遍方式,圖像的后續處理依賴于機器視覺分析。機器視覺以計算機代替人類的視覺,分析拍攝的圖像并從中提取目標信息,得到有價值的結果。在農業領域,對作物生長信息的采集是機器視覺的應用內容之一,在多種作物的長勢監控上取得了良好的效果[14-16]。
油菜是我國重要的經濟作物之一,分布范圍廣,種植面積大。在部分地區,大規模種植的油菜不僅是油料的主要來源,春季盛開的油菜花還可以作為景觀帶動旅游產業的發展。因此,加強對油菜長勢的調查研究,有利于推動種植區域的經濟發展。本文提出了基于機器視覺的油菜長勢調查方法,利用無人機搭載數碼相機進行航拍,通過機器視覺提取油菜圖像中的特征參數,然后分析特征參數與生物量之間的相關性,評價油菜的長勢,為準確預測產量提供依據。
油菜的長勢調查設備包括無人機和機器視覺設備。無人機平臺采用大疆精靈的Phantom4型智能無人機,飛行時間最長可達30min,定位采用GPS/GLONASS兩種模式。無人機的無線通訊為WIFI形式,傳輸的速度可達Class10,最遠距離為7km。無線通訊裝置以PID控制方法和通訊接口將數據傳送回地面站,具有很好的實時性。飛行控制裝置為方向舵、副翼舵和升降舵,通過舵面的偏轉調整無人機飛行方向和姿態,從而對航線進行控制,并賦予無人機紅外避障的能力。地面站安裝大疆DJIGS Pro專業版軟件,可以規劃無人機的飛行路線,實時控制飛行速度、高度和姿態,以及啟動相機拍攝。
機器視覺的關鍵設備是尼康COOLPIX P60型數碼相機,安裝在無人機的下方,豎直向下拍攝。相機產生JPEG格式的圖像,寬度和高度上的像素分別為3 264和2 448點,經過天創UB570型圖像采集卡轉換為數字信號后用于視覺分析處理。核心計算機安裝在地面站,選用聯想揚天M4700型臺式電腦。計算機的主要配置為Intel i6中央處理器和10GB的DDR4內存,具有較好的運行速度和穩定性。計算機安裝Windows10操作系統和MatLab工具箱,用于進行圖像視覺分析和處理結果的存儲。
拍攝圖像時無人機處于飛行的狀態,因此相機與油菜田之間存在相對運動。另外,拍攝的瞬間相機在垂直高度上可能也有變化,或者拍攝角度出現微小的偏轉,這些因素都會導致油菜圖像中含有噪音,影響成像質量和后續分析的效果。針對這些問題,本文通過對圖像進行5×5的中值濾波平滑處理來消除噪音,得到油菜田的原始圖像,如圖1(a)所示。在拍攝的范圍內,相機正下方的位置與相機之間的距離最短,其它越遠的位置與相機之間的距離越長,不同位置上的物體在圖像中存在畸變。圖像的畸變會降低后續特征提取的準確性,本文參考紀壽文等的方法矯正了圖像的畸變。原始圖像由壟溝和油菜組成,壟溝包括綠色的雜草和灰色的土壤。前人的研究表明:RGB顏色坐標系統中的3個分量對目標和土壤背景的分離效果較差,本文采用HSI顏色坐標系統對圖像進行處理,同時克服了光照變化和植物陰影的影響。在HSI顏色空間下以H分量進行灰度化,得到的結果如圖1(b)所示。
圖像中的油菜處于蕾薹期和開花期,蕾薹期的油菜為綠色,始花期為綠色夾雜黃色,盛花期為黃色。這里以生物量作為反映油菜長勢的指標,從蕾薹期至盛花期生物量逐漸增大。因此,從圖像中提取與生育時期相關的特征參數,便可以對油菜的長勢進行評價。圖像中的顏色種類較少,RGB色彩空間難以全面反映特征參數。HSI色彩空間中的3個分量分別為色調H、飽和度S和亮度I,不同生育時期油菜之間的色調H變化最大。因此,將H分量做最大類間方差分析后進行二值化處理,從直方圖中獲得最大閾值T1、T2和T3。閾值分割得到壟溝、蕾薹期油菜、始花期油菜和盛花期油菜的圖像區域,分別用黑色、深灰色、淺灰色和白色表示,如圖1(c)所示。

(a) 原始圖像 (b) 灰度化圖像 (c) 特征提取圖像圖1 圖像的機器視覺分析Fig.1 Analysis of computer-vision to picture
2018年,對油菜長勢調查方法的效果進行試驗,試驗地點位于新疆維吾爾自治區石河子市南山新區,地處天山北麓中段的洪水沖積平原,屬于溫帶大陸性氣候,日照充沛,降雨量較少。試驗區域的土壤為灰潮土,種植的油菜品種為新油2號,按照常規的方式栽培。在油菜從蕾薹期向開花期過渡時采集田間影像,采集時間為上午9:00-11:00,天氣多云無風,具有較好的能見度。無人機的飛行高度100m,速度5m/s。
采集田間影像進行機器視覺分析,同時在試驗區內均勻選擇30個1m×1m的取樣點。收獲取樣點內的油菜地上部分,烘干水分至恒重,然后稱量實際的生物量;將機器視覺圖像中各取樣點的H值除以2π,使H值處于0~1之間,用于建立生物量的計算模型;30個取樣點隨機選擇15個作為建模樣本,剩余15個作為檢驗樣本。
H值與建模樣本生物量實測值的回歸分析如圖2所示。其回歸方程為y=3344.7x+934.76,決定系數為0.864,達到了極顯著的相關性,表明機器視覺圖像中的H值可以用于估算油菜的生物量。

圖2 生物量與H值的相關性Fig.2 The correlation between H value and biomass
用上述的模型估算15個檢驗樣本的生物量,并與實測值進行擬合分析,結果如圖3所示。

圖3 檢驗樣本生物量估算值與實測值的擬合Fig.3 Regression of estimated and measured biomass values (samples for validation)
回歸方程的決定系數和斜率都接近1,分別為0.940和0.960,表明建立的模型可以準確計算油菜的生物量,能夠用于油菜長勢的評價。
提出了基于機器視覺的油菜長勢調查方法,利用無人機搭載數碼相機拍攝油菜田的圖像。機器視覺分析軟件對圖像依次進行預處理和HIS色彩空間中的灰度化,然后通過閾值分割將圖像劃分為壟溝和不同生育時期的油菜。以H分量值作為與生物量相關的特征參數,用于反映油菜的長勢。在實際的應用中,H值與建模樣本的實測生物量之間存在極顯著的相關性,以此建立的模型可以準確計算油菜的生物量,能夠用于油菜長勢的評價,為產量預測提供依據。