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

    云南省測繪成果向2000國家大地坐標(biāo)系轉(zhuǎn)換的方法

    2016-12-20 09:59:07周峻松李石華李文華彭雙云
    測繪通報 2016年11期
    關(guān)鍵詞:柵格殘差西安

    周峻松,李石華,,李文華,彭雙云

    (1. 云南省基礎(chǔ)地理信息中心,云南 昆明 650034; 2. 云南師范大學(xué)旅游與地理科學(xué)學(xué)院,云南 昆明 650500)

    ?

    云南省測繪成果向2000國家大地坐標(biāo)系轉(zhuǎn)換的方法

    周峻松1,李石華1,2,李文華1,彭雙云2

    (1. 云南省基礎(chǔ)地理信息中心,云南 昆明 650034; 2. 云南師范大學(xué)旅游與地理科學(xué)學(xué)院,云南 昆明 650500)

    以國家測繪地理信息局對2000國家大地坐標(biāo)系(CGCS2000)推廣使用的要求為指導(dǎo),根據(jù)坐標(biāo)系建立的原理,結(jié)合云南省現(xiàn)有測繪成果實際,利用布爾莎七參數(shù)模型和柵格模型法,開展了云南省省級基礎(chǔ)測繪成果由1954北京坐標(biāo)系、1980西安坐標(biāo)系向2000國家大地坐標(biāo)系轉(zhuǎn)換方法研究。結(jié)果表明:①兩種方法的殘差中誤差均在其反轉(zhuǎn)時的殘余誤差范圍內(nèi),證明了七參數(shù)法和柵格模型法是可靠的;②從對兩種方法的檢驗結(jié)果來看,模型法比柵格法的改算點位殘差?。虎蹡鸥衲P透恼捎糜诟呔纫蟮淖鴺?biāo)改算,而七參數(shù)法可用于對精度要求不高或需要保持舊成果整體精度的坐標(biāo)轉(zhuǎn)換。

    1954北京坐標(biāo)系;1980西安坐標(biāo)系;2000國家大地坐標(biāo)系;坐標(biāo)轉(zhuǎn)換;布爾莎模型;柵格模型

    此前,關(guān)于啟用2000國家大地坐標(biāo)系(CGCS2000)的提議已獲國務(wù)院批準(zhǔn),國家測繪地理信息局規(guī)定:自2008年7月1日啟用CGCS2000,并用8~10年的時間完成現(xiàn)行國家大地坐標(biāo)系向CGCS2000的過渡和轉(zhuǎn)換。國家測繪地理信息局提供了CGCS2000的技術(shù)參數(shù),并對新舊坐標(biāo)系的轉(zhuǎn)換和使用作出了說明。由于1980西安坐標(biāo)系已采用20多年,大量的測繪成果都采用的是1980西安坐標(biāo)系,甚至是1954北京坐標(biāo)系[1],因此面臨著大量的坐標(biāo)轉(zhuǎn)換問題。眾多學(xué)者就其轉(zhuǎn)換方法開展了深入研究與試驗[2-5],其研究思路與方法可為后續(xù)研究提供借鑒與參考。

    根據(jù)國家測繪地理信息局《關(guān)于加快2000國家大地坐標(biāo)系推廣使用的通知》(國測國發(fā)〔2013〕11號)的要求,尚未開展此項工作的省份,應(yīng)盡快完成各自省份地理信息成果的梳理、省級成果轉(zhuǎn)換,推動市縣級成果轉(zhuǎn)換。為此,云南省測繪地理信息局將該項工作作為2014年度基礎(chǔ)測繪任務(wù)之一,由云南省基礎(chǔ)地理信息中心承擔(dān)并完成“云南省基礎(chǔ)測繪成果2000國家大地坐標(biāo)系轉(zhuǎn)換”任務(wù)。本文分析了現(xiàn)有坐標(biāo)系轉(zhuǎn)換方法,在調(diào)研其他省份的工作經(jīng)驗與充分分析相關(guān)文獻研究結(jié)論的基礎(chǔ)上,結(jié)合本省基礎(chǔ)測繪成果和地域?qū)嶋H,有必要就省級測繪成果坐標(biāo)系轉(zhuǎn)換方法展開研究,以期為1954北京坐標(biāo)系、1980西安坐標(biāo)系向CGCS2000轉(zhuǎn)換方法提供案例。

    一、資料利用與轉(zhuǎn)換方法

    1. 資料利用

    本研究中充分利用了國家、省級、地方現(xiàn)有的測繪成果。具體包括:

    1) 云南省內(nèi)的天文大地網(wǎng)與高精度GPS2000網(wǎng)聯(lián)合平差地面網(wǎng)點CGCS2000成果2168個點,但經(jīng)過與現(xiàn)有的舊成果進行人工比對、檢核和甄別,預(yù)選2100個點作為公共已知點使用。

    2) 為改善已知點的圖形分布結(jié)構(gòu)、提高轉(zhuǎn)換精度,從西安大地數(shù)據(jù)處理中心補充收集到省外周邊428個點CGCS2000坐標(biāo)系成果、省內(nèi)68個點1980西安坐標(biāo)系與1954北京坐標(biāo)系的成果。經(jīng)過與現(xiàn)有的成果比對,將343個點(省內(nèi)68個點、省外275個點)作為預(yù)選公共已知點使用。

    3) 似大地水準(zhǔn)面精化成果(CGCS2000坐標(biāo)系),獲取高程異常值,將正常高換算為CGCS2000坐標(biāo)系大地高。

    4) 利用《1980西安大地坐標(biāo)系統(tǒng)成果圖集》,采集其中的1954北京坐標(biāo)系、1980西安坐標(biāo)系高程異常圖等值線,用于高程異常值空間分布的內(nèi)插,獲取高程異常值。

    5) 國家測繪地理信息局提供的“1980西安坐標(biāo)系到CGCS2000坐標(biāo)系空間格網(wǎng)改正量”,主要用于結(jié)果比較及檢驗。

    6) 云南省1980西安坐標(biāo)系、2000國家大地坐標(biāo)系兩套成果的GPS B級點31個點、GPS C級點1311個點,作為結(jié)果的檢驗。

    7) 以似大地水準(zhǔn)面精化工作時聯(lián)測轉(zhuǎn)換的成果40個點作為結(jié)果檢驗。

    8) 利用少量國家Ⅱ等點的控制成果。聯(lián)測了安寧、玉溪交接區(qū)域的高魯山Ⅱ、雞背山Ⅱ、大黑山(二)Ⅱ,對送鬼山ⅢⅡ進行了聯(lián)測檢驗工作。

    2. 轉(zhuǎn)換方法

    利用布爾莎模型系統(tǒng)轉(zhuǎn)換后,在已知點上還有殘留誤差,且各已知點的模型誤差表現(xiàn)出明顯的區(qū)域差異性。而這種區(qū)域差異性期望用一種二維空間場模型來描述,常用的場模型有6種(如圖1所示)。

    圖1 二維空間場模型

    這些場模型的核心是利用已知點進行空間插值??臻g內(nèi)插可分為幾何方法、統(tǒng)計方法、空間統(tǒng)計方法、函數(shù)方法、隨機模擬方法、物理模型模擬方法和綜合方法[6]。其中幾何方法主要是基于“地理學(xué)第一定律”的基本假設(shè)[7],即鄰近的區(qū)域比距離遠的區(qū)域更相似。最常用的幾何方法有泰森多邊形(最近距離法)和反距離加權(quán)方法。反距離加權(quán)法是最常用的空間內(nèi)插方法之一,認(rèn)為與未采樣點距離最近的若干個點對未采樣點值的貢獻最大,其貢獻與距離成反比可表示為

    (1)

    式中,Z為待插值的柵格點;Zi為第i個已知點的改正值;di為擬合點到第i個已知點的大圓距離;n為用于擬合插值的已知點的數(shù)目。

    用該方法計算出每一個點空間分布的模型殘差改正值,最后形成兩個二維分量上的改正量柵格集,在本研究中將其定義為柵格改正法。

    移動擬合法的計算模型中,其定權(quán)用的是空間直線距離d,而在實際操作時用的是地理坐標(biāo),其計算的距離與空間直線距離有一定差異。但由于其作用范圍只是附近的點,比較地理坐標(biāo)距離計算的權(quán)值與空間直線距離計算的權(quán)值發(fā)現(xiàn),其權(quán)值的比重基本一致,而計算工作量得到大大簡化,故用地理坐標(biāo)代替計算。

    二、計算過程與結(jié)果分析

    1. 布爾莎模型七參數(shù)解算

    在求解七參數(shù)時,先將大地坐標(biāo)轉(zhuǎn)換到空間直角坐標(biāo)系下,再求解七參數(shù),流程如圖2所示。

    圖2 計算流程

    (1) 大地坐標(biāo)與空間直角坐標(biāo)系之間的轉(zhuǎn)換

    大地坐標(biāo)(B、L、H)轉(zhuǎn)換為空間直角坐標(biāo)(X、Y、Z)的公式為

    (2)

    空間直角坐標(biāo)(X、Y、Z)轉(zhuǎn)換為大地坐標(biāo)(B、L、H)的公式為

    (3)

    式(3)用迭代法計算。

    (2) 七參數(shù)轉(zhuǎn)換

    七參數(shù)轉(zhuǎn)換模型采用布爾莎模型,它在兩空間直角坐標(biāo)系之間存在嚴(yán)密的轉(zhuǎn)換關(guān)系,其形式非常簡單,物理意義也非常明晰,即

    (4)

    式中,(X1,Y1,Z1)為原坐標(biāo)系下三維坐標(biāo);(X2,Y2,Z2)為新坐標(biāo)系下三維坐標(biāo);ΔX、ΔY、ΔZ、εX、εY、εZ、m為7個參數(shù)。

    通過這樣的過程改進,不但計算結(jié)果更為精確,7個參數(shù)的物理含義也很明晰。

    (3) 轉(zhuǎn)換步驟

    1) 利用GPS C級點計算的1980西安坐標(biāo)系/1954北京坐標(biāo)系→CGCS2000坐標(biāo)系的參數(shù),將參與1980西安坐標(biāo)系/1954北京坐標(biāo)系天文大地網(wǎng)聯(lián)合平差的地面網(wǎng)點概略改算到2000國家大地坐標(biāo),然后根據(jù)點名、點位(B,L),對照CGCS2000系成果查找重合點。以此為基礎(chǔ),利用新舊成果的正常高信息,檢查正常高是否發(fā)生大的變化,以此來確認(rèn)該標(biāo)石是否發(fā)生變動,進一步甄選重合點。

    2) 將重合點坐標(biāo)換算為大地坐標(biāo)及空間直角坐標(biāo)。

    3) 使用篩選出的全部重合點求取轉(zhuǎn)換參數(shù)(三維轉(zhuǎn)換),再利用該轉(zhuǎn)換參數(shù)回代,求得已知點的坐標(biāo),計算已知點的殘差。

    4) 分析殘差值大于2~3倍殘差中誤差的點的可靠性,檢查點位和成果的正確性,重新確定重合點。

    5) 再重復(fù)步驟3)和4),直到回代殘差到一定范圍內(nèi)為止,最終確定參與求解參數(shù)的已知點。

    根據(jù)上述步驟,可計算出1980西安坐標(biāo)系/1954北京坐標(biāo)系→CGCS2000坐標(biāo)系的參數(shù)(見表1)。

    表1 七參數(shù)及中誤差統(tǒng)計

    2. 柵格法改正模型

    (1) 改算步驟

    ①甄選已知點,剔出與周圍互差較大的異常點;②確定柵格改正量的空間格網(wǎng)大??;③根據(jù)已知點分布及空間插值特點,劃定柵格改正量的有效范圍;④計算經(jīng)度B、緯度L兩個分量的柵格改正量;⑤結(jié)果檢驗。

    (2) 改算過程

    1) 根據(jù)已知點改正量,利用TIN計算已知點構(gòu)成的每一個三角區(qū)的變化梯度,得出已知點經(jīng)度、緯度改正量的差分值,在不影響使用精度的前提下確定格網(wǎng)大小。其中,緯度的殘差梯度最大值0.002 45°、經(jīng)度的殘差梯度最大值0.003 72°,對應(yīng)到平面距離上,每1000 m分別對應(yīng)64.9、42.8 mm。

    2) 經(jīng)反復(fù)試算得出,改正量的格網(wǎng)大小為100 m時,在改正量變化最大的區(qū)域,由規(guī)則格網(wǎng)存儲引起的誤差:在X方向<6.49 mm、在Y方向<4.28 mm。因此,本省改正量的格網(wǎng)大小定為100 m左右(實際執(zhí)行為3″)。

    3)雖然收集了云南省周邊一定緩沖區(qū)域的已知點成果,但云南地處邊境,國外部分沒有可用的已知點。并且采用的IDW方法要對已知點的殘差改正量進行空間插值,該方法具有一定外推特性,若外推距離遠,則精度低。結(jié)合現(xiàn)有已知點及本省1∶1萬地形圖生產(chǎn)的需要,確定了柵格改正量的有效值范圍。

    (3) 柵格改正量的計算

    利用布爾莎模型系統(tǒng)轉(zhuǎn)換后,在已知點上還有殘留誤差,可對這些殘差再次進行改正。利用式(2)計算各柵格點的殘差改正量,再利用模型計算系統(tǒng)改正量,兩者疊加(系統(tǒng)+殘差)形成最后成果。

    3. 精度評定

    模型坐標(biāo)轉(zhuǎn)換的精度是通過求取轉(zhuǎn)換參數(shù)的重合點的殘差來體現(xiàn)的。坐標(biāo)轉(zhuǎn)換精度估計依據(jù)以下公式計算

    1) 殘差V=重合點轉(zhuǎn)換坐標(biāo)-重合點已知坐標(biāo)。

    2) 空間坐標(biāo)X與平面坐標(biāo)x的殘差中誤差為

    (5)

    3) 空間坐標(biāo)Y與平面坐標(biāo)y的殘差中誤差為

    (6)

    4) 空間坐標(biāo)Z與大地高H的殘差中誤差為

    (7)

    (8)

    此外,對轉(zhuǎn)換成果還采用了外業(yè)聯(lián)測檢測和與GPSC級點成果進行對比檢驗的方式開展了精度檢核工作。

    (1) 七參數(shù)模型計算精度評價

    此次驗證工作,使用的點共有2443個,其中已知點1222個、檢驗點1221個。為驗證“模型殘差(已知點殘差中誤差)”評定指標(biāo)的可靠性,將篩選出一半的重合點用作已知點計算模型參數(shù),而另一半重合點用作檢驗點進行改算精度檢驗。七參數(shù)模型的誤差見表2。最后,再與全部作為已知點的模型殘差作比較(見表3)。

    表2 七參數(shù)模型誤差統(tǒng)計

    表3 模型殘差與檢驗中誤差的比較 m

    從表3可以看出,1222個點與2443個點的模型殘差與1221個點的檢驗中誤差基本相當(dāng)。再從檢驗中誤差的大小、空間分布(圖略)上看,兩者也基本一致。

    依據(jù)式(6)—式(8)可計算出,1980西安坐標(biāo)系到CGCS2000改算模型殘差:RMS dS=0.283 69 m,dS<0.3 m的占73.33%,最小值為0.012 002 m,最大值為0.88 m,在貢山縣境內(nèi)(麻紀(jì)娃Ⅱ);1954北京坐標(biāo)系到CGCS2000改算的模型殘差:RMS dS=0.718 46 m,dS<0.5 m的占61.2%,最小值為0.011 4 m,最大值為3.124 m,在貢山縣境內(nèi)(卡瓦嘿Ⅱ)。轉(zhuǎn)換模型殘差區(qū)間統(tǒng)計見表4。

    從表4中可以看出,1980西安坐標(biāo)系到CGCS2000改算模型殘差和1954北京坐標(biāo)系到CGCS2000改算的模型殘差差異較大,各點位差值范圍的統(tǒng)計比例差異明顯。1980西安坐標(biāo)系到CGCS2000改算模型殘差的點位誤差差值范圍集中在[0~0.15]、[0.15~0.3],占73.33%;而1954北京坐標(biāo)系到CGCS2000改算模型殘差的點位誤差差值范圍集中在[0.3~0.5]、[0.5~0.883 6],占63.67%。

    表4 轉(zhuǎn)換模型殘差區(qū)間統(tǒng)計

    (2) 柵格模型計算精度評價

    按照常理,改正量在空間分布上不應(yīng)有大的突變;另外,按照計算改正量公式(式(2))的特點,在已知點上應(yīng)有殘差為0的特性。

    ① 圖面檢驗

    對1980西安坐標(biāo)系到CGSC2000的經(jīng)度(dL)、緯度(dB),分別用色彩表示(圖略),色彩過渡均勻,沒有突變現(xiàn)象。

    對1954北京坐標(biāo)系到CGCS2000的經(jīng)度(dL)、緯度(dB),分別用色彩表示(圖略),除西北角色彩稍有不均勻外,其余色彩過渡均勻,沒有突變現(xiàn)象。

    ② 已知點殘差檢驗

    1980西安坐標(biāo)系已知點殘差檢驗,對在有效范圍內(nèi)的2374個已知點,讀取柵格改正量、轉(zhuǎn)換到CGCS2000后,與國家新成果比較,差異均為0,說明1980西安坐標(biāo)系柵格改正量正確。

    1954北京坐標(biāo)系已知點殘差檢驗,對在有效范圍內(nèi)的2369個已知點,讀取柵格改正量、轉(zhuǎn)換為CGSC2000后,與國家新成果比較,差異均為0,說明1954北京坐標(biāo)系柵格改正量正確。

    對1980西安坐標(biāo)系轉(zhuǎn)換為CGCS2000檢驗:柵格法計算2392個已知點的平面位置殘差,結(jié)果均為0;對1954北京坐標(biāo)系轉(zhuǎn)換為CGCS2000檢驗:計算2369個已知點的平面位置殘差,結(jié)果均為0。說明柵格法坐標(biāo)改正量計算正確。

    (3) 其他檢驗

    利用國家測繪地理信息局提供的云南省1311個1980西安坐標(biāo)系、CGCS2000的GPS C級點兩套成果,通過篩選128個重合點,利用布爾莎七參數(shù)坐標(biāo)轉(zhuǎn)換模型解算CGCS2000成果。經(jīng)檢驗,模型法轉(zhuǎn)換平面位置殘差中誤差為0.110 m;柵格法轉(zhuǎn)換平面位置殘差中誤差為0.287 m。

    三、結(jié) 論

    本文根據(jù)云南省基礎(chǔ)測繪成果坐標(biāo)轉(zhuǎn)換的工作需求,基于坐標(biāo)系建立的原理,對其轉(zhuǎn)換方法進行討論與分析,并對結(jié)果進行驗證,結(jié)論與建議如下:

    1) 利用七參數(shù)方法和基于七參數(shù)的柵格模型改正方法對云南省的基礎(chǔ)測繪成果數(shù)據(jù)進行CGCS2000大地坐標(biāo)改算,從概算精度的評價結(jié)果來看,兩種方法是可靠的,均可使用。

    2) 從兩種方法的檢驗結(jié)果來看,模型法比柵格法的改算點位殘差要小,但這不能說明模型法的結(jié)果更好,只能說明模型法與國家測繪地理信息局的轉(zhuǎn)換法方法更接近。柵格改正法可用于高精度要求的坐標(biāo)改算。對精度要求不高或需要保持舊成果整體精度的場合,只需利用七參數(shù)方法進行改正即可。

    3) 在高等級控制點較稀少的地區(qū),可酌情考慮選用三、四等控制點成果作為已知點,以保證圖形結(jié)構(gòu)。

    4) 由于空間坐標(biāo)轉(zhuǎn)換需要用到大地高,而GPS測量可直接獲取大地高,CGCS2000的新成果就直接提供了大地高。但1954北京坐標(biāo)系、1980西安坐標(biāo)系不是真三維空間的坐標(biāo),其提供的正常高成果是通過水準(zhǔn)、三角高程測量得出的,與坐標(biāo)成果是相互獨立的,大地高需要間接推算獲取。王文利等的研究結(jié)果表明[8],重合點的大地高誤差對布爾莎參數(shù)模型轉(zhuǎn)換的平面結(jié)果影響較小,大地高誤差在10 m以內(nèi)時,基本不影響轉(zhuǎn)換的平面結(jié)果。而轉(zhuǎn)換點的大地高誤差對布爾莎七參數(shù)模型轉(zhuǎn)換的平面結(jié)果影響非常小。轉(zhuǎn)換點大地高誤差≤100 m時,可用正常高代替大地高進行轉(zhuǎn)換,且轉(zhuǎn)換結(jié)果不受影響。因此,在本次研究中,1954北京坐標(biāo)系、1980西安坐標(biāo)系的高程異常值可通過《1980西安大地坐標(biāo)系統(tǒng)成果圖集》獲取,且精度不受影響。

    5) DEM、DOM數(shù)據(jù)由于生產(chǎn)單位提交的原始數(shù)據(jù)部分缺少坐標(biāo)信息,沒有頭文件或頭文件格式不對的均轉(zhuǎn)換失敗,須補充頭文件,手工將頭文件里的坐標(biāo)值加改正量平移。

    6) 靠近中央子午線的圖幅,DLG數(shù)據(jù)中的部分要素難以獲取改正量,需重新設(shè)置帶號進行轉(zhuǎn)換。

    [1] 鐘業(yè)勛,童新華,王龍波.從1980西安坐標(biāo)系到2000國家大地坐標(biāo)系的坐標(biāo)變換[J].海洋測繪,2010,30(1):1-3.

    [2] 過家春,趙秀俠,吳艷蘭.空間直角坐標(biāo)與大地坐標(biāo)轉(zhuǎn)換的拉格朗日反演方法[J].測繪學(xué)報,2014,43(10):998-1004.

    [3] 王海棟,柴洪洲,宋國大. 大地坐標(biāo)系換代對海圖的影響分析[J].測繪通報,2009(1):6-9.

    [4] 曾安敏. 基于擬合推估的1980西安坐標(biāo)系到2000國家坐標(biāo)系的變換[J].大地測量與地球動力學(xué),2008,28(5):125-128.

    [5] 呂志平,魏子卿,李軍,等. 我國CGCS2000高精度坐標(biāo)轉(zhuǎn)換格網(wǎng)模型的建立[J]. 測繪學(xué)報,2013,42(6):797-791.

    [6] 李新,程國棟,盧玲.空間內(nèi)插方法比較[J].地球科學(xué)進展,2000,15(3):260-265.

    [7] 李小文,曹春香,常超一.地理學(xué)第一定律與時空鄰近度的提出[J].自然雜志,2007(2):69-71.

    [8] 王文利,郭春喜,程傳錄.大地高誤差對Bursa七參數(shù)平面轉(zhuǎn)換精度的影響[J].測繪科學(xué),2011,36(5):37-38.

    Study on Coordinate System Conversion Methods from BJ54,XA80 to CGCS2000 of Surveying and Mapping Archives in Yunnan Province

    ZHOU Junsong,LI Shihua,LI Wenhua,PENG Shuangyun

    2016-07-19

    國家自然科學(xué)基金(41561086);國家測繪地理信息局地理國情監(jiān)測示范項目(測國土函〔2014〕35號)

    周峻松(1968—),男,高級工程師,主要從事地圖制圖、遙感與地理信息工程和應(yīng)用項目管理等工作。E-mail: 78121467@qq.com

    李石華。E-mail:lsh8010@163.com

    周峻松,李石華,李文華,等.云南省測繪成果向2000國家大地坐標(biāo)系轉(zhuǎn)換的方法[J].測繪通報,2016(11):80-84.

    10.13474/j.cnki.11-2246.2016.0371.

    P258

    B

    0494-0911(2016)11-0080-05

    猜你喜歡
    柵格殘差西安
    基于雙向GRU與殘差擬合的車輛跟馳建模
    西安2021
    Oh 西安
    基于鄰域柵格篩選的點云邊緣點提取方法*
    基于殘差學(xué)習(xí)的自適應(yīng)無人機目標(biāo)跟蹤算法
    基于遞歸殘差網(wǎng)絡(luò)的圖像超分辨率重建
    《西安人的歌》突如其來?
    追根溯源 回到西安
    平穩(wěn)自相關(guān)過程的殘差累積和控制圖
    河南科技(2015年8期)2015-03-11 16:23:52
    不同剖面形狀的柵格壁對柵格翼氣動特性的影響
    亚洲欧美精品综合久久99| 亚洲国产欧美日韩在线播放| 老司机深夜福利视频在线观看| 国产日韩一区二区三区精品不卡| av网站在线播放免费| 欧美一级毛片孕妇| 国产97色在线日韩免费| av福利片在线| 在线观看日韩欧美| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久久久久久久久大奶| 99国产综合亚洲精品| 国内毛片毛片毛片毛片毛片| 香蕉丝袜av| 91av网站免费观看| 国内毛片毛片毛片毛片毛片| 51午夜福利影视在线观看| 亚洲av第一区精品v没综合| 久久中文字幕人妻熟女| 国产有黄有色有爽视频| 中文字幕另类日韩欧美亚洲嫩草| 成人免费观看视频高清| 视频区欧美日本亚洲| 久久久精品欧美日韩精品| 精品国产乱码久久久久久男人| 18禁美女被吸乳视频| 国产不卡一卡二| 男人舔女人的私密视频| 欧美中文综合在线视频| 国产精品乱码一区二三区的特点 | 亚洲国产精品sss在线观看 | 精品国产一区二区久久| 国产蜜桃级精品一区二区三区| 黄色毛片三级朝国网站| 国产亚洲欧美98| 久久精品人人爽人人爽视色| 国产三级在线视频| 啦啦啦免费观看视频1| 丁香欧美五月| 女人高潮潮喷娇喘18禁视频| 91字幕亚洲| 天天躁夜夜躁狠狠躁躁| 自线自在国产av| 久久人人97超碰香蕉20202| 国产欧美日韩一区二区三区在线| 国产无遮挡羞羞视频在线观看| 999精品在线视频| 黄色成人免费大全| 亚洲成人免费电影在线观看| 国产无遮挡羞羞视频在线观看| 黄色 视频免费看| 亚洲精品国产区一区二| 身体一侧抽搐| 亚洲国产欧美日韩在线播放| 国产无遮挡羞羞视频在线观看| 久久久久久免费高清国产稀缺| 欧美日韩av久久| 一级黄色大片毛片| 777久久人妻少妇嫩草av网站| 欧美精品啪啪一区二区三区| 99国产精品一区二区三区| 国产免费男女视频| 国产精品乱码一区二三区的特点 | 91在线观看av| 国产伦一二天堂av在线观看| 热re99久久国产66热| 看黄色毛片网站| 99久久99久久久精品蜜桃| 欧美乱码精品一区二区三区| 欧美日韩国产mv在线观看视频| 国产精品国产高清国产av| 国产精品1区2区在线观看.| 少妇 在线观看| 国产无遮挡羞羞视频在线观看| 日韩免费av在线播放| 婷婷精品国产亚洲av在线| 在线观看www视频免费| 每晚都被弄得嗷嗷叫到高潮| 久久人人97超碰香蕉20202| 色哟哟哟哟哟哟| 18禁观看日本| 国产精品影院久久| 天天影视国产精品| 一区福利在线观看| av视频免费观看在线观看| 99热只有精品国产| 国产成人影院久久av| 香蕉国产在线看| 三级毛片av免费| 精品日产1卡2卡| 熟女少妇亚洲综合色aaa.| 91老司机精品| 精品高清国产在线一区| 国产免费男女视频| 亚洲专区国产一区二区| 老司机在亚洲福利影院| 精品一区二区三区视频在线观看免费 | 男人操女人黄网站| 久99久视频精品免费| 日本撒尿小便嘘嘘汇集6| 成人亚洲精品av一区二区 | 丁香欧美五月| 岛国视频午夜一区免费看| 99在线人妻在线中文字幕| 国产一区二区激情短视频| 午夜亚洲福利在线播放| 亚洲精华国产精华精| 99热国产这里只有精品6| 老司机午夜福利在线观看视频| 9色porny在线观看| 亚洲午夜理论影院| 少妇 在线观看| ponron亚洲| 国产又爽黄色视频| 久久精品影院6| 性色av乱码一区二区三区2| 中文欧美无线码| 女人精品久久久久毛片| 亚洲精品国产色婷婷电影| 99久久99久久久精品蜜桃| 男男h啪啪无遮挡| 成人精品一区二区免费| 久久精品亚洲精品国产色婷小说| 在线观看66精品国产| 久久人人97超碰香蕉20202| 国产主播在线观看一区二区| 日韩精品免费视频一区二区三区| 欧美黑人欧美精品刺激| 俄罗斯特黄特色一大片| 69av精品久久久久久| 麻豆国产av国片精品| 久久精品国产99精品国产亚洲性色 | 嫩草影视91久久| 一a级毛片在线观看| 免费高清视频大片| 国产有黄有色有爽视频| a在线观看视频网站| 性少妇av在线| 91大片在线观看| 久久久久久久久中文| 国产精品98久久久久久宅男小说| 国产欧美日韩一区二区三区在线| 国产真人三级小视频在线观看| 欧美激情 高清一区二区三区| 成在线人永久免费视频| 日韩人妻精品一区2区三区| bbb黄色大片| 精品午夜福利视频在线观看一区| 欧美久久黑人一区二区| 无限看片的www在线观看| 亚洲精品中文字幕在线视频| 无遮挡黄片免费观看| 亚洲中文av在线| 男人的好看免费观看在线视频 | 精品欧美一区二区三区在线| xxxhd国产人妻xxx| 激情在线观看视频在线高清| 天天躁夜夜躁狠狠躁躁| 久久精品国产清高在天天线| 啦啦啦 在线观看视频| 亚洲国产精品一区二区三区在线| 国产欧美日韩综合在线一区二区| 国产精品一区二区免费欧美| 欧美日韩瑟瑟在线播放| 国产欧美日韩一区二区三区在线| 久久天堂一区二区三区四区| 黄频高清免费视频| 欧美日韩视频精品一区| 高清av免费在线| 久久精品国产综合久久久| 欧美亚洲日本最大视频资源| 免费少妇av软件| 久久精品人人爽人人爽视色| 久久性视频一级片| 亚洲第一欧美日韩一区二区三区| 9色porny在线观看| 一本大道久久a久久精品| 老熟妇仑乱视频hdxx| 欧美黑人精品巨大| 欧美激情高清一区二区三区| 性少妇av在线| 日日摸夜夜添夜夜添小说| 女同久久另类99精品国产91| 无遮挡黄片免费观看| 亚洲中文av在线| 国产成人影院久久av| 亚洲第一欧美日韩一区二区三区| 成人精品一区二区免费| 一进一出抽搐gif免费好疼 | 大码成人一级视频| 咕卡用的链子| 午夜视频精品福利| 另类亚洲欧美激情| 在线观看一区二区三区| 久久久国产精品麻豆| 成年人黄色毛片网站| 国产亚洲欧美精品永久| 美女福利国产在线| 欧美日韩黄片免| av视频免费观看在线观看| 高清黄色对白视频在线免费看| 日本wwww免费看| 国产av一区在线观看免费| 精品国内亚洲2022精品成人| 久久国产精品影院| 久久狼人影院| 亚洲欧美一区二区三区黑人| 无人区码免费观看不卡| 夜夜躁狠狠躁天天躁| 又紧又爽又黄一区二区| 免费看a级黄色片| 91成年电影在线观看| 99香蕉大伊视频| 久久久久久久久免费视频了| 亚洲成人免费电影在线观看| 亚洲精品久久成人aⅴ小说| 亚洲av第一区精品v没综合| 亚洲免费av在线视频| 久久国产亚洲av麻豆专区| 在线观看www视频免费| 日韩欧美一区二区三区在线观看| 免费观看人在逋| 久久这里只有精品19| 精品国内亚洲2022精品成人| 精品国产超薄肉色丝袜足j| 国产免费av片在线观看野外av| 日韩人妻精品一区2区三区| 国产亚洲欧美精品永久| 757午夜福利合集在线观看| 免费看a级黄色片| 日本免费a在线| 香蕉国产在线看| 欧美另类亚洲清纯唯美| 不卡av一区二区三区| 波多野结衣av一区二区av| 欧美日韩亚洲综合一区二区三区_| 国产成人av教育| 亚洲 欧美一区二区三区| 亚洲人成网站在线播放欧美日韩| 丰满的人妻完整版| 色精品久久人妻99蜜桃| www国产在线视频色| 午夜福利影视在线免费观看| 女人高潮潮喷娇喘18禁视频| 一个人免费在线观看的高清视频| 日本wwww免费看| 91九色精品人成在线观看| 亚洲伊人色综图| 90打野战视频偷拍视频| 国产精品秋霞免费鲁丝片| 啦啦啦在线免费观看视频4| 老熟妇仑乱视频hdxx| 天堂√8在线中文| 麻豆国产av国片精品| 免费一级毛片在线播放高清视频 | 国产成人欧美| 亚洲熟妇熟女久久| 99久久国产精品久久久| 国产99久久九九免费精品| 欧美日韩黄片免| 五月开心婷婷网| 久久久国产成人免费| 中文字幕av电影在线播放| 午夜精品国产一区二区电影| 久久久国产精品麻豆| 十分钟在线观看高清视频www| 无限看片的www在线观看| 国产xxxxx性猛交| 久久亚洲真实| 最近最新中文字幕大全电影3 | 男女下面插进去视频免费观看| 制服人妻中文乱码| 精品国产超薄肉色丝袜足j| 国产精品日韩av在线免费观看 | 一进一出抽搐动态| 国产激情欧美一区二区| 国产精品影院久久| 国内毛片毛片毛片毛片毛片| 91麻豆av在线| 国产乱人伦免费视频| 80岁老熟妇乱子伦牲交| 国产精品久久久久久人妻精品电影| 高清av免费在线| 黄色成人免费大全| 深夜精品福利| 视频区图区小说| 无限看片的www在线观看| 日本 av在线| 久久久久九九精品影院| 午夜免费鲁丝| 日本免费一区二区三区高清不卡 | 精品免费久久久久久久清纯| 欧美日韩亚洲国产一区二区在线观看| 久久国产精品影院| 亚洲精品国产精品久久久不卡| 热99re8久久精品国产| 色播在线永久视频| 国产三级在线视频| 久久99一区二区三区| xxx96com| 婷婷精品国产亚洲av在线| 日韩欧美免费精品| 欧洲精品卡2卡3卡4卡5卡区| 日本黄色视频三级网站网址| 啦啦啦 在线观看视频| 在线av久久热| av免费在线观看网站| 欧美黄色片欧美黄色片| bbb黄色大片| 欧美国产精品va在线观看不卡| 国产主播在线观看一区二区| 久久精品亚洲av国产电影网| 亚洲精品av麻豆狂野| 国产野战对白在线观看| 久久亚洲真实| 精品久久久久久久久久免费视频 | 免费少妇av软件| 高清黄色对白视频在线免费看| 久久午夜综合久久蜜桃| 欧美亚洲日本最大视频资源| 国产一区二区三区在线臀色熟女 | 丰满的人妻完整版| 日韩有码中文字幕| 午夜福利影视在线免费观看| 久热这里只有精品99| 国产极品粉嫩免费观看在线| 亚洲自偷自拍图片 自拍| 亚洲成人国产一区在线观看| 国产成人av教育| 色老头精品视频在线观看| 国产精品一区二区三区四区久久 | 日韩欧美国产一区二区入口| 精品福利永久在线观看| av国产精品久久久久影院| 一a级毛片在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品乱码一区二三区的特点 | 91av网站免费观看| av天堂久久9| 精品国产一区二区三区四区第35| 国产亚洲欧美98| 99国产极品粉嫩在线观看| 欧美日韩一级在线毛片| 成人国产一区最新在线观看| 操出白浆在线播放| 国产欧美日韩一区二区三区在线| 国产成人av教育| 久久久久久人人人人人| 国产成人一区二区三区免费视频网站| 亚洲伊人色综图| 亚洲人成伊人成综合网2020| 日本五十路高清| 国产极品粉嫩免费观看在线| 亚洲国产欧美日韩在线播放| 在线观看日韩欧美| 亚洲av成人不卡在线观看播放网| 精品无人区乱码1区二区| 18禁观看日本| 免费在线观看亚洲国产| 无遮挡黄片免费观看| 日本撒尿小便嘘嘘汇集6| 亚洲精品在线观看二区| 国产成人一区二区三区免费视频网站| 老司机靠b影院| 午夜精品在线福利| 女人被躁到高潮嗷嗷叫费观| 亚洲人成电影免费在线| 纯流量卡能插随身wifi吗| 桃红色精品国产亚洲av| 女同久久另类99精品国产91| 香蕉国产在线看| 999精品在线视频| 波多野结衣一区麻豆| 老汉色∧v一级毛片| 国产色视频综合| 欧美精品啪啪一区二区三区| 国产精品爽爽va在线观看网站 | 午夜福利免费观看在线| 看片在线看免费视频| 19禁男女啪啪无遮挡网站| 啦啦啦在线免费观看视频4| 国产亚洲欧美在线一区二区| 黑人操中国人逼视频| 免费高清在线观看日韩| 亚洲自偷自拍图片 自拍| 亚洲精品中文字幕一二三四区| 欧美精品啪啪一区二区三区| 久久天躁狠狠躁夜夜2o2o| 国产精品99久久99久久久不卡| 一级,二级,三级黄色视频| 亚洲欧美精品综合一区二区三区| 桃红色精品国产亚洲av| 欧美日韩亚洲高清精品| 久久这里只有精品19| 国产亚洲精品第一综合不卡| 大码成人一级视频| 少妇被粗大的猛进出69影院| 国产片内射在线| 丰满的人妻完整版| 黄色a级毛片大全视频| 久久精品国产清高在天天线| 日韩大尺度精品在线看网址 | 欧美日韩一级在线毛片| 女人高潮潮喷娇喘18禁视频| av网站在线播放免费| 国产av一区二区精品久久| 伦理电影免费视频| 国产欧美日韩精品亚洲av| 一级毛片女人18水好多| 淫秽高清视频在线观看| 热99国产精品久久久久久7| 免费人成视频x8x8入口观看| 亚洲成人免费电影在线观看| 在线观看免费午夜福利视频| 国产三级在线视频| 欧洲精品卡2卡3卡4卡5卡区| 人人妻人人爽人人添夜夜欢视频| 男男h啪啪无遮挡| 久久中文字幕一级| 色婷婷久久久亚洲欧美| 男女之事视频高清在线观看| 一本大道久久a久久精品| 一夜夜www| 超碰成人久久| 在线看a的网站| 又黄又粗又硬又大视频| a级毛片在线看网站| 亚洲成人免费电影在线观看| 免费在线观看黄色视频的| 乱人伦中国视频| 国产成人欧美在线观看| 80岁老熟妇乱子伦牲交| 免费在线观看日本一区| 久久中文字幕一级| 亚洲精品一卡2卡三卡4卡5卡| 真人一进一出gif抽搐免费| 黑人巨大精品欧美一区二区mp4| 久久九九热精品免费| 成人av一区二区三区在线看| 两人在一起打扑克的视频| 亚洲国产毛片av蜜桃av| 天堂影院成人在线观看| 日韩三级视频一区二区三区| 国产伦一二天堂av在线观看| 中亚洲国语对白在线视频| 精品国产亚洲在线| 午夜老司机福利片| 美女午夜性视频免费| 中出人妻视频一区二区| 国产野战对白在线观看| 国产麻豆69| 国产激情久久老熟女| 777久久人妻少妇嫩草av网站| 大陆偷拍与自拍| 国产精品野战在线观看 | 久久这里只有精品19| 丝袜在线中文字幕| 亚洲一卡2卡3卡4卡5卡精品中文| 免费在线观看影片大全网站| 中文字幕人妻熟女乱码| 欧美亚洲日本最大视频资源| 别揉我奶头~嗯~啊~动态视频| 色哟哟哟哟哟哟| 一区在线观看完整版| 女同久久另类99精品国产91| 日韩av在线大香蕉| 女性被躁到高潮视频| 欧美激情高清一区二区三区| 日日夜夜操网爽| 三上悠亚av全集在线观看| 桃红色精品国产亚洲av| 一边摸一边抽搐一进一出视频| 亚洲九九香蕉| 国产亚洲精品综合一区在线观看 | 午夜久久久在线观看| 亚洲精品一区av在线观看| 99久久国产精品久久久| 日韩免费av在线播放| 久99久视频精品免费| 亚洲成人久久性| 后天国语完整版免费观看| 亚洲熟女毛片儿| 这个男人来自地球电影免费观看| 老司机在亚洲福利影院| 亚洲aⅴ乱码一区二区在线播放 | 99re在线观看精品视频| 性少妇av在线| av天堂在线播放| 性欧美人与动物交配| avwww免费| 免费在线观看日本一区| 亚洲精品中文字幕一二三四区| 中文亚洲av片在线观看爽| 国产精品秋霞免费鲁丝片| 男人舔女人下体高潮全视频| 一级毛片女人18水好多| 亚洲七黄色美女视频| 成人特级黄色片久久久久久久| 亚洲人成伊人成综合网2020| 国产免费av片在线观看野外av| 一级毛片精品| 99国产精品一区二区蜜桃av| 丰满饥渴人妻一区二区三| 色婷婷av一区二区三区视频| 一进一出好大好爽视频| 亚洲av片天天在线观看| 亚洲一区二区三区欧美精品| 国产深夜福利视频在线观看| 夫妻午夜视频| 色综合站精品国产| 97人妻天天添夜夜摸| 在线观看免费午夜福利视频| avwww免费| 啦啦啦在线免费观看视频4| 日日干狠狠操夜夜爽| 亚洲精品中文字幕在线视频| 久久精品国产99精品国产亚洲性色 | 纯流量卡能插随身wifi吗| 91麻豆av在线| 99在线人妻在线中文字幕| 欧美精品啪啪一区二区三区| 99在线视频只有这里精品首页| 亚洲自拍偷在线| 亚洲情色 制服丝袜| 欧美日韩亚洲综合一区二区三区_| 国产精品亚洲一级av第二区| av片东京热男人的天堂| 亚洲成人国产一区在线观看| 一本综合久久免费| 欧美精品一区二区免费开放| 亚洲国产精品sss在线观看 | 国产蜜桃级精品一区二区三区| 美女福利国产在线| 18美女黄网站色大片免费观看| 黑人巨大精品欧美一区二区蜜桃| 在线播放国产精品三级| 午夜福利在线免费观看网站| 麻豆成人av在线观看| aaaaa片日本免费| 国产一卡二卡三卡精品| 成人永久免费在线观看视频| 91国产中文字幕| 亚洲国产精品sss在线观看 | 亚洲精品一二三| 99久久99久久久精品蜜桃| 这个男人来自地球电影免费观看| 欧美日本亚洲视频在线播放| 国产在线精品亚洲第一网站| 久久精品国产亚洲av高清一级| 老司机靠b影院| 校园春色视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 真人做人爱边吃奶动态| 18美女黄网站色大片免费观看| 99国产极品粉嫩在线观看| 91大片在线观看| 少妇裸体淫交视频免费看高清 | 久久精品91无色码中文字幕| 日本五十路高清| 9色porny在线观看| 日日干狠狠操夜夜爽| 三上悠亚av全集在线观看| 校园春色视频在线观看| 淫秽高清视频在线观看| 中文亚洲av片在线观看爽| 国产成人欧美在线观看| 99精品欧美一区二区三区四区| 在线国产一区二区在线| 亚洲人成电影观看| 天天影视国产精品| 黄色毛片三级朝国网站| 国产黄色免费在线视频| 亚洲欧美日韩无卡精品| 日韩大尺度精品在线看网址 | 黑人猛操日本美女一级片| 国产精品自产拍在线观看55亚洲| 亚洲人成网站在线播放欧美日韩| 国产单亲对白刺激| 国产一区二区三区综合在线观看| 国产精品一区二区精品视频观看| 97超级碰碰碰精品色视频在线观看| 精品乱码久久久久久99久播| 欧美精品啪啪一区二区三区| 午夜老司机福利片| 亚洲av美国av| 欧美不卡视频在线免费观看 | 少妇被粗大的猛进出69影院| 香蕉久久夜色| 久9热在线精品视频| 国产av一区二区精品久久| 欧美日韩国产mv在线观看视频| 午夜精品在线福利| 一区二区日韩欧美中文字幕| 精品久久久精品久久久| 久久久国产成人免费| 午夜福利一区二区在线看| 黑人操中国人逼视频| 后天国语完整版免费观看| 91成人精品电影| 啦啦啦在线免费观看视频4| 久久精品成人免费网站| 亚洲欧美日韩无卡精品| 国产1区2区3区精品| 国产黄色免费在线视频| 天天影视国产精品| 男男h啪啪无遮挡| 久久九九热精品免费| 国产主播在线观看一区二区| 亚洲一区二区三区不卡视频| av福利片在线| 亚洲成人国产一区在线观看| 又黄又爽又免费观看的视频| aaaaa片日本免费| 亚洲精品美女久久av网站| 亚洲五月色婷婷综合|