基于分區(qū)域特征點(diǎn)聚類的秧苗行中心線提取

廖 娟 汪 鷂 尹俊楠 張 順 劉 路 朱德泉

(安徽農(nóng)業(yè)大學(xué)工學(xué)院， 合肥 230036)

0 引言

作為智能農(nóng)業(yè)機(jī)械的重要組成部分，農(nóng)業(yè)機(jī)械自動(dòng)導(dǎo)航技術(shù)備受人們關(guān)注，并已廣泛應(yīng)用于農(nóng)作物種植、植保、收獲等方面[1-3]。為適應(yīng)現(xiàn)代化農(nóng)業(yè)的發(fā)展需求，水稻插秧機(jī)的自主行走成為農(nóng)業(yè)機(jī)械智能化的重要研究內(nèi)容，在水稻插秧機(jī)自主行走技術(shù)中，自動(dòng)導(dǎo)航技術(shù)是其重要組成部分[4-5]。目前，農(nóng)業(yè)機(jī)械自動(dòng)導(dǎo)航領(lǐng)域的研究主要集中在機(jī)器視覺技術(shù)和衛(wèi)星定位(GPS)技術(shù)兩種方式上[6]，其中視覺導(dǎo)航具有成本低、范圍廣、適應(yīng)性強(qiáng)等特點(diǎn)，已成為農(nóng)業(yè)機(jī)械導(dǎo)航的研究熱點(diǎn)。

根據(jù)我國水稻按照近似直線的行種植特點(diǎn)，插秧機(jī)視覺導(dǎo)航可通過提取秧苗行中心線，為確定插秧機(jī)的相對位置提供依據(jù)[7]。關(guān)于作物行中心線的提取，國內(nèi)外學(xué)者提出了多種方法，其中應(yīng)用較多的是基于Hough變換的作物行中心線檢測方法，該方法結(jié)合色彩特征分割出圖像中的作物或壟溝像素，通過Hough變換擬合這些像素，以提取作物行中心線[8-10]。但標(biāo)準(zhǔn)的Hough算法往往存在峰值檢測難、計(jì)算復(fù)雜度高等缺點(diǎn)，難以滿足農(nóng)田作業(yè)實(shí)時(shí)處理的要求，且農(nóng)田環(huán)境中光照、雜草等因素會(huì)直接影響作物綠色特征的識(shí)別[11]，這給準(zhǔn)確提取作物行中心線帶了一定困難。

近年來，研究者開始關(guān)注改進(jìn)Hough變換的應(yīng)用。文獻(xiàn)[12]提取秧苗特征點(diǎn)，用近鄰聚類法對特征點(diǎn)進(jìn)行聚類，并采用基于已知點(diǎn)的Hough變換[13]提取秧苗行中心線。文獻(xiàn)[14]結(jié)合過綠特征和中值濾波算法分割作物行像素，通過隨機(jī)Hough變換[15]擬合出作物行中心線。盡管改進(jìn)的Hough算法在一定程度上降低了內(nèi)存消耗，提高了線提取速度，但在噪聲大的農(nóng)田場景下，作物行中心線檢測的準(zhǔn)確性和實(shí)時(shí)性仍有待提高。

除Hough及其相應(yīng)的改進(jìn)算法，最小二乘法也常用于作物行檢測[16-17]。相比于Hough變換，最小二乘法實(shí)時(shí)性好，但其準(zhǔn)確性易受噪聲影響，作物行中心線的魯棒性不佳。針對最小二乘法對噪聲敏感等問題，文獻(xiàn)[18-19]基于顏色特征分割出綠色像素，利用線性回歸算法提取作物行，克服了傳統(tǒng)Hough變換的不足，但該方法需已知作物行間距。文獻(xiàn)[20]假設(shè)作物行間距相等，采用多窗口移動(dòng)方法對特征點(diǎn)進(jìn)行聚類，獲取代表作物行的特征點(diǎn)，擬合特征點(diǎn)提取作物行中心線，該方法能有效去除雜草噪聲的影響，但行間距相等的限定使其適用范圍受限。文獻(xiàn)[21]采用垂直投影法提取水稻稻田圖像秧苗行定位點(diǎn)，利用定位點(diǎn)的鄰近關(guān)系對定位點(diǎn)分類，擬合定位點(diǎn)獲取秧苗行中心線，但垂直投影法定位性較差，如不對定位點(diǎn)進(jìn)行篩選，少量誤差較大的定位點(diǎn)便會(huì)對最終的直線擬合效果產(chǎn)生較大的影響。

