肖菲菲, 蔣蘋, 胡文武, 廖榮華, 張丹慧, 金生

(湖南農(nóng)業(yè)大學機電工程學院, 長沙 410128)

目前，精準農(nóng)業(yè)是農(nóng)業(yè)智能化發(fā)展的研究熱點之一，其中田間自動導航技術(shù)是完成精準農(nóng)業(yè)的關(guān)鍵環(huán)節(jié)[1]。農(nóng)業(yè)機械的自動導航技術(shù)現(xiàn)已被廣泛應(yīng)用于噴藥、施肥和除草等農(nóng)業(yè)作業(yè)過程[2-3]，該項技術(shù)不僅可以降低勞動強度和農(nóng)業(yè)生產(chǎn)成本，還可提高農(nóng)機作業(yè)效率。而在農(nóng)業(yè)機械的自動導航過程中，農(nóng)作物的苗帶識別和定位是核心技術(shù)之一，在很大程度上影響農(nóng)業(yè)機械的導航效果和作業(yè)的精準度[4-7]。苗帶是指根據(jù)農(nóng)作物種植方式，在同一方向上農(nóng)作物連續(xù)性的集合，在作物生長初期，苗帶為條狀；在作物生長后期，作物封行后，苗帶為塊狀。苗帶識別是指通過某種方法區(qū)分農(nóng)作物不同生長狀態(tài)下的苗帶區(qū)域與相鄰苗帶間的行間間隙所在位置。通過苗帶識別，在可行走區(qū)域規(guī)劃行走路徑，使機械行走過程中避免碾壓苗帶，實現(xiàn)田間管理機械自主低傷苗作業(yè)，對促進精準農(nóng)業(yè)發(fā)展有重要意義[8]。

當農(nóng)業(yè)機械行駛在環(huán)境復雜的農(nóng)田中，為了避免障礙物和壓傷農(nóng)作物，需要通過傳感器進行感知識別，當前常用于苗帶檢測識別的傳感器有視覺傳感器和激光雷達等[9-14]。王姍姍等[15]通過視覺傳感器檢測水稻秧苗，采用基于特征點鄰域 Hough變換算法完成水稻苗中心線提取，獲得平均識別準確率為92%，其識別率較好，但響應(yīng)速度慢，只適用于速度低于1.5 m·s-1。張博立等[16]利用機器視覺提取玉米苗帶中心線，采用動態(tài)網(wǎng)格與分區(qū)域聚類算法相結(jié)合進行識別，獲得平均玉米苗帶識別準確率為93.4%，有效地處理了雜草對識別效果的影響，但由于圖像相關(guān)處理問題，導致中心線檢測失敗,且目前僅適用于識別生長初期玉米作物行。姜國權(quán)等[17]通過粒子群優(yōu)化算法獲取麥田作物行的特征點，利用最小二乘法進行線性擬合，運行速度相對傳統(tǒng)Hough算法有很大提高，但并不適用于處理雜草較多的圖像。Hu等[18]利用二維伺服平臺，通過激光雷達掃描獲得水稻苗帶點云信息，設(shè)計基于斜率虛化均值迭代聚類算法，采用斜率虛化處理點云數(shù)據(jù)，通過區(qū)域迭代聚類獲得苗帶特征，利用最小二乘法進行苗帶軌跡線性化擬合，經(jīng)試驗獲得橫向最大誤差為45 mm，中位點偏差為7 mm，但未考慮激光掃描偏移問題，且存在重復運行算法時出現(xiàn)聚類失敗的情況。

由以上研究可知，視覺傳感器成本較低，可實時獲取作物行相對位置，因此廣泛用于作物苗帶識別，但傳感器易受外界環(huán)境因素的影響，且數(shù)據(jù)處理相對繁雜，部分系統(tǒng)并不能滿足真實的農(nóng)田環(huán)境。而激光雷達具有易操作、分辨率高、抗干擾性強等特點，更宜用于農(nóng)田環(huán)境，但目前大多研究集中于室內(nèi)環(huán)境或已知結(jié)構(gòu)的室外環(huán)境，對于水田的作物識別研究相對較少，且并未考慮掃描偏移對識別結(jié)果的影響。

