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

    機(jī)器人系統(tǒng)中小麥病害識別與施藥算法研究

    2017-11-15 09:09:40刁智華刁春迎魏玉泉毋媛媛錢曉亮劉玉翠
    江蘇農(nóng)業(yè)科學(xué) 2017年17期
    關(guān)鍵詞:機(jī)器人

    刁智華+刁春迎+魏玉泉+毋媛媛+錢曉亮+劉玉翠

    摘要:為解決傳統(tǒng)施藥機(jī)械農(nóng)藥利用率低且在施藥過程中不能對作物病害進(jìn)行檢測等問題,根據(jù)常見農(nóng)作物小麥的生長規(guī)律,提出了一種適用于施藥機(jī)器人系統(tǒng)的病害識別與施藥算法。研究方法能使機(jī)器人根據(jù)在行進(jìn)過程中拍攝的作物圖片識別出作物的病害種類并根據(jù)受害面積得出受害程度;通過改進(jìn)的Hough變換算法識別出作物行;將病害信息和作物行信息作為精準(zhǔn)施藥依據(jù),在不改變噴藥壓力和流量的基礎(chǔ)上,通過中間控制器調(diào)節(jié)機(jī)器人噴桿行走速度以作物行中心線為軌跡進(jìn)行變量施藥。試驗結(jié)果表明,機(jī)器人系統(tǒng)算法工作穩(wěn)定,具有良好的穩(wěn)定性和較高的實用性。

    關(guān)鍵詞:精準(zhǔn)施藥;機(jī)器人;病害識別;作物行

    中圖分類號: TP242文獻(xiàn)標(biāo)志碼: A文章編號:1002-1302(2017)17-0192-04

    收稿日期:2017-03-09

    基金項目:國家自然科學(xué)基金(編號:61501407),河南省科技廳科技攻關(guān)項目(編號:162102110118);河南省高等學(xué)校重點科研項目(編號:15A413006);河南省高等學(xué)校青年骨干教師培養(yǎng)計劃(編號:2016GGJS-088);鄭州輕工業(yè)學(xué)院研究生科技創(chuàng)新基金資助項目(編號:2016040)。

    作者簡介:刁智華(1982—),男,河南鄭州人,博士,副教授,主要從事圖像處理及模式識別的有關(guān)研究。E-mail:diaozhua@163.com。我國是一個農(nóng)業(yè)大國,據(jù)相關(guān)研究數(shù)據(jù)統(tǒng)計,每年的農(nóng)作物感染病害面積多達(dá)2億hm2,造成巨大的經(jīng)濟(jì)損失。每年用于病害防治的農(nóng)藥約25萬t(有效利用率僅為20%~30%),真正施在病害部位的僅占0.02%左右,其余的農(nóng)藥都散布在農(nóng)田、水域、空氣中,造成嚴(yán)重的水土資源污染[1-3]。農(nóng)業(yè)操作如病害檢測、噴灑農(nóng)藥等十分繁瑣,大多數(shù)情況下作物病害識別是人工進(jìn)行的,通過人工測量或顯微鏡,但是這需要有經(jīng)驗的專業(yè)人士,而且有些病害的癥狀非常相似,通過肉眼很難分辨,難以達(dá)到精準(zhǔn)的效果[4-7];傳統(tǒng)施藥時會采取過多的防護(hù)措施,如穿防水的衣服、戴口罩和手套等,這些措施可以避免農(nóng)藥對他們造成傷害,但也不能完全避免[8-10]。在這樣的情況下,使用施藥機(jī)器人便成了最好的解決方案,因此對機(jī)器人系統(tǒng)下作物病害識別及施藥算法的研究具有重大的意義。

    自20世紀(jì)80年代發(fā)達(dá)國家開始研究農(nóng)業(yè)機(jī)器人后,相繼出現(xiàn)了采摘、嫁接、施藥機(jī)器人,但早期的施藥機(jī)器人各方面性能低、智能化程度不高,沒有針對性的病害檢測系統(tǒng)和精準(zhǔn)施藥系統(tǒng),后來由于人工智能、機(jī)器視覺等新技術(shù)的發(fā)展及應(yīng)用,有效促進(jìn)了施藥機(jī)器人的自動化和精準(zhǔn)化研發(fā)[11-13]。在國外,Arima等提出了一種多操作機(jī)器人作物信息獲取及精準(zhǔn)施藥的方法,達(dá)到了傳感、控制和數(shù)據(jù)管理3種功能的融合。Kozlowski等提出了一種機(jī)器人系統(tǒng)中精準(zhǔn)施藥自動控制方法,解決了復(fù)雜的精準(zhǔn)施藥問題。在國內(nèi),耿長興等設(shè)計了機(jī)器人系統(tǒng)下PLC控制的精準(zhǔn)施藥方法,在3種病情等級上對機(jī)器人進(jìn)行了試驗,藥液節(jié)省率達(dá)到60%以上。彭吳琦等使用MATLAB對小麥葉部圖像依次進(jìn)行圖像轉(zhuǎn)換反映病害分布情況,根據(jù)病害覆蓋率識別出作物病害程度,為本研究的病害識別方法的研究提供了可行性依據(jù)。目前,國內(nèi)外對精準(zhǔn)施藥機(jī)器人的研究尚處于研究試驗階段,環(huán)境的復(fù)雜性、作業(yè)對象形態(tài)的多義性、光照條件的不可控制性等都給機(jī)器人系統(tǒng)中信息獲取和決策實施帶來很大的困難。本研究以生長早中期的小麥常見病害為噴霧研究對象,利用機(jī)器視覺來獲取小麥的病害信息,以病害的分布及受害面積作為施藥依據(jù),最終實現(xiàn)以作物行中心線為軌跡的精準(zhǔn)施藥。

    1材料與方法

    1.1系統(tǒng)原理

    施藥機(jī)器人系統(tǒng)中小麥病害識別與施藥算法的研究由圖像采集部分、信息識別部分以及精準(zhǔn)施藥控制部分組成。圖像采集部分是通過工業(yè)相機(jī)對農(nóng)作物信息進(jìn)行采集;信息識別部分是通過hough變換提取出作物行中心線以及通過分析采集圖片上病斑的顏色、形狀、紋理等特征識別出病害的種類,進(jìn)而根據(jù)病斑所占面積得出受害程度,在得到施藥決策后,精準(zhǔn)施藥控制部分通過控制噴霧桿行走速度進(jìn)行變量噴霧。

    1.2試驗規(guī)劃

    小麥種植密度大而集中,一般按照等間距行種植,在生長早期植株間距大,葉片分布均勻,無遮擋現(xiàn)象,但很少出現(xiàn)病害信息;生長中期小麥植株間距縮小,葉片呈層次分布,行距分明,病害信息明顯易查,開始出現(xiàn)遮擋現(xiàn)象;而生長晚期小麥植株間距最小,葉片遮擋嚴(yán)重,病害信息雜亂不易檢測,基本無行距。鑒于以上分析,生長早中期的小麥行距分明,葉片遮擋率小,便于施藥,病害的癥狀明顯,易于區(qū)分,所以研究選擇生長早中期的小麥進(jìn)行試驗。本研究在河南農(nóng)業(yè)大學(xué)毛莊農(nóng)業(yè)科研實驗基地進(jìn)行,利用試驗基地現(xiàn)存條件,結(jié)合小麥種植模式及葉面積分布構(gòu)造了適用于機(jī)器人系統(tǒng)的病害識別與施藥算法。

    2系統(tǒng)設(shè)計

    2.1圖像采集與預(yù)處理

    為實現(xiàn)精準(zhǔn)施藥首先要準(zhǔn)確獲取作物的病害信息,即對作物病害部位進(jìn)行準(zhǔn)確的衡量。作物在感染病害后葉部形狀、顏色、紋理會出現(xiàn)變化,且感染不同的病害會出現(xiàn)不同的癥狀,系統(tǒng)根據(jù)這些癥狀的不同區(qū)分病害的種類。系統(tǒng)采用西安艾菲特光電技術(shù)有限公司生產(chǎn)的MV-VDM033SC型號的工業(yè)攝像機(jī),分辨率為640×480,輸出為RGB彩色圖像,相機(jī)在試驗中每隔2 s采集一次圖像,為保證圖像處理的精確度每次采集2張圖片,一張近距離圖像,一張遠(yuǎn)距離圖像,分別用于病害識別和作物行檢測。

    2.1.1灰度變換采集的圖像由于環(huán)境、光照度等的影響會存在噪聲及其他物體的干擾,因此要對圖像進(jìn)行一定的預(yù)處理,包括灰度變換、中值濾波等。采集的圖片中,小麥的正常葉片及陰影部分呈現(xiàn)綠色,而感染病害的葉片呈現(xiàn)異于綠色的白色、黃白色等,可以利用這一特性使用超綠顏色特征2G-R-B 將感染病害的小麥葉部病斑提取出來。超綠算法表達(dá)式為式(1)所示。endprint

    exg=0,2G-R-B<0

    255,2G-R-B>255

    2G-R-B,0≤2G-R-B≤255。(1)

    圖1-a、圖1-c為普通灰度化圖像,圖1-b、圖1-d為超綠灰度化圖像。通過分析可知,超綠灰度化后的圖1-b更能突出作物行的走勢,可清晰地觀察到作物行的形態(tài)特征;圖1-d更加突出了小麥葉片的病斑位置、大小、形狀等形態(tài)特征,減少了冗余的復(fù)雜背景信息。

    2.1.2中值濾波為對灰度化后的圖像除去噪聲對病斑及作物行提取的影響,采用3×3的模板對采集圖像進(jìn)行中值濾波,即對模板的9個灰度值進(jìn)行從小到大或從大到小的排序,選擇中間的灰度值來替代處理的灰度值。

    2.2圖像分割

    對預(yù)處理后的圖像進(jìn)行圖像分割是識別病害必不可少的一部分,分割的效果直接影響分析的結(jié)果,因此選擇的分割方法對病害識別有很大的影響。由于圖像中非病害部位與病害部位存在較大的顏色差異,閾值分割因簡單快速是解決此類問題的最優(yōu)方法。為減少閾值分割參與數(shù)據(jù),提高圖像處理速度,本研究選擇一種改進(jìn)的適用于本次圖像分割的自動選擇閾值的最大類間方差比閾值分割方法,具體算法如下:

    (1)找出整個圖像像素最多點對應(yīng)的灰度級K,并求出K之前的灰度級像素總數(shù)N以及所有像素的灰度均值μ和分布概率p0,p1,p2,…,pk;

    μ=1N∑mi=1∑nj=1f(i,j)。(2)

    式中:f(i,j)為點(i,j)的像素灰度值。

    (2)給定一個非零初始閾值T=T0,將圖像分為C1、C2兩類;

    (3)計算灰度均值μ1、μ2和兩類方差δ12、δ22;

    μi=1NCi∑(x,y)∈Cif(x,y),(i=1,2);(3)

    δi2=∑(x,y)∈Ci[f(x,y)-μi]2,(i=1,2)。(4)

    (4)計算C1、C2的發(fā)生概率P1、P2;

    P1=∑Ti=0pi,P2=1-P1。(5)

    (5)計算類內(nèi)方差δm2和類間方差δb2;

    δm2=P1·δ12+P2·δ22;(6)

    δb2=P1·(μ1-μ)2+P2·(μ2-μ)2。(7)

    令T=δb2/δm2,重復(fù)步驟2~5,第二次運算時取T0=1/2(μ1+μ2),當(dāng)且僅當(dāng)T=Tmax時的T為最佳閾值。本研究選擇的白粉病病斑分割圖片經(jīng)計算后得出閾值為200時分割效果為最佳,分割病斑后的圖像效果如圖2所示。

    對分割后的病害區(qū)域作面積統(tǒng)計,根據(jù)受害面積比例設(shè)置4種不同的受害程度,結(jié)果如表1所示。

    2.3病害識別

    病害識別中顏色是最明顯、直觀的特征,根據(jù)顏色特征能夠分辨出一部分病害,紋理和形狀特征可作為分辨類別的輔助依據(jù),通過分析病害的顏色、紋理及形狀特征能夠分析出病害的種類。本研究根據(jù)病斑顏色特征的顏色矩和標(biāo)準(zhǔn)顏色直方圖,紋理特征的角二階矩(ASM)、對比度(CON)及相關(guān)度(COR)和形狀特征的不變矩(μp)、偏心率(E)的計算方法[14],對小麥幾種常見病害進(jìn)行了識別,病害識別結(jié)果如表2所示。

    由于圖像分割后的病害特征明顯,識別率達(dá)到85%以上;但光照條件或相似病害的特征參數(shù)數(shù)據(jù)重疊存在一定的影響,因此病害識別存在少量的誤差。

    2.4作物行識別

    作物行的背景信息復(fù)雜,直接提取具有較大的困難,且不能保證提取的精度,因此要進(jìn)行一定的圖像處理步驟。二值化是識別作物行的第一步,不僅簡化了作物行信息,而且可以減少圖像處理時間;二值化后進(jìn)行細(xì)化處理突出作物行骨架,在提取作物行時只保留中間2行,最后通過Hough變換擬合直線檢測出作物行的中心線。

    2.4.1二值化為了播種和收割等的方便,大多數(shù)的農(nóng)作物都是以等間距行的形式種植,這也為農(nóng)作物的施藥提供了方便。提取作物行信息首先要對采集的作物行圖片進(jìn)行二值化處理,采用膨脹和腐蝕算法消除二值化后的圖像噪聲,本研究采用3×3結(jié)構(gòu)元素進(jìn)行圖像膨脹,采用5×1的結(jié)構(gòu)元素進(jìn)行圖像腐蝕,二值化后圖像如圖3-a所示。

    2.4.2提取作物行二值化后的圖像已經(jīng)得到作物行的大致輪廓,為清晰地提取出作物行的骨架信息,對圖像采取細(xì)化處理,提取作物行中心骨骼,減少冗余的部分突出形狀特點,細(xì)化后圖像如圖3-b所示??紤]到作物行的間距一般是相等的,只要對準(zhǔn)其中2行,其他行按等距排列就可以實現(xiàn)所有行的對準(zhǔn),一般選取圖像的中心2行為參考行。

    2.4.3改進(jìn)的Hough變換Hough變換具有較好的抗干擾性和魯棒性,是用于直線檢測的一種常用方法,但傳統(tǒng)的Hough變換具有運算量大、斷線錯連、運算時間長等缺陷,在快速運算中失去了本來的應(yīng)用意義。本研究提出一種改進(jìn)的基于定向處理的快速Hough變換,減少了傳統(tǒng)Hough變換巨大的運算量,降低了計算復(fù)雜度。

    Hough變換將作物行擬合為直線,其基本思想是利用點線的對偶性把圖像中的直線檢測問題轉(zhuǎn)換為對參數(shù)空間中點的檢測問題,通過參數(shù)空間中點的累加統(tǒng)計完成直線檢測。提高Hough變換的速度有2種方法,一是減少參與的點數(shù),二是縮小計算的空間,本文提出的基于定向處理的快速Hough變換是這2種方法的結(jié)合。為了減少計算的復(fù)雜度采用直線極坐標(biāo)方程表示,如式(8)所示。

    ρ=xcosθ+ysinθ。(8)

    式中:ρ是原點到直線的距離,θ是x軸與直線的夾角。計算時首先把整體圖像量化為幾個子區(qū)域,在細(xì)化圖像中近直線段上隨機(jī)選取兩點(xm,ym)、(xn,yn)確定一條直線,參數(shù)如式(9)、(10)所示。

    θ=arctanym-ynxm-xn;(9)

    ρ=xmcosθ+ymsinθ。(10)

    把每段近直線中所有像素點組成一集合A,隨機(jī)取集合A中兩點求斜率,如不滿足近直線要求,重新選擇。由于只提取作物行圖像的中間2行,因此可以把θ的取值范圍壓縮為30°~150°,根據(jù)ρ=x2+y2sin(θ+)計算ρ的值,其中θ=sin-1(x/x2+y2)。最后根據(jù)ρ和θ的值對累加數(shù)組進(jìn)行累加,得到共線點的個數(shù),由ρ和θ的取值范圍及它們的分辨率可以確定累加器的大小,從而檢測直線。Hough變換后的圖像如圖3-d所示。endprint

    2.5施藥決策

    通過對采集的圖像進(jìn)行灰度化轉(zhuǎn)換、中值濾波、分割、作物行提取等一系列處理后,將得到的小麥病害種類、分布、受病害程度以及作物行信息作為施藥依據(jù),根據(jù)識別出的病害類型選擇相應(yīng)的藥液,4種病害等級相應(yīng)地設(shè)置4種不同的施藥速度,使其保持在0.15~0.45 m/s的變化范圍內(nèi),在不改變噴藥壓力(0.8 MPa)和噴頭流量(6.5 L/min)的情況下,只調(diào)節(jié)機(jī)器人噴桿行走的速度實現(xiàn)自動變量施藥。為使精準(zhǔn)施藥機(jī)器人在一次噴霧中能夠有效覆蓋盡量多的小麥植株,

    研究采取可移動六噴嘴式噴霧系統(tǒng),每個噴嘴位于作物行中心線的正上方向下垂直噴霧,以噴嘴距離工作區(qū)域為300 mm時對應(yīng)120 mm×120 mm的噴幅為標(biāo)準(zhǔn)。實際施藥時由于相機(jī)與作物行底部中心存在偏差,因此可根據(jù)偏差信息進(jìn)行適量的調(diào)整。

    3算法實現(xiàn)形式

    本研究的機(jī)器人系統(tǒng)中病害識別與精準(zhǔn)施藥的實現(xiàn)形式如圖4所示,首先根據(jù)相機(jī)獲取小麥圖像,對圖像進(jìn)行預(yù)處理;然后進(jìn)行作物行檢測及病害檢測,根據(jù)分割病斑后圖像計算出感染病害面積,繼而劃分病害等級及施藥等級;最后將病害及作物行信息發(fā)送給精準(zhǔn)施藥控制系統(tǒng),由精準(zhǔn)施藥控制系統(tǒng)給控制閥發(fā)送命令控制藥箱的啟閉并控制壓力泵給藥箱施加恒等的壓力,藥液流出時通過流量計檢測藥液流出量的多少,并把信息傳遞給精準(zhǔn)施藥控制系統(tǒng)進(jìn)行信息負(fù)反饋,從而進(jìn)行相應(yīng)的調(diào)整,通過流量計實時監(jiān)測的噴頭進(jìn)行噴藥,實現(xiàn)了噴頭的自動對準(zhǔn)和藥液流量的自動控制。

    4結(jié)論

    本研究提出的機(jī)器人系統(tǒng)中病害識別與施藥算法有效地減少了農(nóng)藥的浪費,通過機(jī)器視覺來獲取病害信息,采用基于小麥葉片顏色、形狀及紋理信息的圖像處理方法實現(xiàn)了小麥病害的精準(zhǔn)識別,并根據(jù)受害面積建立了病害病情等級評價標(biāo)準(zhǔn),為后續(xù)的精準(zhǔn)施藥系統(tǒng)提供了施藥依據(jù)。通過VC++6.0的開發(fā)語言實現(xiàn)小麥常見病害的分割、識別及提取作物行等算法以及對精準(zhǔn)施藥系統(tǒng)的控制,一幅圖像300 ms左右的處理時間滿足了實際的需求,在精度和準(zhǔn)確度上都達(dá)到了較好的效果,且使用軟件處理極大地節(jié)約了成本,為精準(zhǔn)施藥機(jī)器人的廣泛使用提供了可能。但是,研究過程中由于光照、葉片遮擋等的影響,獲取的病害信息仍然存在一定的誤差,從采集圖像到處理結(jié)束中間算法步驟稍多,因此如何準(zhǔn)確地獲取病害信息、提高圖像處理時間仍然是需要解決的問題。

    參考文獻(xiàn):

    [1]趙春江,陳天恩,陳立平,等. 遷飛性害蟲精準(zhǔn)施藥決策分析方法研究[J]. 農(nóng)業(yè)工程學(xué)報,2008,24(增刊2):149-155.

    [2]傅錫敏,薛新宇. 基于我國施藥技術(shù)與裝備現(xiàn)狀的發(fā)展思路[J]. 中國農(nóng)機(jī)化,2008(6):72-76.

    [3]You G,Wang X,Yang S,et al. Precision spraying system of crops disease stress based on acoustic emission[J]. Journal of Computers,2011,6(4):635-642.

    [4]邵慶,張楠,路陽. 小麥病害圖像識別處理及形狀特征提取研究[J]. 農(nóng)機(jī)化研究,2013(8):35-37.

    [5]刁智華,王歡,宋寅卯,等. 復(fù)雜背景下棉花病葉害螨圖像分割方法[J]. 農(nóng)業(yè)工程學(xué)報,2013,29(5):147-152.

    [6]王獻(xiàn)鋒,張善文,王震,等. 基于葉片圖像和環(huán)境信息的黃瓜病害識別方法[J]. 農(nóng)業(yè)工程學(xué)報,2014,30(14):148-153.

    [7]王翔宇,溫皓杰,李鑫星,等. 農(nóng)業(yè)主要病害檢測與預(yù)警技術(shù)研究進(jìn)展分析[J]. 農(nóng)業(yè)機(jī)械學(xué)報,2016,47(9):266-277.

    [8]劉雪美,李揚,李明,等. 噴桿噴霧機(jī)精確對靶施藥系統(tǒng)設(shè)計與試驗[J]. 農(nóng)業(yè)機(jī)械學(xué)報,2016,47(3):37-44.

    [9]Dammer K. Real-time variable-rate herbicide application for weed control in carrots[J]. Weed Res,2016,56(3):237-246.

    [10]Tackenberg M,Volkmar C,Dammer K. Sensor-based variable-rate fungicide application in winter wheat[J]. Pest Manag Sci,2016,72(10):1888-1896.

    [11]Aishwarya B V,Archana G,Umayal C. Agriculture robotic vehicle based pesticide sprayer with efficiency optimization[C]. Technological Innovation in ICT for Agriculture and Rural Development,2015:59-65.

    [12]Neudecker S,Bruns C,Gerbers R,et al. A new robotic spray technology for generative manufacturing of complex concrete structures without formwork[J]. Procedia Cirp,2016,43:333-338.

    [13]刁智華,趙明珍,宋寅卯,等. 基于機(jī)器視覺的玉米精準(zhǔn)施藥系統(tǒng)作物行識別算法及系統(tǒng)實現(xiàn)[J]. 農(nóng)業(yè)工程學(xué)報,2015,31(7):47-52.

    [14]劉濤,仲曉春,孫成明,等. 基于計算機(jī)視覺的水稻葉部病害識別研究[J]. 中國農(nóng)業(yè)科學(xué),2014,47(4):664-674.何玲,陳長喜. 生豬屠宰自動化監(jiān)管系統(tǒng)的物聯(lián)網(wǎng)構(gòu)建[J]. 江蘇農(nóng)業(yè)科學(xué),2017,45(17):196-199.endprint

    猜你喜歡
    機(jī)器人
    “壞”機(jī)器人“大叛亂”
    機(jī)器人
    超酷機(jī)器人
    瘋狂機(jī)器人
    機(jī)器人崛起
    少先隊活動(2021年5期)2021-07-22 09:00:00
    機(jī)器人,讓未來走近你
    金橋(2019年10期)2019-08-13 07:15:16
    擁抱未來——聚焦機(jī)器人
    大社會(2016年3期)2016-05-04 03:41:07
    機(jī)器人來幫你
    認(rèn)識機(jī)器人
    機(jī)器人來啦
    欧美变态另类bdsm刘玥| 久久久久久久久久人人人人人人| 亚洲精品aⅴ在线观看| 久久久久精品久久久久真实原创| 青春草亚洲视频在线观看| 国产淫片久久久久久久久| 成年av动漫网址| 秋霞伦理黄片| 天堂中文最新版在线下载| 看非洲黑人一级黄片| 搡老乐熟女国产| 水蜜桃什么品种好| 午夜视频国产福利| 国产女主播在线喷水免费视频网站| 十分钟在线观看高清视频www | 美女国产视频在线观看| 嫩草影院入口| 日日啪夜夜爽| 成人亚洲精品一区在线观看 | 久久久久国产精品人妻一区二区| 国产又色又爽无遮挡免| 男人舔奶头视频| 色哟哟·www| 国产精品熟女久久久久浪| 尤物成人国产欧美一区二区三区| 欧美日韩一区二区视频在线观看视频在线| 色婷婷av一区二区三区视频| 最后的刺客免费高清国语| 丝瓜视频免费看黄片| 大香蕉久久网| 看非洲黑人一级黄片| 国产黄片视频在线免费观看| 中文字幕亚洲精品专区| 国产成人精品福利久久| 日韩精品有码人妻一区| 亚洲色图综合在线观看| 免费黄网站久久成人精品| 人妻制服诱惑在线中文字幕| 亚洲欧美成人综合另类久久久| 啦啦啦中文免费视频观看日本| 成人二区视频| h日本视频在线播放| 久久久久久九九精品二区国产| 一区二区三区精品91| 在线观看三级黄色| 亚洲av电影在线观看一区二区三区| 亚洲人与动物交配视频| 97超碰精品成人国产| 亚洲精品亚洲一区二区| 91久久精品国产一区二区三区| 国产人妻一区二区三区在| 99热这里只有是精品50| 日本黄大片高清| 自拍欧美九色日韩亚洲蝌蚪91 | 国产色爽女视频免费观看| 久久久久久久久久人人人人人人| 国产精品一及| 一本一本综合久久| 亚洲人成网站高清观看| av网站免费在线观看视频| 高清av免费在线| 国内精品宾馆在线| 久久婷婷青草| 只有这里有精品99| 国产av国产精品国产| 国产精品蜜桃在线观看| 大陆偷拍与自拍| 美女国产视频在线观看| 国产老妇伦熟女老妇高清| 亚洲三级黄色毛片| 亚洲精品乱码久久久久久按摩| 免费av中文字幕在线| 在线免费观看不下载黄p国产| 国产av一区二区精品久久 | 王馨瑶露胸无遮挡在线观看| 尤物成人国产欧美一区二区三区| 亚洲第一区二区三区不卡| 国产欧美另类精品又又久久亚洲欧美| 水蜜桃什么品种好| 在线精品无人区一区二区三 | 国产成人a区在线观看| 大话2 男鬼变身卡| 丰满乱子伦码专区| 一级av片app| 免费大片黄手机在线观看| 亚洲成人中文字幕在线播放| av国产久精品久网站免费入址| 涩涩av久久男人的天堂| 22中文网久久字幕| 欧美性感艳星| 婷婷色av中文字幕| 成人高潮视频无遮挡免费网站| 成人黄色视频免费在线看| 久久97久久精品| 蜜桃久久精品国产亚洲av| 波野结衣二区三区在线| 成人国产av品久久久| 久久久久视频综合| 丝袜脚勾引网站| 午夜免费男女啪啪视频观看| 精华霜和精华液先用哪个| 在线观看免费高清a一片| 高清在线视频一区二区三区| 肉色欧美久久久久久久蜜桃| 亚洲精品国产成人久久av| 热re99久久精品国产66热6| 色5月婷婷丁香| 精品人妻视频免费看| 啦啦啦视频在线资源免费观看| 亚洲精品乱码久久久v下载方式| 国产av码专区亚洲av| 国产欧美另类精品又又久久亚洲欧美| 一级黄片播放器| 国产伦在线观看视频一区| 精品人妻视频免费看| 啦啦啦视频在线资源免费观看| 高清毛片免费看| 精品国产一区二区三区久久久樱花 | 午夜福利网站1000一区二区三区| 欧美少妇被猛烈插入视频| 成年人午夜在线观看视频| 久久精品熟女亚洲av麻豆精品| 在线观看av片永久免费下载| 黄色怎么调成土黄色| 成年av动漫网址| 欧美三级亚洲精品| 99九九线精品视频在线观看视频| 国产在线一区二区三区精| 日本-黄色视频高清免费观看| 精品久久久久久久久av| 嘟嘟电影网在线观看| 亚洲色图av天堂| 亚洲精品乱久久久久久| 秋霞在线观看毛片| 老司机影院毛片| 亚洲国产欧美在线一区| 国产精品久久久久久精品古装| 三级经典国产精品| 高清av免费在线| 日本黄色片子视频| 国产精品一区二区在线观看99| 一本—道久久a久久精品蜜桃钙片| 国产精品久久久久久精品电影小说 | 亚洲av.av天堂| 一本一本综合久久| 一个人看视频在线观看www免费| 视频区图区小说| 免费观看的影片在线观看| 成人午夜精彩视频在线观看| 日韩不卡一区二区三区视频在线| 欧美精品亚洲一区二区| 国产精品99久久久久久久久| 精品午夜福利在线看| 秋霞在线观看毛片| 大陆偷拍与自拍| 一级a做视频免费观看| 少妇熟女欧美另类| 波野结衣二区三区在线| 久久ye,这里只有精品| 女的被弄到高潮叫床怎么办| 夜夜看夜夜爽夜夜摸| 97在线视频观看| 嘟嘟电影网在线观看| 国产有黄有色有爽视频| 深爱激情五月婷婷| 观看免费一级毛片| 日日啪夜夜撸| 欧美成人一区二区免费高清观看| 国产伦精品一区二区三区四那| 国产精品无大码| 亚洲高清免费不卡视频| 国产淫语在线视频| 免费播放大片免费观看视频在线观看| 乱系列少妇在线播放| 亚洲欧美日韩东京热| 免费不卡的大黄色大毛片视频在线观看| 狂野欧美激情性bbbbbb| 欧美区成人在线视频| 亚洲av福利一区| 性色av一级| 看十八女毛片水多多多| 国产精品欧美亚洲77777| 晚上一个人看的免费电影| 国产高清不卡午夜福利| 国国产精品蜜臀av免费| 99久久精品国产国产毛片| 中文字幕久久专区| 国产黄色免费在线视频| 搡老乐熟女国产| 亚洲图色成人| 亚洲真实伦在线观看| 一区二区三区免费毛片| 国产成人一区二区在线| 色吧在线观看| 韩国高清视频一区二区三区| 日本-黄色视频高清免费观看| 欧美一区二区亚洲| 夜夜爽夜夜爽视频| 国产伦精品一区二区三区视频9| 精品一区二区三区视频在线| 少妇高潮的动态图| 激情五月婷婷亚洲| 欧美一级a爱片免费观看看| 亚洲精品日韩在线中文字幕| 美女福利国产在线 | 欧美日韩视频精品一区| 国语对白做爰xxxⅹ性视频网站| 日韩一区二区视频免费看| 观看av在线不卡| 最后的刺客免费高清国语| 日韩欧美一区视频在线观看 | 亚洲综合色惰| 国产视频内射| 国产av码专区亚洲av| 国产高清国产精品国产三级 | 菩萨蛮人人尽说江南好唐韦庄| 啦啦啦中文免费视频观看日本| 国产老妇伦熟女老妇高清| 欧美zozozo另类| 卡戴珊不雅视频在线播放| 51国产日韩欧美| 在线观看免费视频网站a站| 十八禁网站网址无遮挡 | 99久久精品一区二区三区| 大陆偷拍与自拍| 国产精品秋霞免费鲁丝片| 国产黄片视频在线免费观看| 国产乱来视频区| 日韩在线高清观看一区二区三区| 男女边摸边吃奶| 丰满乱子伦码专区| 狂野欧美白嫩少妇大欣赏| 午夜免费观看性视频| 中文天堂在线官网| 午夜激情久久久久久久| 国产人妻一区二区三区在| 国产精品国产av在线观看| 男人舔奶头视频| 深爱激情五月婷婷| 18禁动态无遮挡网站| 女人久久www免费人成看片| 日韩不卡一区二区三区视频在线| av国产久精品久网站免费入址| 涩涩av久久男人的天堂| 我的老师免费观看完整版| 多毛熟女@视频| 久久精品久久久久久久性| 不卡视频在线观看欧美| 国产 精品1| 联通29元200g的流量卡| 欧美精品一区二区大全| 黄片wwwwww| 人妻少妇偷人精品九色| 国产成人精品久久久久久| 91久久精品国产一区二区成人| 国产精品人妻久久久久久| av免费在线看不卡| 老熟女久久久| 久久精品久久久久久噜噜老黄| 在线免费观看不下载黄p国产| 日本免费在线观看一区| 亚洲成人手机| 99re6热这里在线精品视频| 午夜激情久久久久久久| 中文字幕人妻熟人妻熟丝袜美| 亚洲欧美日韩无卡精品| 久久99热6这里只有精品| 1000部很黄的大片| 日本一二三区视频观看| 99视频精品全部免费 在线| 91在线精品国自产拍蜜月| 久久久欧美国产精品| 一本—道久久a久久精品蜜桃钙片| 亚洲国产成人一精品久久久| 亚洲欧美日韩无卡精品| 国产精品国产三级国产专区5o| 亚洲欧美精品自产自拍| 91久久精品国产一区二区成人| 91精品一卡2卡3卡4卡| 久久人人爽人人爽人人片va| 国产 一区 欧美 日韩| 久久精品久久久久久噜噜老黄| 国产乱来视频区| 国产高潮美女av| 18禁在线无遮挡免费观看视频| 亚洲国产精品999| 国产精品三级大全| 97在线人人人人妻| 久久精品久久久久久噜噜老黄| 久久99热6这里只有精品| 亚洲av欧美aⅴ国产| 亚洲va在线va天堂va国产| www.av在线官网国产| 激情 狠狠 欧美| 一级毛片aaaaaa免费看小| 日日撸夜夜添| 中文在线观看免费www的网站| 黄片wwwwww| 欧美丝袜亚洲另类| 亚洲精品久久久久久婷婷小说| 最后的刺客免费高清国语| 色综合色国产| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产精品999| 欧美少妇被猛烈插入视频| 国产色婷婷99| 在线观看三级黄色| 男人添女人高潮全过程视频| 一级片'在线观看视频| av在线app专区| 在线观看一区二区三区| 午夜精品国产一区二区电影| 亚洲美女黄色视频免费看| 亚洲精品第二区| 欧美精品人与动牲交sv欧美| 国产 精品1| av在线蜜桃| 春色校园在线视频观看| 日本-黄色视频高清免费观看| 欧美xxxx黑人xx丫x性爽| 欧美激情极品国产一区二区三区 | 国产在线免费精品| 久久 成人 亚洲| 国产一区二区三区综合在线观看 | 国产av精品麻豆| 国产精品久久久久久av不卡| 成人免费观看视频高清| 亚洲aⅴ乱码一区二区在线播放| 国产探花极品一区二区| 国产高清有码在线观看视频| 久久久久国产网址| 免费观看无遮挡的男女| 久久久久久久久久久免费av| 18禁在线无遮挡免费观看视频| 色综合色国产| 亚洲成人手机| 99视频精品全部免费 在线| 大片免费播放器 马上看| 女人久久www免费人成看片| 一区二区三区四区激情视频| 伊人久久精品亚洲午夜| 啦啦啦中文免费视频观看日本| 国产精品99久久久久久久久| 久久精品久久精品一区二区三区| 欧美激情极品国产一区二区三区 | av.在线天堂| 国产免费一区二区三区四区乱码| 欧美区成人在线视频| 国模一区二区三区四区视频| 肉色欧美久久久久久久蜜桃| 女性被躁到高潮视频| 国产在线免费精品| 欧美xxxx黑人xx丫x性爽| 黄色日韩在线| 少妇猛男粗大的猛烈进出视频| 十八禁网站网址无遮挡 | 中文字幕亚洲精品专区| 一区二区av电影网| 老女人水多毛片| 亚洲综合色惰| 老女人水多毛片| 丰满人妻一区二区三区视频av| av在线app专区| 99久久中文字幕三级久久日本| videos熟女内射| 99热网站在线观看| 国产成人a区在线观看| 黑人高潮一二区| 成人美女网站在线观看视频| 欧美精品国产亚洲| 赤兔流量卡办理| 欧美另类一区| 亚洲国产精品999| 天堂8中文在线网| 亚洲,欧美,日韩| 黑人猛操日本美女一级片| 青春草视频在线免费观看| 久久人人爽人人爽人人片va| 国产午夜精品一二区理论片| 中国国产av一级| 熟女av电影| 免费大片18禁| 九草在线视频观看| 能在线免费看毛片的网站| 欧美日韩综合久久久久久| 日韩免费高清中文字幕av| 日本-黄色视频高清免费观看| 亚洲av免费高清在线观看| 欧美xxⅹ黑人| 国产成人aa在线观看| 这个男人来自地球电影免费观看 | 国产精品久久久久成人av| 欧美日韩一区二区视频在线观看视频在线| 国产老妇伦熟女老妇高清| 午夜福利高清视频| 久久久久性生活片| 成人黄色视频免费在线看| 欧美极品一区二区三区四区| 久久99热这里只有精品18| 青春草亚洲视频在线观看| 一级毛片aaaaaa免费看小| 亚洲国产欧美人成| 中国美白少妇内射xxxbb| 一级毛片久久久久久久久女| 免费大片黄手机在线观看| 久久 成人 亚洲| 亚洲欧美日韩东京热| 国产高清国产精品国产三级 | 亚洲成人中文字幕在线播放| 国产亚洲最大av| 最近2019中文字幕mv第一页| 亚洲久久久国产精品| 97超视频在线观看视频| 国产精品久久久久久精品古装| 亚洲av国产av综合av卡| 精品人妻熟女av久视频| 久久久久久久久久久免费av| 超碰av人人做人人爽久久| videos熟女内射| 日韩中字成人| 97精品久久久久久久久久精品| 亚洲电影在线观看av| 免费观看无遮挡的男女| 成人午夜精彩视频在线观看| 少妇人妻久久综合中文| 狂野欧美激情性bbbbbb| 黄色一级大片看看| 国产精品精品国产色婷婷| 日韩视频在线欧美| 99热网站在线观看| 18禁动态无遮挡网站| 婷婷色麻豆天堂久久| 国产亚洲精品久久久com| 日本黄色日本黄色录像| 久久久精品94久久精品| 丝瓜视频免费看黄片| 久久国产乱子免费精品| 欧美精品亚洲一区二区| 熟女电影av网| 国产免费视频播放在线视频| 日韩成人伦理影院| 日本猛色少妇xxxxx猛交久久| 精品人妻一区二区三区麻豆| 亚洲成人中文字幕在线播放| 成人二区视频| 免费高清在线观看视频在线观看| a 毛片基地| 久久久久久九九精品二区国产| 国产精品熟女久久久久浪| 你懂的网址亚洲精品在线观看| 五月天丁香电影| 蜜臀久久99精品久久宅男| 特大巨黑吊av在线直播| 各种免费的搞黄视频| 赤兔流量卡办理| 国产黄片视频在线免费观看| 在线天堂最新版资源| 国产美女午夜福利| 大话2 男鬼变身卡| 偷拍熟女少妇极品色| 久久久久久久久久久免费av| 久久久久网色| 日本av免费视频播放| 午夜精品国产一区二区电影| 亚洲自偷自拍三级| 亚洲国产精品999| 久久综合国产亚洲精品| 久久亚洲国产成人精品v| 亚洲av电影在线观看一区二区三区| 日日摸夜夜添夜夜添av毛片| 成人18禁高潮啪啪吃奶动态图 | 另类亚洲欧美激情| 亚洲自偷自拍三级| 直男gayav资源| 少妇被粗大猛烈的视频| 超碰97精品在线观看| 国产黄色免费在线视频| 精品人妻熟女av久视频| 午夜福利网站1000一区二区三区| 婷婷色综合www| 国产av一区二区精品久久 | 亚洲精品日本国产第一区| 肉色欧美久久久久久久蜜桃| 97热精品久久久久久| 亚洲色图av天堂| 女人十人毛片免费观看3o分钟| 中国美白少妇内射xxxbb| 久久ye,这里只有精品| 小蜜桃在线观看免费完整版高清| xxx大片免费视频| 国产精品麻豆人妻色哟哟久久| 亚洲伊人久久精品综合| 免费在线观看成人毛片| 这个男人来自地球电影免费观看 | 不卡视频在线观看欧美| 我要看日韩黄色一级片| 午夜视频国产福利| 国产精品欧美亚洲77777| a级一级毛片免费在线观看| 啦啦啦在线观看免费高清www| 久久毛片免费看一区二区三区| 久久精品久久精品一区二区三区| 美女脱内裤让男人舔精品视频| 能在线免费看毛片的网站| 成人亚洲欧美一区二区av| 国产人妻一区二区三区在| 亚洲精品日韩av片在线观看| 亚洲av不卡在线观看| 日本午夜av视频| 涩涩av久久男人的天堂| 国产精品精品国产色婷婷| 免费av不卡在线播放| 国产 精品1| 欧美bdsm另类| 国产欧美日韩精品一区二区| 久久6这里有精品| 身体一侧抽搐| 国产v大片淫在线免费观看| 免费不卡的大黄色大毛片视频在线观看| 国产精品不卡视频一区二区| 性色avwww在线观看| 国产极品天堂在线| 中文在线观看免费www的网站| 成人影院久久| 大陆偷拍与自拍| 人妻系列 视频| 日韩制服骚丝袜av| 欧美变态另类bdsm刘玥| 亚洲高清免费不卡视频| 日日啪夜夜撸| 国产深夜福利视频在线观看| 精品一区在线观看国产| 久久国产精品大桥未久av | 国产伦精品一区二区三区四那| 一二三四中文在线观看免费高清| 久久热精品热| 亚洲美女搞黄在线观看| 国产亚洲最大av| 亚洲成人av在线免费| 免费在线观看成人毛片| 色吧在线观看| 精品少妇黑人巨大在线播放| 中文乱码字字幕精品一区二区三区| 免费观看a级毛片全部| 欧美97在线视频| 91在线精品国自产拍蜜月| 国产女主播在线喷水免费视频网站| 欧美丝袜亚洲另类| 五月伊人婷婷丁香| 熟妇人妻不卡中文字幕| 亚州av有码| 午夜免费男女啪啪视频观看| 久久6这里有精品| 亚洲电影在线观看av| a 毛片基地| 91aial.com中文字幕在线观看| 99热这里只有精品一区| 亚洲精品色激情综合| 久久久欧美国产精品| 青春草亚洲视频在线观看| 伦精品一区二区三区| 91在线精品国自产拍蜜月| 国产黄片美女视频| 我要看黄色一级片免费的| 22中文网久久字幕| 国产精品一二三区在线看| 成人黄色视频免费在线看| 亚洲av在线观看美女高潮| 亚洲一级一片aⅴ在线观看| 联通29元200g的流量卡| 久久久久人妻精品一区果冻| 久久鲁丝午夜福利片| 一二三四中文在线观看免费高清| 插逼视频在线观看| 国产精品一区www在线观看| 日日啪夜夜撸| 狠狠精品人妻久久久久久综合| 国产真实伦视频高清在线观看| 欧美高清性xxxxhd video| 丰满人妻一区二区三区视频av| 欧美97在线视频| 美女xxoo啪啪120秒动态图| 久久久久久久久久久免费av| 一本一本综合久久| 2022亚洲国产成人精品| 中文字幕制服av| 国产欧美另类精品又又久久亚洲欧美| 国产精品国产av在线观看| 一区二区三区四区激情视频| 插逼视频在线观看| 狂野欧美激情性bbbbbb| 日本色播在线视频| 国产亚洲av片在线观看秒播厂| 国产亚洲欧美精品永久| 久久av网站| 大片免费播放器 马上看| 成人毛片60女人毛片免费| 成年免费大片在线观看| 日本wwww免费看| 久久精品熟女亚洲av麻豆精品| 深爱激情五月婷婷| 免费av不卡在线播放| 丝袜脚勾引网站| 日本午夜av视频| 日日啪夜夜撸| 午夜精品国产一区二区电影| 日韩大片免费观看网站| 国产伦在线观看视频一区| 简卡轻食公司| 十分钟在线观看高清视频www | 日本黄色日本黄色录像| 久久久久国产精品人妻一区二区| 爱豆传媒免费全集在线观看| 五月开心婷婷网| 热re99久久精品国产66热6|