針對上述問題，為提高秧苗行中心線提取算法的實(shí)時(shí)性和精度，本文提出一種基于分區(qū)域特征點(diǎn)聚類的秧苗行中心線提取方法。在對圖像進(jìn)行預(yù)處理后，基于分區(qū)域的思想獲取秧苗行候選特征點(diǎn)，并結(jié)合垂直投影方法確定秧苗行數(shù)以及起始點(diǎn)，利用特征點(diǎn)的鄰近關(guān)系對特征點(diǎn)聚類，根據(jù)秧苗行的特征點(diǎn)分布，確定秧苗行的趨勢線，通過趨勢線與特征點(diǎn)的位置關(guān)系，篩選特征點(diǎn)，去除噪聲點(diǎn)，并用最小二乘法擬合，獲取秧苗行中心線。

1 秧苗行像素點(diǎn)分析

傳統(tǒng)的作物行中心線提取方法，多是基于作物像素的顏色、亮度等特征將作物從土壤背景中提取出來，采用Hough變換、最小二乘法等方法對分割后的二值圖像進(jìn)行直線擬合。但分割出的一系列作物像素點(diǎn)還包含許多噪聲和冗余的像素點(diǎn)，會(huì)對行中心線的識(shí)別造成干擾，降低算法的運(yùn)行速度。

圖1 二值圖像的垂直投影曲線Fig.1 Vertical projection curves of binary image

圖1所示為秧苗圖像經(jīng)二值化后對應(yīng)的垂直投影曲線，其中，秧苗行線所在位置出現(xiàn)投影值變化的峰值，從土壤背景到秧苗行，垂直投影值向上躍變，從秧苗行到土壤背景，垂直投影值向下躍變。但在秧苗圖像中由于透視投影作用，相鄰兩條平行的秧苗行會(huì)交于圖像內(nèi)部或外部一點(diǎn)，造成部分秧苗行像素出現(xiàn)跨行顯示，如圖1中紅線標(biāo)記的區(qū)域。第1行中的部分秧苗像素點(diǎn)出現(xiàn)在第2行，會(huì)被誤檢為第2行的秧苗像素點(diǎn)，同時(shí)，田間場景常伴有天氣、光線變化等，這些干擾因素給秧苗分割帶來困難，產(chǎn)生誤分割，從而影響秧苗行中心線的正確提取。此外，二值圖像中包含大量秧苗像素點(diǎn)，如果直接對整幅圖像進(jìn)行行中心線提取，會(huì)增加算法的計(jì)算開銷。

從圖1的二值圖像中隨機(jī)提取部分秧苗像素點(diǎn)來表示秧苗行特征信息，如圖2所示。圖中紅點(diǎn)為隨機(jī)提取的特征點(diǎn)，統(tǒng)計(jì)圖中每行白色像素點(diǎn)和紅色特征點(diǎn)的個(gè)數(shù)，結(jié)果如表1所示。采用運(yùn)算簡單的最小二乘法分別對二值圖像和隨機(jī)特征點(diǎn)進(jìn)行秧苗行中心線提取，得到擬合直線的耗時(shí)分別為0.379、0.295 s，行中心線提取結(jié)果如圖3所示。由處理結(jié)果可以看出，基于秧苗行特征點(diǎn)提取秧苗行中心線可消除冗余的秧苗像素點(diǎn)對檢測精度的影響，又可提高檢測速度，但秧苗行特征點(diǎn)選取準(zhǔn)確與否直接影響提取秧苗行中心線準(zhǔn)確性?；诖耍?本文提出一種基于分區(qū)域特征點(diǎn)聚類的秧苗行中心線提取方法。