綜上，本文提出了基于偏移補償模型的極值點聚類苗帶識別算法，重點研究激光掃描偏移對識別結(jié)果的影響。首次提出了苗帶識別的相關(guān)定義，根據(jù)獲得的苗帶點云線特征，采用極值檢測算法獲取聚類要素點，根據(jù)苗帶特征，分析激光測距儀的掃描特點，建立偏移補償模型，對獲取的聚類要素點進行偏移矯正，最后通過最小二乘法進行線性擬合，獲得苗帶軌跡信息。提高了軌跡重合精度，增強了算法運行效率，獲得的苗帶軌跡參數(shù)可為田間管理機械低傷苗直線行駛提供理論依據(jù)。

1 材料與方法

1.1 系統(tǒng)構(gòu)成及工作原理

導航定位控制系統(tǒng)由植保機行駛控制單元(electronic control unit)和苗帶識別系統(tǒng)構(gòu)成，其中苗帶識別系統(tǒng)包含激光雷達、舵機調(diào)平單元與慣性測量模塊，如圖1所示。

注：1—施肥斗; 2—發(fā)動機; 3—苗帶識別系統(tǒng)； 4—植保機行車控制單元；5—藥箱；6—噴桿。

激光雷達安裝于舵機平臺上，慣性測量模塊采集姿態(tài)數(shù)據(jù)，舵機調(diào)平單元通過采集慣性測量模塊的數(shù)據(jù)，控制舵機進行俯仰、翻滾、旋轉(zhuǎn)動作，進而實現(xiàn)激光雷達底座的俯仰角、偏航角、橫滾角可控調(diào)節(jié)，并通過CAN接口與植保機ECU單元進行數(shù)據(jù)交換。ECU單元采集激光雷達掃描的點云數(shù)據(jù)與舵機平臺角度數(shù)據(jù)后，通過處理獲得水稻苗帶信息，并識別可行駛區(qū)間，控制田間管理機械行駛于苗帶行間間隙中，減少作業(yè)過程中對水稻的損害。

1.2 苗帶識別系統(tǒng)硬件與數(shù)據(jù)采集軟件

水稻苗帶識別系統(tǒng)的基本結(jié)構(gòu)(圖2)由三部分構(gòu)成，分別為激光雷達、二維舵機調(diào)平單元和支撐平臺。本研究選取日本HOKUYO的UTM-30LX二維激光雷達，測量范圍為100～30 000 mm，100～10 000 mm的測距誤差為±30 mm，掃描角度范圍-135°～135°，掃描周期40 Hz，角度分辨率為0.25°。舵機調(diào)平單元采用達盛舵機RDS3115，實現(xiàn)X、Y軸的水平調(diào)節(jié)控制，在履帶底盤底部安裝旋轉(zhuǎn)調(diào)節(jié)舵機，實現(xiàn)X、Y軸的旋轉(zhuǎn)調(diào)節(jié)，完成激光雷達俯仰角、橫滾角、航向角的調(diào)節(jié)控制。為實現(xiàn)舵機角度初始化校準和激光雷達中位調(diào)整，采用核心處理器STM32F103配合MPU6050六軸姿態(tài)傳感器，實現(xiàn)舵機角度調(diào)節(jié)。

注：1—控制單元; 2—支撐平臺; 3—旋轉(zhuǎn)調(diào)節(jié)舵機； 4—激光雷達；5—俯仰調(diào)節(jié)舵機；6—翻滾調(diào)節(jié)舵機。

激光雷達通過軟件UrgBenri Standard V1.8.1進行數(shù)據(jù)采集，數(shù)據(jù)存儲格式為ubh，使用Lenovo CPU1.60GHz在Win10環(huán)境下運用Matlab R2016b進行數(shù)據(jù)分析與算法研究。

1.3 點云數(shù)據(jù)獲取與苗帶軌跡提取分析

1.3.1點云數(shù)據(jù)獲取為確定激光束發(fā)出位置，確保激光雷達采集數(shù)據(jù)的準確性，本研究在數(shù)據(jù)采集前通過最小二乘平差法求得真實距離與激光雷達所測距離之間的數(shù)學關(guān)系[19-20]對測量儀進行標定，以消除激光雷達自身掃描誤差。苗帶數(shù)據(jù)采集于2020年7月湖南農(nóng)科院試驗田中，此時水稻處于移栽后15 d，當天天氣預報風速為2 km·h-1，環(huán)境溫度32 ℃，目測水稻葉片處于靜止狀態(tài)。采用卷尺測量水稻高度、行株距和激光雷達距離水田平面的高度等參數(shù)，采用5組連續(xù)15株水稻苗信息求均值的方法，獲得水稻高度和行距株距。通過測試計算，得到水稻苗與水面高度平均值為360 mm，平均行距為220 mm，激光雷達掃描中心點到水田水平面的垂直距離為960 mm。在實際稻田環(huán)境數(shù)據(jù)采集過程中，二維舵機調(diào)平單元的俯仰掃描初始角度為9°，步進角度0.5°，終止角度17°，激光雷達掃描角度范圍為-16°～16°，掃描寬度對應(yīng)5行苗帶，以此試驗條件獲取苗帶點云數(shù)據(jù)。

