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

    基于車載4線激光雷達的前方道路可行駛區(qū)域檢測*

    2016-03-30 08:16:44段建民李龍杰鄭凱華
    汽車技術(shù) 2016年2期
    關(guān)鍵詞:激光雷達障礙物聚類

    段建民 李龍杰 鄭凱華

    (北京工業(yè)大學(xué)北京市交通工程重點實驗室,北京100124)

    ?

    基于車載4線激光雷達的前方道路可行駛區(qū)域檢測*

    段建民李龍杰鄭凱華

    (北京工業(yè)大學(xué)北京市交通工程重點實驗室,北京100124)

    【摘要】通過分析大量激光雷達數(shù)據(jù),提出了改進的DBSCAN算法對雷達數(shù)據(jù)進行聚類,以使智能汽車獲得前方道路和障礙物信息。根據(jù)不同的密度參數(shù),多次調(diào)用該算法完成多密度聚類。結(jié)合提出的道路邊沿、路面和障礙物等信息提取方法,在結(jié)構(gòu)化或半結(jié)構(gòu)化的城市道路中對提出的方法進行實車試驗,結(jié)果表明該算法可實時準確的提取智能汽車的可行駛區(qū)域信息。

    1 前言

    智能汽車利用可通行區(qū)域信息來保障其在道路區(qū)域內(nèi)安全行駛,防止側(cè)翻和碰撞等危險狀況發(fā)生。車輛可通行區(qū)域的主要信息有邊界的位置信息(路沿、路口)和道路內(nèi)的障礙物信息等。用于環(huán)境信息檢測的傳感器主要有攝像頭和激光雷達等。由于攝像頭采集數(shù)據(jù)時受環(huán)境影響較大,因此本文研究其中一種受環(huán)境因素影響較小的車載傳感器——激光雷達。激光雷達具有全天候、遠距離、不受光線影響等優(yōu)點,可以得到汽車前方環(huán)境的高精度距離和方位數(shù)據(jù)構(gòu)成的深度信息。

    目前使用激光雷達提取道路信息的研究成果可參考文獻[1]~文獻[4],從中可以看出,對雷達掃描數(shù)據(jù)進行處理的方法主要使用聚類分析,其被廣泛運用于數(shù)據(jù)挖掘技術(shù)中,是一種常用的模式識別方法。常用的聚類方法有基于劃分的聚類、基于層次聚類、基于密度的聚類、基于網(wǎng)格的聚類、基于神經(jīng)網(wǎng)絡(luò)的聚類、基于統(tǒng)計學(xué)的聚類和高維可視化數(shù)據(jù)聚類等。K-均值算法一般實現(xiàn)比較簡單,適合用來處理大量數(shù)據(jù)甚至是高維數(shù)據(jù),但是K-均值算法具有需要預(yù)先給定K值和對外部噪聲敏感等缺點[4]。ISODATA算法在這兩個方面作了改進,但是同時也增加了許多限制條件,而這些限制條件在實際應(yīng)用中難以滿足。根據(jù)激光雷達的數(shù)據(jù)特點,當激光雷達掃描到物體表面時數(shù)據(jù)點會呈現(xiàn)密集分布特性,因此基于密度的方法更適合處理雷達數(shù)據(jù)。本文改進了基于密度的DBSCAN聚類算法,運用一維核密度估計和加權(quán)歐式距離來解決參數(shù)Eps和Minpts的選取問題,改進搜索范圍來提高運行速度,再提出可行駛區(qū)域方法。試驗結(jié)果表明,該算法可以很好的提取道路信息和得到可通行區(qū)域信息。

    2 激光雷達數(shù)據(jù)特點和坐標轉(zhuǎn)換

    2.1激光雷達數(shù)據(jù)特點分析

    本文采用德國生產(chǎn)的IBEO-LUX-2010型號4線激光雷達,它是一種三維的4層激光雷達,利用多層技術(shù)可以獲得高度信息,與單線激光雷達相比具有數(shù)據(jù)量大、檢測精度高、探測距離遠、范圍廣等優(yōu)點。雖然與64線激光雷達相比還有一定差距,但是對于智能汽車而言,4線激光雷達性價比較高,在數(shù)據(jù)量和實時性方面能很好的滿足實際需求。該型號雷達垂直視場角度為2.4°,4條掃描線在垂直方向以不同的俯仰角同時在水平方向上構(gòu)成4個掃描面,相鄰的掃描面之間夾角為0.8°,水平視場寬,最大掃描角度范圍約為110°,掃描頻率為12.5 Hz,水平分辨率為0.125°,距離分辨率為4 cm,測量范圍為0.3~200 m,性能穩(wěn)定,具有全天候的工作特性,其具體參數(shù)如表1所列。

    表1 IBEO-LUX-2010型激光雷達技術(shù)參數(shù)

    IBEO-LUX-2010每次掃描采集到的數(shù)據(jù)包含周圍物體的距離、位置和掃描層數(shù)等信息,通過以太網(wǎng)進行數(shù)據(jù)信息的傳遞,以1幀數(shù)據(jù)為單位進行數(shù)據(jù)處理。

    在給定的初始化條件下,通過對采集的單幀數(shù)據(jù)進行分析,得到如下特點:

    a.激光雷達數(shù)據(jù)為三維數(shù)據(jù),掃描點自左向右順序排列,掃描方式為周期掃描,每幀的返回數(shù)據(jù)點數(shù)目相差不大;

    b.激光雷達返回的掃描數(shù)據(jù)時刻隨著行駛環(huán)境的變化而變化,因為需要根據(jù)智能汽車周圍環(huán)境進行聚類,所以聚類的數(shù)目不確定,形狀也是任意的,這就要求在實際聚類時聚類算法能根據(jù)數(shù)據(jù)特性對聚類數(shù)目有自適應(yīng)性;

    c激光雷達掃描到道路路面上的點在x、y軸方向上基本為均勻緊密分布,距離很??;掃描到路沿(一般道路方向與車輛行駛方向同向)上的點為同一層中相鄰的點,x軸方向分布稀疏、y軸方向分布特別緊密,局部呈直線形;

    d.激光雷達掃描到前方物體上的點,在相同掃描角度方向上,隨著障礙物距離的增大,返回的相鄰數(shù)據(jù)點間距變大,呈現(xiàn)正相關(guān)性;等距離情況時,隨著掃描到障礙物上的點與y軸夾角(0°~90°)的增大而增大,同樣也呈現(xiàn)正相關(guān)性。

    根據(jù)雷達數(shù)據(jù)的掃描和分布特性,本文對道路信息提取時,主要使用最低的兩層數(shù)據(jù)來提取車輛前方的路面寬度和車輛所處位置等信息,而最高和次高層數(shù)據(jù)則用來探測較遠處的路面、丁字路口和阻擋物等信息,將4層數(shù)據(jù)融合處理后的結(jié)果與實際場景、高程數(shù)據(jù)信息等多特征比對,以保證提取的信息可靠。

    2.2坐標轉(zhuǎn)換

    從車輛的安全角度出發(fā),為保證激光雷達對前方道路觀察的全面和準確,將4線激光雷達安裝在智能汽車正前方的中間位置,距離地面高度h為0.846 m,安裝示意圖如圖1所示。

    圖1 激光雷達安裝示意圖

    車體坐標系是以車輛后軸中心為坐標原點,前進方向為xv軸,水平向左為yv軸,豎直向上為zv軸。建立臨時坐標系以激光雷達在地面的投影中心為坐標原點OT,水平向前為xT軸,水平向左為yT軸,垂直于xoy平面豎直向上為zT軸。激光雷達坐標系以激光雷達的中心為坐標原點,取第1層藍層為xoy平面,水平向前為xL軸,水平向左為yL軸,垂直于藍層豎直向上為zL軸,如圖2所示。為保證在水平地面時能有兩條掃描線掃描到路面,以準確對道路信息進行提取,故安裝雷達時向下傾斜α>2.4°,如圖2所示?;谥悄芷嚨膽?yīng)用需要將傳感器坐標轉(zhuǎn)換到車體坐標系中,IBEO-LUX-2010激光雷達的垂直掃描精度為固定值θ0=0.8°,所以本文使用一種分層轉(zhuǎn)換4線激光雷達坐標的方法[5,6],關(guān)系式為:

    式中,ρ是雷達原點到掃描點的距離值;θ是激光雷達掃描過的角度值;i分別代表4層激光雷達的層號;ρ1(θ=0° )是智能汽車在水平地面上第1層掃描線掃描角度θ=0°時雷達坐標原點到地面距離值;R3×3為安裝旋轉(zhuǎn)矩陣;T1×3為平移矩陣,代表從臨時坐標系的原點到車輛坐標系的坐標原點的平移量,可由實際測量得到。

    圖2 雷達數(shù)據(jù)的坐標轉(zhuǎn)換

    由于雷達安裝固定后與車體的姿態(tài)時刻一致,因此利用車輛的位姿傳感器(如陀螺儀等)測量車輛姿態(tài),得到旋轉(zhuǎn)矩陣R′,利用式(1)對雷達采集到的數(shù)據(jù)進行預(yù)處理,可以很好的處理車輛俯仰及道路結(jié)構(gòu)導(dǎo)致的激光雷達相對俯仰角起伏變化情況。

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

    3.1DBSCAN算法

    DBSCAN[7]聚類算法是一種經(jīng)典的基于密度的聚類算法。該算法首先遍歷數(shù)據(jù)集中的每個數(shù)據(jù)對象,若該數(shù)據(jù)對象的半徑鄰域內(nèi)包含的對象個數(shù)不少于給定的閾值,則將該數(shù)據(jù)對象標記為核心對象點,迭代尋找其附近密度可達的數(shù)據(jù)對象,標記為同類號形成一個簇,否則將剩余的點標記為噪聲點。DBSCAN算法的優(yōu)點是能夠發(fā)現(xiàn)任意形狀的簇,且對噪聲數(shù)據(jù)點不敏感。該算法需要預(yù)先給定領(lǐng)域搜索半徑Eps和領(lǐng)域內(nèi)成為核心對象的最小數(shù)據(jù)個數(shù)MinPts。傳統(tǒng)的DBSCAN算法描述為[8]:

    a.檢查輸入的數(shù)據(jù)集合D中未被檢查過的數(shù)據(jù)點p,假如p未被處理過(未標記為某一類或歸為噪聲),檢查p的Eps領(lǐng)域內(nèi)數(shù)據(jù)點個數(shù)是否大于MinPts,若是,p標記為核心對象,且建立新簇C,將其領(lǐng)域內(nèi)的直接密度可達點歸入此簇。

    b.檢查C中所有未被處理的數(shù)據(jù)點p′的Eps領(lǐng)域,若包含不少于MinPts個數(shù)據(jù)點,則將其領(lǐng)域中未處理的數(shù)據(jù)點歸入簇C中。

    c.重復(fù)步驟b直至沒有新的對象加入簇C中,然后進行下一步。

    d.重復(fù)上述3步驟,形成最大密度相連對象的集合,中間包含有一些密度可達對象集合的合并等操作。直到所有數(shù)據(jù)都被處理完畢,最后輸出帶有簇號或噪聲標記的數(shù)據(jù)集D′。

    3.2一維核密度估計方法

    基于前述激光雷達數(shù)據(jù)特征分析,本文選用一維核密度估計方法對一幀數(shù)據(jù)的數(shù)據(jù)特性進行初步統(tǒng)計,以提取需要的數(shù)據(jù)集。核密度估計是由Rosenblatt[9]和Parzen[10]等人提出的一種非參數(shù)估計方法,該方法可以不依賴于數(shù)據(jù)分布的先驗知識,從數(shù)據(jù)本身出發(fā)得到數(shù)據(jù)的分布特性。將激光雷達數(shù)據(jù)點的某一維度看作是一組服從獨立同分布的隨機變量x1,x2,……,xn,假設(shè)隨機變量服從同一概率密度函數(shù)f(x ),那么定義:

    核函數(shù)都滿足Mercer定理,很容易確定。常用的核函數(shù)有Uniform、Triangular、Epanechnikov和Gaussian等,這些單個核函數(shù)與坐標軸形成的面積為1。本文采用常用的Epanechnikov核函數(shù),其表達式如下:

    該核密度函數(shù)不僅形式簡單,而且在所有核函數(shù)中其效率也是接近最優(yōu)的。在核密度估計中窗寬的選擇也很重要,若窗寬選擇過小,容易出現(xiàn)波動,呈現(xiàn)不規(guī)則形態(tài),可能會掩蓋f(x)的重要特性;若窗寬選擇過大則經(jīng)壓縮后使得平均作用突出,會使(x)過于平穩(wěn),減小靈敏性。結(jié)合經(jīng)驗和理論分析,運用交叉驗證策略選擇合適的窗寬h來平衡這兩種效應(yīng)[11],本文選擇的窗寬為:

    式中,σ?為待估計數(shù)據(jù)標準差的無偏估計。

    由于相鄰多幀間數(shù)據(jù)相似度較大,結(jié)合當前幾幀待估計數(shù)據(jù)的σ?和n的中間值來進行計算,以防止出現(xiàn)突變值。經(jīng)試驗驗證所得的h能滿足實際應(yīng)用要求。

    3.3加權(quán)歐式距離

    加權(quán)歐式距離可以看作是對簡單歐式距離方法的一種改進,可以利用權(quán)重值將數(shù)據(jù)各個維度的重要性區(qū)別開來。設(shè)n維向量p1(x11,x12,……,x1n)、p2(x21,x22,……,x2n),加權(quán)歐式距離的表達公式如下:

    式中,當λk=1時,為歐式距離;當λk≠1時,為加權(quán)歐式距離;當n=2時,dij為平面內(nèi)兩點的距離;當n=3時,dij為空間內(nèi)兩點間的距離。

    對雷達數(shù)據(jù)進行核密度估計后提取各個峰值區(qū)域的點集,計算其橫、縱坐標的標準差值,對標準差較大的維度賦予較大的權(quán)重值。為了計算方便,將λ2進行歸一化處理,λ值的確定可根據(jù)當前幀核密度估計的峰值區(qū)域點集的橫縱坐標標準差求得,若σ′≠0,則;否則σ′=ε(很?。?。從而可以對掃描到路沿的點的x、y賦予不同權(quán)重值,而對于掃描到障礙物上的點,也可以采用此類辦法進行更精確的聚類。

    3.4最小二乘法的擬合道路邊界

    根據(jù)雷達的掃描數(shù)據(jù)特點,掃描到道路邊沿上的點一般都是呈線性分布,因此本文采用最小二乘法對聚類產(chǎn)生的數(shù)據(jù)點進行擬合,得出被掃描的道路邊沿表面輪廓。

    最小二乘法常用于直線段y=kx+b擬合,其k、b的表達式為:

    式中,n表示待擬合的點的個數(shù);(xi,yi) 為待擬合點的坐標值;k代表斜率;b代表截距。

    4 道路信息的提取方法

    4.1改進的DBSCAN聚類算法

    DBSCAN聚類方法能夠在存在噪聲的空間中發(fā)現(xiàn)任意數(shù)目和任意形狀的簇,但該聚類算法還存在以下不足[12,13]:

    a.聚類時必須預(yù)先手工對Eps(半徑鄰域)和MinPts(核心對象鄰域內(nèi)的最少點數(shù))兩個參數(shù)進行賦值,由于DBSCAN的聚類效果受這兩個參數(shù)準確性的影響嚴重,因此根據(jù)待處理數(shù)據(jù)的特征對這兩個參數(shù)值的自適應(yīng)確定至關(guān)重要;

    b.通用的DBSCAN算法僅對這兩個參數(shù)進行一次賦值,只能對一個全局的單一密度聚類,而在處理密度多變的數(shù)據(jù)集時,往往得不到理想的聚類效果;

    c.聚類算法在進行確立核心對象運算時,通常采用遍歷法,且隨著數(shù)據(jù)的遍歷核心對象不斷增多,會導(dǎo)致運算速度較慢且內(nèi)存耗費很大。

    針對DBSCAN算法的缺點結(jié)合激光雷達數(shù)據(jù)的特點,本文提出了以下改進措施:

    a.針對參數(shù)敏感且無法適用于多密度數(shù)據(jù)集聚類的缺點,運用一種改進的基于一維估計分析的多密度聚類算法[13]來對激光雷達數(shù)據(jù)進行聚類。即對一幀激光雷達數(shù)據(jù)進行一維核密度估計分析,估計激光雷達數(shù)據(jù)的特征分布,結(jié)合加權(quán)歐式距離,得到估計的峰值點附近的數(shù)據(jù)密度和估計的4-dist值(4最近鄰點的距離),對4- dist值從低到高進行排序[sort0,sort1,……,sortP];在聚類時,先聚類出密度大的數(shù)據(jù),再依次聚類出密度小的數(shù)據(jù),從而作為DBSCAN的輸入?yún)?shù)進行多密度聚類,針對不同障礙物數(shù)據(jù)會有不同的密度輸入和處理提取方法。

    b.為了降低算法的搜索復(fù)雜度,根據(jù)激光雷達掃描特性,運用采用峰值近鄰域掃描+分步處理的方法來代替常規(guī)DBSCAN方法中的遍歷法。現(xiàn)將KDE核密度估計的值高于一定密度值的區(qū)域[作為待處理區(qū)域,ζ表示核密度估計的鄰域范圍[12]。分步處理是將4線激光雷達分為兩部分,最低和次低的兩層由于直接掃描到路面上,主要用來提取橫向道路(如邊界等)信息,而最高和次高層由于掃描射程較遠,主要用來提取路面前方縱向道路(如障礙物等)信息(必要時可利用局部4層數(shù)據(jù))。這樣兩層結(jié)合處理,既可以大大降低由于單層激光雷達處理時的虛警和誤報,又可以更好的探測車輛前方足夠大的行駛區(qū)域,同時降低了算法復(fù)雜度。

    4.2道路邊沿信息提取算法

    首先對1層雷達數(shù)據(jù)的橫坐標y進行一維核密度估計分析,尋找密度估計的峰值點,作為DBSCAN算法的核心對象,再對峰值點附近滿足式(7)條件的密度數(shù)據(jù)點區(qū)域進行提取:

    這樣可以最大限度的保留掃描到路沿上的點,同時將一些散亂點數(shù)據(jù)濾除,如掃到的花草等。將提取出的目標區(qū)域數(shù)據(jù)點按密度從高到低進行排序,以聚類出密度大的簇,再依次對剩余的數(shù)據(jù)進行聚類,以把不同密度的數(shù)據(jù)點聚成不同的類別。而在對目標區(qū)域聚類前,首先結(jié)合加權(quán)歐式距離計算峰值點的4-dist近鄰點距離,作為DBSCAN的Eps參數(shù)輸入;然后將自適應(yīng)產(chǎn)生的參數(shù)輸入到改進的DBSCAN算法中,從而將線形分布的路沿數(shù)據(jù)點聚成若干簇,再運用最小二乘法擬合出道路邊沿線若干條。同理,對次低層數(shù)據(jù)進行同樣的分析操作,得到若干條道路邊沿線,保留距離車輛最近的左右兩邊各一條道路邊沿線。利用得到的高程信息z值對提取出的路沿點高度(一般為20 cm左右)特征進行驗證,再結(jié)合場景圖像實際對比和試驗驗證本方法可以很好的提取出道路邊沿線。

    4.3可行駛區(qū)域的檢測

    根據(jù)前述提取出的路沿信息,可以橫向界定可行駛區(qū)域范圍,為車輛橫向控制提供依據(jù),同時需要進一步在此范圍內(nèi)對車輛縱向行駛區(qū)域進行檢測。其過程可分為道路路面的檢測提取和對路面上凹凸障礙物信息的提取兩部分。為了減小計算復(fù)雜度,根據(jù)路面上掃描點之間緊密連續(xù)排列的數(shù)據(jù)特點,將預(yù)處理后的雷達三維數(shù)據(jù)在xoy、yoz的二維平面中進行投影。在yoz平面內(nèi),結(jié)合高程信息,可以將雷達各掃描層掃描到路面點部分提取出來。在一定范圍內(nèi),路面具有一定的連續(xù)性,因此其高度也應(yīng)該在一定范圍內(nèi)變化。設(shè)當前車輛所處位置高度為he,由于地面起伏不平造成的局部高度差極值為τ,則激光雷達所掃描到的部分道路高度[2]應(yīng)滿足[he-τ,he+τ]。由于路面上的點彼此間距相差不大,且連續(xù)密集分布,在xoy平面內(nèi),可以進一步將剩余路面點滿足式(8)條件標記為路面點:

    式中,i,j為同一掃描層相鄰的點的序號;α為比例系數(shù);θs為雷達各層的水平分辨率。

    在濾除道路以外和路面的掃描點數(shù)據(jù)后,剩余的道路內(nèi)數(shù)據(jù)點聚類后主要分為兩種情況:凹障礙物,具有負高度值,如坑洼地區(qū)等,對智能汽車的行駛有一定安全隱患;凸障礙物,具有正高度值,如行人、車輛、路口的路肩或墻壁等。當這兩類障礙物出現(xiàn)在道路范圍內(nèi)后,需要對它們的信息進行提取。因為雷達掃描方式是自左向右,在掃描過程中只掃到同角度且距離雷達近的物體表面,所以對道路內(nèi)最高兩層(必要時局部4層)激光雷達的數(shù)據(jù)y值運用改進后的DBSCAN算法,對道路內(nèi)的障礙物進行不同密度的聚類。在yoz平面利用高程值提取出凸凹障礙物的最大高(深)度信息。通過障礙物高(深)度信息可以判斷出此障礙物所屬的類別,從而可以得到道路區(qū)域內(nèi)的障礙物屬性信息表示object= {distance,center_x,center_y,object_width,angle}。

    前方道路被分割成若干部分,通過計算其與左右路沿的距離,可通行區(qū)域表示如下:

    式中,i代表左、右路沿;ζ表示車輛通行的寬度安全系數(shù)。

    5 實車試驗與結(jié)果分析

    5.1道路路沿、路面檢測和提取試驗

    為了驗證本文所提算法檢測路沿和路面的效果,在以下道路場景中進行試驗。試驗平臺為BJUT-IV自主研發(fā)的智能汽車平臺,裝備有型號為IBEO-LUX-2010 的4線激光雷達,前向安裝俯仰角角度約為2.88°,橫滾角約為-1.1°,正常行駛時保證在水平地面上有兩條掃描線掃描到路沿上。試驗場景1是校園內(nèi)一條雙向單車道的結(jié)構(gòu)化道路;試驗場景2是一條沒有標準車道和道路邊沿的道路。試驗場景1和場景2進行道路邊沿提取和路面提取。

    圖3a所示為校園內(nèi)一條結(jié)構(gòu)化的雙向單車道道路,具有標準高度的路沿作為道路分界,道路兩邊有人行道,遠近道路寬度近似一致。圖3b是對原始一幀雷達數(shù)據(jù)的第0、1層進行一維(y坐標值)核密度估計(KDE)結(jié)果,可以看出掃描到路沿的點的密度估計值較大。由于人行道左側(cè)墻壁的干擾最低和次低層利用最小二乘法共擬合出5條路沿線(圖3c中的橢圓形標示路沿點),將距離較近的進行合并,再提取出距離車體位置最近的左右兩條路沿線(直線標示)。通過將提取出的數(shù)據(jù)點投影到y(tǒng)oz平面,進一步通過高度信息對比,對所提取的路沿進行確認,高度約為0.2 m,符合一般結(jié)構(gòu)化道路的路沿特征。

    圖4所示為校園內(nèi)一條半結(jié)構(gòu)化的道路,一側(cè)具有標準的道路邊界,而另一側(cè)為具有一定高度的圍墻,并且距離車輛較近的道路寬度較窄,而較遠處較寬(道路中部凹陷是因為中部有排水設(shè)施)。同樣進行圖3的操作,可以得到兩條路沿擬合線。

    圖3 結(jié)構(gòu)化城市道路試驗結(jié)果

    圖4 半結(jié)構(gòu)化城市道路試驗結(jié)果

    在提取的道路邊沿范圍內(nèi)對路面點進行提取,由于路面點高度在參考平面附近,且掃描到路面上的點為同一層連續(xù)緊密排列的點,因此先利用高度條件提取出部分路沿點集,再利用式(7)來將符合特征的點歸為路面點集中。從對結(jié)構(gòu)化和半結(jié)構(gòu)化的道路路面提取試驗可以證明,所提出的路面提取方法可以很好的提取路面點,如圖3c和4c中的路面點標示。再與圖3d、4d的高度信息作對比,高度在參考水平面附近,符合一般道路的路面特征。

    5.2道路可行駛區(qū)域信息檢測與提取

    在對道路橫向邊沿和路面進行提取后,進一步對道路內(nèi)剩余的4層雷達數(shù)據(jù)點提取道路縱向信息,進而得出智能汽車的可行駛區(qū)域。本次試驗分為兩種情形的信息提取,第一是對前方道路阻斷情形檢測(如丁字路口);第二是對道路中凹凸障礙物信息的提取,得到更具體的可行駛區(qū)域。

    對路面行駛區(qū)域的界定,包括橫向和縱向行駛區(qū)域界定。圖5是一個丁字路口情形。利用本文提出的算法提取出兩條道路邊沿線,再分別提取道路寬度、車輛中心距離左右道路邊沿的距離和前方路口距離信息。圖5表明,當前道路寬度為3.71 m,實際測量為3.80 m,誤差為9 cm,車輛距離左側(cè)路沿距離為2.65 m,右側(cè)為1.07 m,當前車輛距離前方丁字路口邊界距離為39.81 m。圖6是一個直道情形,當前道路寬度為10.56 m,實際測量為10.51 m,誤差僅為5 cm,車輛距離左側(cè)路沿距離為5.97m,右側(cè)為4.61 m,當前行駛方向無路口。

    圖5 丁字路口檢測

    由于最高和次高層掃描到較遠地方和入射角度的原因,一般沒有路面回波,而當前方出現(xiàn)障礙物時才會在路面范圍內(nèi)具有一維投影且具有一定的高度值。利用上述特性對凸障礙物進行檢測。圖6檢測到前方有兩個行人、1輛車和1個垃圾桶共4個障礙物,具體信息為{distance,center_x,center_y,width,angle},在圖6的下方顯示,并用方框標識。

    圖6 直路和凸障礙物檢測

    對于凹陷區(qū)域,主要采用提取路面區(qū)域的斷點區(qū)域,利用最低兩層的掃描點數(shù)據(jù)找到凹障礙物,如圖7左側(cè)圖顯示為路面中一個凹陷的坑洼,對應(yīng)圖7右側(cè)圖激光雷達檢測到黑色橢圓所標記區(qū)域,凹障礙物寬度為1.22 m,深度為-0.11 m,坐標位置為(20.23,0.62),距離車體坐標原點為20.24 m,相對車輛前進方向(x軸)夾角為1.75°。

    圖7 凹障礙物檢測

    6 結(jié)束語

    根據(jù)4線激光雷達掃描特性,提出了一種改進的DBSCAN算法。該算法通過一維核密度估計方法對雷達的原始數(shù)據(jù)進行統(tǒng)計分析,對不同維度進行密度估計,再結(jié)合4-dist提取峰值附近的加權(quán)歐式距離值,從而解決了原始DBSCAN聚類方法全局單一參數(shù)問題,能夠根據(jù)數(shù)據(jù)特征自適應(yīng)的對參數(shù)Eps和MinPts進行調(diào)整,同時解決了原聚類算法只能對單一密度區(qū)域進行聚類的問題。根據(jù)道路中不同物體的特性,提出了可行駛區(qū)域信息提取方法,并且在規(guī)則和半規(guī)則化城市道路中的實車試驗中得到驗證。同時還可以根據(jù)得到的可行駛區(qū)域信息進行建立局部地圖,結(jié)合智能汽車全局的路徑規(guī)劃進行導(dǎo)航工作,對道路內(nèi)的障礙物進行識別、跟蹤,進一步完善車輛周圍的局部地圖,也可以作為未來的研究方向。

    參考文獻

    1史鵬波.基于單線激光雷達的道路特征檢測:[學(xué)位論文].南京:南京理工大學(xué), 2013.

    2袁夏,趙春霞.一種激光雷達可通行區(qū)域提取算法.兵工學(xué)報,2010(12): 1702~1707.

    3周俊靜,段建民,楊光祖.基于雷達測距的車輛識別與跟蹤方法.汽車工程,2014(11): 1415~1420.

    4吳維一,劉大學(xué),戴斌.一種處理激光雷達數(shù)據(jù)的聚類分析方法.計算機仿真,2007(08):236~240.

    5萬忠濤.基于激光雷達的道路與障礙檢測研究:[學(xué)位論文].長沙:國防科學(xué)技術(shù)大學(xué), 2010.

    6楊象軍.基于四線激光雷達的道路檢測與跟蹤:[學(xué)位論文].杭州:浙江大學(xué),2013.

    7 Martin Ester, Hans-Peter Kriegel, Jorg Sander,Xiaowei Xu. A density- based algorithm for discovering clusters in large spatial databases with noise.The 2nd International Confer?ence On Knowledge Discovery and Data Ming,1996: 226~231.

    8李雙慶,慕升弟.一種改進的DBSCAN算法及其應(yīng)用.計算機工程與應(yīng)用,2014(08): 72~76.

    9 Rosenblatt M. Remarks on some nonparametric estimates of a density function. Annals of Mathematical Statistic,1956,27 (6): 832~837.

    10ParzenE. On estimation of a probability density function and mode.Annals of Mathematical Statistics,1962,33 (8): 1065~1076.

    11 Christoph Heinzand, Bernhard Seeger. Cluster Kernels: Re?source- Aware Kernel Density Estimators over Streaming Data. IEEE Transactions on Knowledge and Data Engineer?ing, 2008(20):884~885.

    12關(guān)超華,陳泳丹,陳慧巖,龔建偉.基于改進DBSCAN算法的激光雷達車輛探測方法.北京理工大學(xué)學(xué)報,2010 (06): 732~736.

    13錢美旋,葉東毅.利用一維投影分析的無參數(shù)多密度聚類算法.小型微型計算機系統(tǒng),2013(08):1866~1871.

    (責(zé)任編輯簾青)

    修改稿收到日期為2015年9月1日。

    Preceding Drivable Area Detection Based on Four-layer Laser Radar

    Duan Jianmin, Li Longjie,Zheng Kaihua
    (Beijing Key Laboratory of Traffic Engineering, Beijing University of Technology, Beijing 100124)

    【Abstract】By analyzing a large number of laser radar’s data, a modified Density Based Spatial Clustering of Applications with Noise (DBSCAN) algorithm is proposed to obtain the information on roads and obstacles for intelligent vehicle. According to different density parameters, this algorithm is called to complete multiple density clustering. Also, methods for extracting the information on road edge, road surface and obstacles etc., are proposed, which are tested on vehicle on structured and semi- structured urban roads. The results show that the methods can real- time extract the accurate information on drivable area for intelligent vehicle.

    Key words:Four- layer Laser Radar, DBSCAN Algorithm, Weighted Euclidean Distance, Drivable Area

    中圖分類號:U463.6

    文獻標識碼:A

    文章編號:1000-3703(2016)02-0055-07

    *基金項目:北京市屬高等學(xué)校人才強教計劃資助項目(038000543115025)。

    主題詞:4線激光雷達DBSCAN算法加權(quán)歐式距離可行駛區(qū)域

    猜你喜歡
    激光雷達障礙物聚類
    手持激光雷達應(yīng)用解決方案
    北京測繪(2022年5期)2022-11-22 06:57:43
    法雷奧第二代SCALA?激光雷達
    汽車觀察(2021年8期)2021-09-01 10:12:41
    高低翻越
    SelTrac?CBTC系統(tǒng)中非通信障礙物的設(shè)計和處理
    基于激光雷達通信的地面特征識別技術(shù)
    基于激光雷達的多旋翼無人機室內(nèi)定位與避障研究
    電子制作(2018年16期)2018-09-26 03:27:00
    基于DBSACN聚類算法的XML文檔聚類
    電子測試(2017年15期)2017-12-18 07:19:27
    基于改進的遺傳算法的模糊聚類算法
    一種層次初始的聚類個數(shù)自適應(yīng)的聚類方法研究
    自適應(yīng)確定K-means算法的聚類數(shù):以遙感圖像聚類為例
    久久香蕉激情| 黄色丝袜av网址大全| 亚洲成人国产一区在线观看| 韩国精品一区二区三区| 国产一区二区三区视频了| 午夜91福利影院| 亚洲一区高清亚洲精品| 国产成人精品久久二区二区91| 久久国产乱子伦精品免费另类| 757午夜福利合集在线观看| 人妻丰满熟妇av一区二区三区 | 19禁男女啪啪无遮挡网站| 人人妻,人人澡人人爽秒播| 99久久99久久久精品蜜桃| 午夜成年电影在线免费观看| 首页视频小说图片口味搜索| 精品一区二区三区四区五区乱码| 狠狠狠狠99中文字幕| 午夜影院日韩av| 中文字幕av电影在线播放| 欧美成狂野欧美在线观看| 啦啦啦 在线观看视频| 中文字幕制服av| 欧美最黄视频在线播放免费 | 涩涩av久久男人的天堂| 国产精品久久久久成人av| 亚洲免费av在线视频| 久久精品aⅴ一区二区三区四区| 日韩中文字幕欧美一区二区| 一边摸一边抽搐一进一出视频| 亚洲午夜理论影院| 咕卡用的链子| 99久久99久久久精品蜜桃| 欧美乱色亚洲激情| 亚洲一区二区三区不卡视频| av国产精品久久久久影院| 一级片'在线观看视频| 精品久久久久久久毛片微露脸| a级片在线免费高清观看视频| 天天影视国产精品| 精品亚洲成国产av| 久久久久国产精品人妻aⅴ院 | 男人舔女人的私密视频| 亚洲精品一卡2卡三卡4卡5卡| 日日夜夜操网爽| 国产男女内射视频| 99热网站在线观看| 午夜福利在线观看吧| 俄罗斯特黄特色一大片| 久热这里只有精品99| 老熟女久久久| 三上悠亚av全集在线观看| 欧美精品亚洲一区二区| 国产精品久久久久成人av| 午夜成年电影在线免费观看| 中国美女看黄片| 欧美日本中文国产一区发布| av线在线观看网站| 亚洲一区二区三区不卡视频| 美女扒开内裤让男人捅视频| 欧美 亚洲 国产 日韩一| 色婷婷av一区二区三区视频| 满18在线观看网站| 欧美老熟妇乱子伦牲交| 18禁观看日本| 在线十欧美十亚洲十日本专区| 男女高潮啪啪啪动态图| 国产1区2区3区精品| 99久久99久久久精品蜜桃| 757午夜福利合集在线观看| 99国产精品99久久久久| 国产精品 欧美亚洲| 欧美丝袜亚洲另类 | 久热爱精品视频在线9| 久久精品亚洲熟妇少妇任你| 91精品国产国语对白视频| 欧美人与性动交α欧美精品济南到| 满18在线观看网站| 国产伦人伦偷精品视频| 99精国产麻豆久久婷婷| 人妻一区二区av| 精品国产美女av久久久久小说| 亚洲少妇的诱惑av| 欧美日韩亚洲高清精品| 飞空精品影院首页| 久久精品亚洲熟妇少妇任你| 天天躁日日躁夜夜躁夜夜| 91九色精品人成在线观看| 丁香六月欧美| 美女扒开内裤让男人捅视频| 成人三级做爰电影| 99热国产这里只有精品6| 亚洲av欧美aⅴ国产| 啪啪无遮挡十八禁网站| 中亚洲国语对白在线视频| 色精品久久人妻99蜜桃| 一级片免费观看大全| 亚洲视频免费观看视频| 韩国精品一区二区三区| 一级a爱视频在线免费观看| 脱女人内裤的视频| 高清av免费在线| 欧美日韩一级在线毛片| 国产精品久久久久久人妻精品电影| 亚洲熟妇熟女久久| 免费黄频网站在线观看国产| 黄网站色视频无遮挡免费观看| 91在线观看av| 久久ye,这里只有精品| 无人区码免费观看不卡| 视频在线观看一区二区三区| 国产免费现黄频在线看| 丰满的人妻完整版| av天堂在线播放| 日韩中文字幕欧美一区二区| av电影中文网址| 国产99久久九九免费精品| 欧美日韩亚洲国产一区二区在线观看 | www日本在线高清视频| 成年人午夜在线观看视频| 九色亚洲精品在线播放| 超碰97精品在线观看| 新久久久久国产一级毛片| 国产一区在线观看成人免费| 国产蜜桃级精品一区二区三区 | 中文字幕色久视频| 国产男女内射视频| 国产亚洲欧美98| 视频在线观看一区二区三区| 如日韩欧美国产精品一区二区三区| 一区在线观看完整版| 性少妇av在线| 国产精品美女特级片免费视频播放器 | 窝窝影院91人妻| 亚洲欧美日韩高清在线视频| 99久久综合精品五月天人人| 久久久久精品国产欧美久久久| 丰满迷人的少妇在线观看| 免费久久久久久久精品成人欧美视频| 国产精品免费大片| 男女高潮啪啪啪动态图| 黑人欧美特级aaaaaa片| 日日夜夜操网爽| 色老头精品视频在线观看| √禁漫天堂资源中文www| aaaaa片日本免费| 黄片大片在线免费观看| 亚洲欧美一区二区三区黑人| 成人免费观看视频高清| 亚洲男人天堂网一区| 国产主播在线观看一区二区| 精品久久久久久,| 自线自在国产av| 国产精品久久久久久精品古装| 精品福利观看| 国产精品久久久久成人av| 真人做人爱边吃奶动态| 久久影院123| 熟女少妇亚洲综合色aaa.| 极品人妻少妇av视频| 成人黄色视频免费在线看| 丝袜人妻中文字幕| 亚洲色图 男人天堂 中文字幕| 亚洲av欧美aⅴ国产| 午夜福利在线观看吧| 久9热在线精品视频| 在线av久久热| 成人18禁高潮啪啪吃奶动态图| 亚洲欧洲精品一区二区精品久久久| 欧美国产精品一级二级三级| 亚洲色图综合在线观看| 国产成+人综合+亚洲专区| 亚洲午夜精品一区,二区,三区| 变态另类成人亚洲欧美熟女 | 日本欧美视频一区| 国产欧美日韩一区二区三| 91精品三级在线观看| 中文字幕人妻丝袜制服| 视频在线观看一区二区三区| 大型av网站在线播放| 久久久久国内视频| 国产成人av教育| 国产精品美女特级片免费视频播放器 | 久久国产乱子伦精品免费另类| 国产av一区二区精品久久| 99re在线观看精品视频| 久久中文看片网| 91成年电影在线观看| av超薄肉色丝袜交足视频| 久9热在线精品视频| 午夜久久久在线观看| www.精华液| 午夜成年电影在线免费观看| 老司机福利观看| 中出人妻视频一区二区| 国产单亲对白刺激| 久久久国产成人精品二区 | 丝袜人妻中文字幕| 久久精品成人免费网站| 亚洲人成电影观看| 男人的好看免费观看在线视频 | 国产激情欧美一区二区| 啦啦啦视频在线资源免费观看| 丁香六月欧美| 亚洲专区国产一区二区| 99国产极品粉嫩在线观看| 精品国产美女av久久久久小说| 俄罗斯特黄特色一大片| 桃红色精品国产亚洲av| 高清av免费在线| 热re99久久国产66热| 一本综合久久免费| 9热在线视频观看99| 老熟妇仑乱视频hdxx| 在线观看免费视频网站a站| 国产欧美亚洲国产| 1024视频免费在线观看| 亚洲精华国产精华精| 视频在线观看一区二区三区| aaaaa片日本免费| 最近最新中文字幕大全电影3 | 国产伦人伦偷精品视频| 91av网站免费观看| 久久九九热精品免费| 国产97色在线日韩免费| 在线av久久热| 9色porny在线观看| 成年版毛片免费区| 精品一区二区三区视频在线观看免费 | 国产一卡二卡三卡精品| 国产精品免费一区二区三区在线 | 午夜福利在线免费观看网站| 亚洲熟妇中文字幕五十中出 | 美女午夜性视频免费| 两个人免费观看高清视频| 人人妻人人澡人人爽人人夜夜| 91av网站免费观看| 在线观看日韩欧美| av一本久久久久| 免费在线观看视频国产中文字幕亚洲| 国产精品免费视频内射| 在线视频色国产色| 90打野战视频偷拍视频| 男人的好看免费观看在线视频 | 久热爱精品视频在线9| 日韩欧美国产一区二区入口| 亚洲第一青青草原| 69av精品久久久久久| 久久久国产成人精品二区 | 国产欧美日韩一区二区精品| 777米奇影视久久| 国产精品永久免费网站| 美女福利国产在线| 侵犯人妻中文字幕一二三四区| 精品亚洲成国产av| 国产一区二区三区综合在线观看| 狠狠婷婷综合久久久久久88av| 在线观看午夜福利视频| 久9热在线精品视频| www.精华液| 精品国产国语对白av| 国内毛片毛片毛片毛片毛片| 老司机午夜福利在线观看视频| 18禁裸乳无遮挡动漫免费视频| www.自偷自拍.com| 欧美日本中文国产一区发布| 欧美老熟妇乱子伦牲交| 宅男免费午夜| 国产蜜桃级精品一区二区三区 | 91av网站免费观看| 91麻豆精品激情在线观看国产 | 日韩制服丝袜自拍偷拍| 亚洲国产毛片av蜜桃av| av视频免费观看在线观看| 亚洲精品美女久久久久99蜜臀| 久久精品国产清高在天天线| 国产三级黄色录像| 老司机影院毛片| cao死你这个sao货| 亚洲精品国产精品久久久不卡| 国产免费现黄频在线看| 亚洲精品粉嫩美女一区| 成年动漫av网址| 日韩欧美三级三区| av有码第一页| 最近最新中文字幕大全免费视频| a级毛片在线看网站| 一进一出抽搐动态| 女人精品久久久久毛片| 国产精品欧美亚洲77777| av国产精品久久久久影院| 欧美日韩国产mv在线观看视频| 精品国产国语对白av| xxxhd国产人妻xxx| 欧美国产精品va在线观看不卡| 老熟妇仑乱视频hdxx| 国产一区在线观看成人免费| 日韩欧美一区二区三区在线观看 | 18禁美女被吸乳视频| 亚洲av日韩在线播放| 桃红色精品国产亚洲av| 国产av又大| 午夜久久久在线观看| 少妇猛男粗大的猛烈进出视频| 超碰成人久久| 妹子高潮喷水视频| 亚洲av成人不卡在线观看播放网| 天天躁狠狠躁夜夜躁狠狠躁| 国产91精品成人一区二区三区| 90打野战视频偷拍视频| 老鸭窝网址在线观看| 99热国产这里只有精品6| e午夜精品久久久久久久| 成人精品一区二区免费| 老汉色∧v一级毛片| 国产精品影院久久| 99re在线观看精品视频| 少妇被粗大的猛进出69影院| 久久国产精品大桥未久av| 最新的欧美精品一区二区| 老司机午夜十八禁免费视频| 国产精品一区二区在线不卡| 少妇猛男粗大的猛烈进出视频| 美女国产高潮福利片在线看| 午夜精品国产一区二区电影| 国产精品香港三级国产av潘金莲| 久久久久精品人妻al黑| 久久久久久久久久久久大奶| 色婷婷av一区二区三区视频| 多毛熟女@视频| 美女福利国产在线| 亚洲成人免费电影在线观看| 欧美av亚洲av综合av国产av| 99精品欧美一区二区三区四区| 91av网站免费观看| 一进一出抽搐动态| 伦理电影免费视频| 亚洲av成人一区二区三| 亚洲成人免费电影在线观看| 午夜福利欧美成人| 一个人免费在线观看的高清视频| 1024视频免费在线观看| 色婷婷av一区二区三区视频| 窝窝影院91人妻| 一级黄色大片毛片| 亚洲 国产 在线| 99re在线观看精品视频| 国产精品.久久久| 又大又爽又粗| 久久久国产成人精品二区 | 香蕉国产在线看| 亚洲av熟女| 视频区欧美日本亚洲| 18在线观看网站| 久久久久久亚洲精品国产蜜桃av| 久久久久久久久免费视频了| av不卡在线播放| 久久精品国产a三级三级三级| 欧美日本中文国产一区发布| 久久精品亚洲av国产电影网| 亚洲国产中文字幕在线视频| 亚洲成人手机| 日韩中文字幕欧美一区二区| 国产又色又爽无遮挡免费看| 夜夜爽天天搞| 欧美日韩亚洲国产一区二区在线观看 | 国产不卡av网站在线观看| 久久人妻熟女aⅴ| 精品无人区乱码1区二区| 热re99久久国产66热| 午夜精品国产一区二区电影| 日本精品一区二区三区蜜桃| 亚洲性夜色夜夜综合| 国产av又大| 国产日韩一区二区三区精品不卡| 在线观看舔阴道视频| av中文乱码字幕在线| 人人妻人人澡人人爽人人夜夜| 免费不卡黄色视频| 一级毛片精品| 香蕉久久夜色| 亚洲伊人色综图| 久久精品人人爽人人爽视色| 中国美女看黄片| 女人被狂操c到高潮| 亚洲av片天天在线观看| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品美女久久av网站| 欧美日韩成人在线一区二区| 熟女少妇亚洲综合色aaa.| 在线十欧美十亚洲十日本专区| 涩涩av久久男人的天堂| 国产午夜精品久久久久久| 电影成人av| 国产男女内射视频| 久久精品亚洲精品国产色婷小说| 国产成人系列免费观看| 欧美激情高清一区二区三区| aaaaa片日本免费| 国产成人精品久久二区二区免费| 麻豆乱淫一区二区| 在线观看免费视频网站a站| 在线永久观看黄色视频| 91麻豆精品激情在线观看国产 | 91字幕亚洲| 午夜免费成人在线视频| 飞空精品影院首页| 曰老女人黄片| 日韩有码中文字幕| 免费观看a级毛片全部| 每晚都被弄得嗷嗷叫到高潮| 亚洲五月色婷婷综合| 日本黄色视频三级网站网址 | 99re6热这里在线精品视频| 深夜精品福利| 亚洲aⅴ乱码一区二区在线播放 | 亚洲av片天天在线观看| 精品人妻在线不人妻| av福利片在线| 精品久久久精品久久久| a级片在线免费高清观看视频| 国产精品秋霞免费鲁丝片| 欧美日韩精品网址| 久久久水蜜桃国产精品网| 大码成人一级视频| 国产亚洲一区二区精品| 色94色欧美一区二区| www.熟女人妻精品国产| 国产深夜福利视频在线观看| 亚洲av第一区精品v没综合| 日日夜夜操网爽| 日韩欧美一区视频在线观看| 欧美精品一区二区免费开放| 亚洲情色 制服丝袜| 日本五十路高清| 日韩人妻精品一区2区三区| 中文字幕高清在线视频| 色在线成人网| 亚洲黑人精品在线| 99久久国产精品久久久| 国产精品久久久人人做人人爽| 久久香蕉国产精品| 大香蕉久久网| 黑人猛操日本美女一级片| 亚洲精品国产区一区二| 18禁美女被吸乳视频| 淫妇啪啪啪对白视频| 欧美在线黄色| 亚洲色图 男人天堂 中文字幕| av电影中文网址| a级片在线免费高清观看视频| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品亚洲精品国产色婷小说| 亚洲熟妇熟女久久| 欧美另类亚洲清纯唯美| 午夜影院日韩av| 久久午夜亚洲精品久久| 午夜福利影视在线免费观看| 亚洲熟妇中文字幕五十中出 | 日韩成人在线观看一区二区三区| 国产精华一区二区三区| 欧美性长视频在线观看| 十八禁网站免费在线| 亚洲第一欧美日韩一区二区三区| 老熟妇仑乱视频hdxx| 黄片大片在线免费观看| 欧美性长视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 两性夫妻黄色片| 欧美日韩av久久| 成人免费观看视频高清| 久久久久精品国产欧美久久久| 国产成人系列免费观看| 免费在线观看黄色视频的| 久久国产精品影院| 久久久国产精品麻豆| 亚洲在线自拍视频| 99re在线观看精品视频| 亚洲成a人片在线一区二区| 丰满饥渴人妻一区二区三| 久久人妻熟女aⅴ| 亚洲午夜理论影院| 日本黄色日本黄色录像| av一本久久久久| 精品亚洲成国产av| 男人操女人黄网站| av免费在线观看网站| 90打野战视频偷拍视频| 脱女人内裤的视频| 宅男免费午夜| 亚洲第一欧美日韩一区二区三区| 嫩草影视91久久| 18在线观看网站| 757午夜福利合集在线观看| 久久精品人人爽人人爽视色| 亚洲精品久久午夜乱码| 久久精品国产亚洲av香蕉五月 | 亚洲精品国产区一区二| 一边摸一边抽搐一进一小说 | 久久人妻熟女aⅴ| 一级毛片高清免费大全| 久久中文看片网| 伦理电影免费视频| 午夜福利乱码中文字幕| 可以免费在线观看a视频的电影网站| 女人被躁到高潮嗷嗷叫费观| 超色免费av| 一级黄色大片毛片| 麻豆av在线久日| 下体分泌物呈黄色| 国产成人精品久久二区二区91| 国产成人精品久久二区二区免费| www.999成人在线观看| 正在播放国产对白刺激| 村上凉子中文字幕在线| 性少妇av在线| 成人国产一区最新在线观看| 如日韩欧美国产精品一区二区三区| 午夜两性在线视频| 亚洲五月天丁香| 超碰成人久久| av有码第一页| 国产成人免费无遮挡视频| 亚洲成av片中文字幕在线观看| 国产高清国产精品国产三级| 亚洲第一av免费看| 国产麻豆69| 成年人午夜在线观看视频| 免费观看人在逋| 在线观看舔阴道视频| 90打野战视频偷拍视频| 亚洲av片天天在线观看| 国产不卡一卡二| 69av精品久久久久久| 高清av免费在线| 欧美激情久久久久久爽电影 | 欧美 日韩 精品 国产| 国产av一区二区精品久久| 久久久久久久久久久久大奶| xxxhd国产人妻xxx| 水蜜桃什么品种好| 精品亚洲成国产av| 黄色女人牲交| 久久久久久久久久久久大奶| 女人被狂操c到高潮| 老司机亚洲免费影院| 久久精品熟女亚洲av麻豆精品| 亚洲国产欧美一区二区综合| 色综合欧美亚洲国产小说| 岛国在线观看网站| 国产成人系列免费观看| 少妇裸体淫交视频免费看高清 | 亚洲伊人色综图| 久久精品人人爽人人爽视色| 视频区图区小说| 在线观看舔阴道视频| 午夜影院日韩av| 国产成人影院久久av| 久久影院123| 999久久久国产精品视频| 黄色视频,在线免费观看| 夜夜夜夜夜久久久久| 99国产精品一区二区三区| 99re在线观看精品视频| 最近最新免费中文字幕在线| 91麻豆av在线| 精品少妇一区二区三区视频日本电影| 丰满的人妻完整版| 久久精品亚洲av国产电影网| 成人精品一区二区免费| 午夜免费成人在线视频| cao死你这个sao货| 人妻丰满熟妇av一区二区三区 | 亚洲精品一卡2卡三卡4卡5卡| 一区二区日韩欧美中文字幕| 免费女性裸体啪啪无遮挡网站| 免费在线观看亚洲国产| 色老头精品视频在线观看| 涩涩av久久男人的天堂| 桃红色精品国产亚洲av| 制服人妻中文乱码| 一级a爱视频在线免费观看| 成年动漫av网址| 亚洲一区二区三区欧美精品| 色播在线永久视频| 一区二区日韩欧美中文字幕| 18禁裸乳无遮挡免费网站照片 | 老司机深夜福利视频在线观看| 黑人巨大精品欧美一区二区蜜桃| 老司机影院毛片| 久久九九热精品免费| 免费不卡黄色视频| 亚洲avbb在线观看| xxxhd国产人妻xxx| 99久久99久久久精品蜜桃| 777久久人妻少妇嫩草av网站| 日韩中文字幕欧美一区二区| 午夜影院日韩av| 久久久久精品人妻al黑| 一夜夜www| 一区二区三区精品91| 国产av又大| 三级毛片av免费| √禁漫天堂资源中文www| 脱女人内裤的视频| 欧美日韩精品网址| 亚洲 国产 在线| 国产蜜桃级精品一区二区三区 | 黄色视频不卡| 久久99一区二区三区| 久9热在线精品视频| 国产成人免费观看mmmm| 亚洲在线自拍视频| 一本大道久久a久久精品| 精品福利观看|