圖2 隨機(jī)選取的特征點(diǎn)圖Fig.2 Random selection of feature points

個(gè)

圖4 圖像預(yù)處理Fig.4 Image preprocessing

圖3 秧苗行中心線提取結(jié)果Fig.3 Seedling row lines extraction

2 秧苗行中心線提取

本文提出的基于分區(qū)域特征點(diǎn)聚類的秧苗行中心線提取方法主要包括圖像預(yù)處理、秧苗行特征點(diǎn)獲取以及秧苗行中心線擬合，首先，對采集的秧苗圖像進(jìn)行預(yù)處理，獲取秧苗的二值圖像；然后，對二值圖像進(jìn)行分區(qū)域特征點(diǎn)聚類，提取每一秧苗行的特征點(diǎn)，并利用最小二乘法擬合秧苗行的特征點(diǎn)，獲取秧苗行中心線。

2.1 圖像預(yù)處理

圖像預(yù)處理包括秧苗圖像灰度化和圖像二值化處理。秧苗圖像中秧苗與非秧苗背景的綠色分量差別較大，基于此，采用超綠特征提取方法[22]進(jìn)行秧苗圖像灰度化，超綠特征提取公式為

I=2G-R-B

(1)

式中R、G、B——彩色圖像中紅色、綠色、藍(lán)色分量

I——轉(zhuǎn)換后灰度圖像的亮度值

將圖4a所示原始圖像轉(zhuǎn)換為灰度圖像，如圖4b所示，有效提高綠色秧苗和非綠色背景的對比度，有利于非綠色背景的去除。

彩色圖像經(jīng)超綠特征法轉(zhuǎn)換為灰度圖像后，仍不能直接提取圖像中秧苗像素坐標(biāo)，因此需要對圖像進(jìn)行二值化處理，分割出秧苗像素和背景像素。對圖4b所示灰度圖像進(jìn)行最大類間方差法(簡稱Otsu)分割[23]，獲得二值圖像，如圖4c所示，圖像中白色像素點(diǎn)為分割的秧苗像素點(diǎn)，黑色為背景。

2.2 候選特征點(diǎn)提取

為減小后期圖像處理的計(jì)算量，聚類前先在二值圖像中分區(qū)域提取秧苗行特征點(diǎn)，具體步驟如下：

(1)將圖4c所示二值圖像等分為上下兩部分，對于上半部分再次等分為n份，下半部分等分為m份，根據(jù)相機(jī)成像原理，實(shí)際獲取到的圖像往往存在遠(yuǎn)小近大的特性，因此，取n>m。

(2)將m+n個(gè)水平子帶豎直等分為h份，形成如圖5所示若干個(gè)矩形塊，其中圖像的寬為W，高為H，所得的矩形寬皆為W/h，上半部分矩形(圖5中矩形a)的高為H/(2n), 下半部分矩形(圖5中矩形b)的高為H/(2m)。

(3)從下到上、從左到右掃描每一矩形塊并統(tǒng)計(jì)矩形內(nèi)白色秧苗像素?cái)?shù)，以圖5中矩形b為例，若矩形內(nèi)秧苗像素點(diǎn)數(shù)大于預(yù)先設(shè)置的閾值Th，計(jì)算所有秧苗像素點(diǎn)坐標(biāo)pb(xj,yj)平均值pavg_b(x,y)，計(jì)算公式為

(2)

式中xj、yj——矩形b中第j個(gè)秧苗像素點(diǎn)坐標(biāo)

N——矩形b中秧苗像素點(diǎn)總數(shù)

將pavg_b(x,y)作為候選特征點(diǎn)坐標(biāo)，存入候選特征點(diǎn)集Q中；否則，跳至下一矩形塊進(jìn)行相關(guān)判斷，提取的候選秧苗特征點(diǎn)如圖6所示。