1.3.2坐標轉(zhuǎn)化激光雷達測距顯示參數(shù)主要為當前激光雷達掃描的起止角、步長、掃描線數(shù)，導出數(shù)據(jù)為距離與角度的極坐標數(shù)據(jù)，系統(tǒng)參數(shù)還包括當前激光雷達的俯仰角度數(shù)據(jù)。圖3為苗帶識別系統(tǒng)工作時的相關(guān)參數(shù)，以水田水平面為參考，以激光雷達激光旋轉(zhuǎn)中心為原點，垂直方向為Z軸，激光雷達初始光軸與機具直行重合方向為Y軸，垂直前進方向為X軸，激光雷達俯仰角度θ，激光測距長度為L，對應(yīng)激光線掃描線角度為γ，Y軸對應(yīng)激光雷達中心位置偏移初始角度為γ0。為方便數(shù)據(jù)分析，對激光極坐標數(shù)據(jù)進行三維坐標轉(zhuǎn)換，基于上述坐標轉(zhuǎn)換模型獲得被測點相對激光雷達當前位置的相對直角坐標表達式(1)。

圖3 掃描系統(tǒng)相關(guān)參數(shù)

(1)

1.3.3苗帶軌跡獲取方法為獲得最終苗帶軌跡，本研究系統(tǒng)將激光雷達和慣性測量模塊采集的數(shù)據(jù)融合，獲得被測目標點的位置坐標，運用Matlab R2016b進行數(shù)據(jù)分析與算法研究。根據(jù)農(nóng)業(yè)機械作業(yè)環(huán)境和激光雷達點云數(shù)據(jù)特征，確定本文數(shù)據(jù)處理流程為數(shù)據(jù)預處理、點云極值點提取、苗帶特征點聚類和軌跡獲取與評價四個部分。數(shù)據(jù)預處理主要包括點云數(shù)據(jù)格式轉(zhuǎn)換、數(shù)據(jù)去噪和數(shù)據(jù)坐標轉(zhuǎn)換等步驟；點云極值點提取主要實現(xiàn)苗帶特征點的提??；偏移補償均值聚類主要進行苗帶特征點聚類；軌跡獲取與評價主要完成聚類點的線性擬合，并對獲取的軌跡線進行有效性評價。

1.4 基于偏移補償模型的極值點聚類

1.4.1苗帶模型建立由于田間管理機械直線行駛僅需要獲取以水平面為參考的平面坐標，因此，對三維數(shù)據(jù)平面化，獲得基于激光雷達為原點的XY坐標定位數(shù)據(jù)。

(2)

為了方便計算苗帶位置，獲得苗帶函數(shù)，根據(jù)激光雷達掃描測距方法建立理想苗帶模型，如圖4所示，長條框為苗帶，激光測距儀距離植株距離b，苗帶寬度d1,相鄰苗帶間距為d2。為驗證模型的正確性，使用5個紙箱進行模擬試驗，中心間距500 mm，UrgBenri Standard V1.8.1采集圖像。由圖像對比可知 ，建立的模型與模擬實驗獲得的點云線相對應(yīng)，圖像在X軸方向都具有周期性，且為具有規(guī)律性的波浪峰，模型結(jié)果成立。

圖4 數(shù)據(jù)采集模型

從數(shù)據(jù)采集模型(圖4)可知，紅色線為激光雷達采集的點云數(shù)據(jù)線，圖像在X方向具有周期性，周期為d1+d2，且圖像關(guān)于Y軸對稱，根據(jù)模型圖和已知參數(shù)計算模型函數(shù)，函數(shù)共分為三段，計算得出函數(shù)關(guān)系式(3)。

(3)

