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

    基于多傳感器信息融合的未知環(huán)境下移動機器人的地圖創(chuàng)建

    2021-06-07 06:10:30唐達明戴慶瑜
    陜西科技大學學報 2021年3期
    關(guān)鍵詞:深度融合優(yōu)化

    李 艷, 唐達明, 戴慶瑜

    (1.陜西科技大學 電氣與控制工程學院, 陜西 西安 710021; 2.陜西農(nóng)產(chǎn)品加工技術(shù)研究院, 陜西 西安 710021)

    0 引言

    未知環(huán)境中,機器人的地圖創(chuàng)建與自主定位密切相關(guān),SLAM(Simultaneous Localization and Mapping)將二者相結(jié)合,通過已經(jīng)獲得的環(huán)境地圖信息完成自主定位,并根據(jù)定位結(jié)果完成對環(huán)境地圖的更新[1].

    SLAM問題是指機器人在未知的環(huán)境下,僅依靠自身傳感器的觀測數(shù)據(jù)zt和控制輸入量ut的情況下確定機器人的大致位置,并且通過觀測信息構(gòu)建出大致的環(huán)境地圖,同時用構(gòu)建出來的地圖繼續(xù)對機器人的位置進行估計.如圖1所示,在機器人的移動過程中,通過一系列的控制輸入量u1∶t={u1,u2,…,ut}和觀測數(shù)據(jù)z1∶t={z1,z2,…,zt},對自身的運動軌跡x1∶t={x1,x2,…,xt}及環(huán)境地圖m進行準確的估計.

    圖1 SLAM問題示意圖

    SLAM問題用數(shù)學語言描述為求解如下概率分布:

    Bel(xt)=p(xt,m|u1∶t,z1∶t)

    或Bel(x0∶t)=p(x0∶t,m|u1∶t,z1∶t)

    (1)

    當前,求解SLAM問題的方法可大致分為兩大類:一類是基于濾波的方法,另一類則是基于圖優(yōu)化的方法.其中,基于濾波器的方法由于其算法簡單而被大量應用于SLAM中[2-4],但隨著研究的不斷深入,SLAM的應用范圍也越來越大,當使用基于濾波的方法時,機器人每次移動只考慮當前狀態(tài),在建圖時誤差會不斷累積,導致其最后構(gòu)建出的地圖與實際環(huán)境有差別,使其難以應用于較大環(huán)境下的地圖創(chuàng)建.而基于圖優(yōu)化的方法在估計機器人運動軌跡時,會用到傳感器的所有觀測信息,在建圖的過程中把對地圖特征的描述轉(zhuǎn)化為對運動軌跡的約束,提高了建圖的效率與準確性,使其適用于較大環(huán)境下的地圖創(chuàng)建.

    目前,對基于圖優(yōu)化的建圖算法研究越來越多,文獻[5]使用激光雷達完成地圖的構(gòu)建,將基于圖優(yōu)化的方法應用到地圖構(gòu)建中提高了建圖的效率;文獻[6]使用雙目相機完成地圖的構(gòu)建,將多層次松弛的方法應用在圖優(yōu)化的地圖構(gòu)建中,很好地解決了因算法形成閉環(huán)而無法更新地圖的問題;文獻[7]使用深度相機完成地圖的構(gòu)建,并在建圖時應用梯度下降方法,提高了SLAM的魯棒性.目前,對于SLAM的研究主要集中在對地圖構(gòu)建算法的研究[2-7],而對于傳感器數(shù)據(jù)處理方面的研究較少,并且以上方法都是基于單一的距離傳感器或視覺傳感器,都存在著很多缺陷,例如,二維激光雷達所獲取的信息會因為安裝位置不同而受到影響,視覺傳感器獲取的信息會因為環(huán)境的變化而受到影響,因此導致實驗效果仍然不佳.

    針對移動機器人在未知環(huán)境下利用單一傳感器進行地圖構(gòu)建時所出現(xiàn)的問題,本文提出利用基于圖優(yōu)化的算法構(gòu)建地圖,并在建圖時將移動機器人的激光與視覺信息進行融合,使得構(gòu)建出的地圖與實際環(huán)境更加接近.

    1 傳感器數(shù)據(jù)處理

    1.1 激光雷達數(shù)據(jù)處理

    激光雷達通過旋轉(zhuǎn)的方式掃描外界物體獲得信息[8],并以極坐標的形式表示為Zl=[dl,θl]T,其中,θl為激光雷達掃描時,目標點和激光雷達掃描起始位置之間的夾角,并且0 °<θl<360 °,dl是目標到激光雷達的直線距離,將原始掃描數(shù)據(jù)轉(zhuǎn)換至直角坐標,則表示形式為Zl=[x,y]T,如圖2所示.

    圖2 激光雷達掃描測量示意圖

    利用式(2),可將掃描得到的數(shù)據(jù)點轉(zhuǎn)換至直角坐標系中.

    (2)

    處理激光雷達數(shù)據(jù)時,可以將周圍環(huán)境的信息抽象的表示為點、線等特征[8].在掃描過的環(huán)境信息中可以得到線性與非線性的數(shù)據(jù),若是線性的數(shù)據(jù)則直接化為直線進行處理;若是非線性的數(shù)據(jù)則需要將這些數(shù)據(jù)變換為折線特征(即直線加角點),同時采用迭代的方法,對非線性區(qū)域中的角點進行尋找,直到掃描區(qū)域內(nèi)不存在角點,由公式(3)可以計算出角點之間的距離.這個過程不斷往復循環(huán),直到最后可以形成2D環(huán)境地圖.

    db=|(xi-x1)sinδ+(yi-y1)cosδ|

    (3)

    式(3)中:db為角點間的距離,(x1,y1)為機器人當前的位置坐標,(xi,yi)為掃描到的目標點的坐標,δ為角點之間形成的連線與X軸的夾角.

    1.2 深度相機數(shù)據(jù)處理

    深度相機的數(shù)據(jù)處理主要是將深度相機獲取到的圖像數(shù)據(jù)轉(zhuǎn)換為模擬激光數(shù)據(jù)[9].深度相機能夠通過相機模型將獲取到的深度圖像數(shù)據(jù)轉(zhuǎn)換成為模擬的二維激光雷達數(shù)據(jù),實現(xiàn)將三維數(shù)據(jù)壓縮為二維數(shù)據(jù)的功能,從而完成傳感器之間的數(shù)據(jù)融合.

    深度圖像數(shù)據(jù)轉(zhuǎn)換模擬激光數(shù)據(jù)的步驟如下:

    (1)計算圖像坐標系下的點m(u,v)在世界坐標系下的位置坐標M(XW,YW,ZW) .深度數(shù)據(jù)轉(zhuǎn)激光數(shù)據(jù)的原理如圖3所示.

    圖3 深度轉(zhuǎn)激光原理圖

    利用相機模型能夠完成世界坐標系與圖像坐標系之間的轉(zhuǎn)換[10],如式(4)所示.

    (4)

    式(4)中:(u,v)為圖像坐標系中的某一坐標點,(u0,v0)為圖像坐標系中心點,(xw,yw,zw)表示該點在世界坐標系下的坐標,其中zc表示目標到相機的距離.R、T為攝像機的外參數(shù),用來表示旋轉(zhuǎn)向量和平移向量;f為內(nèi)參數(shù),通過攝像機標定得到.

    由于世界坐標系的原點與相機的原點重合,因此沒有平移和旋轉(zhuǎn)向量,所以外參數(shù)矩陣R、T可以設(shè)置為

    (5)

    由于在相機坐標系和世界坐標系下對于同一個目標的深度值相同,即ZC=ZW.因此公式可進一步簡化為:

    (6)

    三維空間中任意一點的坐標,都能夠從以上的變換矩陣公式推導得出,因此,圖像點[u,v]T到世界坐標點[XW,YW,ZW]T的變換公式可以表示為:

    (7)

    (2)計算直線AO和DO的夾角AOD,計算公式如式(8)所示.

    θ=arctan(x/z)

    (8)

    (3)將角AOD投影到對應的激光掃描范圍內(nèi).已知激光雷達的掃描范圍為[α,β],將該范圍內(nèi)的光束分為N份,則可用數(shù)組laser[N]來表示[α,β]內(nèi)的所有激光數(shù)據(jù).因此點M投影到數(shù)組laser中的索引值n可用式(9)計算.

    n=(θ-α)/((β-α)/N)=N(θ-α)/(β-α)

    (9)

    laser[n]的值為M在x軸上投影的點D到相機光心O的距離r,即

    (10)

    因此,以激光雷達數(shù)據(jù)格式為標準對深度數(shù)據(jù)進行轉(zhuǎn)換時,通過將深度圖像在水平范圍內(nèi)取整數(shù)角度上的最小深度值作為激光雷達距離數(shù)據(jù),相當于取三維深度數(shù)據(jù)的水平投影,并加上對應角度值,這樣就完成了深度數(shù)據(jù)轉(zhuǎn)換為激光雷達數(shù)據(jù)的過程,再將數(shù)據(jù)保存到激光傳感器獲取的數(shù)組中,完成了深度相機關(guān)于前方障礙物距離信息的獲取.Ranges[]數(shù)組體現(xiàn)了橫向每個障礙點到相機的距離,完成數(shù)據(jù)獲取后,將Ranges[]數(shù)組中關(guān)于同一角度的較大的數(shù)據(jù)都用小數(shù)據(jù)來替換,這樣就完成了對于障礙物在高度上的壓縮,形成了二維平面地圖.

    2 傳感器數(shù)據(jù)融合

    當單獨使用激光雷達或深度相機進行地圖構(gòu)建時都會存在一定的誤差,因此本文先采用擴展卡爾曼濾波將激光雷達的數(shù)據(jù)與深度相機的數(shù)據(jù)進行融合,再利用貝葉斯估計的方法將深度相機與激光雷達分別生成的地圖進行二次融合,來提高建圖的準確性.

    2.1 基于擴展卡爾曼濾波的數(shù)據(jù)融合

    在本文設(shè)計的方案中,利用擴展卡爾曼濾波器融合傳感器的數(shù)據(jù)時,因為激光雷達的測量精度相較于深度相機更高,所以將激光雷達的測量值作為觀測值,將深度相機的測量值作為預測值,這樣能保證融合效果更好.令Xt=[xt,yt,θt]T作為系統(tǒng)的狀態(tài)變量,擴展卡爾曼濾波數(shù)據(jù)融合的過程如下:

    (1)預測階段

    (11)

    (12)

    (2)觀測階段

    將激光雷達的測量值ZS作為系統(tǒng)測量值,

    ZS=Ht+1Xt+1

    (13)

    式(13)中:Ht+1為激光雷達的觀測矩陣.

    (3)觀測的預測階段

    將深度相機的測量值ZV作為預測的觀測值,

    (14)

    式(14)中:Ht+1為深度相機觀測矩陣,且與激光雷達的觀測矩陣相同.

    (4)數(shù)據(jù)關(guān)聯(lián)階段

    在對激光雷達和深度相機進行數(shù)據(jù)融合之前,還需要解決兩個傳感器的數(shù)據(jù)關(guān)聯(lián)問題,目的是辨別兩個傳感器獲取的信息是不是源自于外界環(huán)境中的同一個物體[11].這是數(shù)據(jù)融合的前提,如果數(shù)據(jù)關(guān)聯(lián)出錯會導致對于環(huán)境目標的獲取出現(xiàn)錯誤,最后無法構(gòu)建出完整的地圖.

    根據(jù)激光雷達與深度相機獲取的觀測信息ZS,ZV之間的Mahalanobis距離[12],來判斷兩個傳感器獲取到的數(shù)據(jù)的關(guān)聯(lián)性,當距離小于設(shè)定閾值χ2時,即則說明關(guān)聯(lián)性較高.

    (15)

    τt=ZS-ZV

    (16)

    式(15)、(16)中:d=rank(τij),α為可信度,τt為測量的位置誤差.當ZS,ZV中有多個特征相匹配時,則選取距離最近的兩個特征作為最佳匹配,當完成了數(shù)據(jù)關(guān)聯(lián)后就能夠完成狀態(tài)的更新.

    (5)更新階段

    首先計算卡爾曼增益Kg,利用多傳感器的測量值修正之前的系統(tǒng)預測值,完成激光雷達采集的數(shù)據(jù)和深度相機采集的數(shù)據(jù)之間的融合,

    (17)

    再利用卡爾曼增益計算傳感器數(shù)據(jù)融合的最優(yōu)估計值,

    (18)

    最后利用卡爾曼增益完成對協(xié)方差矩陣的更新,

    (19)

    2.2 基于貝葉斯估計的地圖融合

    在通過擴展卡爾曼濾波完成對于傳感器的數(shù)據(jù)融合后,如果激光雷達與深度相機觀測到的信息沒有完成數(shù)據(jù)關(guān)聯(lián),那么系統(tǒng)狀態(tài)不會完成更新,之后兩個傳感器的數(shù)據(jù)也不會進行融合,環(huán)境中的某些障礙物就不會被檢測出來,進而導致地圖構(gòu)建出現(xiàn)誤差.因此在對兩個傳感器獲取的數(shù)據(jù)完成一次融合后,還應該在分別生成局部地圖之后,將兩個傳感器分別生成的局部地圖進行二次融合,以保證構(gòu)建的地圖的準確性.

    本文利用貝葉斯估計的方法將激光雷達和深度相機獲得的局部地圖完成二次融合.貝葉斯估計通過將傳感器采集的環(huán)境信息轉(zhuǎn)換成后驗概率分布函數(shù)[13],并通過貝葉斯條件概率公式更新后驗概率分布值,從而實現(xiàn)對多傳感器的信息融合[14].

    貝葉斯估計可以描述為:通過已知的向量Z,估計未知的m維狀態(tài)向量X并完成對于后驗概率的計算,已知t時刻的概率為p=(xt),且k組觀測數(shù)據(jù)為Zk=(Z1,…,Zk),先驗分布公式如下:

    (20)

    式(20)中:p(zt|xt)是基于傳感器觀測模型的似然函數(shù);p(xt|zt-1)是基于系統(tǒng)模型得到的前一時刻的先驗分布;p(zt|zt-1)表示當前時刻與上一時刻之間的狀態(tài)變化.

    在建圖的初始階段,由于環(huán)境信息未知,因此地圖中存在障礙物和無障礙物的先驗概率初始化均為0.5.隨著建圖的進行,利用傳感器不斷掃描周圍環(huán)境,通過多傳感器獲得的環(huán)境信息并利用貝葉斯估計得到傳感器信息融合的公式如下所示:

    (21)

    具體融合規(guī)則如表1所示,柵格地圖屬性分為三類:占據(jù)(O)、未占據(jù)(E)以及未知(U).

    表1 激光和視覺傳感器融合規(guī)則

    通過本文提出的方法得到的傳感器數(shù)據(jù)融合構(gòu)建地圖的流程如圖4所示.

    圖4 數(shù)據(jù)融合構(gòu)建地圖流程

    3 移動機器人的地圖構(gòu)建

    完成對傳感器數(shù)據(jù)處理之后,可以利用融合后的傳感器數(shù)據(jù),采用基于圖優(yōu)化的算法完成地圖的構(gòu)建.基于圖優(yōu)化的地圖構(gòu)建算法由圖構(gòu)建和圖優(yōu)化構(gòu)成,如圖5所示.其中圖構(gòu)建包括相鄰兩幀傳感器數(shù)據(jù)之間的匹配與回環(huán)檢測,圖優(yōu)化包括對于運動姿態(tài)的優(yōu)化[15].

    圖5 基于圖優(yōu)化的SLAM算法框架

    如圖6所示,基于圖優(yōu)化的方法通過對機器人的位姿和地圖節(jié)點之間的最小均方誤差的求解來完成對SLAM問題的求解:將各個時刻的位姿xt及地圖m中的路標點(P1,P2,…,Pn) 看成一個個待優(yōu)化的節(jié)點,將相鄰位姿節(jié)點間的邊作為一個約束條件,將位姿節(jié)點和路標節(jié)點間的邊作為另一個約束條件[16],因此,建圖的過程主要是通過傳感器對環(huán)境的掃描匹配來建立位姿節(jié)點和路標節(jié)點之間的約束關(guān)系,當滿足這個約束關(guān)系時,找到最佳的位姿節(jié)點xt和地圖m值,同時讓總體的誤差平方和最小,以此來完成地圖的構(gòu)建[17].

    圖6 基于圖優(yōu)化模型的SLAM

    在機器人仿真軟件 中搭建環(huán)境模型,并分別對基于濾波和基于圖優(yōu)化的地圖構(gòu)建算法進行仿真.其中,如圖7(a)所示為仿真環(huán)境圖,圖7(b)為濾波算法創(chuàng)建的環(huán)境地圖,圖7(c)為圖優(yōu)化算法創(chuàng)建的環(huán)境地圖.通過對比圖7(b)和圖7(c)可以發(fā)現(xiàn),基于濾波的建圖算法在構(gòu)建地圖的過程中,由于數(shù)據(jù)不能及時更新,導致誤差不斷被累積[18],因此構(gòu)建出的地圖與實際地圖相差較大.而基于圖優(yōu)化算法構(gòu)建的地圖,通過對傳感器的連續(xù)幀數(shù)據(jù)的提取,完成對位姿軌跡的確定[19],并構(gòu)建出位姿圖,同時將得到的位姿圖用回環(huán)檢測進行校正,最后通過位姿優(yōu)化消除建圖中產(chǎn)生的誤差,得到與仿真環(huán)境基本一致的地圖.

    (a) GAZEBO仿真環(huán)境圖

    (b) 濾波算法創(chuàng)建的地圖

    (c) 圖優(yōu)化算法創(chuàng)建的地圖圖7 實際地圖及算法仿真對比圖

    4 實驗結(jié)果及分析

    本文利用一臺履帶式的移動機器人進行實驗,如圖8所示.系統(tǒng)硬件框圖如圖9所示,系統(tǒng)由上位機和下位機兩部分構(gòu)成,其中下位機系統(tǒng)由STM32、電機驅(qū)動電路、慣性測量單元組成,上位機系統(tǒng)由樹莓派、深度相機和激光雷達組成,通過樹莓派的WIFI模塊可以遠程連接PC端,并且通過PC端完成對機器人的遠程控制,最后,環(huán)境地圖的構(gòu)建結(jié)果將在PC端的Rviz中顯示[20].

    實驗環(huán)境如圖10(a)和圖11(a)所示,分別是大小為6.0 mm×0.8 m的簡單走廊環(huán)境和大小為6.0 m×5.0 m的較復雜室內(nèi)環(huán)境,其中室內(nèi)環(huán)境包含桌椅、箱子等障礙物.在實際應用時,由于二維激光雷達獲取的信息會受到安裝位置的影響,深度相機獲取的信息會受到環(huán)境光線的影響,因此都會出現(xiàn)因為漏檢障礙物而導致構(gòu)建的環(huán)境地圖精度不高的情況,所以通過設(shè)計這兩種環(huán)境進行地圖構(gòu)建實驗,將兩個傳感器獲取的信息相融合,來驗證本文提出的方法的可靠性.

    圖8 移動機器人模型

    圖9 系統(tǒng)硬件框圖

    圖10(b)和圖11(b)為僅使用激光雷達對環(huán)境進行建圖.可以看出,當只使用激光雷達進行建圖時,隨著時間的推移,建圖誤差越來越大,構(gòu)建出的地圖與實際環(huán)境不一致,并且由于受激光雷達安裝位置的影響,有些障礙物如椅子(B、C、D、E)的某些部分沒有在掃描平面內(nèi),因此導致掃描獲得的地圖信息不夠全面,造成建圖精度低、噪點多的問題.

    圖10(c)和圖11(c)為僅使用深度相機對環(huán)境進行建圖.當只有深度相機完成地圖構(gòu)建時,雖然深度相機相對于激光傳感器能夠獲取到環(huán)境的三維信息,對于障礙物的檢測更加準確,但是深度相機的拍攝精度容易受到環(huán)境光線的影響,因此在建圖的準確性上較差.

    將深度相機與激光雷達的數(shù)據(jù)融合以后,構(gòu)建的地圖如圖10(d)和圖11(d)所示.通過對比可以發(fā)現(xiàn),將兩個傳感器獲取的信息融合之后,有效地彌補了單一傳感器掃描環(huán)境的局限性,克服了復雜環(huán)境的干擾,使構(gòu)建的地圖更加完整,構(gòu)建的地圖與真實環(huán)境更接近,同時地圖上的噪點也更少.

    (a)走廊環(huán)境

    (b)激光雷達構(gòu)建地圖

    (c)深度相機構(gòu)建地圖

    (d)數(shù)據(jù)融合后構(gòu)建地圖圖10 走廊環(huán)境下構(gòu)建的地圖

    (c)深度相機構(gòu)建地圖

    (d)數(shù)據(jù)融合后構(gòu)建地圖圖11 室內(nèi)環(huán)境下構(gòu)建的地圖

    移動機器人對于環(huán)境地圖的構(gòu)建精度測試可以轉(zhuǎn)換為機器人在環(huán)境中移動時對于運動軌跡的精度測試[21].機器人在不同環(huán)境下的運動軌跡的對比圖,如圖12(a)和圖13(a)所示.圖12(b)和圖13(b)所示為數(shù)據(jù)融合前后機器人運動軌跡的均方根誤差對比,可以看出,數(shù)據(jù)融合后的機器人運動軌跡明顯優(yōu)于融合之前,并且走廊環(huán)境下的均方根誤差在3 cm左右,室內(nèi)環(huán)境下的均方根誤差在8 cm左右,均小于數(shù)據(jù)融合之前的誤差,同時曲線波動幅度更小且更穩(wěn)定,表明了數(shù)據(jù)融合之后構(gòu)建的地圖更精確,更符合真實環(huán)境,從而驗證了本文算法的準確性.

    (a)走廊環(huán)境下不同軌跡對比

    (b)數(shù)據(jù)融合前后的均方根誤差圖12 走廊環(huán)境下運動軌跡和誤差對比

    (a) 室內(nèi)環(huán)境下不同軌跡對比

    (b) 數(shù)據(jù)融合前后的均方根誤差圖13 室內(nèi)環(huán)境下運動軌跡和誤差對比

    5 結(jié)論

    本文提出了一種將激光雷達和深度相機信息融合的地圖構(gòu)建方法,先利用擴展卡爾曼濾波器將移動機器人獲取的深度視覺數(shù)據(jù)與激光雷達的數(shù)據(jù)進行融合,再利用貝葉斯估計的方法將激光雷達與深度相機分別建立的局部地圖進行二次融合,并在構(gòu)建地圖時采用基于圖優(yōu)化的地圖構(gòu)建算法,在完全未知的環(huán)境中,完成對環(huán)境地圖的準確構(gòu)建.實驗結(jié)果表明,利用本文提出的多傳感器信息融合的方法,獲得的環(huán)境信息更加全面,提高了地圖構(gòu)建的可靠性和準確性,并且有效提高了移動機器人在SLAM時的魯棒性.

    猜你喜歡
    深度融合優(yōu)化
    超限高層建筑結(jié)構(gòu)設(shè)計與優(yōu)化思考
    村企黨建聯(lián)建融合共贏
    融合菜
    民用建筑防煙排煙設(shè)計優(yōu)化探討
    關(guān)于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    從創(chuàng)新出發(fā),與高考數(shù)列相遇、融合
    深度理解一元一次方程
    《融合》
    深度觀察
    久久精品国产自在天天线| 成年人黄色毛片网站| aaaaa片日本免费| 亚洲在线自拍视频| 成人特级黄色片久久久久久久| 免费在线观看影片大全网站| 精品国产三级普通话版| 一区二区三区免费毛片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 波多野结衣巨乳人妻| 亚洲中文字幕日韩| 91麻豆精品激情在线观看国产| 深爱激情五月婷婷| 免费在线观看亚洲国产| 久久精品国产综合久久久| 日本黄色视频三级网站网址| 他把我摸到了高潮在线观看| 一边摸一边抽搐一进一小说| 日本三级黄在线观看| 日本熟妇午夜| 色在线成人网| 亚洲性夜色夜夜综合| 熟女人妻精品中文字幕| 麻豆一二三区av精品| 国产99白浆流出| 色吧在线观看| 午夜精品在线福利| 国产激情欧美一区二区| 国产精品久久久久久亚洲av鲁大| 丰满人妻一区二区三区视频av | 免费一级毛片在线播放高清视频| 国内精品美女久久久久久| av欧美777| 69人妻影院| 国产成人欧美在线观看| 99久久九九国产精品国产免费| 男人舔奶头视频| 琪琪午夜伦伦电影理论片6080| 男人的好看免费观看在线视频| 亚洲一区高清亚洲精品| 成人永久免费在线观看视频| 日韩 欧美 亚洲 中文字幕| 深爱激情五月婷婷| 精品久久久久久成人av| 国产视频一区二区在线看| 国产精品av视频在线免费观看| 偷拍熟女少妇极品色| 99热6这里只有精品| 国产激情偷乱视频一区二区| 精品人妻1区二区| 欧美性猛交黑人性爽| 18禁黄网站禁片午夜丰满| 天堂av国产一区二区熟女人妻| 成人特级黄色片久久久久久久| 国产主播在线观看一区二区| 18禁裸乳无遮挡免费网站照片| 欧美日韩亚洲国产一区二区在线观看| 国产毛片a区久久久久| eeuss影院久久| 两个人视频免费观看高清| 亚洲专区国产一区二区| 亚洲自拍偷在线| 欧美一区二区亚洲| av片东京热男人的天堂| 国产精品一区二区三区四区免费观看 | 久久久久久久久久黄片| 国产精品美女特级片免费视频播放器| 成人精品一区二区免费| 亚洲av成人不卡在线观看播放网| 一个人免费在线观看的高清视频| 很黄的视频免费| 色综合婷婷激情| 亚洲精品在线观看二区| 日韩av在线大香蕉| 黑人欧美特级aaaaaa片| 伊人久久大香线蕉亚洲五| 女生性感内裤真人,穿戴方法视频| 久久久色成人| 国产精品野战在线观看| 久久精品人妻少妇| 亚洲七黄色美女视频| 99久久综合精品五月天人人| 国产精品亚洲一级av第二区| 在线免费观看的www视频| 日本一二三区视频观看| 免费在线观看日本一区| 色在线成人网| 免费看十八禁软件| 两个人看的免费小视频| 伊人久久大香线蕉亚洲五| 久久久精品欧美日韩精品| 99久久精品国产亚洲精品| 韩国av一区二区三区四区| 午夜福利18| 性色av乱码一区二区三区2| 欧美日韩亚洲国产一区二区在线观看| 亚洲av二区三区四区| www.色视频.com| 俺也久久电影网| 人妻丰满熟妇av一区二区三区| 久久香蕉精品热| 一二三四社区在线视频社区8| 性色avwww在线观看| 亚洲在线观看片| 一进一出抽搐gif免费好疼| 国产精品日韩av在线免费观看| 成人av在线播放网站| 少妇的逼水好多| av天堂在线播放| 欧美一区二区亚洲| 制服人妻中文乱码| 亚洲av成人不卡在线观看播放网| 一个人免费在线观看电影| 婷婷丁香在线五月| 午夜两性在线视频| 亚洲内射少妇av| 亚洲人与动物交配视频| 伊人久久精品亚洲午夜| 成人国产一区最新在线观看| 天美传媒精品一区二区| 亚洲天堂国产精品一区在线| 我的老师免费观看完整版| 色综合站精品国产| 又爽又黄无遮挡网站| 亚洲av五月六月丁香网| 免费电影在线观看免费观看| av女优亚洲男人天堂| 中文字幕高清在线视频| 国产一区在线观看成人免费| 久久午夜亚洲精品久久| 欧美乱码精品一区二区三区| 亚洲精品成人久久久久久| 男女视频在线观看网站免费| 窝窝影院91人妻| 91字幕亚洲| 欧美一区二区亚洲| 国产色爽女视频免费观看| 久久久久久九九精品二区国产| 波多野结衣高清作品| 99热6这里只有精品| 欧美精品啪啪一区二区三区| 亚洲18禁久久av| 国产麻豆成人av免费视频| 韩国av一区二区三区四区| avwww免费| 午夜免费激情av| 美女被艹到高潮喷水动态| 亚洲真实伦在线观看| 日韩av在线大香蕉| 麻豆一二三区av精品| 亚洲国产精品sss在线观看| 日本 欧美在线| 亚洲在线观看片| 国产av在哪里看| 丝袜美腿在线中文| 国产激情偷乱视频一区二区| 久久精品国产亚洲av香蕉五月| 久久精品91无色码中文字幕| 午夜视频国产福利| 少妇人妻一区二区三区视频| 一本一本综合久久| 在线播放无遮挡| 欧美bdsm另类| 女警被强在线播放| 亚洲av二区三区四区| АⅤ资源中文在线天堂| 少妇的丰满在线观看| av天堂中文字幕网| 极品教师在线免费播放| 精品人妻一区二区三区麻豆 | 亚洲av第一区精品v没综合| 国内久久婷婷六月综合欲色啪| 九九在线视频观看精品| 18禁在线播放成人免费| 观看免费一级毛片| av国产免费在线观看| 日韩高清综合在线| 18+在线观看网站| avwww免费| tocl精华| 久久久久国产精品人妻aⅴ院| 国产免费av片在线观看野外av| 高清毛片免费观看视频网站| 内射极品少妇av片p| 18禁国产床啪视频网站| 国内精品久久久久久久电影| 欧美三级亚洲精品| 亚洲aⅴ乱码一区二区在线播放| 国产av一区在线观看免费| 老熟妇仑乱视频hdxx| 国产精品爽爽va在线观看网站| 国产v大片淫在线免费观看| 国产午夜福利久久久久久| 免费看美女性在线毛片视频| 人妻丰满熟妇av一区二区三区| 此物有八面人人有两片| 国产亚洲av嫩草精品影院| 免费av观看视频| 热99在线观看视频| 成年免费大片在线观看| 国内精品美女久久久久久| 在线观看美女被高潮喷水网站 | 国产精品一区二区免费欧美| av天堂在线播放| 色精品久久人妻99蜜桃| 午夜日韩欧美国产| 少妇丰满av| 亚洲av免费高清在线观看| 久久精品国产亚洲av涩爱 | 午夜老司机福利剧场| 中文字幕av在线有码专区| 婷婷六月久久综合丁香| 免费看日本二区| 国产极品精品免费视频能看的| a级毛片a级免费在线| 欧美中文综合在线视频| 亚洲精品一区av在线观看| 午夜老司机福利剧场| 韩国av一区二区三区四区| 免费电影在线观看免费观看| 久久久国产成人精品二区| 亚洲 欧美 日韩 在线 免费| 免费人成视频x8x8入口观看| a级一级毛片免费在线观看| 在线观看一区二区三区| 亚洲人与动物交配视频| 日本成人三级电影网站| 亚洲av第一区精品v没综合| 一级毛片女人18水好多| 久久精品国产清高在天天线| 成熟少妇高潮喷水视频| 日韩欧美国产一区二区入口| 精品欧美国产一区二区三| 日本在线视频免费播放| 好男人电影高清在线观看| 国产精品 国内视频| 波多野结衣巨乳人妻| 一进一出好大好爽视频| 特大巨黑吊av在线直播| 久久欧美精品欧美久久欧美| 好男人在线观看高清免费视频| 亚洲午夜理论影院| 五月玫瑰六月丁香| 老司机福利观看| 精品人妻偷拍中文字幕| 男人的好看免费观看在线视频| 亚洲午夜理论影院| 国产精品三级大全| 制服人妻中文乱码| 国产精品av视频在线免费观看| 精品一区二区三区视频在线 | 桃红色精品国产亚洲av| 亚洲欧美日韩东京热| 少妇的逼好多水| 蜜桃久久精品国产亚洲av| 神马国产精品三级电影在线观看| 69av精品久久久久久| 波多野结衣巨乳人妻| 国产精品免费一区二区三区在线| 国产精品亚洲av一区麻豆| 国内精品美女久久久久久| 久久6这里有精品| 亚洲中文日韩欧美视频| 有码 亚洲区| 99久国产av精品| 国产免费一级a男人的天堂| 99久久无色码亚洲精品果冻| 欧美绝顶高潮抽搐喷水| 国产精品三级大全| 看免费av毛片| 中文字幕人妻丝袜一区二区| 亚洲男人的天堂狠狠| 国产av在哪里看| 久久久色成人| 精品一区二区三区视频在线观看免费| 久久6这里有精品| 欧美高清成人免费视频www| 97超视频在线观看视频| 精品99又大又爽又粗少妇毛片 | 九九在线视频观看精品| 中亚洲国语对白在线视频| 国产熟女xx| АⅤ资源中文在线天堂| 天堂网av新在线| 热99在线观看视频| 婷婷丁香在线五月| 一级a爱片免费观看的视频| 国产麻豆成人av免费视频| 亚洲国产高清在线一区二区三| 欧美一区二区亚洲| 国产亚洲精品av在线| 哪里可以看免费的av片| 亚洲在线观看片| 日韩欧美一区二区三区在线观看| 欧美丝袜亚洲另类 | 国产探花极品一区二区| 超碰av人人做人人爽久久 | 在线观看一区二区三区| 欧美成人免费av一区二区三区| 麻豆久久精品国产亚洲av| 啪啪无遮挡十八禁网站| 色老头精品视频在线观看| 欧美性猛交╳xxx乱大交人| 男女那种视频在线观看| 欧美黄色片欧美黄色片| 久久国产精品人妻蜜桃| 国产真人三级小视频在线观看| 成人午夜高清在线视频| 91九色精品人成在线观看| 久久久久久久久大av| 午夜精品久久久久久毛片777| 搡老熟女国产l中国老女人| 国产成人福利小说| 99热精品在线国产| 亚洲精品粉嫩美女一区| 国产激情欧美一区二区| 国产精品三级大全| 午夜视频国产福利| 久久久国产精品麻豆| 人妻丰满熟妇av一区二区三区| 国产高清激情床上av| 免费看日本二区| 最好的美女福利视频网| 两人在一起打扑克的视频| 国产精品综合久久久久久久免费| 日韩欧美精品免费久久 | 亚洲无线观看免费| 黄色视频,在线免费观看| 脱女人内裤的视频| АⅤ资源中文在线天堂| 国产av麻豆久久久久久久| 深夜精品福利| 国产又黄又爽又无遮挡在线| 欧美黄色片欧美黄色片| 真实男女啪啪啪动态图| 亚洲,欧美精品.| 级片在线观看| 久久精品国产清高在天天线| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 久久这里只有精品中国| 校园春色视频在线观看| 少妇丰满av| 人妻久久中文字幕网| 亚洲乱码一区二区免费版| 免费在线观看成人毛片| 亚洲乱码一区二区免费版| aaaaa片日本免费| 此物有八面人人有两片| 精品无人区乱码1区二区| 国产乱人视频| 午夜激情欧美在线| 精品人妻1区二区| 精品无人区乱码1区二区| 又粗又爽又猛毛片免费看| 国产真实乱freesex| 国产高潮美女av| 国产淫片久久久久久久久 | 日本黄大片高清| 久久久久九九精品影院| 国产极品精品免费视频能看的| 久久久久久久精品吃奶| 黄色成人免费大全| 久久久久国内视频| 欧美三级亚洲精品| 亚洲国产色片| 久久久久久人人人人人| 黄色片一级片一级黄色片| 精品电影一区二区在线| 狂野欧美白嫩少妇大欣赏| 日韩人妻高清精品专区| 国产日本99.免费观看| 男人舔奶头视频| 国内精品久久久久久久电影| 日韩欧美精品v在线| 少妇的丰满在线观看| 美女大奶头视频| 蜜桃亚洲精品一区二区三区| 国产精品香港三级国产av潘金莲| 特大巨黑吊av在线直播| 黄色女人牲交| 国产高清videossex| 欧美午夜高清在线| 精品久久久久久久久久免费视频| 99热只有精品国产| 99热精品在线国产| 国产一区二区三区在线臀色熟女| 国产亚洲精品综合一区在线观看| 国产精品爽爽va在线观看网站| 成人欧美大片| 听说在线观看完整版免费高清| 可以在线观看毛片的网站| 国产亚洲欧美在线一区二区| 少妇的逼水好多| 最近视频中文字幕2019在线8| 99精品在免费线老司机午夜| 国产成人福利小说| 韩国av一区二区三区四区| 超碰av人人做人人爽久久 | 19禁男女啪啪无遮挡网站| 成人一区二区视频在线观看| 淫秽高清视频在线观看| 国产亚洲精品一区二区www| 日韩欧美免费精品| 国产毛片a区久久久久| 嫩草影院入口| 精品国产超薄肉色丝袜足j| 1024手机看黄色片| 成人国产一区最新在线观看| 女人被狂操c到高潮| av国产免费在线观看| 内地一区二区视频在线| 18禁国产床啪视频网站| 日日夜夜操网爽| a级毛片a级免费在线| 亚洲无线观看免费| 久久亚洲真实| 夜夜看夜夜爽夜夜摸| 亚洲国产色片| 亚洲乱码一区二区免费版| 欧美成人免费av一区二区三区| 亚洲乱码一区二区免费版| 特大巨黑吊av在线直播| 国产精品一及| 窝窝影院91人妻| 一个人免费在线观看的高清视频| 亚洲成av人片在线播放无| 噜噜噜噜噜久久久久久91| 内射极品少妇av片p| 亚洲黑人精品在线| 亚洲人成网站在线播放欧美日韩| 欧美乱妇无乱码| 精品电影一区二区在线| www国产在线视频色| 一a级毛片在线观看| 国产av一区在线观看免费| 五月玫瑰六月丁香| 亚洲午夜理论影院| 亚洲七黄色美女视频| 国产精品99久久99久久久不卡| 国产精品久久视频播放| 悠悠久久av| 国内精品久久久久久久电影| 老司机午夜福利在线观看视频| 国产精品影院久久| 蜜桃亚洲精品一区二区三区| 99久久精品热视频| 免费人成在线观看视频色| 国产成人福利小说| 欧美绝顶高潮抽搐喷水| 一边摸一边抽搐一进一小说| 国产精品综合久久久久久久免费| 国产成人av教育| 国产野战对白在线观看| 欧洲精品卡2卡3卡4卡5卡区| 中文字幕人妻丝袜一区二区| 韩国av一区二区三区四区| 欧美zozozo另类| 全区人妻精品视频| 非洲黑人性xxxx精品又粗又长| 亚洲精品在线观看二区| 蜜桃久久精品国产亚洲av| 亚洲av美国av| 久久久久久久久中文| а√天堂www在线а√下载| 波野结衣二区三区在线 | 欧美不卡视频在线免费观看| 最后的刺客免费高清国语| 亚洲 欧美 日韩 在线 免费| 免费在线观看亚洲国产| 国产高清激情床上av| 亚洲欧美日韩高清在线视频| 中文字幕熟女人妻在线| 国产精品98久久久久久宅男小说| 在线观看午夜福利视频| 日韩人妻高清精品专区| 色av中文字幕| 老汉色av国产亚洲站长工具| 国产伦在线观看视频一区| 嫩草影院入口| 日韩免费av在线播放| 欧美丝袜亚洲另类 | 亚洲美女黄片视频| 日本三级黄在线观看| 天堂√8在线中文| 国产精品亚洲美女久久久| 久久精品91无色码中文字幕| 日韩人妻高清精品专区| 一卡2卡三卡四卡精品乱码亚洲| 草草在线视频免费看| 国产成人a区在线观看| 午夜福利免费观看在线| 99国产精品一区二区三区| 一本久久中文字幕| 国产精品影院久久| 国产精品永久免费网站| 九色成人免费人妻av| 又黄又粗又硬又大视频| 俺也久久电影网| 日韩欧美国产一区二区入口| 日本三级黄在线观看| 国产成人福利小说| 国产精品国产高清国产av| 午夜精品在线福利| 欧美av亚洲av综合av国产av| 日本一本二区三区精品| 国产精品女同一区二区软件 | 色综合婷婷激情| 欧美zozozo另类| 精品一区二区三区人妻视频| 久久久国产成人精品二区| 国产精品免费一区二区三区在线| 免费在线观看日本一区| 亚洲欧美一区二区三区黑人| 成人特级黄色片久久久久久久| 熟女人妻精品中文字幕| 一区福利在线观看| 欧美最黄视频在线播放免费| 亚洲欧美一区二区三区黑人| 国产一区二区三区视频了| 成人鲁丝片一二三区免费| 国产精品,欧美在线| 色综合婷婷激情| 精品久久久久久久末码| 波多野结衣高清无吗| 国产久久久一区二区三区| 久久久久久国产a免费观看| 99久久九九国产精品国产免费| 久久性视频一级片| 无遮挡黄片免费观看| 国产高清三级在线| 午夜亚洲福利在线播放| 动漫黄色视频在线观看| 亚洲精品456在线播放app | 国内揄拍国产精品人妻在线| 国产精品野战在线观看| 男女视频在线观看网站免费| 97碰自拍视频| 日本在线视频免费播放| 嫩草影院精品99| 女同久久另类99精品国产91| 97超级碰碰碰精品色视频在线观看| 男女下面进入的视频免费午夜| xxx96com| 超碰av人人做人人爽久久 | 两个人视频免费观看高清| 我的老师免费观看完整版| 特大巨黑吊av在线直播| 麻豆国产97在线/欧美| 国产精品国产高清国产av| 亚洲人成电影免费在线| 国产精品香港三级国产av潘金莲| 夜夜躁狠狠躁天天躁| 一夜夜www| 国产免费一级a男人的天堂| 成熟少妇高潮喷水视频| 久久精品综合一区二区三区| 一级黄色大片毛片| 在线观看免费午夜福利视频| 国产精品av视频在线免费观看| 欧美一级a爱片免费观看看| 两个人视频免费观看高清| 老熟妇仑乱视频hdxx| 国产精品久久电影中文字幕| 露出奶头的视频| 日韩中文字幕欧美一区二区| 国产成人欧美在线观看| 国产国拍精品亚洲av在线观看 | 精品久久久久久久毛片微露脸| 国内精品一区二区在线观看| 久99久视频精品免费| 国产伦精品一区二区三区视频9 | 欧美午夜高清在线| 国产亚洲精品一区二区www| 给我免费播放毛片高清在线观看| 亚洲成人精品中文字幕电影| 国产成人av激情在线播放| 欧美成人性av电影在线观看| 国产爱豆传媒在线观看| 精品久久久久久久久久久久久| 久久久久久久午夜电影| 一级黄色大片毛片| 欧美丝袜亚洲另类 | 日本 av在线| 一区二区三区国产精品乱码| 俺也久久电影网| 亚洲精品成人久久久久久| 亚洲av中文字字幕乱码综合| 99riav亚洲国产免费| 亚洲欧美日韩无卡精品| 国产精品野战在线观看| 日韩高清综合在线| 精华霜和精华液先用哪个| 一级毛片高清免费大全| 男插女下体视频免费在线播放| 人人妻人人澡欧美一区二区| 精品一区二区三区av网在线观看| 我要搜黄色片| 日韩免费av在线播放| 亚洲天堂国产精品一区在线| 极品教师在线免费播放| 99久久久亚洲精品蜜臀av| 国内久久婷婷六月综合欲色啪| 床上黄色一级片| 在线观看免费视频日本深夜| 国产精品亚洲av一区麻豆| 亚洲内射少妇av| 久久国产精品影院| 国产精品久久久久久久久免 | 狠狠狠狠99中文字幕| 国内少妇人妻偷人精品xxx网站| 欧美中文综合在线视频| 18禁在线播放成人免费| 成人av一区二区三区在线看| 级片在线观看| 成人午夜高清在线视频|