圖5 區(qū)域分割矩形示意圖Fig.5 Schematic of area division

圖6 候選特征點(diǎn)圖Fig.6 Initial feature points image

2.3 特征點(diǎn)聚類

秧苗行中心線提取可以看作是對秧苗行特征點(diǎn)進(jìn)行直線擬合，為確定秧苗行的精確位置，需判斷當(dāng)前圖像中有效的秧苗行數(shù)，并利用每行秧苗特征點(diǎn)間的鄰近關(guān)系對特征點(diǎn)進(jìn)行準(zhǔn)確聚類。具體如下：

(1)對圖5中圖像頂端首條水平子帶進(jìn)行垂直投影，計(jì)算出該水平子帶的每一列中白色像素點(diǎn)數(shù)s(j)，并求出子帶的平均白色像素?cái)?shù)savg。

(2)根據(jù)s(j)與savg的關(guān)系，確定升降點(diǎn)位置[24]。如果s(j-1)savg>s(j+1)，記點(diǎn)j為下降點(diǎn)pd。

(3)計(jì)算升降點(diǎn)pu和pd之間的橫向距離d(k)和距離平均值davg，若d(k)

(4)設(shè)第i行起始點(diǎn)的坐標(biāo)為pri(x,y)，在圖7中定義以pri(x,y)為圓心，R為半徑的圓CR(pri(x,y))，如圖7a所示，圖中黑色實(shí)心圓點(diǎn)代表候選特征點(diǎn)集Q中特征點(diǎn)，如果pi(xj,yj)在圓CR(ai(xj,yj))內(nèi)，則表明點(diǎn)pi(xj,yj)屬于pri(x,y)所在的秧苗行，將pi(xj,yj)點(diǎn)作為當(dāng)前點(diǎn)并把該點(diǎn)坐標(biāo)存入特征點(diǎn)子集Qi中，以pi(xj,yj)為圓心構(gòu)建上述圓，繼續(xù)掃描其他特征點(diǎn)。

(5)采用步驟(4)中方法掃描該行所有特征點(diǎn)，直至當(dāng)前圓內(nèi)除圓心外不存在其他特征點(diǎn)時(shí)，結(jié)束該行特征點(diǎn)掃描，輸出該行特征點(diǎn)子集Qi。

(6)重復(fù)步驟(4)、(5)，將候選特征點(diǎn)集Q中所有特征點(diǎn)均歸類，圖6中候選秧苗行特征點(diǎn)的聚類結(jié)果如圖7b所示，不同顏色代表不同秧苗行的特征點(diǎn)聚類結(jié)果。

圖7 特征點(diǎn)聚類Fig.7 Feature points clustering

2.4 中心線提取

秧苗行特征點(diǎn)選取準(zhǔn)確與否直接影響提取秧苗行中心線準(zhǔn)確性，但候選特征點(diǎn)會(huì)因田間環(huán)境復(fù)雜而存在噪聲點(diǎn)，且遠(yuǎn)離秧苗行中心線的噪聲點(diǎn)對檢測結(jié)果影響較大。由于秧苗為行種植模式，每行有條趨勢線[25]，通過趨勢線與特征點(diǎn)的位置關(guān)系，設(shè)置合適的距離閾值，對已聚類的各行秧苗行特征點(diǎn)進(jìn)行篩選，去除離趨勢線遠(yuǎn)的特征點(diǎn)。

圖8 特征點(diǎn)篩選示意圖Fig.8 Schematic of feature point filtering

基于趨勢線篩選特征點(diǎn)的示意圖如圖8所示，其中趨勢線方程為y=ax+b，設(shè)置合適的距離閾值c，可以確定兩條與趨勢線平行的直線，兩直線構(gòu)成的條形區(qū)域以外的特征點(diǎn)為遠(yuǎn)距離點(diǎn)，如圖8中紅色標(biāo)記的特征點(diǎn)，將被移出特征點(diǎn)集，從而獲取最終的秧苗行特征點(diǎn)集，具體算法流程如下：