1.4.2建立偏移補償模型根據(jù)建立的模型可知，由于苗帶的物理遮擋，導致部分位置的點云數(shù)據(jù)無法正確采集。在掃描中間苗帶時(圖5)，掃描線經(jīng)過Q點后，右側(cè)面的點云數(shù)據(jù)由于遮擋無法采集，使得在QN線上的點云數(shù)據(jù)產(chǎn)生了偏移，所以需要建立偏移補償模型進行矯正。在激光掃描范圍內(nèi)，以原點為中心左右兩側(cè)最高峰為大波峰，最低峰稱為小波峰，Y小于小波峰以下的數(shù)據(jù)相對較為準確，偏移量較小，可通過聚類算法直接獲得苗帶聚類點；Y在小波峰和大波峰之間的數(shù)據(jù)偏移較為嚴重，需對獲得的苗帶聚類點進行偏移矯正?；谝陨蠁栴}，建立偏移補償模型，如圖5所示，在縱向小波峰到大波峰范圍(GA≤Y≤NA)內(nèi)，獲得偏移量ΔX的矯正范圍。

圖5 偏移補償模型

在EM線上，通過苗帶模型函數(shù)，在點E、M處分別得到如下公式。

(4)

(5)

(6)

在QN線上，通過苗帶模型函數(shù)，在點H、N處分別得到如下公式。

t×tanγ1=GA=MB

(7)

(8)

(9)

則苗帶間隙最大偏移量計算如下。

(10)

因此，偏移量范圍為0≤ΔX,基于此建立的偏移模型在同一波中，取低于小波峰的點云偏移量趨向0，在小波峰與大波峰之間的點云偏移量為ΔX，獲得矯正模型(圖6)，表達如式(11)所示。

圖6 矯正模型

(11)

1.4.3聚類要素點的獲取根據(jù)建立的模型(圖4)可知，每條點云數(shù)據(jù)線在X方向都存在規(guī)律性的波浪峰，而波峰波谷所在位置與苗帶和苗帶行間間隙相關(guān)，然而波峰波谷位置又對應(yīng)函數(shù)中的極值所在位置，因此，可通過計算函數(shù)的極值所在位置得到與苗帶位置相關(guān)信息。本研究采用極大值檢測算法進行極值檢測，尋找苗帶位置，由式(12)確定。

s=find{diff[sign(diff(w))]<0}+1

(12)

式中，diff指差分運算，sign為符號函數(shù)，find表示查找點云數(shù)據(jù)線中滿足條件的數(shù)據(jù)位置[21]。

當采樣點滿足式(12)時，其點云數(shù)據(jù)在鄰域范圍內(nèi)呈先增加后減小的趨勢，即局部極大值點，獲得波峰位置；將原數(shù)據(jù)的縱坐標值取反，再利用函數(shù)查找波谷位置。因農(nóng)田環(huán)境復雜，水稻葉片具有下垂交疊等特征，且此特征也可作為判別苗帶位置相關(guān)因素，為避免信息缺失，所以閾值范圍取最小值到最大值，使用Matlab 繪制多角度點云數(shù)據(jù)線的波峰波谷圖，進而獲得苗帶聚類特征信息。

1.4.4改進K-means均值聚類均值聚類算法的基本思想是以空間中k個點為中心聚類，對最接近他們的對象進行歸類，通過迭代的方法不斷調(diào)整各分類中心位置，直到收斂[22]。聚類算法中聚類數(shù)量與初始聚類中心的選擇會直接影響聚類結(jié)果，通過改進K-means聚類算法提高苗帶點云數(shù)據(jù)的聚類效果。通過矯正函數(shù)獲得最佳聚類個數(shù)K,具體流程如圖7。初始中心的選擇通過密度公式(13)確定，在r∈(0,d1+d2)范圍內(nèi)，定義對應(yīng)運算區(qū)域的聚類點位置X，計算xi內(nèi)的點密度ρ(xi),當密度達到最大時，在對應(yīng)范圍內(nèi)采用歐式距離度量方式選取距離xi最遠的點作為第(i+1)類初始聚類中心xi+1，并且除去xi、xi+1，以此類推，重復上述步驟再獲得下一個聚類中心，不斷迭代，直至獲得第K個聚類中心結(jié)束。

圖7 獲取聚類點數(shù)流程

ρ(xi)={x∈X|dist(xi,x)≤r}

(13)

1.4.5苗帶軌跡處理為便于獲取苗帶軌跡，獲得更多用于線性擬合的聚類點，同時降低數(shù)據(jù)采樣維度，采用縱向等區(qū)間劃分增加數(shù)據(jù)獲取量，獲得更為準確的軌跡參數(shù)，對Y軸向范圍進行n等分區(qū)間聚類。通過分層聚類后，獲得分段范圍內(nèi)的聚類點。

(14)

式中，k為單次聚類數(shù)量，(x,y)為聚類點對應(yīng)的坐標位置,n為Y軸等分處理的次數(shù)，從而獲得nk個聚類點。則對應(yīng)苗帶行間間隙點為(xv,yv),矩陣為O(n)。

