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

    提取標(biāo)記點(diǎn)中心在子孔徑拼接檢測中的應(yīng)用

    2014-02-02 14:36:33隋永新楊懷江
    中國光學(xué) 2014年5期
    關(guān)鍵詞:面形孔徑光學(xué)

    張 敏,隋永新,楊懷江

    (1.中國科學(xué)院長春光學(xué)精密機(jī)械與物理研究所應(yīng)用光學(xué)國家重點(diǎn)實(shí)驗(yàn)室,吉林長春130033;2.中國科學(xué)院大學(xué),北京100049)

    提取標(biāo)記點(diǎn)中心在子孔徑拼接檢測中的應(yīng)用

    張 敏1,2,隋永新1*,楊懷江1

    (1.中國科學(xué)院長春光學(xué)精密機(jī)械與物理研究所應(yīng)用光學(xué)國家重點(diǎn)實(shí)驗(yàn)室,吉林長春130033;2.中國科學(xué)院大學(xué),北京100049)

    子孔徑拼接干涉儀中子孔徑定位精度難以在大行程范圍內(nèi)得到保證,為此本文提出了基于提取標(biāo)記點(diǎn)中心定位子孔徑的拼接方法。以標(biāo)記點(diǎn)的中心坐標(biāo)為標(biāo)記點(diǎn)坐標(biāo),根據(jù)標(biāo)記點(diǎn)在兩子孔徑局部坐標(biāo)系下的坐標(biāo)計(jì)算兩子孔徑之間的坐標(biāo)變換,將所有子孔徑數(shù)據(jù)坐標(biāo)變換到統(tǒng)一坐標(biāo)系下,利用機(jī)械誤差補(bǔ)償算法拼接出全口徑面形。在搭建的拼接檢測系統(tǒng)上實(shí)現(xiàn)了外徑468mm的平面鏡拋光過程和最終的全口徑面形測量,加工過程中的測量結(jié)果為面形誤差修正提供了準(zhǔn)確的數(shù)據(jù),保證了最終全口徑面形誤差RMS快速收斂到35 nm。實(shí)驗(yàn)證明,基于提取標(biāo)記點(diǎn)中心的子孔徑拼接檢測能放寬對機(jī)械定位精度的要求,有效檢測大口徑光學(xué)元件面形。

    光學(xué)檢測;干涉儀;子孔徑拼接算法;標(biāo)記點(diǎn)中心

    1 引言

    子孔徑拼接干涉檢測技術(shù)是由美國Arizona光學(xué)中心的Kim與Wyant在1981年首次提出的[1],他們實(shí)現(xiàn)了用小口徑的平面反射鏡陣列來代替大口徑的平面反射鏡,從而完成拋物面鏡的自準(zhǔn)直檢驗(yàn),并提出可將子孔徑拼接測量用于非球面的檢測,在此之后國內(nèi)外很多科技工作者都對子孔徑測量技術(shù)進(jìn)行了研究[2]。在國內(nèi),國防科技大學(xué)在子孔徑迭代算法和高數(shù)值孔徑的子孔徑拼接檢測上碩果累累[3-4],西安交通大學(xué)、中國科學(xué)院光電技術(shù)研究所、長春光學(xué)精密機(jī)械與物理研究所等單位[5-8]的科研人員在子孔徑拼接檢測大口徑光學(xué)元件方向上均有不少研究。最近幾年一些國內(nèi)外的子孔徑拼接技術(shù)研究者提出了新的子孔徑定位方法,如圖像處理中的搜索方法[6,9,10]。其中Roland Maurer[9]所在的工作組采用標(biāo)記點(diǎn)輔助搜索的子孔徑定位方法,找出使兩個(gè)子孔徑相關(guān)系數(shù)最大的位置,將其作為子孔徑的定位,但該方法僅在x、y方向的平移定位中有效。

    子孔徑拼接檢測較大口徑的光學(xué)元件需要干涉儀或待測元件進(jìn)行大行程的移動(dòng),移動(dòng)距離越大定位精度就越低,其次還有一些不具備定位條件的檢測情況。在相機(jī)標(biāo)定、三維重建等三維形貌測量等領(lǐng)域,科研工作者廣泛采用人工標(biāo)記點(diǎn)作為重要的定位方法,其中易于識別、定位精度高的圓形標(biāo)記點(diǎn)最為常用[11]。所以本文采用圓形標(biāo)記點(diǎn)定位的方法,提出了基于自動(dòng)提取標(biāo)記點(diǎn)中心的拼接算法,解決測量大口徑平面鏡定位精度較低的問題,避免了為提高機(jī)構(gòu)運(yùn)動(dòng)準(zhǔn)確度而增加成本。本文通過外口徑為468 mm、內(nèi)口徑為202 mm的環(huán)形平面鏡鏡面拋光階段的面形檢測,驗(yàn)證了該方法的可行性以及實(shí)用性。

    2 基于標(biāo)記點(diǎn)定位的子孔徑拼接過程

    2.1 子孔徑拼接過程

    子孔徑拼接檢測過程如圖1,首先進(jìn)行子孔徑劃分,按干涉儀視場、待測鏡大小以及子孔徑兩兩之間重疊面積大小,確定子孔徑的個(gè)數(shù)以及分布。在重疊區(qū)域做標(biāo)記點(diǎn)之后,干涉檢測各子孔徑面形信息。對各子孔徑數(shù)據(jù)邊緣檢測,檢測出標(biāo)記點(diǎn)邊緣,根據(jù)標(biāo)記點(diǎn)邊緣擬合出相應(yīng)的橢圓方程,以橢圓中心坐標(biāo)為標(biāo)記點(diǎn)的中心坐標(biāo),計(jì)算出相鄰兩子孔徑之間的坐標(biāo)變換。最后,將各子孔徑的坐標(biāo)系轉(zhuǎn)換到全局坐標(biāo)系下,采用合適的子孔徑拼接算法拼接出全口徑的面形信息。

    2.2 提取標(biāo)記點(diǎn)的中心位置

    在完成各子孔徑的測量后,首先對各帶有標(biāo)記點(diǎn)子孔徑數(shù)據(jù)進(jìn)行Canny邊緣檢測和對邊緣檢測結(jié)果進(jìn)行亞像素定位。Canny邊緣檢測算法對圖像進(jìn)行平滑處理的是二維高斯函數(shù)[12]:

    式中,σ為高斯函數(shù)參數(shù),由它控制平滑處理的平滑程度[13]。當(dāng)σ較大時(shí)算法定位精度較低,但信噪比較高,σ較小時(shí)算法情況正好相反。所以選取高斯函數(shù)參數(shù)時(shí)得根據(jù)需要選擇最佳值。本文中σ取圖像灰度值的均方根值,并在雙閾值方法中將圖像的平均灰度值作為高閾值選取的參數(shù)[14]。Canny邊緣檢測后可以得到一幅邊緣的二值圖像,如圖2(a),但圖中邊緣有背景產(chǎn)生的雜亂邊緣。因此,在對圖2(a)中的邊緣數(shù)據(jù)施加邊緣像素?cái)?shù)、邊緣閉合、擬合誤差等約束[15],可以得到只包含標(biāo)記點(diǎn)邊緣的圖像,如圖2(b)所示。

    從圖2(b)中可看出邊緣仍有毛刺,所以仍需要將邊緣附近的像素分解以得到標(biāo)記點(diǎn)邊緣的精確坐標(biāo),來準(zhǔn)確地?cái)M合橢圓方程,這就是邊緣亞像素定位。由于標(biāo)記點(diǎn)邊緣是階躍邊緣,所以標(biāo)記點(diǎn)邊緣的準(zhǔn)確坐標(biāo)是使鄰近邊緣像素點(diǎn)擬合出的曲面的梯度方向二階導(dǎo)數(shù)等于零的點(diǎn)。為了同時(shí)兼顧準(zhǔn)確性以及效率,本文采用三次多項(xiàng)式擬合邊緣[16],對每一個(gè)邊緣點(diǎn)位置做亞像素修正,整體修正效果如圖2(c)所示。

    各子孔徑中的圓形標(biāo)記點(diǎn)的成像為平面橢圓,所以對各子孔徑數(shù)據(jù)邊緣檢測,亞像素邊緣定位后,可根據(jù)標(biāo)記點(diǎn)邊緣數(shù)據(jù)利用最小二乘擬合出相應(yīng)的橢圓方程,用橢圓方程的系數(shù)計(jì)算出相應(yīng)的橢圓中心坐標(biāo),將橢圓中心坐標(biāo)作為標(biāo)記點(diǎn)坐標(biāo)來計(jì)算兩子孔徑之間的坐標(biāo)變換。在二維平面坐標(biāo)系中,橢圓利用圓錐曲線方程的代數(shù)形式表示如下[15]:

    為提高計(jì)算的精度,采用亞像素邊緣檢測得到的所有離散點(diǎn)的坐標(biāo)數(shù)據(jù)組成超定方程組,通過最小二乘法擬合確定A、B、C、D、E、F的值。橢圓中心坐標(biāo)為:

    當(dāng)標(biāo)記點(diǎn)所在圖像質(zhì)量較好的時(shí)候,亞像素邊緣定位得到的標(biāo)記點(diǎn)中心坐標(biāo)精度高于0.01 pixel,在圖像質(zhì)量比較差時(shí),該方法仍有0.02 pixel的精度[16]。

    3 坐標(biāo)變換及全局拼接的數(shù)學(xué)模型

    3.1 坐標(biāo)變換

    若有兩個(gè)空間A和B,對于空間中點(diǎn)和矢量的剛體變換,主要包括平移(tx,ty)、旋轉(zhuǎn)(θ)。點(diǎn)(x,y)在空間B中,則其在空間A中的坐標(biāo)(u,υ)可表示為:

    若兩相鄰子孔徑,子孔徑a的3個(gè)標(biāo)記點(diǎn)坐標(biāo)為(x1,y1)、(x2,y2)、(x3,y3),子孔徑b中對應(yīng)的標(biāo)記點(diǎn)坐標(biāo)為(x′1,y′1)、(x′2,y′2)、(x′3,y′3),可以通過式(5)計(jì)算出子孔徑b相對于子孔徑a的相對旋轉(zhuǎn)角θ。通過3個(gè)標(biāo)記點(diǎn)算出3個(gè)相對選轉(zhuǎn)角度θ,最后取均值θavg,并可通過式(4)計(jì)算出tx、ty,對子孔徑b所有坐標(biāo)點(diǎn)(x,y)進(jìn)行式(4)表示的坐標(biāo)變換,得到在子孔徑a坐標(biāo)系下的坐標(biāo)值。

    3.2 全局拼接的子孔徑拼接方法

    干涉儀檢測子孔徑過程中,干涉儀與待測鏡之間的相對運(yùn)動(dòng)會(huì)引入機(jī)械運(yùn)動(dòng)誤差以及調(diào)整誤差,這將在子孔徑測量結(jié)果中引入平移(piston)和傾斜(tilt)等誤差,而且會(huì)使實(shí)際檢測的子孔徑所在位置與規(guī)劃的子孔徑位置有偏差。因此相鄰兩子孔徑重疊區(qū)相位并不一致。機(jī)械定位誤差補(bǔ)償算法在坐標(biāo)變換后子孔徑相位分布wi(x,y)(i =1,2,...,N-2)中引入調(diào)整誤差補(bǔ)償項(xiàng)、機(jī)械定位誤差補(bǔ)償項(xiàng),采用對所有子孔徑同時(shí)優(yōu)化的最小二乘法擬合出各子孔徑的補(bǔ)償系數(shù)來補(bǔ)償機(jī)械誤差。

    若N為子孔徑個(gè)數(shù),w0(x,y)為基準(zhǔn)子孔徑相位分布。引入調(diào)整誤差補(bǔ)償項(xiàng)、機(jī)械定位誤差

    補(bǔ)償項(xiàng)的相位分布可記為[17]:

    式中,P是z方向的平移系數(shù),Tx、Ty分別是x、y方向子孔徑的傾斜系數(shù),Δx、Δy分別是x、y方向子孔徑定位誤差系數(shù),各系數(shù)采用最小二乘法擬合,即:

    其中

    4 實(shí)驗(yàn)驗(yàn)證

    結(jié)合本實(shí)驗(yàn)室進(jìn)行的外直徑468 mm、內(nèi)直徑202 mm的環(huán)形平面鏡鏡面的拋光加工過程,在各次加工迭代過程中以及最終的面形檢測中,應(yīng)用已搭建的拼接裝置和上述算法對平面鏡進(jìn)行子孔徑拼接測量實(shí)驗(yàn)。如圖3所示,搭建的拼接檢測系統(tǒng)中檢測平臺可進(jìn)行二維(傾斜和俯仰)微調(diào),并在平臺上放置可360°連續(xù)旋轉(zhuǎn)的旋轉(zhuǎn)臺;測量各個(gè)子孔徑的干涉儀為Zygo干涉儀,參考面(TF)直徑為150 mm。

    子孔徑規(guī)劃如圖4,共20個(gè)子孔徑。平面鏡放置在旋轉(zhuǎn)臺上后進(jìn)行調(diào)整,使平面鏡中心在旋轉(zhuǎn)臺的旋轉(zhuǎn)軸上。測量第一個(gè)子孔徑后,利用MetroPro軟件中Fiducials功能標(biāo)記其干涉條紋圖中的標(biāo)記點(diǎn),并將它們作為其他子孔徑中標(biāo)記點(diǎn)的初始位置。檢測完所有子孔徑,圖5為其中的一子孔徑的干涉條紋圖。接下來確定每個(gè)子孔徑數(shù)據(jù)中標(biāo)記點(diǎn)的位置,在每個(gè)標(biāo)記點(diǎn)初始位置上下左右各20個(gè)像素的范圍內(nèi)進(jìn)行邊緣檢測,其中一個(gè)標(biāo)記點(diǎn)邊緣檢測結(jié)果如圖2所示。

    各子孔徑數(shù)據(jù)按檢測平臺自帶的定位的拼接結(jié)果如圖6(a)所示,拼接結(jié)果中有拼接痕跡,并且最后一個(gè)子孔徑與第一個(gè)子孔徑重疊處明顯錯(cuò)開。將子孔徑測量結(jié)果和邊緣檢測結(jié)果代入到拼接算法中,得到的環(huán)形平面鏡初始面形如圖6(b),PV值為2 663 nm,RMS值為214 nm,拼接結(jié)果均利用Data Fill功能填補(bǔ)標(biāo)記點(diǎn)造成的空白數(shù)據(jù)。盡管子孔徑的數(shù)目較多,圖6(b)仍顯示拼接得到的面形平滑連續(xù);圖6(b)與6(a)相比幾乎看不到拼接痕跡,說明按標(biāo)記點(diǎn)定位的拼接結(jié)果較為準(zhǔn)確,進(jìn)一步驗(yàn)證了在檢測平臺機(jī)械定位精度不高或自身機(jī)械運(yùn)動(dòng)定位較困難的情況下標(biāo)記點(diǎn)定位的實(shí)用性。

    將數(shù)據(jù)提交給數(shù)控加工系統(tǒng),以此對鏡面進(jìn)行誤差修正。第一輪加工和第五輪加工后面形檢測結(jié)果分別如圖7、8所示,經(jīng)過五輪加工后RMS值達(dá)35 nm,并且中、低頻誤差明顯減少。上述拼接結(jié)果均除去加工造成的邊緣效應(yīng)。將基于提取標(biāo)記點(diǎn)中心的拼接檢測數(shù)據(jù)作為每次加工的依據(jù),最終得到較好的加工結(jié)果,這間接表明了該方法的有效性,同時(shí)在實(shí)踐中驗(yàn)證了該方法的可行性。

    5 結(jié)論

    本文利用圓形標(biāo)記點(diǎn)中心位置來計(jì)算相鄰子孔徑之間坐標(biāo)變換的方法,采用全局優(yōu)化的機(jī)械誤差補(bǔ)償算法拼接統(tǒng)一坐標(biāo)系后的各子孔徑數(shù)據(jù),為環(huán)形平面鏡拋光過程提供了準(zhǔn)確的面形數(shù)據(jù),最終指導(dǎo)加工出了面形RMS值為35 nm的光學(xué)鏡面。實(shí)驗(yàn)驗(yàn)證了圓形標(biāo)記點(diǎn)定位和機(jī)械誤差補(bǔ)償算法相結(jié)合的拼接檢測方法的可靠性,為大口徑平面光學(xué)元件的檢測提供了技術(shù)基礎(chǔ),同時(shí)也是子孔徑拼接干涉檢測方法在生產(chǎn)實(shí)踐中應(yīng)用的經(jīng)驗(yàn),但是根據(jù)標(biāo)記點(diǎn)中心坐標(biāo)得到的坐標(biāo)變換矩陣的精度仍需進(jìn)一步分析。

    [1]KIM C J,WYANT JC.Subaperture test of a large flat on a fast spherical surface[J].Opt.Soc.Am.,1981,71:15-87.

    [2]侯溪,伍凡,楊力,等.子孔徑拼接干涉測試技術(shù)現(xiàn)狀及發(fā)展趨勢[J].光學(xué)與光電技術(shù),2005,3(3):50-53. HOU X,WU F,YANG L,et al..Status and development t rend of sub-aperture stitching interferometric testing technique[J].Opt.Optoelectronic Technology,2005,3(3):50-53.(in Chinese)

    [3]CHEN SH Y,DAIY F.Error reductions for stitching test of large optical flats[J].Opt.Laser Technology,2012,44(5):1543-1550.

    [4]CHEN SH Y,LIAOW L.Self-calibrated subaperture stitching testof hyper-hemispheres using latitude and longitude coordinates[J].Optical Society America,2012,51(17):3817-3825.

    [5]張鵬飛,趙宏,周翔,等.采用立體視覺實(shí)現(xiàn)子孔徑拼接測量的工件定位[J].光學(xué)精密工程,2010,18(2):503-511. ZHANG P F,ZHAO H,ZHOU X,et al..Work-piece localization in sub-aperture stitching test based on stereo vision[J]. Opt.Precision Eng.,2010,18(2):503-511.(in Chinese).

    [6]汪利華,吳時(shí)彬,任戈,等.子孔徑拼接檢測光學(xué)系統(tǒng)波前機(jī)械定位誤差補(bǔ)償算法[J].光學(xué)學(xué)報(bào),2012,32(1):113-118. WANG L H,WU SH B,REN G,et al..Location error compensation algorithm formeasure optical system wave front[J]. Acta Optica Sinica,2012,32(1):113-118.(in Chinese)

    [7]王孝坤,王麗輝,鄧偉杰,等.用非零位補(bǔ)償法檢測大口徑非球面反射鏡[J].光學(xué)精密工程,2011,19(3):520-528. WANG X K,WANG L H,DENGW J,etal.Measurement of large asphericmirrors by non-null testing[J].Opt.Precision Eng.,2011,19(3):520-528.(in Chinese).

    [8]陳一巍,王飛,王高文,等.基于變換的子孔徑拼接新算法[J].光學(xué)學(xué)報(bào),2013,33(9):0912004. CHEN YW,WANG F,WANGGW,etal..New sub-aperture stitching algorithm based on transformation[J].Acta Optica Sinica,2013,33(9):0912004.(in Chinese)

    [9]ROLAND M,FLORIAN S,CHRISTIAN V,et al..Physicalmarker based stitching process of circular and non-circular interferograms[J].SPIE,2011,8083:80830Q.

    [10]LIN PC,CHEN Y A,CHANG H SH,et al..Aberration compensation and position scanning of a subaperture stitching algorithm[J].SPIE,2012,8494:84940L.

    [11]許策,馬杰,趙全明,等.改進(jìn)的圓形標(biāo)志亞像素級中心檢測方法[J].計(jì)算機(jī)工程,2013,39(1):217-220. XU C,MA J,ZHAO Q M,et al..Improved method of sub-pixel centre detection in circlemark[J].Computer Engineering,2013,39(1):217-220.(in Chinese)

    [12]CANNY JF.A computational approach to edge detection[J].IEEE Trans,Pattern Analysis and Machine Intelligence, 1986,8(6):679-698.

    [13]YANG J,ZHANG Q.An improved Canny edge detection algorithm combined with steering kernel regression[J].SPIE, 2013,8768:87687L.

    [14]薛麗霞,李濤,王佐成.一種自適應(yīng)的Canny邊緣檢測算法[J].計(jì)算機(jī)應(yīng)用研究,2010,27(90):3588-3590. XUE L X,LI T,WANG Z CH.Adaptive Canny edge detection algorithm[J].Appl.Res.Computers,2010,27(90):3588-3590.(in Chinese)

    [15]殷永凱,劉曉利,李阿蒙,等.圓形標(biāo)志點(diǎn)的亞像素定位及其應(yīng)用[J].紅外與激光工程,2008,37(增刊):47-50. YIN Y K,LIU X L,LIA M,et al..Sub-pixel location of circle target and its application[J].Infrared and Laser Engineering,2008,37(Supp.):47-50.(in Chinese)

    [16]HARALICK R M.Digital step edges from zero crossing of second directional derivatives[J].IEEE Trans,Pattern Analysis and Machine Intelligence,1984,6(1):58-68.

    [17]DONALD G,GREG F,PAUL M.Method for selfcalibrated subaperture stitching for surface figure measurement:USA, 69566572B2[P].2005.

    張敏(1988-),女,河南杞縣人,博士研究生,2010年于吉林大學(xué)獲得學(xué)士學(xué)位,主要從事光學(xué)檢測方面的研究。E-mail:zhangminxiaowei@163.com

    楊懷江(1966-),男,遼寧丹東人,博士,研究員,博士生導(dǎo)師,主要從事紫外光刻技術(shù)、光學(xué)信息融合、網(wǎng)絡(luò)信息安全等方面的研究。E-mail:yanghj@ sklao.ac.cn

    隋永新(1970-),男,吉林長春人,博士,研究員,博士生導(dǎo)師,主要從事超高精度精密光學(xué)元件檢測、光學(xué)信息融合等方面的研究。E-mail:suiyx@sklao. ac.cn

    Subaperture stitching interferometry based on detection ofmarker center

    ZHANG Min1,2,SUIYong-xin1*,YANG Huai-jiang1
    (1.State Key Laboratory of Applied Optics,Changchun Institute of Optic, Fine Mechanics and Physics,Chinese Academy of Sciences,Changchun 130033,China;2.Uniυersity of Chinese Academy of Sciences,Beijing 100049,China)*Corresponding author,E-mail:suiyx@sklao.ac.cn

    The stitching interferometer systemswith larger relativemovementwill show a significantly lower positioning accuracy of subapertures.As a consequence a stitchingmethod based on detection of artificial circularmark center to find the necessary translation between two neighborly subapertures is implemented.Firstly, we take coordinates ofmark centers as themarks′coordinates by which the translation is computed.Then all the subaperture data are unified into the same reference by homogeneous coordinate transformation and the full aperture phase are stitched by usingmechanical system error compensation algorithm.A subaperture stitching process for a 468 mm flatmirror was carried out including surface accuracy tests during the polishing.In this process,subaperture stitching test offered the surface data precisely for polishing,which ensured the surfaceerror converged quickly to a final RMS of 35 nm.The experimental results show that themethod relaxes the precision requirement for subaperture location and can get the full aperture phase for large optical element correctly.

    opticalmeasurement;interferometer;subaperture stitching algorithm;marker center

    O436.1;TG84

    A

    10.3788/CO.20140705.0830

    2095-1531(2014)05-0830-07

    2014-04-15;

    2014-07-22

    國家科技重大專項(xiàng)基金資助項(xiàng)目(No.2009ZX02205)

    猜你喜歡
    面形孔徑光學(xué)
    反射鏡面形三坐標(biāo)白光掃描檢測精度研究
    滑輪組的裝配
    光學(xué)常見考題逐個(gè)擊破
    使用最小二乘迭代相移方法測量透明元件*
    不同滲透率巖芯孔徑分布與可動(dòng)流體研究
    超光譜儀光柵面形對光譜性能影響的仿真分析
    分布式孔徑相參合成雷達(dá)技術(shù)
    基于子孔徑斜率離散采樣的波前重構(gòu)
    采樣點(diǎn)分布對基于面形斜率徑向基模型的自由曲面擬合精度的影響
    光學(xué)遙感壓縮成像技術(shù)
    亚洲男人的天堂狠狠| av天堂在线播放| 久久天堂一区二区三区四区| 男女视频在线观看网站免费 | 最好的美女福利视频网| 欧美黄色片欧美黄色片| 久久性视频一级片| 老司机午夜福利在线观看视频| 欧美高清成人免费视频www| 亚洲中文字幕一区二区三区有码在线看 | 18禁观看日本| 91成年电影在线观看| 别揉我奶头~嗯~啊~动态视频| 久久久久精品国产欧美久久久| www日本在线高清视频| 亚洲美女黄片视频| 久久久久久国产a免费观看| 香蕉国产在线看| 日本在线视频免费播放| 国产免费男女视频| 亚洲乱码一区二区免费版| 欧美一区二区国产精品久久精品 | 黄色成人免费大全| 精华霜和精华液先用哪个| 哪里可以看免费的av片| 男人舔女人下体高潮全视频| 中文字幕最新亚洲高清| 亚洲欧美精品综合久久99| 国产单亲对白刺激| 天堂影院成人在线观看| 啦啦啦观看免费观看视频高清| 麻豆一二三区av精品| 不卡av一区二区三区| 亚洲国产高清在线一区二区三| 国产精品久久久久久亚洲av鲁大| 五月伊人婷婷丁香| 亚洲av美国av| 黄频高清免费视频| 一本久久中文字幕| 一进一出抽搐动态| 久久人妻av系列| 国产视频一区二区在线看| 国产av一区二区精品久久| 亚洲黑人精品在线| 91麻豆av在线| av福利片在线观看| 草草在线视频免费看| 69av精品久久久久久| 成熟少妇高潮喷水视频| 国产精品一区二区三区四区免费观看 | 丰满人妻熟妇乱又伦精品不卡| 又黄又爽又免费观看的视频| 丰满的人妻完整版| a在线观看视频网站| 欧美大码av| 亚洲熟女毛片儿| 久久中文看片网| 十八禁人妻一区二区| 国产成人啪精品午夜网站| 日韩精品中文字幕看吧| 久久精品夜夜夜夜夜久久蜜豆 | 搡老熟女国产l中国老女人| 久久欧美精品欧美久久欧美| 久久精品综合一区二区三区| 国产单亲对白刺激| 岛国在线观看网站| 国产一区二区在线观看日韩 | 久久这里只有精品中国| 亚洲国产精品成人综合色| 精品午夜福利视频在线观看一区| av免费在线观看网站| 99久久精品热视频| 18禁裸乳无遮挡免费网站照片| 高清在线国产一区| 午夜日韩欧美国产| www国产在线视频色| 这个男人来自地球电影免费观看| 久久天躁狠狠躁夜夜2o2o| 色播亚洲综合网| 两个人视频免费观看高清| 51午夜福利影视在线观看| 1024视频免费在线观看| 亚洲人成电影免费在线| 人妻久久中文字幕网| 男男h啪啪无遮挡| 在线a可以看的网站| 男人舔奶头视频| 他把我摸到了高潮在线观看| 极品教师在线免费播放| 毛片女人毛片| 露出奶头的视频| 欧美成人性av电影在线观看| 欧美成狂野欧美在线观看| 91大片在线观看| 在线观看www视频免费| 亚洲18禁久久av| 一级作爱视频免费观看| 看黄色毛片网站| 国产av在哪里看| 免费观看人在逋| 97超级碰碰碰精品色视频在线观看| 国产午夜精品论理片| www日本黄色视频网| 亚洲aⅴ乱码一区二区在线播放 | 亚洲,欧美精品.| 久久九九热精品免费| 人妻丰满熟妇av一区二区三区| 国产高清视频在线播放一区| 欧美性猛交黑人性爽| 免费在线观看亚洲国产| 国内精品一区二区在线观看| 五月玫瑰六月丁香| 午夜激情福利司机影院| 国产成+人综合+亚洲专区| 国产真人三级小视频在线观看| 久久亚洲真实| 看片在线看免费视频| 亚洲精品在线观看二区| 搡老妇女老女人老熟妇| 日本三级黄在线观看| 久久香蕉精品热| 亚洲av成人av| 黄色片一级片一级黄色片| 成人永久免费在线观看视频| 国产91精品成人一区二区三区| 麻豆成人午夜福利视频| 欧美日韩国产亚洲二区| 久久精品夜夜夜夜夜久久蜜豆 | 老熟妇乱子伦视频在线观看| 午夜a级毛片| 午夜福利视频1000在线观看| 国产午夜精品论理片| 一区福利在线观看| 久久久久九九精品影院| 日日爽夜夜爽网站| 少妇人妻一区二区三区视频| 给我免费播放毛片高清在线观看| 搡老熟女国产l中国老女人| 国产一区二区在线观看日韩 | 国产高清视频在线观看网站| 可以免费在线观看a视频的电影网站| 精品福利观看| 亚洲国产欧美人成| 狂野欧美激情性xxxx| 日日夜夜操网爽| 欧美日本视频| 日韩大码丰满熟妇| 午夜福利视频1000在线观看| 成人亚洲精品av一区二区| netflix在线观看网站| ponron亚洲| 99热只有精品国产| 夜夜爽天天搞| 久久精品国产99精品国产亚洲性色| 丁香六月欧美| 日本黄大片高清| 日韩欧美一区二区三区在线观看| а√天堂www在线а√下载| 国内精品久久久久久久电影| 日日爽夜夜爽网站| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品美女久久av网站| 午夜精品在线福利| 久久热在线av| 神马国产精品三级电影在线观看 | 久久久久免费精品人妻一区二区| 在线观看一区二区三区| 麻豆国产av国片精品| 搡老熟女国产l中国老女人| 婷婷六月久久综合丁香| tocl精华| 国产精品一区二区三区四区免费观看 | 国产亚洲av高清不卡| 亚洲精品粉嫩美女一区| 精品午夜福利视频在线观看一区| 亚洲 国产 在线| 国产成人精品久久二区二区91| 精品第一国产精品| 国产精品影院久久| 亚洲国产欧美一区二区综合| 日日摸夜夜添夜夜添小说| 一级a爱片免费观看的视频| 十八禁人妻一区二区| 国产成人精品久久二区二区免费| 国产精品免费一区二区三区在线| 国产麻豆成人av免费视频| 手机成人av网站| 十八禁网站免费在线| xxxwww97欧美| 成人特级黄色片久久久久久久| 中文字幕熟女人妻在线| 黄色丝袜av网址大全| www.自偷自拍.com| 免费看十八禁软件| 成人精品一区二区免费| 国产精品久久久人人做人人爽| 国产麻豆成人av免费视频| 亚洲av五月六月丁香网| 国产午夜精品久久久久久| 亚洲第一欧美日韩一区二区三区| 天天躁夜夜躁狠狠躁躁| 亚洲国产欧美一区二区综合| 特级一级黄色大片| 天堂√8在线中文| 婷婷丁香在线五月| 国内少妇人妻偷人精品xxx网站 | 欧美久久黑人一区二区| av中文乱码字幕在线| 国产成人影院久久av| 成人欧美大片| 免费av毛片视频| 亚洲欧美激情综合另类| 精品国产亚洲在线| 51午夜福利影视在线观看| 高潮久久久久久久久久久不卡| 一级a爱片免费观看的视频| 亚洲欧美日韩高清专用| 国产精品国产高清国产av| 九色成人免费人妻av| 久久久久久亚洲精品国产蜜桃av| 国内精品一区二区在线观看| 午夜福利成人在线免费观看| 99精品在免费线老司机午夜| 午夜视频精品福利| 在线免费观看的www视频| 免费在线观看视频国产中文字幕亚洲| 18禁观看日本| 叶爱在线成人免费视频播放| 欧美中文日本在线观看视频| 中文字幕高清在线视频| 免费在线观看影片大全网站| 久久国产乱子伦精品免费另类| 88av欧美| 母亲3免费完整高清在线观看| aaaaa片日本免费| 亚洲午夜理论影院| 国产av在哪里看| 免费观看人在逋| 99热这里只有精品一区 | 精品午夜福利视频在线观看一区| 岛国在线观看网站| 久久精品国产亚洲av香蕉五月| av超薄肉色丝袜交足视频| 国产人伦9x9x在线观看| 久久久久久免费高清国产稀缺| АⅤ资源中文在线天堂| 中文字幕精品亚洲无线码一区| 熟妇人妻久久中文字幕3abv| 国产av不卡久久| 亚洲国产中文字幕在线视频| 欧美av亚洲av综合av国产av| 久久久久性生活片| 老汉色∧v一级毛片| 99久久无色码亚洲精品果冻| 超碰成人久久| 丝袜美腿诱惑在线| 啦啦啦观看免费观看视频高清| 国产主播在线观看一区二区| 国产区一区二久久| 妹子高潮喷水视频| 亚洲国产欧美一区二区综合| 岛国在线观看网站| 久久中文字幕一级| 亚洲av中文字字幕乱码综合| 国产高清视频在线观看网站| 亚洲国产高清在线一区二区三| 日本免费a在线| 免费观看精品视频网站| 国产又色又爽无遮挡免费看| 欧美黄色淫秽网站| 身体一侧抽搐| 淫秽高清视频在线观看| 欧美人与性动交α欧美精品济南到| 亚洲国产精品sss在线观看| 午夜福利视频1000在线观看| 久久久久久亚洲精品国产蜜桃av| 中文亚洲av片在线观看爽| 天天添夜夜摸| 色综合欧美亚洲国产小说| 国产乱人伦免费视频| 草草在线视频免费看| 国产成人av教育| 久久人妻av系列| 亚洲美女黄片视频| 精品免费久久久久久久清纯| 亚洲av电影不卡..在线观看| 999久久久精品免费观看国产| 日韩欧美三级三区| 日韩大码丰满熟妇| 国产伦在线观看视频一区| 国产精品久久久人人做人人爽| 99riav亚洲国产免费| 精品久久久久久久久久久久久| 成人国产一区最新在线观看| 久久天躁狠狠躁夜夜2o2o| 黄色毛片三级朝国网站| 丁香欧美五月| a在线观看视频网站| 亚洲国产中文字幕在线视频| 国产一级毛片七仙女欲春2| 久久久久久大精品| 欧美一级毛片孕妇| 成人av一区二区三区在线看| 免费看美女性在线毛片视频| 国产成人欧美在线观看| www日本在线高清视频| 亚洲男人天堂网一区| 国产成人一区二区三区免费视频网站| 超碰成人久久| 精品国产美女av久久久久小说| 亚洲成av人片免费观看| 午夜日韩欧美国产| 日韩精品免费视频一区二区三区| 成年人黄色毛片网站| 久久精品国产综合久久久| 国产又色又爽无遮挡免费看| 亚洲成av人片免费观看| 9191精品国产免费久久| 国产精品av视频在线免费观看| av福利片在线| 欧美精品啪啪一区二区三区| 99国产极品粉嫩在线观看| 日本免费a在线| 国产亚洲精品av在线| av在线播放免费不卡| 久久性视频一级片| 精品一区二区三区视频在线观看免费| 两性夫妻黄色片| 国产精品免费视频内射| a级毛片在线看网站| 成人国产综合亚洲| 免费人成视频x8x8入口观看| 婷婷精品国产亚洲av| 成人午夜高清在线视频| 岛国在线观看网站| cao死你这个sao货| 激情在线观看视频在线高清| 欧美日韩乱码在线| 熟妇人妻久久中文字幕3abv| 亚洲欧美一区二区三区黑人| 亚洲午夜精品一区,二区,三区| 性色av乱码一区二区三区2| 亚洲成av人片在线播放无| 欧美日韩中文字幕国产精品一区二区三区| 中文字幕高清在线视频| or卡值多少钱| 国产成+人综合+亚洲专区| 国产精品av久久久久免费| 国产成人精品久久二区二区免费| or卡值多少钱| 一本大道久久a久久精品| 99国产综合亚洲精品| 国产亚洲精品av在线| 久久久精品欧美日韩精品| 中文字幕高清在线视频| 国产精品自产拍在线观看55亚洲| 又黄又粗又硬又大视频| 亚洲av第一区精品v没综合| 亚洲自偷自拍图片 自拍| 午夜亚洲福利在线播放| 又黄又爽又免费观看的视频| av超薄肉色丝袜交足视频| 午夜激情av网站| 啪啪无遮挡十八禁网站| 午夜激情av网站| 中亚洲国语对白在线视频| 黄色成人免费大全| 桃红色精品国产亚洲av| 老汉色∧v一级毛片| 最近在线观看免费完整版| 日本黄大片高清| 99在线视频只有这里精品首页| 免费在线观看影片大全网站| 午夜老司机福利片| 2021天堂中文幕一二区在线观| 免费在线观看亚洲国产| 日韩欧美精品v在线| 免费看日本二区| www日本黄色视频网| 亚洲av电影不卡..在线观看| 他把我摸到了高潮在线观看| 亚洲在线自拍视频| 19禁男女啪啪无遮挡网站| 免费一级毛片在线播放高清视频| 听说在线观看完整版免费高清| 在线国产一区二区在线| 在线观看美女被高潮喷水网站 | av中文乱码字幕在线| 97超级碰碰碰精品色视频在线观看| 国产精品影院久久| www国产在线视频色| 欧美在线黄色| 国产精品永久免费网站| 亚洲精品色激情综合| 精品国产乱子伦一区二区三区| 哪里可以看免费的av片| 在线国产一区二区在线| 国产黄色小视频在线观看| 欧美激情久久久久久爽电影| 欧美成人一区二区免费高清观看 | 男女床上黄色一级片免费看| 天堂av国产一区二区熟女人妻 | 看免费av毛片| 久久久国产精品麻豆| 久久精品国产综合久久久| 91av网站免费观看| 亚洲午夜精品一区,二区,三区| 亚洲国产欧美一区二区综合| 精品国产超薄肉色丝袜足j| 国产熟女午夜一区二区三区| 国产精品久久久久久亚洲av鲁大| 精品久久蜜臀av无| 丰满的人妻完整版| 一级毛片女人18水好多| 91成年电影在线观看| 久久精品亚洲精品国产色婷小说| 日韩av在线大香蕉| 一边摸一边做爽爽视频免费| 一级作爱视频免费观看| 国产亚洲精品第一综合不卡| 精品久久久久久成人av| 国产成年人精品一区二区| 亚洲精品av麻豆狂野| 嫩草影院精品99| 日韩国内少妇激情av| 久久亚洲真实| 亚洲一区二区三区色噜噜| 一进一出抽搐动态| 国产在线观看jvid| 午夜久久久久精精品| 天堂动漫精品| netflix在线观看网站| 一夜夜www| 成年版毛片免费区| 国产精华一区二区三区| 此物有八面人人有两片| 午夜福利成人在线免费观看| 美女午夜性视频免费| 成人亚洲精品av一区二区| 久久天堂一区二区三区四区| 丁香六月欧美| 亚洲欧美一区二区三区黑人| 亚洲人成网站在线播放欧美日韩| 在线观看日韩欧美| 亚洲乱码一区二区免费版| 美女扒开内裤让男人捅视频| 亚洲男人天堂网一区| 亚洲av中文字字幕乱码综合| 日本a在线网址| 又紧又爽又黄一区二区| 日韩有码中文字幕| 成人三级黄色视频| 日韩欧美国产在线观看| 欧美日韩中文字幕国产精品一区二区三区| 日本熟妇午夜| 久久性视频一级片| av中文乱码字幕在线| 欧美另类亚洲清纯唯美| 男女床上黄色一级片免费看| 久久精品国产综合久久久| 69av精品久久久久久| 亚洲国产高清在线一区二区三| 亚洲成人久久爱视频| 少妇人妻一区二区三区视频| 国产又黄又爽又无遮挡在线| 久久久国产精品麻豆| 99国产精品一区二区三区| 国模一区二区三区四区视频 | 亚洲av日韩精品久久久久久密| 久久久久免费精品人妻一区二区| 亚洲精华国产精华精| 人妻夜夜爽99麻豆av| av片东京热男人的天堂| 黄色女人牲交| 久久 成人 亚洲| 一个人免费在线观看电影 | 久久久久久久久久黄片| 国产亚洲精品综合一区在线观看 | 老汉色∧v一级毛片| 两人在一起打扑克的视频| 久久精品aⅴ一区二区三区四区| 久久久精品大字幕| 久久中文字幕一级| 亚洲人成网站在线播放欧美日韩| 99热6这里只有精品| 丁香欧美五月| 淫妇啪啪啪对白视频| 俄罗斯特黄特色一大片| 黑人巨大精品欧美一区二区mp4| 国产真实乱freesex| 成人三级黄色视频| 日本一本二区三区精品| 国产精品av久久久久免费| 一二三四在线观看免费中文在| 欧美日韩一级在线毛片| 国产不卡一卡二| 亚洲人成网站高清观看| 久久人妻福利社区极品人妻图片| 久久香蕉精品热| 精品一区二区三区四区五区乱码| 男女做爰动态图高潮gif福利片| 亚洲 欧美一区二区三区| 少妇裸体淫交视频免费看高清 | 欧美成人性av电影在线观看| 亚洲精品中文字幕在线视频| 久久久国产成人精品二区| 国产精品久久久久久精品电影| 看片在线看免费视频| cao死你这个sao货| 久久久久久九九精品二区国产 | 亚洲av成人av| 亚洲最大成人中文| 国产欧美日韩一区二区精品| 五月伊人婷婷丁香| 啦啦啦韩国在线观看视频| 亚洲av第一区精品v没综合| 精品一区二区三区视频在线观看免费| 好男人在线观看高清免费视频| 悠悠久久av| 在线看三级毛片| 色综合亚洲欧美另类图片| 国产成人av激情在线播放| 亚洲精品美女久久av网站| 免费一级毛片在线播放高清视频| 亚洲国产精品成人综合色| 一进一出抽搐gif免费好疼| 成熟少妇高潮喷水视频| 床上黄色一级片| 国产精品国产高清国产av| 国产真实乱freesex| 午夜福利欧美成人| 国产精品野战在线观看| 国产精品 国内视频| 嫩草影视91久久| 高清在线国产一区| 听说在线观看完整版免费高清| 国产视频一区二区在线看| 欧美色欧美亚洲另类二区| 18禁国产床啪视频网站| netflix在线观看网站| 丝袜人妻中文字幕| 亚洲真实伦在线观看| 正在播放国产对白刺激| 黄色视频,在线免费观看| 在线a可以看的网站| 99久久精品国产亚洲精品| 18禁国产床啪视频网站| 特级一级黄色大片| 久9热在线精品视频| 99re在线观看精品视频| 天堂av国产一区二区熟女人妻 | 757午夜福利合集在线观看| 后天国语完整版免费观看| 国产一区在线观看成人免费| 国产av麻豆久久久久久久| 一级片免费观看大全| 亚洲精品在线观看二区| 久久热在线av| 狠狠狠狠99中文字幕| 国产精品综合久久久久久久免费| 亚洲国产精品久久男人天堂| 91成年电影在线观看| 国产亚洲精品久久久久久毛片| 欧美3d第一页| 国产高清视频在线观看网站| 婷婷亚洲欧美| 免费在线观看亚洲国产| 岛国视频午夜一区免费看| 国产精品乱码一区二三区的特点| 一级作爱视频免费观看| 亚洲美女视频黄频| 国产熟女xx| 成人手机av| 亚洲国产看品久久| 国产激情偷乱视频一区二区| 久久久久国内视频| 国产成人系列免费观看| 丁香六月欧美| a在线观看视频网站| 国产在线观看jvid| 国产1区2区3区精品| 久久久久国内视频| 美女大奶头视频| 天堂av国产一区二区熟女人妻 | 老熟妇乱子伦视频在线观看| 久久这里只有精品19| 亚洲激情在线av| 老司机深夜福利视频在线观看| 后天国语完整版免费观看| 久久亚洲精品不卡| 香蕉久久夜色| 国产真实乱freesex| 国产精品久久久久久精品电影| 国产成+人综合+亚洲专区| 在线观看日韩欧美| 久久天堂一区二区三区四区| 99国产精品一区二区三区| 女警被强在线播放| 丰满人妻一区二区三区视频av | 欧美国产日韩亚洲一区| 91av网站免费观看| 欧美一区二区精品小视频在线| 久久久久精品国产欧美久久久| 91av网站免费观看| 日本免费一区二区三区高清不卡| 午夜久久久久精精品| 最近最新中文字幕大全免费视频| 此物有八面人人有两片| 黄色a级毛片大全视频| 国产精品久久久久久人妻精品电影| www.自偷自拍.com| 欧美日韩亚洲综合一区二区三区_| 麻豆久久精品国产亚洲av| 亚洲第一电影网av| 日本免费a在线|