(1)連接第i行秧苗特征點(diǎn)集Qi中起點(diǎn)與終點(diǎn)，得到第i行秧苗的趨勢線Li，趨勢線方程為y=aix+bi。

(2)統(tǒng)計(jì)第i行特征點(diǎn)集Qi中特征點(diǎn)的總數(shù)t，計(jì)算特征點(diǎn)集Qi中每一個(gè)特征點(diǎn)到趨勢線Li的距離和距離閾值ci。

(3)

式中i——秧苗行數(shù)

j——Qi的第j個(gè)特征點(diǎn)

di,j——第i行第j個(gè)特征點(diǎn)到趨勢線的距離

(3)逐點(diǎn)掃描第i行特征點(diǎn)，若di,j>ci，則將該特征點(diǎn)從Qi中去除，反之不處理，直至該行所有特征點(diǎn)掃描完畢。

(4)重復(fù)步驟(1)～(3)直至圖像中所有行特征點(diǎn)篩選完畢，輸出最終的特征點(diǎn)集Qi。

圖9 特征點(diǎn)篩選結(jié)果Fig.9 Result of feature point filtering

圖9是基于趨勢線與特征點(diǎn)位置篩選后秧苗行特征點(diǎn)圖，對比圖7b，可見遠(yuǎn)距離的特征點(diǎn)被有效去除。最小二乘法通過最小化誤差的平方和尋找數(shù)據(jù)的最佳函數(shù)匹配，尤其對數(shù)據(jù)量較少時(shí)，擬合精度高、速度快。因此，本文采用最小二乘法對圖9中每一行特征點(diǎn)進(jìn)行直線擬合，擬合結(jié)果如圖10所示，圖中紅色的直線即為提取的秧苗行中心線。綜上所述，基于分區(qū)域特征點(diǎn)聚類的秧苗行中心線提取方法流程如圖11所示。

圖10 秧苗行中心線檢測結(jié)果Fig.10 Result of seedling rows detection

圖11 基于分區(qū)域特征點(diǎn)聚類的秧苗行中心線提取流程圖Fig.11 Flow chart of seedling row centerlines detection based on sub-regional feature points clustering

3 實(shí)驗(yàn)結(jié)果與分析

為判斷秧苗行中心線提取的準(zhǔn)確性，人工擬合出秧苗行中心線作為衡量基準(zhǔn)，以判斷算法提取出的秧苗行中心線的準(zhǔn)確率。人工擬合是在實(shí)驗(yàn)圖像上通過手工繪制秧苗行直線并取其平均值作為秧苗行的參考直線，人工擬合的秧苗行中心線與實(shí)際提取的秧苗行中心線之間的夾角為誤差角，該誤差角用來衡量算法獲取行線的準(zhǔn)確率，為滿足農(nóng)機(jī)農(nóng)藝要求，當(dāng)誤差角大于5°，則視提取的秧苗行線無效，以此計(jì)算識(shí)別率。

于2018年6月在安徽農(nóng)業(yè)大學(xué)郭河實(shí)驗(yàn)基地采集秧苗圖像，圖像自動(dòng)采集裝置如圖12所示，主要包括便攜式計(jì)算機(jī)、帶有USB接口的MV-CA050-20UM/C型相機(jī)和插秧機(jī)移動(dòng)平臺(tái)，相機(jī)安裝在插秧機(jī)移動(dòng)平臺(tái)的前方，距離地面高度為1.5 m，相機(jī)中心線相對移動(dòng)平臺(tái)縱軸中心線向下傾斜15°，獲得的圖像通過USB線傳輸?shù)接?jì)算機(jī)，并保存為2 560像素×1 920像素JPG格式的24位彩色圖像。為減小圖像處理的計(jì)算量，實(shí)驗(yàn)中圖像分辨率縮小為320像素×237像素，所用計(jì)算機(jī)的CPU為Intel Core i7-7700HQ，主頻為2.8 GHz，內(nèi)存8 GB，操作系統(tǒng)為Windows 10 64位系統(tǒng)，圖像處理所用的軟件是Microsoft Visual Studio 2017和OpenCV 3.0。