(15)

(16)

(17)

(18)

基于偏移模型進行矯正獲得聚類矯正后的矩陣T。

(19)

從而獲得對應(yīng)Y軸向的s組聚類點矩陣M(s)。

(20)

對上述矩陣坐標點進行線性擬合獲取對應(yīng)苗帶行間間隙的軌跡方程。對聚類點坐標矩陣通過最小二乘法進行線性擬合，獲得Ak和Bk的值，從而獲取苗帶對應(yīng)的理論軌跡。

(21)

(22)

對應(yīng)苗帶間隙軌跡表達式如下。

y=Akx+Bk

(23)

1.5 苗帶軌跡評價

參考文獻[23]的評價指標，在獲得理論苗帶軌跡后，與對應(yīng)水稻苗帶真實軌跡進行對比，采用測量軌跡與真實軌跡平行度和中心偏差距離兩個指標評價算法的平行度(式24)和中心偏差(式25)。

Err_max=|xk-x′k|max,yk∈(a～b)

(24)

Err_mid=|xk-x′k|max,yk=(a+b)/2

(25)

式中，xk為對應(yīng)第k行x的真實值，x′k為對應(yīng)第k行x的測量值。

2 結(jié)果與分析

2.1 測量儀標定誤差分析

通過5組不同距離的數(shù)據(jù)，對激光測距儀進行標定，標定后實際距離與激光雷達測量結(jié)果見表1?？梢钥闯觯S著激光雷達所測距離由0.645增加到14.389 m,數(shù)據(jù)標定前與標定后的誤差也逐漸增加，但標定后的誤差明顯低于標定前的。當被測物體距激光雷達14.389 m時，標定后的誤差為0.011 m，相對標定前的誤差0.037 m明顯減少。說明測距儀使用前有必要進行標定。同樣，完成二維舵機平臺與激光雷達的聯(lián)合標定。

表1 雷達標定前后測量值誤差

2.2 試驗數(shù)據(jù)模型構(gòu)建

農(nóng)田環(huán)境中，作物通常采用彼此平行的種植方式[24]，采集試驗區(qū)域?qū)?yīng)有效測量范圍的水平投影數(shù)據(jù)。以激光雷達為坐標原點，為消除側(cè)偏導致的方位差，激光雷達前視方向與中間列位置基本平行，為滿足前視行駛需要，設(shè)定前視水平投影距離為3 030 mm，前視水平投影探測有效范圍為2 000～3 030 mm，有效長度為1 030 mm，左右水平投影探測有效寬度為-500～+500 mm。測試環(huán)境的真實苗帶軌跡如式(26)。

(26)

則測試環(huán)境的可行駛苗帶行間間隙軌跡如式(27)。

(27)

2.3 苗帶點云數(shù)據(jù)分析

2.3.1苗帶三維點云顯示及分析通過Matlab R2016b對激光測距儀采集的點云數(shù)據(jù)進行預處理，繪制苗帶軌跡三維立體圖，如圖8所示。

圖8 MATLAB轉(zhuǎn)換后的點云數(shù)據(jù)

根據(jù)田間采集的實況與三維點云圖外觀對比分析可知，對應(yīng)的苗帶信息清晰有規(guī)律，且相鄰苗帶行間間隙的初始軌跡鮮明。箭頭標明了苗帶行間間隙軌跡所在位置，但受葉片傾斜折彎和交叉下垂以及激光雷達掃描方式的影響，行間間隙寬度較窄，且隨著前視距離的增大間隙軌跡不易識別。多角度掃描軌跡線在靠近X軸有五個突出的波谷，為激光雷達掃描的第一行水稻苗所呈現(xiàn)出的形狀，且每根點云數(shù)據(jù)線都存在這種規(guī)律性的波浪峰狀特征。

2.3.2點云極值分析通過對每根點云線進行極值檢測，獲得波峰波谷數(shù)據(jù)點位置，獲取聚類要素點，如圖9所示。

圖9 波峰波谷數(shù)據(jù)點

注：圖中不同顏色的圈圈代表不同的分類，*代表此類型的聚類點。

表2 最終聚類中心

2.3.4苗帶軌跡提取通過聚類后進而獲得水平矩陣R(n)。

(28)

(29)

則對應(yīng)苗帶行間間隙點為(xv,yv),矩陣為O(n)。

(30)

(31)

(32)

對應(yīng)獲得水平聚類點矩陣Q。

