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

    基于局部熵的點(diǎn)云精簡(jiǎn)算法

    2020-03-03 13:20:44董嘉敏田華
    現(xiàn)代電子技術(shù) 2020年1期
    關(guān)鍵詞:點(diǎn)云數(shù)據(jù)means算法信息熵

    董嘉敏 田華

    摘 ?要: 針對(duì)目前流行的三維物體激光掃描儀獲取的點(diǎn)云數(shù)據(jù)量大,冗余度高等問(wèn)題,提出一種基于信息熵的點(diǎn)云精簡(jiǎn)算法。首先,定義數(shù)據(jù)點(diǎn)的曲率、點(diǎn)到鄰域點(diǎn)重心的距離、點(diǎn)到鄰域點(diǎn)的平均距離的倒數(shù),三者乘積為權(quán)值積;然后,使用K?means聚類算法劃分點(diǎn)云數(shù)據(jù),根據(jù)類內(nèi)估計(jì)曲率差值區(qū)分特征區(qū)域與非特征區(qū)域;最后,針對(duì)特征區(qū)域,利用提出的精簡(jiǎn)方法精簡(jiǎn)點(diǎn)云。實(shí)驗(yàn)結(jié)果表明,該方法計(jì)算相對(duì)簡(jiǎn)單,能夠有效避免孔洞現(xiàn)象,同時(shí),更好地保留了點(diǎn)云數(shù)據(jù)的原始物理特征。

    關(guān)鍵詞: 點(diǎn)云精簡(jiǎn); 信息熵; K?means算法; 特征區(qū)域區(qū)分; 點(diǎn)云數(shù)據(jù); 曲率估計(jì)

    中圖分類號(hào): TN911?34; TP391.9 ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)識(shí)碼: A ? ? ? ? ? ? ? ? ? ?文章編號(hào): 1004?373X(2020)01?0020?04

    A point cloud reduction algorithm based on local entropy

    DONG Jiamin, TIAN Hua

    Abstract: A point cloud reduction algorithm based on information entropy is proposed to deal with the large amount of point cloud data and high redundancy produced by the popular three?dimensional laser scanner. Firstly, the product of curvature of the data points, the distance from the points to the center of gravity of the neighborhood points, and the reciprocal of the average distance from the points to the neighborhood points is defined as the product of weight; secondly, the K?means clustering algorithm is used to classify the point cloud data, and the characteristic region and the non?characteristic region are distinguished according to the estimated curvature difference; finally, for the characteristic region, the proposed reduction algorithm is used to reduce the point cloud. The experimental results show that the calculation process of the proposed algorithm is relatively simple, which can effectively avoid the hole phenomenon, and preserve the original physical characteristics of the point cloud data better.

    Keywords: point cloud reduction; information entropy; K?means algorithm; characteristic region distinction; point cloud data; curvature estimation

    0 ?引 ?言

    三維重建技術(shù)一直是計(jì)算機(jī)視覺(jué)、逆向工程學(xué)和計(jì)算機(jī)圖形學(xué)等熱門研究領(lǐng)域的研究重點(diǎn)[1]。關(guān)于點(diǎn)云技術(shù)三維重建,已經(jīng)有不少的研究工作和成果[2?3]。隨著FreeScanX3等三維物體激光掃描儀不斷推陳出新,點(diǎn)云數(shù)據(jù)的獲得變得更加方便,同時(shí)點(diǎn)云數(shù)據(jù)的精度也非常高。然而,如此密集的點(diǎn)云數(shù)據(jù)中卻存在著大量的冗余數(shù)據(jù),這些冗余數(shù)據(jù)極大地增加了點(diǎn)云數(shù)據(jù)的存儲(chǔ)開銷和計(jì)算開銷,從而影響到曲面擬合和模型生成等重要后續(xù)工作的效率。綜上所述,點(diǎn)云數(shù)據(jù)的精簡(jiǎn)工作對(duì)于物體的三維重建和繪制工作是非常重要的[4]。

    近年來(lái),許多學(xué)者對(duì)點(diǎn)云精簡(jiǎn)的研究做出了貢獻(xiàn)并取得了很大的成功。根據(jù)對(duì)點(diǎn)云數(shù)據(jù)精簡(jiǎn)的原則,現(xiàn)有的點(diǎn)云精簡(jiǎn)方法可以分為兩類:基于網(wǎng)格的簡(jiǎn)化和基于點(diǎn)的簡(jiǎn)化?;诰W(wǎng)格的簡(jiǎn)化方法主要通過(guò)構(gòu)造多邊形網(wǎng)格(三角形網(wǎng)格或四邊形網(wǎng)格),依據(jù)一些規(guī)則來(lái)判定這些網(wǎng)格的重要性,刪除冗余網(wǎng)格的邊來(lái)實(shí)現(xiàn)精簡(jiǎn)。這些方法的不足之處在于網(wǎng)格的生成需要花費(fèi)大量的時(shí)間和存儲(chǔ)。相比較這類方法而言,基于點(diǎn)的簡(jiǎn)化方法直接對(duì)點(diǎn)云數(shù)據(jù)進(jìn)行精簡(jiǎn)而不用構(gòu)建任何網(wǎng)格。文獻(xiàn)[5]采用K均值聚類算法在空間域中將相似點(diǎn)聚集在一起,并使用最大法向量偏差作為聚類分散的度量,將聚集的點(diǎn)集劃分為特征域中的一系列子聚類,提出一種新的劃分和細(xì)分方法,實(shí)現(xiàn)對(duì)點(diǎn)云的精簡(jiǎn);文獻(xiàn)[6]針對(duì)散亂點(diǎn)云簡(jiǎn)化中易丟失幾何特征及潛在曲面形狀信息的問(wèn)題,采用基于泊松分布的區(qū)域生長(zhǎng)法自適應(yīng)檢測(cè)特征點(diǎn),通過(guò)設(shè)定不同的聚類閾值,運(yùn)用不同的簡(jiǎn)化策略簡(jiǎn)化點(diǎn)云數(shù)據(jù);文獻(xiàn)[7]運(yùn)用多判別參數(shù)混合的方法首先識(shí)別特征點(diǎn)并保留,然后對(duì)剩余點(diǎn)云數(shù)據(jù)K?means聚類,根據(jù)類內(nèi)最大曲率差作為細(xì)分標(biāo)準(zhǔn),最終完成特征保持的點(diǎn)云數(shù)據(jù)精簡(jiǎn)。

    這些不同的基于點(diǎn)的簡(jiǎn)化方法有一個(gè)共同點(diǎn),就是先評(píng)估點(diǎn)云數(shù)據(jù)中每個(gè)點(diǎn)的重要性,通過(guò)刪除不重要的點(diǎn)完成點(diǎn)云數(shù)據(jù)的精簡(jiǎn),因此評(píng)估參數(shù)就顯得尤為重要。本文提出一種基于信息熵的點(diǎn)云精簡(jiǎn)算法,針對(duì)以往最小二乘法擬合曲面計(jì)算曲率過(guò)于復(fù)雜耗時(shí),提出一種新的曲率估計(jì)方法,并定義數(shù)據(jù)點(diǎn)的曲率、點(diǎn)到鄰域點(diǎn)重心的距離、點(diǎn)到鄰域點(diǎn)的平均距離的倒數(shù),三者乘積為權(quán)值積,以權(quán)值積衡量點(diǎn)的重要性,降低曲率估計(jì)誤差影響,然后結(jié)合K?means聚類算法和權(quán)值積的信息熵保留特征點(diǎn),完成點(diǎn)云數(shù)據(jù)的精簡(jiǎn)。

    1 ?權(quán)值積的計(jì)算

    定義曲率[H]、點(diǎn)到鄰域點(diǎn)重心的距離[S]、點(diǎn)到鄰域點(diǎn)的平均距離的倒數(shù)[1A]的乘積為權(quán)值積[w],即:

    [w=H?S?1A] (1)

    式中:[1A]與[S]可通過(guò)歐氏距離的簡(jiǎn)單計(jì)算獲得,此處不再贅述。

    眾所周知,曲率是曲面的基本信息,反映曲面的局部幾何特征。因此根據(jù)曲率的大小可以反映面的彎曲程度。在三維歐氏空間中,一般使用高斯曲率來(lái)反映某點(diǎn)的彎曲程度。一般實(shí)現(xiàn)方法是將此點(diǎn)和鄰域內(nèi)的點(diǎn)擬合成曲面,再求取曲面的主曲率和高斯曲率。設(shè)[P]點(diǎn)為點(diǎn)云中的一點(diǎn),使用K?D樹算法尋得[P]點(diǎn)的8近鄰[(P1,P2,…,P8)]。

    1.1 ?傳統(tǒng)最小二乘法求解曲率[8]

    建立二次曲面的參數(shù)方程:

    [ru,v=i,j=02qijuivj] (2)

    令:

    [ ?a=[a00 a01 a02 a10 a11 a12 a20 a21 a22] ?b=[b00 ?b01 ?b02 ?b10 ?b11 ?b12 ?b20 ?b21 ?b22] ?c=[c00 ?c01 ?c02 ?c10 ?c11 ?c12 ?c20 ?c21 ?c22] ?l=[u0v0 ?u0v1 ?u0v2 ?u1v0 ?u1v1 ?u1v2 ?u2v0 ?u2v1 ?u2v2] ?Q=[a ?b ?c]] (3)

    式(2)可改寫為:

    [r(u,v)=[x(u,v)y(u,v)z(u,v)]=[lTalTblTc]=lTQ] ? (4)

    引入兩個(gè)矩陣:

    [B=x0y0z0x1?y1?z1?xkykzk, ? ? ? l=lT0lT1?lTk] (5)

    待測(cè)點(diǎn)[P]的坐標(biāo)為[(x0,y0,z0)],這里[k=8]。逼近曲面和數(shù)據(jù)點(diǎn)的誤差函數(shù)為:

    [Ω=lQ-B] ? (6)

    運(yùn)用最小二乘法,推導(dǎo)出系數(shù)矩陣[Q]為:

    [Q=(lTl)-1lTB] (7)

    從而求出曲面[r(u,v)]的偏導(dǎo)數(shù)[ru],[rv],[ruu],[ruv],[rvv],曲面的單位法矢[n]為:

    [n=ru×rvru×rv] (8)

    曲面的第1、第2基本量為:

    [E=r2u,F(xiàn)=rurv,G=r2,L=nruu,M=nruv,N=nrvv] (9)

    求出高斯曲率[K]和平均曲率[H]為:

    [K=LN×l2EG-F2, ? ?H=LG-2MF+NE2(EG-F2)] (10)

    由于三維掃描的數(shù)據(jù)量往往很大,因此對(duì)每個(gè)點(diǎn)擬合曲面將會(huì)耗費(fèi)大量的內(nèi)存和時(shí)間,所以本文根據(jù)點(diǎn)云間的三角拓?fù)潢P(guān)系提出一種新的曲率估計(jì)方法。

    1.2 ?新的曲率估計(jì)方法

    如圖2所示,首先找到待測(cè)點(diǎn)[P]的8近鄰[P1~P8],根據(jù)這9個(gè)點(diǎn)最小二乘擬合確定平面[S],向量[α]表示平面[S]的法向量。估計(jì)曲率的計(jì)算過(guò)程如下:

    1) 首先連接[P]與其余8個(gè)點(diǎn),分別得到向量:[PP1],[PP2],[PP3],[PP4],[PP5],[PP6],[PP7]和[PP8]。

    2) 向量[α]與向量[PP1]的余弦值為:

    [cosθ1=α×PP1α·PP1] ?(11)

    圖1中,實(shí)心點(diǎn)為待測(cè)點(diǎn),空心點(diǎn)為鄰近點(diǎn)。角[θ]的余弦值的大小反映了兩個(gè)向量夾角的大小,由圖1可知,待測(cè)點(diǎn)附近的凹凸情況可由余弦值的絕對(duì)值來(lái)說(shuō)明,若余弦值的絕對(duì)值越大,說(shuō)明兩向量的夾角越大或越小,而這兩種情況都說(shuō)明[P]點(diǎn)附近的曲面凹凸程度越大,是特征區(qū)域。根據(jù)式(11),同樣可得[cos θ2],[cos θ3],[cosθ4],[cos θ5],[cos θ6],[cos θ7]和[cos θ8]。令曲率估計(jì)值為:

    [H=i=18cos θi] (12)

    [H]的大小反映了點(diǎn)云數(shù)據(jù)在待測(cè)點(diǎn)[P]附近的凹凸程度。[H]越大,表明在[P]點(diǎn)附近凹凸程度越大,[H]越小表示在[P]點(diǎn)附近的凹凸程度越小。權(quán)值積的另兩個(gè)參數(shù),[S]越大說(shuō)明[P]點(diǎn)附近幾何特征越明顯;[1A]越大說(shuō)明[P]點(diǎn)附近點(diǎn)云密度較為稀疏,應(yīng)該保存較多的點(diǎn)以保護(hù)原始點(diǎn)云特征。

    完整的待測(cè)點(diǎn)曲率估計(jì)示意圖如圖2所示。

    2 ?本文算法

    2.1 ?K?means聚類算法

    K?means聚類算法是典型的基于歐氏距離的聚類算法。在點(diǎn)云的分割過(guò)程中,能夠保持穩(wěn)定和快速的優(yōu)點(diǎn)。若原始點(diǎn)集是位于[N]維空間中的點(diǎn)集,表示為[Pi],[i∈{1,2,…,m}]。首先找到[K]個(gè)聚類中心[{C1,C2,…,Ck}],然后通過(guò)迭代計(jì)算移動(dòng)[K]個(gè)聚類中心,最終使得[D]值最?。?/p>

    [D=argmini=1kj=1uPj-Ci] (13)

    式中:[u]表示以[Ci]為聚類中心的簇中的點(diǎn)的個(gè)數(shù);[Pj]表示屬于以[Ci]為聚類中心的簇中的點(diǎn)。由于[K?D]樹方法的特殊性,點(diǎn)云依照[K?D]樹劃分后會(huì)非常均勻,即可以由[K?D]樹的葉節(jié)點(diǎn)決定聚類重心,該算法計(jì)算過(guò)程如下:

    1) 選定聚類[K]值的大小,[K?D]樹劃分點(diǎn)云數(shù)據(jù),直到葉節(jié)點(diǎn)的數(shù)目最接近[K]時(shí),初始的聚類中心選擇[K?D]樹的葉節(jié)點(diǎn);

    2) 根據(jù)點(diǎn)到聚類中心的歐氏距離將點(diǎn)云劃分為[K]個(gè)簇;

    3) 計(jì)算各個(gè)簇的重心代替原來(lái)簇中心;

    4) 按照新的聚類中心重新聚類;

    5) 重復(fù)步驟3)和步驟4),直到簇的中心不再發(fā)生變化,得到[K]個(gè)簇。

    2.2 ?信息熵理論

    信息熵是一種非常重要的系統(tǒng)科學(xué)理論,是系統(tǒng)的一個(gè)狀態(tài)函數(shù),其含義非常豐富,最初由著名科學(xué)家香農(nóng)(Shannon)于1948年在其《A Mathematical Theory of Communication》一文中提出,又稱為Shannon熵,用來(lái)度量通信中的不確定性和信息量。Shannon給出了信息熵的定量表述[9]:

    [H(X)=-Ki=1nPilog Pi] ?(14)

    式中:[K]為某個(gè)固定正常數(shù);[Pi]是事件[Xi]出現(xiàn)的概率,滿足條件[0≤Pi(i=1,2,…,n)]和[i=1nPi=1]。

    根據(jù)信息熵理論,本文將信息熵引入點(diǎn)云精簡(jiǎn)算法之中,在數(shù)據(jù)點(diǎn)權(quán)值積的特征表達(dá)上,信息熵理論可用于將某一數(shù)據(jù)點(diǎn)權(quán)值積表征為對(duì)某類特征的隸屬程度。對(duì)于點(diǎn)云中一點(diǎn)[P],它的信息熵可由下式計(jì)算:

    [I(P)=-Pwlog Pw-i=1kPwi] ?(15)

    其中:

    [Pw=ww+i=1kwi, ?Pwi=wiw+i=1kwi] (16)

    式中:[Pw]是點(diǎn)[P]的權(quán)值積概率;[Pwi]是第[i]個(gè)相鄰點(diǎn)的權(quán)值積概率。較大[I(P)]值意味著該點(diǎn)位于凹陷和凸起明顯變化的區(qū)域中,如果移除此點(diǎn),則局部幾何圖形將發(fā)生變化,因此應(yīng)該保留信息熵較大值的點(diǎn)云數(shù)據(jù)。

    2.3 ?保持特征的精簡(jiǎn)算法流程

    設(shè)聚類后得到[n]個(gè)簇([S1]~[Sn]),針對(duì)每一個(gè)簇[Si],給定權(quán)值[λ]用來(lái)定義點(diǎn)云特征保持度,將權(quán)值積之差與[λ]比較,若[λ]越小則特征保持效果越好,[λ]越大則特征保持越小但精簡(jiǎn)程度更高。算法步驟流程如圖3所示,具體步驟如下:

    1) 計(jì)算簇中每個(gè)點(diǎn)的權(quán)值積[wi],其中,估計(jì)曲率最大值為[Hmax],最小值為[Hmin];

    2) 若簇中[Hmax-Hmin<λ],則認(rèn)為此簇處于比較平坦的區(qū)域,即用簇中心代替此簇;

    3) 若簇中[Hmax-Hmin>λ],則判斷此簇為特征區(qū)域,計(jì)算此簇中各點(diǎn)權(quán)值積的信息熵,求得平均值,信息熵大于平均值的點(diǎn)予以保留。

    通過(guò)以上算法即可實(shí)現(xiàn)對(duì)點(diǎn)云數(shù)據(jù)的精簡(jiǎn),相對(duì)平滑區(qū)域,使用簇中心替代整個(gè)簇,不會(huì)造成幾何特征的明顯變化;特征區(qū)域由于使用信息熵值來(lái)衡量點(diǎn)的重要性,保留了最重要的點(diǎn),因此保留了原始點(diǎn)云的幾何特征。

    3 ?實(shí)驗(yàn)結(jié)果及分析

    本文使用Matlab對(duì)算法進(jìn)行實(shí)現(xiàn)。實(shí)驗(yàn)用到的點(diǎn)云數(shù)據(jù)為斯坦福大學(xué)建立的3D點(diǎn)云數(shù)據(jù)庫(kù)中的斯坦福兔子Buuny和Dragon點(diǎn)云數(shù)據(jù)模型,格式均為PLY。本文算法有兩個(gè)變量需要自行確定,初始聚類數(shù)目[k]和權(quán)值[λ]。本次實(shí)驗(yàn)取初始聚類數(shù)[k]為2 048,[λ]分別取0.6,0.1。

    圖4所示為原始點(diǎn)云由Matlab的呈現(xiàn)結(jié)果,其中圖4a)為Buuny模型,共有35 947個(gè)點(diǎn),圖4b)為有10 248個(gè)點(diǎn)的Dragon模型。為了對(duì)點(diǎn)云原始數(shù)據(jù)和精簡(jiǎn)后數(shù)據(jù)進(jìn)行定量分析,本文采用Geomagic Studio軟件實(shí)現(xiàn)對(duì)點(diǎn)云數(shù)據(jù)的重構(gòu)[10],依此衡量精簡(jiǎn)結(jié)果的優(yōu)劣。圖5為原始點(diǎn)云數(shù)據(jù)的重構(gòu)結(jié)果。

    圖6a)和圖7a)為[λ=]0.1時(shí)的精簡(jiǎn)結(jié)果,圖6b)和圖7b)為[λ=]0.6時(shí)的精簡(jiǎn)結(jié)果。Buuny模型精簡(jiǎn)后點(diǎn)的個(gè)數(shù)為7 271,3 571,Dragon模型精簡(jiǎn)后點(diǎn)的個(gè)數(shù)為177 416,8 654。

    精簡(jiǎn)后的模型重建結(jié)果如圖8所示。與原始模型進(jìn)行比較,在[λ=]0.1時(shí)可以較好地保留其特征信息;在[λ=]0.6時(shí),由于精簡(jiǎn)度過(guò)高,會(huì)丟失部分特征信息。

    4 ?結(jié) ?語(yǔ)

    針對(duì)最小二乘法擬合曲面計(jì)算數(shù)據(jù)點(diǎn)曲率過(guò)于復(fù)雜和費(fèi)時(shí),本文提出一種新的點(diǎn)云數(shù)據(jù)點(diǎn)曲率估計(jì)方法,并將信息熵的理論引入點(diǎn)云精簡(jiǎn)的過(guò)程中,更好地保留了特征點(diǎn)。經(jīng)過(guò)實(shí)驗(yàn)證明,該算法精簡(jiǎn)效果較好,保留了原始數(shù)據(jù)的基本幾何特征,精簡(jiǎn)后的點(diǎn)云數(shù)據(jù)用于重建,效果良好。

    注:本文通訊作者為田華。

    參考文獻(xiàn)

    [1] THANOU D, CHOU P A, FROSSARD P. Graph?based compression of dynamic 3D point cloud sequences [J]. IEEE tran?sactions on image processing, 2016, 25(4): 1765?1778.

    [2] M BERGER, A TAGLIASACCHI, L SEVERSKY, et al. State of the art in surface reconstruction from point cloud [J]. Eurographics star reports, 2014, 1(1): 161?185.

    [3] NAGAI Y, OHTAKE Y, SUZUKI H, et al. Tomographic surface reconstruction from point cloud [J]. Computers & graphics, 2015, 46: 55?63.

    [4] SANCHEZ G, LEAL E, LEAL N. A linear programming approach for 3D point cloud simplification [J]. Iaeng international journal of computer science, 2017, 44(1): 60?67.

    [5] SHI B Q, LIANG J, LIU Q. Adaptive simplification of point cloud using k?means clustering [J]. Computer?aided design, 2011, 43(8): 910?922.

    [6] ZHANG Y, GENG G, WEI X, et al. A statistical approach for extraction of feature lines from point clouds [J]. Computers & graphics, 2016, 56: 31?45.

    [7] 陳龍,蔡勇,張建生.自適應(yīng)K?means聚類的散亂點(diǎn)云精簡(jiǎn)[J].中國(guó)圖象圖形學(xué)報(bào),2017,22(8):1089?1097.

    [8] 周煜,張萬(wàn)兵,杜發(fā)榮,等.散亂點(diǎn)云數(shù)據(jù)的曲率精簡(jiǎn)算法[J].北京理工大學(xué)學(xué)報(bào),2010,30(7):785?789.

    [9] 姜茸,廖鴻志,楊明.信息熵在軟件領(lǐng)域中的應(yīng)用研究現(xiàn)狀[J].自動(dòng)化技術(shù)與應(yīng)用,2015,34(4):1?6.

    [10] 郭培閃,杜黎明.運(yùn)用Geomagic Studio實(shí)現(xiàn)點(diǎn)云數(shù)據(jù)的曲面重建及誤差分析[J].地理信息世界,2015,22(1):57?60.

    作者簡(jiǎn)介:董嘉敏(1993—),男,山西運(yùn)城人,碩士研究生,研究方向?yàn)樘摂M現(xiàn)實(shí)。

    田 ?華(1983—),女,山西晉城人,博士,碩士生導(dǎo)師,主要研究方向?yàn)橛?jì)算材料。

    猜你喜歡
    點(diǎn)云數(shù)據(jù)means算法信息熵
    基于信息熵可信度的測(cè)試點(diǎn)選擇方法研究
    基于信息熵的實(shí)驗(yàn)教學(xué)量化研究
    一種基于信息熵的雷達(dá)動(dòng)態(tài)自適應(yīng)選擇跟蹤方法
    校園典型建筑物點(diǎn)云數(shù)據(jù)三維建模方法
    三維激光掃描技術(shù)在異型建筑玻璃幕墻中的應(yīng)用研究
    一種改進(jìn)的點(diǎn)云數(shù)據(jù)組合精簡(jiǎn)算法
    基于K—Means聚類算法入侵檢測(cè)系統(tǒng)研究
    基于Weka的Apriori算法在原油產(chǎn)量預(yù)測(cè)中的應(yīng)用
    基于HSI顏色空間的小麥粉精度自動(dòng)識(shí)別研究
    基于聚類的Web日志挖掘
    亚洲欧美日韩高清专用| 90打野战视频偷拍视频| 欧美最新免费一区二区三区 | 高潮久久久久久久久久久不卡| 亚洲av成人精品一区久久| 久久久久久久久中文| 国产欧美日韩一区二区三| 男女午夜视频在线观看| 成人av一区二区三区在线看| 少妇的逼好多水| 一个人免费在线观看电影| 中国美女看黄片| 青草久久国产| 特大巨黑吊av在线直播| 丁香六月欧美| 亚洲中文字幕日韩| 两性午夜刺激爽爽歪歪视频在线观看| 啦啦啦免费观看视频1| 美女黄网站色视频| 久久久精品欧美日韩精品| 麻豆国产av国片精品| 国产麻豆成人av免费视频| 91麻豆av在线| 最新美女视频免费是黄的| 精品福利观看| 男人舔奶头视频| 久久精品国产自在天天线| 亚洲狠狠婷婷综合久久图片| 最后的刺客免费高清国语| 欧美日韩福利视频一区二区| 色尼玛亚洲综合影院| 99久久九九国产精品国产免费| 亚洲一区二区三区色噜噜| 亚洲自拍偷在线| 国产私拍福利视频在线观看| 日日干狠狠操夜夜爽| 日本一二三区视频观看| 老司机午夜福利在线观看视频| 中文字幕av在线有码专区| 99久久九九国产精品国产免费| 12—13女人毛片做爰片一| 精华霜和精华液先用哪个| 俄罗斯特黄特色一大片| 久久亚洲真实| 在线看三级毛片| 国产精品亚洲av一区麻豆| 9191精品国产免费久久| 人妻夜夜爽99麻豆av| eeuss影院久久| а√天堂www在线а√下载| 男插女下体视频免费在线播放| 免费在线观看亚洲国产| 高清在线国产一区| 亚洲成人免费电影在线观看| 亚洲精品日韩av片在线观看 | 美女cb高潮喷水在线观看| 在线观看免费午夜福利视频| 女人被狂操c到高潮| 国内精品久久久久久久电影| 欧美乱色亚洲激情| 欧美一区二区精品小视频在线| 久久久久久久久久黄片| 天天躁日日操中文字幕| 国产视频内射| 国产精品一区二区免费欧美| 三级国产精品欧美在线观看| 久久婷婷人人爽人人干人人爱| 国产一区二区亚洲精品在线观看| 天堂动漫精品| 桃红色精品国产亚洲av| 黄色丝袜av网址大全| 免费观看精品视频网站| 嫩草影视91久久| 午夜a级毛片| 国产一区二区激情短视频| 精品熟女少妇八av免费久了| 国产精品一区二区免费欧美| 夜夜看夜夜爽夜夜摸| 欧美又色又爽又黄视频| 成人午夜高清在线视频| 老司机午夜十八禁免费视频| 国产色婷婷99| 亚洲精品在线观看二区| 国产精品免费一区二区三区在线| 一夜夜www| 免费搜索国产男女视频| 在线免费观看不下载黄p国产 | 在线观看美女被高潮喷水网站 | 久久精品国产自在天天线| 亚洲av熟女| 亚洲成人久久爱视频| 亚洲成人久久爱视频| 18+在线观看网站| 全区人妻精品视频| 男人舔女人下体高潮全视频| 国产精品久久久人人做人人爽| 国内毛片毛片毛片毛片毛片| 亚洲五月婷婷丁香| 精品久久久久久久久久免费视频| 亚洲精品一卡2卡三卡4卡5卡| 身体一侧抽搐| 极品教师在线免费播放| 18禁黄网站禁片午夜丰满| 亚洲第一电影网av| 我的老师免费观看完整版| 最新中文字幕久久久久| 中文字幕熟女人妻在线| 好男人在线观看高清免费视频| 成人一区二区视频在线观看| 国产极品精品免费视频能看的| 一a级毛片在线观看| 女人十人毛片免费观看3o分钟| 国产成人av教育| 最新美女视频免费是黄的| 婷婷亚洲欧美| 波多野结衣高清无吗| 日本精品一区二区三区蜜桃| 国产成人影院久久av| 国产麻豆成人av免费视频| 欧美日韩国产亚洲二区| 伊人久久大香线蕉亚洲五| 日韩欧美精品v在线| 18禁美女被吸乳视频| 精品午夜福利视频在线观看一区| 欧美绝顶高潮抽搐喷水| 免费看a级黄色片| 天美传媒精品一区二区| 国产精品久久久久久久久免 | 法律面前人人平等表现在哪些方面| 精品福利观看| 国产主播在线观看一区二区| 久久精品91无色码中文字幕| 法律面前人人平等表现在哪些方面| 亚洲va日本ⅴa欧美va伊人久久| 国产乱人视频| 亚洲精品色激情综合| 99精品久久久久人妻精品| av黄色大香蕉| 日韩精品中文字幕看吧| 他把我摸到了高潮在线观看| 欧美在线一区亚洲| 欧美中文综合在线视频| 动漫黄色视频在线观看| 国产欧美日韩一区二区三| 亚洲aⅴ乱码一区二区在线播放| 国产成人av教育| 免费看光身美女| 又粗又爽又猛毛片免费看| 88av欧美| 色视频www国产| 中文字幕高清在线视频| 日韩欧美国产一区二区入口| 日日干狠狠操夜夜爽| 亚洲精品456在线播放app | tocl精华| 99久久综合精品五月天人人| 一区二区三区激情视频| 99精品欧美一区二区三区四区| 欧美激情在线99| 国产又黄又爽又无遮挡在线| 欧美三级亚洲精品| 九色成人免费人妻av| 十八禁人妻一区二区| 午夜久久久久精精品| 激情在线观看视频在线高清| 色在线成人网| 欧美成狂野欧美在线观看| 99久久成人亚洲精品观看| 国产精品99久久久久久久久| 草草在线视频免费看| 亚洲欧美一区二区三区黑人| 婷婷丁香在线五月| 少妇的逼水好多| 日本免费一区二区三区高清不卡| 色综合站精品国产| 午夜精品久久久久久毛片777| 久久精品亚洲精品国产色婷小说| 久久精品国产自在天天线| 日日干狠狠操夜夜爽| 国产在线精品亚洲第一网站| 国内精品美女久久久久久| 黄片小视频在线播放| 狂野欧美白嫩少妇大欣赏| 日日夜夜操网爽| 久久精品国产综合久久久| 日本黄大片高清| 身体一侧抽搐| 少妇熟女aⅴ在线视频| 精品一区二区三区视频在线观看免费| 精品免费久久久久久久清纯| 国产精品日韩av在线免费观看| 一本综合久久免费| 免费观看人在逋| 国产一区二区激情短视频| 深爱激情五月婷婷| 日韩亚洲欧美综合| 三级毛片av免费| 级片在线观看| 韩国av一区二区三区四区| 欧美精品啪啪一区二区三区| 国产极品精品免费视频能看的| 久久精品人妻少妇| 国产男靠女视频免费网站| 草草在线视频免费看| 久久99热这里只有精品18| 嫩草影院精品99| 亚洲狠狠婷婷综合久久图片| 国产精品爽爽va在线观看网站| 免费观看人在逋| 在线观看日韩欧美| 最新在线观看一区二区三区| 久久久久精品国产欧美久久久| 啦啦啦观看免费观看视频高清| 搡老熟女国产l中国老女人| 熟女电影av网| 91av网一区二区| 成人午夜高清在线视频| 成人av在线播放网站| 内地一区二区视频在线| 亚洲午夜理论影院| 亚洲国产欧美网| 亚洲专区国产一区二区| 国产精品久久久久久精品电影| 一夜夜www| 一进一出抽搐动态| 国产久久久一区二区三区| 欧美高清成人免费视频www| 日韩免费av在线播放| 18禁黄网站禁片免费观看直播| 久久99热这里只有精品18| 国产一区二区三区在线臀色熟女| 91在线精品国自产拍蜜月 | 啪啪无遮挡十八禁网站| 亚洲欧美日韩卡通动漫| 最后的刺客免费高清国语| 人妻久久中文字幕网| 欧美中文综合在线视频| 少妇的逼水好多| 欧美成狂野欧美在线观看| 一进一出抽搐动态| 亚洲第一欧美日韩一区二区三区| 精品人妻一区二区三区麻豆 | 手机成人av网站| av天堂中文字幕网| 成年女人看的毛片在线观看| 在线观看av片永久免费下载| 久久国产乱子伦精品免费另类| 91麻豆精品激情在线观看国产| 亚洲av免费在线观看| 欧美性猛交黑人性爽| 俄罗斯特黄特色一大片| 免费人成视频x8x8入口观看| 欧美xxxx黑人xx丫x性爽| 成人性生交大片免费视频hd| 最新在线观看一区二区三区| 亚洲avbb在线观看| 日韩欧美在线二视频| 国产在线精品亚洲第一网站| 亚洲国产日韩欧美精品在线观看 | 亚洲国产精品999在线| 亚洲18禁久久av| 精品免费久久久久久久清纯| 一进一出抽搐动态| 亚洲精品乱码久久久v下载方式 | 亚洲精品日韩av片在线观看 | 成人三级黄色视频| 日韩欧美免费精品| 蜜桃久久精品国产亚洲av| 免费看a级黄色片| 久久6这里有精品| 在线观看免费午夜福利视频| 熟女电影av网| 极品教师在线免费播放| 国产单亲对白刺激| 日本五十路高清| 欧美激情在线99| 色综合婷婷激情| 国产精品香港三级国产av潘金莲| 久久午夜亚洲精品久久| 九色国产91popny在线| 女人高潮潮喷娇喘18禁视频| 亚洲精品久久国产高清桃花| 亚洲专区国产一区二区| 哪里可以看免费的av片| 最后的刺客免费高清国语| 少妇熟女aⅴ在线视频| 99热只有精品国产| 校园春色视频在线观看| 日本黄大片高清| 国产精品自产拍在线观看55亚洲| 亚洲在线自拍视频| 日本与韩国留学比较| 女人十人毛片免费观看3o分钟| 亚洲成av人片在线播放无| www日本黄色视频网| 青草久久国产| 国产亚洲精品一区二区www| 成人国产一区最新在线观看| 69av精品久久久久久| 日韩欧美国产在线观看| 精品国产超薄肉色丝袜足j| 亚洲黑人精品在线| 可以在线观看毛片的网站| 国产精品久久久人人做人人爽| 久99久视频精品免费| 好看av亚洲va欧美ⅴa在| 亚洲国产精品sss在线观看| 午夜福利在线观看吧| 波多野结衣巨乳人妻| 免费高清视频大片| xxxwww97欧美| 日本黄色视频三级网站网址| 精品久久久久久久末码| 色综合婷婷激情| 2021天堂中文幕一二区在线观| 国产成人aa在线观看| 熟妇人妻久久中文字幕3abv| 香蕉av资源在线| 亚洲精品乱码久久久v下载方式 | 一级黄片播放器| 久久精品国产亚洲av涩爱 | 12—13女人毛片做爰片一| 九九热线精品视视频播放| 亚洲av成人精品一区久久| 女人高潮潮喷娇喘18禁视频| 禁无遮挡网站| 欧美三级亚洲精品| 最近最新中文字幕大全免费视频| 内地一区二区视频在线| 最新中文字幕久久久久| 国产中年淑女户外野战色| 三级国产精品欧美在线观看| 91在线精品国自产拍蜜月| 亚洲天堂国产精品一区在线| 成人国产麻豆网| 国产老妇伦熟女老妇高清| 亚洲国产日韩欧美精品在线观看| 91久久精品国产一区二区三区| 免费黄网站久久成人精品| 国产午夜福利久久久久久| 国语对白做爰xxxⅹ性视频网站| 国产精品人妻久久久久久| eeuss影院久久| 蜜桃亚洲精品一区二区三区| 日韩欧美 国产精品| 日本一二三区视频观看| 天天一区二区日本电影三级| 久久99热这里只频精品6学生| 干丝袜人妻中文字幕| 男人舔奶头视频| 51国产日韩欧美| 大话2 男鬼变身卡| 好男人视频免费观看在线| 三级经典国产精品| 成人无遮挡网站| 免费少妇av软件| 黑人高潮一二区| 国产高清三级在线| 少妇丰满av| 乱系列少妇在线播放| 波野结衣二区三区在线| 日韩强制内射视频| 久久精品夜色国产| 国产三级在线视频| 男人爽女人下面视频在线观看| 国产一区二区三区综合在线观看 | 一个人看视频在线观看www免费| 91午夜精品亚洲一区二区三区| 亚洲自偷自拍三级| 六月丁香七月| 亚洲精品456在线播放app| 欧美日韩一区二区视频在线观看视频在线 | 国产老妇伦熟女老妇高清| 99热这里只有是精品在线观看| 观看免费一级毛片| 免费观看精品视频网站| 精品酒店卫生间| 男的添女的下面高潮视频| 精品久久国产蜜桃| 国产在视频线精品| 日韩在线高清观看一区二区三区| 久久久久国产网址| 肉色欧美久久久久久久蜜桃 | 久久久久久久久久久丰满| 午夜免费男女啪啪视频观看| 肉色欧美久久久久久久蜜桃 | videossex国产| 成人美女网站在线观看视频| 人妻夜夜爽99麻豆av| xxx大片免费视频| 一级毛片我不卡| 尤物成人国产欧美一区二区三区| h日本视频在线播放| 久久精品人妻少妇| 人妻制服诱惑在线中文字幕| freevideosex欧美| 色播亚洲综合网| 亚洲国产成人一精品久久久| 国产高清不卡午夜福利| 丰满人妻一区二区三区视频av| 欧美zozozo另类| 国产精品久久久久久久久免| 国产伦一二天堂av在线观看| 欧美bdsm另类| 欧美日韩国产mv在线观看视频 | 久热久热在线精品观看| 欧美高清成人免费视频www| 国产亚洲5aaaaa淫片| 精品一区在线观看国产| 在线免费十八禁| 亚洲人成网站在线观看播放| freevideosex欧美| 国产国拍精品亚洲av在线观看| 一区二区三区免费毛片| 久久久久久久亚洲中文字幕| 日本免费在线观看一区| 亚州av有码| 国产黄色视频一区二区在线观看| 国产美女午夜福利| 久久久a久久爽久久v久久| 99视频精品全部免费 在线| 亚洲av成人av| 久久久久网色| 十八禁国产超污无遮挡网站| 亚洲成色77777| 亚洲av男天堂| 老司机影院成人| 午夜福利成人在线免费观看| 秋霞伦理黄片| 五月伊人婷婷丁香| 嫩草影院精品99| 少妇熟女欧美另类| 日日摸夜夜添夜夜添av毛片| 国产一区二区三区av在线| 国产精品女同一区二区软件| 亚洲精品国产av蜜桃| 精品久久久久久电影网| 欧美 日韩 精品 国产| 在线观看一区二区三区| 又爽又黄无遮挡网站| 亚洲电影在线观看av| 成年人午夜在线观看视频 | 少妇被粗大猛烈的视频| 中国国产av一级| 久久久久久久亚洲中文字幕| 色综合亚洲欧美另类图片| 爱豆传媒免费全集在线观看| 久久久久久久久中文| 高清在线视频一区二区三区| 简卡轻食公司| 精品久久久久久久人妻蜜臀av| 国内少妇人妻偷人精品xxx网站| 99久久中文字幕三级久久日本| 一个人观看的视频www高清免费观看| 国产午夜精品论理片| 菩萨蛮人人尽说江南好唐韦庄| 高清欧美精品videossex| 久久久精品免费免费高清| 一区二区三区高清视频在线| 国产一区亚洲一区在线观看| 国产精品久久久久久av不卡| 国产伦在线观看视频一区| 搡女人真爽免费视频火全软件| 国产精品久久久久久av不卡| 天堂网av新在线| 国产精品福利在线免费观看| 久久6这里有精品| 国产亚洲av片在线观看秒播厂 | 成人亚洲精品一区在线观看 | 黄色欧美视频在线观看| 麻豆久久精品国产亚洲av| 高清av免费在线| 一区二区三区免费毛片| 伊人久久国产一区二区| 欧美日韩在线观看h| 天天躁日日操中文字幕| 一级毛片aaaaaa免费看小| 91aial.com中文字幕在线观看| 亚洲精品影视一区二区三区av| 亚洲av在线观看美女高潮| 国产免费一级a男人的天堂| 一级片'在线观看视频| 成人性生交大片免费视频hd| 肉色欧美久久久久久久蜜桃 | 亚洲精品中文字幕在线视频 | 狂野欧美激情性xxxx在线观看| 国产精品一二三区在线看| 欧美最新免费一区二区三区| 国产国拍精品亚洲av在线观看| 久久久久性生活片| 日韩,欧美,国产一区二区三区| 在线免费观看的www视频| 狂野欧美白嫩少妇大欣赏| 天天一区二区日本电影三级| 美女黄网站色视频| 51国产日韩欧美| 伊人久久精品亚洲午夜| 亚洲精品成人久久久久久| 精品久久久久久久末码| av专区在线播放| 亚洲精品成人久久久久久| 午夜精品在线福利| 国产又色又爽无遮挡免| 亚洲精品视频女| 日本免费a在线| 色尼玛亚洲综合影院| 少妇丰满av| 日韩亚洲欧美综合| 亚洲性久久影院| 色综合亚洲欧美另类图片| 人妻少妇偷人精品九色| 国产精品不卡视频一区二区| 久久国产乱子免费精品| 国产亚洲一区二区精品| 网址你懂的国产日韩在线| 亚洲一级一片aⅴ在线观看| 国产精品国产三级国产专区5o| 国产白丝娇喘喷水9色精品| 色哟哟·www| 国产69精品久久久久777片| 成人高潮视频无遮挡免费网站| 中国国产av一级| 又黄又爽又刺激的免费视频.| 亚洲乱码一区二区免费版| 免费观看的影片在线观看| 亚洲在线自拍视频| 亚洲精品视频女| 久99久视频精品免费| 女人久久www免费人成看片| 一级a做视频免费观看| 精品99又大又爽又粗少妇毛片| 亚洲精品视频女| 一二三四中文在线观看免费高清| 成年女人在线观看亚洲视频 | 欧美成人午夜免费资源| www.色视频.com| 久久精品国产亚洲av天美| 日韩伦理黄色片| 别揉我奶头 嗯啊视频| 亚洲av男天堂| 老女人水多毛片| 亚洲自偷自拍三级| 亚洲丝袜综合中文字幕| 亚洲最大成人av| 久久久久免费精品人妻一区二区| 亚洲av一区综合| 免费大片黄手机在线观看| 久久6这里有精品| 伦理电影大哥的女人| 美女内射精品一级片tv| 国产白丝娇喘喷水9色精品| 欧美精品国产亚洲| av在线老鸭窝| 美女黄网站色视频| 天堂网av新在线| 久久久久久国产a免费观看| 久久久欧美国产精品| 亚洲最大成人av| av在线亚洲专区| 亚洲aⅴ乱码一区二区在线播放| 欧美一级a爱片免费观看看| 亚洲国产精品sss在线观看| 在线 av 中文字幕| 国产伦在线观看视频一区| 热99在线观看视频| av黄色大香蕉| 91aial.com中文字幕在线观看| 一本久久精品| 日韩国内少妇激情av| 久久人人爽人人爽人人片va| 免费在线观看成人毛片| 亚洲国产高清在线一区二区三| 舔av片在线| 国产色婷婷99| 久久精品夜夜夜夜夜久久蜜豆| 国产亚洲精品久久久com| 小蜜桃在线观看免费完整版高清| 亚洲最大成人av| 亚洲在线自拍视频| 亚洲欧洲国产日韩| 青春草国产在线视频| 深夜a级毛片| 51国产日韩欧美| 亚洲欧美一区二区三区黑人 | 国产 一区 欧美 日韩| 欧美xxⅹ黑人| 亚洲精品久久午夜乱码| 免费无遮挡裸体视频| 一级毛片电影观看| 亚洲人成网站高清观看| 日本午夜av视频| 99热全是精品| 日本色播在线视频| 日韩 亚洲 欧美在线| 国产精品一区二区在线观看99 | 精品不卡国产一区二区三区| 3wmmmm亚洲av在线观看| 亚洲高清免费不卡视频| 国产午夜精品一二区理论片| 国产精品精品国产色婷婷| 国产 一区 欧美 日韩| 黄片wwwwww| 久久精品国产亚洲网站| 美女脱内裤让男人舔精品视频| 美女xxoo啪啪120秒动态图| 黄色一级大片看看| 免费看美女性在线毛片视频| 久久久久久久亚洲中文字幕| 少妇被粗大猛烈的视频| 国产亚洲av嫩草精品影院| 三级男女做爰猛烈吃奶摸视频| 69av精品久久久久久| 日韩成人伦理影院| 中文精品一卡2卡3卡4更新| 久久久久久久久久久丰满|