圖12 圖像采集裝置Fig.12 Equipment of image acquisition1.便攜式計(jì)算機(jī) 2.工業(yè)相機(jī)

共采集了3種不同天氣下若干數(shù)量的圖像，實(shí)驗(yàn)中以300幅圖像進(jìn)行秧苗行中心線檢測測試，其中晴天天氣下秧苗圖像120幅，多云天氣下秧苗圖像90幅，陰天天氣下秧苗圖像90幅。將本文算法分別與標(biāo)準(zhǔn)Hough變換和隨機(jī)Hough變換進(jìn)行比較和分析，標(biāo)準(zhǔn)Hough變換和隨機(jī)Hough變換擬合秧苗行中心線的大致步驟如下: 采用2G-R-B提取秧苗超綠特征圖像，并進(jìn)行 Otsu二值分割，以及膨脹、腐蝕等形態(tài)學(xué)操作，將秧苗從土壤背景中分割出來，采用Hough變換或隨機(jī)Hough變換對分割后的秧苗二值圖像進(jìn)行直線擬合。其中，標(biāo)準(zhǔn)Hough變換算法中參數(shù)θ∈[-90°,90°]、ρ∈[-398,398]， 且θ的量化精度為1°；隨機(jī)Hough變換算法中最短直線的投票閾值為3，本文算法中參數(shù)m=10、n=20、h=50、Th=25、R=20。表2為3種算法的識(shí)別率、準(zhǔn)確率和耗時(shí)數(shù)據(jù)。

表2 不同算法的識(shí)別率、準(zhǔn)確率和耗時(shí)的比較Tab.2 Comparison of recognition rate, accuracy rate and time-consuming

由于標(biāo)準(zhǔn)Hough變換是對全局的二值圖像進(jìn)行參數(shù)空間的映射，圖像信息量大，降低了檢測速度，且向外分散生長的秧苗葉片和噪聲像素被檢測為秧苗像素，導(dǎo)致秧苗行中心線檢測存在偏差。相比標(biāo)準(zhǔn)Hough，隨機(jī)Hough變換采用多對一的映射，降低了運(yùn)算量，提高了檢測精度，但在處理復(fù)雜圖像時(shí)隨機(jī)采樣會(huì)引入無效采樣，從而降低算法的性能。而本文算法在二值圖像的基礎(chǔ)上通過識(shí)別、聚類秧苗行特征點(diǎn)來對秧苗行進(jìn)行定位，克服了最小二乘法受干擾點(diǎn)影響的不足，又充分利用了最小二乘法耗時(shí)低于Hough變換的優(yōu)勢，通過表2可知，該算法可以成功檢測出300幅秧苗圖像中所有秧苗行線，識(shí)別率達(dá)95.2%，且具有耗時(shí)少、準(zhǔn)確率高的優(yōu)勢。

圖13a為選取的3幅不同光照條件下水稻秧苗圖像，從左至右分別為晴天天氣下水稻秧苗圖像、多云天氣下水稻秧苗圖像、陰天天氣下水稻秧苗圖像。圖13b為經(jīng)過2.1節(jié)和2.2節(jié)中算法處理后提取候選特征點(diǎn)圖像，由圖可以看出，提取出的特征點(diǎn)可以很好地代表秧苗行的趨勢。用本文提出的算法及標(biāo)準(zhǔn)Hough變換算法、隨機(jī)Hough變換算法對圖13b所示的特征點(diǎn)圖進(jìn)行秧苗行中心線檢測。

