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

    基于QTM的導(dǎo)航星庫劃分方法

    2016-12-15 03:49:16趙啟坤嚴(yán)小軍
    導(dǎo)航與控制 2016年6期
    關(guān)鍵詞:天球星圖經(jīng)緯度

    趙啟坤,嚴(yán)小軍

    (北京航天控制儀器研究所,北京100039)

    基于QTM的導(dǎo)航星庫劃分方法

    趙啟坤,嚴(yán)小軍

    (北京航天控制儀器研究所,北京100039)

    為了提高星敏感器觀測星的檢索速度,提出了一種基于球面四元三角網(wǎng)(Quaternary Triangular Mesh,QTM)的導(dǎo)航星庫劃分方法。該方法首先利用天球的內(nèi)接正八面體把全天球劃分成8個子天區(qū),然后在每個子天區(qū)中進(jìn)行與經(jīng)緯度相關(guān)的剖分和編碼,最后掃描恒星星表,把每顆導(dǎo)航星信息劃歸到相應(yīng)的子天區(qū)內(nèi)存儲。利用此星庫進(jìn)行的局部天球星圖識別仿真實驗表明該星庫具有很高的檢索效率。

    星敏感器;導(dǎo)航星星庫;球面四元三角網(wǎng);星圖識別

    0 引言

    在星敏感器星圖識別中,合適的導(dǎo)航星庫組織結(jié)構(gòu)可以快速從導(dǎo)航星表中提取出需要檢索的導(dǎo)航星信息,以提高星圖識別效率。在組織導(dǎo)航星庫時一般采取的方法是把天球按一定的規(guī)則劃分成若干子區(qū)域,與某一子區(qū)域有關(guān)的識別模式用相應(yīng)的數(shù)據(jù)文件分別存儲或用數(shù)據(jù)文件的特定部分存儲。然而目前主流的劃分方法還存在兩個不足:一方面各個子區(qū)域的覆蓋面積的差異會導(dǎo)致各個子區(qū)域內(nèi)導(dǎo)航星數(shù)目存在差別,導(dǎo)航星信息在各個子塊中的搜索效率有較大差異;另一方面各個子區(qū)域缺乏高效的索引機制,赤經(jīng)、赤緯與子區(qū)域地址碼的轉(zhuǎn)換比較繁瑣,子區(qū)域的鄰近索引不能直接利用地址碼進(jìn)行。

    本文首先采用基于球面四元三角網(wǎng)的空間剖分方法對導(dǎo)航星庫進(jìn)行剖分,以保證剖分出的各個子塊形狀和面積基本相同;然后對各個子塊進(jìn)行與經(jīng)緯度相關(guān)的編碼,提高子塊的檢索效率和鄰近搜索效率;最后利用局部天區(qū)星圖識別算法對模擬的150幅星圖進(jìn)行識別,驗證本方法的性能。

    1 常用的導(dǎo)航星庫劃分方法

    現(xiàn)有的導(dǎo)航星庫劃分方法主要有赤緯帶法[1]、圓錐法[2]、球矩陣法[3]、內(nèi)接正六面體法等方法[4]。各種劃分方法和優(yōu)缺點如表1所示。

    表1 現(xiàn)有的導(dǎo)航星庫劃分方法Table 1 Star catalog partition methods

    2 基于球面四元三角網(wǎng)(QTM)的導(dǎo)航星庫劃分方法

    本文提出了一種基于球面四元三角網(wǎng)的導(dǎo)航星庫劃分方法,本方法主要包括以下幾個步驟:天球剖分、利用QTM對剖分后的天區(qū)進(jìn)行編碼、利用QTM碼在全天球范圍內(nèi)進(jìn)行鄰近檢索。

    2.1 天球剖分方法

    1)用天球的內(nèi)接正八面體把球面均勻的分成8個區(qū)域S1~S8[5],每個子區(qū)域稱為一個八分體,如圖1所示。正八面體的6個頂點在天球坐標(biāo)系中的坐標(biāo)如表2所示。

    圖1 內(nèi)接正八面體剖分球面Fig.1 Spherical facet based on octahedron

    表2 正八面體頂點在天球坐標(biāo)系中的坐標(biāo)Table 2 Vertex coordinates ofregular octahedron in celestial coordinate system

    正八面體各個邊在天球上的投影與赤道、0°經(jīng)線、180°經(jīng)線、東經(jīng)90°線、西經(jīng)90°線重合,可以方便地確定天球上的某一點在哪個八分體上。

    圖2 三角形剖分示意圖Fig.2 Triangle subdivisions

    2)對于S1~S88個八分體,按照圖2所示的方法將其劃分為4N個小區(qū)域。整個天球就被劃分為8×4N個小區(qū)域。

    如圖3所示,對球面三角形進(jìn)行遞歸剖分的方法分為兩種,一種是取大圓弧的中點(圖3中的M1點)作為新剖分三角形的一個頂點,這種方法稱為大弧平分法;另一種方法是取經(jīng)緯線的中點(圖3中的M2點)作為新剖分三角形的一個頂點,這種方法被稱為經(jīng)緯度平分法。大弧平分法得到的球面三角形比較均勻,形狀、大小也最接近,缺點是計算比較復(fù)雜。經(jīng)緯度平分法雖然在三角形的相似性上有一些損失,但是保證了數(shù)據(jù)與經(jīng)緯度關(guān)聯(lián)的特性,有利于數(shù)據(jù)的編碼和索引。在導(dǎo)航星庫中星的坐標(biāo)都是與經(jīng)緯度相關(guān)聯(lián)的,因此為了快速索引導(dǎo)航星,本文采用經(jīng)緯度平分法對球面三角形的各個面進(jìn)行細(xì)分。

    圖3 大弧平分法和經(jīng)緯度平分法Fig.3 The methods of triangulation

    2.2 利用QTM對天球進(jìn)行編碼

    在利用QTM對天球編碼時,一個QTM位置碼由一個八分碼(0~7)和若干個四分碼(0~3)組成。在第k個剖分層次,球面三角形A的編碼可以表示為:a0a1a2…ak。其中,a0是八分碼,a1~ak是k個四分碼。

    八分碼a0代表8個八分體之一,這8個八分體的編碼如圖4所示。

    圖4 內(nèi)接正八面體剖分球面的初始編碼Fig.4 Initial partition and coding of spherical facet based on octahedron

    四分碼ak(k≥1)可以通過當(dāng)前區(qū)域在其父區(qū)域中的相對位置來編碼:0表示當(dāng)前區(qū)域在其父區(qū)域的中間,該三角形區(qū)域稱為中間三角形;1表示當(dāng)前區(qū)域在其父區(qū)域的上面(或者下面),該三角形區(qū)域稱為頂邊三角形;2表示當(dāng)前區(qū)域在其父區(qū)域的左邊,該三角形區(qū)域稱為左邊三角形;3表示當(dāng)前區(qū)域在其父區(qū)域的右邊,該三角形區(qū)域稱為右邊三角形,如圖5(a)所示。圖5(b)是對為0的八分體作為初始球面三角形進(jìn)行2次剖分后的QTM編碼示意圖。其余7個區(qū)域的剖分與0區(qū)域類似。

    圖5 QTM編碼示意圖Fig.5 Subdivision and QTM coding

    2.3 QTM地址碼與赤經(jīng)赤緯的相互轉(zhuǎn)換

    天球上的每一個星點都由赤經(jīng)、赤緯來確定位置,因此赤經(jīng)赤緯與QTM編碼之間的相互轉(zhuǎn)換是完成星點快速搜索的關(guān)鍵。

    在現(xiàn)有的QTM編碼與赤經(jīng)、赤緯的轉(zhuǎn)換算法中,ZOT(Zenithal Ortho Triangular)投影法轉(zhuǎn)換速度快,但生成的編碼缺乏方向性;ETP(Equal?Triangles Projection)投影法生成的編碼具有固定的方向性,但轉(zhuǎn)換速度慢。趙學(xué)勝等提出的行列逼近法(Cavalcade Approach Method,CAM)[6]是一種從赤經(jīng)、赤緯到QTM碼的快速轉(zhuǎn)換方法。該方法基于經(jīng)緯度平分法,根據(jù)位置區(qū)域的行列以及地址碼的方向進(jìn)行層次遞歸,轉(zhuǎn)換速度較快,因此本文選擇行列逼近法進(jìn)行QTM地址碼與赤經(jīng)、赤緯的相互轉(zhuǎn)換。

    2.4 基于QTM編碼的鄰近檢索

    在QTM格網(wǎng)中鄰近三角形的定義為:具有公共邊的為邊鄰近Edge?adjacent(三角形);只具有公共頂點的為角鄰近(也稱為頂點鄰近,Vertex?adja?cent)三角形。

    (1)邊鄰近三角形的搜索

    在QTM格網(wǎng)中所有三角形具有3個邊鄰近:左鄰近三角形L(u)、右鄰近三角形R(u)和頂(上或下)鄰近三角形T(u),每種鄰近又可以按照鄰近的三角形是否在同一個八分體內(nèi)分為兩種情況。

    下面以左鄰近搜索為例介紹邊鄰近三角形的搜索算法,如圖6所示。

    圖6 對八分體剖分兩次的左鄰域搜索示意圖Fig.6 Searching sketch map of octant splited two times on the left neighborhood

    如果左鄰近三角形在同一個八分體內(nèi)部:

    1)如果當(dāng)前三角形區(qū)域地址碼的最后一位是0,也即當(dāng)前三角形是其父三角形的中間三角形,其左鄰近三角形為其父三角形的左邊三角形;

    2)如果當(dāng)前三角形區(qū)域地址碼的最后一位是3,也即當(dāng)前三角形是其父三角形的右邊三角形,其左鄰近三角形為其父三角形的中間三角形;

    3)如果當(dāng)前三角形區(qū)域地址碼的最后一位是1,也即當(dāng)前三角形是其父三角形的頂邊三角形,其左鄰近三角形為其父三角形的左鄰近三角形的右邊三角形;

    4)如果當(dāng)前三角形區(qū)域地址碼的最后一位是2,也即當(dāng)前三角形是其父三角形的左邊三角形,其左鄰近三角形為其父三角形的左鄰近三角形的頂邊三角形。

    如果左鄰近三角形不在同一個八分體內(nèi)部:

    1)計算當(dāng)前八分體左側(cè)的八分體的八分碼:

    2)所有為2的四分碼變?yōu)?,其余四分碼保持不變:

    這樣就得到了三角形左鄰近三角形的QTM編碼,算法的流程圖如圖7所示。

    對于右鄰近和頂鄰近搜索有類似的流程,在此不再贅述。

    (2)角鄰近三角形的搜索

    三角形在八分體的最上、最左、最右側(cè)時,角鄰近三角形有7個,如圖8(a)所示;在其他位置時,角鄰近三角形有9個,如圖8(b)所示。

    對角鄰近三角形搜索可以轉(zhuǎn)換為邊鄰近三角形的搜索。以圖8(b)為例,編號為1的三角形與當(dāng)前三角形(U)的關(guān)系可以表示為:

    式中,U1代表編號為1的三角形,U代表當(dāng)前三角形,代表當(dāng)前三角形的頂鄰近三角形,代表當(dāng)前三角形的左鄰近三角形。

    (3)對一個三角形的鄰近搜索包括邊鄰近搜索和角鄰近搜索

    由圖7和圖8可知,三角形的邊鄰近搜索只是簡單的二進(jìn)制操作,而角鄰近搜索可以轉(zhuǎn)換為7個(或9個)邊鄰近搜索。因此,鄰近搜索的效率非常高。已知一個天區(qū)的QTM編碼時,可以通過鄰近搜索算法迅速計算出與該三角形相鄰的所有三角形的QTM編碼。

    圖8 角鄰近三角形Fig.8 Vertex?adjacent triangles

    3 基于QTM劃分的導(dǎo)航星庫在星圖識別中的應(yīng)用與分析

    以10°×10°的視場為例,利用QTM法劃分導(dǎo)航星庫,并基于劃分后的星庫做局部天區(qū)星圖識別,以驗證該劃分方法的性能。

    對每一個八分體Si利用經(jīng)緯度平分法細(xì)分3次,全天球區(qū)域被劃分成了512(8×43)個小三角形天區(qū);掃描導(dǎo)航星庫,將每顆導(dǎo)航星都劃歸到相應(yīng)的子天區(qū);按照QTM編碼的方式為每一個天區(qū)建立如表3所示的索引表。

    表3 天區(qū)索引表Table 3 Index of sub?sky

    對其中一個八分體進(jìn)行統(tǒng)計分析,球面三角形頂點之間夾角的最小值為11.25°,頂點與對邊中點夾角的最小值為7.94°,如表4所示。因此,任意一個視場角小于15.88°×15.88°的視場都可以利用該視場中心所在的三角形天區(qū)及其12鄰近(或10鄰近)天區(qū)完全覆蓋。

    表4 八分體內(nèi)三角形頂點夾角大小的統(tǒng)計分析Table 4 Statistic and analysis for the angles in octant

    在載體運動角速度比較小時,利用上一時刻視軸的所在天區(qū)及其鄰近天區(qū)中包含的導(dǎo)航信息組成實時導(dǎo)航星庫;當(dāng)載體運動角速度比較大時,利用之前幾個時刻的視軸方向和運動信息估計當(dāng)前時刻的視軸方向,組建實時導(dǎo)航星庫。利用這種方法把星圖識別的搜索范圍縮減為整個導(dǎo)航星庫的13/512≈2.5%。

    采用QTM劃分導(dǎo)航星庫后,星圖識別的流程如下:

    1)第一幀星圖的視軸方向由全天球星圖識別計算得出(或者由程序給定);

    2)根據(jù)上一時刻的視軸方向或者之前幾個時刻的視軸方向構(gòu)建(或更新)實時星庫;

    3)利用實時星庫,進(jìn)行局部天區(qū)星圖識別;

    4)如果局部天區(qū)星圖識別未能成功識別則轉(zhuǎn)入全天球星圖識別。

    下面以三角形識別為例,對視軸為BoreSight1(Ra=282.09°,Dec=-28.77°,Roll=40.06°,10°×10°圓形視場,極限星等為6.0等)的視場分別采用全天球星圖識別法與QTM劃分導(dǎo)航星庫的局部天區(qū)星圖識別法進(jìn)行識別以比較兩種方法在識別效率上的差異,識別中星角距的匹配誤差設(shè)置為0.001rad(0.057°)。在進(jìn)行局部天區(qū)星圖識別之前,首先對視軸為BoreSight0(Ra=280.60°,Dec =-28.95°,Roll=38.94°)的視場進(jìn)行全天球星圖識別以建立實時星庫。

    在BoreSight1視場中共有10顆導(dǎo)航星,選取其中最亮的3顆待識別星(S1、S2、S3)組成的3個待匹配星對和待匹配星三角形進(jìn)行匹配,待匹配星對和待匹配星三角形的匹配結(jié)果如表5所示。

    表5 對BoreSight1視場中最亮的三顆星進(jìn)行識別的結(jié)果Table 5 Identify result of the three brightest stars in BoreSight1

    由表5可以看出,采用QTM對導(dǎo)航星庫進(jìn)行劃分后,一方面大大減少了與每個待匹配星對相匹配的星對數(shù)量;另一方面減少了待匹配三角形出現(xiàn)冗余匹配的可能性,算法的運行時間得以減少。

    為了評價該星庫的性能,在計算機上進(jìn)行了仿真測試。測試時假設(shè)星敏感器初始姿態(tài)為Bore?Sight0(Ra=280.60°,Dec=-28.95°,Roll= 38.94°),第i時刻姿態(tài)為:

    式中,j取0、1、2,分別代表視軸的Ra、Dec、Roll 3個方向;rand是范圍為(-1°,1°)的隨機數(shù)。

    測試總共進(jìn)行了150步,每一步都會根據(jù)當(dāng)前的姿態(tài)產(chǎn)生一幅星圖模擬星敏感器采集到的觀測星圖,測試結(jié)果如圖9所示。

    圖9 150次星圖識別所用時間Fig.9 Elapsed time of star identification

    除去第1次識別采用全天球星圖識別外,其余149次均用QTM分區(qū)的局部天區(qū)星圖識別完成,證明了QTM分區(qū)、編碼及鄰域搜索算法的正確性。采用QTM分區(qū)的局部天區(qū)星圖識別每次識別時間最短為0.34ms,最長為8.5ms,平均識別時間為1.9ms,且有84%的匹配時間小于2.0ms,相對于全天球星圖識別來說識別時間少了1個數(shù)量級。以文獻(xiàn)[4]為例,采用內(nèi)接正六面體劃分導(dǎo)航星表的改進(jìn)三角形星圖識別方法做一次全天球星圖識別需要的平均時間為8.4ms。

    4 結(jié)論

    本文把球面四元三角網(wǎng)(QTM)引入到導(dǎo)航星庫的劃分中,提出了一種基于經(jīng)緯度平分法的QTM導(dǎo)航星庫劃分方法,并采用行列逼近算法對各個天區(qū)進(jìn)行編碼?;谶@種劃分方法構(gòu)造的導(dǎo)航星庫繼承了QTM劃分均勻、編碼方向一致、鄰域搜索高效,經(jīng)緯度平分法與經(jīng)緯度高度相關(guān),行列逼近法編碼快速等特點。基于該星庫進(jìn)行的局部天區(qū)星圖識別仿真試驗也驗證了這種方法的高檢索效率。

    [1]Bone J W.On?orbit star processing using multi?star star trackers[J].Proceedings of the SPIE,1994:6?14.

    [2]Ju G,Kim H,Pollock T,et al.Digistar:a low?cost micro star tracker[C].AIAA Space Technology Conference&Exposition,Albuquerque,1999.

    [3]Roelof W H,Van B.True?sky demonstration of an auton?omous star tracker[J].Proceedings of the SPIE,1994:204?216.

    [4]張廣軍.星圖識別[M].北京:國防工業(yè)出版社,2011.ZHANG Guang?jun.Staridentification[M].Beijing:Na?tional Defense Industry Press,2011.

    [5]Randall D A,RinglerTD,HeikesRP.Climate modeling with spherical geodesic grids[J].Computing in Science and Engineering,2002,4(5):32?40.

    [6]趙學(xué)勝,陳軍.QTM地址碼與經(jīng)緯度坐標(biāo)的快速轉(zhuǎn)換算法[J].測繪學(xué)報,2003,32(3):272?277.ZHAO Xue?sheng,CHEN Jun.Fast translation algorithm between QTM code and longitude/latitude coordination[J].Acta Geodaeticaet Cartographica Sinica,2003,32(3):272?277.

    [7]白建軍,孫文彬,趙學(xué)勝.基于QTM的WGS?84托球面層次剖分及其特點分析[J].測繪學(xué)報,2011,40(2):243?248.BAI Jian?jun,SUN Wen?bin,ZHAO Xue?sheng.Charac?teranalysis and hierarchical partition of WGS?84 ellipsoidal facet based on QTM[J].Acta Geodaetica et Cartographica Sinica,2011,40(2):243?248.

    [8]趙學(xué)勝,侯妙樂,白建軍.全球離散格網(wǎng)的空間數(shù)學(xué)建模[M].北京:測繪出版社,2007.ZHAOXue?sheng,HOUMiao?le,BAIJian?jun.Spatialdata modeling based on global discrete grid[M].Beijing:Surveying and Mapping Press,2007.

    [9]孔祥元,郭際明,劉宗全.大地測量學(xué)基礎(chǔ)[M].武漢:武漢大學(xué)出版社,2001.KONG Xiang?yuan,GUO Ji?ming,LIU Zong?quan.Foun?dationofgeodesy[M].Wuhan:WuhanUniversity Press,2001.

    A Method for Guide Star Catalog Partition Based on QTM

    ZHAO Qi?kun,YAN Xiao?jun
    (Beijing Institute of Aerospace Control Devices,Beijing 100039)

    In order to improve the efficiency of star sensor during retrieving the guide star,a method for star catalog partition which based on quaternary triangular mesh was proposed.The approach first utilizes the regular inscribed regular octahedron in the celestial sphere and divided the sphere into eight sub regions,then split and encoded in each sub regions associated with the latitude and longitude,finally scanning the guide star database,each navigation information was storage in the corresponding sub region.The star identification used this star catalog shows that the designed star catalog has high retrieval efficiency.

    star sensor;guide star catalog;quaternary triangular mesh(QTM);star identification

    U666.1

    A

    1674?5558(2016)03?01223

    10.3969/j.issn.1674?5558.2016.06.014

    2016?01?04

    趙啟坤,男,碩士,導(dǎo)航、制導(dǎo)與控制專業(yè),研究方向為星光組合導(dǎo)航技術(shù)。

    猜你喜歡
    天球星圖經(jīng)緯度
    星圖上非線性分?jǐn)?shù)階微分方程邊值問題解的存在唯一性
    乾隆款景泰藍(lán)花開富貴 加座獸足天球瓶
    收藏界(2019年3期)2019-10-10 03:16:30
    詩意聯(lián)結(jié) 水漾星圖——上海龍湖·星圖美學(xué)展示中心
    天球瓶史話
    收藏界(2018年4期)2018-10-12 00:57:20
    基于三角形周長的暗星全天球自主快速識別
    自制中學(xué)實驗操作型經(jīng)緯測量儀
    澳洲位移大,需調(diào)經(jīng)緯度
    一種利用太陽影子定位的數(shù)學(xué)模型
    天文測量仿真器模擬星圖精度分析
    《地球概論》課程中“天球坐標(biāo)”的教學(xué)探究
    久久天躁狠狠躁夜夜2o2o| 色综合站精品国产| 免费无遮挡裸体视频| 久久久久久久亚洲中文字幕| 精品99又大又爽又粗少妇毛片 | 桃红色精品国产亚洲av| 97超视频在线观看视频| 在线观看av片永久免费下载| 99九九线精品视频在线观看视频| 久久久久性生活片| 日本a在线网址| 免费电影在线观看免费观看| 此物有八面人人有两片| 成人av一区二区三区在线看| 亚洲午夜理论影院| 性欧美人与动物交配| 日本免费a在线| 麻豆久久精品国产亚洲av| 欧美国产日韩亚洲一区| 久久午夜福利片| 精品乱码久久久久久99久播| 中文亚洲av片在线观看爽| 级片在线观看| 99在线视频只有这里精品首页| 国产91精品成人一区二区三区| 亚洲成av人片在线播放无| 丰满乱子伦码专区| 久久久久久久久大av| 精品一区二区免费观看| 亚洲精品在线观看二区| 国产精品综合久久久久久久免费| 亚洲av二区三区四区| 欧美+日韩+精品| 午夜久久久久精精品| 精品人妻熟女av久视频| 国产一区二区三区在线臀色熟女| 3wmmmm亚洲av在线观看| 国产亚洲精品久久久久久毛片| 免费一级毛片在线播放高清视频| 免费看av在线观看网站| 午夜精品在线福利| 亚洲aⅴ乱码一区二区在线播放| 国产精品爽爽va在线观看网站| 91av网一区二区| 我的老师免费观看完整版| 国产精品综合久久久久久久免费| 亚洲欧美日韩高清专用| 中文字幕人妻熟人妻熟丝袜美| 国产高清激情床上av| 成人毛片a级毛片在线播放| 中文字幕高清在线视频| 欧美3d第一页| 国产成人影院久久av| 99riav亚洲国产免费| 日韩欧美一区二区三区在线观看| 精品一区二区三区视频在线观看免费| 亚洲在线自拍视频| 小蜜桃在线观看免费完整版高清| 人妻制服诱惑在线中文字幕| 中文资源天堂在线| 成人一区二区视频在线观看| 日韩一区二区视频免费看| 欧美潮喷喷水| 真人做人爱边吃奶动态| 亚洲欧美日韩高清在线视频| 国产伦在线观看视频一区| ponron亚洲| 男人舔女人下体高潮全视频| 国产伦人伦偷精品视频| 精品免费久久久久久久清纯| 成人av在线播放网站| 国产精品久久电影中文字幕| 亚洲国产日韩欧美精品在线观看| 欧美日韩黄片免| www.色视频.com| 综合色av麻豆| 国产成人影院久久av| 国产久久久一区二区三区| 在线国产一区二区在线| 国产视频一区二区在线看| 婷婷精品国产亚洲av| av.在线天堂| 在线天堂最新版资源| 日韩国内少妇激情av| 亚洲av中文av极速乱 | 深夜a级毛片| 一级黄色大片毛片| 精品国内亚洲2022精品成人| 如何舔出高潮| 色av中文字幕| 女人被狂操c到高潮| 久9热在线精品视频| 一边摸一边抽搐一进一小说| 一区二区三区四区激情视频 | 国产大屁股一区二区在线视频| av女优亚洲男人天堂| 99热只有精品国产| 麻豆久久精品国产亚洲av| 久久人妻av系列| 97超级碰碰碰精品色视频在线观看| 亚洲欧美日韩卡通动漫| 久久久久久国产a免费观看| 国产乱人视频| av福利片在线观看| 日韩欧美在线乱码| 国内久久婷婷六月综合欲色啪| 男女视频在线观看网站免费| 久久人人精品亚洲av| 一区二区三区高清视频在线| 欧美一区二区精品小视频在线| 变态另类成人亚洲欧美熟女| 国产熟女欧美一区二区| av福利片在线观看| 国产精品国产高清国产av| 日本一本二区三区精品| 男人狂女人下面高潮的视频| 能在线免费观看的黄片| 久9热在线精品视频| 99久久精品国产国产毛片| 神马国产精品三级电影在线观看| 国产在线男女| 日本黄色视频三级网站网址| 午夜免费成人在线视频| 日韩精品有码人妻一区| 综合色av麻豆| 天堂动漫精品| 国产真实伦视频高清在线观看 | 欧美精品国产亚洲| 国产老妇女一区| 联通29元200g的流量卡| 久久久久九九精品影院| 欧美激情在线99| 免费大片18禁| 69av精品久久久久久| 变态另类丝袜制服| 亚洲av免费在线观看| 国产主播在线观看一区二区| 久久精品国产清高在天天线| 国产高清不卡午夜福利| 少妇的逼水好多| 国产久久久一区二区三区| 日韩亚洲欧美综合| 啪啪无遮挡十八禁网站| 男人的好看免费观看在线视频| 中国美女看黄片| 免费观看的影片在线观看| 99精品久久久久人妻精品| 免费观看的影片在线观看| 亚洲精品久久国产高清桃花| 欧美黑人巨大hd| 99久久无色码亚洲精品果冻| 最近中文字幕高清免费大全6 | 欧美黑人巨大hd| 99久久精品国产国产毛片| 成年人黄色毛片网站| 久久九九热精品免费| 97热精品久久久久久| 桃红色精品国产亚洲av| 午夜福利视频1000在线观看| 乱系列少妇在线播放| a在线观看视频网站| 久9热在线精品视频| 午夜爱爱视频在线播放| 亚洲国产色片| 无人区码免费观看不卡| 麻豆av噜噜一区二区三区| 亚洲成a人片在线一区二区| 中文字幕人妻熟人妻熟丝袜美| 午夜免费男女啪啪视频观看 | 天堂网av新在线| 国产精品野战在线观看| 熟妇人妻久久中文字幕3abv| 少妇高潮的动态图| 日韩中文字幕欧美一区二区| 身体一侧抽搐| 又粗又爽又猛毛片免费看| 精品久久久久久久久av| 亚洲国产欧洲综合997久久,| 最近视频中文字幕2019在线8| 欧美成人a在线观看| 97热精品久久久久久| 一进一出抽搐动态| 免费看av在线观看网站| 最近最新中文字幕大全电影3| 欧美最黄视频在线播放免费| 亚洲性久久影院| 日本成人三级电影网站| 一区二区三区高清视频在线| 欧美日本亚洲视频在线播放| 国产精品免费一区二区三区在线| 成年女人永久免费观看视频| 女人十人毛片免费观看3o分钟| 国产精品嫩草影院av在线观看 | 中文字幕人妻熟人妻熟丝袜美| 尾随美女入室| 成人高潮视频无遮挡免费网站| 色哟哟哟哟哟哟| 精华霜和精华液先用哪个| 精品一区二区免费观看| 久久久久久久午夜电影| 精品久久国产蜜桃| 嫩草影院新地址| av女优亚洲男人天堂| 精品99又大又爽又粗少妇毛片 | 精品午夜福利在线看| 丰满人妻一区二区三区视频av| 国产一区二区三区在线臀色熟女| 精品久久久久久久末码| 色精品久久人妻99蜜桃| 国产精品伦人一区二区| 日韩欧美国产在线观看| 成人美女网站在线观看视频| 中文资源天堂在线| 亚洲人成网站在线播放欧美日韩| 亚洲av中文av极速乱 | 在现免费观看毛片| 中出人妻视频一区二区| 高清在线国产一区| 欧美激情国产日韩精品一区| 小说图片视频综合网站| 午夜免费成人在线视频| 日本精品一区二区三区蜜桃| 男女边吃奶边做爰视频| 成人午夜高清在线视频| 久99久视频精品免费| 一卡2卡三卡四卡精品乱码亚洲| 国产精品野战在线观看| 人妻制服诱惑在线中文字幕| 成年女人看的毛片在线观看| 97人妻精品一区二区三区麻豆| 亚洲熟妇中文字幕五十中出| 亚洲精品粉嫩美女一区| 免费观看人在逋| 18禁裸乳无遮挡免费网站照片| 如何舔出高潮| 亚洲在线自拍视频| 久久精品国产亚洲av涩爱 | 乱人视频在线观看| 日韩av在线大香蕉| 国产人妻一区二区三区在| 中文字幕精品亚洲无线码一区| 欧美日韩亚洲国产一区二区在线观看| 九九爱精品视频在线观看| 联通29元200g的流量卡| 在线观看66精品国产| av在线观看视频网站免费| 免费高清视频大片| 99久久无色码亚洲精品果冻| 热99re8久久精品国产| 老女人水多毛片| 搞女人的毛片| 国产一级毛片七仙女欲春2| 国产一区二区在线观看日韩| 噜噜噜噜噜久久久久久91| 免费观看的影片在线观看| 国产v大片淫在线免费观看| 黄色欧美视频在线观看| 日本黄色视频三级网站网址| 无遮挡黄片免费观看| 久久久久久伊人网av| 成年女人永久免费观看视频| 99久久精品国产国产毛片| 成年女人毛片免费观看观看9| 大又大粗又爽又黄少妇毛片口| 国产午夜精品论理片| 国产女主播在线喷水免费视频网站 | 午夜福利在线观看免费完整高清在 | 精品久久久久久久久av| 搡老妇女老女人老熟妇| 午夜激情欧美在线| 真实男女啪啪啪动态图| 99久久九九国产精品国产免费| 最新中文字幕久久久久| 最近中文字幕高清免费大全6 | 日日夜夜操网爽| 国产午夜福利久久久久久| 亚洲熟妇中文字幕五十中出| 成人午夜高清在线视频| 免费人成视频x8x8入口观看| 亚洲国产色片| 最近视频中文字幕2019在线8| 99久久中文字幕三级久久日本| 精品久久久久久久久av| 日韩av在线大香蕉| 亚洲成人久久爱视频| 精品午夜福利视频在线观看一区| 精品久久久噜噜| 国产白丝娇喘喷水9色精品| 在线观看av片永久免费下载| 亚洲精品日韩av片在线观看| 色哟哟·www| 欧美丝袜亚洲另类 | 搡老岳熟女国产| 久久婷婷人人爽人人干人人爱| 久久久精品欧美日韩精品| 日本 欧美在线| 在现免费观看毛片| 精品国内亚洲2022精品成人| 久久欧美精品欧美久久欧美| 亚洲av免费在线观看| 97热精品久久久久久| 成人性生交大片免费视频hd| 国产成人a区在线观看| 国产精品久久久久久久电影| 3wmmmm亚洲av在线观看| a级毛片免费高清观看在线播放| 国产精品美女特级片免费视频播放器| 亚洲午夜理论影院| 成人永久免费在线观看视频| 国产高清视频在线播放一区| 最好的美女福利视频网| 九九久久精品国产亚洲av麻豆| 成人无遮挡网站| 久久精品国产亚洲av天美| 免费看光身美女| 3wmmmm亚洲av在线观看| 中国美女看黄片| 亚洲三级黄色毛片| 三级国产精品欧美在线观看| 嫩草影院精品99| 亚洲男人的天堂狠狠| 亚洲色图av天堂| 波多野结衣高清作品| 能在线免费观看的黄片| 欧美高清成人免费视频www| 一级黄色大片毛片| 亚洲精品日韩av片在线观看| 丰满乱子伦码专区| 男人和女人高潮做爰伦理| 久久人妻av系列| 亚洲av二区三区四区| 人人妻人人澡欧美一区二区| 天堂√8在线中文| 毛片一级片免费看久久久久 | 日日摸夜夜添夜夜添av毛片 | 亚洲国产精品sss在线观看| 两个人视频免费观看高清| 成人二区视频| 伦精品一区二区三区| 最好的美女福利视频网| 一边摸一边抽搐一进一小说| 日本三级黄在线观看| 免费av不卡在线播放| 精品一区二区免费观看| 欧美成人一区二区免费高清观看| 91麻豆av在线| 成人高潮视频无遮挡免费网站| 国产 一区精品| 欧美日韩乱码在线| xxxwww97欧美| 韩国av一区二区三区四区| 国产男靠女视频免费网站| 亚洲狠狠婷婷综合久久图片| 亚洲成人久久爱视频| 麻豆av噜噜一区二区三区| 亚洲在线自拍视频| 动漫黄色视频在线观看| 免费av毛片视频| 免费高清视频大片| 国产精品久久久久久av不卡| 51国产日韩欧美| 精品日产1卡2卡| 亚洲无线观看免费| 国产精品无大码| 我要看日韩黄色一级片| 久久久国产成人精品二区| 中国美白少妇内射xxxbb| eeuss影院久久| 亚洲av第一区精品v没综合| 成人亚洲精品av一区二区| 欧美区成人在线视频| 99在线人妻在线中文字幕| 亚洲成a人片在线一区二区| 色综合色国产| 99热只有精品国产| 精品无人区乱码1区二区| 黄色配什么色好看| 精品午夜福利在线看| 中文字幕精品亚洲无线码一区| 亚洲欧美日韩东京热| 国产精品,欧美在线| 简卡轻食公司| 日韩,欧美,国产一区二区三区 | 亚洲第一区二区三区不卡| 国产精品精品国产色婷婷| 午夜精品一区二区三区免费看| 中文亚洲av片在线观看爽| 精品午夜福利在线看| 国产美女午夜福利| 免费观看在线日韩| 亚洲精品一卡2卡三卡4卡5卡| 国产真实伦视频高清在线观看 | 久久精品国产鲁丝片午夜精品 | 一区二区三区免费毛片| 老女人水多毛片| 国产免费一级a男人的天堂| 精品久久国产蜜桃| 成人性生交大片免费视频hd| 午夜日韩欧美国产| 成人av一区二区三区在线看| 亚洲无线在线观看| a在线观看视频网站| 国产精品电影一区二区三区| 亚洲 国产 在线| 免费看美女性在线毛片视频| 亚洲欧美精品综合久久99| 天天一区二区日本电影三级| 午夜a级毛片| 午夜精品在线福利| 啦啦啦韩国在线观看视频| 久久热精品热| 亚洲av日韩精品久久久久久密| 国产成人av教育| 中文资源天堂在线| 美女大奶头视频| 欧美在线一区亚洲| 99久久成人亚洲精品观看| 99视频精品全部免费 在线| 久久久久久久久久黄片| 国产人妻一区二区三区在| 十八禁国产超污无遮挡网站| 身体一侧抽搐| 亚洲国产欧美人成| 成熟少妇高潮喷水视频| 国产午夜福利久久久久久| 久久久久久久午夜电影| 老司机午夜福利在线观看视频| 婷婷精品国产亚洲av| 亚洲内射少妇av| 九九爱精品视频在线观看| 97碰自拍视频| 一个人观看的视频www高清免费观看| 国产精品一区二区免费欧美| 午夜免费成人在线视频| 免费人成视频x8x8入口观看| 男女那种视频在线观看| 又粗又爽又猛毛片免费看| 搡老妇女老女人老熟妇| 国内久久婷婷六月综合欲色啪| 日本-黄色视频高清免费观看| 男女边吃奶边做爰视频| 91精品国产九色| 亚洲一区二区三区色噜噜| 我要搜黄色片| 国产探花在线观看一区二区| 精华霜和精华液先用哪个| 久99久视频精品免费| 亚洲专区国产一区二区| 欧美性猛交╳xxx乱大交人| 少妇人妻精品综合一区二区 | 偷拍熟女少妇极品色| 别揉我奶头~嗯~啊~动态视频| 在线播放无遮挡| 久久精品夜夜夜夜夜久久蜜豆| 国产精华一区二区三区| 午夜日韩欧美国产| 亚洲av免费在线观看| 日本在线视频免费播放| 国产乱人视频| 高清毛片免费观看视频网站| 亚洲性久久影院| 亚洲自偷自拍三级| 亚洲av不卡在线观看| 久久欧美精品欧美久久欧美| 嫁个100分男人电影在线观看| 亚洲最大成人中文| 亚洲熟妇中文字幕五十中出| 99久久成人亚洲精品观看| 亚洲真实伦在线观看| 日本a在线网址| 国产激情偷乱视频一区二区| 丰满的人妻完整版| 精品无人区乱码1区二区| 国产黄片美女视频| 亚洲狠狠婷婷综合久久图片| 99热这里只有是精品50| 久久亚洲精品不卡| 欧美3d第一页| 欧美日本亚洲视频在线播放| 国内精品久久久久久久电影| av.在线天堂| 亚洲精品影视一区二区三区av| 久久久久久久久久成人| 少妇裸体淫交视频免费看高清| 午夜视频国产福利| 欧美精品国产亚洲| 成人欧美大片| 欧美黑人欧美精品刺激| 亚洲第一电影网av| 3wmmmm亚洲av在线观看| 亚洲最大成人手机在线| 在线观看舔阴道视频| 久久精品国产亚洲av香蕉五月| 一级黄片播放器| 全区人妻精品视频| 久久午夜福利片| 欧美激情久久久久久爽电影| 亚洲狠狠婷婷综合久久图片| 天堂影院成人在线观看| 久久国产乱子免费精品| 国产精品美女特级片免费视频播放器| 麻豆成人av在线观看| 日日干狠狠操夜夜爽| 亚洲午夜理论影院| 国产精品免费一区二区三区在线| 老司机午夜福利在线观看视频| 国产乱人视频| 久久久久国内视频| 老熟妇乱子伦视频在线观看| 丰满乱子伦码专区| 精品久久久久久久久久久久久| 97碰自拍视频| 18禁裸乳无遮挡免费网站照片| 亚洲无线观看免费| 国产精品人妻久久久影院| 国产精品久久久久久久久免| 成人av一区二区三区在线看| 夜夜爽天天搞| 波多野结衣巨乳人妻| 免费av毛片视频| 女人被狂操c到高潮| 身体一侧抽搐| 夜夜爽天天搞| 国产午夜精品论理片| 国产精品人妻久久久影院| 午夜精品久久久久久毛片777| 国产精品亚洲美女久久久| 亚州av有码| 亚洲av电影不卡..在线观看| 成人av一区二区三区在线看| 亚洲五月天丁香| 欧美+日韩+精品| 少妇裸体淫交视频免费看高清| 欧美+日韩+精品| 九色成人免费人妻av| 天堂网av新在线| 一级a爱片免费观看的视频| 欧美丝袜亚洲另类 | 日本免费a在线| 中文字幕高清在线视频| 美女高潮喷水抽搐中文字幕| 天堂影院成人在线观看| 成人国产麻豆网| 国产v大片淫在线免费观看| av国产免费在线观看| 成人亚洲精品av一区二区| 国产一区二区在线观看日韩| 麻豆av噜噜一区二区三区| 久久精品人妻少妇| 非洲黑人性xxxx精品又粗又长| 国产高清视频在线播放一区| 在现免费观看毛片| 国产亚洲精品久久久久久毛片| 亚洲黑人精品在线| 亚洲,欧美,日韩| 国产精品自产拍在线观看55亚洲| 在线国产一区二区在线| 亚洲精品粉嫩美女一区| 91麻豆av在线| 午夜福利在线在线| 欧美色欧美亚洲另类二区| 国产三级在线视频| 很黄的视频免费| 自拍偷自拍亚洲精品老妇| 麻豆精品久久久久久蜜桃| 无遮挡黄片免费观看| 免费看日本二区| 91av网一区二区| 如何舔出高潮| 午夜影院日韩av| 成年女人毛片免费观看观看9| 男人舔奶头视频| 午夜精品在线福利| 成人av一区二区三区在线看| 无人区码免费观看不卡| 国产av在哪里看| 成人二区视频| 老女人水多毛片| 亚洲综合色惰| 日日摸夜夜添夜夜添小说| 日日干狠狠操夜夜爽| 久久精品91蜜桃| 国产麻豆成人av免费视频| 亚洲国产精品sss在线观看| 久久精品国产鲁丝片午夜精品 | 久久精品国产亚洲网站| av中文乱码字幕在线| 九九爱精品视频在线观看| 亚洲av成人av| 免费看光身美女| 久久午夜亚洲精品久久| 久久精品人妻少妇| 亚洲久久久久久中文字幕| 精品乱码久久久久久99久播| 国产欧美日韩精品一区二区| 免费观看人在逋| 变态另类成人亚洲欧美熟女| 高清毛片免费观看视频网站| 久久久久久九九精品二区国产| 欧美成人免费av一区二区三区| 国产精品人妻久久久影院| 欧美日韩国产亚洲二区| 亚洲va在线va天堂va国产| 波多野结衣高清无吗| 亚洲美女搞黄在线观看 | 久久久久久大精品| 色综合色国产| 亚洲国产色片| 国产午夜精品久久久久久一区二区三区 | 国产精品一区二区三区四区久久| 国产亚洲精品综合一区在线观看| 亚洲精品影视一区二区三区av| 国产精品综合久久久久久久免费| 性欧美人与动物交配| 国产精品国产三级国产av玫瑰|