• <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
    不同剖面形狀的柵格壁對柵格翼氣動特性的影響
    亚洲精品乱码久久久久久按摩| 成人无遮挡网站| 亚洲欧美日韩另类电影网站 | 一级av片app| 成年女人在线观看亚洲视频 | 狂野欧美激情性bbbbbb| 精品视频人人做人人爽| 亚洲av欧美aⅴ国产| 最近最新中文字幕免费大全7| 欧美激情在线99| 韩国高清视频一区二区三区| 国产欧美日韩精品一区二区| 人人妻人人爽人人添夜夜欢视频 | 婷婷色综合大香蕉| 亚洲精品乱码久久久久久按摩| 亚洲av免费在线观看| 一区二区三区四区激情视频| av女优亚洲男人天堂| 看十八女毛片水多多多| 亚洲av二区三区四区| 寂寞人妻少妇视频99o| 国产极品天堂在线| 免费播放大片免费观看视频在线观看| 免费不卡的大黄色大毛片视频在线观看| 国国产精品蜜臀av免费| 日本av手机在线免费观看| 国语对白做爰xxxⅹ性视频网站| 人妻少妇偷人精品九色| 极品教师在线视频| 中文字幕人妻熟人妻熟丝袜美| 久久精品熟女亚洲av麻豆精品| 人妻一区二区av| 亚洲国产日韩一区二区| 精品熟女少妇av免费看| 日韩欧美精品v在线| 高清日韩中文字幕在线| 国产一区亚洲一区在线观看| 18禁在线播放成人免费| 亚洲四区av| 免费看av在线观看网站| 成人综合一区亚洲| 九草在线视频观看| 一级毛片黄色毛片免费观看视频| 成年女人看的毛片在线观看| 2022亚洲国产成人精品| 九九久久精品国产亚洲av麻豆| 麻豆国产97在线/欧美| 欧美另类一区| 听说在线观看完整版免费高清| 国产免费一区二区三区四区乱码| 亚洲精品成人av观看孕妇| 欧美极品一区二区三区四区| 欧美高清成人免费视频www| 人体艺术视频欧美日本| 在线观看av片永久免费下载| 国产永久视频网站| 国产免费又黄又爽又色| 天堂网av新在线| 亚洲人成网站在线观看播放| 久久久亚洲精品成人影院| 啦啦啦啦在线视频资源| 久久久久国产精品人妻一区二区| 狂野欧美白嫩少妇大欣赏| 免费黄网站久久成人精品| 国产精品99久久久久久久久| 寂寞人妻少妇视频99o| 日本熟妇午夜| 久久久久久久久久成人| 大陆偷拍与自拍| 2022亚洲国产成人精品| 免费播放大片免费观看视频在线观看| av又黄又爽大尺度在线免费看| 色视频在线一区二区三区| 国产精品熟女久久久久浪| 国产色爽女视频免费观看| 日本黄大片高清| 国产国拍精品亚洲av在线观看| 国产精品嫩草影院av在线观看| 成人美女网站在线观看视频| 国产精品国产三级专区第一集| 久久ye,这里只有精品| 少妇熟女欧美另类| 天天一区二区日本电影三级| 免费电影在线观看免费观看| 日韩 亚洲 欧美在线| 秋霞伦理黄片| 亚洲熟女精品中文字幕| 一区二区三区乱码不卡18| 丝袜喷水一区| 精品人妻偷拍中文字幕| 日本爱情动作片www.在线观看| 中文字幕亚洲精品专区| 国产成人精品久久久久久| 一级片'在线观看视频| 在线免费十八禁| 久久久精品免费免费高清| 国产免费一级a男人的天堂| 成人国产麻豆网| 国产精品国产三级国产av玫瑰| 欧美日韩一区二区视频在线观看视频在线 | 久久久久久国产a免费观看| 亚洲国产av新网站| 黄片wwwwww| 夫妻午夜视频| 能在线免费看毛片的网站| 亚洲精品第二区| 国产在线男女| 22中文网久久字幕| 亚洲天堂av无毛| 我的女老师完整版在线观看| av又黄又爽大尺度在线免费看| 免费大片18禁| 97热精品久久久久久| 青春草国产在线视频| 成人美女网站在线观看视频| 人妻制服诱惑在线中文字幕| 我的老师免费观看完整版| 欧美日韩精品成人综合77777| 日本猛色少妇xxxxx猛交久久| 卡戴珊不雅视频在线播放| 狂野欧美激情性bbbbbb| tube8黄色片| www.av在线官网国产| 国产精品福利在线免费观看| 日韩 亚洲 欧美在线| 可以在线观看毛片的网站| 男人舔奶头视频| 免费大片18禁| av在线老鸭窝| 欧美日韩在线观看h| 男男h啪啪无遮挡| 国产大屁股一区二区在线视频| 国内精品美女久久久久久| 国产亚洲最大av| 超碰97精品在线观看| 日韩成人伦理影院| 亚洲精品乱久久久久久| 91精品国产九色| 欧美区成人在线视频| 国产av不卡久久| 人妻系列 视频| 成人高潮视频无遮挡免费网站| 国产伦精品一区二区三区视频9| 亚洲欧美清纯卡通| 人妻一区二区av| 精华霜和精华液先用哪个| 嫩草影院精品99| 偷拍熟女少妇极品色| 国产爱豆传媒在线观看| 亚洲怡红院男人天堂| 美女xxoo啪啪120秒动态图| 亚洲国产色片| 国产 一区 欧美 日韩| 99久久精品一区二区三区| av女优亚洲男人天堂| 亚洲精品国产av蜜桃| 国产精品一区二区三区四区免费观看| 欧美激情在线99| 噜噜噜噜噜久久久久久91| 高清欧美精品videossex| 高清毛片免费看| 国产白丝娇喘喷水9色精品| eeuss影院久久| 少妇裸体淫交视频免费看高清| 久久久久性生活片| 国产精品麻豆人妻色哟哟久久| 亚洲色图av天堂| 亚洲精品自拍成人| 最近最新中文字幕免费大全7| 91在线精品国自产拍蜜月| 国产黄色视频一区二区在线观看| 丝瓜视频免费看黄片| 久久亚洲国产成人精品v| 国内精品美女久久久久久| 久久久久精品久久久久真实原创| 搡老乐熟女国产| 国产亚洲精品久久久com| 欧美高清性xxxxhd video| 久久久久久久久久成人| 国产乱人视频| 国产精品国产三级国产av玫瑰| 国产成人a∨麻豆精品| 国产精品蜜桃在线观看| 在线播放无遮挡| 1000部很黄的大片| 久久久久久伊人网av| av国产精品久久久久影院| 黄色视频在线播放观看不卡| 日韩电影二区| 日日啪夜夜撸| 青青草视频在线视频观看| 国产视频内射| 久久精品国产亚洲av涩爱| 亚洲欧美精品专区久久| av在线播放精品| 亚洲精品乱码久久久v下载方式| 人妻一区二区av| 亚洲激情五月婷婷啪啪| 大码成人一级视频| 王馨瑶露胸无遮挡在线观看| 在线亚洲精品国产二区图片欧美 | 可以在线观看毛片的网站| 人体艺术视频欧美日本| 久久精品久久久久久噜噜老黄| 舔av片在线| 精品人妻熟女av久视频| 中文字幕久久专区| 国产男女超爽视频在线观看| av免费观看日本| 我要看日韩黄色一级片| 青青草视频在线视频观看| 嫩草影院精品99| 乱系列少妇在线播放| 日韩电影二区| 亚洲精品aⅴ在线观看| 美女国产视频在线观看| 午夜免费鲁丝| 人妻系列 视频| 国产在线男女| 欧美精品人与动牲交sv欧美| 麻豆国产97在线/欧美| 人妻少妇偷人精品九色| 男女无遮挡免费网站观看| 午夜福利在线观看免费完整高清在| 高清日韩中文字幕在线| 69av精品久久久久久| 晚上一个人看的免费电影| 久久鲁丝午夜福利片| 日韩精品有码人妻一区| 国产免费一区二区三区四区乱码| 精品视频人人做人人爽| 蜜臀久久99精品久久宅男| 免费观看无遮挡的男女| 亚洲精品一区蜜桃| 亚洲激情五月婷婷啪啪| 亚洲色图综合在线观看| 亚洲欧洲国产日韩| av国产精品久久久久影院| 国产欧美日韩一区二区三区在线 | 亚洲美女视频黄频| 亚洲精品国产av成人精品| 美女被艹到高潮喷水动态| 成人黄色视频免费在线看| 免费观看无遮挡的男女| 成年av动漫网址| 夫妻性生交免费视频一级片| www.色视频.com| 国模一区二区三区四区视频| 亚洲天堂国产精品一区在线| 免费观看在线日韩| 777米奇影视久久| 亚洲人与动物交配视频| 亚洲欧美成人精品一区二区| 欧美日韩在线观看h| 日日摸夜夜添夜夜爱| videos熟女内射| 身体一侧抽搐| 国产亚洲最大av| 熟女人妻精品中文字幕| 高清视频免费观看一区二区| 3wmmmm亚洲av在线观看| 妹子高潮喷水视频| 国产精品国产三级专区第一集| 99久国产av精品国产电影| 欧美激情极品国产一区二区三区| 久久午夜综合久久蜜桃| 另类亚洲欧美激情| 一边亲一边摸免费视频| 欧美av亚洲av综合av国产av | 欧美av亚洲av综合av国产av | 精品视频人人做人人爽| xxx大片免费视频| 亚洲图色成人| 亚洲天堂av无毛| 少妇人妻精品综合一区二区| 老鸭窝网址在线观看| 亚洲综合精品二区| 亚洲成色77777| 久久综合国产亚洲精品| 另类亚洲欧美激情| 一级毛片电影观看| 欧美黑人精品巨大| 免费久久久久久久精品成人欧美视频| 日本欧美视频一区| 国产精品人妻久久久影院| 黑人巨大精品欧美一区二区蜜桃| 国产熟女午夜一区二区三区| 日韩中文字幕欧美一区二区 | 国产激情久久老熟女| 亚洲av在线观看美女高潮| 国产精品一国产av| 一区二区三区激情视频| 欧美亚洲 丝袜 人妻 在线| 国产精品嫩草影院av在线观看| 亚洲av电影在线进入| 青青草视频在线视频观看| 午夜福利一区二区在线看| 各种免费的搞黄视频| 男人舔女人的私密视频| 免费在线观看黄色视频的| 精品国产超薄肉色丝袜足j| 欧美精品一区二区大全| 久久久久精品性色| 国产成人精品久久二区二区91 | 国产成人精品福利久久| av国产精品久久久久影院| 欧美激情 高清一区二区三区| 欧美中文综合在线视频| 欧美黑人精品巨大| 亚洲国产中文字幕在线视频| 亚洲精品久久成人aⅴ小说| 18禁裸乳无遮挡动漫免费视频| 亚洲图色成人| 日本vs欧美在线观看视频| 晚上一个人看的免费电影| 街头女战士在线观看网站| 水蜜桃什么品种好| av网站免费在线观看视频| 日本欧美国产在线视频| 男女免费视频国产| 国产片特级美女逼逼视频| 你懂的网址亚洲精品在线观看| 大片免费播放器 马上看| 国产在线免费精品| 国产免费视频播放在线视频| 午夜福利免费观看在线| 男女高潮啪啪啪动态图| 亚洲专区中文字幕在线 | 王馨瑶露胸无遮挡在线观看| 国产精品香港三级国产av潘金莲 | 国产成人免费观看mmmm| 国产黄色视频一区二区在线观看| 成人亚洲精品一区在线观看| 一本—道久久a久久精品蜜桃钙片| 最近的中文字幕免费完整| 19禁男女啪啪无遮挡网站| 777米奇影视久久| 超碰97精品在线观看| 国产熟女欧美一区二区| 夫妻午夜视频| 久久人人97超碰香蕉20202| 亚洲精品aⅴ在线观看| 啦啦啦在线免费观看视频4| 国产又爽黄色视频| 免费观看a级毛片全部| 老司机影院成人| 不卡视频在线观看欧美| 久久久久国产一级毛片高清牌| 黄网站色视频无遮挡免费观看| 国产精品av久久久久免费| 国产精品女同一区二区软件| 国产日韩欧美在线精品| 波多野结衣一区麻豆| 欧美人与性动交α欧美软件| avwww免费| 久久毛片免费看一区二区三区| 午夜老司机福利片| 97精品久久久久久久久久精品| 成年美女黄网站色视频大全免费| av一本久久久久| 丁香六月天网| 又大又爽又粗| 夜夜骑夜夜射夜夜干| 亚洲欧美成人精品一区二区| 国产精品秋霞免费鲁丝片| 久热爱精品视频在线9| 国产一区二区激情短视频 | 国产不卡av网站在线观看| 99热网站在线观看| 青春草国产在线视频| 免费观看av网站的网址| 久久热在线av| 国产精品久久久av美女十八| 少妇被粗大猛烈的视频| 国产精品99久久99久久久不卡 | netflix在线观看网站| 国产成人午夜福利电影在线观看| 男女高潮啪啪啪动态图| 国产午夜精品一二区理论片| 久久影院123| 黄片小视频在线播放| 黄片播放在线免费| av有码第一页| 午夜福利影视在线免费观看| 永久免费av网站大全| 777久久人妻少妇嫩草av网站| 久久免费观看电影| 热re99久久精品国产66热6| 亚洲色图 男人天堂 中文字幕| 美女视频免费永久观看网站| 久久久久精品久久久久真实原创| av在线app专区| 女人精品久久久久毛片| 国产黄色免费在线视频| 久久久久久免费高清国产稀缺| 制服诱惑二区| 国产爽快片一区二区三区| bbb黄色大片| 高清在线视频一区二区三区| 各种免费的搞黄视频| 日韩 欧美 亚洲 中文字幕| 亚洲国产中文字幕在线视频| 久久亚洲国产成人精品v| 黄片播放在线免费| 中文字幕人妻丝袜一区二区 | 狠狠婷婷综合久久久久久88av| 久久精品亚洲熟妇少妇任你| 1024视频免费在线观看| 欧美亚洲日本最大视频资源| 视频区图区小说| 97在线人人人人妻| 在线观看三级黄色| 久久久欧美国产精品| 欧美日韩亚洲国产一区二区在线观看 | 国产精品国产三级国产专区5o| 只有这里有精品99| 久久久久久人妻| 国产精品秋霞免费鲁丝片| 国产亚洲av片在线观看秒播厂| 自拍欧美九色日韩亚洲蝌蚪91| 国产淫语在线视频| 久久影院123| 在线天堂中文资源库| 人人妻人人添人人爽欧美一区卜| av有码第一页| 美国免费a级毛片| 国产av精品麻豆| 国产精品熟女久久久久浪| 亚洲一码二码三码区别大吗| 夜夜骑夜夜射夜夜干| 美女高潮到喷水免费观看| 99久久99久久久精品蜜桃| 捣出白浆h1v1| 人人妻人人添人人爽欧美一区卜| 国产精品二区激情视频| 精品国产国语对白av| 亚洲av日韩在线播放| 亚洲欧美日韩另类电影网站| 亚洲一卡2卡3卡4卡5卡精品中文| 精品亚洲成国产av| √禁漫天堂资源中文www| 一本一本久久a久久精品综合妖精| 最近中文字幕高清免费大全6| 亚洲av电影在线进入| 日本欧美国产在线视频| av在线观看视频网站免费| 亚洲国产欧美日韩在线播放| 另类精品久久| 伊人亚洲综合成人网| 麻豆精品久久久久久蜜桃| 最黄视频免费看| 国产亚洲av片在线观看秒播厂| 免费久久久久久久精品成人欧美视频| 人人妻人人爽人人添夜夜欢视频| 亚洲成人免费av在线播放| 各种免费的搞黄视频| 少妇被粗大的猛进出69影院| 看免费av毛片| 国产av码专区亚洲av| 国产精品嫩草影院av在线观看| 婷婷色麻豆天堂久久| 如何舔出高潮| 精品久久蜜臀av无| 亚洲,欧美精品.| 国产成人a∨麻豆精品| 欧美xxⅹ黑人| 久久久精品区二区三区| 日韩一卡2卡3卡4卡2021年| 日本欧美视频一区| 人人澡人人妻人| 午夜福利在线免费观看网站| 午夜激情av网站| 成人黄色视频免费在线看| 国产在线一区二区三区精| 大码成人一级视频| 亚洲国产精品999| 国产成人免费无遮挡视频| 大陆偷拍与自拍| 国产精品99久久99久久久不卡 | 久久精品国产a三级三级三级| 亚洲综合色网址| 18禁国产床啪视频网站| 午夜日韩欧美国产| 母亲3免费完整高清在线观看| 久久久精品区二区三区| 亚洲欧美成人综合另类久久久| 午夜影院在线不卡| 国产高清国产精品国产三级| 女性生殖器流出的白浆| a级毛片黄视频| 国产极品天堂在线| 999久久久国产精品视频| 亚洲av福利一区| 高清视频免费观看一区二区| 国产在视频线精品| 国产精品偷伦视频观看了| 十八禁高潮呻吟视频| 国产精品欧美亚洲77777| www.av在线官网国产| 亚洲四区av| 麻豆av在线久日| 亚洲欧美日韩另类电影网站| 午夜日本视频在线| 欧美日韩福利视频一区二区| 日韩 亚洲 欧美在线| 永久免费av网站大全| 亚洲自偷自拍图片 自拍| 天天添夜夜摸| 91国产中文字幕| 婷婷色综合www| 欧美最新免费一区二区三区| 久久国产亚洲av麻豆专区| 黄色视频不卡| 日本av手机在线免费观看| 日本黄色日本黄色录像| 国精品久久久久久国模美| 精品酒店卫生间| 91精品伊人久久大香线蕉| 国产精品一区二区在线观看99| 男男h啪啪无遮挡| 大陆偷拍与自拍| 亚洲精品国产av成人精品| 丝袜美腿诱惑在线| 亚洲综合精品二区| 尾随美女入室| 成人毛片60女人毛片免费| 精品久久久精品久久久| 成人亚洲欧美一区二区av| 日韩免费高清中文字幕av| 国产精品一区二区在线观看99| 日日摸夜夜添夜夜爱| 国产一区二区激情短视频 | 欧美黄色片欧美黄色片| 综合色丁香网| 老汉色av国产亚洲站长工具| 高清不卡的av网站| 波野结衣二区三区在线| 久久精品国产综合久久久| 大陆偷拍与自拍| 国产精品久久久久成人av| 纵有疾风起免费观看全集完整版| 狂野欧美激情性bbbbbb| av在线观看视频网站免费| 欧美亚洲日本最大视频资源| 久久精品久久久久久噜噜老黄| 亚洲精品一区蜜桃| 中国三级夫妇交换| 久久国产精品男人的天堂亚洲| 成人漫画全彩无遮挡| 一级毛片黄色毛片免费观看视频| 亚洲av综合色区一区| 国产男人的电影天堂91| 大香蕉久久网| 少妇人妻精品综合一区二区| 七月丁香在线播放| 亚洲成av片中文字幕在线观看| 日本欧美国产在线视频| 少妇被粗大猛烈的视频| 欧美变态另类bdsm刘玥| 午夜91福利影院| 青草久久国产| 午夜影院在线不卡| 黄色毛片三级朝国网站| 亚洲av中文av极速乱| 亚洲国产精品国产精品| 国产日韩一区二区三区精品不卡| 中文字幕人妻丝袜制服| 999精品在线视频| 亚洲av福利一区| 国产精品.久久久| 秋霞在线观看毛片| 成人国语在线视频| 啦啦啦 在线观看视频| 亚洲欧美中文字幕日韩二区| 亚洲av成人精品一二三区| 亚洲精华国产精华液的使用体验| 中文字幕精品免费在线观看视频| 欧美人与性动交α欧美精品济南到| 精品一品国产午夜福利视频| 精品国产国语对白av| 欧美精品一区二区免费开放| 国产精品香港三级国产av潘金莲 | 精品亚洲成a人片在线观看| 9色porny在线观看| 精品免费久久久久久久清纯 | 大话2 男鬼变身卡| 亚洲国产成人一精品久久久| svipshipincom国产片| 美女脱内裤让男人舔精品视频| 亚洲av综合色区一区| 国产激情久久老熟女| 国产高清不卡午夜福利| 亚洲精品国产区一区二| 美女扒开内裤让男人捅视频| 国产免费现黄频在线看| 性色av一级| 国产精品免费大片| 精品国产一区二区三区久久久樱花| 999精品在线视频| 免费av中文字幕在线| 免费观看人在逋| 久久久久久久国产电影| 成人18禁高潮啪啪吃奶动态图| 在线看a的网站| 亚洲精品一区蜜桃| 天美传媒精品一区二区| 色精品久久人妻99蜜桃| 亚洲伊人色综图| 天美传媒精品一区二区| 男女边摸边吃奶| 国产日韩欧美在线精品| 午夜日韩欧美国产| 欧美日韩视频精品一区| 另类亚洲欧美激情| 免费日韩欧美在线观看|