圖13c為采用標(biāo)準(zhǔn)Hough變換算法對圖13b中特征點(diǎn)圖像進(jìn)行秧苗行中心線檢測的結(jié)果，紅色直線是實(shí)際檢測的行中心線，藍(lán)色直線是人工擬合的參考線。由圖中實(shí)驗(yàn)檢測結(jié)果與參考直線的夾角可以看出，在秧苗行周圍背景噪聲較少時(shí)，標(biāo)準(zhǔn)Hough變換可以成功檢測出秧苗行，但是當(dāng)秧苗行周圍背景噪聲過多和光照不均時(shí)，特征點(diǎn)較為分散，此時(shí)標(biāo)準(zhǔn)Hough的檢測結(jié)果將出現(xiàn)較大偏差，甚至?xí)霈F(xiàn)秧苗行的漏檢測和誤檢測現(xiàn)象，其平均角度偏差為2.847°。

圖13d為基于特征點(diǎn)圖像采用隨機(jī)Hough變換進(jìn)行秧苗行中心線提取的結(jié)果，由圖中實(shí)驗(yàn)檢測結(jié)果與參考直線的夾角可以看出，隨機(jī)Hough變換算法在漏檢測和誤檢測方面較Hough變換有很大改善，但當(dāng)秧苗行周圍泥土過多和光照不均時(shí)，檢測結(jié)果依然存在較大偏差，其平均角度偏差為1.803°。

針對以上2種算法對于特征點(diǎn)分散的圖像檢測結(jié)果不夠準(zhǔn)確的問題，本文對候選特征點(diǎn)進(jìn)行進(jìn)一步篩選，如圖13e所示，篩選后的特征點(diǎn)更好地描述了秧苗行中心線的位置。圖13f顯示了本文算法的秧苗行中心線檢測結(jié)果，其平均角度偏差為0.794°，對于秧苗行周圍有明顯土壤雜物，以及光照不均的圖像，本文算法仍然可以成功檢測出秧苗行中心線。

圖13 不同算法檢測結(jié)果的比較Fig.13 Comparisons of detection lines by different algorithms

綜合上述秧苗行中心線檢測效果分析，所提算法的準(zhǔn)確性和檢測速度高于標(biāo)準(zhǔn)Hough變換算法和隨機(jī)Hough變換算法。但本文選取的實(shí)驗(yàn)對象是水田秧苗，沒有驗(yàn)證算法的作物行檢測的適應(yīng)性，且本文算法處理一幅320像素×237像素的彩色圖像平均耗時(shí)為0.467 s，不能滿足實(shí)際攝像機(jī)幀率為25 f/s的需要。因此在后續(xù)的研究中需要開展算法對不同類型農(nóng)作物和作物不同生長階段的適應(yīng)性研究，并改進(jìn)、優(yōu)化算法，提高算法的適應(yīng)性和運(yùn)算速度。

4 結(jié)論

(1)采用2G-R-B特征因子和Otsu法對彩色圖像進(jìn)行預(yù)處理，基于秧苗像素點(diǎn)的位置分布，進(jìn)行分區(qū)域的秧苗行特征點(diǎn)提取和聚類，通過最小二乘法擬合秧苗行中心線，可以有效解決秧苗行周圍背景噪聲過多和光照不均的問題，具有較好的魯棒性。

(2)通過對不同天氣下的300幅圖像進(jìn)行測試，利用本文算法提取秧苗行中心線的準(zhǔn)確率為95.6%、識(shí)別率為95.2%，高于標(biāo)準(zhǔn)Hough變換和隨機(jī)Hough變換算法，同時(shí)，算法處理速率也有明顯提升。

(3)采用各算法對隨機(jī)選取的不同天氣下3幅秧苗圖像提取秧苗行中心線，并與人工擬合秧苗行中心線對比，結(jié)果顯示，標(biāo)準(zhǔn)Hough變換、隨機(jī)Hough變換和本文算法的平均偏差角度分別為2.847°、1.803°和0.794°，表明本文算法能夠準(zhǔn)確提取秧苗行中心線，可靠性較高。