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

    基于內(nèi)表面特征的碗狀碎塊匹配方法

    2022-10-12 05:58:50習(xí)俊通朱興龍
    機械設(shè)計與制造 2022年10期
    關(guān)鍵詞:碎塊曲率高斯

    孫 進,丁 煜,習(xí)俊通,朱興龍

    (1.揚州大學(xué)機械工程學(xué)院,江蘇 揚州 225100;2.上海交通大學(xué)機械系統(tǒng)與振動國家重點實驗室,上海 200240)

    1 引言

    人類燦爛的文明和悠久的歷史使得河流、湖泊、海洋下隱藏著眾多的文化遺跡,其中以碗、瓶、罐、壺、盤等為器型的瓷器在水下文物研究中占有重要地位。備受矚目的宋代沉船“南海Ⅰ號”近期結(jié)束清理工作,出水的十八多萬件文物中以碗狀器型存在的瓷器為多數(shù)[1]。這些水下文物在水底經(jīng)過長時間的侵蝕和顛簸,大多數(shù)變得破碎殘缺。隨著信息技術(shù)的不斷發(fā)展,利用計算機的虛擬復(fù)原技術(shù)來輔助文物碎塊的實體拼合逐漸成為可能。為便于文物數(shù)字化表征,基于厚度值大致將破損物體分為兩類:厚度小于1mm的如紙幣、油畫和地圖等,稱為碎片;厚度大于1mm的如碗、壁畫和兵馬俑等,稱為碎塊??紤]到水下文物中大部分是以碗狀器型存在的物件,故以碗狀碎塊為研究對象。

    傳統(tǒng)的碎塊匹配方法主要圍繞斷裂曲面展開研究。研究方向上大體可分為三類:第一類是基于碎塊斷裂曲面特征的重組匹配。文獻[2]提出利用碎塊表面的積分不變量刻畫斷裂曲面的彎曲能量,并采用向前搜索算法以及表面一致性約束完成碎塊斷裂曲面的匹配。文獻[3]直接根據(jù)斷裂曲面上所有頂點的三維坐標(biāo)和曲率信息,采用基于隨機采樣算法和分級二叉樹進行碎塊的配對。目前,文獻[4]提出了確定碎塊匹配對的四個準(zhǔn)則,即兩個匹配碎塊的間隙體積、碎塊在匹配位置上的重疊性、兩個接觸面的平均曲率和接觸弧的長度,該準(zhǔn)則也得到了廣泛應(yīng)用。在此基礎(chǔ)上,文獻[5]先是根據(jù)斷裂曲面的顯著性區(qū)域?qū)崿F(xiàn)碎塊的粗匹配,然后通過加入高斯概率模型、角度約束和動態(tài)迭代系數(shù)的方式來改進ICP算法完成精匹配,提高了碎塊匹配的正確率和對齊精度。但是斷裂曲面的面特征信息過多,直接運算會造成算法的效率低下。一些專家學(xué)者受到碎片拼接方法[6-7]的啟發(fā),提出基于碎塊輪廓曲線特征的重組匹配。文獻[8]引進一種新的描述子表示碎塊輪廓曲線上的特征點和曲線,然后使用快速傅里葉變換完成碎塊的兩兩匹配,其中特征點的描述子是由曲率計算得到的,曲線的描述子用傅里葉級數(shù)表示。文獻[9]針對碎塊表面的特征輪廓線,提出了一種使用曲線矩的不變量和創(chuàng)建在碎塊邊界區(qū)域的錨點來實現(xiàn)碎塊匹配和對齊的方法,并通過距離誤差度量來估算特征輪廓線對齊的準(zhǔn)確性。文獻[10]識別出砂粒碎塊特征曲線網(wǎng)絡(luò)中所有簡單的無弦周期,然后采用Hausdorff距離法和四點同余集算法對潛在的無弦周期進行匹配,從而完成對斷裂砂粒碎塊的重組匹配。又由于基于輪廓曲線的拼合方案被限定在碎塊形狀互補的理想條件下,任何幾何外形上微小的缺失都會對算法產(chǎn)生很大的影響,基于自適應(yīng)交互式的匹配方法應(yīng)運而生。文獻[11]提出一種基于自適應(yīng)鄰域的多特征融合文物碎塊自動拼接方法,該方法基于自適應(yīng)鄰域機制融合碎塊斷裂面及原曲面上的多種特征,將斷裂部位輪廓線的匹配轉(zhuǎn)化為輪廓帶的匹配,通過自適應(yīng)鄰域的匹配反饋碎塊的匹配關(guān)系。

    上述的碎塊匹配算法比較單一,僅圍繞斷裂曲面的幾何特征進行。倘若碎塊的斷裂曲面有部分幾何特征遺失,那么容易在重組過程中產(chǎn)生錯誤的匹配;而且在曲面分割過程中,斷裂曲面的“過分割”現(xiàn)象嚴(yán)重,需要研究者在提取斷裂曲面特征前進行曲面合并,增加了重組算法的運行時間。相較于斷裂曲面可能存在的問題,碗狀碎塊的內(nèi)表面不僅具有和斷裂曲面一樣的幾何信息,還擁有比斷裂曲面更豐富的顏色信息和紋理信息,同時內(nèi)表面的識別方法不需要對斷裂曲面分割后產(chǎn)生的眾多曲面進行合并,只需要正確區(qū)分出內(nèi)表面即可,可以巧妙地避開過分割帶來的問題。另外,當(dāng)現(xiàn)有文獻在討論碎塊匹配算法時,往往默認待匹配的一對碎塊是可以拼接的,并沒有展開碎塊篩選方法的研究,而從碎塊集中篩選出最佳匹配的一對碎塊是計算機實現(xiàn)自動重組的必經(jīng)之路。受碎片拼接算法[12]的啟發(fā),考慮到碎塊的表面都是由復(fù)雜的曲面構(gòu)成的,不同的曲面都具有自身的特殊性,這里將這些曲面的幾何特征和顏色特征通過高斯映射在擴展高斯球[13]中表現(xiàn)出來。鑒于此,提出一種基于內(nèi)表面特征的碗狀碎塊匹配方法,用于協(xié)助計算機自動篩選出兩個最相似的待匹配碎塊。

    2 方法

    2.1 碗狀碎塊的內(nèi)表面提取

    在一般的碎塊模型中,主要存在兩種類型的曲面:斷裂曲面和非斷裂曲面,通??梢愿鶕?jù)曲面的粗糙程度[5]來識別斷裂曲面。相比于非斷裂曲面,斷裂曲面上空間點的法向量變化程度更大,表面更凹凸不平。而對于碗狀碎塊的點云數(shù)據(jù)模型,采用區(qū)域生長算法[14]經(jīng)過曲面分割后,如圖1所示??梢詫⒎菙嗔亚婕毞譃榈酌妗⒃急砻?、外表面和內(nèi)表面四種類型,其中原始表面指的是區(qū)域生長算法在外表面和內(nèi)表面交界處過度分割的曲面,這部分曲面表面積較小,僅存在外表面和內(nèi)表面相鄰的部分碗狀碎塊中。

    圖1 碗狀碎塊的曲面分類Fig.1 Surface Classification of Bowl-Shaped Broken Pieces

    這里采用點云占比定性分析碗狀碎塊各分割曲面的性質(zhì),這主要是源于工程實踐中的發(fā)現(xiàn):相比于其他類型的碎塊,碗狀碎塊的厚度較薄,導(dǎo)致其外表面和內(nèi)表面的表面積要大于其他曲面,也造成外表面和內(nèi)表面的點云占比明顯大于其他曲面。點云占比可以描述碎塊中各分割曲面的大小,一般采用曲面的空間點個數(shù)在碎塊點云中所占的百分比來表示某一曲面的點云占比數(shù)值。通常根據(jù)各個曲面的點云占比,可以從碗狀碎塊的眾多分割曲面中識別出外表面和內(nèi)表面。

    考慮到碗狀物體在平時生活中的實際應(yīng)用,外表面常用于抓取接觸,內(nèi)表面要便于清洗使用,而且一般情況下外表面還會設(shè)計一些凸出的紋飾特征,因此內(nèi)表面一般比外表面更加光滑平整。數(shù)學(xué)表征為:相比于外表面,內(nèi)表面上某一空間點的法向量與其周圍點的法向量變化程度較小。為了便于定量表征,這里引入平滑度的概念。假設(shè)曲面S由n個點構(gòu)成,即S={p1,p2,…,pn},其中第i個點pi(1 ≤i≤n)的鄰域點為pj(1 ≤j≤k),那么曲面S的平滑度ε可以表示為:

    式中:φi,j—pi的法向量與其鄰域點pj的法向量的夾角。

    ε的大小可以定量地反映曲面S的平滑度,曲面中相鄰點之間的法向量變化愈小,則說明曲面的平滑度愈大,因此碗狀碎塊內(nèi)表面的平滑度數(shù)值要大于外表面的平滑度數(shù)值,這樣基于點云占比和平滑度能從碗狀碎塊的眾多分割曲面中識別出內(nèi)表面。接著采用KD樹[15]建立內(nèi)表面上各空間點的位置關(guān)系,然后進行K鄰域搜索并求出各點的法向量和曲率,最后依據(jù)這些信息進行內(nèi)表面邊緣點的檢測和提取。

    2.2 碗狀碎塊內(nèi)表面的特征提取

    2.2.1 內(nèi)表面邊緣點的特征圓曲率計算

    內(nèi)表面邊緣點所構(gòu)成的輪廓線可以被看成是由一系列凹凸的平滑曲線組成,這些曲線凹凸部分的形狀類似拋物線。對于一段類拋物線,采用弧長、短邊與特征圓半徑之比、長邊與特征圓半徑之比和弦長與特征圓半徑之比來進行描述。P1、P2和P3均為碎塊內(nèi)表面的邊緣點,P1P3的弧長用s表示,O是由P1、P2和P3確定的特征圓的圓心,該特征圓的半徑用r表示,如圖2所示。類拋物線的幾何描述符可以表示為(s,a,b,c),其中a表示長邊與特征圓半徑之比即a=| |P2P1/r,b表示短邊與特征圓半徑之比即b=|P2P3|/r,c表示弦長與特征圓半徑之比即c=|P1P3|/r。

    圖2 類拋物線模型Fig.2 Likeness Parabolic Model

    在類拋物線模型中,根據(jù)圓周角定理、扇形弧長公式和余弦定理可以計算出特征圓半徑r,又因為特征圓曲率在數(shù)值上等于r的倒數(shù),則內(nèi)表面上第n個邊緣點的特征圓曲率gn可以表示為:

    2.2.2 改進后的擴展高斯球

    曲線和曲面的幾何特征可通過高斯映射[16]來表示。高斯映射是將曲線或曲面上所有點的法矢先進行單位化,并將法矢的起點平移至同一個端點上,則曲線的各法矢端點落在半徑為1 的單位圓上,曲面的各法矢端點落在單位球面上。法矢端點在圓上的投影點所構(gòu)成的圖像即為高斯圖,在球上的投影點所構(gòu)成的圖像即為高斯球。傳統(tǒng)擴展高斯球[17]的建立依賴于法向量和曲率,只能表征曲面的幾何特征,這里為同時表征碗狀碎塊內(nèi)表面的顏色特征和幾何特征,改進擴展高斯球的建立方法。

    首先單位化內(nèi)表面邊緣點的顏色矢量。為了讓內(nèi)表面邊緣點的顏色矢量能均勻分布在各個卦限內(nèi),將RGB顏色空間模型的重心平移至坐標(biāo)系的原點,如圖3所示。那么單位化后的邊緣點顏色矢量γn可表示為:

    圖3 RGB顏色空間模型Fig.3 RGB Color Space Model

    式中:x=r-255/2,y=g-255/2,z=b-255/2,且(r,g,b)—內(nèi)表面邊緣點的RGB 顏色值;ε,μ和φ—x軸,y軸和z軸上的單位向量。

    然后附加特征圓曲率信息。將單位化后的顏色矢量起點平移至原點,則其端點落在高斯球的球面上,再將特征圓曲率信息附加至高斯球的端點便形成擴展高斯球。改進后的擴展高斯球,邊緣點顏色矢量的起點落在擴展高斯球的球心上,端點落在球面上,每個顏色矢量都帶有該點的特征圓曲率。

    2.3 基于相似度的碗狀碎塊篩選

    建立擴展高斯球的目的是為了尋找兩塊相似的待匹配碎塊之間的對應(yīng)點,由于碎塊內(nèi)表面輪廓曲線的擴展高斯球都是在球坐標(biāo)系下,所以使尋找對應(yīng)點成為可能。相似的兩塊待匹配碎塊之間由于RGB顏色特征和特征圓曲率的計算方法相同,其擴展高斯球上必有一部分形狀相似。在球坐標(biāo)系下,計算碎塊內(nèi)表面邊緣點的顏色矢量與平面XOY的夾角,并按照夾角的大小沿著球半徑將擴展高斯球平均分割成q個空間區(qū)域,令gM()為碎塊內(nèi)表面輪廓曲線M在第x個空間區(qū)域qx內(nèi)第i個邊緣點的特征圓曲率,其中i=1,2,…,m,同樣令gN()為碎塊內(nèi)表面輪廓曲線N在第x個空間區(qū)域qx內(nèi)第j個邊緣點的特征圓曲率,其中j=1,2,…,n,若將M和N在空間區(qū)域qx內(nèi)具有相同特征圓曲率的邊緣點個數(shù)記為f(qx),則f(qx)可用如下公式表示:

    式中:gM(qx),gN(qx)—M和N在第x個空間區(qū)域qx內(nèi)的邊緣點個數(shù);q—沿著球半徑將擴展高斯球分割得到的空間區(qū)域數(shù)量。相似度的數(shù)值越大,則代表兩塊待匹配的碗狀碎塊越相似。

    2.4 匹配策略

    結(jié)合基于內(nèi)表面輪廓特征的粗匹配方法和基于ICP算法的精匹配方法[5],提出一種基于內(nèi)表面特征的碗狀碎塊重組策略,如圖4所示。

    圖4 碗狀碎塊的整體匹配策略Fig.4 Whole Reassembly Strategy of Bowl-Shaped Broken Pieces

    具體方法步驟如下:

    (1)將待匹配的碎塊分別標(biāo)記為block1,block2,…,blockn,狀態(tài)設(shè)置為未訪問,然后提取這些未訪問碎塊的內(nèi)表面及其輪廓曲線,以碎塊blocki(1 ≤i≤n)為基準(zhǔn)模型,計算其他碎塊內(nèi)表面輪廓曲線與碎塊blocki內(nèi)表面輪廓曲線的相似度;

    (3)若拼接誤差在允許范圍內(nèi),則選擇ICP算法完成精匹配,并將碎塊blocki的狀態(tài)設(shè)置為已訪問,然后將更新輪廓后的碎塊設(shè)置為基準(zhǔn),回到(1)繼續(xù)計算相似度;若拼接誤差不在允許范圍內(nèi),則回到(2),選擇s值第二大的碎塊與基準(zhǔn)碎塊進行拼接,直至拼接誤差在允許范圍內(nèi);若拼接誤差始終不在允許范圍內(nèi),則將碎塊blocki剔除出待匹配的碎塊集合,重新選擇基準(zhǔn)碎塊進行拼接;

    (4)每次精匹配結(jié)束后更新碎塊的輪廓,直至完成碎塊的整體拼合,待所有碎塊均標(biāo)記為已訪問狀態(tài)時,拼接結(jié)束。

    3 實驗

    3.1 實驗結(jié)果

    實驗的軟件架構(gòu)為Visual Studio 2010,VTK 5.8.0(Visualization Toolkit,視覺化工具函式庫)和PCL1.7.2(Point Cloud Library,點云庫),實驗在Intel Core i3-4160 CPU@3.60Hz 處理器,4.0GB內(nèi)存的PC上運行。碗狀碎塊的三維點云數(shù)據(jù)由Einscan-S三維激光掃描儀采集,首先建立碗狀碎塊的擴展高斯球模型,如圖5所示。高斯球上點的矢量方向代表了邊緣點的顏色特征,擴展高斯球上點的向量模代表了邊緣點的特征圓曲率,此時擴展高斯球能完整地表征出碗狀碎塊內(nèi)表面的幾何特征和顏色特征。

    圖5 擴展高斯球的建立過程Fig.5 Establishment Process of the Extended Gauss Sphere

    接著本實驗將碎塊01作為基準(zhǔn)碎塊,剩余的碎塊作為待匹配碎塊,進行相似度的交叉計算,如圖6所示。相似度計算結(jié)果表明與碎塊01最相似的是碎塊10。根據(jù)2.4節(jié)中所提出的匹配策略,對篩選出碎塊01和碎塊10進行拼接匹配。在一對碎塊的重組完成后進行輪廓更新,篩選出下一個最佳待匹配碎塊,然后進行下一步的拼接匹配,最后直至完成碎塊集中所有碗狀碎塊的整體拼合。

    3.2 實驗討論

    傳統(tǒng)的擴展高斯球模型建立方法依賴于點云的法向量和曲率值[16-18],對于幾何特征明顯的碗狀碎塊可以很好地篩選出最佳待匹配碎塊。但是對于區(qū)域相對平坦的碗狀碎塊,單純依靠幾何特征的方法容易產(chǎn)生錯誤篩選。如圖7所示,傳統(tǒng)的擴展高斯球建立方法在最佳待匹配模型的篩選過程中存在明顯的錯誤。碗狀碎塊集B中碎塊01的最佳待匹配碎塊應(yīng)為碎塊02,而相似度計算結(jié)果卻表明碎塊01的最佳待匹配碎塊是碎塊06。在同樣的情況下,如圖8所示?;赗GB顏色特征和特征圓曲率的擴展高斯球模型建立方法,可以正確地完成碎塊01的最佳待匹配碎塊的篩選。

    圖8 我們方法在碎塊集B的正確篩選Fig.8 Correct Screening Result in Broken Pieces Set B with Our Method

    為了進一步驗證以RGB顏色特征與特征圓曲率為基礎(chǔ)建立的擴展高斯球模型更適合于最佳待匹配碎塊的篩選,分別使用我們方法(即RGB顏色特征與特征圓曲率相結(jié)合的方法),初始方法(即HSV顏色特征與曲率相結(jié)合的方法)和傳統(tǒng)方法(即法向量與曲率相結(jié)合的方法[16-18])對碎塊集A(圖6所示的碎塊)和碎塊集B(圖7和圖8所示的碎塊)進行重組匹配。如表1所示,相似度計算結(jié)果包括篩選準(zhǔn)確率和算法運行時間。計算結(jié)果表明結(jié)合RGB顏色特征與特征圓曲率建立擴展高斯球模型來進行最佳待匹配碎塊的篩選,碎塊篩選的平均準(zhǔn)確率可達到93.87%。與初始方法相比,我們方法的計算準(zhǔn)確率比其高5.78%,平均運算速度提高了20.64%;與傳統(tǒng)方法相比,我們方法的計算準(zhǔn)確率比其高6.60%,平均運算速度提高了16.08%。因此基于RGB 顏色特征與特征圓曲率建立擴展高斯球模型的方法更適合于碗狀碎塊最佳待匹配碎塊的篩選,在不同碗狀碎塊數(shù)據(jù)上的實驗表明,該方法相比于其他擴展高斯球模型的建立方法更加優(yōu)越。

    表1 相似度交叉計算結(jié)果的比較Tab.1 Comparison of Similarity Calculation Result

    4 結(jié)束語

    傳統(tǒng)的碎塊匹配方法通過研究斷裂曲面的幾何性質(zhì)來達到目的,容易在曲面分割時產(chǎn)生過分割現(xiàn)象,從而增加了拼接過程的時間,為了避免這一情況,這里結(jié)合點云占比和平滑度提取碗狀碎塊的內(nèi)表面,并將碗狀碎塊的內(nèi)表面作為研究對象。此外,為了幫助計算機自動識別兩個最相似的待匹配碎塊,建立擴展高斯球模型將碎塊內(nèi)表面邊緣點的幾何特征和顏色特征融合起來,并通過構(gòu)造類拋物線模型描述內(nèi)表面邊緣點的幾何特征,然后定義相似度度量函數(shù)來描述擴展高斯球模型的全局特征,最后選擇相似度最大的一對碎塊進行碎塊匹配,從而完成對碗狀碎塊的整體拼合。

    實驗表明結(jié)合RGB顏色特征與特征圓曲率建立擴展高斯球模型來進行相似度計算,相似度的平均計算準(zhǔn)確率可達到93.87%,相比于傳統(tǒng)的擴展高斯球模型建立方法,結(jié)合RGB顏色特征與特征圓曲率的方法在相似度計算準(zhǔn)確率和平均運算速度上分別提高了6.60%和16.08%。在接下來的工作中,研究者不會僅僅局限于提高碗狀碎塊的匹配準(zhǔn)確率和拼接效率,將進一步展開對不同形狀碎塊的相似度計算方法和匹配策略的研究,從而為建立水下文物虛擬拼接復(fù)原的工程化體系框架作鋪墊。

    猜你喜歡
    碎塊曲率高斯
    小高斯的大發(fā)現(xiàn)
    大曲率沉管安裝關(guān)鍵技術(shù)研究
    一類雙曲平均曲率流的對稱與整體解
    希臘遺址
    溫度效應(yīng)對隧道圍巖影響的試驗研究
    半正迷向曲率的四維Shrinking Gradient Ricci Solitons
    天才數(shù)學(xué)家——高斯
    淺析槍擊鋼化玻璃矩形碎塊特征
    有限域上高斯正規(guī)基的一個注記
    Esn+1中具有至多兩個不同主曲率的2-調(diào)和超曲面
    噜噜噜噜噜久久久久久91| 男女啪啪激烈高潮av片| 国产真实伦视频高清在线观看| 国产一区有黄有色的免费视频| 久久久色成人| 成人毛片a级毛片在线播放| 全区人妻精品视频| 97在线人人人人妻| 国产av国产精品国产| 精品人妻偷拍中文字幕| 免费大片黄手机在线观看| 成人美女网站在线观看视频| 亚洲国产色片| 久久久久久久精品精品| 免费av毛片视频| freevideosex欧美| 国产精品.久久久| 欧美一级a爱片免费观看看| 男人狂女人下面高潮的视频| 午夜福利视频1000在线观看| 中国三级夫妇交换| 日韩欧美精品v在线| 国内精品美女久久久久久| 十八禁网站网址无遮挡 | 午夜老司机福利剧场| 舔av片在线| 精品久久久久久久久亚洲| 成人毛片60女人毛片免费| 久久亚洲国产成人精品v| 亚洲精品一区蜜桃| 51国产日韩欧美| 国产免费视频播放在线视频| 日韩三级伦理在线观看| 久久国产乱子免费精品| 欧美精品人与动牲交sv欧美| 99久久九九国产精品国产免费| 成年女人在线观看亚洲视频 | 全区人妻精品视频| 国产精品麻豆人妻色哟哟久久| 久久韩国三级中文字幕| 日韩av免费高清视频| 日本wwww免费看| 一级片'在线观看视频| 中文天堂在线官网| 国产精品国产av在线观看| 最近最新中文字幕大全电影3| 亚洲精品日本国产第一区| 国产精品人妻久久久影院| 国产熟女欧美一区二区| 日韩视频在线欧美| 亚洲人与动物交配视频| 97在线视频观看| 亚洲精品,欧美精品| 一个人看的www免费观看视频| 熟女av电影| 成年女人在线观看亚洲视频 | 国产一区二区在线观看日韩| 亚洲精品日韩在线中文字幕| 成年女人看的毛片在线观看| 大香蕉久久网| 亚洲性久久影院| 欧美变态另类bdsm刘玥| 欧美一级a爱片免费观看看| 97在线人人人人妻| 中文字幕免费在线视频6| 啦啦啦啦在线视频资源| 国产黄色视频一区二区在线观看| 春色校园在线视频观看| 精品久久久久久久末码| 最近中文字幕高清免费大全6| 国产精品三级大全| 中文字幕亚洲精品专区| 超碰av人人做人人爽久久| 亚洲国产精品专区欧美| 中文在线观看免费www的网站| 亚洲精品中文字幕在线视频 | 亚洲人成网站高清观看| 一个人观看的视频www高清免费观看| 美女内射精品一级片tv| 国产成人freesex在线| 大又大粗又爽又黄少妇毛片口| 极品教师在线视频| 校园人妻丝袜中文字幕| 亚洲国产欧美人成| 国产精品不卡视频一区二区| 在线观看免费高清a一片| 欧美日本视频| 好男人在线观看高清免费视频| 国产 一区精品| 人人妻人人看人人澡| 夜夜爽夜夜爽视频| 乱系列少妇在线播放| 亚洲高清免费不卡视频| 国产精品福利在线免费观看| 不卡视频在线观看欧美| 国产欧美日韩精品一区二区| 国产 精品1| 青青草视频在线视频观看| 亚洲精品,欧美精品| av天堂中文字幕网| 99久久精品一区二区三区| 人妻系列 视频| 日韩成人av中文字幕在线观看| 777米奇影视久久| 最近最新中文字幕免费大全7| 免费电影在线观看免费观看| 精品酒店卫生间| 国产爽快片一区二区三区| 日本av手机在线免费观看| 国产精品福利在线免费观看| 亚洲精品一二三| 人妻系列 视频| 久久影院123| 精品久久久久久久久亚洲| .国产精品久久| 女的被弄到高潮叫床怎么办| 性插视频无遮挡在线免费观看| 国产一区二区在线观看日韩| 人人妻人人澡人人爽人人夜夜| 国产高清三级在线| 好男人视频免费观看在线| 最近最新中文字幕免费大全7| 亚洲怡红院男人天堂| 七月丁香在线播放| 亚洲精品日韩av片在线观看| 亚洲精品,欧美精品| 男女边摸边吃奶| 在线看a的网站| 亚洲电影在线观看av| 欧美最新免费一区二区三区| 大香蕉久久网| 亚洲av中文字字幕乱码综合| 又爽又黄a免费视频| 亚洲欧美日韩另类电影网站 | 日韩欧美精品免费久久| 高清av免费在线| 国产精品福利在线免费观看| 久久鲁丝午夜福利片| 久久久久网色| 国产成人精品婷婷| 亚洲精品亚洲一区二区| 人人妻人人看人人澡| 国产精品久久久久久精品古装| 在线观看国产h片| 97人妻精品一区二区三区麻豆| 大香蕉久久网| 一级爰片在线观看| 麻豆乱淫一区二区| 免费大片18禁| 精品久久久久久久久亚洲| 国产精品一区二区三区四区免费观看| 国产综合懂色| 乱系列少妇在线播放| 天美传媒精品一区二区| 国产黄色免费在线视频| 色综合色国产| 大又大粗又爽又黄少妇毛片口| 日韩 亚洲 欧美在线| 亚洲在线观看片| 色哟哟·www| 只有这里有精品99| 天天一区二区日本电影三级| 精品99又大又爽又粗少妇毛片| 嫩草影院入口| 亚洲国产欧美在线一区| 久久99精品国语久久久| 成年版毛片免费区| 午夜免费观看性视频| 亚洲不卡免费看| 99热6这里只有精品| 一区二区三区乱码不卡18| 狂野欧美激情性xxxx在线观看| 亚洲av成人精品一二三区| 久久久久久九九精品二区国产| 黄色日韩在线| 亚洲性久久影院| 欧美一级a爱片免费观看看| 高清在线视频一区二区三区| 少妇熟女欧美另类| 99久国产av精品国产电影| 国产免费视频播放在线视频| 搡老乐熟女国产| 国产一区亚洲一区在线观看| 黄片无遮挡物在线观看| 交换朋友夫妻互换小说| 国产精品一区二区三区四区免费观看| 国产精品人妻久久久久久| 午夜激情福利司机影院| 99久久精品国产国产毛片| 久久久久久久精品精品| 欧美少妇被猛烈插入视频| 国产精品一及| 精华霜和精华液先用哪个| 91午夜精品亚洲一区二区三区| 精品国产露脸久久av麻豆| 亚洲最大成人中文| 赤兔流量卡办理| 伊人久久国产一区二区| 99热这里只有是精品在线观看| 水蜜桃什么品种好| 亚洲国产色片| 久久久色成人| 高清欧美精品videossex| 欧美zozozo另类| 99九九线精品视频在线观看视频| 国产国拍精品亚洲av在线观看| 99久久中文字幕三级久久日本| 久久精品国产自在天天线| 国产精品偷伦视频观看了| 婷婷色麻豆天堂久久| 日韩人妻高清精品专区| 永久免费av网站大全| 欧美精品人与动牲交sv欧美| 婷婷色综合大香蕉| 免费观看无遮挡的男女| 精品少妇久久久久久888优播| 97热精品久久久久久| 一级黄片播放器| 欧美高清性xxxxhd video| 青春草视频在线免费观看| 日韩,欧美,国产一区二区三区| 欧美三级亚洲精品| 国产午夜精品久久久久久一区二区三区| 99久久精品热视频| 久久精品综合一区二区三区| 69av精品久久久久久| 国产人妻一区二区三区在| 狠狠精品人妻久久久久久综合| tube8黄色片| 亚洲国产精品999| 欧美区成人在线视频| 一区二区三区精品91| 亚洲国产精品专区欧美| 中国国产av一级| 久久午夜福利片| 秋霞在线观看毛片| 国产 一区 欧美 日韩| 伦精品一区二区三区| 寂寞人妻少妇视频99o| 国产 一区精品| 高清av免费在线| 国产精品久久久久久av不卡| 亚洲精品第二区| 国产 一区 欧美 日韩| 一个人看视频在线观看www免费| 国内精品美女久久久久久| 少妇被粗大猛烈的视频| 精品久久国产蜜桃| av国产久精品久网站免费入址| 亚洲精品成人av观看孕妇| 老女人水多毛片| 日日摸夜夜添夜夜添av毛片| 插逼视频在线观看| 日本与韩国留学比较| 国产成人精品婷婷| 男男h啪啪无遮挡| 51国产日韩欧美| 一级a做视频免费观看| 黄色怎么调成土黄色| 在线免费十八禁| 亚洲精品一区蜜桃| av在线亚洲专区| 国产成人精品久久久久久| av免费观看日本| 精品一区在线观看国产| 日韩欧美精品免费久久| 久热久热在线精品观看| 联通29元200g的流量卡| .国产精品久久| 深夜a级毛片| 成人特级av手机在线观看| 美女cb高潮喷水在线观看| 国产精品伦人一区二区| 久热这里只有精品99| 亚洲成人av在线免费| 五月伊人婷婷丁香| 亚洲最大成人av| 国产老妇伦熟女老妇高清| 亚洲精品乱码久久久v下载方式| 尾随美女入室| 久久久久久九九精品二区国产| 99热网站在线观看| 国产高清有码在线观看视频| 在线 av 中文字幕| 尤物成人国产欧美一区二区三区| 不卡视频在线观看欧美| 永久免费av网站大全| 久久久久久伊人网av| 亚洲,一卡二卡三卡| 在线免费观看不下载黄p国产| 亚洲av在线观看美女高潮| 免费看光身美女| 2022亚洲国产成人精品| 国产黄a三级三级三级人| 国产在视频线精品| 内地一区二区视频在线| 内射极品少妇av片p| 亚洲精品视频女| 国产男女内射视频| 2021少妇久久久久久久久久久| 久久久久久久久久久丰满| 亚洲欧美日韩东京热| 街头女战士在线观看网站| 91久久精品国产一区二区三区| 久久久久网色| 丝袜喷水一区| 三级男女做爰猛烈吃奶摸视频| 久久国产乱子免费精品| 中国国产av一级| 亚洲精品成人av观看孕妇| 美女被艹到高潮喷水动态| 高清午夜精品一区二区三区| 日韩制服骚丝袜av| 免费观看无遮挡的男女| 如何舔出高潮| 国产探花极品一区二区| 国产黄色免费在线视频| 国国产精品蜜臀av免费| 国产淫片久久久久久久久| 成人美女网站在线观看视频| 高清午夜精品一区二区三区| 日韩制服骚丝袜av| 少妇猛男粗大的猛烈进出视频 | 制服丝袜香蕉在线| 丰满人妻一区二区三区视频av| 王馨瑶露胸无遮挡在线观看| 亚洲久久久久久中文字幕| 九草在线视频观看| 久久99蜜桃精品久久| 偷拍熟女少妇极品色| 三级国产精品片| 在线播放无遮挡| 久久99精品国语久久久| 嫩草影院精品99| 国产精品国产三级国产专区5o| 18+在线观看网站| 亚洲精品自拍成人| 男女国产视频网站| 水蜜桃什么品种好| 午夜免费男女啪啪视频观看| 亚洲精品自拍成人| 国产精品爽爽va在线观看网站| 少妇 在线观看| 18+在线观看网站| 国产亚洲5aaaaa淫片| 国产精品女同一区二区软件| 干丝袜人妻中文字幕| 亚洲欧美日韩东京热| 精品国产乱码久久久久久小说| 久久久久久久精品精品| 18禁动态无遮挡网站| 白带黄色成豆腐渣| 久久久成人免费电影| 国产在视频线精品| 国产免费又黄又爽又色| 边亲边吃奶的免费视频| 亚洲内射少妇av| 久久久久国产精品人妻一区二区| 免费人成在线观看视频色| 国产爱豆传媒在线观看| 亚洲内射少妇av| 少妇人妻精品综合一区二区| 亚洲国产精品999| 日本色播在线视频| 免费观看性生交大片5| 熟妇人妻不卡中文字幕| 99久久九九国产精品国产免费| 精品国产三级普通话版| 18禁在线无遮挡免费观看视频| 老女人水多毛片| 国产精品一区二区性色av| 高清av免费在线| 51国产日韩欧美| 91精品一卡2卡3卡4卡| 亚洲精品自拍成人| 国产精品精品国产色婷婷| 久热久热在线精品观看| 99久久九九国产精品国产免费| 国产成人91sexporn| 热99国产精品久久久久久7| a级一级毛片免费在线观看| 只有这里有精品99| 免费观看无遮挡的男女| av在线app专区| 日本黄大片高清| 精品一区在线观看国产| 精品久久久久久电影网| 精品久久久久久久末码| 亚洲在线观看片| 亚洲一区二区三区欧美精品 | 精品久久久久久久人妻蜜臀av| 欧美亚洲 丝袜 人妻 在线| av在线播放精品| 亚洲天堂国产精品一区在线| 成人毛片60女人毛片免费| 欧美另类一区| 国产一级毛片在线| 网址你懂的国产日韩在线| 国产欧美另类精品又又久久亚洲欧美| 水蜜桃什么品种好| 高清视频免费观看一区二区| 精华霜和精华液先用哪个| 国模一区二区三区四区视频| 国产色婷婷99| 亚洲精品,欧美精品| 亚洲精品影视一区二区三区av| 国产日韩欧美在线精品| 汤姆久久久久久久影院中文字幕| 国产av国产精品国产| 涩涩av久久男人的天堂| 交换朋友夫妻互换小说| 一级毛片久久久久久久久女| 欧美成人午夜免费资源| 国产精品一区二区性色av| 欧美性感艳星| 99久久人妻综合| av网站免费在线观看视频| 久久久久久久精品精品| 老女人水多毛片| 婷婷色麻豆天堂久久| 久久韩国三级中文字幕| 日韩视频在线欧美| 日韩欧美精品v在线| 精品久久久精品久久久| 小蜜桃在线观看免费完整版高清| 精华霜和精华液先用哪个| 日本一二三区视频观看| 亚洲欧美成人精品一区二区| 日韩不卡一区二区三区视频在线| 国产精品嫩草影院av在线观看| 91狼人影院| 99久久精品热视频| 国产精品99久久99久久久不卡 | 亚洲国产成人一精品久久久| 欧美xxⅹ黑人| 国产综合精华液| 男女那种视频在线观看| 亚洲av国产av综合av卡| 好男人视频免费观看在线| a级一级毛片免费在线观看| 日韩伦理黄色片| 在线免费十八禁| 国产毛片a区久久久久| 国产av不卡久久| 欧美日韩亚洲高清精品| av在线蜜桃| 亚洲色图综合在线观看| 国产淫语在线视频| 日韩伦理黄色片| 国产精品一区www在线观看| 伊人久久精品亚洲午夜| 国产亚洲91精品色在线| 如何舔出高潮| 91久久精品电影网| 午夜福利视频精品| 赤兔流量卡办理| 精品一区二区免费观看| 久久精品国产亚洲网站| 国产熟女欧美一区二区| 一级毛片久久久久久久久女| tube8黄色片| 亚洲精品日韩av片在线观看| 老司机影院成人| 大香蕉久久网| 国产亚洲最大av| 国产综合懂色| 久久国产乱子免费精品| 久久久久网色| 伦理电影大哥的女人| 男女国产视频网站| 涩涩av久久男人的天堂| 深夜a级毛片| 国产精品一区www在线观看| 亚洲av二区三区四区| 国产精品精品国产色婷婷| 久久久精品94久久精品| 久久99热6这里只有精品| 精品一区在线观看国产| 老师上课跳d突然被开到最大视频| av在线观看视频网站免费| 男插女下体视频免费在线播放| 国产伦在线观看视频一区| 中文精品一卡2卡3卡4更新| 国产欧美日韩一区二区三区在线 | 黄片wwwwww| 亚洲成人久久爱视频| 色网站视频免费| 69人妻影院| 欧美高清成人免费视频www| 亚洲一级一片aⅴ在线观看| 如何舔出高潮| 午夜福利在线在线| 日韩亚洲欧美综合| 亚洲成色77777| 视频中文字幕在线观看| 亚洲丝袜综合中文字幕| 亚洲av男天堂| 香蕉精品网在线| 久久久亚洲精品成人影院| 男插女下体视频免费在线播放| 人妻制服诱惑在线中文字幕| av一本久久久久| 欧美一级a爱片免费观看看| 九九久久精品国产亚洲av麻豆| 成人鲁丝片一二三区免费| 国内少妇人妻偷人精品xxx网站| 性插视频无遮挡在线免费观看| 高清视频免费观看一区二区| 97在线视频观看| 国产色婷婷99| 一级毛片 在线播放| 国产欧美日韩一区二区三区在线 | 日韩制服骚丝袜av| 香蕉精品网在线| 免费看a级黄色片| 国产精品一区www在线观看| 一级毛片黄色毛片免费观看视频| 亚洲aⅴ乱码一区二区在线播放| 色5月婷婷丁香| 一级毛片电影观看| 亚洲三级黄色毛片| 精品视频人人做人人爽| 成人国产av品久久久| av在线亚洲专区| 夫妻午夜视频| 国产探花在线观看一区二区| 成人国产av品久久久| 国产有黄有色有爽视频| 毛片一级片免费看久久久久| 国产精品99久久久久久久久| av一本久久久久| 欧美潮喷喷水| 欧美精品国产亚洲| 国产综合懂色| 国产成人91sexporn| 亚洲美女搞黄在线观看| 婷婷色综合大香蕉| 久久99热6这里只有精品| 亚洲一区二区三区欧美精品 | 天堂中文最新版在线下载 | 别揉我奶头 嗯啊视频| 少妇被粗大猛烈的视频| 国产一级毛片在线| 亚洲成人精品中文字幕电影| 国产乱来视频区| 日产精品乱码卡一卡2卡三| 久久精品夜色国产| 日韩人妻高清精品专区| 男人爽女人下面视频在线观看| 成年女人看的毛片在线观看| 亚洲va在线va天堂va国产| 国产高清不卡午夜福利| 男男h啪啪无遮挡| 久久久久久久久久久免费av| 色播亚洲综合网| 看十八女毛片水多多多| 国产免费视频播放在线视频| 欧美激情久久久久久爽电影| 日本av手机在线免费观看| 搞女人的毛片| 最近中文字幕高清免费大全6| 精品少妇黑人巨大在线播放| av在线天堂中文字幕| 91精品一卡2卡3卡4卡| 青春草亚洲视频在线观看| 国产伦理片在线播放av一区| 我要看日韩黄色一级片| 久久这里有精品视频免费| 大片电影免费在线观看免费| 男女无遮挡免费网站观看| 日韩av免费高清视频| 日韩一区二区视频免费看| 一级a做视频免费观看| 国产在线男女| 国产精品爽爽va在线观看网站| 97超碰精品成人国产| 成人漫画全彩无遮挡| 中国美白少妇内射xxxbb| 激情五月婷婷亚洲| 国产精品不卡视频一区二区| 欧美一区二区亚洲| 国产一区二区三区综合在线观看 | 国内少妇人妻偷人精品xxx网站| 久久久色成人| 成人毛片60女人毛片免费| 人妻少妇偷人精品九色| 午夜日本视频在线| 韩国高清视频一区二区三区| 亚洲自拍偷在线| 亚洲av不卡在线观看| 日韩伦理黄色片| 亚洲av在线观看美女高潮| 中国三级夫妇交换| 欧美97在线视频| 91久久精品国产一区二区三区| 色5月婷婷丁香| av在线播放精品| av天堂中文字幕网| 禁无遮挡网站| 中文天堂在线官网| a级毛色黄片| av国产免费在线观看| 别揉我奶头 嗯啊视频| 成人国产av品久久久| 身体一侧抽搐| 最近的中文字幕免费完整| 国产精品偷伦视频观看了| 18禁动态无遮挡网站| 99热网站在线观看| 国产色婷婷99| 久久ye,这里只有精品| 搞女人的毛片| 欧美成人a在线观看| 日本-黄色视频高清免费观看| 国产精品一区www在线观看|