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

    三維數(shù)字圖像相關(guān)法的初值估計和仿真圖研究

    2017-04-05 07:53:08程百順毛建國梁中漢
    關(guān)鍵詞:散斑數(shù)字圖像初值

    程百順,毛建國,沈 峘,梁中漢,王 瑩

    (南京航空航天大學(xué) 能源與動力學(xué)院,南京 210016)

    三維數(shù)字圖像相關(guān)法的初值估計和仿真圖研究

    程百順,毛建國,沈 峘,梁中漢,王 瑩

    (南京航空航天大學(xué) 能源與動力學(xué)院,南京 210016)

    三維數(shù)字圖像相關(guān)法主要用來測量物體表面的三維形貌和變形,其關(guān)鍵技術(shù)主要包括相機(jī)標(biāo)定、圖像匹配、三維重建和三維位移場和應(yīng)變場的計算等。針對現(xiàn)有的初值估計方法在計算速度和精度上的局限性,提出基于sift的形函數(shù)和線性插值的初值估計方法,有效地提高了計算的速度和精度。傳統(tǒng)的三維算法精度分析只能借助于高精度的儀器測量對比,為此提出了基于反向映射法的三維仿真散斑圖分析方法。仿真結(jié)果表明:在小變形情況下z方向的誤差在2×10-4%左右,面內(nèi)位移測量誤差在8 μm以內(nèi)。驗證了三維仿真散斑圖的有效性和三維初值估計策略的可行性。

    三維數(shù)字圖像相關(guān)法;初值估計;仿真散斑圖

    在現(xiàn)代的科學(xué)研究、工程實踐及工業(yè)生產(chǎn)中,獲取材料的位移和應(yīng)變信息并對材料的變形及力學(xué)性能進(jìn)行深入研究具有重要的意義。傳統(tǒng)的接觸式測量方法和光學(xué)測量方法[1]由于自身的局限性不能滿足現(xiàn)在的測量要求。隨著計算機(jī)技術(shù)和相機(jī)產(chǎn)業(yè)的發(fā)展,在力學(xué)測量領(lǐng)域出現(xiàn)了新的測量方法,其中數(shù)字圖像相關(guān)法[2-3]是一種基于現(xiàn)代數(shù)字圖像和分析技術(shù)的新型光測技術(shù),通過分析變形前后物體表面的數(shù)字圖像獲得被測物體的變形信息,已經(jīng)成為現(xiàn)代光測力學(xué)領(lǐng)域中一個重要的方法。

    數(shù)字圖像相關(guān)法是由日本學(xué)者Yamaguchi[4]與美國學(xué)者Peters[5]等在20世紀(jì)80年代幾乎同時獨立提出的。在數(shù)字圖像相關(guān)法中相關(guān)匹配運算一直是該算法的核心,其重要的進(jìn)展是1989年Bruck[6]等提出的Newton-Rapshon優(yōu)化算法,該方法不僅精度高,而且只要變形初值估計合理,算法會迅速收斂,現(xiàn)在的大多數(shù)研究都采用此種算法。但是二維數(shù)字圖像相關(guān)法只能測量平面物體的面內(nèi)位移,所以引進(jìn)三維數(shù)字圖像相關(guān)法。在實驗力學(xué)領(lǐng)域,Luo和Sutton等[7]利用2個攝像機(jī)從不同角度對被測物體表面進(jìn)行成像,首先對雙目立體視覺模型進(jìn)行標(biāo)定獲得2個攝像機(jī)的內(nèi)外參數(shù),然后直接利用二維數(shù)字圖像相關(guān)法中的相關(guān)匹配算法得到左右兩幅圖像中對應(yīng)點的視差,從各點的視差數(shù)據(jù)和預(yù)先獲得的標(biāo)定參數(shù)恢復(fù)物體表面的三維形貌,從而得到全場的位移分布。

    隨著現(xiàn)代數(shù)學(xué)和光學(xué)技術(shù)的發(fā)展,越來越多的國內(nèi)外學(xué)者在數(shù)字圖像相關(guān)法的基礎(chǔ)理論[8-9]、匹配精度[10]和計算精度[11]上作出巨大的貢獻(xiàn)。其中初值估計是三維數(shù)字圖像相關(guān)法(3D-DIC)中的基礎(chǔ)環(huán)節(jié),目前主要方法有整像素估計法[12]和種子點法[13-14],但其在計算精度和速度上有缺陷。本文基于sift[15]的初值估計方法,得到精度較高的匹配點,然后利用形函數(shù)和線性插值得到所有關(guān)心點的估計點。本文在二維數(shù)字圖像相關(guān)法的散斑技術(shù)[16]基礎(chǔ)上,提出建立三維仿真散斑圖的方法,解決現(xiàn)有的利用仿真實驗驗證算法的問題。

    1 三維數(shù)字圖像相關(guān)法

    三維數(shù)字圖像相關(guān)法的主要技術(shù)包括相機(jī)標(biāo)定、圖像匹配、位移場和應(yīng)變場的計算,其中圖像匹配包括粗匹配和細(xì)匹配。粗匹配是指感興趣點的初值估計,細(xì)匹配是指相關(guān)函數(shù)的迭代收斂計算。應(yīng)變場的計算采用將文獻(xiàn)[17]中的三維應(yīng)變計算和潘兵最小二乘[18]計算相結(jié)合的方法。

    1.1 相機(jī)標(biāo)定

    相機(jī)標(biāo)定是3D-DIC的基礎(chǔ)。本文采用的張正友[19]的標(biāo)定方法。利用采集平面標(biāo)定板上的標(biāo)記點在空間中的不同位置反向計算出相機(jī)的內(nèi)參數(shù)和外參數(shù)。其中的小孔成像模型公式為:

    (1)

    式中:s為任意的比例系數(shù);R和t為外部參數(shù),分別為從世界坐標(biāo)系到攝像機(jī)坐標(biāo)系的旋轉(zhuǎn)矩陣和平移向量;A為攝像機(jī)的內(nèi)部參數(shù)。由于簡單的針孔模型不能準(zhǔn)確描述鏡頭的成像關(guān)系,本文采用tsai兩步標(biāo)定法[20],公式如下:

    (2)

    其中:k1、k2分別為1階、2階徑向畸變系數(shù);r為像素點到主點(u,v)的距離。

    1.2 圖像匹配

    圖像匹配現(xiàn)在的主流方法是基于NR迭代的相關(guān)函數(shù)方法。物體表面在發(fā)生變形,左相機(jī)所采的系列圖像相當(dāng)于在二維平面上的變形,而對應(yīng)的右相機(jī)所采圖像就是因為視差所導(dǎo)致的類變形圖像,如圖1所示。

    圖1 物體變形前后示意圖

    圖1中:p點是計算中心點;(x,y)是p點的圖像坐標(biāo);(x′,y′)是變形后的對應(yīng)點。(x,y)就是所要計算的位移。采用歸一化最小平方距離相關(guān)函數(shù)評價子區(qū)的相似度。由于圖像整體存在變形,因此采用形函數(shù)來描述這種變形,用于同一組相機(jī)所采圖像的匹配和兩個相機(jī)之間的匹配。

    (1)

    其中:Δx=x′-x0,Δy=y′-y0;ξ0、η0是參考圖像子區(qū)中心點P(x0,y0)的視差;ξx、ξy、ηx、ηy為圖像子區(qū)內(nèi)視差的1階導(dǎo)數(shù);ξxx、ξxy、ξyy、ηxx、ηxy、ηyy為圖像子區(qū)內(nèi)視差的2階導(dǎo)數(shù)。

    2 初值估計

    2.1 初值估計基本理論

    初值估計在數(shù)學(xué)上的意義是給上述形函數(shù)初值的過程,在實驗中是要賦予計算的點的初始位置。由于sift可能不會提取所有的特征點,此時應(yīng)用形函數(shù)和線性插值可以解決此類問題。網(wǎng)格點只是從所有的像素點中提取的一部分,而且是整像素的網(wǎng)格點,此時就要提取精度較高的sift匹配點,同時去除誤匹配點。

    圖2中:p點是參考圖像中的計算點;p′是p點在目標(biāo)圖像中的初值估計點;p1,…,p5是以p點為中心,寬W、高H范圍內(nèi)的sift提取的特征點;p1′,…,p5′是p1,…,p5在目標(biāo)圖像中的對應(yīng)sift特征點。

    圖2 設(shè)定區(qū)域內(nèi)sift匹配結(jié)果

    1) 若在計算點周圍存在大于k(k=3)個sift特征點,那么就認(rèn)為p點可以被估計。

    2) 計算sift特征點到p0點的距離d。sift的匹配有一定的不確定性,取di為中值的點作為計算參考點。

    3) 利用線性插值公式

    (2)

    計算得到p0點的估計點p′0 的坐標(biāo)(xp0′,yp0′)。但并不是每個計算點都能滿足sift估計的要求,所以引進(jìn)形函數(shù)初值估計。在將已經(jīng)估計出來的計算點進(jìn)行N-R迭代計算,若是收斂成功,則利用式(3)將其周圍某一點的初值計算出來了。若某一估值點計算不收斂,就會導(dǎo)致其某一側(cè)的值不能通過形函數(shù)來估計,此時引入線性插值式(4)來計算,如圖3所示。其中,“+”為N-R計算收斂點,“△”為計算不收斂的點,“o”是利用線性插值計算的點。

    圖3

    2.2 實驗分析

    現(xiàn)在常用的初值估計方法有2種:一種是基于人機(jī)交互的種子點法;一種是結(jié)合模板匹配的整像素法?;谌藱C(jī)交互種子點法是在目標(biāo)圖像中,利用人工選取點作為初值估計點,然后進(jìn)行匹配計算,通過這一點的收斂結(jié)果估計其周圍點的初值,然后幅散開來,完成整張圖像的匹配工作。結(jié)合模板匹配的整像素估計是利用參考圖像子區(qū)和目標(biāo)圖像子區(qū)進(jìn)行模板匹配,找到初值估計點,然后進(jìn)行匹配計算,這是一種由粗匹配到細(xì)匹配的過程。

    本實驗是利用仿真實驗完成,硬件配置如下:采用Intel(R) Core(TM) i7-6700 CPU,內(nèi)存為16G。軟件采用Matlab。仿真圖為旋轉(zhuǎn)5°、15°、30°的散斑圖。其中計算點為336個,網(wǎng)格間距20個像素,計算子區(qū)半徑20個像素。如圖4所示,其中第1列是二維仿真圖,第2列是三維仿真散斑圖的左相機(jī)圖像,第3列是三維仿真散斑圖的右相機(jī)圖像,生成圖像具體參數(shù)參考本文第4節(jié)。

    實驗結(jié)果如表1、2所示,其中sift的插值方法在計算時間上比種子點法和整像素法要少,而且在大旋轉(zhuǎn)的情況下,種子點法和整像素法會因為無法完成初值估計導(dǎo)致計算不能收斂,而基于sift的插值計算方法使得每個點計算的都收斂。

    種子點法的缺陷是:即使在選對第一個點作為種子點的情況下,若該點不收斂那么無法進(jìn)行進(jìn)一步的計算,而且是靠人工選取,時間較長。整像素法初值估計主要用基于模版匹配的方法,其抗旋轉(zhuǎn)能力較差,所以在大變形的情況下不能滿足計算要求。

    旋轉(zhuǎn)角度角度/(°)51530sift法收斂點數(shù)336336301計算時間/s19.3023.71137.54種子點法收斂點數(shù)336336不收斂計算時間/s37.0078.45不收斂整像素法收斂點數(shù)336336不收斂計算時間/s63.16111.12不收斂

    表2 右相機(jī)圖片計算結(jié)果

    3 基于反向映射法的三維仿真散斑圖

    3.1 散斑圖生成原理

    散斑圖生成公式如下:

    (3)

    式中:k表示散斑顆粒的數(shù)目;(xk,yk)為散斑顆粒的中心位置,隨機(jī)產(chǎn)生;R為散斑顆粒的半徑;Ik為隨機(jī)生成的散斑中心的亮度,范圍為0~255, 圖像大小為M×N。當(dāng)生成變形圖像時,散斑顆粒的中心點發(fā)生偏移,滿足式(6)中的變形模式。

    (4)

    其中(u0,v0,ux,vx,uy,vy)為變形參數(shù)。由(xk,yk)可求得(xk′,yk′),代入式(4)可得到變形散斑圖的灰度分布。

    3.2 三維反向映射法

    本文生成三維仿真圖的方法如下:

    物體是一個平面,ZW=1 200 mm,物體沿x反向平移0.2個像素:

    1) 設(shè)定相機(jī)標(biāo)定參數(shù),主要是設(shè)定每個相機(jī)的內(nèi)參數(shù)A和左右相機(jī)之間的外參數(shù)(R|T),設(shè)定左相機(jī)的相機(jī)坐標(biāo)為世界坐標(biāo)。

    2) 利用式(5)生成散斑圖,作為左相機(jī)的第1張圖像及參考圖像,利用反向映射法生成沿x方向平移0.2個像素的目標(biāo)圖像,作為左相機(jī)的第2幅圖像。

    3) 利用先前的相機(jī)內(nèi)外參數(shù)和左相機(jī)第1幅圖中的像素點的位置(UL,VL),結(jié)合式(1),得到左相機(jī)中的像素點在右相機(jī)圖像中的坐標(biāo)位置(UR,VR),將左相機(jī)中對應(yīng)像素點的灰度值賦給右相機(jī)圖像中的對應(yīng)點。

    4 實驗驗證

    參考圖像為隨機(jī)生成的散斑圖像,分辨率為700像素×700 像素,散斑顆粒數(shù)為19 600,散斑顆粒半徑為3。假設(shè)物體是一個ZW=1 200 mm的平面。其中左右相機(jī)的內(nèi)外參數(shù)為:

    4.1 位移測試

    沿x向生成間隔為0.05像素的平移圖像,共20張,在三維空間里對應(yīng)的實際位移從0到0.12 mm,間隔為0.006 mm。實驗結(jié)果如圖5~7所示。

    圖5 x方向的位移場誤差

    圖6 y反向位移場誤差

    如圖5和6所示,模擬0~0.12 mm的位移測量中測量誤差在8 μm以內(nèi)。如圖7所示,在z方向的距離測量誤差在1 μm以內(nèi),相對誤差為2×10-4%。結(jié)果表明:在小變形的情況下,三維仿真散斑圖的生成是有效的,而且精度較高,滿足測試三維形貌重建的算法測試標(biāo)準(zhǔn)。

    圖7 z方向距離測量誤差

    4.1 應(yīng)變測試

    沿x方向有應(yīng)變2 000 με。參考圖像和目標(biāo)圖像如圖8所示。(a)和(b)是左右相機(jī)采集的第1對圖像,(c)和(d)左右相機(jī)采集的第2對圖像。

    圖8 參考圖像和目標(biāo)圖像

    如圖9所示,x方向的應(yīng)變場計算結(jié)果為1 960 με,方差為58 με。如圖10所示,z反向的距離測量誤差均值為0.2 μm,最大值為2 μm。因此,三維仿真散斑圖可用于三維重建和應(yīng)變的精度分析,驗證應(yīng)變和重建算法的可行性。

    圖9 x方向應(yīng)變場

    圖10 z方向尺寸測量誤差

    5 結(jié)束語

    根據(jù)實驗結(jié)果可知:本文采用的三維仿真散斑圖可以有效驗證立體匹配、三維重建和應(yīng)變算法的精度和有效性,同時本文的三維初值估計及其插值算法滿足三維數(shù)字圖像相關(guān)法的精度要求。

    [1] 金觀昌.計算機(jī)輔助光學(xué)測量[M].北京:清華大學(xué)出版社,2007.

    JIN Guanchang.Computer aided optical measurement[M].Beijing:Tsinghua University Press,2007.

    [2] 王懷文,亢一瀾,謝和平.數(shù)字散斑相關(guān)方法與應(yīng)用研究進(jìn)展[J].力學(xué)進(jìn)展,2005(2):195-203.

    WANG Huaiwen,KANG Yilan,XIE Heping.Advance in digital speckle correlation method and its application[J].Advances in Mechanics,2005(2):195-203.

    [3] HUNG P C,VOLOSHIN A S.In-plane strain measurement by digital image correlation[J].Journal of the Brazilian Society of Mechanical Sciences & Engineering,2003,25(3):215-221.

    [4] YAMAGUCHI I.A laser-speckle strain gauge[J].Phys E:Sci Instrum.1981,14(11):1270-1273.

    [5] PETERS W H,RANSON W F.Digital imaging techniques on experimental stress analysis[J].Optical Engineering 1982,21(3):427-431.

    [6] BRUCK H A,MCNEILL S R,SUTTON M A,et al.Digital image correlation using Newton-Raphson method of partial differential correction[J].Experimental Mechanics,1989,29(29):261-267.

    [7] LUO P F,CHAO Y J,SUTTON M A.Application of stereo vision to three-dimensional deformation analyses in fracture experiments[J].Optical Engineering,1994,33(3):990.

    [8] SUTTON M A,TURNER J L,BRUCK H A,et al.Full-field representation of discretely sampled surface deformation for displacement and strain analysis[J].Experimental Mechanics,1991,31(31):168-177.

    [9] ORTEU J J.3-D computer vision in experimental mechanics[J].Optics & Lasers in Engineering,2009,47(3/4):282-291.

    [10]SCHREIER H,ORTEU J J,SUTTON M A.Image Correlation for Shape,Motion and Deformation Measurements[M].[S.l.]:Basic Concepts,Theory and Applications,2009:565-600.

    [11]PAN B,YUAN J,XIA Y.Strain field denoising for digital image correlation using a regularized cost-function[J].Optics & Lasers in Engineering,2015,65:9-17.

    [12]GAO Y,CHENG T,SU Y,et al.High-efficiency and high-accuracy digital image correlation for three-dimensional measurement[J].Optics & Lasers in Engineering,2015,65:73-80.

    [13]PAN B,ASUNDI A,XIE H,et al.Digital image correlation using iterative least squares and pointwise least squares for displacement field and strain field measurements[J].Optics & Lasers in Engineering,2009,47(7/8):865-874.

    [14]HELM J D,MCNEILL S R,SUTTON M A.Improved three‐dimensional image correlation for surface displacement measurement[J].Proceedings of SPIE-The International Society for Optical Engineering,1994,35(7):1911-1920.

    [15]ZHOU Y,CHEN Y Q.Feature matching for automated and reliable initialization in three-dimensional digital image correlation[J].Optics & Lasers in Engineering,2013,51(3):213-223.

    [16]MAO Jianguo,ZHANG Peize,SHEN Huan,et al,Reverse Mapping for Generating Simulated Deformed Speckle Patterns[J].Journal of Optoelectronics.Laser,2015(4):2433-2439.

    [17]高越.三維數(shù)字圖像相關(guān)法的關(guān)鍵技術(shù)及應(yīng)用研究[D].合肥:中國科學(xué)技術(shù)大學(xué),2014.

    GAO Yue.Research on key technologies and applications of three-dimensional digital image correlation[D].Hefei:University of science and Technology of China,2014.

    [18]潘兵,謝惠民.數(shù)字圖像相關(guān)中基于位移場局部最小二乘擬合的全場應(yīng)變測量[J].光學(xué)學(xué)報,2007,27(11):1980-1986.

    PAN Bing,XIE Huimin.Full-Filed strain measurement based on Least-Square fitting of displacement for digital image correlation[J].Acta Optica Sinica,2007,27(11):1980-1986.

    [19]ZHANG.A flexible new technique for camera calibration[J].IEEE Transactions on Pattern Analysis & Machine Intelligence,2000,22(11):1330-1334.

    [20]TSAI R.A versatile camera calibration technique for high-accuracy 3D machine vision metrology using off-the-shelf TV cameras and lenses[J].Robotics & Automation IEEE Journal,1987,3(4):323-344.

    (責(zé)任編輯 楊黎麗)

    Study on the Initialization for Feature Matching and Simulation Speckles in Three-Dimensional Digital Image Correlation

    CHENG Bai-shun, MAO Jian-guo, SHEN Huan, LIANG Zhong-han, WANG Ying

    (Energy and Power College, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China)

    Three-dimensional digital image correlation (3D-DIC) method has been increasingly employed to measure 3D surface profile and displacement. The key technologies applied in this method include stereo camera calibration, digital image correlation, reconstruction, and 3D displacement/strain field computation. This paper mainly studies on an initial guess technique based on SIFT matching and linear interpolation to improve the accuracy of 3D-DIC. Additionally, a new method based on reverse mapping is proposed to generate the 3D simulated speckles for testing the accuracy of the initial guess and digital image correlation technique. The numerical simulations reveal that the error in Z-direction is about 2×10-4% of the object distance. Test and analysis of the in-plane displacement error is within 8 μm under small displacement. Results show that the 3D simulated speckles are useful and effective for the verification of 3D-DIC method and show that the initial guess technique can achieve high accuracy.Key words: three-dimensional digital image correlation; initial guess technique; 3D simulated speckles

    2016-11-03 基金項目:航空科學(xué)基金資助項目(111111)

    程百順(1991—),江蘇徐州人,碩士研究生,主要從事機(jī)器視覺,數(shù)字圖像相關(guān)法在力學(xué)方向的研究,E-mail:cbs_1012@163.com。

    程百順,毛建國,沈峘,等.三維數(shù)字圖像相關(guān)法的初值估計和仿真圖研究[J].重慶理工大學(xué)學(xué)報(自然科學(xué)),2017(3):116-122.

    format:CHENG Bai-shun, MAO Jian-guo, SHEN Huan, et al.Study on the Initialization for Feature Matching and Simulation Speckles in Three-Dimensional Digital Image Correlation[J].Journal of Chongqing University of Technology(Natural Science),2017(3):116-122.

    10.3969/j.issn.1674-8425(z).2017.03.017

    TP391.41

    A

    1674-8425(2017)03-0116-07

    猜你喜歡
    散斑數(shù)字圖像初值
    具非定常數(shù)初值的全變差方程解的漸近性
    激光顯示中的彩色散斑測量研究
    一種適用于平動點周期軌道初值計算的簡化路徑搜索修正法
    激光投影顯示散斑抑制方法研究
    三維擬線性波方程的小初值光滑解
    ARGUS-100 藝術(shù)品鑒證數(shù)字圖像比對系統(tǒng)
    用于檢驗散斑協(xié)方差矩陣估計性能的白化度評價方法
    基于塊效應(yīng)測度的JPEG數(shù)字圖像盲取證
    電鏡成像方式對數(shù)字散斑相關(guān)方法結(jié)果的影響
    數(shù)字圖像修復(fù)在圖像壓縮上的應(yīng)用
    欧美久久黑人一区二区| 日韩 欧美 亚洲 中文字幕| 蜜桃在线观看..| 精品久久久久久久毛片微露脸 | 亚洲国产av新网站| 真人做人爱边吃奶动态| 国产97色在线日韩免费| 99热国产这里只有精品6| 国产激情久久老熟女| 色播在线永久视频| 人人妻,人人澡人人爽秒播 | 性少妇av在线| 精品国产国语对白av| 婷婷成人精品国产| 十八禁网站网址无遮挡| 午夜av观看不卡| 多毛熟女@视频| 中文字幕精品免费在线观看视频| 国产日韩欧美亚洲二区| 90打野战视频偷拍视频| 亚洲国产av新网站| 人妻人人澡人人爽人人| 极品少妇高潮喷水抽搐| 肉色欧美久久久久久久蜜桃| 在线av久久热| av不卡在线播放| 18禁观看日本| 90打野战视频偷拍视频| 女人精品久久久久毛片| 一本大道久久a久久精品| 国产深夜福利视频在线观看| 夫妻性生交免费视频一级片| e午夜精品久久久久久久| 高清不卡的av网站| av线在线观看网站| 欧美激情高清一区二区三区| 日本vs欧美在线观看视频| 国产1区2区3区精品| 妹子高潮喷水视频| 黄色视频不卡| 亚洲欧美成人综合另类久久久| 精品少妇久久久久久888优播| 秋霞在线观看毛片| 亚洲,欧美精品.| 日韩一本色道免费dvd| 亚洲精品国产一区二区精华液| 高清欧美精品videossex| 极品少妇高潮喷水抽搐| 美女脱内裤让男人舔精品视频| 女警被强在线播放| 99久久99久久久精品蜜桃| 亚洲国产欧美日韩在线播放| 欧美+亚洲+日韩+国产| 在线观看人妻少妇| 美女脱内裤让男人舔精品视频| e午夜精品久久久久久久| 99久久99久久久精品蜜桃| 黑人巨大精品欧美一区二区蜜桃| 天天添夜夜摸| 爱豆传媒免费全集在线观看| 久久久久网色| av视频免费观看在线观看| 波野结衣二区三区在线| 欧美黄色淫秽网站| 亚洲午夜精品一区,二区,三区| tube8黄色片| 美女脱内裤让男人舔精品视频| 日韩一卡2卡3卡4卡2021年| 99久久99久久久精品蜜桃| 黑人巨大精品欧美一区二区蜜桃| 国产伦人伦偷精品视频| 美女扒开内裤让男人捅视频| 日韩一区二区三区影片| 啦啦啦 在线观看视频| 欧美日韩一级在线毛片| 咕卡用的链子| 国产免费又黄又爽又色| 欧美在线一区亚洲| 青春草亚洲视频在线观看| 久热爱精品视频在线9| 成人午夜精彩视频在线观看| 免费人妻精品一区二区三区视频| 日日夜夜操网爽| 美女视频免费永久观看网站| 女人久久www免费人成看片| 久久精品久久久久久噜噜老黄| 久久青草综合色| 999精品在线视频| 国产欧美日韩精品亚洲av| 国产精品一国产av| xxxhd国产人妻xxx| 欧美成人精品欧美一级黄| 欧美日本中文国产一区发布| 操美女的视频在线观看| 下体分泌物呈黄色| 人妻一区二区av| 国产精品成人在线| 国产片内射在线| 午夜激情久久久久久久| 日韩视频在线欧美| 亚洲男人天堂网一区| 精品久久久久久久毛片微露脸 | 99久久人妻综合| tube8黄色片| 国产成人精品久久久久久| 亚洲欧洲国产日韩| av不卡在线播放| avwww免费| 一级片免费观看大全| 多毛熟女@视频| 亚洲欧美日韩另类电影网站| 黄片小视频在线播放| 这个男人来自地球电影免费观看| 国产精品三级大全| 女人被躁到高潮嗷嗷叫费观| 人体艺术视频欧美日本| 国产成人啪精品午夜网站| 日韩 亚洲 欧美在线| 亚洲国产精品国产精品| 国产在线视频一区二区| 亚洲伊人久久精品综合| xxxhd国产人妻xxx| 嫩草影视91久久| 国产成人精品在线电影| 久久久国产一区二区| 日韩 亚洲 欧美在线| 亚洲第一青青草原| 老司机午夜十八禁免费视频| 热99国产精品久久久久久7| 国产一区有黄有色的免费视频| 满18在线观看网站| 久久久久久久久免费视频了| 久久精品久久久久久久性| 国产精品久久久久久精品古装| 高清av免费在线| 两个人看的免费小视频| 欧美人与性动交α欧美软件| 一二三四在线观看免费中文在| 午夜福利乱码中文字幕| 最黄视频免费看| 国产免费一区二区三区四区乱码| 午夜免费观看性视频| 亚洲五月婷婷丁香| 国产免费视频播放在线视频| 亚洲少妇的诱惑av| 一级,二级,三级黄色视频| √禁漫天堂资源中文www| 欧美成人午夜精品| 叶爱在线成人免费视频播放| 久久热在线av| 乱人伦中国视频| 最近最新中文字幕大全免费视频 | 国产真人三级小视频在线观看| 母亲3免费完整高清在线观看| 亚洲国产精品一区三区| av天堂在线播放| 久久人妻熟女aⅴ| 午夜久久久在线观看| 国产日韩欧美视频二区| 丝袜美腿诱惑在线| 亚洲av电影在线观看一区二区三区| 悠悠久久av| 99精国产麻豆久久婷婷| 国产视频首页在线观看| 乱人伦中国视频| 最近最新中文字幕大全免费视频 | 国产亚洲欧美在线一区二区| 午夜激情久久久久久久| 每晚都被弄得嗷嗷叫到高潮| 久久中文字幕一级| 黄色a级毛片大全视频| 欧美日韩亚洲高清精品| 亚洲欧美精品自产自拍| 久久99热这里只频精品6学生| av线在线观看网站| 亚洲精品久久久久久婷婷小说| 啦啦啦在线免费观看视频4| 老熟女久久久| 欧美xxⅹ黑人| 精品久久久久久电影网| 美女大奶头黄色视频| 久久国产精品大桥未久av| 一区在线观看完整版| 久久国产精品人妻蜜桃| 超色免费av| 99热国产这里只有精品6| 久久热在线av| 老司机深夜福利视频在线观看 | 婷婷色综合www| 日韩视频在线欧美| 91国产中文字幕| 中文字幕高清在线视频| 亚洲欧美色中文字幕在线| 丝袜人妻中文字幕| 好男人视频免费观看在线| 美女午夜性视频免费| 丝袜人妻中文字幕| 91九色精品人成在线观看| 1024香蕉在线观看| 成人三级做爰电影| 在线亚洲精品国产二区图片欧美| 免费女性裸体啪啪无遮挡网站| 国产精品一区二区在线不卡| 搡老岳熟女国产| 色精品久久人妻99蜜桃| 男女边吃奶边做爰视频| 母亲3免费完整高清在线观看| 91国产中文字幕| 欧美黄色淫秽网站| 91九色精品人成在线观看| 老汉色∧v一级毛片| 亚洲人成网站在线观看播放| 丝袜美腿诱惑在线| 欧美日韩福利视频一区二区| 国产精品欧美亚洲77777| 国产日韩欧美视频二区| 亚洲欧美一区二区三区黑人| 黄色毛片三级朝国网站| 午夜福利一区二区在线看| 日韩伦理黄色片| 国产精品熟女久久久久浪| 欧美黄色片欧美黄色片| 国产男女内射视频| av天堂在线播放| 免费高清在线观看视频在线观看| 欧美精品啪啪一区二区三区 | 欧美黄色片欧美黄色片| 欧美中文综合在线视频| videos熟女内射| 两性夫妻黄色片| 久久精品亚洲熟妇少妇任你| 99re6热这里在线精品视频| 婷婷丁香在线五月| 纯流量卡能插随身wifi吗| 纵有疾风起免费观看全集完整版| 国产成人av激情在线播放| 老司机靠b影院| 国产精品久久久av美女十八| 久久精品国产亚洲av涩爱| 日韩大码丰满熟妇| 欧美国产精品一级二级三级| 欧美 亚洲 国产 日韩一| 在线亚洲精品国产二区图片欧美| 久久中文字幕一级| 男男h啪啪无遮挡| 热99久久久久精品小说推荐| 91麻豆av在线| 免费黄频网站在线观看国产| 精品一区二区三区四区五区乱码 | 亚洲人成77777在线视频| av一本久久久久| av线在线观看网站| 欧美成狂野欧美在线观看| 高清欧美精品videossex| h视频一区二区三区| 日韩av不卡免费在线播放| 91麻豆av在线| 国产精品久久久久久精品古装| 天天操日日干夜夜撸| 国产免费视频播放在线视频| 日本欧美视频一区| 国产成人一区二区三区免费视频网站 | 成人18禁高潮啪啪吃奶动态图| 午夜免费成人在线视频| 国产精品 国内视频| 最近手机中文字幕大全| 大香蕉久久网| 欧美日韩亚洲国产一区二区在线观看 | 1024香蕉在线观看| 日韩伦理黄色片| 午夜两性在线视频| 国产精品九九99| 精品一区在线观看国产| 国产成人系列免费观看| 欧美日韩成人在线一区二区| 亚洲精品第二区| 大话2 男鬼变身卡| av一本久久久久| 国产在线观看jvid| 久久精品久久久久久噜噜老黄| 最近手机中文字幕大全| 久久国产精品大桥未久av| 欧美日韩视频高清一区二区三区二| 嫩草影视91久久| 香蕉国产在线看| 午夜福利视频精品| 国产欧美亚洲国产| 50天的宝宝边吃奶边哭怎么回事| 国产野战对白在线观看| 99久久精品国产亚洲精品| 亚洲精品中文字幕在线视频| 亚洲欧美一区二区三区黑人| 亚洲自偷自拍图片 自拍| 国产精品.久久久| 亚洲成国产人片在线观看| 欧美日韩亚洲高清精品| 婷婷色综合大香蕉| 亚洲专区国产一区二区| 一级,二级,三级黄色视频| 欧美精品啪啪一区二区三区 | 成年人免费黄色播放视频| 成人亚洲欧美一区二区av| 成人国产av品久久久| 午夜免费男女啪啪视频观看| 精品熟女少妇八av免费久了| 国产精品久久久久久精品古装| cao死你这个sao货| 真人做人爱边吃奶动态| 成人免费观看视频高清| 人人澡人人妻人| 精品亚洲乱码少妇综合久久| 国产成人一区二区三区免费视频网站 | 亚洲欧美激情在线| 午夜免费男女啪啪视频观看| 亚洲精品国产色婷婷电影| 美女中出高潮动态图| 亚洲精品日韩在线中文字幕| 国产在线视频一区二区| 69精品国产乱码久久久| 亚洲精品国产区一区二| 亚洲精品久久久久久婷婷小说| 国产成人a∨麻豆精品| 中文字幕人妻丝袜一区二区| 欧美日韩亚洲国产一区二区在线观看 | 久久久久久久久久久久大奶| 男女床上黄色一级片免费看| 深夜精品福利| 在线看a的网站| 十分钟在线观看高清视频www| 日日夜夜操网爽| 丝袜喷水一区| 亚洲国产日韩一区二区| 99精品久久久久人妻精品| 日本vs欧美在线观看视频| 啦啦啦视频在线资源免费观看| 成年人午夜在线观看视频| 99久久人妻综合| 午夜福利免费观看在线| av网站在线播放免费| 国产黄色视频一区二区在线观看| 亚洲av美国av| 香蕉国产在线看| 一区二区三区乱码不卡18| 亚洲国产精品国产精品| av网站免费在线观看视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品一二三区在线看| 国产高清不卡午夜福利| 可以免费在线观看a视频的电影网站| 亚洲精品成人av观看孕妇| 18禁国产床啪视频网站| √禁漫天堂资源中文www| 亚洲专区中文字幕在线| 老司机午夜十八禁免费视频| 欧美日韩黄片免| 成年人黄色毛片网站| 妹子高潮喷水视频| 美女大奶头黄色视频| a 毛片基地| 精品少妇黑人巨大在线播放| 免费一级毛片在线播放高清视频 | 国产一区二区 视频在线| 国产欧美日韩一区二区三 | 可以免费在线观看a视频的电影网站| 亚洲成人免费电影在线观看 | 飞空精品影院首页| 男女下面插进去视频免费观看| 国产老妇伦熟女老妇高清| 久久人妻熟女aⅴ| 欧美日韩视频高清一区二区三区二| 欧美黄色片欧美黄色片| 欧美日韩视频高清一区二区三区二| 国产精品久久久av美女十八| 日韩制服骚丝袜av| 欧美亚洲日本最大视频资源| 视频区图区小说| 99久久人妻综合| 1024视频免费在线观看| 韩国高清视频一区二区三区| 欧美成狂野欧美在线观看| 婷婷色综合www| 51午夜福利影视在线观看| 国产精品久久久久成人av| 亚洲,欧美,日韩| 久久综合国产亚洲精品| 国产一区二区激情短视频 | 热99国产精品久久久久久7| 少妇被粗大的猛进出69影院| 两人在一起打扑克的视频| 无遮挡黄片免费观看| 嫁个100分男人电影在线观看 | 91国产中文字幕| 又黄又粗又硬又大视频| 亚洲一码二码三码区别大吗| 国产黄色视频一区二区在线观看| 国产亚洲精品久久久久5区| 91九色精品人成在线观看| 狠狠婷婷综合久久久久久88av| 国产成人精品久久二区二区免费| 日韩,欧美,国产一区二区三区| 久久99热这里只频精品6学生| 久久久久精品人妻al黑| 午夜免费鲁丝| 国产在线视频一区二区| 午夜福利视频精品| 人妻一区二区av| av国产精品久久久久影院| www.av在线官网国产| av天堂久久9| 国产精品二区激情视频| 啦啦啦视频在线资源免费观看| 蜜桃在线观看..| 国产成人av激情在线播放| www.av在线官网国产| 在线观看国产h片| 在线精品无人区一区二区三| 丝袜美足系列| 在线观看人妻少妇| 自拍欧美九色日韩亚洲蝌蚪91| 久久狼人影院| 桃花免费在线播放| 男人添女人高潮全过程视频| 久久久久精品国产欧美久久久 | 两个人免费观看高清视频| 欧美少妇被猛烈插入视频| 精品亚洲成a人片在线观看| 一本色道久久久久久精品综合| 免费黄频网站在线观看国产| 两人在一起打扑克的视频| 欧美在线黄色| 夫妻午夜视频| 国产日韩欧美在线精品| 在线观看国产h片| 国产欧美亚洲国产| 高潮久久久久久久久久久不卡| 国产欧美亚洲国产| 亚洲欧美清纯卡通| 日韩熟女老妇一区二区性免费视频| 一级毛片我不卡| 在线看a的网站| 久久 成人 亚洲| 捣出白浆h1v1| 日韩免费高清中文字幕av| 欧美在线一区亚洲| 老司机影院成人| 日韩大码丰满熟妇| 91字幕亚洲| 最新的欧美精品一区二区| 超碰97精品在线观看| 欧美人与性动交α欧美精品济南到| 日韩 亚洲 欧美在线| 久久性视频一级片| 三上悠亚av全集在线观看| 免费不卡黄色视频| 国产在线免费精品| 久久ye,这里只有精品| 人人妻人人爽人人添夜夜欢视频| 搡老乐熟女国产| 国产欧美亚洲国产| 99国产精品一区二区三区| 国产男女超爽视频在线观看| 欧美 亚洲 国产 日韩一| 麻豆av在线久日| 国产精品久久久久久精品古装| 又大又爽又粗| 亚洲图色成人| 午夜激情久久久久久久| 丝袜喷水一区| 在线观看国产h片| 看十八女毛片水多多多| 男女国产视频网站| 成年动漫av网址| 精品第一国产精品| 国产成人精品在线电影| 精品少妇内射三级| 亚洲人成网站在线观看播放| 制服人妻中文乱码| 国产在线视频一区二区| 高清不卡的av网站| 久久99一区二区三区| 一边摸一边做爽爽视频免费| 精品免费久久久久久久清纯 | 久久久国产精品麻豆| 女人爽到高潮嗷嗷叫在线视频| 国产一区二区激情短视频 | 亚洲图色成人| 国产精品一国产av| 国产精品麻豆人妻色哟哟久久| 国产精品欧美亚洲77777| 久久狼人影院| 两个人看的免费小视频| 99国产综合亚洲精品| 亚洲精品一卡2卡三卡4卡5卡 | 一级片免费观看大全| www日本在线高清视频| 国产成人影院久久av| 久久久久久久国产电影| 亚洲美女黄色视频免费看| 国产在线一区二区三区精| 国产精品三级大全| www.999成人在线观看| 黄片播放在线免费| 叶爱在线成人免费视频播放| 一边摸一边做爽爽视频免费| 超色免费av| 午夜免费成人在线视频| 久久ye,这里只有精品| 国产亚洲欧美在线一区二区| 国产成人精品久久二区二区免费| 久久免费观看电影| 国产爽快片一区二区三区| 午夜av观看不卡| 亚洲国产毛片av蜜桃av| 制服人妻中文乱码| 国产男女超爽视频在线观看| 又黄又粗又硬又大视频| 亚洲美女黄色视频免费看| 午夜福利视频在线观看免费| 操出白浆在线播放| 美女中出高潮动态图| 久久精品久久精品一区二区三区| 黄频高清免费视频| 亚洲一区中文字幕在线| 一区在线观看完整版| 搡老岳熟女国产| 久久久久久久久免费视频了| videosex国产| 侵犯人妻中文字幕一二三四区| 久久这里只有精品19| 亚洲精品国产区一区二| 操美女的视频在线观看| 免费在线观看完整版高清| 国产在视频线精品| 免费观看人在逋| 免费看av在线观看网站| 91麻豆av在线| 19禁男女啪啪无遮挡网站| 少妇猛男粗大的猛烈进出视频| 国产女主播在线喷水免费视频网站| 欧美日韩亚洲高清精品| 大片电影免费在线观看免费| 亚洲中文av在线| 色播在线永久视频| 久久精品aⅴ一区二区三区四区| 国产日韩欧美在线精品| 国产国语露脸激情在线看| 久久精品亚洲熟妇少妇任你| 男女免费视频国产| av在线播放精品| 夫妻性生交免费视频一级片| 91精品国产国语对白视频| 亚洲精品一二三| 欧美日韩视频高清一区二区三区二| 国产亚洲精品久久久久5区| 一本大道久久a久久精品| 国产一区二区三区av在线| 男女午夜视频在线观看| av天堂久久9| 亚洲久久久国产精品| 国产精品久久久久成人av| 久9热在线精品视频| 国产精品熟女久久久久浪| 后天国语完整版免费观看| 亚洲熟女毛片儿| 午夜av观看不卡| 日韩电影二区| 黄色怎么调成土黄色| 美女脱内裤让男人舔精品视频| 啦啦啦啦在线视频资源| 精品久久蜜臀av无| 一区二区日韩欧美中文字幕| 午夜91福利影院| 免费人妻精品一区二区三区视频| 亚洲国产欧美日韩在线播放| 多毛熟女@视频| 国产av精品麻豆| 国产精品一二三区在线看| 别揉我奶头~嗯~啊~动态视频 | 乱人伦中国视频| 制服人妻中文乱码| 亚洲欧美一区二区三区黑人| 999久久久国产精品视频| 亚洲欧洲国产日韩| 亚洲,欧美精品.| 五月天丁香电影| 又大又黄又爽视频免费| 国产精品麻豆人妻色哟哟久久| av在线播放精品| 麻豆av在线久日| 少妇 在线观看| 91国产中文字幕| 精品视频人人做人人爽| 国产免费一区二区三区四区乱码| 欧美日韩亚洲综合一区二区三区_| 999精品在线视频| 国产精品 欧美亚洲| 国产男女超爽视频在线观看| 一本大道久久a久久精品| 在线观看免费午夜福利视频| 午夜福利一区二区在线看| 一本大道久久a久久精品| 亚洲欧美精品自产自拍| 嫩草影视91久久| 亚洲国产欧美在线一区| 成人午夜精彩视频在线观看| 国产精品二区激情视频| 精品一区二区三区四区五区乱码 | 亚洲情色 制服丝袜| 一二三四社区在线视频社区8| 久久午夜综合久久蜜桃| 一级毛片黄色毛片免费观看视频| 成年人午夜在线观看视频| 日本vs欧美在线观看视频| 成年av动漫网址| 亚洲av电影在线观看一区二区三区|