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

    基于VB的坐標轉換程序的開發(fā)

    2013-03-06 02:11:59田繼成
    城市勘測 2013年4期
    關鍵詞:源代碼橢球高斯

    田繼成

    (大同市勘察測繪院,山西大同 037000)

    基于VB的坐標轉換程序的開發(fā)

    田繼成?

    (大同市勘察測繪院,山西大同 037000)

    網(wǎng)上免費的坐標轉換軟件很多,功能各異,其中COORD功能強大、使用廣泛,是坐標轉換軟件中的杰出代表。盡管如此,COORD在有些情況下仍不能實現(xiàn)一步轉換。因此,有必要開發(fā)一種能夠實現(xiàn)一步轉換的坐標轉換程序。本文以VB 6.0為開發(fā)環(huán)境,重點闡述了開發(fā)坐標轉換程序的思路和算法。程序投入運行以來,極大地提高了坐標轉換工作的效率,也有利于坐標轉換成果的標準化。

    坐標轉換;程序設計;模塊化;VB 6.0

    1 引 言

    坐標轉換問題在測量中經(jīng)常遇到,計算過程較復雜,一般使用軟件計算。網(wǎng)上免費的坐標轉換軟件很多,功能各異,其中COORD功能強大、使用廣泛,是坐標轉換軟件中的杰出代表。

    盡管如此,COORD在有些情況下仍不能實現(xiàn)一步轉換,如:轉換過程中需要進行兩次不同自定義的高斯正算和反算時。由于大同市獨立坐標系采用非國家標準3°帶中央子午線投影,故大同獨立坐標系轉西安80坐標系即屬于這種情況。因此,有必要開發(fā)一種能夠實現(xiàn)一步轉換的坐標轉換程序。

    為了使開發(fā)的程序精確可靠、功能強大、易于使用,并縮短程序開發(fā)周期,開發(fā)過程中筆者閱讀了大量坐標轉換文獻,對比和借鑒了大量坐標轉換軟件,測試分析了大量坐標轉換源代碼,并在此基礎上對源代碼進行了改進和利用。

    本文以VB 6.0為開發(fā)環(huán)境,重點闡述了開發(fā)坐標轉換程序的思路和算法。程序投入運行以來,極大地提高了大同市勘察測繪院坐標轉換工作的效率,也有利于坐標轉換成果的標準化。

    2 功能及界面設計

    程序應實現(xiàn):

    (1)一步坐標轉換

    程序實行模塊化,通過對坐標轉換子模塊的任意組合,實現(xiàn)一步坐標轉換。

    (2)實時坐標轉換

    當輸入或改變源坐標時,實時計算當前行目標坐標。改變源坐標系或目標坐標系時,實時全部重新計算目標坐標。改變保留小數(shù)位數(shù)或角度格式后,實時刷新數(shù)據(jù)。

    (3)轉換參數(shù)的安全

    轉換參數(shù)加密內(nèi)置在程序中,既方便使用和防止誤操作,又避免了使用明碼參數(shù)帶來的易泄露和傳播等不安全問題。

    為防止對程序進行逆向工程,程序采用P碼編譯并加殼處理。通過第三方軟件實現(xiàn)程序的防拷貝。通過設定登錄密碼、綁定微機硬件、使用次數(shù)限制、時間炸彈等措施,實現(xiàn)程序非授權禁用并自毀。

    程序界面如圖1所示:

    圖1 坐標轉換

    3 關鍵技術及實現(xiàn)

    3.1 模塊化坐標轉換

    本程序坐標轉換分為專用坐標轉換和通用坐標轉換兩部分。專用坐標轉換中預置了大同市常用的坐標轉換參數(shù),通用坐標轉換需用戶定制后才可使用。

    將坐標轉換模塊劃分4大類:四參數(shù)轉換、七參數(shù)轉換、大地坐標系與空間直角坐標系互轉、高斯正反算。

    (1)專用坐標轉換:支持6種坐標(大同地方坐標系、大同獨立坐標系、北京54坐標系、西安80坐標系、WGS84坐標系、CGCS2000坐標系)互轉。將6種坐標互轉過程按4大類劃分為28個模塊,編程時指定每兩種坐標轉換時對應的模塊組合,即可實現(xiàn)專用坐標轉換。

    以大同獨立坐標系轉西安80坐標系為例,需要6步:

    ①四參數(shù)轉換:大同獨立坐標轉WGS84橢球的高斯投影平面坐標;

    ②高斯反算:WGS84橢球的高斯投影平面坐標轉WGS84大地坐標;

    ③大地轉空間:WGS84大地坐標轉WGS84空間直角坐標;

    ④七參數(shù)轉換:WGS84空間直角坐標轉西安80空間直角坐標;

    ⑤空間轉大地:西安80空間直角坐標轉西安80大地坐標;

    ⑥高斯正算:西安80大地坐標轉西安80平面坐標。

    (注:WGS84橢球的高斯投影平面坐標,具體為高斯投影3°帶,中央子午線???°??′,投影面高程????m。)

    (2)通用坐標轉換:按4大類劃分為6個模塊:四參數(shù)轉換、七參數(shù)轉換、大地轉空間、空間轉大地、高斯正算、高斯反算。使用時用戶自定義模塊組合,即可實現(xiàn)通用坐標轉換。添加模塊后,需設置模塊對應參數(shù),如:四參數(shù)轉換需輸入四參數(shù)、高斯正算需輸入投影參數(shù)等。全部設置好后,可保存為坐標轉換配置文件,下次使用時選擇即可。坐標轉換配置文件為加密格式。

    圖2 定制轉換方案

    3.2 高斯投影代碼的改進

    網(wǎng)上收集了很多高斯投影源代碼,經(jīng)測試或多或少存在不足,如:代碼存在漏錯、代碼格式混亂、支持橢球較少、不支持自定義高斯投影、不支持橢球變換、轉換精度不夠等。

    經(jīng)比較,選定一款較好的高斯投影源代碼,在此基礎上進行改進,實現(xiàn)了:

    (1)增加支持的橢球

    源代碼只支持北京54坐標系的參考和西安80坐標系的參考橢球,在此基礎上,增加了對WGS84和CGCS 2000橢球的支持。

    (2)提高轉換精度

    經(jīng)對比發(fā)現(xiàn),源代碼中用于計算X(B)的高階項系數(shù)只有A1、A2、A3、A4,通過增加A5、A6,提高了轉換精度。

    (3)支持橢球變換

    橢球變換主要有橢球平移法、橢球膨脹法、橢球變形法等,其中橢球膨脹法最為簡潔,為工程應用之首選??紤]到坐標轉換軟件COORD采用橢球膨脹法,為了使用COORD求解的坐標轉換參數(shù),故本程序橢球變換也采用橢球膨脹法。

    橢球膨脹法公式:

    式中:Hm—該地區(qū)平均大地高;Bm—該地區(qū)的平均緯度。

    改進后的高斯投影代碼如下(限于篇幅,部分省略):

    '設置橢球參數(shù)

    Public Sub SetEllipsoidPara(ByVal CoordID%)

    Select Case CoordID

    Case 1 '北京54 a=6378245# f=298.3

    Case 2 '西安80 a=6378140# f=298.257

    Case 3 'WGS84 a=6378137# f=298.257223563

    Case 4 'CGCS2000 a=6378137# f=298.257222101

    End Select

    End Sub

    '設置坐標系統(tǒng)參數(shù)

    Public Sub SetCoordPara(ByVal CoordID%,ByVal Hm#,By-Val Bm#)

    Dim TA#,TB#,TC#,TD#,TE#,TF#

    Dim Rm#,Rn#

    '設置橢球參數(shù)

    Call SetEllipsoidPara(CoordID) '橢球變換(橢球膨脹法)

    e2=1-((f-1)/f)?((f-1)/f)

    Rm=a?Sqr(1-e2)/(1-e2?Sin(Bm)^2)

    Rn=Rm+Hm

    a=Rn?(1-e2?Sin(Bm)^2)/Sqr(1-e2)

    '計算用于計算X的橢球參數(shù)

    e2=1-((f-1)/f)?((f-1)/f)

    e12=(f/(f-1))?(f/(f-1))-1

    TA=a?(1-e2)?(1+3?e2/4+45?e2^2/64+175?e2^3/256+11025?e2^4/16384)

    TB=a?(1-e2)?(3?e2/4+15?e2^2/16+525?e2^3/ 512+2205?e2^4/2048)

    TC=a?(1-e2)?(15?e2^2/64+105?e2^3/256+2205?e2^4/4096+10395?e2^5/16384)

    TD=a?(1-e2)?(35?e2^3/512+315?e2^4/2048+ 31185?e2^5/131072)

    TE=a?(1-e2)?(315?e2^4/16384+3465?e2^5/ 65536)

    TF=a?(1-e2)?(693?e2^5/131072)

    A1=TA

    A2=-TB/2

    A3=TC/4

    A4=-TD/6

    A5=TE/8

    A6=-TF/10

    End Sub

    '高斯正算

    Public Sub BL2XY(ByRef B#,ByRef L#,ByVal L0#,ByVal Y0#,ByVal Hm#,ByVal Bm#,ByVal Coord ID%)

    ………………

    x=A1?B+A2?Sin(2?B)+A3?Sin(4?B)+A4?Sin (6?B)+A5?Sin(8?B)+A6?Sin(10?B)

    ………………End Sub

    '高斯反算

    Public Sub XY2BL(ByRef x#,ByRef y#,ByVal L0#,ByVal Y0#,ByVal Hm#,ByVal Bm#,ByVal CoordID%)

    ………………

    DoWhile True

    preB0=B0

    B0=B0?pi/180

    B0=(x-(A2?Sin(2?B0)+A3?Sin(4?B0)+A4?Sin(6?B0)+A5?Sin(8?B0)+A6?Sin(10?B0)))/A1?180/pi

    eta=Abs(B0-preB0)

    If eta<1E-16 Then Exit Do

    Loop

    ………………

    End Sub

    3.3 角度轉度分秒函數(shù)的改進

    度轉度分秒的難點在于容易出現(xiàn)60″。出現(xiàn)60″的原因:VB雙精度變量最多存儲16位(不包括正負號),秒單獨存儲時保留小數(shù)較多,此時的秒小于60″。秒和度分組合后,由于雙精度最大存儲限制導致秒四舍五入,此時的秒有可能等于60″。

    避免出現(xiàn)60″的方法:在秒和度分組合前,對秒進行四舍五入,并判斷其是否為60″。實現(xiàn)保留最多小數(shù)位數(shù)的計算公式:小數(shù)位數(shù)=16-度的位數(shù)-分秒的位數(shù)4-小數(shù)點位數(shù)1。

    通過對源代碼的修改,實現(xiàn)了避免出現(xiàn)60″、保留最多小數(shù)位數(shù)、支持負角度。

    修改后的角度轉度分秒代碼如下:

    Public Function Deg2Dms(ByVal Deg As Double)As Double

    ………………

    Deg=Abs(Deg) DD=Int(Deg)

    MM=Int((Deg-DD)?60#)

    SS=((Deg-DD)?60#-MM)?60#

    SecDec=16-Len(CStr(DD))-5

    If Round(SS,SecDec)=60 Then

    SS=0:MM=MM+1

    End If

    If MM=60 Then MM=0:DD=DD+1

    End If

    Deg2Dms=Sign?(DD+MM/100+SS/10000)

    End Function

    3.4 支持“一轉多”和“多轉一”

    支持源坐標同時轉換多種目標坐標,簡稱“一轉多”。支持源坐標按多種源坐標系轉換一種目標坐標,簡稱“多轉一”。

    “一轉多”用于同時轉換兩種或兩種以上坐標,避免了多次轉換的麻煩,而且成果可生成在一張報表中。

    “多轉一”用于源坐標系的輔助確定,有時無法確定源坐標系,但知道具體位置時,可按不同坐標系轉換后,與地形圖匹配確定源坐標系,再進行坐標的轉換。

    3.5 文件的導入導出

    程序除支持常用的導入導出格式外,還支持對用戶鍵入格式的解析,用戶鍵入的格式要求以逗號分隔。

    用戶鍵入的導入格式,可以是6種關鍵字(序號、點名、代碼、X、Y、H)的自由組合。

    用戶鍵入的導出格式,可以是9種關鍵字(序號、點名、代碼、前X、前Y、前H、后X、后Y、后H)的自由組合。其中,“前”代表轉換前的坐標,“后”代表轉換后的坐標。

    此外還支持輸出Word格式的坐標轉換報表,表格的列寬和字號隨坐標系名稱、坐標系類型、小數(shù)位數(shù)、角度格式等智能調整。

    3.6 生成經(jīng)緯網(wǎng)格

    經(jīng)緯網(wǎng)格可用于專題圖的制作、柵格圖的糾正拼接等。較復雜的經(jīng)緯網(wǎng)格一般使用ArcGIS生成,為了能夠直接生成符合要求的、DWG格式的經(jīng)緯網(wǎng)格,專門開發(fā)了生成經(jīng)緯網(wǎng)格功能。

    本功能支持多種橢球、支持自定義高斯投影、支持注記網(wǎng)格經(jīng)緯度,可自定義網(wǎng)格范圍、經(jīng)線緯線間隔、設置網(wǎng)格線樣式等。

    圖3 繪制經(jīng)緯風格

    4 程序測試

    筆者對程序的各個模塊進行了大量獨立測試和聯(lián)合測試,并用大同市C級GPS網(wǎng)多套坐標成果對程序進行了全面驗證。大同市C級GPS網(wǎng)覆蓋面積約1 500 km2,東西跨度57 km,南北跨度36 km。用該成果驗證轉換精度,具有很好的代表性。

    作者編制了較完整程序測試報告,并制作了轉換精度分布圖。測試結果表明,程序穩(wěn)定可靠、轉換精度較高。由于報告內(nèi)容較多,限于篇幅,僅以大同獨立坐標系轉西安80坐標系為例簡要介紹。

    測試方法如下:

    (1)轉換精度:輸入源坐標,轉換得到目標坐標,將目標坐標與控制點坐標進行比較。

    (2)逆轉相符性:輸入源坐標,轉換得到目標坐標,再將目標坐標轉換為源坐標,將轉換得到的源坐標與輸入的源坐標進行比較。

    (3)源坐標高程對平面精度的影響:輸入源坐標時,高程分別輸入0、1 000、大地高,將3種不同高程時得到的目標坐標進行比較。

    (4)COORD轉換精度:在COORD中輸入與本程序一樣的坐標轉換參數(shù),將COORD得到目標坐標與控制點坐標進行比較。

    (5)本程序與COORD比較:在COORD中輸入與本程序一樣的坐標轉換參數(shù),將COORD得到目標坐標與本程序得到目標坐標進行比較,如表1所示。

    坐標轉換測試項目及精度 表1

    5 結 語

    投入使用以來,程序在不斷地改進和完善,目前仍有諸多不足之處,如不支持坐標轉換參數(shù)的求解,不支持高程的轉換等。今后筆者會適時增加這方面的功能。

    通過此次開發(fā),筆者掌握了全套源代碼,通過定制可服務不同需要,對于提高工作效率具有重要意義。此外,通過實踐加深了對坐標轉換原理和方法的認識。

    [1] 王順興,顏美代.精解Word 2000 VBA與范例解析[M].北京:北京大學出版社,2001.

    [2] 張晉西.Visual Basic與AutoCAD二次開發(fā)[M].北京:清華大學出版社,2002.

    [3] 何光渝.Visual Basic常用數(shù)值算法集[M].北京:科學出版社,2002.

    [4] 許謖.Visual Basic應用與開發(fā)案例教程[M].北京:清華大學出版社,2005.

    [5] 李懷明,駱原,王育新.Visual Basic 6.0中文版參考詳解[M].北京:清華大學出版社,1999.

    [6] 陳涌泉.Visual Basic第三方控件大全[M].北京:中國電力出版社,200.

    The Program Development of Coordinate Transformation Based on VB

    Tian Jicheng
    (Datong Institute of Prospecting and Mapping,Datong 037000,China)

    There's lots of free coordinate transformation software with different functions.As an outstanding representative of the software,COORD is powerful and widely used.But,in some cases,COORD still can not achieve one-step conversion.Therefore,it is necessary to develop a coordinate transformation program to achieve one-step conversion. Based on VB 6.0,this paper focusing on the thoughtand algorithms of development coordinate transformation programs.It greatly improves thework efficiency of the coordinate transformation since put into operation,and conducive to the standardization of coordinate transformation results.

    coordinate transformation;program development;modularization;VB 6.0

    1672-8262(2013)04-71-04

    P226+.3

    B

    2013—03—09

    田繼成(1965—),男,高級工程師,注冊規(guī)劃師,主要從事生產(chǎn)管理和技術研究工作。

    猜你喜歡
    源代碼橢球高斯
    小高斯的大發(fā)現(xiàn)
    人工智能下復雜軟件源代碼缺陷精準校正
    計算機仿真(2023年8期)2023-09-20 11:23:42
    獨立坐標系橢球變換與坐標換算
    橢球槽宏程序編制及其Vericut仿真
    智能制造(2021年4期)2021-11-04 08:54:44
    基于TXL的源代碼插樁技術研究
    天才數(shù)學家——高斯
    軟件源代碼非公知性司法鑒定方法探析
    橢球精加工軌跡及程序設計
    基于外定界橢球集員估計的純方位目標跟蹤
    揭秘龍湖產(chǎn)品“源代碼”
    黑丝袜美女国产一区| 九九爱精品视频在线观看| 黑人巨大精品欧美一区二区蜜桃| 男女床上黄色一级片免费看| 亚洲激情五月婷婷啪啪| 一级毛片电影观看| 免费日韩欧美在线观看| 亚洲专区中文字幕在线 | av免费观看日本| www.熟女人妻精品国产| 亚洲四区av| 国产一级毛片在线| 在线天堂中文资源库| 一区二区三区精品91| 国产亚洲午夜精品一区二区久久| 国产 精品1| 免费少妇av软件| 最近中文字幕2019免费版| av福利片在线| 日本av手机在线免费观看| 秋霞在线观看毛片| 国产精品二区激情视频| 欧美日韩视频精品一区| 观看美女的网站| 一边摸一边做爽爽视频免费| 尾随美女入室| 欧美精品一区二区免费开放| 亚洲精品av麻豆狂野| 啦啦啦啦在线视频资源| 不卡视频在线观看欧美| 亚洲精品,欧美精品| www.熟女人妻精品国产| 国产亚洲一区二区精品| 亚洲在久久综合| 国产极品天堂在线| www.熟女人妻精品国产| 亚洲成人一二三区av| a级片在线免费高清观看视频| 亚洲国产欧美一区二区综合| 国产片内射在线| 老司机在亚洲福利影院| 久热这里只有精品99| 久久精品国产亚洲av高清一级| 十分钟在线观看高清视频www| 亚洲婷婷狠狠爱综合网| www.自偷自拍.com| 亚洲人成77777在线视频| 亚洲美女黄色视频免费看| 在线观看国产h片| 一二三四在线观看免费中文在| 男女无遮挡免费网站观看| 亚洲七黄色美女视频| 中文字幕色久视频| 丝袜脚勾引网站| 美女中出高潮动态图| 亚洲欧美一区二区三区黑人| 国产高清国产精品国产三级| 国产无遮挡羞羞视频在线观看| 亚洲,欧美精品.| 日韩大码丰满熟妇| 丰满迷人的少妇在线观看| 在线免费观看不下载黄p国产| 日韩一区二区视频免费看| 国产成人精品久久二区二区91 | 久久精品国产a三级三级三级| 啦啦啦 在线观看视频| 国产一区二区三区av在线| 成人18禁高潮啪啪吃奶动态图| 亚洲五月色婷婷综合| 国产精品久久久人人做人人爽| 亚洲国产欧美一区二区综合| 91老司机精品| 亚洲美女搞黄在线观看| 国产精品av久久久久免费| 日本一区二区免费在线视频| 亚洲婷婷狠狠爱综合网| 男人添女人高潮全过程视频| 这个男人来自地球电影免费观看 | 夫妻性生交免费视频一级片| av网站在线播放免费| 女性被躁到高潮视频| 国产色婷婷99| 男的添女的下面高潮视频| 国产激情久久老熟女| 国产精品久久久久久精品电影小说| 日韩 欧美 亚洲 中文字幕| 美女大奶头黄色视频| 欧美国产精品一级二级三级| 日韩一区二区三区影片| 国产探花极品一区二区| 国产成人91sexporn| 女人爽到高潮嗷嗷叫在线视频| 91成人精品电影| 国产有黄有色有爽视频| 多毛熟女@视频| 韩国av在线不卡| 免费在线观看完整版高清| 亚洲色图 男人天堂 中文字幕| 看免费av毛片| 亚洲成人免费av在线播放| 你懂的网址亚洲精品在线观看| 成人国语在线视频| 国产片特级美女逼逼视频| 成人午夜精彩视频在线观看| 日本色播在线视频| 美国免费a级毛片| 国产一区有黄有色的免费视频| 交换朋友夫妻互换小说| 两个人看的免费小视频| 精品亚洲乱码少妇综合久久| 久久久亚洲精品成人影院| 亚洲精品国产色婷婷电影| 日韩一区二区三区影片| 老司机影院毛片| 最近最新中文字幕免费大全7| 久久久久网色| 亚洲在久久综合| 亚洲av福利一区| 欧美国产精品va在线观看不卡| 美女中出高潮动态图| 日韩大码丰满熟妇| 午夜91福利影院| 一个人免费看片子| 国产极品天堂在线| 青春草亚洲视频在线观看| 国产日韩欧美视频二区| 亚洲av日韩精品久久久久久密 | 国产亚洲最大av| 国产亚洲一区二区精品| 国产高清不卡午夜福利| 国产视频首页在线观看| 亚洲美女搞黄在线观看| 久久久亚洲精品成人影院| 人人妻人人添人人爽欧美一区卜| 丁香六月天网| 中文欧美无线码| 丝袜美足系列| 一区二区日韩欧美中文字幕| 人人妻人人澡人人看| 国语对白做爰xxxⅹ性视频网站| 国产精品香港三级国产av潘金莲 | 中文字幕最新亚洲高清| 极品少妇高潮喷水抽搐| 视频在线观看一区二区三区| 在线观看免费高清a一片| 91成人精品电影| 电影成人av| 看免费av毛片| h视频一区二区三区| 超色免费av| 久久99热这里只频精品6学生| 亚洲国产精品国产精品| 十八禁人妻一区二区| 啦啦啦中文免费视频观看日本| 丁香六月天网| 如日韩欧美国产精品一区二区三区| 王馨瑶露胸无遮挡在线观看| 国产又色又爽无遮挡免| 亚洲国产av新网站| 秋霞在线观看毛片| 欧美日韩一区二区视频在线观看视频在线| 一级毛片电影观看| 男女午夜视频在线观看| 嫩草影院入口| 欧美亚洲日本最大视频资源| a级片在线免费高清观看视频| 亚洲av中文av极速乱| 9色porny在线观看| 国产日韩欧美视频二区| 人人妻人人爽人人添夜夜欢视频| 观看av在线不卡| 久久青草综合色| www.自偷自拍.com| 男女之事视频高清在线观看 | 亚洲国产欧美一区二区综合| 欧美中文综合在线视频| 99热全是精品| 日本欧美视频一区| 国产免费福利视频在线观看| 人人澡人人妻人| 亚洲av在线观看美女高潮| 成人毛片60女人毛片免费| 精品一区二区三区av网在线观看 | 最近中文字幕高清免费大全6| 韩国av在线不卡| 天天影视国产精品| 深夜精品福利| 亚洲av日韩在线播放| 如日韩欧美国产精品一区二区三区| 亚洲国产中文字幕在线视频| 老司机深夜福利视频在线观看 | 亚洲av电影在线观看一区二区三区| 亚洲欧美清纯卡通| 悠悠久久av| 人妻一区二区av| 麻豆av在线久日| 国产精品免费大片| 国产老妇伦熟女老妇高清| 精品视频人人做人人爽| 人成视频在线观看免费观看| 国产成人av激情在线播放| 丝袜喷水一区| 大香蕉久久成人网| 久久久久网色| 午夜福利在线免费观看网站| 国产精品女同一区二区软件| 高清黄色对白视频在线免费看| 亚洲国产精品一区三区| 伊人久久大香线蕉亚洲五| 欧美黑人精品巨大| 欧美久久黑人一区二区| 日本av手机在线免费观看| 狂野欧美激情性bbbbbb| 亚洲精品一二三| 高清黄色对白视频在线免费看| 国产极品粉嫩免费观看在线| 日韩成人av中文字幕在线观看| 亚洲色图综合在线观看| 高清欧美精品videossex| 乱人伦中国视频| 波野结衣二区三区在线| 日本一区二区免费在线视频| 在线观看免费高清a一片| 天天躁夜夜躁狠狠躁躁| 黄色视频不卡| 在线观看免费午夜福利视频| 精品少妇内射三级| 久久久久久人妻| 女性生殖器流出的白浆| 少妇被粗大猛烈的视频| av电影中文网址| 精品一品国产午夜福利视频| 亚洲av国产av综合av卡| 中国国产av一级| 晚上一个人看的免费电影| 亚洲精品aⅴ在线观看| 久久97久久精品| 久久久精品免费免费高清| 精品一区二区免费观看| 韩国高清视频一区二区三区| 亚洲av欧美aⅴ国产| 亚洲婷婷狠狠爱综合网| 岛国毛片在线播放| 一区二区三区精品91| 日本爱情动作片www.在线观看| 国产淫语在线视频| 国产一区亚洲一区在线观看| 婷婷成人精品国产| 电影成人av| 欧美乱码精品一区二区三区| 亚洲久久久国产精品| av网站在线播放免费| 伊人久久大香线蕉亚洲五| 日韩大码丰满熟妇| 欧美精品亚洲一区二区| 男女之事视频高清在线观看 | 一级片'在线观看视频| 婷婷色av中文字幕| 亚洲av国产av综合av卡| 黄色视频在线播放观看不卡| 最近最新中文字幕大全免费视频 | 黄色 视频免费看| 国产成人精品无人区| 啦啦啦在线免费观看视频4| 久热这里只有精品99| 十分钟在线观看高清视频www| 飞空精品影院首页| 午夜日本视频在线| 亚洲欧美精品综合一区二区三区| 亚洲欧美一区二区三区黑人| 一本久久精品| 岛国毛片在线播放| 日本爱情动作片www.在线观看| 日韩欧美一区视频在线观看| 精品午夜福利在线看| 丝瓜视频免费看黄片| 免费高清在线观看视频在线观看| 欧美日韩av久久| 欧美日韩亚洲综合一区二区三区_| 啦啦啦中文免费视频观看日本| 少妇被粗大的猛进出69影院| 国产精品欧美亚洲77777| 在线观看人妻少妇| 国产一区二区在线观看av| 性色av一级| 天天操日日干夜夜撸| 美国免费a级毛片| 九草在线视频观看| 中文字幕最新亚洲高清| 国产淫语在线视频| av网站在线播放免费| 涩涩av久久男人的天堂| 99热国产这里只有精品6| 日韩中文字幕欧美一区二区 | 亚洲欧洲精品一区二区精品久久久 | 久久久久久久大尺度免费视频| 国产黄色免费在线视频| 97人妻天天添夜夜摸| 人体艺术视频欧美日本| 好男人视频免费观看在线| 看免费成人av毛片| 综合色丁香网| 另类亚洲欧美激情| 国产亚洲最大av| 男女无遮挡免费网站观看| 一级毛片 在线播放| 国产爽快片一区二区三区| 91国产中文字幕| a 毛片基地| 国产成人精品久久久久久| 欧美日韩国产mv在线观看视频| 大香蕉久久网| 亚洲av电影在线观看一区二区三区| 一区二区日韩欧美中文字幕| 免费av中文字幕在线| av在线老鸭窝| 国产精品秋霞免费鲁丝片| 亚洲精品国产区一区二| 免费观看a级毛片全部| 国产亚洲av高清不卡| 男女下面插进去视频免费观看| 一本大道久久a久久精品| 免费高清在线观看日韩| 99久久综合免费| 日韩,欧美,国产一区二区三区| 黑人欧美特级aaaaaa片| 尾随美女入室| 精品一区二区三区av网在线观看 | 国产精品三级大全| 午夜福利一区二区在线看| 一级片'在线观看视频| 国产精品久久久av美女十八| 一本色道久久久久久精品综合| 高清在线视频一区二区三区| 国产一区二区在线观看av| 成年av动漫网址| 在线观看免费午夜福利视频| 日韩制服丝袜自拍偷拍| 欧美日韩一级在线毛片| 九九爱精品视频在线观看| 日本色播在线视频| 久久久久久久国产电影| 亚洲国产av新网站| 丁香六月天网| 久久久久国产一级毛片高清牌| 少妇精品久久久久久久| 成年女人毛片免费观看观看9 | 亚洲成av片中文字幕在线观看| 久久久久人妻精品一区果冻| 国产一卡二卡三卡精品 | 国产精品久久久av美女十八| 亚洲成国产人片在线观看| 国产野战对白在线观看| 精品午夜福利在线看| 国产日韩欧美在线精品| 亚洲国产中文字幕在线视频| 午夜免费男女啪啪视频观看| av一本久久久久| 伊人久久国产一区二区| 中文天堂在线官网| 欧美日韩福利视频一区二区| 午夜福利网站1000一区二区三区| 亚洲情色 制服丝袜| 另类精品久久| 高清黄色对白视频在线免费看| av国产久精品久网站免费入址| 国产又色又爽无遮挡免| 99热国产这里只有精品6| 老司机深夜福利视频在线观看 | 中文字幕最新亚洲高清| 91成人精品电影| 电影成人av| 啦啦啦在线免费观看视频4| 亚洲欧美色中文字幕在线| 国产精品成人在线| 日本欧美国产在线视频| 中文字幕人妻丝袜一区二区 | 一级爰片在线观看| 国产精品无大码| 下体分泌物呈黄色| 九草在线视频观看| 国产精品免费大片| 亚洲国产精品国产精品| 18在线观看网站| 人妻人人澡人人爽人人| 久久av网站| 精品久久蜜臀av无| 女人被躁到高潮嗷嗷叫费观| a级毛片在线看网站| 一级毛片我不卡| 精品人妻熟女毛片av久久网站| 久久久国产一区二区| 99久久人妻综合| 一区二区三区精品91| 在线 av 中文字幕| 久久国产精品大桥未久av| 久久综合国产亚洲精品| 9191精品国产免费久久| 丰满少妇做爰视频| 成人免费观看视频高清| 搡老乐熟女国产| 亚洲成人免费av在线播放| 一区二区三区精品91| 国产在视频线精品| 国产精品.久久久| 高清不卡的av网站| 亚洲国产精品999| 国产精品免费大片| 国产日韩欧美亚洲二区| 亚洲精品自拍成人| 精品久久久久久电影网| a 毛片基地| av.在线天堂| 精品国产国语对白av| 看十八女毛片水多多多| 啦啦啦视频在线资源免费观看| 亚洲久久久国产精品| 国产精品女同一区二区软件| 亚洲精华国产精华液的使用体验| 久久久久国产一级毛片高清牌| tube8黄色片| 一级毛片 在线播放| 三上悠亚av全集在线观看| av免费观看日本| 国产精品一二三区在线看| 老鸭窝网址在线观看| 我要看黄色一级片免费的| 国产欧美亚洲国产| 久久精品久久精品一区二区三区| 亚洲精品一二三| 精品一区在线观看国产| 青春草视频在线免费观看| 日韩大片免费观看网站| 18禁观看日本| 免费观看a级毛片全部| 激情视频va一区二区三区| 午夜91福利影院| 男女午夜视频在线观看| 久久婷婷青草| 精品国产露脸久久av麻豆| 91精品伊人久久大香线蕉| 亚洲国产看品久久| 自拍欧美九色日韩亚洲蝌蚪91| 丝袜喷水一区| 一级毛片 在线播放| 黄色 视频免费看| 欧美国产精品一级二级三级| 国产男女超爽视频在线观看| e午夜精品久久久久久久| 国产毛片在线视频| 国产在视频线精品| 国产极品粉嫩免费观看在线| 老汉色∧v一级毛片| 日韩av不卡免费在线播放| 免费在线观看视频国产中文字幕亚洲 | 精品一区二区三区四区五区乱码 | 色婷婷久久久亚洲欧美| 精品国产露脸久久av麻豆| 精品国产国语对白av| 欧美最新免费一区二区三区| 亚洲,一卡二卡三卡| 国产99久久九九免费精品| 19禁男女啪啪无遮挡网站| 大码成人一级视频| 女人高潮潮喷娇喘18禁视频| 99久久精品国产亚洲精品| 人人妻,人人澡人人爽秒播 | 亚洲,欧美精品.| 搡老乐熟女国产| a级片在线免费高清观看视频| 国产精品免费大片| 欧美97在线视频| 中国三级夫妇交换| 狠狠精品人妻久久久久久综合| 一级,二级,三级黄色视频| 夜夜骑夜夜射夜夜干| 午夜福利在线免费观看网站| 亚洲 欧美一区二区三区| 考比视频在线观看| 日本猛色少妇xxxxx猛交久久| 精品国产超薄肉色丝袜足j| a级片在线免费高清观看视频| 久久久国产一区二区| 国产亚洲一区二区精品| 最近最新中文字幕大全免费视频 | 国产精品女同一区二区软件| 国产午夜精品一二区理论片| 精品午夜福利在线看| 国产精品久久久av美女十八| 免费在线观看黄色视频的| 久久久久网色| 岛国毛片在线播放| 综合色丁香网| 男女床上黄色一级片免费看| 一边亲一边摸免费视频| 亚洲第一av免费看| 菩萨蛮人人尽说江南好唐韦庄| 精品少妇一区二区三区视频日本电影 | 久久精品熟女亚洲av麻豆精品| 国产极品天堂在线| av国产久精品久网站免费入址| 亚洲成色77777| 亚洲精品久久午夜乱码| 赤兔流量卡办理| 国产亚洲午夜精品一区二区久久| 日韩一区二区视频免费看| 色综合欧美亚洲国产小说| 亚洲图色成人| 精品福利永久在线观看| av不卡在线播放| 免费高清在线观看日韩| 国语对白做爰xxxⅹ性视频网站| 性高湖久久久久久久久免费观看| 亚洲色图综合在线观看| 亚洲,欧美精品.| 可以免费在线观看a视频的电影网站 | 久热这里只有精品99| 成人三级做爰电影| 国产精品麻豆人妻色哟哟久久| 日韩大片免费观看网站| 男女午夜视频在线观看| 精品一区二区免费观看| 91成人精品电影| 午夜激情久久久久久久| 国产精品香港三级国产av潘金莲 | 午夜激情久久久久久久| 国产黄色免费在线视频| 国产日韩欧美视频二区| 亚洲av国产av综合av卡| 亚洲欧美激情在线| 国产一区二区三区av在线| 婷婷色av中文字幕| 如何舔出高潮| bbb黄色大片| 国产乱人偷精品视频| 国产成人欧美在线观看 | 天美传媒精品一区二区| 男人爽女人下面视频在线观看| 国产免费福利视频在线观看| 超碰成人久久| 国产精品久久久久久人妻精品电影 | 精品一区二区免费观看| 国产亚洲av片在线观看秒播厂| 爱豆传媒免费全集在线观看| 天天影视国产精品| 一边亲一边摸免费视频| 男女无遮挡免费网站观看| 午夜激情久久久久久久| 熟女av电影| 啦啦啦中文免费视频观看日本| 大陆偷拍与自拍| 亚洲欧美中文字幕日韩二区| 亚洲美女黄色视频免费看| 亚洲精品一区蜜桃| 2018国产大陆天天弄谢| 热99久久久久精品小说推荐| 黑人猛操日本美女一级片| 国产精品二区激情视频| 精品卡一卡二卡四卡免费| 丝袜美腿诱惑在线| 国产精品免费视频内射| 日韩大码丰满熟妇| 国产成人欧美| 美女主播在线视频| 久久人人97超碰香蕉20202| 国产一区二区在线观看av| 无遮挡黄片免费观看| 一级毛片 在线播放| 十八禁高潮呻吟视频| 亚洲av男天堂| 国产精品香港三级国产av潘金莲 | 国产精品 国内视频| 少妇 在线观看| 日本午夜av视频| 亚洲国产毛片av蜜桃av| 欧美精品av麻豆av| 国产成人精品福利久久| 女人爽到高潮嗷嗷叫在线视频| 叶爱在线成人免费视频播放| 最近手机中文字幕大全| 久久人妻熟女aⅴ| 欧美日韩精品网址| 亚洲精品久久久久久婷婷小说| 老汉色∧v一级毛片| 精品国产一区二区久久| 欧美精品一区二区免费开放| 亚洲精品国产av成人精品| 亚洲第一区二区三区不卡| 亚洲熟女毛片儿| 国产精品无大码| 在线观看www视频免费| 国产麻豆69| 一区二区日韩欧美中文字幕| 国产成人av激情在线播放| 叶爱在线成人免费视频播放| 欧美成人午夜精品| 国产精品秋霞免费鲁丝片| 久久久精品免费免费高清| 中文字幕av电影在线播放| 在线天堂最新版资源| 中文乱码字字幕精品一区二区三区| 欧美精品人与动牲交sv欧美| 国产一级毛片在线| 女人精品久久久久毛片| 午夜影院在线不卡| www.精华液| 国产探花极品一区二区| 91精品伊人久久大香线蕉| 免费女性裸体啪啪无遮挡网站| 亚洲精品第二区| 成人黄色视频免费在线看| 精品人妻一区二区三区麻豆| 久久久久视频综合| 国产午夜精品一二区理论片| 啦啦啦在线免费观看视频4|