• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于機(jī)器視覺的農(nóng)田地頭邊界線檢測(cè)方法

    2020-07-07 06:10:48楊鵬樹孟志軍
    關(guān)鍵詞:邊界線波谷波峰

    王 僑 劉 卉 楊鵬樹, 孟志軍

    (1.國(guó)家農(nóng)業(yè)智能裝備工程技術(shù)研究中心,北京100097;2.首都師范大學(xué)信息工程學(xué)院,北京100048)

    0 引言

    農(nóng)業(yè)生產(chǎn)機(jī)器人化是解決農(nóng)業(yè)勞動(dòng)力短缺、促進(jìn)農(nóng)業(yè)機(jī)械精準(zhǔn)化和智能化、提高農(nóng)業(yè)生產(chǎn)效率和經(jīng)濟(jì)效益的關(guān)鍵舉措。農(nóng)田作業(yè)機(jī)器人化是農(nóng)業(yè)生產(chǎn)機(jī)器人化的重要組成部分,非結(jié)構(gòu)農(nóng)田作業(yè)環(huán)境中農(nóng)機(jī)全程自主導(dǎo)航行走技術(shù)是當(dāng)前農(nóng)田作業(yè)由自動(dòng)化過渡至機(jī)器人化的關(guān)鍵技術(shù)之一。

    實(shí)現(xiàn)農(nóng)機(jī)田間作業(yè)過程中的全程自主導(dǎo)航行走,主要包括實(shí)現(xiàn)其田內(nèi)正常作業(yè)時(shí)的自主導(dǎo)航行走[1-4]、田內(nèi)遇障時(shí)的自主導(dǎo)航避障[5-8]以及地頭處的自主導(dǎo)航轉(zhuǎn)彎[9-12]三部分。隨著農(nóng)田場(chǎng)景的變化,田內(nèi)障礙物通常表現(xiàn)出隨機(jī)性、不確定性和非必然存在性,故針對(duì)農(nóng)機(jī)田內(nèi)遇障時(shí)的自主導(dǎo)航避障技術(shù)研究暫不具有迫切性。而針對(duì)農(nóng)機(jī)田內(nèi)正常作業(yè)時(shí)的自主導(dǎo)航行走,目前基于全球衛(wèi)星導(dǎo)航定位系統(tǒng)(Global navigation satellite system,GNSS)[13]的田內(nèi)直線作業(yè)導(dǎo)航技術(shù)已較成熟,相關(guān)產(chǎn)品已實(shí)現(xiàn)批量銷售。但是,農(nóng)機(jī)抵達(dá)地頭處的轉(zhuǎn)向掉頭過程,目前仍需人工操作來實(shí)現(xiàn)。因此,研究地頭處的自主導(dǎo)航轉(zhuǎn)彎技術(shù),是目前實(shí)現(xiàn)農(nóng)機(jī)田間作業(yè)過程中全程自主導(dǎo)航行走最為緊迫的任務(wù)。

    在非結(jié)構(gòu)化的農(nóng)田作業(yè)環(huán)境中,為實(shí)現(xiàn)農(nóng)機(jī)在地頭處的自主導(dǎo)航轉(zhuǎn)彎,首先要感知地頭的空間位置信息。目前,國(guó)內(nèi)外針對(duì)農(nóng)機(jī)在地頭處的轉(zhuǎn)彎路徑規(guī)劃和轉(zhuǎn)向控制進(jìn)行了較多研究,但是針對(duì)地頭空間位置信息檢測(cè)方面的研究卻較少,且相關(guān)研究主要基于機(jī)器視覺感知技術(shù)來判斷是否到達(dá)地頭[12,14-15]或完成對(duì)地頭邊界線的定位檢測(cè)[16-20]。地頭邊界線是描述地頭空間位置的關(guān)鍵參數(shù),在短距離范圍內(nèi)通常近似呈線性特征,在農(nóng)機(jī)田間作業(yè)過程中,實(shí)時(shí)探測(cè)地頭邊界線的位置,可為后續(xù)規(guī)劃地頭轉(zhuǎn)彎路徑、實(shí)施自主轉(zhuǎn)向掉頭提供關(guān)鍵的參考線或基準(zhǔn)線。文獻(xiàn)[16 -17]中分別針對(duì)土質(zhì)田埂和水泥田埂,研究了水田中田埂線的視覺檢測(cè)方法,主要依據(jù)秧苗插種之前水田內(nèi)外邊緣紋理特征的差異來獲取對(duì)應(yīng)分界線,檢測(cè)準(zhǔn)確率在98%以上,但該方法并不適用于旱田環(huán)境。對(duì)于旱田環(huán)境中的地頭邊界線,已有研究[18-20]主要依據(jù)農(nóng)田內(nèi)外的亮度或顏色差異,通過尋找跳變行確定其位置,所獲得的地頭邊界線在圖像中往往呈水平直線狀。很顯然,相關(guān)方法僅適用于檢測(cè)地頭邊界與農(nóng)機(jī)作業(yè)方向近似垂直(此時(shí)地頭邊界近似呈水平延伸狀)、較規(guī)整的地頭場(chǎng)景,對(duì)于實(shí)際農(nóng)田環(huán)境中較常見的非垂直、呈傾斜延伸狀、非規(guī)整地頭邊界類型并不適用,具有較大的應(yīng)用局限性。并且所檢測(cè)的邊界線也無法保證農(nóng)機(jī)在其內(nèi)實(shí)施自主轉(zhuǎn)向掉頭操作時(shí)安全可靠。目前,關(guān)于地頭邊界線的視覺檢測(cè)算法所涉及的地頭場(chǎng)景單一、特殊,算法的魯棒性較差,無法適應(yīng)實(shí)際農(nóng)田中復(fù)雜多變的地頭環(huán)境,且所檢測(cè)的地頭邊界線并不能直接用于農(nóng)機(jī)在地頭處掉頭轉(zhuǎn)向的基準(zhǔn)線。

    本文面向農(nóng)機(jī)實(shí)現(xiàn)地頭處自主導(dǎo)航轉(zhuǎn)彎的信息感知需求,基于機(jī)器視覺技術(shù),針對(duì)耕播時(shí)期的幾種典型地頭場(chǎng)景,研究非規(guī)整地頭邊界線的檢測(cè)方法,為農(nóng)機(jī)在地頭處實(shí)現(xiàn)自主轉(zhuǎn)向掉頭提供可靠的轉(zhuǎn)向基準(zhǔn)線,為實(shí)現(xiàn)農(nóng)機(jī)田間作業(yè)過程中全程自主導(dǎo)航行走提供技術(shù)支撐。

    1 農(nóng)田地頭圖像樣本采集

    本質(zhì)上,地頭屬于農(nóng)機(jī)田間作業(yè)過程中一種必然存在的障礙。每個(gè)自然農(nóng)田地塊均存在4 條以上邊界,通常由防護(hù)林、田埂、農(nóng)田道路、溝渠等包圍農(nóng)田地塊形成。其中,阻礙農(nóng)機(jī)正常前行作業(yè)或中斷農(nóng)機(jī)連續(xù)作業(yè)的田塊邊界(即壟行末端的田塊邊界)視為地頭邊界所在位置。本研究所檢測(cè)的地頭邊界線,主要產(chǎn)生于自然田塊邊界位置。農(nóng)機(jī)田間作業(yè)過程中,需在抵達(dá)地頭邊界線之前完成轉(zhuǎn)向掉頭,同時(shí)需在開始轉(zhuǎn)向掉頭之前完成地頭邊界線的探測(cè)。

    非結(jié)構(gòu)農(nóng)田環(huán)境中,地頭類型呈現(xiàn)復(fù)雜多樣性,本研究選取耕播時(shí)期的3 類典型地頭場(chǎng)景進(jìn)行地頭邊界線檢測(cè)算法的研究:①植被田埂/農(nóng)田,如灌叢、雜草、樹木等植被田埂。②泥土路/農(nóng)田。③水泥路/農(nóng)田。

    選用USB 接口的TXY_616_1080P 型高清彩色數(shù)碼攝像機(jī),輸出分辨率為800 像素×600 像素。于2018 年10 月在小湯山精準(zhǔn)農(nóng)業(yè)示范基地,針對(duì)3 類共9 種農(nóng)田場(chǎng)景,由田內(nèi)至田外以20 f/s 的速率連續(xù)采集獲取視頻圖像,圖1 為采集的3 類共9 種場(chǎng)景下地頭出現(xiàn)時(shí)的圖像樣本,圖2 為相應(yīng)地頭出現(xiàn)前的田內(nèi)場(chǎng)景圖像樣本。

    圖1 不同地頭場(chǎng)景圖像樣本Fig.1 Image samples of different headland scenes

    2 地頭邊界線視覺檢測(cè)算法

    圖2 圖1 中不同地頭場(chǎng)景下的田內(nèi)場(chǎng)景圖像樣本Fig.2 Image samples of field scenes corresponding to headland scenes in Fig.1

    面向非結(jié)構(gòu)復(fù)雜農(nóng)田環(huán)境,結(jié)合農(nóng)機(jī)田間動(dòng)態(tài)作業(yè)過程以及農(nóng)機(jī)在地頭處安全可靠地實(shí)施自主轉(zhuǎn)向掉頭的實(shí)際應(yīng)用需求,基于Matlab 平臺(tái)進(jìn)行非規(guī)整地頭邊界線視覺檢測(cè)算法的研究,完成以下3 項(xiàng)檢測(cè)任務(wù):①首要檢測(cè)任務(wù):判斷地頭出現(xiàn)與否。②根本檢測(cè)任務(wù):獲取非規(guī)整地頭邊界的主體延伸方位線。③關(guān)鍵檢測(cè)任務(wù):獲取地頭轉(zhuǎn)向基準(zhǔn)線。

    2.1 地頭出現(xiàn)與否判斷

    農(nóng)機(jī)田內(nèi)動(dòng)態(tài)作業(yè)過程中,需實(shí)時(shí)判斷其前方視野范圍內(nèi)是否出現(xiàn)地頭,只有在確定地頭出現(xiàn)的前提下,才需進(jìn)行地頭邊界線的檢測(cè)。尤其是當(dāng)田外場(chǎng)景較明顯地出現(xiàn)在圖像中時(shí),再檢測(cè)地頭邊界線,較為穩(wěn)妥。

    同一地塊中地頭還未出現(xiàn)時(shí)的田內(nèi)場(chǎng)景以及地頭出現(xiàn)后的地頭場(chǎng)景的行灰度平均值(即水平方向像素灰度平均值)分布如圖3 所示,每組圖中,左側(cè)為地頭圖像的行灰度平均值,右側(cè)為田內(nèi)圖像的行灰度平均值,對(duì)于田內(nèi)圖像,其分布圖整體變化較為平穩(wěn),而對(duì)于地頭圖像,對(duì)比地頭邊界兩側(cè)區(qū)域,其行灰度平均值通常存在較明顯的跳變。因此,基于該跳變特征判斷地頭是否出現(xiàn)。

    圖3 不同農(nóng)田場(chǎng)景下的地頭圖像和田內(nèi)圖像的行灰度平均值變化曲線Fig.3 Curves of row gray average values of different headland images and their corresponding field images

    具體檢測(cè)步驟如下:

    (1)彩色圖像灰度化。根據(jù)計(jì)算式

    式中 R(i,j)、G(i,j)、B(i,j)——原彩色圖像中(i,j)處像素點(diǎn)的R、G、B 分量值

    f(i,j)——綠色分量增強(qiáng)后的(i,j)處像素點(diǎn)的灰度

    對(duì)彩色圖像進(jìn)行灰度化。

    (2)計(jì)算行灰度平均值的分散程度。從上往下逐行掃描灰度圖,第j 行像素的灰度平均值為

    所有像素的灰度平均值為

    正向分布偏差為

    負(fù)向分布偏差為

    式中 xsize——圖像寬度 ysize——圖像高度

    (3)判斷地頭是否出現(xiàn)。若正向分布偏差D+>Dth(判斷閾值),或者負(fù)向分布偏差D - >Dth,則表明出現(xiàn)跳變特征,即認(rèn)為地頭出現(xiàn),繼續(xù)執(zhí)行后續(xù)檢測(cè)步驟。否則認(rèn)為地頭未出現(xiàn),結(jié)束當(dāng)前幀圖像的檢測(cè)。

    其中,地頭判斷閾值Dth主要取決于田內(nèi)場(chǎng)景圖像中上述行灰度平均值數(shù)據(jù)的正、負(fù)向分布偏差的分布范圍,尤其是分布范圍的上限值。

    另外,定義跳變方向標(biāo)志量f1,初始化為0。在判斷出現(xiàn)地頭的前提下,若D+≥D -,則表明地頭邊界處表現(xiàn)為正向跳變特征,令f1=1,否則認(rèn)為表現(xiàn)為負(fù)向跳變特征,令f1= -1。

    2.2 非規(guī)整地頭邊界的主體延伸方位線獲取

    非結(jié)構(gòu)復(fù)雜農(nóng)田環(huán)境中,地頭邊界通常呈非規(guī)則曲線性,但在短距離范圍內(nèi)可近似擬合為直線,擬合直線的方位即可視為當(dāng)前段地頭邊界的主體延伸方位。實(shí)際農(nóng)田環(huán)境中,地頭邊界的主體延伸方向不垂直壟行方向的情況較常見,在圖4 所示地頭場(chǎng)景圖像中,黑色箭頭、黑色曲線和橙色直線分別示意壟行方向、實(shí)際地頭邊界和曲形地頭邊界的主體延伸方位。此時(shí),采用求取跳變行的方法[18-20]來確定地頭邊界在圖像中的位置,僅能獲得一條紅色虛線所示的位于地頭邊界附近的水平行線,該行線并不能真實(shí)地展現(xiàn)地頭邊界的主體延伸方位。

    圖4 非規(guī)整地頭邊界線獲取示意圖Fig.4 Schematic of acquisition of irregular headland boundary line

    本研究擬將圖像沿水平方向平均分成8 個(gè)子處理區(qū)域,如圖4 所示,分別針對(duì)各子處理區(qū)域求取跳變位置處的特征點(diǎn)(圖4 中紅色小圓圈),最后再針對(duì)這些特征點(diǎn)進(jìn)行線性擬合,獲取實(shí)際非規(guī)整地頭邊界的主體延伸方位線(圖4 中橙色直線)。具體求取過程如下:

    (1)分區(qū)求取跳變位置處特征點(diǎn)。將灰度圖沿圖像寬度方向均分為8 個(gè)子處理區(qū)域,第m 個(gè)子處理區(qū)域沿寬度方向的像素區(qū)間范圍為

    式中 xsm、xem——第m 個(gè)子處理區(qū)域沿圖像寬度方向的起始、終止像素位置

    定義標(biāo)志變量f2,初始化為0。從第1 個(gè)子處理區(qū)域開始,從左往右依次針對(duì)各子處理區(qū)域進(jìn)行檢測(cè):

    步驟1:參照式(2)求取當(dāng)前子區(qū)域中每一行像素的灰度平均值,存儲(chǔ)于數(shù)組Fa中?;趕mooth 函數(shù)利用局部加權(quán)回歸方法對(duì)Fa中數(shù)據(jù)進(jìn)行平滑處理,平移窗口寬度設(shè)置為60 像素,平滑后的數(shù)據(jù)存儲(chǔ)于數(shù)組Fs中。以Fs數(shù)據(jù)作為縱坐標(biāo),各數(shù)據(jù)所對(duì)應(yīng)的行數(shù)作為橫坐標(biāo),繪制行灰度平均值平滑曲線,尋找該曲線上的所有波峰點(diǎn)和波谷點(diǎn),存儲(chǔ)其位置坐標(biāo)。

    步驟2:結(jié)合當(dāng)前地頭邊界跳變方向的不同(即f1的取值),確定位于跳變位置處的波峰點(diǎn)和波谷點(diǎn)方法如下:

    若f1=1,則針對(duì)上述所有波峰點(diǎn),按照其行數(shù)(即其橫坐標(biāo))從大到小的順序,對(duì)其進(jìn)行排列,并基于排序后的波峰點(diǎn)的縱坐標(biāo)數(shù)據(jù)來確定位于跳變位置處的波峰點(diǎn),從第2 個(gè)數(shù)據(jù)開始,按序計(jì)算當(dāng)前數(shù)據(jù)的離群度,計(jì)算式為

    式中 yi、yi-1——排序后的第i 和第i-1 個(gè)波峰點(diǎn)或波谷點(diǎn)的縱坐標(biāo),i∈[2,N],N表示波峰點(diǎn)或波谷點(diǎn)的總數(shù)

    直至尋找到首個(gè)滿足不等式R >Rth(判斷閾值)的數(shù)據(jù)為止,則該數(shù)據(jù)所對(duì)應(yīng)的波峰點(diǎn)即為跳變波峰點(diǎn),記其位置坐標(biāo)為(k,yk),同時(shí)本文將排于其之前的相鄰波峰點(diǎn)稱之為預(yù)跳變點(diǎn)(即其靠近田內(nèi)一側(cè)的相鄰波峰點(diǎn),亦即行數(shù)較大的相鄰波峰點(diǎn)),記其坐標(biāo)位置為(p,yp)。同時(shí),在步驟1 中的行灰度平均值平滑曲線上,尋找點(diǎn)(k,yk)右側(cè)最鄰近的波谷點(diǎn),稱之為跳前波谷點(diǎn),記其位置坐標(biāo)為(h,yh)。

    若f1= -1,則針對(duì)上述所有波谷點(diǎn),按照同樣的方式,尋找滿足條件的跳變波谷點(diǎn),記其位置坐標(biāo)為(h,yh),同樣將排于其之前的相鄰波谷點(diǎn)稱之為預(yù)跳變點(diǎn),記其位置坐標(biāo)為(p,yp)。同樣地,在步驟1中的行灰度平均值平滑曲線上,尋找點(diǎn)(h,yh)右側(cè)最鄰近的波峰點(diǎn),稱之為跳前波峰點(diǎn),記其位置坐標(biāo)為(k,yk)。

    其中,跳變特征判斷閾值Rth主要取決于田內(nèi)場(chǎng)景中的上述對(duì)應(yīng)曲線下的波峰和波谷點(diǎn)的按序離群程度,即其所對(duì)應(yīng)的|R|的分布范圍的上限值。

    若不存在滿足R >Rth的跳變波峰點(diǎn)或跳變波谷點(diǎn),則認(rèn)為當(dāng)前子處理區(qū)域中無跳變特征點(diǎn),結(jié)束當(dāng)前幀圖像的檢測(cè)。

    步驟3:基于跳變波峰點(diǎn)和跳前波谷點(diǎn)或跳變波谷點(diǎn)和跳前波峰點(diǎn)的位置坐標(biāo),確定當(dāng)前子處理區(qū)域中跳變特征點(diǎn)的像素坐標(biāo)位置。取跳變特征點(diǎn)所在行數(shù)為(k+h)/2,跳變特征點(diǎn)所在列為當(dāng)前子處理區(qū)域的中心列。若當(dāng)前為第m 個(gè)子處理區(qū)域,則當(dāng)前子處理區(qū)域中跳變特征點(diǎn)的像素坐標(biāo)為((xsm+xem)/2,(k+h)/2)。

    另外,定義標(biāo)志變量f2,初始化為0。當(dāng)6 個(gè)以上子處理區(qū)域均存在跳變特征點(diǎn)時(shí),令f2=1,同時(shí)記有效子處理區(qū)域個(gè)數(shù)(即存在跳變特征點(diǎn)的子處理區(qū)域個(gè)數(shù),亦即跳變特征點(diǎn)個(gè)數(shù))為M,否則令f2= -1。

    (2)線性擬合跳變特征點(diǎn),獲取主體延伸方位線。若f2=1,則針對(duì)上述求取的跳變特征點(diǎn),基于穩(wěn)健回歸法[21-22]進(jìn)行線性擬合,獲取擬合直線,從而獲得地頭邊界的主體延伸方位線。

    2.3 地頭轉(zhuǎn)向基準(zhǔn)線獲取

    從實(shí)際應(yīng)用的角度看,上述主體延伸方位線并不能直接用作農(nóng)機(jī)在當(dāng)前地頭處轉(zhuǎn)向掉頭的基準(zhǔn)線。如圖4 所示,當(dāng)農(nóng)機(jī)以橙色直線作為轉(zhuǎn)向基準(zhǔn)線(即要求農(nóng)機(jī)在抵達(dá)該轉(zhuǎn)向基準(zhǔn)線之前完成轉(zhuǎn)向掉頭)時(shí),A 點(diǎn)處凸向該基準(zhǔn)線以內(nèi)的田外場(chǎng)景(如田埂或道路等)很可能會(huì)對(duì)農(nóng)機(jī)當(dāng)前的正常轉(zhuǎn)向掉頭造成阻礙,導(dǎo)致轉(zhuǎn)向掉頭失敗,甚至產(chǎn)生碰撞等安全事故。

    為此,本研究擬將上述主體延伸方位線向下平行移動(dòng),使之經(jīng)過甚至跨越最遠(yuǎn)的邊界位置(以下簡(jiǎn)稱安全位置,圖4 中的A 點(diǎn)處),以此獲取在當(dāng)前地頭處安全轉(zhuǎn)向掉頭的基準(zhǔn)線,如圖4 中綠色直線所示。由此,該轉(zhuǎn)向基準(zhǔn)線既順沿實(shí)際地頭邊界的主體延伸方向,又將當(dāng)前所有凸向田內(nèi)的外部區(qū)域分割在外,可保證農(nóng)機(jī)在該轉(zhuǎn)向基準(zhǔn)線之內(nèi)實(shí)施自主轉(zhuǎn)向掉頭操作時(shí)安全可靠。操作步驟如下:

    (1)在前述M 個(gè)預(yù)跳變點(diǎn)中,尋找對(duì)應(yīng)行數(shù)最大(即橫坐標(biāo)值最大)的預(yù)跳變點(diǎn),并記該預(yù)跳變點(diǎn)所屬的子處理區(qū)域?yàn)殛P(guān)鍵子處理區(qū)域。

    (2)基于關(guān)鍵子處理區(qū)域的行灰度平均值平滑曲線,求取曲線上預(yù)跳變點(diǎn)左側(cè)數(shù)據(jù)的平均值和標(biāo)準(zhǔn)偏差D0。

    (3)在關(guān)鍵子處理區(qū)域內(nèi),以主體延伸方位線的初始位置作為起始位置,將其逐行往下平移,同時(shí)逐行計(jì)算線上像素的灰度平均值,并存儲(chǔ)于數(shù)組FL中。同樣基于smooth 函數(shù),設(shè)置平移窗口寬度為60像素,利用局部加權(quán)回歸方法對(duì)FL中數(shù)據(jù)進(jìn)行平滑處理,將平滑后的數(shù)據(jù)存儲(chǔ)于數(shù)組FLs中。以FLs中數(shù)據(jù)作為縱坐標(biāo),以對(duì)應(yīng)的平移行數(shù)作為橫坐標(biāo),繪制灰度平均值平滑曲線。

    (4)再次結(jié)合當(dāng)前地頭邊界跳變方向的不同,分別采取以下方法將主體延伸方位線平移至安全位置處:

    當(dāng)f1= -1 時(shí),在上述平滑曲線上,從左往右按序?qū)ふ沂讉€(gè)滿足縱坐標(biāo)值大于-D0的點(diǎn),該點(diǎn)所對(duì)應(yīng)的平移行數(shù)即可認(rèn)為是原主體延伸方位線抵達(dá)安全位置處所需向下平移的最小行數(shù),本文將滿足以上平移條件的平移方案稱之為平移途徑1。

    當(dāng)f1=1 時(shí),根據(jù)D -與Dth-關(guān)系,采取平移方法如下:①若D -≤Dth-,則首先在上述平滑曲線上,尋找到第1 個(gè)波谷點(diǎn),然后在該波谷點(diǎn)之前的平滑曲線上尋找到縱坐標(biāo)值小于且最接近+ D0的點(diǎn),該點(diǎn)所對(duì)應(yīng)的平移行數(shù),即可認(rèn)為是原主體延伸方位線抵達(dá)安全位置處所需向下平移的最小行數(shù),本文將滿足以上平移條件的平移方案稱之為平移途徑2。②若D - >Dth-,則在上述平滑曲線上,首先尋找到第1 個(gè)波谷點(diǎn),然后以第1 個(gè)波谷點(diǎn)為起點(diǎn),向后按序?qū)ふ业绞讉€(gè)滿足縱坐標(biāo)值大于的點(diǎn),同樣地,該點(diǎn)所對(duì)應(yīng)的平移行數(shù),即可認(rèn)為是原主體延伸方位線抵達(dá)安全位置處所需向下平移的最小行數(shù),本文將滿足以上平移條件的平移方案稱之為平移途徑3。

    依據(jù)上述3 種平移途徑所獲取的最小平移行數(shù),將主體延伸方位線平移到位后,即可獲得當(dāng)前地頭處轉(zhuǎn)向掉頭的基準(zhǔn)線,即為本研究最終所要求取的用于農(nóng)機(jī)自主導(dǎo)航的地頭邊界線。

    2.4 總體檢測(cè)流程

    總體檢測(cè)流程如圖5 所示。

    3 試驗(yàn)與結(jié)果分析

    3.1 閾值確定

    從前述3 類共9 種農(nóng)田場(chǎng)景視頻圖像中,截取各場(chǎng)景下的田內(nèi)場(chǎng)景視頻段,從第1 幀開始,以每20 幀選取一幀的方式,從截取的每段田內(nèi)場(chǎng)景視頻中,各選取50 幀圖像。

    圖5 總體檢測(cè)流程圖Fig.5 Flow chart of detection method

    (1)地頭判斷閾值

    分別針對(duì)上述9 種田內(nèi)場(chǎng)景下的各50 幀圖像樣本,計(jì)算每幀圖像的行灰度平均值數(shù)據(jù)間的正向分布偏差DIn+和負(fù)向分布偏差DIn-,并統(tǒng)計(jì)這50 組DIn+和DIn-的分布區(qū)間,設(shè)兩分布區(qū)間上限值中的最大值為Dmax,最后根據(jù)Dmax確定當(dāng)前場(chǎng)景下的地頭判斷閾值Dth。事實(shí)上,本研究設(shè)置Dth=1.3Dmax,通過預(yù)留一定的余量以便地頭較明顯地出現(xiàn)在圖像上方后再確認(rèn)地頭出現(xiàn)并進(jìn)行后續(xù)檢測(cè),以此提高地頭邊界線的可見性,便于后續(xù)對(duì)其進(jìn)行準(zhǔn)確檢測(cè)。由此,當(dāng)圖像中僅冒出少部分地頭時(shí)本研究視之為地頭未出現(xiàn)的情況。9 種田內(nèi)場(chǎng)景下的DIn+和DIn-分布區(qū)間、Dmax以及對(duì)應(yīng)的Dth見表1。

    表1 地頭判斷閾值和跳變特征判斷閾值Tab.1 Determination of headland judgement threshold and jumping feature judgement threshold

    農(nóng)機(jī)田間作業(yè)過程中,在作業(yè)前期,通常不會(huì)遇見地頭,所采集的圖像一般為田內(nèi)場(chǎng)景圖像,由此在實(shí)際應(yīng)用中,可基于前期連續(xù)采集的圖像樣本集確定地頭判斷閾值。

    (2)跳變特征判斷閾值

    針對(duì)上述9 種田內(nèi)場(chǎng)景下的各50 幀圖像樣本,計(jì)算每幀圖像的行灰度平均值平滑曲線上波峰點(diǎn)間和波谷點(diǎn)間的按序離群程度|R|的取值范圍,其中f1取1。并統(tǒng)計(jì)這100 組取值范圍的分布區(qū)間,設(shè)分布區(qū)間上限值為Rmax,根據(jù)Rmax確定當(dāng)前場(chǎng)景下的跳變特征判斷閾值Rth。本研究設(shè)置Rth=1.3Rmax,同樣通過預(yù)留一定的余量以避免凸起的農(nóng)田道路其側(cè)面非均質(zhì)的小段斷層陰影對(duì)田內(nèi)外灰度跳變特征的判斷造成負(fù)向干擾,如圖1c(場(chǎng)景9)所示的水泥地頭場(chǎng)景,從其行灰度平均值分布(圖3d)可知,水泥道路相對(duì)田內(nèi)存在明顯的灰度正向跳變,而其側(cè)面小段的斷層陰影相對(duì)田內(nèi)存在局部灰度負(fù)向跳變。由此,針對(duì)植被田埂類地頭場(chǎng)景和非凸起的道路類地頭場(chǎng)景,通過獲取田埂或道路與田內(nèi)分界處的灰度跳變特征來確定地頭邊界的主體延伸方位,而針對(duì)凸起的道路類地頭場(chǎng)景,基于道路與其側(cè)面斷層區(qū)分界處的灰度跳變特征來確定地頭邊界的主體延伸方位。9 種田內(nèi)場(chǎng)景下的Rmax以及對(duì)應(yīng)的Rth見表1。

    在實(shí)際應(yīng)用中,同樣地,可基于作業(yè)前期連續(xù)采集的田內(nèi)圖像樣本集來確定跳變特征判斷閾值。

    3.2 地頭出現(xiàn)與否判斷結(jié)果

    針對(duì)圖1 中的3 類共9 種地頭場(chǎng)景圖像樣本,計(jì)算獲得其正向分布偏差D+和負(fù)向分布偏差D -依次為:17.63、43.83;20.07、46.54;50.95、15.04;32.11、24.95;36.65、17.23;20.70、12.12;25.00、19.43;54.88、25.04;20.69、18.49。結(jié)合表1 可知,圖1a 中場(chǎng)景1、2 的D -,圖1a 中場(chǎng)景3 的D+,圖1b中場(chǎng)景4、5 的D+,以及圖1c 中場(chǎng)景7、8 的D+,均小于對(duì)應(yīng)場(chǎng)景下的閾值Dth,由此判斷上述圖像均為地頭出現(xiàn)時(shí)的場(chǎng)景,符合實(shí)際情況,均判斷正確。對(duì)于圖1b 中場(chǎng)景6 和圖1c 中場(chǎng)景9,由于圖像中僅有少部分地頭冒出,導(dǎo)致其D+雖然大于對(duì)應(yīng)場(chǎng)景下的Dmax,但卻均小于對(duì)應(yīng)場(chǎng)景下設(shè)置的閾值Dth,由此均判斷為地頭未出現(xiàn)時(shí)的場(chǎng)景,此時(shí)停止當(dāng)前幀圖像的檢測(cè)。而當(dāng)這2 種場(chǎng)景下對(duì)應(yīng)的地頭較明顯地出現(xiàn)在圖像上方后,如圖6 所示,分別為場(chǎng)景6、9各自所對(duì)應(yīng)場(chǎng)景下地頭出現(xiàn)較明顯后的圖像樣本,計(jì)算獲得其D+和D -依次為:38.32、16.09;42.33、24.10。此時(shí)圖6a 和圖6b 的D+均小于對(duì)應(yīng)場(chǎng)景下的閾值Dth,由此均判斷為地頭出現(xiàn)時(shí)的場(chǎng)景,此時(shí)執(zhí)行后續(xù)檢測(cè)步驟。由此可見,本研究的地頭判斷閾值設(shè)置合理,判斷結(jié)果符合預(yù)期效果。

    圖6 出現(xiàn)較明顯地頭時(shí)的圖像樣本Fig.6 Image samples of headland scenes after obvious headland appeared at top of image

    同樣地,針對(duì)圖2 中的3 類共9 種田內(nèi)場(chǎng)景圖像樣本,計(jì)算獲得其正向分布偏差D+和負(fù)向分布偏差D - 依次為:8.58、10.33;10.22、13.52;6.71、6.04;18.60、12.59;8.07、10.57;10.14、10.02;11.00、15.23;11.47、13.35;9.42、9.65。結(jié)合表1可知,以上各值均小于對(duì)應(yīng)場(chǎng)景下的閾值Dth,由此判斷上述圖像均為地頭未出現(xiàn)時(shí)的場(chǎng)景,符合實(shí)際情況,均判斷正確。

    3.3 主體延伸方位線檢測(cè)結(jié)果

    (1)跳變特征點(diǎn)求取結(jié)果

    以圖6b 為例,其灰度圖如圖7a 所示,其跳變方向標(biāo)志量f1=1,將其灰度圖沿水平方向平均分成8個(gè)子處理區(qū)域,圖7a 中紅色方框區(qū)域表示首個(gè)子處理區(qū)域。該子處理區(qū)域所對(duì)應(yīng)的行灰度平均值見圖7b 中藍(lán)色曲線。該平滑曲線上共計(jì)有14 個(gè)波峰點(diǎn)和14 個(gè)波谷點(diǎn)。按照行數(shù)從大到小的順序,14 個(gè)波峰點(diǎn)的位置坐標(biāo)依次為:(546,148)、(530,147)、(483,130)、(469,133)、(413,134)、(403,130)、(336,145)、(274,145)、(235,138)、(197,147)、(182,148)、(135,118)、(47,213)、(29,215)。從以上第2 個(gè)波峰點(diǎn)開始,依次計(jì)算并獲得其縱坐標(biāo)值的按序離群度R 分別為:-1.29%、-11.15%、2.34%、0.65%、- 2.68%、10.54%、- 0.21%、-5.01%、6.44%、0.80%、- 21.36%、68.47%、1.46%。

    結(jié)合表1 可知,第13 個(gè)波峰點(diǎn)的R 值大于當(dāng)前的跳變特征判斷閾值19.44%,由此確定波峰點(diǎn)(47,213)為跳變波峰點(diǎn),波峰點(diǎn)(135,118)為預(yù)跳變波峰點(diǎn)。另外,結(jié)合波谷點(diǎn)的位置坐標(biāo),確定跳變前波谷點(diǎn)的位置坐標(biāo)為(104,92)。由此,基于跳變波峰點(diǎn)和跳變前波谷點(diǎn)位置坐標(biāo)確定首個(gè)子處理區(qū)域中跳變特征點(diǎn)的像素坐標(biāo)為(50,76)。同理依次獲得剩余7 個(gè)子處理區(qū)域中跳變特征點(diǎn)的像素坐標(biāo)為(150,67)、(250,65)、(350,70)、(450,66)、(550,75)、(650,67)、(750,80)。

    如圖7b 所示,紅色和綠色小圓圈中心分別表示所獲取的波谷點(diǎn)、波峰點(diǎn)位置,藍(lán)色、青綠色和黑色小圓圈中心分別表示所獲取的跳變波峰點(diǎn)、跳變前波谷點(diǎn)和預(yù)跳變點(diǎn)。8 個(gè)跳變特征點(diǎn)在原圖中的位置如圖7c 中紅色圓圈中心所示,由圖可知,檢測(cè)獲取的跳變特征點(diǎn)定位較準(zhǔn)確。

    圖7 圖6b 地頭場(chǎng)景圖下的跳變特征點(diǎn)檢測(cè)過程示意圖Fig.7 Schematics of detection process of jump feature points for Fig.6b

    (2)主體延伸方位線擬合結(jié)果

    上述8 個(gè)跳變特征點(diǎn)線性擬合結(jié)果如圖8i 中黃色直線所示,該直線即為圖6b 中地頭場(chǎng)景的主體延伸方位線。另外,圖1a、1b 中場(chǎng)景4、5 和圖6a、1c中場(chǎng)景7、8 的地頭主體延伸方位線檢測(cè)結(jié)果如圖8中黃色直線所示。其中,圖1c 中場(chǎng)景8 和圖6b 為凸起的道路類地頭場(chǎng)景,由檢測(cè)結(jié)果可知,所獲取的主體延伸方位線較好地?cái)M合了道路與其側(cè)面斷層區(qū)之間的分界位置,而對(duì)于剩余的地頭場(chǎng)景圖像,所獲取的主體延伸方位線也較好地?cái)M合了田內(nèi)與田外區(qū)域之間的分界位置,檢測(cè)結(jié)果均較準(zhǔn)確。

    3.4 轉(zhuǎn)向基準(zhǔn)線檢測(cè)結(jié)果

    同樣地,以圖6b 中地頭場(chǎng)景圖像為例,8 個(gè)預(yù)跳變點(diǎn)中,第2 和第5 個(gè)子處理區(qū)域的預(yù)跳變點(diǎn)所對(duì)應(yīng)的行數(shù)最大,兩者任選其一作為關(guān)鍵子區(qū)域。選取第5 個(gè)子處理區(qū)域作為關(guān)鍵子區(qū)域。在該關(guān)鍵子區(qū)域的行灰度平均值平滑曲線上,計(jì)算獲得預(yù)跳變點(diǎn)左側(cè)數(shù)據(jù)的平均值為137.34。在關(guān)鍵子區(qū)域內(nèi),將主體延伸方位線逐行向下平移,并逐行計(jì)算平移后的線上像素灰度平均值,如圖9a 中黃色曲線所示,獲得其平滑曲線如圖中藍(lán)色曲線所示。另外,圖中紅色小圓圈中心表示平滑曲線上的波谷點(diǎn)。鑒于圖6b 中地頭場(chǎng)景圖像下f1=1,且D - =24.10 >Dth= 23.39,由此依據(jù)平移途徑3 進(jìn)行平移。由圖9a 可知,藍(lán)色小圓圈中心所示位置處,其縱坐標(biāo)值為139.71,為平滑曲線上第1 個(gè)波谷點(diǎn)之后首個(gè)大于的值,所對(duì)應(yīng)的平移行數(shù)為112 行,故原主體延伸方位線抵達(dá)安全位置處所需向下平移的最小行數(shù)為112 行。

    另外,對(duì)于圖1a 中地頭場(chǎng)景1 和圖1b 中地頭場(chǎng)景5,其f1、D -、Dth、關(guān)鍵子區(qū)域、D0依次為-1、43.83、18.31、第5 個(gè)子處理區(qū)域、150.82、16.76,1、17.23、21.55、第5 個(gè)子處理區(qū)域、137.11、18.26。同樣地,在各自的關(guān)鍵子區(qū)域內(nèi),計(jì)算并獲得各自的主體延伸方位線逐行向下平移過程中其線上灰度平均值數(shù)據(jù)的平滑曲線,分別如圖9b、9c 所示。結(jié)合平滑曲線,確定上述3 類地頭場(chǎng)景圖像下所需采取的平移途徑依次為平移途徑1 和平移途徑2。由圖9 可知,以上兩類地頭場(chǎng)景圖像下,原主體延伸方位線抵達(dá)安全位置處所需向下平移的最小行數(shù)分別為42、19 行。

    圖9 主體延伸方位線平移行數(shù)確定示意圖Fig.9 Determination of row number of downward translation of main-body extended azimuth line

    按照上述最小平移行數(shù)將圖6b 場(chǎng)景9、圖1a中場(chǎng)景1、圖1b 中場(chǎng)景5 的主體延伸方位線分別平移至安全位置處,由此獲得相應(yīng)地頭處轉(zhuǎn)向掉頭的基準(zhǔn)線,分別如圖8i、8a、8e 中綠色直線所示。

    此外,剩余地頭場(chǎng)景圖像的轉(zhuǎn)向基準(zhǔn)線獲取結(jié)果如圖8 中綠色直線所示。由圖8 可知,檢測(cè)獲取的轉(zhuǎn)向基準(zhǔn)線,可保證農(nóng)機(jī)在該轉(zhuǎn)向基準(zhǔn)線以內(nèi)完成掉頭轉(zhuǎn)向操作時(shí)安全可靠,即為本研究所要求取的用于農(nóng)機(jī)自主導(dǎo)航的地頭邊界線。

    3.5 總體檢測(cè)結(jié)果

    從前述采集的3 類共9 種農(nóng)田場(chǎng)景視頻中,針對(duì)每種場(chǎng)景,從中截取地頭出現(xiàn)前的田內(nèi)連續(xù)圖像幀共計(jì)2 000 幀,并截取地頭出現(xiàn)較明顯后的連續(xù)圖像幀共計(jì)100 幀。針對(duì)每種場(chǎng)景下的2 100 幀圖像,逐幀檢測(cè)地頭是否出現(xiàn),并判斷檢測(cè)結(jié)果的準(zhǔn)確性,同時(shí)針對(duì)每種場(chǎng)景下的上述100 幀圖像,逐幀進(jìn)行地頭邊界線檢測(cè),僅當(dāng)能檢測(cè)出地頭邊界線且所檢測(cè)出的地頭邊界線可保證農(nóng)機(jī)在其內(nèi)實(shí)施自主轉(zhuǎn)向掉頭操作安全可靠時(shí),才視當(dāng)前幀地頭邊界線檢測(cè)正確,同時(shí)統(tǒng)計(jì)單幀圖像檢測(cè)處理的最長(zhǎng)時(shí)間,試驗(yàn)結(jié)果如表2 所示。

    表2 檢測(cè)結(jié)果Tab.2 Test results

    由表2 可知,地頭出現(xiàn)與否判斷準(zhǔn)確率不低于96%,地頭邊界線檢測(cè)準(zhǔn)確率不低于92%,基于Matlab 平臺(tái)單幀圖像處理時(shí)間不高于0.52 s。試驗(yàn)結(jié)果表明,本算法可準(zhǔn)確地判斷是否抵達(dá)地頭,并快速檢測(cè)出用于農(nóng)機(jī)自主導(dǎo)航的地頭邊界線,所檢測(cè)出的地頭邊界線可保證農(nóng)機(jī)在該轉(zhuǎn)向基準(zhǔn)線以內(nèi)完成轉(zhuǎn)向掉頭操作時(shí)安全可靠。

    4 結(jié)論

    (1)基于農(nóng)田內(nèi)外像素灰度值的跳變特征來判斷地頭是否出現(xiàn)。建立并計(jì)算正向分布偏差和負(fù)向分布偏差兩個(gè)度量,分別用以描述圖像行灰度平均值數(shù)據(jù)的正向和負(fù)向分散程度,當(dāng)兩度量之一大于判斷閾值(即數(shù)據(jù)分布較為分散)時(shí),即可認(rèn)為出現(xiàn)跳變特征,判斷地頭出現(xiàn)。該判斷方法主要適用于農(nóng)田內(nèi)外具有灰度跳變特性且田內(nèi)勻質(zhì)分布的農(nóng)田場(chǎng)景,檢測(cè)準(zhǔn)確率不低于96%。當(dāng)田內(nèi)存在大塊異物或雜質(zhì),或存在其他導(dǎo)致田內(nèi)局部區(qū)域灰度差異較大的影響因素(如陰影、灌溉不均勻等)時(shí),則會(huì)造成判斷準(zhǔn)確性下降,同時(shí)也將影響后續(xù)檢測(cè)。

    (2)將圖像沿水平方向平均分成8 個(gè)子處理區(qū)域,針對(duì)各子處理區(qū)域求取行灰度平均值分布圖,基于局部加權(quán)回歸法對(duì)該分布圖進(jìn)行平滑處理。針對(duì)存在正向(或負(fù)向)跳變特征的場(chǎng)景,基于平滑曲線上的波峰點(diǎn)(或波谷點(diǎn)),引入按序離群度,用以描述按行數(shù)大小倒序排列后的波峰點(diǎn)或波谷點(diǎn)其所對(duì)應(yīng)的行灰度平均值的按序離群程度,將首個(gè)離群度較大的波峰點(diǎn)或波谷點(diǎn)視為跳變波峰點(diǎn)或跳變波谷點(diǎn)。基于跳變波峰點(diǎn)和跳變前波谷點(diǎn)或跳變波谷點(diǎn)和跳變前波峰點(diǎn)位置坐標(biāo),確定跳變特征點(diǎn)的像素坐標(biāo)位置。最后,基于穩(wěn)健回歸法對(duì)跳變特征點(diǎn)進(jìn)行線性擬合,獲取實(shí)際非規(guī)整地頭邊界的主體延伸方位線。試驗(yàn)結(jié)果表明,所獲取的主體延伸方位線較好地?cái)M合了田內(nèi)與田外區(qū)域之間的分界位置,檢測(cè)結(jié)果較為準(zhǔn)確。

    (3)根據(jù)預(yù)跳變點(diǎn)所對(duì)應(yīng)的行數(shù)確定田外區(qū)域凸向田內(nèi)程度最大的關(guān)鍵子區(qū)域,基于該子區(qū)域的行灰度平均值擬合曲線,求取預(yù)跳變點(diǎn)左側(cè)數(shù)據(jù)的平均值和標(biāo)準(zhǔn)差,確定田內(nèi)像素的灰度分布特征。將主體延伸方位線向下平行移動(dòng),當(dāng)其線上像素的灰度平均值接近于田內(nèi)像素的灰度分布特征時(shí),認(rèn)為抵達(dá)安全位置處,由此獲取農(nóng)機(jī)在當(dāng)前地頭處安全轉(zhuǎn)向掉頭的邊界線。試驗(yàn)結(jié)果表明,3 類共9 種農(nóng)田場(chǎng)景下,地頭邊界線檢測(cè)準(zhǔn)確率不低于92%。本研究可為農(nóng)機(jī)在地頭處實(shí)施自主導(dǎo)航轉(zhuǎn)彎提供較為準(zhǔn)確、可靠的地頭信息感知技術(shù)支持。

    猜你喜歡
    邊界線波谷波峰
    弟弟尿床了
    板厚與波高對(duì)波紋鋼管涵受力性能影響分析
    梅緣稻
    作用于直立堤墻與樁柱的波峰高度分析計(jì)算
    “邊界線”風(fēng)波
    “邊界線”風(fēng)波
    神奇的邊界線:一不留神就出國(guó)
    智慧少年(2017年8期)2018-01-10 21:39:12
    兒童標(biāo)準(zhǔn)12導(dǎo)聯(lián)T波峰末間期的分析
    基于音節(jié)時(shí)間長(zhǎng)度高斯擬合的漢語音節(jié)切分方法
    Dynamic Loads and Wake Prediction for Large Wind Turbines Based on Free Wake Method
    嫁个100分男人电影在线观看| 亚洲精品一卡2卡三卡4卡5卡| 最好的美女福利视频网| 露出奶头的视频| 在线观看免费午夜福利视频| 夜夜躁狠狠躁天天躁| 国产黄片美女视频| 国产伦精品一区二区三区视频9 | 免费在线观看日本一区| 成人午夜高清在线视频| 在线观看午夜福利视频| 热99re8久久精品国产| 国产视频内射| 国产69精品久久久久777片 | 一级黄色大片毛片| 亚洲 欧美 日韩 在线 免费| 亚洲欧美日韩高清在线视频| 色综合欧美亚洲国产小说| av片东京热男人的天堂| 国产成人啪精品午夜网站| 国产免费男女视频| av国产免费在线观看| 我要搜黄色片| x7x7x7水蜜桃| 国产精品电影一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | 噜噜噜噜噜久久久久久91| 麻豆国产av国片精品| 国产 一区 欧美 日韩| 亚洲色图av天堂| 啦啦啦免费观看视频1| 国产黄色小视频在线观看| 97超级碰碰碰精品色视频在线观看| 日日摸夜夜添夜夜添小说| 精品国产超薄肉色丝袜足j| 国产精品免费一区二区三区在线| 日韩高清综合在线| 精品一区二区三区av网在线观看| 在线观看免费视频日本深夜| 91麻豆精品激情在线观看国产| 亚洲黑人精品在线| 亚洲男人的天堂狠狠| 在线十欧美十亚洲十日本专区| 此物有八面人人有两片| 90打野战视频偷拍视频| 国产69精品久久久久777片 | 亚洲欧美日韩高清专用| 男女那种视频在线观看| 免费av不卡在线播放| 午夜两性在线视频| 亚洲av成人av| 中文字幕精品亚洲无线码一区| 日日夜夜操网爽| 999精品在线视频| 国产成年人精品一区二区| 桃红色精品国产亚洲av| 国产探花在线观看一区二区| 丰满人妻一区二区三区视频av | 国产爱豆传媒在线观看| 久久久水蜜桃国产精品网| 高清毛片免费观看视频网站| 99热精品在线国产| 欧美又色又爽又黄视频| 18禁美女被吸乳视频| 中出人妻视频一区二区| 国产精品98久久久久久宅男小说| 亚洲成人精品中文字幕电影| 久久中文字幕一级| 精品日产1卡2卡| 国产精品乱码一区二三区的特点| 国产伦精品一区二区三区视频9 | 久久久国产成人精品二区| 好看av亚洲va欧美ⅴa在| 午夜福利高清视频| 18禁黄网站禁片免费观看直播| 国产高清激情床上av| 亚洲片人在线观看| 女生性感内裤真人,穿戴方法视频| 搡老妇女老女人老熟妇| 麻豆av在线久日| 最新在线观看一区二区三区| 免费在线观看日本一区| 亚洲 欧美 日韩 在线 免费| 18禁观看日本| 日本与韩国留学比较| 久久精品夜夜夜夜夜久久蜜豆| 在线免费观看的www视频| 午夜福利免费观看在线| 国产精品 国内视频| 99精品久久久久人妻精品| 国产亚洲欧美98| 亚洲天堂国产精品一区在线| 丁香六月欧美| 亚洲专区字幕在线| 一个人看的www免费观看视频| 免费在线观看影片大全网站| 无遮挡黄片免费观看| 亚洲国产看品久久| 久久久水蜜桃国产精品网| 亚洲aⅴ乱码一区二区在线播放| 99在线视频只有这里精品首页| 国内精品久久久久精免费| 中文字幕人成人乱码亚洲影| 久久中文看片网| 香蕉国产在线看| 啦啦啦韩国在线观看视频| 91字幕亚洲| 噜噜噜噜噜久久久久久91| 性欧美人与动物交配| 国产毛片a区久久久久| 久久精品夜夜夜夜夜久久蜜豆| 老鸭窝网址在线观看| 午夜久久久久精精品| 国产午夜福利久久久久久| 久久久水蜜桃国产精品网| 免费在线观看亚洲国产| 国产黄片美女视频| 亚洲精品国产精品久久久不卡| 操出白浆在线播放| 曰老女人黄片| 亚洲av成人av| 欧美一级毛片孕妇| 动漫黄色视频在线观看| 极品教师在线免费播放| 亚洲天堂国产精品一区在线| av天堂在线播放| 男女之事视频高清在线观看| 精品久久久久久久久久久久久| 国产午夜福利久久久久久| 欧美日韩中文字幕国产精品一区二区三区| 特大巨黑吊av在线直播| 999久久久精品免费观看国产| 99在线人妻在线中文字幕| av片东京热男人的天堂| 在线a可以看的网站| 99热只有精品国产| 成年免费大片在线观看| 在线观看免费视频日本深夜| 亚洲中文av在线| 首页视频小说图片口味搜索| 亚洲 欧美 日韩 在线 免费| 午夜亚洲福利在线播放| 99热这里只有是精品50| 无遮挡黄片免费观看| 男女视频在线观看网站免费| 国产1区2区3区精品| 日本与韩国留学比较| 99久久无色码亚洲精品果冻| 成人欧美大片| 午夜a级毛片| 久久久色成人| 欧美精品啪啪一区二区三区| 色吧在线观看| 亚洲人成伊人成综合网2020| 亚洲成人免费电影在线观看| 叶爱在线成人免费视频播放| 观看免费一级毛片| 无遮挡黄片免费观看| 中国美女看黄片| 99精品在免费线老司机午夜| 在线视频色国产色| av女优亚洲男人天堂 | 亚洲欧美日韩卡通动漫| 国产伦精品一区二区三区四那| 亚洲美女黄片视频| 两性夫妻黄色片| 日韩欧美在线乱码| 99久久精品热视频| 最新在线观看一区二区三区| 色哟哟哟哟哟哟| 欧美性猛交╳xxx乱大交人| 法律面前人人平等表现在哪些方面| 亚洲一区二区三区色噜噜| 黄色视频,在线免费观看| 亚洲男人的天堂狠狠| 人妻夜夜爽99麻豆av| 听说在线观看完整版免费高清| 精品久久久久久成人av| 欧美+亚洲+日韩+国产| 国产成人精品无人区| 日本一本二区三区精品| 国产又色又爽无遮挡免费看| bbb黄色大片| 成年版毛片免费区| or卡值多少钱| 亚洲欧美日韩高清在线视频| 亚洲美女视频黄频| 亚洲国产色片| 白带黄色成豆腐渣| 国产精品自产拍在线观看55亚洲| 精品日产1卡2卡| 黄色 视频免费看| 我要搜黄色片| 看片在线看免费视频| 国产一区在线观看成人免费| 99久国产av精品| 国产精品乱码一区二三区的特点| 在线观看一区二区三区| 免费人成视频x8x8入口观看| 国产精品永久免费网站| 小说图片视频综合网站| 岛国视频午夜一区免费看| 成人av在线播放网站| 9191精品国产免费久久| 亚洲国产精品sss在线观看| 欧美zozozo另类| 日本在线视频免费播放| 亚洲成人免费电影在线观看| АⅤ资源中文在线天堂| 99在线视频只有这里精品首页| 一级毛片女人18水好多| 精品人妻1区二区| 成人无遮挡网站| 99久久99久久久精品蜜桃| 久久久成人免费电影| 狂野欧美白嫩少妇大欣赏| 手机成人av网站| 国产三级黄色录像| 国产精品久久电影中文字幕| 男女午夜视频在线观看| 午夜福利视频1000在线观看| 又黄又粗又硬又大视频| 女生性感内裤真人,穿戴方法视频| 亚洲国产中文字幕在线视频| 女人高潮潮喷娇喘18禁视频| 又大又爽又粗| 国产探花在线观看一区二区| 精品国产亚洲在线| 国产精品久久久人人做人人爽| 免费观看人在逋| 欧美日韩亚洲国产一区二区在线观看| 亚洲午夜精品一区,二区,三区| 国产蜜桃级精品一区二区三区| 日日干狠狠操夜夜爽| 成年女人毛片免费观看观看9| 亚洲国产精品合色在线| 欧美在线黄色| 精品久久久久久久末码| 99久久精品国产亚洲精品| 精品一区二区三区av网在线观看| 999精品在线视频| 亚洲七黄色美女视频| 免费观看精品视频网站| 一进一出好大好爽视频| 亚洲一区高清亚洲精品| 波多野结衣高清无吗| 嫁个100分男人电影在线观看| 免费在线观看视频国产中文字幕亚洲| 视频区欧美日本亚洲| 亚洲第一电影网av| av国产免费在线观看| 熟女人妻精品中文字幕| 亚洲精品在线美女| 亚洲精品久久国产高清桃花| 手机成人av网站| 国产精品日韩av在线免费观看| 在线播放国产精品三级| 啦啦啦免费观看视频1| bbb黄色大片| 国产视频一区二区在线看| 欧美日韩乱码在线| 波多野结衣巨乳人妻| 久久婷婷人人爽人人干人人爱| 精品无人区乱码1区二区| 成年女人看的毛片在线观看| 成人无遮挡网站| 久久精品亚洲精品国产色婷小说| 久久久久久国产a免费观看| 国产人伦9x9x在线观看| 在线观看一区二区三区| 夜夜夜夜夜久久久久| 亚洲中文字幕一区二区三区有码在线看 | 国产综合懂色| 精品国产三级普通话版| e午夜精品久久久久久久| 亚洲精品在线美女| 高清毛片免费观看视频网站| 三级男女做爰猛烈吃奶摸视频| 欧美性猛交╳xxx乱大交人| 高潮久久久久久久久久久不卡| 男人舔女人下体高潮全视频| 久久精品国产99精品国产亚洲性色| 国产精品综合久久久久久久免费| 禁无遮挡网站| 中文亚洲av片在线观看爽| 国产黄色小视频在线观看| 一级毛片高清免费大全| 老汉色av国产亚洲站长工具| 国产精品永久免费网站| 国产又色又爽无遮挡免费看| 91麻豆av在线| 午夜影院日韩av| 麻豆国产av国片精品| 在线a可以看的网站| 亚洲av第一区精品v没综合| 久久久久久久久中文| 中文字幕人妻丝袜一区二区| xxx96com| 97超视频在线观看视频| 级片在线观看| 九九久久精品国产亚洲av麻豆 | 国内少妇人妻偷人精品xxx网站 | 成人三级黄色视频| 亚洲成a人片在线一区二区| 一进一出抽搐gif免费好疼| 久久久久九九精品影院| 欧美日韩综合久久久久久 | 国产熟女xx| 哪里可以看免费的av片| 国产综合懂色| 国产精品,欧美在线| 美女扒开内裤让男人捅视频| 2021天堂中文幕一二区在线观| 欧美日韩一级在线毛片| 两个人的视频大全免费| 成年人黄色毛片网站| 在线a可以看的网站| 亚洲av熟女| 91九色精品人成在线观看| 国产成人福利小说| 窝窝影院91人妻| 欧美黑人欧美精品刺激| 亚洲在线自拍视频| 在线免费观看的www视频| 国内久久婷婷六月综合欲色啪| 在线视频色国产色| 亚洲第一欧美日韩一区二区三区| 97人妻精品一区二区三区麻豆| 美女 人体艺术 gogo| 身体一侧抽搐| 在线观看日韩欧美| 久久久久久九九精品二区国产| 男插女下体视频免费在线播放| 国产黄片美女视频| 亚洲无线观看免费| 99国产综合亚洲精品| 久久久久久久久中文| 国产三级黄色录像| 国产精品电影一区二区三区| 老汉色av国产亚洲站长工具| 丁香六月欧美| 一个人免费在线观看的高清视频| 欧美乱码精品一区二区三区| 国产精品精品国产色婷婷| 亚洲av中文字字幕乱码综合| 亚洲av日韩精品久久久久久密| 成人欧美大片| 国产精品日韩av在线免费观看| 日日夜夜操网爽| 美女黄网站色视频| 淫妇啪啪啪对白视频| 亚洲激情在线av| 久久伊人香网站| 国内揄拍国产精品人妻在线| 欧洲精品卡2卡3卡4卡5卡区| 在线永久观看黄色视频| 中文字幕最新亚洲高清| 99精品欧美一区二区三区四区| 国内久久婷婷六月综合欲色啪| 久久精品aⅴ一区二区三区四区| 久久久久国产一级毛片高清牌| 国产一区在线观看成人免费| 精品久久久久久成人av| 在线永久观看黄色视频| 日韩欧美国产在线观看| 男女床上黄色一级片免费看| 国产探花在线观看一区二区| 视频区欧美日本亚洲| 天堂av国产一区二区熟女人妻| 国产av不卡久久| 亚洲 欧美一区二区三区| 欧美黑人欧美精品刺激| 丁香六月欧美| 欧美性猛交黑人性爽| 99热这里只有精品一区 | 日本一二三区视频观看| 又黄又粗又硬又大视频| 国产一区二区在线av高清观看| 怎么达到女性高潮| 女生性感内裤真人,穿戴方法视频| 99久久精品一区二区三区| 久久久水蜜桃国产精品网| 91老司机精品| 19禁男女啪啪无遮挡网站| 成人一区二区视频在线观看| 久久久久久大精品| 又黄又粗又硬又大视频| 国产精品久久久久久人妻精品电影| 少妇熟女aⅴ在线视频| 日本a在线网址| 两性午夜刺激爽爽歪歪视频在线观看| 免费搜索国产男女视频| 日韩中文字幕欧美一区二区| 99精品久久久久人妻精品| 999久久久国产精品视频| 欧美日韩一级在线毛片| 亚洲精品一区av在线观看| 亚洲成a人片在线一区二区| 久久久国产欧美日韩av| 国产激情久久老熟女| 日韩欧美一区二区三区在线观看| 中文在线观看免费www的网站| 精品乱码久久久久久99久播| 国产探花在线观看一区二区| 亚洲狠狠婷婷综合久久图片| 久久精品亚洲精品国产色婷小说| 色播亚洲综合网| 日本成人三级电影网站| 欧美乱妇无乱码| 一区二区三区激情视频| av视频在线观看入口| 人人妻人人看人人澡| 婷婷六月久久综合丁香| 无人区码免费观看不卡| 国内精品一区二区在线观看| 国产精品日韩av在线免费观看| 欧洲精品卡2卡3卡4卡5卡区| 午夜福利在线在线| 男人和女人高潮做爰伦理| 久久久久免费精品人妻一区二区| 19禁男女啪啪无遮挡网站| 午夜成年电影在线免费观看| 国产精品98久久久久久宅男小说| 精品久久久久久久毛片微露脸| 九色国产91popny在线| 天天躁狠狠躁夜夜躁狠狠躁| 最近最新中文字幕大全电影3| 国内揄拍国产精品人妻在线| 男女那种视频在线观看| 视频区欧美日本亚洲| 一本综合久久免费| 亚洲五月天丁香| 亚洲国产高清在线一区二区三| 香蕉国产在线看| 久久精品人妻少妇| 欧美一级毛片孕妇| 宅男免费午夜| 天堂av国产一区二区熟女人妻| 亚洲无线在线观看| 看免费av毛片| 欧美av亚洲av综合av国产av| 国产精品香港三级国产av潘金莲| 久久九九热精品免费| 又黄又爽又免费观看的视频| 国产一区二区在线av高清观看| 黄色丝袜av网址大全| 欧美在线黄色| www国产在线视频色| 婷婷精品国产亚洲av在线| 一本一本综合久久| 老司机在亚洲福利影院| 搡老妇女老女人老熟妇| 搡老熟女国产l中国老女人| АⅤ资源中文在线天堂| 伦理电影免费视频| 国产真人三级小视频在线观看| 日本撒尿小便嘘嘘汇集6| 麻豆久久精品国产亚洲av| 中亚洲国语对白在线视频| 又黄又粗又硬又大视频| 亚洲中文av在线| 亚洲欧美日韩高清在线视频| 国产精品一及| 一本精品99久久精品77| 黄色女人牲交| 亚洲国产中文字幕在线视频| 久久中文看片网| 麻豆一二三区av精品| 亚洲欧美一区二区三区黑人| 亚洲电影在线观看av| 精品无人区乱码1区二区| 国产又色又爽无遮挡免费看| 亚洲欧美日韩高清在线视频| 在线观看美女被高潮喷水网站 | x7x7x7水蜜桃| 国产精品av视频在线免费观看| 国产综合懂色| 久久久久久大精品| 成年版毛片免费区| 精品熟女少妇八av免费久了| 国产真人三级小视频在线观看| 九色国产91popny在线| АⅤ资源中文在线天堂| 亚洲专区字幕在线| 日本一二三区视频观看| 男女午夜视频在线观看| 天堂影院成人在线观看| 变态另类丝袜制服| www日本在线高清视频| 丝袜人妻中文字幕| 午夜福利在线在线| 久久久久精品国产欧美久久久| 国产一区二区激情短视频| 韩国av一区二区三区四区| 亚洲五月天丁香| 国内少妇人妻偷人精品xxx网站 | 久久精品影院6| 特级一级黄色大片| 久久九九热精品免费| 2021天堂中文幕一二区在线观| www.熟女人妻精品国产| 国产亚洲欧美在线一区二区| 免费无遮挡裸体视频| 精品欧美国产一区二区三| 亚洲精品一区av在线观看| 高清在线国产一区| 精品久久久久久久毛片微露脸| 婷婷精品国产亚洲av| 在线观看免费午夜福利视频| 国产精品香港三级国产av潘金莲| 窝窝影院91人妻| 亚洲欧美日韩卡通动漫| 成人av一区二区三区在线看| www国产在线视频色| 88av欧美| 亚洲自偷自拍图片 自拍| 又大又爽又粗| 国产三级在线视频| 最近最新中文字幕大全电影3| 中文字幕熟女人妻在线| 又紧又爽又黄一区二区| 国产日本99.免费观看| 少妇丰满av| 国产av麻豆久久久久久久| 日日干狠狠操夜夜爽| 日本黄大片高清| 老司机午夜福利在线观看视频| 五月玫瑰六月丁香| 超碰成人久久| av女优亚洲男人天堂 | 精品国产乱子伦一区二区三区| 国产精品av久久久久免费| 国产精品98久久久久久宅男小说| 别揉我奶头~嗯~啊~动态视频| 女人高潮潮喷娇喘18禁视频| 久久人妻av系列| 日韩欧美免费精品| 成年女人看的毛片在线观看| 狂野欧美激情性xxxx| 亚洲18禁久久av| 欧美乱色亚洲激情| 给我免费播放毛片高清在线观看| 在线观看日韩欧美| 国产精品av视频在线免费观看| 中文字幕高清在线视频| 又大又爽又粗| 精品久久久久久,| 夜夜躁狠狠躁天天躁| 岛国在线免费视频观看| 精品电影一区二区在线| 国产成人啪精品午夜网站| 视频区欧美日本亚洲| 久久国产精品影院| 老司机午夜十八禁免费视频| 国产精品一及| 一个人看视频在线观看www免费 | 精品国产乱子伦一区二区三区| 免费观看人在逋| 夜夜看夜夜爽夜夜摸| 成年女人毛片免费观看观看9| 国产av不卡久久| 国语自产精品视频在线第100页| 三级毛片av免费| 精品国产美女av久久久久小说| 欧美日韩国产亚洲二区| 1024手机看黄色片| 日韩 欧美 亚洲 中文字幕| 国产伦精品一区二区三区四那| 51午夜福利影视在线观看| 不卡一级毛片| 午夜福利成人在线免费观看| 国产亚洲欧美在线一区二区| 亚洲人成网站在线播放欧美日韩| 少妇的丰满在线观看| 亚洲精品粉嫩美女一区| 久久国产精品影院| 久久久国产成人精品二区| 亚洲熟妇熟女久久| 又紧又爽又黄一区二区| www国产在线视频色| 国产午夜精品久久久久久| 一级黄色大片毛片| www国产在线视频色| 国产三级黄色录像| 中文亚洲av片在线观看爽| 91麻豆av在线| 欧美性猛交黑人性爽| 久久香蕉精品热| 少妇的丰满在线观看| 国产一区二区激情短视频| 国产精品一区二区精品视频观看| 国产成人一区二区三区免费视频网站| 宅男免费午夜| 日韩精品中文字幕看吧| 老司机午夜十八禁免费视频| 亚洲国产欧洲综合997久久,| 欧美不卡视频在线免费观看| 97超级碰碰碰精品色视频在线观看| 亚洲国产欧洲综合997久久,| 国产精品一区二区精品视频观看| 国产69精品久久久久777片 | 亚洲最大成人中文| 美女 人体艺术 gogo| 亚洲精品色激情综合| www日本在线高清视频| 99热6这里只有精品| 黄片小视频在线播放| 曰老女人黄片| 免费在线观看成人毛片| 丰满人妻一区二区三区视频av | 99久久久亚洲精品蜜臀av| 99久久成人亚洲精品观看| 亚洲欧美一区二区三区黑人| 日韩欧美国产一区二区入口| 亚洲欧美日韩高清在线视频|