• <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中文字幕| 成人国语在线视频| bbb黄色大片| 日本撒尿小便嘘嘘汇集6| 国产区一区二久久| 两人在一起打扑克的视频| 久久久国产一区二区| 91麻豆av在线| 午夜免费鲁丝| 一本大道久久a久久精品| 一二三四在线观看免费中文在| 老司机靠b影院| 日本五十路高清| 久久中文看片网| 亚洲av日韩在线播放| 日韩有码中文字幕| svipshipincom国产片| 国产免费一区二区三区四区乱码| 女性被躁到高潮视频| 韩国高清视频一区二区三区| 欧美成人午夜精品| 亚洲午夜精品一区,二区,三区| 五月天丁香电影| 日日摸夜夜添夜夜添小说| 美女福利国产在线| 国产免费一区二区三区四区乱码| 一边摸一边抽搐一进一出视频| 久久人人爽人人片av| 久久久久久人人人人人| 91成人精品电影| 国产免费av片在线观看野外av| 美女高潮喷水抽搐中文字幕| 成人免费观看视频高清| 桃红色精品国产亚洲av| 亚洲欧美日韩高清在线视频 | 亚洲国产精品成人久久小说| 午夜影院在线不卡| 欧美日本中文国产一区发布| 少妇裸体淫交视频免费看高清 | 俄罗斯特黄特色一大片| 久久精品国产亚洲av高清一级| 国产精品久久久久久人妻精品电影 | av天堂在线播放| 久久国产精品男人的天堂亚洲| 久久精品aⅴ一区二区三区四区| a级片在线免费高清观看视频| 在线天堂中文资源库| 叶爱在线成人免费视频播放| 麻豆av在线久日| 欧美日韩一级在线毛片| 人妻久久中文字幕网| 男女国产视频网站| 建设人人有责人人尽责人人享有的| 一区二区日韩欧美中文字幕| 午夜福利在线观看吧| 美女扒开内裤让男人捅视频| 麻豆乱淫一区二区| 黄色视频不卡| 亚洲伊人久久精品综合| 欧美性长视频在线观看| 两个人免费观看高清视频| 国产福利在线免费观看视频| 国产高清videossex| 九色亚洲精品在线播放| 深夜精品福利| 成年人黄色毛片网站| 青青草视频在线视频观看| 日韩中文字幕欧美一区二区| 亚洲中文av在线| 精品熟女少妇八av免费久了| 亚洲av成人一区二区三| 国产成人系列免费观看| 一级,二级,三级黄色视频| 中文精品一卡2卡3卡4更新| 热99国产精品久久久久久7| 丁香六月天网| 亚洲欧美一区二区三区黑人| 国产高清国产精品国产三级| 精品国产超薄肉色丝袜足j| 午夜福利,免费看| 后天国语完整版免费观看| 久久99热这里只频精品6学生| 悠悠久久av| 久久久久精品人妻al黑| 老司机午夜十八禁免费视频| 亚洲国产精品成人久久小说| 成年动漫av网址| 人成视频在线观看免费观看| 久久影院123| 日韩中文字幕欧美一区二区| 久久久水蜜桃国产精品网| 国产熟女午夜一区二区三区| 国产av国产精品国产| 亚洲国产欧美日韩在线播放| 国内毛片毛片毛片毛片毛片| 91麻豆av在线| 动漫黄色视频在线观看| 性高湖久久久久久久久免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 男人爽女人下面视频在线观看| 国产野战对白在线观看| 交换朋友夫妻互换小说| 天堂8中文在线网| 侵犯人妻中文字幕一二三四区| 91麻豆精品激情在线观看国产 | 三级男女做爰猛烈吃奶摸视频| 免费av毛片视频| 18美女黄网站色大片免费观看| 国产av一区在线观看免费| 亚洲电影在线观看av| 天堂动漫精品| xxxwww97欧美| 三级男女做爰猛烈吃奶摸视频| 国产成年人精品一区二区| 久9热在线精品视频| 国产激情偷乱视频一区二区| 在线观看一区二区三区| 后天国语完整版免费观看| 精品久久久久久久人妻蜜臀av| 波多野结衣高清无吗| 香蕉丝袜av| 亚洲av电影不卡..在线观看| 99久久久亚洲精品蜜臀av| 国产一区二区激情短视频| x7x7x7水蜜桃| 国产精品99久久99久久久不卡| 国产精品一及| 国内精品久久久久久久电影| 国产乱人伦免费视频| 国产精品免费一区二区三区在线| 国内精品久久久久久久电影| 婷婷亚洲欧美| av有码第一页| 国产激情久久老熟女| 午夜两性在线视频| 欧美日韩亚洲综合一区二区三区_| 精品日产1卡2卡| 麻豆成人av在线观看| 天堂√8在线中文| 非洲黑人性xxxx精品又粗又长| 一本大道久久a久久精品| 久久亚洲精品不卡| 亚洲人成电影免费在线| 老熟妇乱子伦视频在线观看| а√天堂www在线а√下载| 女人被狂操c到高潮| 欧美不卡视频在线免费观看 | 99久久99久久久精品蜜桃| 国产真人三级小视频在线观看| 最近最新免费中文字幕在线| 亚洲中文日韩欧美视频| 色综合亚洲欧美另类图片| 99在线视频只有这里精品首页| 波多野结衣巨乳人妻| 久久人人精品亚洲av| 免费搜索国产男女视频| 50天的宝宝边吃奶边哭怎么回事| 国产三级在线视频| 久久午夜亚洲精品久久| 超碰成人久久| 亚洲熟妇熟女久久| 色综合欧美亚洲国产小说| 欧美成人午夜精品| 男女之事视频高清在线观看| 国产精品电影一区二区三区| 午夜a级毛片| 一本一本综合久久| 99久久99久久久精品蜜桃| 午夜a级毛片| 亚洲欧美日韩高清专用| 老司机在亚洲福利影院| 天天添夜夜摸| 精品一区二区三区av网在线观看| 美女免费视频网站| 久99久视频精品免费| 在线十欧美十亚洲十日本专区| 宅男免费午夜| 久久久久性生活片| 在线观看www视频免费| 露出奶头的视频| 精品日产1卡2卡| 亚洲精品中文字幕在线视频| av视频在线观看入口| 欧美国产日韩亚洲一区| 午夜视频精品福利| av超薄肉色丝袜交足视频| 中文字幕久久专区| 午夜a级毛片| 毛片女人毛片| 久久久久久久久免费视频了| 久久久久免费精品人妻一区二区| 女人被狂操c到高潮| 亚洲精品av麻豆狂野| 亚洲成人免费电影在线观看| 99久久精品热视频| 99久久无色码亚洲精品果冻| 亚洲第一欧美日韩一区二区三区| 在线观看舔阴道视频| 国产精品98久久久久久宅男小说| 国产午夜福利久久久久久| 久久精品夜夜夜夜夜久久蜜豆 | 久久午夜亚洲精品久久| 丰满人妻熟妇乱又伦精品不卡| av福利片在线观看| 美女免费视频网站| 制服丝袜大香蕉在线| 精品一区二区三区四区五区乱码| 可以在线观看毛片的网站| 一边摸一边做爽爽视频免费| avwww免费| 搡老妇女老女人老熟妇| 久久久久久亚洲精品国产蜜桃av| 欧美乱码精品一区二区三区| 18禁黄网站禁片免费观看直播| 午夜福利视频1000在线观看| 国语自产精品视频在线第100页| 久久久久久国产a免费观看| 99热这里只有精品一区 | 亚洲精品一区av在线观看| 国模一区二区三区四区视频 | 亚洲精品一区av在线观看| www.www免费av| 欧美日本视频| 国产精品日韩av在线免费观看| 国产激情久久老熟女| 亚洲 国产 在线| 亚洲精华国产精华精| 色综合婷婷激情| 国产91精品成人一区二区三区| 国产视频内射| 国产高清视频在线播放一区| 91国产中文字幕| 国产成+人综合+亚洲专区| 精品乱码久久久久久99久播| 法律面前人人平等表现在哪些方面| ponron亚洲| 色播亚洲综合网| 午夜日韩欧美国产| 国产欧美日韩精品亚洲av| 美女扒开内裤让男人捅视频| www.自偷自拍.com| 好男人在线观看高清免费视频| 婷婷亚洲欧美| 国内毛片毛片毛片毛片毛片| 国产99久久九九免费精品| 日韩精品免费视频一区二区三区| av片东京热男人的天堂| 精品电影一区二区在线| 欧美乱妇无乱码| 国产精品亚洲美女久久久| 亚洲国产欧洲综合997久久,| 亚洲第一电影网av| 日韩 欧美 亚洲 中文字幕| 亚洲人与动物交配视频| 又紧又爽又黄一区二区| 亚洲九九香蕉| 精品高清国产在线一区| 久久精品影院6| 欧美极品一区二区三区四区| 久久欧美精品欧美久久欧美| 亚洲一卡2卡3卡4卡5卡精品中文| 久久精品91蜜桃| 99国产精品一区二区蜜桃av| 国产精品亚洲美女久久久| 九色成人免费人妻av| 国内久久婷婷六月综合欲色啪| 女人高潮潮喷娇喘18禁视频| 日本一区二区免费在线视频| 人妻夜夜爽99麻豆av| 亚洲欧美一区二区三区黑人| 一本综合久久免费| 舔av片在线| 成人18禁高潮啪啪吃奶动态图| 欧美日本亚洲视频在线播放| 18禁黄网站禁片午夜丰满| 国产亚洲av高清不卡| 叶爱在线成人免费视频播放| 91在线观看av| 亚洲人成77777在线视频| 亚洲欧美日韩高清在线视频| 国产精品久久久人人做人人爽| 国产精品亚洲av一区麻豆| 成人av一区二区三区在线看| 在线观看免费日韩欧美大片| 国产亚洲精品久久久久久毛片| 性色av乱码一区二区三区2| 淫妇啪啪啪对白视频| 亚洲精品粉嫩美女一区| 日本一本二区三区精品| 国产精品久久久人人做人人爽| 亚洲第一欧美日韩一区二区三区| 久久国产乱子伦精品免费另类| 欧美人与性动交α欧美精品济南到| 日本成人三级电影网站| 亚洲国产欧美一区二区综合| 亚洲一区中文字幕在线| 男女之事视频高清在线观看| 在线观看舔阴道视频| 久久久水蜜桃国产精品网| 一本精品99久久精品77| 黄色女人牲交| 精品久久蜜臀av无| 99国产极品粉嫩在线观看| 精品久久久久久久末码| 老司机在亚洲福利影院| 母亲3免费完整高清在线观看| 欧美绝顶高潮抽搐喷水| 免费看美女性在线毛片视频| 9191精品国产免费久久| 天天添夜夜摸| 久久精品国产综合久久久| 亚洲精品粉嫩美女一区| 可以在线观看毛片的网站| 岛国视频午夜一区免费看| 亚洲男人的天堂狠狠| 成人三级黄色视频| 精品欧美一区二区三区在线| 巨乳人妻的诱惑在线观看| 亚洲免费av在线视频| 欧美一区二区精品小视频在线| 青草久久国产| 国产午夜福利久久久久久| 色精品久久人妻99蜜桃| 曰老女人黄片| 亚洲专区字幕在线| 久久这里只有精品中国| 美女扒开内裤让男人捅视频| 中文字幕人成人乱码亚洲影| 人妻久久中文字幕网| 无遮挡黄片免费观看| 国产成人一区二区三区免费视频网站| 伊人久久大香线蕉亚洲五| 久久午夜亚洲精品久久| 悠悠久久av| 精品欧美国产一区二区三| 熟女电影av网| 久久性视频一级片| 999精品在线视频| 日韩大码丰满熟妇| 日韩中文字幕欧美一区二区| 毛片女人毛片| 免费在线观看影片大全网站| 后天国语完整版免费观看| 亚洲熟女毛片儿| 亚洲国产欧美一区二区综合| 国产精品免费视频内射| 国产精品久久视频播放| 麻豆成人av在线观看| 亚洲精华国产精华精| 成人高潮视频无遮挡免费网站| 香蕉丝袜av| netflix在线观看网站| 床上黄色一级片| 啦啦啦观看免费观看视频高清| 精品乱码久久久久久99久播| 久久久精品大字幕| 九色成人免费人妻av| 亚洲 欧美一区二区三区| 天堂动漫精品| 国产激情偷乱视频一区二区| 舔av片在线| 在线观看舔阴道视频| 亚洲性夜色夜夜综合| 国产亚洲精品综合一区在线观看 | 久久国产精品人妻蜜桃| 在线免费观看的www视频| 中文字幕最新亚洲高清| 五月玫瑰六月丁香| 脱女人内裤的视频| 久久久久久大精品| 黑人操中国人逼视频| 欧美性猛交╳xxx乱大交人| 91九色精品人成在线观看| 成在线人永久免费视频| 成人永久免费在线观看视频| 18禁美女被吸乳视频| 在线十欧美十亚洲十日本专区| 狠狠狠狠99中文字幕| www日本在线高清视频| 成人午夜高清在线视频| 免费无遮挡裸体视频| a级毛片在线看网站| 国产精品一区二区精品视频观看| 欧美黑人巨大hd| 精品福利观看| 日韩大码丰满熟妇| 我的老师免费观看完整版| 别揉我奶头~嗯~啊~动态视频| 人妻夜夜爽99麻豆av| 在线观看免费午夜福利视频| 午夜免费成人在线视频| 欧美中文日本在线观看视频| 巨乳人妻的诱惑在线观看| 亚洲avbb在线观看| 日本 av在线| 国产久久久一区二区三区| 免费av毛片视频| 男女那种视频在线观看| 亚洲精品粉嫩美女一区| 韩国av一区二区三区四区| 丝袜美腿诱惑在线| 国产精品乱码一区二三区的特点| 老鸭窝网址在线观看| 变态另类成人亚洲欧美熟女| 这个男人来自地球电影免费观看| 不卡一级毛片| 免费在线观看黄色视频的| 嫩草影视91久久| 桃色一区二区三区在线观看| 国产av在哪里看| 狂野欧美激情性xxxx| 一a级毛片在线观看| 日本三级黄在线观看| 在线永久观看黄色视频| 啦啦啦观看免费观看视频高清| 国产一区二区在线av高清观看| 一级毛片精品| 在线观看美女被高潮喷水网站 | 男女床上黄色一级片免费看| 一区二区三区国产精品乱码| 久久中文字幕一级| 国语自产精品视频在线第100页| 十八禁人妻一区二区| 亚洲精品在线美女| 一进一出抽搐动态| 午夜福利18| 99久久综合精品五月天人人| 91麻豆精品激情在线观看国产| 高潮久久久久久久久久久不卡| 黄色视频不卡| 又黄又粗又硬又大视频| 午夜福利在线观看吧| 亚洲人成网站在线播放欧美日韩| 真人一进一出gif抽搐免费| 不卡一级毛片| 欧美色欧美亚洲另类二区| 巨乳人妻的诱惑在线观看| 日本三级黄在线观看| 久久久久久久久久黄片| 波多野结衣巨乳人妻| 国产亚洲精品久久久久5区| 精品福利观看| 1024手机看黄色片| 国产野战对白在线观看| 老司机午夜福利在线观看视频| 国产男靠女视频免费网站| 色综合欧美亚洲国产小说| 久久天堂一区二区三区四区| 国产av麻豆久久久久久久| 国产精品av久久久久免费| 欧美激情久久久久久爽电影| 日韩大码丰满熟妇| 伊人久久大香线蕉亚洲五| www国产在线视频色| av欧美777| xxx96com| 最近视频中文字幕2019在线8| 亚洲最大成人中文| 在线观看免费日韩欧美大片| 精品电影一区二区在线| 99久久综合精品五月天人人| 久久久久久久久久黄片| 黑人操中国人逼视频| 香蕉丝袜av| av福利片在线观看| 操出白浆在线播放| 国产真人三级小视频在线观看| 中文字幕最新亚洲高清| www.999成人在线观看| 久久久久久久久中文| 亚洲精品中文字幕一二三四区| 变态另类成人亚洲欧美熟女| 亚洲成人久久爱视频| 少妇人妻一区二区三区视频| 中文在线观看免费www的网站 | 老司机深夜福利视频在线观看| 1024视频免费在线观看| 12—13女人毛片做爰片一| 久久久久久亚洲精品国产蜜桃av| 欧美午夜高清在线| 久久久久久久久中文| 亚洲一码二码三码区别大吗| 不卡一级毛片| 欧美3d第一页| 1024香蕉在线观看| 亚洲av片天天在线观看| 国产精品 国内视频| 久久婷婷成人综合色麻豆| 免费在线观看成人毛片| 久久精品aⅴ一区二区三区四区| 精品国产乱子伦一区二区三区| 久久久精品欧美日韩精品| 亚洲色图av天堂| 亚洲色图 男人天堂 中文字幕| 亚洲精品一区av在线观看| 国产精品永久免费网站| 国产亚洲欧美在线一区二区| 老熟妇乱子伦视频在线观看| 成人av在线播放网站| 日韩欧美精品v在线| 国产精品永久免费网站| 亚洲无线在线观看| 婷婷精品国产亚洲av| 亚洲人成电影免费在线| 国产不卡一卡二| 91老司机精品| 91av网站免费观看| 动漫黄色视频在线观看| 首页视频小说图片口味搜索| 国产精品免费一区二区三区在线| 亚洲成av人片在线播放无| 中亚洲国语对白在线视频| 精品电影一区二区在线| 国产精品久久电影中文字幕| 国产一区二区激情短视频| 国内精品一区二区在线观看| 久久午夜亚洲精品久久| 搡老岳熟女国产| 国产成人av激情在线播放| 国产亚洲精品久久久久5区| 亚洲国产高清在线一区二区三| 亚洲成a人片在线一区二区| 他把我摸到了高潮在线观看| 亚洲午夜理论影院| a在线观看视频网站| 亚洲av成人不卡在线观看播放网| 国产av一区在线观看免费| a级毛片a级免费在线| 亚洲精品一卡2卡三卡4卡5卡| 午夜精品久久久久久毛片777| 国产真人三级小视频在线观看| 日韩大尺度精品在线看网址| 亚洲欧美一区二区三区黑人| 亚洲成人久久性| 欧美极品一区二区三区四区| 女生性感内裤真人,穿戴方法视频| 激情在线观看视频在线高清| 一进一出抽搐动态| 少妇被粗大的猛进出69影院| 91在线观看av| 在线观看免费日韩欧美大片| 在线视频色国产色| 看免费av毛片| 亚洲黑人精品在线| 天天躁夜夜躁狠狠躁躁| 亚洲第一电影网av| 国产精品精品国产色婷婷| 全区人妻精品视频| 99久久精品国产亚洲精品| 国产91精品成人一区二区三区| 欧美日韩国产亚洲二区| 成年人黄色毛片网站| 色在线成人网| 一级毛片女人18水好多| 啪啪无遮挡十八禁网站| 亚洲av成人不卡在线观看播放网| 一本精品99久久精品77| 国产亚洲精品久久久久5区| 亚洲一区二区三区色噜噜| 亚洲欧美精品综合久久99| 国产一区二区在线av高清观看| 国语自产精品视频在线第100页| 99热6这里只有精品| 黄色女人牲交| 国产亚洲精品综合一区在线观看 | 亚洲精品在线观看二区| 最新美女视频免费是黄的| 国产黄色小视频在线观看| 动漫黄色视频在线观看| 精品日产1卡2卡| 欧美一区二区国产精品久久精品 | 99国产精品一区二区蜜桃av| 伊人久久大香线蕉亚洲五| 在线播放国产精品三级| 叶爱在线成人免费视频播放| 老汉色∧v一级毛片| 中出人妻视频一区二区| 国产三级黄色录像| 757午夜福利合集在线观看| 欧美人与性动交α欧美精品济南到| 午夜a级毛片| 9191精品国产免费久久| 成人手机av| 中文字幕人成人乱码亚洲影| 亚洲中文日韩欧美视频| 亚洲五月婷婷丁香| 国产又色又爽无遮挡免费看| 中文字幕av在线有码专区| 日韩高清综合在线| 丁香欧美五月| 国产1区2区3区精品| 精品久久久久久久末码| 日韩欧美精品v在线| 亚洲精品av麻豆狂野| 久久久久国内视频| 亚洲成人中文字幕在线播放| 日日摸夜夜添夜夜添小说| 色精品久久人妻99蜜桃| 天天躁狠狠躁夜夜躁狠狠躁| 久久人妻av系列| 黄频高清免费视频| 日本免费一区二区三区高清不卡| 亚洲天堂国产精品一区在线| 51午夜福利影视在线观看| 欧美+亚洲+日韩+国产| 国产精品一区二区免费欧美| 日韩av在线大香蕉| 免费看a级黄色片| 夜夜躁狠狠躁天天躁| 狠狠狠狠99中文字幕| 午夜福利在线在线| 人妻丰满熟妇av一区二区三区| 色精品久久人妻99蜜桃| 国产成人系列免费观看| 亚洲av中文字字幕乱码综合|