• <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ù)
    久久精品91蜜桃| 欧美日韩中文字幕国产精品一区二区三区 | 婷婷丁香在线五月| 国产午夜精品久久久久久| 人人妻人人爽人人添夜夜欢视频| 亚洲精品久久成人aⅴ小说| 精品电影一区二区在线| 久久精品91蜜桃| 激情视频va一区二区三区| 9热在线视频观看99| av电影中文网址| 亚洲精品久久午夜乱码| 级片在线观看| 在线播放国产精品三级| 国产精品永久免费网站| 成人国语在线视频| 欧美中文日本在线观看视频| 国产成人欧美在线观看| 很黄的视频免费| 精品无人区乱码1区二区| 很黄的视频免费| 女人被狂操c到高潮| 十八禁人妻一区二区| 亚洲狠狠婷婷综合久久图片| 亚洲国产精品999在线| 精品久久久久久,| 亚洲精品一区av在线观看| 欧美+亚洲+日韩+国产| 国产精品亚洲一级av第二区| 水蜜桃什么品种好| 9热在线视频观看99| 不卡av一区二区三区| av电影中文网址| 国产精品野战在线观看 | av中文乱码字幕在线| 精品人妻1区二区| 免费在线观看影片大全网站| 999久久久精品免费观看国产| 国产无遮挡羞羞视频在线观看| 在线观看免费视频日本深夜| 欧美精品亚洲一区二区| 激情在线观看视频在线高清| 777久久人妻少妇嫩草av网站| 日本一区二区免费在线视频| 精品国产一区二区久久| 18禁美女被吸乳视频| 1024视频免费在线观看| 天堂中文最新版在线下载| 国产精品久久久久成人av| 日韩大尺度精品在线看网址 | 97超级碰碰碰精品色视频在线观看| 久久久国产成人免费| 一区在线观看完整版| 99热国产这里只有精品6| 桃红色精品国产亚洲av| 一级,二级,三级黄色视频| 欧美激情高清一区二区三区| 亚洲性夜色夜夜综合| 丝袜美腿诱惑在线| 日韩视频一区二区在线观看| av有码第一页| 长腿黑丝高跟| 亚洲中文av在线| 国产精品日韩av在线免费观看 | 精品人妻1区二区| 自线自在国产av| 悠悠久久av| 看片在线看免费视频| 亚洲欧美日韩高清在线视频| 久9热在线精品视频| a在线观看视频网站| 国产黄色免费在线视频| 免费看十八禁软件| 麻豆av在线久日| 超色免费av| 成人手机av| 亚洲人成电影免费在线| 一级片'在线观看视频| 成人18禁在线播放| 久久久久精品国产欧美久久久| 国产精品 欧美亚洲| 亚洲精品国产色婷婷电影| xxx96com| 亚洲av成人一区二区三| 国产成人精品久久二区二区免费| 香蕉久久夜色| 亚洲av成人一区二区三| 国产免费男女视频| 女性被躁到高潮视频| 一二三四在线观看免费中文在| 老司机在亚洲福利影院| 老司机在亚洲福利影院| 午夜视频精品福利| 精品久久久久久,| 国产97色在线日韩免费| 久久人妻熟女aⅴ| 亚洲欧美一区二区三区黑人| 久久久久久久久久久久大奶| 亚洲人成电影观看| 中国美女看黄片| 久久久国产成人免费| 女警被强在线播放| 日本撒尿小便嘘嘘汇集6| 熟女少妇亚洲综合色aaa.| 国产精品永久免费网站| 日韩欧美在线二视频| 看免费av毛片| 最近最新中文字幕大全电影3 | 国产精品98久久久久久宅男小说| 9热在线视频观看99| 国产成人欧美| 99国产极品粉嫩在线观看| 男女下面进入的视频免费午夜 | 久久午夜亚洲精品久久| 别揉我奶头~嗯~啊~动态视频| 美女 人体艺术 gogo| 色综合欧美亚洲国产小说| 一级片'在线观看视频| 久久人妻av系列| xxxhd国产人妻xxx| 在线观看舔阴道视频| 又黄又爽又免费观看的视频| av免费在线观看网站| 亚洲成人免费电影在线观看| 亚洲全国av大片| 亚洲成人精品中文字幕电影 | 免费高清在线观看日韩| 中文字幕人妻丝袜制服| 亚洲av熟女| 男人的好看免费观看在线视频 | 国产av又大| 两性夫妻黄色片| 午夜老司机福利片| 中出人妻视频一区二区| 丰满人妻熟妇乱又伦精品不卡| 久久人妻av系列| 在线观看免费高清a一片| 香蕉久久夜色| 午夜激情av网站| 丝袜在线中文字幕| 亚洲精品一区av在线观看| 午夜a级毛片| 一进一出好大好爽视频| 深夜精品福利| 国产精品爽爽va在线观看网站 | 女人精品久久久久毛片| 中文字幕色久视频| 欧美大码av| 日日爽夜夜爽网站| 国产精品久久久久成人av| 亚洲成人免费电影在线观看| 黑人巨大精品欧美一区二区蜜桃| 啦啦啦免费观看视频1| 国产精品九九99| 亚洲黑人精品在线| 亚洲色图 男人天堂 中文字幕| 在线av久久热| 久久久久精品国产欧美久久久| 国产激情欧美一区二区| 又黄又爽又免费观看的视频| 老司机靠b影院| 精品一区二区三卡| 十八禁网站免费在线| 叶爱在线成人免费视频播放| 免费搜索国产男女视频| 日韩有码中文字幕| 高清av免费在线| 久久人人爽av亚洲精品天堂| 欧美乱码精品一区二区三区| 一边摸一边做爽爽视频免费| 亚洲国产欧美日韩在线播放| 婷婷精品国产亚洲av在线| 午夜福利一区二区在线看| 黄色视频不卡| 男女之事视频高清在线观看| 黄片大片在线免费观看| 亚洲成a人片在线一区二区| 首页视频小说图片口味搜索| 丰满人妻熟妇乱又伦精品不卡| 亚洲午夜理论影院| 天堂俺去俺来也www色官网| 婷婷精品国产亚洲av在线| 国产免费av片在线观看野外av| 亚洲精品国产精品久久久不卡| 免费搜索国产男女视频| 一进一出抽搐动态| 久久人人爽av亚洲精品天堂| 国产精品乱码一区二三区的特点 | 国产视频一区二区在线看| 国产精品影院久久| 午夜福利一区二区在线看| 国产主播在线观看一区二区| 在线观看日韩欧美| 狂野欧美激情性xxxx| 国产精品一区二区免费欧美| 美国免费a级毛片| 亚洲欧美激情在线| 天天影视国产精品| 99在线人妻在线中文字幕| 69av精品久久久久久| 后天国语完整版免费观看| av网站在线播放免费| 久久国产亚洲av麻豆专区| 超碰成人久久| 黄色 视频免费看| 国产av一区二区精品久久| 日日干狠狠操夜夜爽| 99热国产这里只有精品6| 正在播放国产对白刺激| 99国产精品一区二区三区| 亚洲全国av大片| 国产精品免费一区二区三区在线| 在线看a的网站| 亚洲欧美激情综合另类| 国产日韩一区二区三区精品不卡| 亚洲专区字幕在线| 老司机靠b影院| 免费少妇av软件| 99精国产麻豆久久婷婷| 午夜精品久久久久久毛片777| 熟女少妇亚洲综合色aaa.| 久久国产亚洲av麻豆专区| 超碰97精品在线观看| 丰满的人妻完整版| 日韩欧美免费精品| 久99久视频精品免费| 国产高清videossex| 色老头精品视频在线观看| 神马国产精品三级电影在线观看 | 中文字幕人妻丝袜一区二区| 最近最新免费中文字幕在线| 亚洲一区高清亚洲精品| 亚洲 欧美一区二区三区| 午夜激情av网站| 在线视频色国产色| 亚洲av日韩精品久久久久久密| 男人的好看免费观看在线视频 | 高清在线国产一区| 老司机亚洲免费影院| 亚洲成人免费av在线播放| 国产精品日韩av在线免费观看 | 亚洲一码二码三码区别大吗| bbb黄色大片| 我的亚洲天堂| 午夜激情av网站| 在线观看免费视频网站a站| 国产有黄有色有爽视频| 桃色一区二区三区在线观看| 日韩大尺度精品在线看网址 | 中文字幕精品免费在线观看视频| 亚洲成人免费电影在线观看| 亚洲国产看品久久| 亚洲av成人一区二区三| 女人爽到高潮嗷嗷叫在线视频| 天堂俺去俺来也www色官网| 国产又色又爽无遮挡免费看| 90打野战视频偷拍视频| 男人操女人黄网站| 天堂影院成人在线观看| 老司机在亚洲福利影院| 国产精品久久久av美女十八| 色综合站精品国产| 一级作爱视频免费观看| 两个人看的免费小视频| 久久精品亚洲熟妇少妇任你| av中文乱码字幕在线| 亚洲中文av在线| 别揉我奶头~嗯~啊~动态视频| 黄片播放在线免费| 国产精品1区2区在线观看.| 免费av毛片视频| 麻豆av在线久日| 亚洲男人的天堂狠狠| 亚洲午夜精品一区,二区,三区| 亚洲成人免费av在线播放| 精品国产乱码久久久久久男人| 夜夜夜夜夜久久久久| 国产成人精品久久二区二区免费| 女人被躁到高潮嗷嗷叫费观| 黄色视频不卡| 欧美日韩福利视频一区二区| 精品国产乱子伦一区二区三区| 一级a爱视频在线免费观看| 欧美在线黄色| 精品福利永久在线观看| 成人三级做爰电影| 高清欧美精品videossex| 国产在线观看jvid| 交换朋友夫妻互换小说| 夫妻午夜视频| 欧美 亚洲 国产 日韩一| 777久久人妻少妇嫩草av网站| 亚洲片人在线观看| 久久久久久久精品吃奶| 国产欧美日韩综合在线一区二区| 国产成人欧美在线观看| 国产亚洲精品综合一区在线观看 | 成人特级黄色片久久久久久久| av视频免费观看在线观看| 啪啪无遮挡十八禁网站| 亚洲欧美精品综合一区二区三区| 亚洲美女黄片视频| 欧美激情高清一区二区三区| 黄色片一级片一级黄色片| 看黄色毛片网站| 亚洲欧美日韩无卡精品| 少妇裸体淫交视频免费看高清 | a级毛片在线看网站| av网站在线播放免费| 欧美人与性动交α欧美软件| 18禁观看日本| 久久久水蜜桃国产精品网| 日本一区二区免费在线视频| 美女福利国产在线| 伊人久久大香线蕉亚洲五| 午夜福利免费观看在线| 人人妻人人澡人人看| 一级毛片高清免费大全| 国产成人欧美| 亚洲av日韩精品久久久久久密| 如日韩欧美国产精品一区二区三区| 欧美大码av| 麻豆久久精品国产亚洲av | svipshipincom国产片| 又紧又爽又黄一区二区| 久久天堂一区二区三区四区| 桃红色精品国产亚洲av| 黄色怎么调成土黄色| 色哟哟哟哟哟哟| 亚洲成av片中文字幕在线观看| 麻豆国产av国片精品| 精品国内亚洲2022精品成人| 日本免费一区二区三区高清不卡 | 国产一区二区在线av高清观看| 18禁黄网站禁片午夜丰满| 一区二区三区国产精品乱码| 伊人久久大香线蕉亚洲五| av在线天堂中文字幕 | 午夜日韩欧美国产| 91成人精品电影| 久久人妻熟女aⅴ| 亚洲avbb在线观看| 9热在线视频观看99| 国产xxxxx性猛交| av超薄肉色丝袜交足视频| 丰满迷人的少妇在线观看| 不卡一级毛片| 国产一区二区激情短视频| 亚洲人成电影免费在线| netflix在线观看网站| 国产xxxxx性猛交| 这个男人来自地球电影免费观看| 9色porny在线观看| 99久久人妻综合| 久久人妻福利社区极品人妻图片| 免费久久久久久久精品成人欧美视频| 激情视频va一区二区三区| 国产精品九九99| 欧美黄色淫秽网站| 国产精品99久久99久久久不卡| 精品国产美女av久久久久小说| 制服诱惑二区| 在线观看免费高清a一片| 色在线成人网| 黄频高清免费视频| 成人亚洲精品av一区二区 | 欧美日韩乱码在线| 成人亚洲精品av一区二区 | 男女之事视频高清在线观看| √禁漫天堂资源中文www| 妹子高潮喷水视频| av有码第一页| 人妻久久中文字幕网| 动漫黄色视频在线观看| 88av欧美| 国产亚洲欧美98| 亚洲一区中文字幕在线| 亚洲视频免费观看视频| 99精品久久久久人妻精品| 搡老熟女国产l中国老女人| 国产又爽黄色视频| 亚洲精华国产精华精| 高潮久久久久久久久久久不卡| 欧美+亚洲+日韩+国产| 精品免费久久久久久久清纯| 制服人妻中文乱码| 久久久久久人人人人人| 一夜夜www| 亚洲欧洲精品一区二区精品久久久| 午夜视频精品福利| 国产精品自产拍在线观看55亚洲| √禁漫天堂资源中文www| 热re99久久精品国产66热6| 视频区欧美日本亚洲| 中出人妻视频一区二区| 国产视频一区二区在线看| 国产91精品成人一区二区三区| 三级毛片av免费| 欧美不卡视频在线免费观看 | 国产一卡二卡三卡精品| 久久午夜亚洲精品久久| 国产成+人综合+亚洲专区| 男人舔女人下体高潮全视频| 国产色视频综合| 亚洲专区字幕在线| 国产精品免费一区二区三区在线| 99久久国产精品久久久| 亚洲av五月六月丁香网| 黄色视频不卡| 精品久久蜜臀av无| 日本黄色视频三级网站网址| 天天添夜夜摸| 成人18禁高潮啪啪吃奶动态图| 大型黄色视频在线免费观看| 国产成人精品无人区| 看免费av毛片| 国产精品香港三级国产av潘金莲| 午夜精品久久久久久毛片777| 99国产精品一区二区三区| 国产成年人精品一区二区 | 亚洲精品中文字幕在线视频| 一夜夜www| 久久青草综合色| 国产主播在线观看一区二区| 欧美日韩福利视频一区二区| 777久久人妻少妇嫩草av网站| 女性生殖器流出的白浆| 伦理电影免费视频| 黄色 视频免费看| 日本a在线网址| 色尼玛亚洲综合影院| 精品福利观看| 天堂动漫精品| 亚洲精品国产一区二区精华液| 色哟哟哟哟哟哟| 亚洲七黄色美女视频| 午夜免费鲁丝| 精品一品国产午夜福利视频| 久久久国产成人免费| 手机成人av网站| 在线天堂中文资源库| 精品无人区乱码1区二区| 侵犯人妻中文字幕一二三四区| 女生性感内裤真人,穿戴方法视频| 国产激情久久老熟女| 纯流量卡能插随身wifi吗| 免费女性裸体啪啪无遮挡网站| 亚洲精品久久午夜乱码| 天天躁狠狠躁夜夜躁狠狠躁| 男人的好看免费观看在线视频 | 女性被躁到高潮视频| www.999成人在线观看| 午夜成年电影在线免费观看| 亚洲专区中文字幕在线| 亚洲国产欧美一区二区综合| netflix在线观看网站| 日韩欧美三级三区| 亚洲人成伊人成综合网2020| 一级a爱视频在线免费观看| 国产精品av久久久久免费| 9色porny在线观看| 久久国产精品男人的天堂亚洲| 中文字幕av电影在线播放| 亚洲精品久久午夜乱码| 国产不卡一卡二| 视频区图区小说| 亚洲三区欧美一区| 俄罗斯特黄特色一大片| 久久人人爽av亚洲精品天堂| 天天影视国产精品| 成人手机av| 亚洲 欧美一区二区三区| 最好的美女福利视频网| 国产欧美日韩一区二区精品| 又黄又爽又免费观看的视频| 88av欧美| 色尼玛亚洲综合影院| 亚洲成人精品中文字幕电影 | 人妻丰满熟妇av一区二区三区| 美女午夜性视频免费| 精品一区二区三区四区五区乱码| 热re99久久精品国产66热6| 母亲3免费完整高清在线观看| 国产成人精品久久二区二区免费| 男女床上黄色一级片免费看| 中文字幕最新亚洲高清| 搡老岳熟女国产| 国产麻豆69| 老熟妇仑乱视频hdxx| 嫁个100分男人电影在线观看| 午夜激情av网站| 日本免费a在线| 麻豆一二三区av精品| 一二三四社区在线视频社区8| 级片在线观看| 国产xxxxx性猛交| 极品人妻少妇av视频| 欧美日韩黄片免| 午夜免费激情av| 黄色成人免费大全| 中文字幕人妻丝袜制服| 我的亚洲天堂| 午夜免费观看网址| 老司机福利观看| av天堂在线播放| av福利片在线| 国产无遮挡羞羞视频在线观看| 久久精品91无色码中文字幕| 欧美激情久久久久久爽电影 | 国产一区二区激情短视频| 国产成人av激情在线播放| 久久精品国产亚洲av香蕉五月| 婷婷丁香在线五月| 亚洲美女黄片视频| 三级毛片av免费| 日本免费一区二区三区高清不卡 | 欧美精品一区二区免费开放| netflix在线观看网站| 国产免费av片在线观看野外av| 如日韩欧美国产精品一区二区三区| 亚洲 欧美一区二区三区| 97碰自拍视频| 极品教师在线免费播放| e午夜精品久久久久久久| 欧美 亚洲 国产 日韩一| 99久久人妻综合| 久久久久久久午夜电影 | 999精品在线视频| 亚洲国产欧美一区二区综合| 中文字幕精品免费在线观看视频| 久久久精品欧美日韩精品| 丝袜美足系列| 国产免费现黄频在线看| 久久青草综合色| 国产一区二区三区综合在线观看| 精品第一国产精品| 欧美激情 高清一区二区三区| 高清在线国产一区| 国产免费av片在线观看野外av| 美女扒开内裤让男人捅视频| 久久人人精品亚洲av| 国产麻豆69| 日本五十路高清| 色婷婷av一区二区三区视频| 亚洲专区中文字幕在线| 热99re8久久精品国产| 自线自在国产av| 18禁黄网站禁片午夜丰满| 欧美一区二区精品小视频在线| 欧美色视频一区免费| 天堂俺去俺来也www色官网| 亚洲五月婷婷丁香| 成人特级黄色片久久久久久久| 视频在线观看一区二区三区| 侵犯人妻中文字幕一二三四区| 一区二区三区精品91| 精品国产超薄肉色丝袜足j| 久久性视频一级片| 夜夜看夜夜爽夜夜摸 | av网站免费在线观看视频| 亚洲午夜精品一区,二区,三区| 51午夜福利影视在线观看| 午夜91福利影院| 免费高清视频大片| 91麻豆av在线| 一本综合久久免费| www.999成人在线观看| 人妻丰满熟妇av一区二区三区| 国产欧美日韩精品亚洲av| 国产精品电影一区二区三区| 国产一区二区三区综合在线观看| 可以在线观看毛片的网站| 视频区欧美日本亚洲| 在线国产一区二区在线| 黑丝袜美女国产一区| 久久精品国产综合久久久| 欧美日韩乱码在线| 日韩精品中文字幕看吧| 亚洲精华国产精华精| 夜夜看夜夜爽夜夜摸 | 精品日产1卡2卡| 亚洲国产看品久久| 国产亚洲欧美精品永久| 男人舔女人下体高潮全视频| √禁漫天堂资源中文www| 久久亚洲精品不卡| 国产精品久久电影中文字幕| 80岁老熟妇乱子伦牲交| 久久人妻av系列| 天堂√8在线中文| 黄色毛片三级朝国网站| 色综合站精品国产| 欧美乱码精品一区二区三区| 欧美激情高清一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 久久99一区二区三区| 国产精品一区二区精品视频观看| 麻豆av在线久日| 999久久久精品免费观看国产| 波多野结衣一区麻豆| 伊人久久大香线蕉亚洲五| 别揉我奶头~嗯~啊~动态视频| 最好的美女福利视频网| 欧美成人免费av一区二区三区| 婷婷六月久久综合丁香| 中文字幕最新亚洲高清| 男女床上黄色一级片免费看| 日日夜夜操网爽| 日韩视频一区二区在线观看| av电影中文网址| 亚洲国产精品一区二区三区在线| 老司机在亚洲福利影院| 18禁裸乳无遮挡免费网站照片 | 女性生殖器流出的白浆| 午夜91福利影院|