亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        外視場拼接測量系統(tǒng)的視場拼接和交匯測量算法及其實現

        2010-05-11 08:26:15宋立維喬彥峰
        中國光學 2010年3期
        關鍵詞:交匯測量儀視場

        王 旻,宋立維,喬彥峰,余 毅

        (1.中國科學院長春光學精密機械與物理研究所,吉林 長春 130033;2.中國科學院研究生院 北京 100039)

        外視場拼接測量系統(tǒng)的視場拼接和交匯測量算法及其實現

        王 旻1,2,宋立維1,喬彥峰1,余 毅1

        (1.中國科學院長春光學精密機械與物理研究所,吉林 長春 130033;2.中國科學院研究生院 北京 100039)

        基于2臺外視場拼接高速電視測量儀原型樣機(硬件)提出了一種將交匯測量和拼接處理相結合的算法。根據需要對地球質心和光電測量系統(tǒng)建立了5個坐標系,介紹了坐標系的定義及其它們之間的變換過程,給出了目標軌跡交匯測量和視場拼接的實例。用2臺高速電視測量儀拼接的8臺測量相機同時對同一運動目標進行測量,對獲得的測量數據進行交匯和視場拼接處理,結果顯示提出的拼接算法是正確、有效的,可以得到唯一解;地球曲率半徑和地球子午線收斂兩項影響因素得到了完全修正,表明交匯測量算法完全可以推廣到其它光電測控儀器的交匯測量。

        外視場拼接;交匯測量;地球曲率半徑;子午線收斂;誤差修正

        1 引言

        在目標高速測量領域,被測目標速度快、目標多、所占的視場范圍大,應用單臺的高速攝像機不能拍攝到完整的畫面。為了克服單臺相機拍攝的不足,研制了一種將4臺高速相機拼接為1臺高速電視測量儀的方法,實現了對多個目標的完整拍攝。本文基于2臺外視場拼接高速電視測量儀原型樣機(硬件),提出了一種將交匯測量和拼接處理相結合的算法。用2臺高速電視測量儀拼接的8臺測量相機同時對同一運動目標進行測量,對獲得的測量數據進行交匯和視場拼接處理。結果顯示,提出的拼接算法是正確、有效的,可以得到唯一解;地球曲率半徑和地球子午線收斂兩項影響因素得到了完全修正,表明提出的交匯測量算法,完全可以推廣到其它光電測控儀器的交匯測量。

        2 外場拼接原理

        外視場拼接測量系統(tǒng),一般由3個測量站光電測量儀組成,首先通過2臺以上的光電測量儀測得的數據進行交匯計算,獲得每個目標在空間的三維位置信息,然后進行視場拼接[1,2]。采用外視場拼接,光學系統(tǒng)簡單、重量輕、體積小、制造難度低,便于使用,但由于外視場拼接不可避免存在視差,因此進行視場拼接時,計算量大。近代計算機技術的發(fā)展,給外視場拼接提供了可能。如圖1所示的配置,每個測量站光電測量儀上安裝4個測量相機,這4個測量相機上下、左右對稱地配置在“T”字型跟蹤架水平軸外伸兩軸端,組成“田”字型布局配置,即為較典型的外視場拼接測量系統(tǒng)。4臺攝像機視場在外視場呈水平方向“一”字拼接,且采用凝視工作模式。

        圖1 外視場拼接測量系統(tǒng)的外觀圖Fig.1 Outlook of external field butting measuring system

        由于各光電測量儀的4個測量相機在上下、左右方向都存在視差,在目標沒有距離信息的情況下,只依靠單臺光電測量儀的信息不可能完成精密視場拼接,只有應用兩臺光電測量儀所獲得的信息,對同一目標同時進行視場拼接和交匯測量,才能完成精密的視場拼接。

        考慮采用“T”字型跟蹤架[3],“田”字型布局配置的光電測量儀,4個測量相機在上下、左右方向都存在視差,其視場拼接和交匯測量的數學模型具有典型性,本文對這一方案進行了分析,并以一個目標為例,建立了視場拼接和交匯測量的數學模型,在建立數學模型時,考慮了地球曲率半徑和子午線收斂誤差的影響[4,5]。

        3 坐標系定義

        根據實際分析的需要共建立5個坐標系[6~8](見圖2)。各坐標系(右手坐標系)定義如下:

        (1)地球質心坐標系I(i1,i2,i3);

        (2)第i個光電測量儀地平坐標系Si(si1,si2,si3);

        (3)第i個光電測量儀垂直軸坐標系Vi(vi1,vi2,vi3);

        (4)第i個光電測量儀水平軸坐標系Hi(hi1,hi2,hi3);

        (5)第i個光電測量儀、第j個測量相機坐標系Cij(cij1,cij2,cij3)。

        圖2 外視場拼接測量系統(tǒng)在地心質心坐標系中的坐標系統(tǒng)示意圖Fig.2 Coordinates of external field stitchingmeasuring system in earth centroid coordinate

        4 測量相機外視場拼接參數ξij的確定

        用文中4個測量相機的配置方式進行外視場拼接時,不可避免地會造成視差。當拼接參數ξij選擇不當時,在一定的測量距離上會產生測量盲區(qū)(即漏測),或者浪費測量相機的視場資源。為了實現每臺光電測量儀的4個測量相機視場在外場的合理拼接,可以讓中間的兩個(第二個和第三個)測量相機之間,在0.7l1距離處的外視場緊密拼接(確保不露縫),而在兩邊的兩個(第一個和第二個、第三個和第四個)測量相機之間為緊密搭接(即視場邊緣的像元既無重疊又無間隙,條件是第一和第二個,第三和第四個測量相機的節(jié)點,上下對齊)。如圖3所示。

        由圖3所示的幾何關系可知:

        由此可解得ξi2和ξi3,并進一步解得ξi1和ξi4:

        式中l(wèi)1為被測量目標到光電測量儀地平坐標系原點Si的最近距離;ξij為第i個光電測量儀,第j個測量相機光軸,相對第i個光電測量儀水平軸坐標系在方位上的偏轉角;dij1、dij2、dij3為第i個光電測量儀,第j個測量相機望遠鏡頭的節(jié)點,相對第i個光電測量儀水平軸坐標系原點Hi,在3個坐標軸上的位移距離;ωxij、ωyij分別為第i個光電測量儀,第j個測量相機在橫向和豎向的半視場角。

        圖3 測量相機排列示意圖Fig.3 Sketch of cameras ordering

        為了使各測量相機的各視場之間有一定的覆蓋,計算時可將各相機的有效視場角減小一定的量值,用來實現視場的覆蓋。

        5從第i個光電測量儀,第j個測量相機坐標系到第i個光電測量儀水平軸坐標系Hi(hi1,hi2,hi3)的變換

        設第k個目標到第i個光電測量儀,第j個測量相機坐標系Cij原點的距離為rijk,則該目標在第i個光電測量儀,第j個測量相機坐標系Cij中的測量值Cijk為:

        式中,k為被測目標的序號;Δαijk、Δλijk分別為第k個目標,在第i個光電測量儀,第j個測量相機坐標系中的方位角和俯仰角;xijk、yijk分別為第i個光電測量儀,第j個測量相機的面陣CCD上,第k個目標的橫向和豎向脫靶量(mm);f′ij為第i個光電測量儀,第j個測量相機鏡頭的焦距(mm)。

        設第k個目標到第一個光電測量儀,第j個測量相機坐標系C1j原點的距離為r1jk,則該目標在第一個光電測量儀第j個測量相機坐標系C1j中的測量值C1jk為:

        其中,

        設第k個目標經過視場拼接后,在第i個光電測量儀水平軸坐標系Hi(測量相機視場拼接后的組合測量相機坐標系)中,測量值Hik為:

        式中Rik為第k個目標到第i個光電測量儀,水平軸坐標系Hi原點Hoi的距離;Aik、Eik分別為第k個目標在第i個光電測量儀水平軸坐標系中的方位角和俯仰角。

        經如下的坐標轉換過程:

        (1)繞cij3軸旋轉(-ξij)角;

        (2)沿cij1軸平移(-dij1);

        (3)沿cij3軸平移(-dij3);

        (4)沿cij2軸平移(-dij2),

        可以實現第k個目標,在第j個測量相機坐標系Cij中的測量值到第i個光電測量儀水平軸坐標系Hi中測量值的變換。

        同理應用下述的坐標轉換過程:

        (1)繞ci13軸旋轉(-ξi1);

        (2)沿ci11軸平移(-di11);

        (3)沿ci13軸平移(-di13);

        (4)沿ci12軸平移(-di12),

        可以實現第k個目標,在第一個測量相機坐標系Ci1中的測量值,到第一個光電測量儀水平軸坐標系H1中的測量值的變換:

        6從第i個光電測量儀水平軸坐標系到第一個光電測量儀水平軸坐標系的變換,即Hi(hi1、hi2、hi3)到H1(h11、h12、h13)的變換

        經如下的坐標轉換過程:

        (1)從第i個光電測量儀水平軸坐標系Hi到第i個光電測量儀垂直軸坐標系Vi的坐標轉換:①繞hi2軸旋轉-λi角;

        (2)從第i個光電測量儀垂直軸坐標系Vi到第i個光電測量儀地平坐標系Si的坐標轉換:②繞si3軸旋轉-αi角;

        (3)從第i個光電測量儀地平坐標系Si到地球質心坐標系I的坐標轉換:③沿si3軸平移-hi;④繞si2軸旋轉-θi角;

        (4)從地球質心坐標系I到第1個光電測量儀地平坐標系S1的坐標轉換:⑤繞i1軸旋轉δiδ1角;⑥繞i2軸旋轉θ1角;⑦沿S13軸平移h1;

        (5)從第1個光電測量儀地平坐標系S1到第1個光電測量儀垂直軸坐標系V1的坐標轉換:⑧繞s13軸旋轉α1角;

        (6)從第1個光電測量儀垂直軸坐標系V1地平坐標系S1到第1個光電測量儀水平軸坐標系H1的坐標轉換:⑨繞v12軸旋轉λ1角。

        可以實現第k個目標,在第i個光電測量儀水平軸坐標系中的測量值Hik,到第一個光電測量儀地平坐標系H1中的測量值H1k的變換:

        式中:δi、θi、hi為第i個光電測量儀坐標系原點Soi處的赤經、赤緯和相對地球質心的高程(單位km);αi為第i個測量站光電測量儀垂直軸坐標系與地平坐標系(正北方向N)之間的夾角(即真方位角);λi為第i個測量站光電測量儀水平軸坐標系相對垂直軸坐標系的俯仰角。

        7 交匯測量和視場拼接結果

        7.1 交匯測量結果

        由于每臺光電測量儀的4臺測量相機之間都存在左右和上下的視差,僅僅應用測量相機的測量數據不可能完成在自身的光電測量儀中的視場拼接。如式(5)和(8)所示,要完成視場拼接首先應解出距離rijk和r1jk,也就是說首先要對2臺光電測量儀的測量相機的測量數據進行交匯[9~13],求得距離信息rijk和r1jk,即通過聯解方程組(15)、(16)求得rijk和r1jk。其中式(15)由式(12)和式(14)聯立而得;式(16)即式(13)。

        7.2 視場拼接結果

        通過解方程(12)可求出Rik,Eik和Aik,則在第i個光電測量儀中,視場拼接后第k個目標在第i個光電測量儀水平軸坐標系中(即視場拼接后在組合視場中)的測量值Hik為:

        由于:

        式中f′為測量相機視場拼接后的組合測量相機的標準鏡頭焦距(mm);xik,yik分別為在第i個光電測量儀中,經視場拼接后,第k個目標在第i個光電測量儀組合測量相機視場中的橫向和豎向脫靶量(mm)。

        由此可得第k個目標在第i個光電測量儀中,視場拼接后在組合視場中兩個方向上的脫靶量xik和yik為:

        通過解(13)可求出R1k,E1k和A1k,則在第1個光電測量儀中,第k個目標在第1個光電測量儀水平軸坐標系中(即視場拼接后在組合視場中)的測量值H1k為:

        由于:

        可得第k個目標在第1個光電測量儀中經視場拼接后在組合視場中兩個方向上的脫靶量x1k和y1k為:

        8 被測目標軌跡的交匯測量和視場拼接實例

        圖4 測量站1、2的地平坐標系,以及目標軌跡在測量站1地平坐標系中的三維圖Fig.4 Target path in geodetic coordinates atmeasurement station 1,and the geodetic coordinates at measurement station 1 and 2

        如圖4所示,圖中以測量站1的地平坐標系為基準,分別示出了測量站1、2的地平坐標系,以及目標在測量站1地平坐標系中的三維軌跡。

        測量站1和2的地平坐標系參數為:

        δ1=125.35°、θ1=43.88°、h1=6 371 000m;

        δ2=125.362 5°、θ2=43.88°、h2=6 371 005m。

        測量站1和2的電視測量儀水平軸坐標系(即組合視場1和2光軸)在各自地平坐標系中的方位和仰角參數為:α1=-35°、λ1=15°;α2=30°、λ2=15°。

        因為擬測目標區(qū)域與測量站之間的最短距離為500m,也就是說大于500m距離,被測目標可實現無盲區(qū)覆蓋,本文采用l1=500m進行拼接參數設計。

        由已知電視測量儀的下列結構參數:

        d111=d121=d131=d141=0mm;

        d211=d221=d231=d241=0mm;

        d112=d122=d212=d222=430mm;

        d132=d142=d232=d242=-430mm;

        d113=d143=d213=d243=180mm;

        d123=d133=d223=d233=-180mm;

        f′11=f′12=f′13=f′14=94mm;

        f′21=f′22=f′23=f′24=94mm;

        ωxij=4.98°、ωyij=4.98°,

        根據公式(1)、(2)、(3)、(4),得出電視測量儀1和2的其它結構參數為:

        ξ12=ξ22=4.9096°;

        ξ11=ξ21=14.8696°;

        ξ13=ξ23=-4.9096°;

        ξ14=ξ24=-14.8696°。

        圖5和圖6分別為在測量站1和2處的兩個高速電視測量儀的8臺測量相機C11,C12,C13,C14和C21,C22,C23,C24,用500frame/s的采樣頻率,對同一運動目標在視場中獲得的測量值,總測量時間為2s,總采樣數據為1 000組。

        圖7和圖8分別為根據測量站1和2的2個高速電視量儀的8臺測量相機對在視場中獲得的同一運動目標測量值進行交匯測量處理和視場拼接后在電視測量儀中拼接視場中的拼接結果(軌跡)。

        圖5 運動目標在測量站1中C11、C12、C13、C144個測量相機視場中的軌跡Fig.5 Target paths in FOV of C11,C12,C13,C14 at measurement station1

        圖6 運動目標在測量站2中C11、C12、C13、C144個測量相機視場中的軌跡Fig.6 Target paths in FOV of C11、C12、C13、C14 at measurement station 2

        圖7 運動目標在測量站1拼接視場中的軌跡Fig.7 Target path in FOV after field stitching ofmeasurement station 1

        圖8 運動目標在測量站2拼接視場中的軌跡Fig.8 Target path in FOV after field stitching ofmeasurement station 2

        9 結論

        能否實現高速電視測量儀外視場拼接,取決于硬件和拼接算法(軟件)兩個方面。本文在2臺高速電視測量儀原型樣機(硬件)研制成功的基礎上,提出了一種將交匯測量和拼接處理相結合的算法,通過用2臺高速電視測量儀的8臺測量相機測量同一運動目標,并對獲得的測量數據進行交匯和視場拼接處理,得出以下結論:

        (1)拼接算法是正確、有效的,可以得到唯一解。(2)要完成視場拼接,首先要根據交匯測量原理,對數據進行交匯測量處理[9~13],解出目標的距離信息。(3)本文中提供的交匯測量拼接算法,由于引入地球質心坐標系變換,地球曲率半徑和地球子午線收斂角兩項影響因素得到了完全修正(當兩站測量儀在同一緯度上時,兩個站點之間的地球子午線收斂角最大,設兩個站點測量儀在緯度為43.88°,相距1km時,其地球子午線收斂角約為45″,兩站地平面鉛垂線之間的夾角約為32.5″)。(4)本文中的交匯測量算法,完全可以推廣到其它光電測控儀器的交匯測量,如果擬測量的目標越遠,或測量站間布站距離越大,則這種算法的效果和優(yōu)點就更明顯。(5)視場拼接和交匯測量的總誤差,主要來自兩個方面:一是光電測量儀本身的測量誤差;二是布站點處地平坐標系原點位置的測量誤差。外視場拼接和交匯測量算法本身不引入誤差。

        需要注意的是,布站時任意2臺光電測量儀到被測目標區(qū)域之間的夾角不要過小,最好符合正交布站原則,否則由于光電測量儀本身的測量誤差和布站點位置測量誤差的影響,會導致交匯測量誤差增大,甚至造成交匯無解。

        [1]劉罡,高曉東,曹學東.大視場光電測量系統(tǒng)實時測量精度的綜合標定[J].光電工程,2001,28(6):10-13.

        LIU G,GAO X D,CAO X D.A comprehensive calibration of the real-time measurement accuracy for a photoelectrical measuring system with wide field of view[J].Opto-electronic Eng.,2001,28(6):10-13.(in Chinese)

        [2]于曉波,盛磊.靶場測量中多鏡頭大視場視頻圖像的拼接[J].光學精密工程,2008,16(11):2145-2150.

        YU X B,SHENG L.Image stitching of multi-lens with large visual field in range instrumentation[J].Opt.Precision Eng.,2008,16(11):2145-2150.(in Chinese)

        [3]賈峰,李桂芝,南雪飛.靶場T型架光電經緯儀精度檢測方法的研究與應用[J].飛行器測控學報,2007,26(2):24-29.

        JIA F,LIG ZH,NAN X F.An accuracy testmethod for type-T cinetheodolites in instrum entation ranges[J].J.Spacecraft TT&C Technol.,2007,26(2):24-29.(in Chinese)

        [4]何照才,胡保安.光學測量系統(tǒng)[M].北京:國防工業(yè)出版社,2002.

        HE ZH C,HU B A.Optical Measure Systems[M].Beijing:National Defence Industry Press,2002.(in Chinese)

        [5]趙學顏,李迎春.靶場光學測量[M].北京:國防工業(yè)出版社,2001.

        ZHAO X Y,LIY CH.Optical Measure at Shooting Range[M].Beijing:National Defense Press,2001.(in Chinese)

        [6]王家騏.光學儀器總體設計[M].長春:中科院長春光學精密機械與物理研究所,2003.

        WANG JQ.Optical Instrument Collectivity Design[M].Changchun:Changchun Institute of Optics,Fine Mechanics and Physics,Chinese Academy of Sciences,2003.(in Chinese)

        [7]金光,王家騏.利用坐標變換推導經緯儀三軸誤差[J].光學精密工程,1999,7(5):89-94.

        JIN G,WANG JQ.The three-error axis of theodolitewith the utilization of the coordinate to the variation[J].Opt.Precision Eng.,1999,7(5):89-94.(in Chinese)

        [8]王家騏,金光,顏昌翔.機載光電跟蹤測量設備的目標定位誤差分析[J].光學精密工程,2005,13(2):105-116.

        WANG JQ,JIN G,YAN CH X.Orientation error analysis of airborne opto-electric tracking and measuring device[J].Opt.Precision Eng.,2005,13(2):105-116.(in Chinese)

        [9]馮桂蘭,田維堅,屈有山,等.實時視場拼接系統(tǒng)的設計與實現[J].光電工程,2007,34(4):124-127.

        FENG G L,TIAN W J,QU Y SH,et al..Design of real-time video mosaic system[J].Opto-electronic Eng.,2007,34(4):124-127.(in Chinese)

        [10]任建岳,孫斌,張星祥,等.TDICCD交錯拼接的精度檢測[J].光學精密工程,2008,16(10):1852-1857.

        REN JY,SUN B,ZHANG X X,et al..Precision measurement of TDICCD interleaving assembly[J].Opt.Precision Eng.,2008,16(10):1852-1857.(in Chinese)

        [11]蘭海濱,王平,龍騰.圖像拼接中相機鏡頭非線性畸變的校正[J].光學精密工程,2009,17(5):1196-1202.

        LAN H B,WANG P,LONG T.Nonlinear aberration correction of lens in imagemosaic[J].Opt.Precision Eng.,2009,17(5):1196-1202.(in Chinese)

        [12]李喆,丁振良,袁峰.基于共面點的多視覺測量系統(tǒng)的全局標定[J].光學精密工程,2008,16(3):467-472.

        LIZH,DING ZH L,YUAN F.Global calibration method formulti-vision measurement system with coplanar targets[J].Opt.Precision Eng.,2008,16(3):467-472.(in Chinese)

        [13]費業(yè)泰.誤差理論與數據處理[M].北京:機械工業(yè)出版社,2007.

        FEIY T.Error Theory and Data Processing[M].Beijing:China Machine Press,2007.(in Chinese)

        Convergence algorithm of field stitching and intersection measurement for external field stitching measuring system and its im plementation

        WANG Min1,2,SONG Li-wei1,QIAO Yan-feng1,YU Yi1
        (1.Changchun Institute of Optics,Fine Mechanics and Physics,Chinese Academy of Sciences,Changchun 130033,China;2.Graduate University of Chinese Academy of Sciences,Beijing 100039,China)

        A convergence algorithm of field stitching and intersection measurement is proposed based on the development of two prototypes of high-speed video-measuring systemswith external field stitching.Five coordinates are established for the earth centroid and optoelectronical systems.The definations of the coordinates are introduced and coordinate transforms are discussed.The experiment of intersection measurement and field stitching for an identicalmoving target is carried out with the data from eight cameras in two high-speed TV measuring instruments.The result shows that the stitching algorithm is correct and effective,and has obtained an unique solution.Meanwhile,the influences of earth curvature and meridional convergent errors are fullycorreated,which indicates that the intersection measurement algorithm can also be used in the intersection measurement of other optoelectronicalmeasurement and control instruments.

        external field stitching,intersection measurement,earth curvature,meridional convergence;error correction

        1674-2915(2010)03-0229-10

        V557.3;TP391.4

        A

        2010-03-12;

        2010-04-13

        王 旻(1977—),女,吉林長春人,碩士研究生,主要從事光電瞄準與測量的研究。E-mail:work6180@yahoo.com.cn

        喬彥峰(1962—),男,吉林長春人,研究員,博士生導師,主要從事光電測量與控制技術方面的研究。E-mail:qiaoyf@ciomp.ac.cn

        猜你喜歡
        交匯測量儀視場
        星模擬器光學系統(tǒng)視場拼接方法的研究
        中國光學(2021年6期)2021-11-25 07:48:32
        水平度與垂直度精密測量儀
        兩水交匯,涇渭分明
        三角的交匯創(chuàng)新問題
        基于單片機的便捷式LCF測量儀
        電子制作(2019年9期)2019-05-30 09:42:02
        揭秘身高體重測量儀
        聚焦離散型隨機變量的交匯創(chuàng)新
        醫(yī)用內窺鏡矩形視場下入瞳視場角的測試方法研究
        寬電容測量儀的設計
        電子制作(2017年7期)2017-06-05 09:36:14
        輕小型面陣擺掃熱紅外成像系統(tǒng)研究
        97一期涩涩97片久久久久久久| 亚洲AV秘 无码一区二区三区| 日本国主产一区二区三区在线观看| 日本人妻系列中文字幕| 亚洲av中文无码乱人伦在线咪咕| 人妻暴雨中被强制侵犯在线| 大陆一级毛片免费播放| 亚洲av日韩一区二三四五六七 | av免费在线播放视频| 久久久av波多野一区二区| 亚洲最大天堂无码精品区| 熟女少妇av免费观看| 一区二区三区四区亚洲免费 | 女女互揉吃奶揉到高潮视频| 秋霞午夜无码鲁丝片午夜精品| 国产一区二区丁香婷婷| 第一九区另类中文字幕| 极品少妇小泬50pthepon| 亚洲国产成人精品女人久久久| 国产精品亚洲精品日产久久久| 九九久久精品国产免费av| 人妻少妇乱子伦精品| 精品欧美乱子伦一区二区三区 | 精品在线亚洲一区二区三区| 日韩av无码社区一区二区三区| 日韩在线一区二区三区免费视频 | 人妻少妇精品中文字幕专区| 在线精品国产一区二区三区| av无码天堂一区二区三区| 狠狠综合久久av一区二区三区| 国产欧美性成人精品午夜| 99久热re在线精品99 6热视频| 日本一区二区亚洲三区| 青青草免费在线爽视频| 久久久久亚洲av片无码v| 国产乱人伦真实精品视频| 中文字幕人妻精品一区| 三级全黄的视频在线观看| 一区二区日韩国产精品| 成人av在线免费播放| 无码国产精品一区二区免费式芒果 |