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

    基于自適應(yīng)差分進化算法擬合圓的樹干胸徑測量方法

    2018-09-17 06:49:34胡春華李萍萍金成磊
    農(nóng)業(yè)機械學(xué)報 2018年9期
    關(guān)鍵詞:胸徑樹干差分

    胡春華 李萍萍 金成磊 朱 樂

    (1.南京林業(yè)大學(xué)信息科學(xué)技術(shù)學(xué)院, 南京 210037; 2.南京林業(yè)大學(xué)生物與環(huán)境學(xué)院, 南京 210037)

    0 引言

    林業(yè)生產(chǎn)中,森林資源調(diào)查與數(shù)據(jù)采集是一個重要環(huán)節(jié)。樹木胸徑是立木測定的最基本因子之一。測量樹木胸徑可以了解樹木生長狀況,為樹木可視化以及整個林區(qū)的生產(chǎn)控制提供重要數(shù)據(jù)。早期測量胸徑采用卡尺或圍尺人工測量[1],這種測量方式勞動強度大,特別是對于有雜草的森林。隨著電子技術(shù)的不斷發(fā)展,各種測量樹干胸徑方法相繼產(chǎn)生,實現(xiàn)無接觸的樹干胸徑測量[2-16]。有學(xué)者采用數(shù)字圖像方法獲取樹干胸徑,如利用條形碼閱讀器掃描并讀取樹干上帶有樹木胸徑的條形碼帶的數(shù)值[2],獲得樹干胸徑。王雪峰等[3]采用原野服務(wù)器傳回的圖像,并通過立體的視覺技術(shù)重建了遠距離樹木的胸徑。王建利等[4]采用光學(xué)三角形測量理論,搭建一個由線性激光傳感器與CMOS圖像傳感器組成的測量系統(tǒng),通過圖像數(shù)據(jù)與樹干坐標信息獲取樹干表面輪廓圖像,最后通過最小二乘擬合橢圓法求取樹干胸徑。除了直接根據(jù)圖像獲取樹干直徑以外,F(xiàn)ORSMAN等[5]采取多視覺圖像重構(gòu)樹干,獲得樹干點云數(shù)據(jù),采用Gauss-Newton法擬合圓,求取樹干胸徑。黃曉東等[6]采用CCD、傾角傳感器、距離傳感器等設(shè)計了一套自動測量胸徑與樹高的測量儀。近幾年,較多學(xué)者采用地面激光點云數(shù)據(jù)分析樹干參數(shù)[7-11],也有較多學(xué)者采用二維激光數(shù)據(jù)測量樹干直徑[12-15]。也有學(xué)者融合激光傳感器與機器視覺測量單木樹胸徑[16]。機器視覺比較符合人的思維,但是在自然環(huán)境下,容易受到光照等條件的影響。地面激光雷達技術(shù)點云數(shù)據(jù)處理非常費時,二維激光技術(shù)特別是胸徑位置處采用二維激光技術(shù)測量速度較快,且輕便,容易制成測量儀。

    在解決優(yōu)化問題中,差分進化算法(Differential evolution, DE)是一種非常有效的優(yōu)化算法[17]。它是由STORN等[18]于1995年提出的。為解決差分進化算法對參數(shù)敏感的問題,許多研究者提出了參數(shù)自適應(yīng)的差分進化算法。QIN等[19]提出了自適應(yīng)差分進化算法(SaDE)。LIU[20]提出參數(shù)自適應(yīng)差分算法(JADE)優(yōu)化JIT(Just-in-time)與相關(guān)向量機(Relevant vector machine,RVM)。PEUURI等[21]針對DE/rand/1/bin策略算子提出一種自適種群數(shù)的DE算法,找出合適的交叉概率與對應(yīng)的種群數(shù),變異尺度在[0,1)內(nèi)隨機取,該算法對較復(fù)雜的函數(shù)優(yōu)化效果較好。本文提出一種新的非接觸式測量樹干胸徑方法,采用實驗室自制的基于激光傳感器的履帶移動機器人胸徑測量平臺測量樹干胸徑。針對樹干近似圓形,將傳感器采集的樹干表面參數(shù)用圓進行擬合,擬合出圓的半徑,從而獲得樹干胸徑。針對本研究數(shù)據(jù),在PEUURI等[21]研究基礎(chǔ)上,加入自適應(yīng)策略與自適應(yīng)控制參數(shù),實現(xiàn)圓的參數(shù)擬合。

    1 實驗材料

    本文實驗數(shù)據(jù)采用實驗室自制胸徑測量儀測量獲得,其測量平臺示意圖如圖1所示,以Arduino為微處理器設(shè)計測量控制系統(tǒng)。采用長、寬、高為540、280、110 mm 的履帶小車為載體,將控制器以及傳感器固定在智能小車平臺上。平臺上固定一伸縮式鋁合金撐桿,底部用螺絲螺母連接于小車底盤上,高度可以調(diào)整。撐桿頂端帶有三維多角度旋轉(zhuǎn)云臺,可以滿足多角度的測量需求。激光傳感器(VL53L0X型)與舵機用L型連接件連接,采用透明亞克力塑料板作為激光傳感器的固定板,用螺絲固定在L型金屬連接件上。測量平臺采用紅外接收模塊與紅外遙控器控制小車前進、后退、左右轉(zhuǎn)、停止以及采集數(shù)據(jù)等。測量數(shù)據(jù)存放于控制系統(tǒng)的SD卡內(nèi)。采用超聲波傳感器測量樹干表面到測量平臺的距離。

    圖1 測量平臺示意圖Fig.1 Sketch of sample collection platform1.激光傳感器 2.舵機 3.控制線 4.可調(diào)整鋁合金撐桿 5.超聲波傳感器 6.履帶移動機器人平臺 7.電動機 8.控制系統(tǒng) 9.紅外接收模塊 10.紅外遙控器

    測量平臺可以通過遙控器實現(xiàn)遠程控制,控制測量平臺運動與數(shù)據(jù)采集。激光傳感器測量范圍在2 m以內(nèi)。本實驗平臺測量范圍為50~1 000 mm,測量精度為±1.0 mm,舵機轉(zhuǎn)角范圍在0°~180°,環(huán)境工作溫度為-10~45℃。

    2 自適應(yīng)差分進化算法擬合圓

    2.1 擬合圓方程的建立

    由于大部分樹干呈圓臺形,對于同一高度的樹干截面呈圓形,故根據(jù)同一平面上任意3個點能確定圓。為了使擬合圓更精確,每次舵機帶動激光傳感器從0°到180°采集數(shù)據(jù),對同一個方位采集3次。如圖2所示,每次采集數(shù)據(jù)時傳感器自動調(diào)整到零初始位置。初始位置設(shè)置為極坐標起始坐標位置,傳感器中心為坐標原點,每隔15 ms激光傳感器轉(zhuǎn)動1°,采集一次數(shù)據(jù),數(shù)據(jù)單位均為mm。將采集到的數(shù)據(jù)利用K-means算法聚類,保留目標數(shù)據(jù)點。如圖2所示,設(shè)轉(zhuǎn)動角度與初始位置之間的夾角為α。設(shè)擬合圓心坐標為(xc,yc),半徑為R,采集樣本點t時刻的坐標(xt,yt)到圓心距離為dt,其中xt=ρtcosα,yt=ρtsinα,ρt為傳感器到樹干表面的距離,為使樣本點盡可能擬合在圓心為(xc,yc)、半徑為R的圓上,則所有樣本點到圓心距離與半徑差的和應(yīng)最小,即優(yōu)化目標方程為

    (1)

    (2)

    (3)

    式中N——目標樣本點數(shù)θ——優(yōu)化參數(shù)

    圖2 樣本采集原理示意圖Fig.2 Sample collection principal diagram

    圓心坐標與半徑均為未知數(shù),即優(yōu)化參數(shù)θ=(xc,yc,R)。需要擬合優(yōu)化算法對其進行優(yōu)化,基于梯度的Levenberg-Marquardt(LM)優(yōu)化算法、進化優(yōu)化算法、粒子群優(yōu)化算法等均有各自的優(yōu)缺點。本文采用差分進化算法擬合圓。

    2.2 自適應(yīng)差分進化算法

    差分進化算法是基于群體智能理論的優(yōu)化算法,通過群體內(nèi)個體間的合作與競爭產(chǎn)生的群體智能指導(dǎo)優(yōu)化搜索,主要有變異、交叉與選擇3個過程,而算法過程中的種群初始化、種群規(guī)模、進化策略、變異尺度以及交叉概率等均影響優(yōu)化算法的收斂性與魯棒性。因此本文主要對種群規(guī)模、進化策略、進化算子以及交叉概率在進化過程中自適應(yīng)調(diào)整進行研究,其算法過程如下:

    (1)初始化

    設(shè)置初始化種群與代數(shù)以及計算過程中的其他參數(shù)初始值。假設(shè)用G=0,1,…,Gmax表示進化代數(shù),則當前代數(shù)下種群中第i個個體表示為

    Xi,G=(xi,G,1,xi,G,2,…,xi,G,D)
    (i=1,2,…,Np)

    (4)

    式中D——種群個體維數(shù)

    Np——種群個數(shù)

    在種群初始化時,初始的種群需覆蓋整個搜索空間RD,初始化為

    xi,0,j=xL,j+r(xU,j-xL,j)
    (j=1,2,…,D)

    (5)

    式中xU,j——第j維的上界

    xL,j——第j維的下界

    r——[0,1]區(qū)間內(nèi)均勻分布的隨機數(shù)

    種群個數(shù)必須大于進化向量數(shù),目前較多的學(xué)者對種群個數(shù)Np進行研究,種群個數(shù)與種群向量幅度有關(guān),本研究采用PEUURI等[21]提出的種群自適應(yīng)算法

    Np=kDln(xU-xL+e)

    (6)

    式中k——比例系數(shù)

    xL——下界xU——上界

    本文D為3,初始種群數(shù)Np0=30,k由PEUURI 等[21]提出的自適應(yīng)算法獲取。

    (2)差分變異操作

    在當前代數(shù)中,對于每一個個體矢量Xi,G(稱為目標矢量),使用變異算子生成一個新的個體Vi,G(稱為變異矢量)。目前變異策略有多種形式,經(jīng)過多次實驗可知DE/rand/1/bin策略與DE/current-to-best/1策略對圓擬合效果較好,從當前種群中隨機挑選3個個體矢量Xr1,G、Xr2,G與Xr3,G,并選出最優(yōu)Xbest,G,則變異策略為[22]

    (7)

    (8)

    式中Fi——變異尺度因子

    ζ——選擇概率

    λ——自適應(yīng)因子,取0.5

    γ——常數(shù),取0.01

    對差分矢量進行縮放,從而控制搜索步長。較大Fi對全局搜索比較有效,較小的Fi可以加速收斂??梢?,選擇合適的Fi非常重要,一般選擇Fi∈(0.4,0.9),然而對不同的進化代數(shù)Fi不是一個定值,F(xiàn)i隨著進化過程變化而動態(tài)變化。本文選用文獻[23]提出的按照Cauchy分布,其參數(shù)隨進化代數(shù)變化而變化。

    (3)交叉操作

    差分進化算法采用離散交叉算子:二項式交叉和指數(shù)交叉。交叉算子把通過變異操作得到的變異矢量Vi,G與個體目標矢量Xi,G進行離散交叉生成試驗矢量Ui,G,具體的二項式交叉操作為

    (9)

    式中Cri——第i次進化交叉概率因子

    jrand——隨機整數(shù)

    一般Cri∈(0,1),但在Cri∈(0.1,0.8)時更有效。在j維上,若隨機生成數(shù)小于Cri時,試驗矢量繼承變異矢量,反之繼承父代目標矢量。一般而言Cri越大則選擇的交叉位也就越多,意味著試驗矢量更傾向于變異矢量。本文Cri隨著進化過程正態(tài)分布,其正態(tài)分布均值隨進化代數(shù)變化而變化[23]。

    (4)選擇操作

    差分進化算法通過變異操作和交叉操作產(chǎn)生后代群體之后,采用一對一的貪婪篩選算子將子個體與相應(yīng)的父個體進行比較,較優(yōu)者保存到下一代,選擇操作為

    (10)

    (5)優(yōu)化結(jié)束條件判斷

    若不滿足最大迭代次數(shù)1 000,繼續(xù)進化過程,采用式(6)更新種群數(shù),進入變異、交叉與選擇過程。否則結(jié)束進化優(yōu)化過程,輸出最優(yōu)值。

    3 實驗與分析

    圖3 樣本采集場景Fig.3 Sample collection scene

    在林分測定中,胸徑(DBH)指樹干距離地面以上1.3 m處樹干的直徑。本次測量平臺為履帶機器人平臺,樣本采集場景如圖3所示。為驗證機器人行走穩(wěn)定性,在校園與人工林地隨機選擇大小不同的40棵樹進行測量,對每一棵樹進行標記。為驗證本文測量系統(tǒng)測量胸徑的準確性,首先對標記好的樹采用圍尺人工測量出胸徑值作為樹干胸徑實際值,其數(shù)據(jù)統(tǒng)計分析表如表1所示。

    表1 樣本胸徑統(tǒng)計Tab.1 DBH statistics of measured sample mm

    為驗證系統(tǒng)測量的穩(wěn)定性,分別在陰天的早晨、中午、下午與晴天的早晨、中午、下午進行了6次數(shù)據(jù)測量,每次對每個場景(校園與人工林地)有標記的40棵樹以相同順序進行測量記錄數(shù)據(jù),每次測量3個方位,求取平均值作為測量值。實驗數(shù)據(jù)處理器為Intel i5-6500 CPU 3.2 GHz,內(nèi)存8.0 GB,數(shù)據(jù)處理軟件為Matlab 2016b。

    激光傳感器采集的數(shù)據(jù)由樹干表面反射到傳感器的時間與速度計算出來,由于激光傳感器存在一定的發(fā)散性,本研究首先將采集的數(shù)據(jù)經(jīng)過發(fā)散角校正,然后計算樹干表面點的坐標值。分別將不同時間段與不同方位采集的數(shù)據(jù)采用K-means算法分割出目標數(shù)據(jù)。采用文獻[14]提出的弧長法矯正激光傳感器的發(fā)散角,擬合出發(fā)散角二次方程,每次計算時采用角度補償?shù)姆绞綔p少誤差,擬合R2為0.907,均方根誤差(RMSE)為0.626°,其方程為

    Fρ=-0.000 592ρ2+0.064 28ρ-11.1

    (11)

    式中ρ——激光傳感器中心到樹干表面點距離,mm

    Fρ——補償角,(°)

    補償角隨著距離的增加而增大,因此每次計算樹干表面點坐標值之前,將對應(yīng)的角度減去補償角即為補償后角度。

    圖4 二維坐標點擬合圓Fig.4 Two dimensional coordinate point fitting circle

    計算出補償后數(shù)據(jù)對應(yīng)的坐標值,采用2.2節(jié)的自適應(yīng)差分進化算法擬合圓。圖4是經(jīng)過處理后的數(shù)據(jù)點擬合圓效果圖,顯示大部分數(shù)據(jù)點比較集中在圓弧上,極少部分數(shù)據(jù)有一定的偏差,由于樹干表面不是絕對的光滑,因此表面的部分數(shù)據(jù)不在同一圓弧上。但是對于同一個方向的數(shù)據(jù)點,大部分在同一個圓弧上。分別對3個方向的數(shù)據(jù)進行分割與擬合,擬合出最優(yōu)半徑,計算3個方向的平均值即為本文所求樹干胸徑。

    本次實驗分別對基于梯度下降法Levenberg-Marquardt(LM)優(yōu)化算法、非梯度算法粒子群方法(PSO)以及本文提出的自適應(yīng)差分進化算法(ADE)進行了比較分析。對同一組數(shù)據(jù)進行比較,其結(jié)果如表2所示,基于梯度下降方法的LM算法收斂較慢,運行時間較長,PSO算法與ADE算法收斂速度快,而ADE算法最快。擬合圓的直徑較相近,因此均方根誤差RMSE相差較小。從實驗結(jié)果可以看出,ADE算法較優(yōu)。

    表2 不同算法擬合結(jié)果比較Tab.2 Comparison of different algorithms

    將Hough變換算法與ADE算法進行比較。對同一棵樹的6次測量數(shù)據(jù)進行了統(tǒng)計分析,測量誤差在±5 mm之間,求取平均值作為測量值。將實驗測量值分別與實際值進行對比分析,圖5與圖6分別給出了校園內(nèi)與人工林地不同直徑的40棵樹胸徑測量值與實際值之間的關(guān)系。由圖5與圖6可知,本文算法數(shù)據(jù)點較集中。表3給出了采用本文算法與Hough變換算法擬合的測量值與實際值的直線擬合結(jié)果。本文算法對校園內(nèi)數(shù)據(jù)擬合直線斜率為0.993,偏離值為3.782,R2為0.989,RMSE為4.996 mm。而Hough變換算法擬合的直線斜率為1.038,偏離值為-7.981,R2為0.947,RMSE為12.18 mm。對于人工林地,本文擬合算法的RMSE為4.500 mm,Hough變換算法的RMSE為11.97 mm。由實驗分析可知,本文測量胸徑方法受樣地場景影響較小。雖然本文算法擬合圓效果較好,但是實際值與測量值之間存在一定的誤差。引起誤差的因素較多,如樹干實際測量的誤差、傳感器誤差以及傳感器發(fā)散角補償誤差等。在以后的研究工作中,應(yīng)盡量減少測量誤差,提出更有效的發(fā)散角補償算法,使測量值更逼近于真實值。

    4 結(jié)論

    (1)采用實驗室自制的基于激光傳感器的胸徑測量平臺對校園內(nèi)與人工林地不同直徑的40棵樹進行多次測量,每次測量對同一棵樹采集3個方向的數(shù)據(jù),每個方向采集3次,獲得不同方位的測量數(shù)據(jù),并將數(shù)據(jù)保存于SD卡。

    (2)對采集的數(shù)據(jù)采用K-means進行聚類,分

    圖5 校園內(nèi)樹干胸徑測量值與實際值關(guān)系Fig.5 Relationship between actual DBH and measured DBH in campus

    圖6 人工林地樹干胸徑測量值與實際值關(guān)系Fig.6 Relationship between actual DBH and measured DBH in artificial forest

    擬合方法擬合方程RMSE/mmR2ADE算法校園內(nèi)y=0.993x+3.7824.9960.989人工林地y=0.979x+4.9754.5000.993Hough變換算法校園內(nèi)y=1.038x-7.98112.180.947人工林地y=0.973x+8.51211.970.955

    割出背景與目標。采用弧長法計算出補償角與激光傳感器到樹干表面點距離的擬合公式。將采集到的數(shù)據(jù)角度分別減去補償角度,對分割的目標樣本數(shù)據(jù)進行補償。對補償后的數(shù)據(jù),以激光傳感器中心為坐標原點,計算出樹干表面點的橫縱坐標。

    (3)根據(jù)每棵樹的橫縱坐標,分別采用基于梯度的LM優(yōu)化算法、粒子群(PSO)優(yōu)化算法、本文提出的自適應(yīng)差分進化算法以及Hough變換算法擬合圓測量出樹干胸徑,并對測量值與實際值進行了分析與比較。本文算法擬合時間最短,為1.41 s,對校園內(nèi)和人工林測量結(jié)果與實際值比較分析,RMSE分別為4.996 mm與4.500 mm。實驗結(jié)果表明,本文提出的自適應(yīng)差分進化算法收斂快,且比Hough變換擬合圓方法更逼近真實值。

    猜你喜歡
    胸徑樹干差分
    白城山新1號楊育苗密度研究
    武漢5種常見園林綠化樹種胸徑與樹高的相關(guān)性研究
    為什么樹干不是方的?
    軍事文摘(2021年16期)2021-11-05 08:49:16
    數(shù)列與差分
    福建省森林資源監(jiān)測體系抽樣調(diào)查中胸徑測量精度范圍的精準確定
    為什么要在樹干上刷白漿
    為什么要在樹干上刷白漿
    為什么要在樹干上刷一層白漿
    基于差分隱私的大數(shù)據(jù)隱私保護
    相對差分單項測距△DOR
    太空探索(2014年1期)2014-07-10 13:41:50
    夜夜看夜夜爽夜夜摸| 亚洲人成网站在线观看播放| 白带黄色成豆腐渣| 此物有八面人人有两片| 日韩制服骚丝袜av| 日本爱情动作片www.在线观看| 午夜激情欧美在线| 国产 一区 欧美 日韩| 欧美变态另类bdsm刘玥| 两个人的视频大全免费| 亚洲av男天堂| 性插视频无遮挡在线免费观看| 一级av片app| 成人一区二区视频在线观看| 亚洲国产高清在线一区二区三| 久久亚洲精品不卡| 美女xxoo啪啪120秒动态图| 亚洲高清免费不卡视频| 99riav亚洲国产免费| 久久午夜福利片| 一边摸一边抽搐一进一小说| 国产精华一区二区三区| 久久久国产成人免费| 简卡轻食公司| 国产91av在线免费观看| 午夜精品一区二区三区免费看| 真实男女啪啪啪动态图| 国产成人a区在线观看| 国产成人a区在线观看| 欧美潮喷喷水| 一级毛片aaaaaa免费看小| 久久久欧美国产精品| 人妻制服诱惑在线中文字幕| 国内精品一区二区在线观看| 女同久久另类99精品国产91| 日韩欧美 国产精品| 最好的美女福利视频网| 久久亚洲精品不卡| 国产一区二区在线av高清观看| 免费大片18禁| 欧美xxxx黑人xx丫x性爽| 乱系列少妇在线播放| 精品99又大又爽又粗少妇毛片| 亚洲,欧美,日韩| 午夜福利在线在线| 亚洲无线在线观看| 听说在线观看完整版免费高清| 老司机影院成人| 国产高清三级在线| 寂寞人妻少妇视频99o| 美女黄网站色视频| 国产成人影院久久av| 99久久人妻综合| 久久久久久久久久久免费av| 久久精品久久久久久噜噜老黄 | 毛片女人毛片| 精品久久久久久久久久免费视频| 精品熟女少妇av免费看| 免费不卡的大黄色大毛片视频在线观看 | 一进一出抽搐动态| 日韩一区二区视频免费看| 国产免费一级a男人的天堂| 国产精品国产三级国产av玫瑰| 欧美在线一区亚洲| 99热这里只有精品一区| 网址你懂的国产日韩在线| 久久这里只有精品中国| 成人无遮挡网站| 国内久久婷婷六月综合欲色啪| 久久久久久久久久久免费av| 国产精品一及| 99riav亚洲国产免费| av在线老鸭窝| 三级国产精品欧美在线观看| 免费不卡的大黄色大毛片视频在线观看 | 免费av毛片视频| 久久午夜福利片| 中文在线观看免费www的网站| 女同久久另类99精品国产91| 亚洲av免费在线观看| 毛片女人毛片| 国产av麻豆久久久久久久| 一进一出抽搐gif免费好疼| 亚洲精品粉嫩美女一区| 99在线视频只有这里精品首页| 亚洲三级黄色毛片| 午夜福利高清视频| 久久久精品欧美日韩精品| 亚洲一区高清亚洲精品| 国产精品乱码一区二三区的特点| 亚洲欧洲国产日韩| 嫩草影院入口| 一本精品99久久精品77| 国产极品精品免费视频能看的| 国产成年人精品一区二区| 午夜精品在线福利| 国产激情偷乱视频一区二区| 亚洲真实伦在线观看| 美女xxoo啪啪120秒动态图| 99久国产av精品| 九草在线视频观看| 欧美日韩精品成人综合77777| 少妇的逼好多水| 成人无遮挡网站| 毛片女人毛片| 成人特级黄色片久久久久久久| 日本熟妇午夜| 人体艺术视频欧美日本| 色噜噜av男人的天堂激情| 青青草视频在线视频观看| 国产黄色视频一区二区在线观看 | 少妇人妻精品综合一区二区 | 精品久久久久久久久久免费视频| 亚洲精华国产精华液的使用体验 | 毛片一级片免费看久久久久| 免费观看的影片在线观看| 女人被狂操c到高潮| 看片在线看免费视频| 亚洲欧美日韩东京热| 国产成人精品久久久久久| 22中文网久久字幕| 在线观看一区二区三区| 五月玫瑰六月丁香| 女人十人毛片免费观看3o分钟| 国产成人aa在线观看| 欧美3d第一页| 小蜜桃在线观看免费完整版高清| 91狼人影院| 国语自产精品视频在线第100页| 最好的美女福利视频网| 国产亚洲91精品色在线| 69av精品久久久久久| 亚洲欧美日韩高清专用| 搡老妇女老女人老熟妇| 好男人视频免费观看在线| 久久久国产成人免费| 晚上一个人看的免费电影| 一级毛片aaaaaa免费看小| 免费一级毛片在线播放高清视频| 三级男女做爰猛烈吃奶摸视频| av在线蜜桃| 国产成年人精品一区二区| 日韩av在线大香蕉| 久99久视频精品免费| 国产精品伦人一区二区| 成人午夜精彩视频在线观看| 亚洲色图av天堂| 久久久久久九九精品二区国产| 国产不卡一卡二| 麻豆国产97在线/欧美| 亚洲av第一区精品v没综合| 中国美女看黄片| 成人性生交大片免费视频hd| 又爽又黄无遮挡网站| 内地一区二区视频在线| 成人美女网站在线观看视频| 伦理电影大哥的女人| 天天躁夜夜躁狠狠久久av| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品久久久久久久久免| 国产亚洲av片在线观看秒播厂 | 国产高清不卡午夜福利| 日韩人妻高清精品专区| 国产成年人精品一区二区| 色播亚洲综合网| 日韩 亚洲 欧美在线| 亚洲美女搞黄在线观看| 国产久久久一区二区三区| 国产成人精品一,二区 | 成年女人永久免费观看视频| 亚洲七黄色美女视频| 免费人成视频x8x8入口观看| 国产片特级美女逼逼视频| 级片在线观看| 搞女人的毛片| 成年女人看的毛片在线观看| 自拍偷自拍亚洲精品老妇| 九色成人免费人妻av| 国产老妇伦熟女老妇高清| 亚洲久久久久久中文字幕| 赤兔流量卡办理| 两个人的视频大全免费| 看非洲黑人一级黄片| 久久久久久国产a免费观看| 久久久久九九精品影院| 欧美丝袜亚洲另类| 国产伦精品一区二区三区四那| 久久久久久九九精品二区国产| 久久久久久久久大av| 深爱激情五月婷婷| 99久久精品国产国产毛片| 久久午夜亚洲精品久久| 91精品国产九色| 亚洲人成网站在线观看播放| 欧美日韩一区二区视频在线观看视频在线 | 搡老妇女老女人老熟妇| 久久人人爽人人片av| 亚洲在线自拍视频| 国产精品乱码一区二三区的特点| 亚洲人成网站在线播放欧美日韩| 18禁黄网站禁片免费观看直播| 97在线视频观看| 好男人在线观看高清免费视频| 亚洲精品456在线播放app| 人人妻人人看人人澡| 亚洲,欧美,日韩| 日本成人三级电影网站| 久久99热这里只有精品18| 亚洲精品影视一区二区三区av| 亚洲第一电影网av| 精品不卡国产一区二区三区| 欧美bdsm另类| 嫩草影院精品99| 亚洲五月天丁香| 久久人妻av系列| 大香蕉久久网| 五月伊人婷婷丁香| 亚洲av第一区精品v没综合| 中出人妻视频一区二区| 日本成人三级电影网站| 婷婷六月久久综合丁香| 最近中文字幕高清免费大全6| 国产美女午夜福利| 国产91av在线免费观看| 少妇高潮的动态图| 少妇人妻一区二区三区视频| а√天堂www在线а√下载| 成年免费大片在线观看| 国产午夜精品论理片| 又粗又爽又猛毛片免费看| a级毛片a级免费在线| 久久草成人影院| 久久精品91蜜桃| 婷婷色综合大香蕉| 一级二级三级毛片免费看| 欧美丝袜亚洲另类| 欧美日韩国产亚洲二区| 18禁黄网站禁片免费观看直播| 插阴视频在线观看视频| 蜜桃久久精品国产亚洲av| 国产伦在线观看视频一区| 国产三级中文精品| 91久久精品国产一区二区三区| 久久午夜亚洲精品久久| 午夜免费男女啪啪视频观看| 国产极品天堂在线| 综合色av麻豆| 国产片特级美女逼逼视频| 国产精品麻豆人妻色哟哟久久 | 久99久视频精品免费| 午夜视频国产福利| 蜜臀久久99精品久久宅男| 久久久久久大精品| 欧美3d第一页| 变态另类成人亚洲欧美熟女| 菩萨蛮人人尽说江南好唐韦庄 | 日本av手机在线免费观看| 久久精品国产亚洲网站| 国产午夜精品论理片| 国产真实伦视频高清在线观看| 观看美女的网站| 国产亚洲欧美98| 国产伦一二天堂av在线观看| 哪个播放器可以免费观看大片| 国产精品女同一区二区软件| 免费在线观看成人毛片| 国产又黄又爽又无遮挡在线| 一个人看视频在线观看www免费| 日韩中字成人| 丝袜美腿在线中文| 亚洲图色成人| 寂寞人妻少妇视频99o| 少妇高潮的动态图| 99国产极品粉嫩在线观看| 亚洲熟妇中文字幕五十中出| 一级毛片电影观看 | 自拍偷自拍亚洲精品老妇| 色综合站精品国产| 国产亚洲精品久久久com| 国产成人a区在线观看| 国产一级毛片七仙女欲春2| 一边亲一边摸免费视频| 欧美色视频一区免费| av在线老鸭窝| 国产免费一级a男人的天堂| 男女边吃奶边做爰视频| 桃色一区二区三区在线观看| 国产成人freesex在线| 国产高清激情床上av| videossex国产| 久久久a久久爽久久v久久| 边亲边吃奶的免费视频| av在线天堂中文字幕| 欧美成人精品欧美一级黄| .国产精品久久| 国产成人精品婷婷| 熟女人妻精品中文字幕| 一个人看视频在线观看www免费| 草草在线视频免费看| 亚洲av免费高清在线观看| 亚洲国产日韩欧美精品在线观看| 少妇被粗大猛烈的视频| 久久久久九九精品影院| 欧美一级a爱片免费观看看| 哪个播放器可以免费观看大片| 非洲黑人性xxxx精品又粗又长| 亚洲一区二区三区色噜噜| av视频在线观看入口| 尾随美女入室| 九九热线精品视视频播放| 亚洲在线自拍视频| 久久久久久伊人网av| 99热这里只有精品一区| 在线观看av片永久免费下载| 久久国内精品自在自线图片| 九九在线视频观看精品| 国模一区二区三区四区视频| 日韩,欧美,国产一区二区三区 | 欧美丝袜亚洲另类| 国产精品久久久久久亚洲av鲁大| 亚洲国产精品成人久久小说 | a级毛色黄片| 欧美精品一区二区大全| 精品欧美国产一区二区三| 女同久久另类99精品国产91| 国产精品久久视频播放| 国产黄片美女视频| 欧美zozozo另类| 亚洲在久久综合| 中国国产av一级| 欧美变态另类bdsm刘玥| 成人特级av手机在线观看| 国产探花在线观看一区二区| av又黄又爽大尺度在线免费看 | 免费av观看视频| 亚洲国产欧美人成| 99久久精品一区二区三区| 久久精品影院6| 精品久久久久久久久久免费视频| 亚洲自偷自拍三级| 天堂网av新在线| 色吧在线观看| 国产亚洲av片在线观看秒播厂 | 免费观看在线日韩| 高清毛片免费观看视频网站| 在线观看66精品国产| 欧美激情在线99| 久久精品国产亚洲av涩爱 | 欧美日韩综合久久久久久| 夜夜爽天天搞| av福利片在线观看| 日韩欧美精品免费久久| 在线a可以看的网站| 毛片一级片免费看久久久久| 国产精品一区二区在线观看99 | 男女做爰动态图高潮gif福利片| 欧美另类亚洲清纯唯美| 一区福利在线观看| 亚洲久久久久久中文字幕| 色噜噜av男人的天堂激情| 女的被弄到高潮叫床怎么办| 国产不卡一卡二| 午夜福利成人在线免费观看| 热99在线观看视频| 精品久久久久久久久久久久久| 久久亚洲精品不卡| 高清在线视频一区二区三区 | 亚洲精品国产成人久久av| 欧美潮喷喷水| 淫秽高清视频在线观看| 欧美高清成人免费视频www| 欧美一区二区国产精品久久精品| 欧美变态另类bdsm刘玥| 欧美性猛交╳xxx乱大交人| 亚洲国产日韩欧美精品在线观看| 在线观看美女被高潮喷水网站| 伦精品一区二区三区| 中国国产av一级| 国产伦理片在线播放av一区 | 一级毛片aaaaaa免费看小| 床上黄色一级片| 国产伦精品一区二区三区四那| 国产精品三级大全| 日韩在线高清观看一区二区三区| 我要搜黄色片| 国产精品一区二区三区四区免费观看| 亚洲婷婷狠狠爱综合网| 小说图片视频综合网站| 男插女下体视频免费在线播放| 美女脱内裤让男人舔精品视频 | 看十八女毛片水多多多| 一级毛片我不卡| 国产熟女欧美一区二区| 在现免费观看毛片| 久久热精品热| 国产亚洲欧美98| 美女大奶头视频| 人妻久久中文字幕网| 久久精品久久久久久噜噜老黄 | 在线观看免费视频日本深夜| 亚洲av免费高清在线观看| 1000部很黄的大片| 一级黄片播放器| 欧美又色又爽又黄视频| 免费看av在线观看网站| 欧美bdsm另类| 国产亚洲av片在线观看秒播厂 | 久久热精品热| 麻豆av噜噜一区二区三区| 久久99蜜桃精品久久| 久久人人爽人人片av| 亚洲自拍偷在线| 极品教师在线视频| 婷婷六月久久综合丁香| 久久久精品大字幕| 国产视频首页在线观看| 内射极品少妇av片p| 一个人免费在线观看电影| 色综合色国产| 午夜福利视频1000在线观看| 1000部很黄的大片| 99热这里只有是精品在线观看| 国产 一区精品| 午夜精品国产一区二区电影 | 午夜福利成人在线免费观看| 高清日韩中文字幕在线| 99热精品在线国产| 深爱激情五月婷婷| 婷婷精品国产亚洲av| 久久国内精品自在自线图片| 此物有八面人人有两片| 人妻久久中文字幕网| 丰满的人妻完整版| 大又大粗又爽又黄少妇毛片口| 麻豆国产97在线/欧美| 在线播放无遮挡| 亚洲18禁久久av| 精品久久久久久久久久免费视频| 亚洲精品国产成人久久av| 特级一级黄色大片| 亚洲精品久久国产高清桃花| 国内精品美女久久久久久| 熟女人妻精品中文字幕| 97人妻精品一区二区三区麻豆| 日韩欧美精品免费久久| 中文字幕久久专区| 舔av片在线| 亚洲国产精品合色在线| 高清毛片免费观看视频网站| 日韩人妻高清精品专区| 如何舔出高潮| 国产私拍福利视频在线观看| 一级毛片aaaaaa免费看小| 欧美性感艳星| .国产精品久久| 国产老妇伦熟女老妇高清| av天堂在线播放| 欧美另类亚洲清纯唯美| 亚洲国产欧洲综合997久久,| av在线亚洲专区| 免费搜索国产男女视频| 身体一侧抽搐| 亚洲成a人片在线一区二区| 色噜噜av男人的天堂激情| 久久久a久久爽久久v久久| 男女边吃奶边做爰视频| 看免费成人av毛片| 日本欧美国产在线视频| 村上凉子中文字幕在线| 中文字幕人妻熟人妻熟丝袜美| 精品久久国产蜜桃| 久久久久国产网址| 男人和女人高潮做爰伦理| 日产精品乱码卡一卡2卡三| 欧美成人精品欧美一级黄| 99久久中文字幕三级久久日本| 精品国内亚洲2022精品成人| 欧美性感艳星| 日本成人三级电影网站| 最近的中文字幕免费完整| 深夜a级毛片| 午夜久久久久精精品| 变态另类丝袜制服| 成人高潮视频无遮挡免费网站| 日本一二三区视频观看| 人人妻人人澡欧美一区二区| 别揉我奶头 嗯啊视频| 波多野结衣高清作品| 色5月婷婷丁香| 午夜福利在线观看吧| 22中文网久久字幕| 欧美另类亚洲清纯唯美| 99在线人妻在线中文字幕| 国产伦理片在线播放av一区 | 天天一区二区日本电影三级| 国产亚洲精品久久久com| 国产亚洲5aaaaa淫片| 色尼玛亚洲综合影院| 色哟哟·www| 久久精品国产亚洲av涩爱 | 人人妻人人澡欧美一区二区| av黄色大香蕉| 欧美性猛交黑人性爽| 亚洲色图av天堂| 亚洲欧美中文字幕日韩二区| 国产一区二区在线观看日韩| 免费观看精品视频网站| а√天堂www在线а√下载| 天堂中文最新版在线下载 | 男人舔奶头视频| 欧美高清性xxxxhd video| 床上黄色一级片| 亚洲国产精品合色在线| av女优亚洲男人天堂| 亚洲精品成人久久久久久| av在线老鸭窝| 观看美女的网站| 久久精品91蜜桃| 国产69精品久久久久777片| 91麻豆精品激情在线观看国产| 久久久欧美国产精品| 好男人视频免费观看在线| 91午夜精品亚洲一区二区三区| 成年av动漫网址| 国产一区二区激情短视频| 99久久九九国产精品国产免费| 精品欧美国产一区二区三| 哪个播放器可以免费观看大片| 亚洲在线观看片| 一夜夜www| 午夜福利在线观看免费完整高清在 | 欧美日韩综合久久久久久| 秋霞在线观看毛片| 日日撸夜夜添| av在线亚洲专区| 美女cb高潮喷水在线观看| 国产av在哪里看| 国产大屁股一区二区在线视频| or卡值多少钱| 亚洲成人av在线免费| 99热全是精品| 国产精品国产三级国产av玫瑰| 欧美日韩国产亚洲二区| 国内精品宾馆在线| av免费观看日本| 精品国内亚洲2022精品成人| 91久久精品国产一区二区成人| 22中文网久久字幕| 国内久久婷婷六月综合欲色啪| or卡值多少钱| 亚洲国产高清在线一区二区三| 欧美丝袜亚洲另类| 深夜精品福利| 两性午夜刺激爽爽歪歪视频在线观看| 男人和女人高潮做爰伦理| 精品久久久久久久久亚洲| 亚洲丝袜综合中文字幕| 美女内射精品一级片tv| av在线观看视频网站免费| 十八禁国产超污无遮挡网站| 观看美女的网站| 国产免费一级a男人的天堂| 99热这里只有精品一区| 狂野欧美激情性xxxx在线观看| 熟女电影av网| 国产精品久久久久久久久免| 麻豆国产97在线/欧美| 免费av不卡在线播放| 国产视频内射| 国产精品人妻久久久影院| 岛国毛片在线播放| 国产中年淑女户外野战色| 亚洲av.av天堂| 国产免费男女视频| 成人性生交大片免费视频hd| 菩萨蛮人人尽说江南好唐韦庄 | 小说图片视频综合网站| 欧美丝袜亚洲另类| 长腿黑丝高跟| 国产黄片美女视频| 夜夜夜夜夜久久久久| 波多野结衣高清作品| 看十八女毛片水多多多| 国产视频内射| АⅤ资源中文在线天堂| 国语自产精品视频在线第100页| 国产日本99.免费观看| 国产精品人妻久久久影院| 99久久九九国产精品国产免费| 秋霞在线观看毛片| 国产精品人妻久久久影院| 国产精品.久久久| 国产日本99.免费观看| 国产男人的电影天堂91| 国产白丝娇喘喷水9色精品| 国产中年淑女户外野战色| 日韩三级伦理在线观看| 老熟妇乱子伦视频在线观看| a级毛色黄片| 久久精品影院6| 国产高清不卡午夜福利| 亚洲精品456在线播放app| 久久精品综合一区二区三区| 在线天堂最新版资源| 亚洲av电影不卡..在线观看| 狂野欧美激情性xxxx在线观看| 国产一级毛片七仙女欲春2| 高清毛片免费看| 少妇裸体淫交视频免费看高清| 午夜福利高清视频| 国产又黄又爽又无遮挡在线| 欧美zozozo另类| 2021天堂中文幕一二区在线观| 99热精品在线国产| 日韩高清综合在线|