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

    基于Python 設(shè)計的TEQC 數(shù)據(jù)質(zhì)量可視化分析軟件

    2021-08-07 05:48:46吳家杰王挺黃爾雙
    全球定位系統(tǒng) 2021年3期
    關(guān)鍵詞:序列圖檢核頻點(diǎn)

    吳家杰,王挺,黃爾雙

    (江西省煤田地質(zhì)局測繪大隊(duì),南昌 330001)

    0 引 言

    伴隨全球衛(wèi)星導(dǎo)航系統(tǒng)(GNSS)的快速發(fā)展與廣泛應(yīng)用,用戶對GNSS 定位服務(wù)質(zhì)量提出了更高的需求.但GNSS 定位服務(wù)的準(zhǔn)確性、可靠性與可用性會受到GNSS 觀測數(shù)據(jù)質(zhì)量好壞的影響[1-2].通常在GNSS 觀測數(shù)據(jù)使用前,需對觀測數(shù)據(jù)進(jìn)行質(zhì)量評估和預(yù)處理,以消除或減弱各種誤差因素的影響,進(jìn)而提高GNSS 定位精度[3].目前工程應(yīng)用上公認(rèn)最為常用的GNSS 數(shù)據(jù)預(yù)處理軟件,是由美國衛(wèi)星導(dǎo)航系統(tǒng)與地殼變形觀測研究大學(xué)聯(lián)合體(UNAVCO Facility)開發(fā)的TEQC 軟件[4-6].

    TEQC 軟件以質(zhì)量檢核功能為核心,包含數(shù)據(jù)格式轉(zhuǎn)換、編輯等功能,可對GNSS 觀測數(shù)據(jù)進(jìn)行多角度全方位的質(zhì)量評估與預(yù)處理.但TEQC 軟件的可視化功能較為薄弱,不利于對觀測數(shù)據(jù)質(zhì)量的直觀評估,需借助QCVIEW、CF2PS 和QC2SKY 等第三方繪圖軟件查看分析.上述繪圖工具已不再支持TEQC 2013-03-15 之后版本生成的compact3 結(jié)果文件[7].為實(shí)現(xiàn)對更新的TEQC 軟件處理結(jié)果可視化分析,許多學(xué)者進(jìn)行了相關(guān)研究.文獻(xiàn)[7]基于MATLAB GUI開發(fā)了一套適用于compact3 結(jié)果文件的可視化界面軟件(TEQC Plot-View).文獻(xiàn)[8]基于MATLAB 二次開發(fā)編制了teqcplot 部分程序模塊,實(shí)現(xiàn)了對TEQC的圖形顯示與可視化查詢.文獻(xiàn)[3]基于MATLAB 對TEQC 繪圖工具進(jìn)行了二次開發(fā),豐富了其圖形顯示功能,實(shí)現(xiàn)了可視化預(yù)處理分析.文獻(xiàn)[9]基于C#語言開發(fā)了TeqcChart(V1.0)軟件.文獻(xiàn)[10]利用C#開發(fā)基于TEQC 的GPS 數(shù)據(jù)預(yù)處理工具箱.上述TEQC質(zhì)量可視化腳本的改進(jìn),依托于MATLAB 和C#,但MATLAB 運(yùn)行效率低,C#開發(fā)效率低且不利于后期維護(hù)和更新.目前,流行的Python 語言具有免費(fèi)開源、開發(fā)效率高、代碼簡潔等特點(diǎn),在一定程度上彌補(bǔ)了該不足.鑒于此,本文基于Python 語言對TEQC質(zhì)量檢核模塊進(jìn)行可視化軟件設(shè)計,并結(jié)合實(shí)驗(yàn)數(shù)據(jù)對編寫的軟件性能進(jìn)行測試.

    1 TEQC 質(zhì)量檢核

    質(zhì)量檢核模塊是TEQC 軟件的核心功能之一,其質(zhì)量檢核原理主要是利用GNSS 觀測數(shù)據(jù)的偽距和載波相位信息進(jìn)行線性組合,實(shí)現(xiàn)對GNSS 數(shù)據(jù)的信噪比(SNR)、多路徑、電離層延遲和電離層延遲變化率等質(zhì)量檢核指標(biāo)的分析[11].用戶執(zhí)行TEQC 質(zhì)量檢核功能后會生成多個質(zhì)量檢核結(jié)果文件.質(zhì)量檢核結(jié)果文件類型包含質(zhì)量檢核摘要文件、基礎(chǔ)數(shù)據(jù)信息文件、質(zhì)量檢核指標(biāo)結(jié)果文件,其文件數(shù)量取決于衛(wèi)星信號頻點(diǎn)總數(shù).在眾多質(zhì)量檢核文件之中,質(zhì)量檢核摘要文件的內(nèi)容最為豐富,給出了觀測歷元、觀測值統(tǒng)計、觀測失鎖統(tǒng)計、觀測值刪除統(tǒng)計、周跳比和多路徑等數(shù)據(jù)質(zhì)量指標(biāo)值.因此,在進(jìn)行GNSS 數(shù)據(jù)質(zhì)量全面評估時,質(zhì)量檢核摘要文件是不可或缺的部分[12].

    TEQC 軟件實(shí)際上是一個可執(zhí)行文件(exe 文件),需在DOS 環(huán)境中設(shè)置不同參數(shù)指令,才可完成TEQC軟件質(zhì)量檢核功能.質(zhì)量檢核模式依據(jù)質(zhì)量檢核指令中是否添加廣播星歷文件分為完整模式和輕量模式[7],在廣播星歷缺省狀態(tài)下,輕量模式只會生成部分質(zhì)量檢核結(jié)果文件(不包含高度角和方位角等信息文件).一般情況下,完整模式輸出全部的質(zhì)量檢核結(jié)果文件,下文如不特別說明均采用完整模式解算.

    隨著TEQC 軟件的持續(xù)改進(jìn),不同版本的質(zhì)量檢核指令存有微小差異,如更新日期在2013-11-13以上的TEQC 版本,在DOS 環(huán)境中輸入teqc+qc+all+ssv?nav brdm1681.20L abmf1682_e.20O,默認(rèn)僅生成1 個質(zhì)量檢核摘要文件.因此,需完全輸出所有質(zhì)量檢核文件,應(yīng)在+qc 指令后添加+plot 參數(shù),相應(yīng)的質(zhì)量檢核指令如下:

    其中:+qc 為選擇質(zhì)量檢核;+plot 為輸出所有質(zhì)量檢核結(jié)果文件;+all 為添加所有觀測類型數(shù)據(jù);+ssv 為輸出每個衛(wèi)星的多路徑值;?nav 為輸入廣播星歷文件,多個廣播星歷文件可并行輸入.

    完整模式下生成的質(zhì)量檢核結(jié)果文件類型如表1 所示.

    表1 完整模式下質(zhì)量檢核結(jié)果文件類型

    長期以來,TEQC 軟件質(zhì)量檢核模塊只支持GPS 和GLONASS 雙系統(tǒng)數(shù)據(jù),對北斗衛(wèi)星導(dǎo)航系統(tǒng)(BDS)和GALIEO 觀測數(shù)據(jù)不兼容.但隨著TEQC軟件的持續(xù)優(yōu)化,TEQC 2019-02-12 版本已提供BDS和GALIEO 數(shù)據(jù)接口,可以進(jìn)行相應(yīng)系統(tǒng)數(shù)據(jù)的質(zhì)量檢核操作[13].

    2 基于Python 的TEQC質(zhì)量可視化設(shè)計

    鑒于TEQC 質(zhì)量檢核可視化功能相對欠缺,本文利用Python 語言對TEQC 質(zhì)量檢核結(jié)果文件進(jìn)行了界面化設(shè)計與封裝處理,編寫了一款TEQC 質(zhì)量檢核的可視化軟件TPP(TEQC Plot of Pyhone).該軟件包含數(shù)據(jù)導(dǎo)入、質(zhì)量檢核、參數(shù)選擇、繪圖和保存等5 個模塊.其中,繪圖模塊是核心,參數(shù)選擇模塊的參數(shù)選項(xiàng)可根據(jù)用戶需求靈活設(shè)置,涵蓋衛(wèi)星數(shù)目限定與衛(wèi)星選定,可隨機(jī)進(jìn)行組合,查看需求衛(wèi)星的數(shù)據(jù)質(zhì)量檢核指標(biāo)情況,以滿足精密單點(diǎn)定位(PPP)的數(shù)據(jù)質(zhì)量標(biāo)準(zhǔn).TPP 軟件結(jié)構(gòu)具體流程如圖1 所示.

    圖1 TPP 軟件結(jié)構(gòu)圖

    繪圖模塊作為TPP 軟件核心,其可依據(jù)質(zhì)量檢核模塊輸出的結(jié)果文件,分別繪制極坐標(biāo)天空圖、衛(wèi)星觀測序列圖和高度角時間序列圖、參數(shù)時間序列圖以及質(zhì)量檢核指標(biāo)參數(shù)時間序列圖(包括SNR、多路徑和電離層延遲等).但由于單個圖形間相互對比性較差,難以充分顯示不同質(zhì)量檢核指標(biāo)相關(guān)性.因此,TPP 將極坐標(biāo)天空圖、衛(wèi)星觀測序列圖和高度角時間序列圖作為底圖形式疊加SNR、多路徑和電離層延遲等質(zhì)量檢核信息,可進(jìn)一步分析高度角、方位角和觀測時段與檢核質(zhì)量指標(biāo)各因素間的相關(guān)關(guān)系;同時增設(shè)了質(zhì)量摘要信息顯示功能,通過質(zhì)量摘要信息與繪圖模塊的聯(lián)合分析,用戶能夠綜合判別每顆衛(wèi)星的數(shù)據(jù)質(zhì)量好壞,利于后期的數(shù)據(jù)預(yù)處理,剔除出部分質(zhì)量較差的衛(wèi)星觀測信息,從而提高基線解算的質(zhì)量.圖2 為TPP 主體界面.

    圖2 TPP 主體界面

    TPP 還可限制觀測時段,解析不同觀測時段下的衛(wèi)星數(shù)目與質(zhì)量檢核情況,將普遍質(zhì)量偏差的時段進(jìn)行切割處理.其次,TEQC 封裝于TPP 之中,通過獨(dú)立線程在后臺運(yùn)行,減免了用戶牢記TEQC 復(fù)雜檢核指令,減輕了底層數(shù)據(jù)處理重復(fù)性.

    3 TPP 軟件性能測試

    為測試TPP 軟件性能穩(wěn)定性和驗(yàn)證TEQC 針對GALILEO 系統(tǒng)的質(zhì)量分析情況,本文選取了MGEX(Multi-GNSS Experiment)基準(zhǔn)站ABMF 站2020 年6 月16 日全天的GAILIEO 衛(wèi)星觀測數(shù)據(jù),作為軟件性能測試的實(shí)驗(yàn)數(shù)據(jù).并通過TPP 軟件對GAILIEO衛(wèi)星E6 頻點(diǎn)上的觀測數(shù)據(jù)質(zhì)量檢核結(jié)果進(jìn)行可視化測試分析.

    3.1 質(zhì)量檢核摘要顯示測試

    TEQC 的質(zhì)量檢核模式會生成兩種用戶所需要的質(zhì)量檢核結(jié)果(質(zhì)量檢核摘要、質(zhì)量檢核指標(biāo)).質(zhì)量檢核摘要文件主要是記錄衛(wèi)星系統(tǒng)質(zhì)量評估匯總結(jié)果,而質(zhì)量檢核指標(biāo)結(jié)果文件更多地是記錄每顆衛(wèi)星的質(zhì)量檢核指標(biāo)在不同時段、方位角及高度角的詳細(xì)信息.通過兩者的相互補(bǔ)充,可實(shí)現(xiàn)對GNSS 觀測數(shù)據(jù)的總體質(zhì)量評估,質(zhì)量檢核摘要界面如圖3 所示.

    圖3 質(zhì)量檢核摘要顯示界面

    3.2 極坐標(biāo)天空圖測試

    極坐標(biāo)天空圖主要包含兩種衛(wèi)星位置信息因素:1)高度角;2)方位角.依據(jù)兩種因素變化趨勢,可給出衛(wèi)星信號的運(yùn)行軌跡.同時極坐標(biāo)天空圖還可添加質(zhì)量檢核指標(biāo)信息于衛(wèi)星軌跡之上,利于從空間角度解析質(zhì)量檢核指標(biāo)的變化情況.

    圖4 為GALILEO 衛(wèi)星的高度角與方位角的變化情況.由圖4 可知,ABMF 站的觀測數(shù)據(jù)總計接收到22 顆GALILEO 衛(wèi)星數(shù)據(jù),表明ABMF 站GALILEO信號較為良好.

    圖4 GALILEO 衛(wèi)星極坐標(biāo)天空圖

    圖5 給出了GALILEO 衛(wèi)星的E6 頻點(diǎn)上SNR 的變化情況.由圖5 可知,E6 頻點(diǎn)SNR 值基本穩(wěn)定在30~55 dB-Hz,表明ABMF 站的觀測位置條件良好,并未受到太多遮擋物干擾.此外,GALILEO 衛(wèi)星E6頻點(diǎn)低SNR 數(shù)值集中于高度角10°~20°,后續(xù)對GALILEO 觀測數(shù)據(jù)進(jìn)行預(yù)處理時,可通過限制高度角范圍,剔除質(zhì)量較差的觀測信息,以保證定位精度.

    圖5 GALILEO 衛(wèi)星E6 頻點(diǎn)上SNR 極坐標(biāo)天空圖

    3.3 衛(wèi)星觀測序列圖測試

    衛(wèi)星觀測序列圖是依據(jù)不同觀測時段下衛(wèi)星的可見性繪制而成,用戶可了解每顆衛(wèi)星信號的接收與結(jié)束時段信息.通過質(zhì)量檢核指標(biāo)的疊加顯示,質(zhì)量檢核信息被有序的分布在每顆可見衛(wèi)星上,利于從時間角度解析不同衛(wèi)星質(zhì)量檢核指標(biāo)在觀測時段的變化情況.

    圖6 為GALILEO 衛(wèi)星在2020 年6 月16 日的全天衛(wèi)星可見性情況.由圖6 可知,多數(shù)GALILEO 衛(wèi)星的觀測時段大于9 h,僅有少數(shù)可視衛(wèi)星E33、E19和E12 的觀測時段相對較短.

    圖6 GALILEO 衛(wèi)星觀測序列圖

    圖7 為GALILEO 衛(wèi)星 E6 頻點(diǎn)上多路徑值與可見性之間的關(guān)系情況.由圖7 可知多路徑值與可見性并未存在關(guān)聯(lián).但在接收機(jī)開始與停止接收衛(wèi)星信號的前后時段,多路徑值較大.因此,可利用TEQC 軟件的數(shù)據(jù)預(yù)處理功能,將多路徑效應(yīng)嚴(yán)重的時段進(jìn)行剔除,以保證GNSS 數(shù)據(jù)質(zhì)量.

    圖7 GALILEO 衛(wèi)星E6 頻點(diǎn)上多路徑觀測序列圖

    3.4 高度角時間序列圖測試

    高度角時間序列圖是反映的不同歷元下衛(wèi)星高度角的變化情況,用戶可從時空角度去判別質(zhì)量檢核指標(biāo)的變化情況.

    圖8 與圖9 中分別給出了GALILEO 衛(wèi)星高度角和E6 頻點(diǎn)上電離層延遲的變化情況,由圖8 和圖9可知,電離層延遲與高度角存有一定關(guān)聯(lián)性,電離層誤差是由電離層效應(yīng)引起的觀測值誤差,其存在一定正負(fù)值[14].伴隨高度角的增加,電離層延遲逐漸減小,反之亦然.

    圖8 GALILEO 衛(wèi)星高度角時間序列圖

    圖9 GALILEO 衛(wèi)星E6 頻點(diǎn)上電離層延遲隨高度角變化信息圖

    3.5 質(zhì)量檢核指標(biāo)參數(shù)時間序列圖測試

    指標(biāo)參數(shù)時間序列圖是質(zhì)量檢核指標(biāo)隨時間的變化關(guān)系,用于分析各項(xiàng)質(zhì)量指標(biāo)的變化情況.為驗(yàn)證單顆衛(wèi)星質(zhì)量可視化效果,結(jié)合圖6 本文選取了觀測時段較長的E03 衛(wèi)星的觀測數(shù)據(jù)進(jìn)行分析.

    圖10 為E03 衛(wèi)星在E6 頻點(diǎn)上質(zhì)量檢核指標(biāo)隨時間的變化關(guān)系.由圖10 可知,E03 共有04:00—14:00 與15:00—20:00 兩個觀測時段,其對應(yīng)的質(zhì)量檢核指標(biāo)波動存有差異.04:00—14:00 時段僅在開始與結(jié)束時出現(xiàn)較大波動,而15:00—20:00整時段波動性都較大.結(jié)合高度角信息可知高度角越大的區(qū)間,質(zhì)量檢核指標(biāo)波動較小.上述實(shí)驗(yàn)證實(shí)了SNR 與高度角成正比[15-16]、多路徑值與高度角成反比[17]以及電離層延遲與高度角成反比[18]等相關(guān)結(jié)論.

    圖10 E03 衛(wèi)星E6 頻點(diǎn)上質(zhì)量檢核指標(biāo)參數(shù)時間序列圖

    4 結(jié)束語

    針對TEQC 2013-03-15 以后版本生成的質(zhì)量檢核結(jié)果文件,難以采用現(xiàn)有的軟件對其可視化表達(dá),不利于準(zhǔn)確可靠的數(shù)據(jù)預(yù)處理這一問題.本文利用Python語言對TEQC 的質(zhì)量檢核模塊進(jìn)行了可視化設(shè)計,編寫了一款TPP 質(zhì)量可視化軟件,并采用實(shí)測數(shù)據(jù)進(jìn)行軟件性能測試,得出以下三點(diǎn)結(jié)論:

    1) TPP 質(zhì)量可視化軟件可對單顆衛(wèi)星以及多顆衛(wèi)星不同觀測時段下數(shù)據(jù)進(jìn)行質(zhì)量檢核結(jié)果的可視化分析,輸出的可視化圖形較為豐富,能較為直觀反映出該系統(tǒng)下每顆衛(wèi)星的質(zhì)量狀況;

    2) TPP 質(zhì)量可視化軟件設(shè)置了質(zhì)量檢核摘要和質(zhì)量檢核指標(biāo)的雙重判別模式,增強(qiáng)了用戶對于整體與局部衛(wèi)星觀測數(shù)據(jù)質(zhì)量的判斷;

    3) TPP 質(zhì)量可視化軟件結(jié)合TEQC 軟件的文件切割與合并、衛(wèi)星系統(tǒng)的選擇以及特定衛(wèi)星的禁用等功能,可提取有效的GNSS 觀測數(shù)據(jù)信息,為后續(xù)GNSS 高精度數(shù)據(jù)處理做準(zhǔn)備.

    猜你喜歡
    序列圖檢核頻點(diǎn)
    基于4G MR的大數(shù)據(jù)分析優(yōu)化EPSFB語音時延
    基于 ROADS 的面向場景業(yè)務(wù)架構(gòu)建模方法
    從“檢核”到“鷹架”
    基于SPSS序列法的商務(wù)談判實(shí)務(wù)課程混合教學(xué)模式實(shí)證研究
    物流科技(2021年10期)2021-05-12 08:41:06
    垂直荷載木結(jié)構(gòu)大跨屋頂設(shè)計
    應(yīng)用ETDFA生成CBTC聯(lián)鎖軟件形式化模型的方法
    思維游戲
    喜劇世界(2016年24期)2017-01-04 05:06:56
    檢核目錄法的研究與應(yīng)用—以書架設(shè)計為例
    福建省廈門第一中學(xué)黃建通老師:中學(xué)生創(chuàng)新思維課程引入“奧斯本檢核表技法”
    新校長(2016年8期)2016-01-10 06:43:49
    基于測量報告數(shù)據(jù)優(yōu)化CSFB頻點(diǎn)配置的方法
    亚洲伊人久久精品综合| av免费观看日本| 9热在线视频观看99| 日本与韩国留学比较| 午夜影院在线不卡| av国产久精品久网站免费入址| 免费观看a级毛片全部| 国产毛片在线视频| 亚洲人与动物交配视频| 亚洲精品色激情综合| 免费观看在线日韩| 午夜激情av网站| 少妇的丰满在线观看| 亚洲人成77777在线视频| 国产乱来视频区| 亚洲欧洲日产国产| 丝袜人妻中文字幕| 午夜视频国产福利| 国产黄色视频一区二区在线观看| 亚洲av电影在线观看一区二区三区| 日本黄大片高清| 久久久久久人人人人人| 99久久综合免费| 街头女战士在线观看网站| 久久久久久久久久成人| 90打野战视频偷拍视频| 看非洲黑人一级黄片| 在线精品无人区一区二区三| www.熟女人妻精品国产 | av不卡在线播放| 国产精品成人在线| 国产一区二区三区综合在线观看 | 国产精品一二三区在线看| 国产一区二区在线观看av| av视频免费观看在线观看| 久久精品国产a三级三级三级| 热re99久久国产66热| 国产成人精品无人区| 少妇被粗大的猛进出69影院 | 久久久久久久久久久久大奶| 在线观看一区二区三区激情| 97人妻天天添夜夜摸| 丰满迷人的少妇在线观看| 亚洲精品乱码久久久久久按摩| 久久亚洲国产成人精品v| 狂野欧美激情性bbbbbb| 亚洲,一卡二卡三卡| 观看av在线不卡| 18禁国产床啪视频网站| 国产一级毛片在线| 男男h啪啪无遮挡| 午夜av观看不卡| 国产成人精品在线电影| 欧美国产精品va在线观看不卡| 一级片'在线观看视频| 少妇高潮的动态图| 一级爰片在线观看| av不卡在线播放| 久久人人爽人人爽人人片va| 欧美日韩成人在线一区二区| 亚洲国产精品一区二区三区在线| 一区二区三区乱码不卡18| 最黄视频免费看| 涩涩av久久男人的天堂| 亚洲一级一片aⅴ在线观看| 亚洲一码二码三码区别大吗| 精品人妻熟女毛片av久久网站| 久久精品熟女亚洲av麻豆精品| 久久精品国产鲁丝片午夜精品| 伊人久久国产一区二区| 亚洲国产成人一精品久久久| 国产乱来视频区| 亚洲色图综合在线观看| 日本wwww免费看| 国产成人精品在线电影| 欧美亚洲日本最大视频资源| 亚洲精品色激情综合| 国产免费福利视频在线观看| 久久久久久久亚洲中文字幕| 亚洲av免费高清在线观看| 亚洲伊人久久精品综合| 国产av精品麻豆| 国产熟女欧美一区二区| 亚洲精品久久成人aⅴ小说| 欧美日韩av久久| av天堂久久9| 街头女战士在线观看网站| 国产伦理片在线播放av一区| 亚洲av综合色区一区| 一区二区三区乱码不卡18| 爱豆传媒免费全集在线观看| 日韩av不卡免费在线播放| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 91午夜精品亚洲一区二区三区| 国产精品久久久久成人av| 在线看a的网站| 黑人巨大精品欧美一区二区蜜桃 | 亚洲欧美日韩卡通动漫| 日产精品乱码卡一卡2卡三| 国产在线免费精品| 建设人人有责人人尽责人人享有的| 少妇猛男粗大的猛烈进出视频| 成人亚洲欧美一区二区av| 熟妇人妻不卡中文字幕| 国内精品宾馆在线| 久久精品久久精品一区二区三区| 久久久久视频综合| 免费观看av网站的网址| 国产毛片在线视频| 在线免费观看不下载黄p国产| 久久久久久伊人网av| 国产精品无大码| 久久久久久久久久久久大奶| 90打野战视频偷拍视频| 久久久精品免费免费高清| 纯流量卡能插随身wifi吗| av不卡在线播放| 久久99热这里只频精品6学生| 国产精品秋霞免费鲁丝片| 母亲3免费完整高清在线观看 | 三上悠亚av全集在线观看| av国产久精品久网站免费入址| 国产精品欧美亚洲77777| 婷婷色麻豆天堂久久| 国产成人精品在线电影| 水蜜桃什么品种好| 久久韩国三级中文字幕| 中国三级夫妇交换| 街头女战士在线观看网站| 国产一区二区激情短视频 | 黄网站色视频无遮挡免费观看| 欧美成人午夜精品| 国产高清三级在线| 午夜精品国产一区二区电影| 欧美日韩亚洲高清精品| 国产精品蜜桃在线观看| 国产欧美日韩一区二区三区在线| 中文字幕制服av| 欧美成人午夜免费资源| 老熟女久久久| 精品亚洲成国产av| 80岁老熟妇乱子伦牲交| 一级爰片在线观看| 一个人免费看片子| 亚洲欧美色中文字幕在线| 巨乳人妻的诱惑在线观看| 日韩在线高清观看一区二区三区| 男人添女人高潮全过程视频| 亚洲欧美色中文字幕在线| av在线app专区| 亚洲av在线观看美女高潮| 日韩视频在线欧美| 久久国产精品男人的天堂亚洲 | 国产精品秋霞免费鲁丝片| 亚洲精品色激情综合| 高清av免费在线| 天堂8中文在线网| 一级毛片我不卡| 精品久久久精品久久久| 夜夜骑夜夜射夜夜干| 亚洲精品国产av成人精品| 亚洲少妇的诱惑av| 啦啦啦啦在线视频资源| 国产爽快片一区二区三区| 91成人精品电影| 人妻 亚洲 视频| 岛国毛片在线播放| 亚洲欧美精品自产自拍| 大香蕉久久成人网| 26uuu在线亚洲综合色| 最黄视频免费看| 少妇的逼好多水| 在线观看美女被高潮喷水网站| 午夜91福利影院| 国产精品一区www在线观看| 国产精品久久久久成人av| 日韩 亚洲 欧美在线| 亚洲精品国产av蜜桃| 免费少妇av软件| videosex国产| 精品福利永久在线观看| 哪个播放器可以免费观看大片| 成人国产av品久久久| 91aial.com中文字幕在线观看| 最黄视频免费看| 色婷婷av一区二区三区视频| 亚洲国产av新网站| 桃花免费在线播放| 国产亚洲午夜精品一区二区久久| 久久综合国产亚洲精品| 日韩成人av中文字幕在线观看| 18禁在线无遮挡免费观看视频| 一级黄片播放器| 自线自在国产av| 午夜精品国产一区二区电影| 纯流量卡能插随身wifi吗| 女的被弄到高潮叫床怎么办| 午夜福利乱码中文字幕| 国产高清不卡午夜福利| 日韩av不卡免费在线播放| 久久99精品国语久久久| 国产一区二区在线观看日韩| 精品久久蜜臀av无| 人人妻人人添人人爽欧美一区卜| 国精品久久久久久国模美| 亚洲国产色片| 男人舔女人的私密视频| 十八禁高潮呻吟视频| 国产av一区二区精品久久| av有码第一页| 亚洲第一区二区三区不卡| av在线观看视频网站免费| 日韩不卡一区二区三区视频在线| 日韩三级伦理在线观看| 久久午夜综合久久蜜桃| 我要看黄色一级片免费的| 在线观看www视频免费| 国产一区有黄有色的免费视频| 免费在线观看完整版高清| 日韩人妻精品一区2区三区| 久久久久久伊人网av| 国产av国产精品国产| 日韩一区二区三区影片| 免费黄频网站在线观看国产| 人妻一区二区av| 大话2 男鬼变身卡| 美女主播在线视频| 亚洲精品国产色婷婷电影| 国产精品久久久久久av不卡| 妹子高潮喷水视频| 日韩精品免费视频一区二区三区 | 午夜日本视频在线| 精品午夜福利在线看| 国产免费视频播放在线视频| 夜夜骑夜夜射夜夜干| 另类精品久久| 欧美97在线视频| 最近中文字幕2019免费版| 国产国语露脸激情在线看| 国产激情久久老熟女| 一级毛片我不卡| 在线观看www视频免费| 国产又爽黄色视频| 午夜免费观看性视频| 成年人免费黄色播放视频| 寂寞人妻少妇视频99o| 亚洲熟女精品中文字幕| 国产精品三级大全| 精品国产国语对白av| 丰满饥渴人妻一区二区三| 精品亚洲成国产av| 亚洲精品美女久久久久99蜜臀 | 久久国产亚洲av麻豆专区| 99九九在线精品视频| 人妻系列 视频| 晚上一个人看的免费电影| tube8黄色片| 国产白丝娇喘喷水9色精品| 男的添女的下面高潮视频| 日韩制服骚丝袜av| 久久青草综合色| 亚洲精品第二区| 一区二区日韩欧美中文字幕 | 亚洲精品久久久久久婷婷小说| av免费在线看不卡| 久久97久久精品| 青春草国产在线视频| 久久狼人影院| 欧美日韩一区二区视频在线观看视频在线| 超色免费av| 777米奇影视久久| 91精品伊人久久大香线蕉| 婷婷成人精品国产| 日韩欧美一区视频在线观看| 在线看a的网站| 777米奇影视久久| 女性被躁到高潮视频| videos熟女内射| 欧美日本中文国产一区发布| 在线观看免费日韩欧美大片| 日本猛色少妇xxxxx猛交久久| 女性被躁到高潮视频| 青春草视频在线免费观看| 大陆偷拍与自拍| 免费在线观看黄色视频的| 国产视频首页在线观看| 97人妻天天添夜夜摸| 久久久国产欧美日韩av| 色网站视频免费| 日韩伦理黄色片| 一二三四中文在线观看免费高清| 飞空精品影院首页| 91精品三级在线观看| 18在线观看网站| 国产高清国产精品国产三级| 最近2019中文字幕mv第一页| 国产精品久久久av美女十八| 色94色欧美一区二区| 免费播放大片免费观看视频在线观看| 久久99热这里只频精品6学生| 国产亚洲一区二区精品| 亚洲精华国产精华液的使用体验| 日本黄大片高清| 26uuu在线亚洲综合色| 免费大片18禁| 性色av一级| 人人妻人人爽人人添夜夜欢视频| 久久久亚洲精品成人影院| 国产精品三级大全| 少妇人妻精品综合一区二区| 亚洲人成网站在线观看播放| 欧美精品人与动牲交sv欧美| 国产男人的电影天堂91| www日本在线高清视频| h视频一区二区三区| 欧美精品人与动牲交sv欧美| av一本久久久久| 麻豆乱淫一区二区| 人妻一区二区av| 七月丁香在线播放| 亚洲婷婷狠狠爱综合网| 中文字幕亚洲精品专区| 国产欧美亚洲国产| 亚洲国产av影院在线观看| 嫩草影院入口| 内地一区二区视频在线| 波野结衣二区三区在线| 精品久久久久久电影网| 午夜老司机福利剧场| 内地一区二区视频在线| 大香蕉久久网| 国产一区二区在线观看日韩| 亚洲av电影在线观看一区二区三区| 三级国产精品片| 蜜桃在线观看..| 汤姆久久久久久久影院中文字幕| 欧美日韩视频高清一区二区三区二| 熟妇人妻不卡中文字幕| 在线天堂最新版资源| 我的女老师完整版在线观看| 99九九在线精品视频| 又黄又爽又刺激的免费视频.| 一级毛片黄色毛片免费观看视频| 又粗又硬又长又爽又黄的视频| 2021少妇久久久久久久久久久| 精品一区二区三卡| 日本黄大片高清| 男女边吃奶边做爰视频| 熟女人妻精品中文字幕| 久久毛片免费看一区二区三区| 国产极品天堂在线| 久久影院123| 亚洲精品国产av蜜桃| 久久久久久久久久久免费av| 街头女战士在线观看网站| 国产综合精华液| 国产精品国产三级国产专区5o| 精品国产一区二区三区四区第35| 成年人午夜在线观看视频| 丁香六月天网| 亚洲第一区二区三区不卡| 国产片内射在线| 久久99热6这里只有精品| 国产精品.久久久| 久久久久久伊人网av| 亚洲婷婷狠狠爱综合网| 99久久人妻综合| 亚洲av日韩在线播放| 国产高清国产精品国产三级| 婷婷色av中文字幕| 亚洲精品美女久久av网站| 成人国产av品久久久| 亚洲精品国产色婷婷电影| 麻豆精品久久久久久蜜桃| 伦理电影大哥的女人| 成人国产av品久久久| 欧美精品亚洲一区二区| a级片在线免费高清观看视频| 美女主播在线视频| 亚洲精品国产av成人精品| 国产午夜精品一二区理论片| 制服诱惑二区| 国产成人a∨麻豆精品| 乱码一卡2卡4卡精品| 日韩人妻精品一区2区三区| 国产高清国产精品国产三级| 老熟女久久久| 男女下面插进去视频免费观看 | 久久狼人影院| 久久久久久久久久久久大奶| 国产综合精华液| 精品亚洲成a人片在线观看| 最近中文字幕2019免费版| 国国产精品蜜臀av免费| 在线观看www视频免费| 欧美最新免费一区二区三区| 又粗又硬又长又爽又黄的视频| 国产精品久久久久久久久免| 99精国产麻豆久久婷婷| 大片电影免费在线观看免费| 久久人人爽人人片av| 国产白丝娇喘喷水9色精品| 自拍欧美九色日韩亚洲蝌蚪91| 黄网站色视频无遮挡免费观看| 免费av不卡在线播放| 777米奇影视久久| 免费看av在线观看网站| 中文字幕另类日韩欧美亚洲嫩草| 精品一区在线观看国产| 免费黄网站久久成人精品| 男的添女的下面高潮视频| 天美传媒精品一区二区| 精品亚洲成国产av| 汤姆久久久久久久影院中文字幕| 国产成人免费无遮挡视频| 日韩精品免费视频一区二区三区 | 多毛熟女@视频| 自线自在国产av| 日韩av不卡免费在线播放| 亚洲情色 制服丝袜| 亚洲国产精品国产精品| 国产精品国产三级专区第一集| 精品卡一卡二卡四卡免费| 成年动漫av网址| 婷婷色综合大香蕉| 免费日韩欧美在线观看| 在线天堂中文资源库| 高清av免费在线| 中文字幕免费在线视频6| 在线看a的网站| 一区二区三区四区激情视频| 精品一品国产午夜福利视频| 日韩中字成人| 九色亚洲精品在线播放| 国产精品一区二区在线观看99| 免费av中文字幕在线| 丝瓜视频免费看黄片| 欧美日韩成人在线一区二区| 少妇猛男粗大的猛烈进出视频| 男女无遮挡免费网站观看| 免费在线观看黄色视频的| 水蜜桃什么品种好| 免费看光身美女| 国产国拍精品亚洲av在线观看| 亚洲精品国产av蜜桃| 午夜福利网站1000一区二区三区| 黄色毛片三级朝国网站| 国产av一区二区精品久久| 国产精品秋霞免费鲁丝片| 国产精品久久久久成人av| 妹子高潮喷水视频| 国产免费福利视频在线观看| 免费观看无遮挡的男女| 中国美白少妇内射xxxbb| 自拍欧美九色日韩亚洲蝌蚪91| 女的被弄到高潮叫床怎么办| 亚洲av国产av综合av卡| 精品国产国语对白av| 熟女av电影| 国产亚洲欧美精品永久| 熟妇人妻不卡中文字幕| 汤姆久久久久久久影院中文字幕| 婷婷色综合www| 精品酒店卫生间| 丰满迷人的少妇在线观看| 精品视频人人做人人爽| 亚洲av.av天堂| 晚上一个人看的免费电影| 22中文网久久字幕| 中文字幕精品免费在线观看视频 | 精品第一国产精品| 巨乳人妻的诱惑在线观看| 国产极品粉嫩免费观看在线| 麻豆精品久久久久久蜜桃| 精品福利永久在线观看| 国产1区2区3区精品| 人人妻人人爽人人添夜夜欢视频| 男女边吃奶边做爰视频| 免费少妇av软件| 女人精品久久久久毛片| 纯流量卡能插随身wifi吗| 国产黄色免费在线视频| 日本欧美国产在线视频| 国产精品一区二区在线观看99| 黄片播放在线免费| 久久毛片免费看一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 日韩制服丝袜自拍偷拍| 国产精品偷伦视频观看了| 在线 av 中文字幕| 免费av不卡在线播放| 久久97久久精品| 最近的中文字幕免费完整| 国产精品.久久久| 天堂8中文在线网| 久久这里有精品视频免费| 亚洲综合色惰| 中文字幕亚洲精品专区| 超色免费av| 国产黄色视频一区二区在线观看| 欧美激情极品国产一区二区三区 | 美女中出高潮动态图| 亚洲国产av影院在线观看| 热re99久久精品国产66热6| 久久久亚洲精品成人影院| 99久国产av精品国产电影| 国产男女超爽视频在线观看| 人人妻人人添人人爽欧美一区卜| 日韩精品免费视频一区二区三区 | 精品少妇内射三级| 亚洲精品美女久久av网站| 欧美性感艳星| 精品一区二区三卡| www.色视频.com| 狠狠精品人妻久久久久久综合| freevideosex欧美| 97人妻天天添夜夜摸| 免费观看av网站的网址| 卡戴珊不雅视频在线播放| 欧美+日韩+精品| 国产片内射在线| 成人综合一区亚洲| 又大又黄又爽视频免费| 少妇熟女欧美另类| 22中文网久久字幕| 久久久久网色| 成人亚洲欧美一区二区av| 丝袜美足系列| 久久久久久久久久久免费av| av天堂久久9| 久久婷婷青草| 黑人高潮一二区| 91成人精品电影| www.色视频.com| 97在线人人人人妻| 国产成人免费观看mmmm| 狠狠精品人妻久久久久久综合| 亚洲色图 男人天堂 中文字幕 | 亚洲少妇的诱惑av| 成人手机av| 免费看av在线观看网站| 最新的欧美精品一区二区| 亚洲人与动物交配视频| 久久免费观看电影| 美女国产高潮福利片在线看| 交换朋友夫妻互换小说| 国产片特级美女逼逼视频| 一级毛片黄色毛片免费观看视频| 午夜福利影视在线免费观看| 美女脱内裤让男人舔精品视频| 99九九在线精品视频| 精品酒店卫生间| 高清不卡的av网站| 大香蕉久久网| 日日摸夜夜添夜夜爱| av有码第一页| 中文字幕av电影在线播放| 黄色配什么色好看| 免费大片18禁| 伦精品一区二区三区| 欧美少妇被猛烈插入视频| 国产av一区二区精品久久| 国产成人91sexporn| 久久久a久久爽久久v久久| 国产成人免费观看mmmm| 高清在线视频一区二区三区| 观看av在线不卡| 国产无遮挡羞羞视频在线观看| 91国产中文字幕| 香蕉精品网在线| 美女主播在线视频| 18禁观看日本| 精品一区二区免费观看| 涩涩av久久男人的天堂| 亚洲精品乱久久久久久| 亚洲国产看品久久| 亚洲成国产人片在线观看| 午夜福利视频在线观看免费| 女人精品久久久久毛片| 国产一区有黄有色的免费视频| 国产男女内射视频| 成人午夜精彩视频在线观看| 亚洲欧洲精品一区二区精品久久久 | 夜夜骑夜夜射夜夜干| 国产高清不卡午夜福利| 亚洲综合精品二区| 制服丝袜香蕉在线| 午夜激情av网站| 1024视频免费在线观看| 欧美成人午夜精品| 只有这里有精品99| 美女视频免费永久观看网站| 国产极品粉嫩免费观看在线| 97在线人人人人妻| 黑人高潮一二区| 免费人成在线观看视频色| 国产精品欧美亚洲77777| 久久精品aⅴ一区二区三区四区 | 美女脱内裤让男人舔精品视频| 久久午夜福利片| 老司机影院毛片| 国产老妇伦熟女老妇高清| 这个男人来自地球电影免费观看 | 99香蕉大伊视频| 亚洲国产成人一精品久久久| 日韩精品有码人妻一区| 十八禁高潮呻吟视频| 亚洲熟女精品中文字幕| 最近的中文字幕免费完整| 午夜av观看不卡| 亚洲国产欧美日韩在线播放| 国产黄色免费在线视频| 97人妻天天添夜夜摸| 日本午夜av视频| 精品亚洲成a人片在线观看|