根據(jù)實際采集情況，基于偏移模型計算獲得偏移量范圍為0≤ΔX≤174 mm，獲得聚類矯正后的矩陣T。

(33)

通過平面投影處理，獲得如圖11的分段均值聚類分布圖，其中黑點為聚類中心點，紅點為矯正點。

注：圖中不同顏色的圈圈代表不同的分類，*代表此類型的聚類點。

獲取Y軸向?qū)?yīng)的4組聚類點矩陣M1～M4。

(34)

對應(yīng)上述縱向4組苗帶行間間隙聚類點，在Matlab下進行線性擬合處理，獲得軌跡矩陣如式(35)，對應(yīng)軌跡如圖12所示。

圖12 苗帶行間間隙軌跡對比

(35)

2.4 軌跡結(jié)果分析

在y∈(2 000,3 000) 范圍內(nèi)，計算獲得軌跡與實際軌對應(yīng)橫向最大偏差Err_max和中位點偏差Err_mid。

Err_max=[16111614]

(36)

Err_mid=[6461]

(37)

表4分別統(tǒng)計了采用兩種不同算法對苗帶點云數(shù)據(jù)處理后的X軸最大橫向偏差、中位點最大偏差和平均耗時。

表4 算法檢測偏差及檢測時間對比

從表5中數(shù)據(jù)可得出，在前視距離2 000～3 000 mm范圍內(nèi),通過斜率虛化聚類算法對點云數(shù)據(jù)進行處理，獲得軌跡與實際苗帶軌跡X軸橫向最大偏差為45 mm, 中位點最大偏差7 mm，平均耗時1.64 s。本文采用基于偏移補償模型的極值點聚類苗帶識別算法獲得X軸橫向最大偏差控制在16 mm以內(nèi)，中位點最大偏差為6 mm，運算耗時1.50 s，相比斜率虛化聚類算法X軸向最大偏減小了29 mm，中心點最大偏差減小1 mm，耗時減少0.14 s，檢測效果明顯提高，所檢測苗帶行間間隙軌跡與實際軌跡具有較好的重合度，滿足水田管理機械免壓苗行駛的需求(輪寬10 cm)。

3 討論

苗帶軌跡的識別是實現(xiàn)田間自動導航技術(shù)關(guān)鍵?，F(xiàn)有對作物的檢測大多針對旱地作物,如蔬菜、玉米、小麥等，對于水田的作物識別研究相對較少。本文首次提出了苗帶識別的概念，利用二維激光雷達與舵機調(diào)平單元相結(jié)合，以獲得三維掃描功能完成苗帶點云數(shù)據(jù)獲取，進行水稻苗帶識別研究，避免了視覺傳感器易受外界環(huán)境影響，數(shù)據(jù)處理的算法復雜度較高等問題，重點研究激光掃描偏移對識別結(jié)果的影響，為解決激光雷達掃描偏移和提高軌跡識別精度等問題，提出了基于偏移補償模型的極值點聚類苗帶識別算法。根據(jù)獲得的點云線特征，利用極值檢測算法獲得聚類要素點，分析激光雷達的苗帶掃描特征，建立了偏移補償模型，結(jié)合改進的K-means聚類算法進行偏移矯正聚類，最后通過最小二乘進行線性擬合，從而實現(xiàn)苗帶軌跡可行駛區(qū)域的識別。

通過對前視距離為2 000～3 000 mm范圍內(nèi)的苗帶軌跡進行識別，獲得最大橫向偏差低于16 mm,中位值偏差控制在6 mm以內(nèi)。相對斜率虛化聚類算法，考慮了激光掃描偏移等問題，為提高聚類效果，改善原始K-means聚類算法中聚類數(shù)量與初始聚類中心需自行設(shè)定等問題，結(jié)合苗帶特性對其進行了改進，試驗結(jié)果表明，基于偏移補償模型的極值點聚類苗帶識別算法軌跡識別精度更高，平均運算耗時更少，算法魯棒性更強，可為田間管理機械進行低傷苗自主導航行駛提供理論依據(jù)。但本文算法進行了Y軸縱向等分聚類，在選取不同等分數(shù)量時，可能會出現(xiàn)聚類效果不佳現(xiàn)象，后期需進一步改進。文中采集數(shù)據(jù)的試驗田為機插秧的標準化水田，測量對象無缺失，測量現(xiàn)場環(huán)境良好，但實際農(nóng)田環(huán)境中有諸多不確定性因素影響，需要在后期進行進一步的研究。