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

    一種IGrubbs-LWLR的區(qū)域高程異常擬合方法

    2023-10-12 07:48:44劉立龍蒙金龍何廣煥胡鵬程
    無(wú)線電工程 2023年10期
    關(guān)鍵詞:檢核置信區(qū)間訓(xùn)練樣本

    張 炎,劉立龍,蒙金龍,徐 勇,何廣煥,胡鵬程

    (1. 桂林理工大學(xué) 測(cè)繪地理信息學(xué)院,廣西 桂林 541006;2.廣西空間信息與測(cè)繪重點(diǎn)實(shí)驗(yàn)室,廣西 桂林 541006;3.廣西建設(shè)職業(yè)技術(shù)學(xué)院 市政與交通學(xué)院,廣西 南寧 530007)

    0 引言

    現(xiàn)今通過(guò)全球?qū)Ш叫l(wèi)星系統(tǒng)(Global Navigation Satellite System,GNSS)獲取地理信息數(shù)據(jù)已成為空間大地測(cè)量領(lǐng)域的一種重要手段,該技術(shù)所采集的高程值是以WGS-84參考橢球面為基準(zhǔn)的大地高[1-2]。實(shí)際工程中,常常采用以似大地水準(zhǔn)面為基準(zhǔn)的正常高,二者之間存在的差值,即為高程異常值[3]。外野采集的數(shù)據(jù)中往往存在一定的粗差,剔除粗差及選定合適的高程異常預(yù)測(cè)模型的擬合方法是在工程項(xiàng)目中完成測(cè)高工作的重要內(nèi)容[4]。對(duì)于觀測(cè)值中的粗差判別/剔除,目前常用的方法有格拉布斯(Grubbs)法則、t檢驗(yàn)法和穩(wěn)健估計(jì)法等,針對(duì)不同的樣本總量及所存異常值比例,選用不同的方法進(jìn)行粗差識(shí)別及剔除的效果也不盡相同[5-7]。韓紅超[8]利用Grubbs檢驗(yàn)法對(duì)某市的沉降監(jiān)測(cè)數(shù)據(jù)進(jìn)行粗差探測(cè)和剔除,證明了Grubbs法應(yīng)用于大量數(shù)據(jù)中的粗差探測(cè)的可行性,取得了良好的效果。在完成原始測(cè)量數(shù)據(jù)中的粗差剔除后,擬合方法的選擇也極為重要,常用的多面函數(shù)法、樣條曲線擬合法、支持向量機(jī)等在面對(duì)山峰突起且山體分布不規(guī)則的情況時(shí),所建高程異常擬合模型的效果不理想,其精度也難以滿足實(shí)際測(cè)繪生產(chǎn)的需求[9]。

    本文提出一種改進(jìn)格拉布斯法(Improved Grubbs,IGrubbs)結(jié)合局部加權(quán)線性回歸(Local Weighted Linear Regression,LWLR)組合算法來(lái)構(gòu)建區(qū)域高程異常擬合模型。在原Grubbs的基礎(chǔ)上,引入自適應(yīng)迭代,構(gòu)建新的評(píng)判指標(biāo)參數(shù),來(lái)提高Grubbs法則對(duì)測(cè)量數(shù)據(jù)中粗差的判別及剔除效率。通過(guò)LWLR法對(duì)預(yù)處理后的高程異常數(shù)據(jù)進(jìn)行擬合,完成精度較高的區(qū)域高程異常預(yù)測(cè)模型的建立,為以后的測(cè)高工作提供一定的參考。

    1 原理與方法

    1.1 Grubbs法則

    在樣本總體標(biāo)準(zhǔn)差未知的情況下,Grubbs通過(guò)對(duì)服從正態(tài)樣本或接近服從正態(tài)樣本的數(shù)據(jù)進(jìn)行異常判決[10-11]。計(jì)算公式如下:

    (1)

    (2)

    1.2 IGrubbs粗差剔除法

    在對(duì)測(cè)量數(shù)據(jù)中的粗差進(jìn)行判定及識(shí)別時(shí),若樣本方差偏大或偏小,將會(huì)存在對(duì)之后的樣本產(chǎn)生誤判或漏判的可能[12]。特別是在樣本中非零值估計(jì)偏差受異常值影響過(guò)大時(shí),會(huì)導(dǎo)致后續(xù)樣本總量偏低,發(fā)生錯(cuò)誤識(shí)別及判定的機(jī)率將會(huì)提高[13]。為了有效避免誤判及漏判的情況發(fā)生,對(duì)原Grubbs作出改進(jìn),具體如下:

    ① 引入自適應(yīng)迭代。以樣本長(zhǎng)度作為自適應(yīng)迭代參數(shù)進(jìn)行循環(huán)判定,從而確保樣本中所有數(shù)值得到充分檢驗(yàn)。

    ② 設(shè)定粗差剔除完成的指標(biāo)參數(shù)。在剔除粗差的過(guò)程中,以觀測(cè)值殘差與標(biāo)準(zhǔn)差的比值與Grubbs臨界值的差作為粗差剔除完成的指標(biāo)參數(shù)。在對(duì)測(cè)量數(shù)據(jù)的粗差識(shí)別及判定的過(guò)程中,當(dāng)樣本中所有數(shù)據(jù)均已低于對(duì)應(yīng)的Grubbs臨界值時(shí),表明當(dāng)前樣本中已不存在異常值,則退出自適應(yīng)迭代,不再進(jìn)行冗余運(yùn)算。指標(biāo)參數(shù)flag具體如下:

    (3)

    式中:數(shù)據(jù)組U為儲(chǔ)存所有樣本數(shù)據(jù)的動(dòng)態(tài)數(shù)組,Gi為最大偏離值,Gg為Grubbs臨界值與標(biāo)準(zhǔn)差的乘積,L為數(shù)據(jù)組U的樣本總量。

    當(dāng)原始測(cè)量數(shù)據(jù)中存在粗差值時(shí),IGrubbs對(duì)其進(jìn)行識(shí)別及判定的具體步驟如下:

    ① 數(shù)據(jù)初步處理:原始樣本數(shù)據(jù)為存在粗差的所有測(cè)量點(diǎn),構(gòu)成三維(平面坐標(biāo)及高程異常)數(shù)組U,獲取每個(gè)高程異常值所對(duì)應(yīng)的平面坐標(biāo)的位置索引,將U中高程異常值按從小到大的順序排列,隨后計(jì)算U中樣本長(zhǎng)度L0,L0的大小即為自適應(yīng)迭代參數(shù)控制最大迭代次數(shù)T,從而保證樣本中所有數(shù)值得到充分檢驗(yàn)。

    ② 分配IGrubbs臨界值區(qū)間:IGrubbs臨界值區(qū)間由選定的置信區(qū)間及樣本長(zhǎng)度L決定,當(dāng)置信水平α=95.00%時(shí)則選取對(duì)應(yīng)95.00%置信區(qū)間的IGrubbs臨界值進(jìn)行異常值判定,隨著置信水平值的變大,其對(duì)應(yīng)的異常判定嚴(yán)密程度也更為強(qiáng)烈。通常情況下,置信區(qū)間的選取應(yīng)為適中狀態(tài),過(guò)高或過(guò)低均難以達(dá)到最佳的粗差識(shí)別目的。根據(jù)前者的研究經(jīng)驗(yàn),實(shí)驗(yàn)時(shí)在對(duì)高程異常數(shù)據(jù)中的粗差識(shí)別過(guò)程中選擇置信水平為95.00%或97.50%時(shí)置信區(qū)間所對(duì)應(yīng)的IGrubbs臨界值來(lái)作為粗差判定的標(biāo)準(zhǔn)。

    圖1 IGrubb粗差剔除流程Fig.1 The process of gross error elimination based on IGrubbs method

    1.3 LWLR算法

    在對(duì)訓(xùn)練樣本數(shù)據(jù)進(jìn)行預(yù)處理后,利用LWLR法來(lái)構(gòu)建區(qū)域高程異常擬合模型。線性回歸法求的是具有小均方誤差的無(wú)偏估計(jì),回歸曲線受整體影響較大,易形成欠擬合現(xiàn)象[14-15]。與之相比,LWLR法的優(yōu)勢(shì)在于每一次迭代預(yù)測(cè)時(shí)都會(huì)更新權(quán)值,重新確定擬合系數(shù),從而使擬合面逐漸逼近真實(shí)值[16]。在著重考慮到局部點(diǎn)的同時(shí),也使局部以外的整體數(shù)據(jù)參與到擬合曲線中來(lái),更適用于非線性變化的高程異常數(shù)據(jù)的擬合及預(yù)測(cè)。具體內(nèi)容如下:

    若存在m個(gè)高程擬合點(diǎn),(xi,yi)為擬合點(diǎn)平面坐標(biāo),f(xi,yi)為曲面擬合函數(shù),vi為殘差,ξi為對(duì)應(yīng)的高程異常值,則已知點(diǎn)與高程異常值的關(guān)系可表示為:

    (4)

    (5)

    (6)

    通過(guò)測(cè)量獲取的高程異常數(shù)據(jù)集來(lái)構(gòu)建一定范圍內(nèi)連續(xù)變化的坐標(biāo)曲面函數(shù),若選擇的擬合點(diǎn)分布不均勻或未選擇到能夠代表地形起伏的特征點(diǎn),則最終形成的擬合面將與實(shí)際相差較大。因此,在式(6)的基礎(chǔ)上,引入距離定權(quán)的權(quán)重矩陣W,以擬合點(diǎn)與檢核點(diǎn)的距離作為評(píng)判二者相關(guān)程度的標(biāo)準(zhǔn),使擬合函數(shù)能夠較為準(zhǔn)確地反映出周圍地形起伏變化的狀態(tài),從而有效降低擬合點(diǎn)遠(yuǎn)離檢核點(diǎn)對(duì)擬合模型產(chǎn)生的負(fù)面效果[17]。加入權(quán)重后的公式如下:

    (7)

    式中:本文選擇高斯核函數(shù)的權(quán)重作為其權(quán)重系數(shù)。

    (8)

    式中:di為檢核點(diǎn)到擬合點(diǎn)的距離,(xi,yi)為擬合點(diǎn)平面坐標(biāo),(ej,fj)為檢核點(diǎn)平面坐標(biāo),g為檢核點(diǎn)數(shù),k為波長(zhǎng)。

    1.4 擬合模型的構(gòu)建

    在工程測(cè)高工作中,粗差往往是難以避免的,若未能得到及時(shí)處理,將會(huì)對(duì)后期項(xiàng)目中待測(cè)點(diǎn)獲取的高程精度產(chǎn)生嚴(yán)重影響。本文利用IGrubbs-LWLR組合法來(lái)完成粗差的剔除及待測(cè)點(diǎn)的高程預(yù)測(cè)工作,建立較高精度的區(qū)域高程異常擬合模型。具體過(guò)程如下:

    ① 對(duì)訓(xùn)練樣本數(shù)據(jù)進(jìn)行初步處理,讀取數(shù)據(jù)集。

    每名學(xué)生都必須參加出科考核,考核內(nèi)容包括病史詢問(wèn)、體格檢查、病歷書寫及技能操作。除病歷書寫外其他項(xiàng)均由總住院醫(yī)師負(fù)責(zé)考核,病歷書寫由總住院醫(yī)師協(xié)助主治醫(yī)師實(shí)施考核。在考核中,總住院醫(yī)師應(yīng)關(guān)注每名學(xué)生的成績(jī),給出客觀、真實(shí)的評(píng)價(jià),指出學(xué)生存在的問(wèn)題,并給予相應(yīng)指導(dǎo);同時(shí)對(duì)于學(xué)生做得好的方面,給予表?yè)P(yáng)和鼓勵(lì),真正做到以考促學(xué)。4重視實(shí)習(xí)生反饋意見,不斷改進(jìn)教學(xué)方法

    ② 記錄數(shù)據(jù)組所對(duì)應(yīng)位置索引并排序。選取置信區(qū)間,確定Grubbs臨界值區(qū)間。

    ④ 迭代開始,對(duì)待檢點(diǎn)進(jìn)行粗差識(shí)別,調(diào)用格拉布斯臨界值λ(a,n),判斷該點(diǎn)的Gi是否大于λ(a,n)與標(biāo)準(zhǔn)差E的乘積,若Gi大于λ(a,n),將該點(diǎn)視為可疑值并記錄,保存該可疑點(diǎn);若Gi小于λ(a,n),則表明該檢測(cè)點(diǎn)中未含有粗差值,不做記錄及保存。

    ⑤ 對(duì)第一個(gè)待測(cè)點(diǎn)判別完成后,則迭代進(jìn)入下一個(gè)待測(cè)點(diǎn)的粗差識(shí)別工作,返回至步驟④中進(jìn)行。直到迭代次數(shù)達(dá)到最大值時(shí),訓(xùn)練樣本數(shù)據(jù)中的粗差判定完畢,統(tǒng)計(jì)訓(xùn)練樣本中含有粗差的數(shù)值,并對(duì)其進(jìn)行剔除,整理獲取新的訓(xùn)練樣本集。

    ⑥ 對(duì)剔除粗差后的訓(xùn)練數(shù)據(jù)集作出歸一化處理,利用LWLR法來(lái)構(gòu)建區(qū)域高程異常擬合模型。

    ⑦ 輸入預(yù)處理后的訓(xùn)練數(shù)據(jù)與測(cè)試數(shù)據(jù),設(shè)置參數(shù)k初始值及最大迭代次數(shù)Tm。

    ⑩ 判斷當(dāng)前的均方誤差值是否滿足閾值或者達(dá)到最大迭代次數(shù),若滿足條件,迭代停止,記錄并保存當(dāng)前模型的相關(guān)系數(shù)及預(yù)測(cè)值,此時(shí)所構(gòu)建的區(qū)域高程異常擬合模型的精度最佳;若未滿足條件,更新參數(shù)k(k=k+0.05),重復(fù)步驟⑧、⑨,直至達(dá)到限定條件位置,輸出結(jié)果,迭代停止。

    2 實(shí)例分析

    實(shí)驗(yàn)數(shù)據(jù)選定為廣西東北小范圍山地區(qū)域內(nèi)的采集的74個(gè)GNSS水準(zhǔn)重合點(diǎn)(四等水準(zhǔn)測(cè)量要求,部分點(diǎn)含有粗差),用于后期區(qū)域高程異常擬合模型建立及精度檢驗(yàn)。研究區(qū)域內(nèi)點(diǎn)位分布較為均勻,2個(gè)GNSS水準(zhǔn)重合點(diǎn)間隔在2.5 km范圍內(nèi),總區(qū)域覆蓋面積約為300 km2。均勻選取其中的59個(gè)GNSS水準(zhǔn)重合點(diǎn)作為模型的訓(xùn)練樣本集,剩余的15個(gè)重合點(diǎn)作為擬合模型精度的檢核點(diǎn)。

    2.1 IGrubbs剔除粗差

    為確保擬合模型的精確度及可靠性,通過(guò)IGrubbs對(duì)參與模型構(gòu)建的高程異常值進(jìn)行預(yù)處理,剔除存在粗差的數(shù)據(jù)。同時(shí)選用傳統(tǒng)Grubbs作對(duì)比分析,設(shè)置不同的置信區(qū)間,對(duì)59個(gè)建模點(diǎn)進(jìn)行粗差探測(cè)/剔除,2種算法在α=99.5%的高置信水平下包含了低置信區(qū)間剔除的所有可疑點(diǎn)。最終,傳統(tǒng)Grubbs篩選出高程異常值分別為20.553、19.258、20.856 m的3個(gè)可疑點(diǎn),IGrubbs篩選出高程異常值分別為20.553、19.258、19.363、20.856、20.461 m的5個(gè)可疑點(diǎn)。具體的效果對(duì)比如表1所示。

    表1 剔除點(diǎn)的檢核統(tǒng)計(jì)表Tab.1 The check statistics of elimination points 單位:m

    為了檢驗(yàn)2種方法是否存在粗差誤判的情況,選用常規(guī)移去恢復(fù)法的二次曲面對(duì)疑似點(diǎn)周邊的數(shù)個(gè)GNSS水準(zhǔn)重合點(diǎn)進(jìn)行擬合,獲取疑似點(diǎn)的預(yù)測(cè)高程異常值。經(jīng)過(guò)對(duì)比分析發(fā)現(xiàn),34號(hào)點(diǎn)的真實(shí)值與預(yù)測(cè)值的差異較小,可認(rèn)定為誤判,對(duì)其進(jìn)行留存,其余4個(gè)疑似點(diǎn)存在的誤差均可作為粗差進(jìn)行剔除。二者最終的結(jié)果比較如表2所示。

    表2 IGrubbs與Grubbs粗差剔除效果對(duì)比Tab.2 Comparison of gross error elimination effect between IGrubbs and Grubbs

    在實(shí)驗(yàn)過(guò)程中,Grubbs在置信水平為90.0%時(shí)僅識(shí)別出1個(gè)誤差值,3個(gè)誤差值未識(shí)別出;相比于Grubbs,IGrubbs在置信水平為95.0%或97.5%時(shí)能夠較為準(zhǔn)確地識(shí)別出4個(gè)誤差值并進(jìn)行剔除。在另外2種相同情況下,誤差值的識(shí)別及判定效率均要優(yōu)于Grubbs。在利用IGrubbs對(duì)訓(xùn)練樣本數(shù)據(jù)進(jìn)行預(yù)處理后,其剔除粗差前后的高程異常數(shù)據(jù)變化如圖2所示。

    由圖2可知,在通過(guò)IGrubbs進(jìn)行粗差探測(cè)及剔除時(shí),發(fā)現(xiàn)參與后期建模的59個(gè)GNSS水準(zhǔn)重合點(diǎn)中存在4個(gè)偏離正常值的數(shù)值,其中存在的粗差與正常值的最小差值為0.27 m。若未及時(shí)對(duì)其進(jìn)行識(shí)別及剔除,將會(huì)嚴(yán)重影響后期所建擬合模型的預(yù)測(cè)精度。在對(duì)訓(xùn)練樣本數(shù)據(jù)進(jìn)行預(yù)處理后,測(cè)試重合點(diǎn)及檢核點(diǎn)分布如圖3所示。

    圖3 點(diǎn)位分布Fig.3 Distribution of points

    2.2 構(gòu)建擬合模型

    在完成訓(xùn)練樣本數(shù)據(jù)的粗差識(shí)別及剔除后,通過(guò)LWLR法來(lái)建立區(qū)域高程異常擬合模型。建模前,對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行歸一化處理,可提高算法的計(jì)算效率。將模型中的參數(shù)k初始值設(shè)置為0.5,最大迭代次數(shù)為70,通過(guò)迭代計(jì)算完成研究區(qū)域內(nèi)最適高程異常擬合模型的選擇,并預(yù)測(cè)出檢核點(diǎn)的高程異常值。在訓(xùn)練數(shù)據(jù)集及檢核點(diǎn)預(yù)測(cè)數(shù)據(jù)的基礎(chǔ)上,借助Matlab繪制出其相應(yīng)的三維效果圖(如圖4所示)。為近一步比較分析擬合模型的預(yù)測(cè)精度,分別利用傳統(tǒng)多面函數(shù)法、基于移去恢復(fù)法的二次曲面擬合來(lái)建立預(yù)測(cè)模型[19-20],對(duì)3種方法建模的預(yù)測(cè)結(jié)果進(jìn)行了統(tǒng)計(jì),其相應(yīng)的殘差值分布變化如圖5所示。

    圖4 高程異常擬合三維圖Fig.4 Three-dimensional graph of elevation anomaly fitting

    圖5 模型檢核點(diǎn)的殘差對(duì)比Fig.5 Residual comparison of model checkpoints

    由圖5可以看出,傳統(tǒng)多面函數(shù)法建模的檢核點(diǎn)殘差波動(dòng)范圍最為明顯,穩(wěn)定性相對(duì)較差;IGrubbs-LWLR所對(duì)應(yīng)的殘差值的最大值及最小值均要優(yōu)于其他2種方法建立的擬合模型,變化更為趨近于0,擬合效果更為理想。為更加直觀地了解各模型的擬合效果,選擇外符合精度(檢核點(diǎn)的均方根誤差)作為模型精度的判定指標(biāo),具體如表3所示。

    表3 不同擬合模型檢核點(diǎn)的預(yù)測(cè)結(jié)果Tab.3 Prediction results of different fitting model checkpoints

    由表3可知,多面函數(shù)法殘差變化為-3.3~3.2 cm,基于移去恢復(fù)法的二次曲面擬合殘差變化為-2.7~2.9 cm,IGrubbs-LWLR建模的擬合殘差變化為-1.7~2.5 cm。相比之下,IGrubbs-LWLR的擬合殘差變化區(qū)間較小,外符合精度為±1.6 cm,其擬合模型預(yù)測(cè)的精準(zhǔn)度更高。表明針對(duì)小范圍的地形起伏較為明顯的區(qū)域,選用IGrubbs-LWLR來(lái)建立高程異常擬合模型具備一定的時(shí)效性及穩(wěn)定性。

    3 結(jié)論

    為解決測(cè)量數(shù)據(jù)中粗差干擾及高程異常擬合模性構(gòu)建方法選擇的問(wèn)題,本文提出利用IGrubbs-LWLR來(lái)完成區(qū)域高程異常擬合模型的建立,結(jié)論如下:

    ① 在常規(guī)Grubbs法則的基礎(chǔ)上,引入自適應(yīng)迭代及判定粗差提出工作完成的指標(biāo)參數(shù),降低了其原先易發(fā)生漏判及誤判的概率,提高了該算法對(duì)樣本數(shù)劇中粗差識(shí)別及判定剔除的效率;

    ② 在對(duì)數(shù)據(jù)預(yù)處理后,選用LWLR法來(lái)構(gòu)建區(qū)域高程異常擬合預(yù)測(cè)模型,與傳統(tǒng)多面函數(shù)法、基于移去恢復(fù)法的二次曲面擬合法相比,其精度分別提高了38.5%、23.8%。表明IGrubbs-LWLR應(yīng)用于區(qū)域高程異常擬合模型的構(gòu)建具備一定的現(xiàn)實(shí)性。

    當(dāng)前實(shí)驗(yàn)區(qū)域研究范圍有限,地勢(shì)變化較為簡(jiǎn)單,在之后的研究中,會(huì)將該方法應(yīng)用于更為復(fù)雜多變的區(qū)域來(lái)構(gòu)建擬合模型,同時(shí)也會(huì)對(duì)模型中的參數(shù)近一步優(yōu)化提取,提高擬合模型的穩(wěn)定性及時(shí)效性,為今后的測(cè)高工作提供更高的參考價(jià)值。

    猜你喜歡
    檢核置信區(qū)間訓(xùn)練樣本
    從“檢核”到“鷹架”
    定數(shù)截尾場(chǎng)合三參數(shù)pareto分布參數(shù)的最優(yōu)置信區(qū)間
    p-范分布中參數(shù)的置信區(qū)間
    多個(gè)偏正態(tài)總體共同位置參數(shù)的Bootstrap置信區(qū)間
    垂直荷載木結(jié)構(gòu)大跨屋頂設(shè)計(jì)
    人工智能
    列車定位中置信區(qū)間的確定方法
    寬帶光譜成像系統(tǒng)最優(yōu)訓(xùn)練樣本選擇方法研究
    融合原始樣本和虛擬樣本的人臉識(shí)別算法
    基于稀疏重構(gòu)的機(jī)載雷達(dá)訓(xùn)練樣本挑選方法
    亚洲激情五月婷婷啪啪| 精品一区二区三区视频在线| 国产男人的电影天堂91| 久久热精品热| 久久久久国产精品人妻aⅴ院| 一级av片app| 久久天躁狠狠躁夜夜2o2o| 长腿黑丝高跟| 婷婷色综合大香蕉| 精品一区二区三区视频在线观看免费| 亚洲人成网站在线观看播放| 国产 一区 欧美 日韩| 久久久国产成人精品二区| 亚洲国产日韩欧美精品在线观看| 国产精品久久久久久av不卡| 18禁在线播放成人免费| 亚洲精品色激情综合| 菩萨蛮人人尽说江南好唐韦庄 | 在线观看免费视频日本深夜| 久久亚洲国产成人精品v| 国产高清三级在线| 免费黄网站久久成人精品| 欧美色欧美亚洲另类二区| 亚洲国产色片| 99久久无色码亚洲精品果冻| 最近视频中文字幕2019在线8| 欧美精品国产亚洲| 色吧在线观看| 香蕉av资源在线| 国产 一区精品| 中文字幕av成人在线电影| 在线观看免费视频日本深夜| 少妇被粗大猛烈的视频| 日韩欧美 国产精品| 中文亚洲av片在线观看爽| 日韩国内少妇激情av| 99热精品在线国产| 国产大屁股一区二区在线视频| 亚洲欧美精品综合久久99| 日韩亚洲欧美综合| 毛片一级片免费看久久久久| 舔av片在线| 国产一区亚洲一区在线观看| 九九久久精品国产亚洲av麻豆| 亚洲欧美日韩无卡精品| 国产男靠女视频免费网站| 国产高潮美女av| 国产高清视频在线播放一区| 我的老师免费观看完整版| 亚洲中文字幕一区二区三区有码在线看| 国语自产精品视频在线第100页| 国产精品亚洲美女久久久| 特级一级黄色大片| 国产伦在线观看视频一区| 国产精品电影一区二区三区| 欧美高清性xxxxhd video| 日本黄色片子视频| 日韩 亚洲 欧美在线| 国产黄色视频一区二区在线观看 | 人妻制服诱惑在线中文字幕| 中文资源天堂在线| 欧美日本视频| 国内精品久久久久精免费| 黄色欧美视频在线观看| 精品久久久噜噜| 精品日产1卡2卡| 精品久久久久久成人av| 午夜免费激情av| 高清毛片免费观看视频网站| 久久天躁狠狠躁夜夜2o2o| 18禁裸乳无遮挡免费网站照片| 久久韩国三级中文字幕| 一夜夜www| 成人午夜高清在线视频| 最近手机中文字幕大全| 最新在线观看一区二区三区| 国产亚洲欧美98| 亚洲美女搞黄在线观看 | www日本黄色视频网| 国产伦精品一区二区三区视频9| 久久精品夜夜夜夜夜久久蜜豆| 日韩欧美一区二区三区在线观看| 久久久久久久午夜电影| 国产精品久久视频播放| 男人的好看免费观看在线视频| 国产精品久久视频播放| 成人无遮挡网站| 在线免费十八禁| 啦啦啦啦在线视频资源| 成熟少妇高潮喷水视频| 99久久成人亚洲精品观看| 天堂影院成人在线观看| 亚洲在线观看片| 成人欧美大片| 国产精品不卡视频一区二区| 国产精品99久久久久久久久| 久久久久国产精品人妻aⅴ院| 在线观看66精品国产| 日本熟妇午夜| 精品午夜福利视频在线观看一区| 99久国产av精品| 亚洲国产精品国产精品| 久久久久免费精品人妻一区二区| 欧美激情久久久久久爽电影| 亚洲精品乱码久久久v下载方式| 两性午夜刺激爽爽歪歪视频在线观看| 久久亚洲精品不卡| 女人十人毛片免费观看3o分钟| 国产亚洲欧美98| 观看美女的网站| 99久久中文字幕三级久久日本| 国产国拍精品亚洲av在线观看| 国产成人福利小说| 亚洲精品国产成人久久av| 我要看日韩黄色一级片| 99国产精品一区二区蜜桃av| 99久久无色码亚洲精品果冻| 熟女人妻精品中文字幕| 蜜桃久久精品国产亚洲av| 国产成人一区二区在线| 天天躁夜夜躁狠狠久久av| 国产精品,欧美在线| 嫩草影院入口| 在线看三级毛片| 又爽又黄a免费视频| 少妇的逼水好多| 亚洲最大成人手机在线| aaaaa片日本免费| 中文字幕人妻熟人妻熟丝袜美| 内射极品少妇av片p| 日本免费一区二区三区高清不卡| 国产男靠女视频免费网站| 我的女老师完整版在线观看| 成人午夜高清在线视频| 久久久国产成人免费| 成年av动漫网址| a级一级毛片免费在线观看| av福利片在线观看| 国产三级在线视频| 天堂√8在线中文| 亚洲精品国产av成人精品 | 国产精品久久久久久亚洲av鲁大| 久久午夜福利片| 精品久久久噜噜| 日产精品乱码卡一卡2卡三| 日本a在线网址| 午夜久久久久精精品| 91久久精品电影网| 成人午夜高清在线视频| 99久久精品一区二区三区| 亚洲不卡免费看| 乱码一卡2卡4卡精品| 精品一区二区三区人妻视频| 成人漫画全彩无遮挡| 精品99又大又爽又粗少妇毛片| 成人特级黄色片久久久久久久| 国产一区二区激情短视频| 非洲黑人性xxxx精品又粗又长| 五月玫瑰六月丁香| 熟女人妻精品中文字幕| av福利片在线观看| 国产精品女同一区二区软件| 观看免费一级毛片| 97人妻精品一区二区三区麻豆| 男女视频在线观看网站免费| 免费av观看视频| 看非洲黑人一级黄片| 春色校园在线视频观看| 一个人看视频在线观看www免费| 99热网站在线观看| 国产淫片久久久久久久久| 搡老熟女国产l中国老女人| www日本黄色视频网| 午夜福利视频1000在线观看| 久久鲁丝午夜福利片| 久久久久久伊人网av| 成人性生交大片免费视频hd| 三级男女做爰猛烈吃奶摸视频| 亚洲精品粉嫩美女一区| 亚洲av不卡在线观看| 国产熟女欧美一区二区| 一个人免费在线观看电影| 99久久成人亚洲精品观看| 日韩人妻高清精品专区| 毛片一级片免费看久久久久| 麻豆一二三区av精品| 婷婷色综合大香蕉| 精品无人区乱码1区二区| 少妇高潮的动态图| 国产精品人妻久久久影院| 国产精品一二三区在线看| 亚洲av二区三区四区| 亚洲熟妇熟女久久| 欧美国产日韩亚洲一区| 国产精品无大码| 一级毛片我不卡| 日韩欧美在线乱码| 国产精品一区二区三区四区久久| 国产精品久久电影中文字幕| 欧美又色又爽又黄视频| 亚洲四区av| 欧美激情国产日韩精品一区| 国产精品三级大全| 最近最新中文字幕大全电影3| 亚洲精品日韩在线中文字幕 | avwww免费| 国产人妻一区二区三区在| 免费看光身美女| 波多野结衣高清作品| 国产单亲对白刺激| 看十八女毛片水多多多| 国产一区二区亚洲精品在线观看| 尤物成人国产欧美一区二区三区| 亚洲av中文av极速乱| 国模一区二区三区四区视频| 亚洲最大成人中文| 精品一区二区免费观看| 最新在线观看一区二区三区| 午夜福利高清视频| 俺也久久电影网| 少妇的逼水好多| 色综合亚洲欧美另类图片| 韩国av在线不卡| 熟女电影av网| 免费人成视频x8x8入口观看| 欧美+日韩+精品| 精品午夜福利在线看| 蜜桃亚洲精品一区二区三区| 悠悠久久av| 少妇猛男粗大的猛烈进出视频 | 亚洲国产高清在线一区二区三| av在线天堂中文字幕| 乱码一卡2卡4卡精品| 久久人人精品亚洲av| 久久久久九九精品影院| 成年版毛片免费区| 日本免费a在线| 亚洲欧美日韩高清在线视频| 菩萨蛮人人尽说江南好唐韦庄 | 麻豆国产97在线/欧美| 九九爱精品视频在线观看| 亚洲精品亚洲一区二区| 午夜老司机福利剧场| 免费av观看视频| 在线免费观看的www视频| 欧美丝袜亚洲另类| 欧美日韩综合久久久久久| 激情 狠狠 欧美| 亚洲综合色惰| 亚洲av不卡在线观看| 男女啪啪激烈高潮av片| 一本一本综合久久| 国产视频一区二区在线看| 俄罗斯特黄特色一大片| 99热这里只有精品一区| 2021天堂中文幕一二区在线观| 午夜免费激情av| 可以在线观看毛片的网站| 国产不卡一卡二| 国产精品一区www在线观看| 亚洲国产精品合色在线| 噜噜噜噜噜久久久久久91| 又黄又爽又刺激的免费视频.| 国产三级在线视频| 亚洲国产色片| 三级国产精品欧美在线观看| 91久久精品国产一区二区成人| 真人做人爱边吃奶动态| 久久精品国产清高在天天线| 最近手机中文字幕大全| 青春草视频在线免费观看| 久久热精品热| 在线观看美女被高潮喷水网站| 亚洲三级黄色毛片| 一进一出好大好爽视频| 免费看av在线观看网站| 亚洲激情五月婷婷啪啪| 午夜爱爱视频在线播放| 悠悠久久av| 一级毛片我不卡| 久久鲁丝午夜福利片| 亚洲熟妇中文字幕五十中出| 一进一出抽搐动态| 国产视频一区二区在线看| 亚洲七黄色美女视频| 人人妻人人看人人澡| 亚洲自偷自拍三级| 十八禁国产超污无遮挡网站| 99久久中文字幕三级久久日本| 老司机福利观看| 午夜福利在线在线| 午夜激情欧美在线| 国产乱人视频| 亚洲国产精品成人久久小说 | 我的女老师完整版在线观看| 亚洲国产精品国产精品| 九九久久精品国产亚洲av麻豆| 亚洲精品456在线播放app| 男人舔奶头视频| 老司机影院成人| 波野结衣二区三区在线| 亚洲婷婷狠狠爱综合网| 婷婷亚洲欧美| 免费人成视频x8x8入口观看| 国产精品国产高清国产av| av在线天堂中文字幕| 国产av麻豆久久久久久久| 国产成年人精品一区二区| 日本三级黄在线观看| 国产 一区精品| 韩国av在线不卡| 91久久精品电影网| 成人av在线播放网站| 国产黄片美女视频| 哪里可以看免费的av片| 女的被弄到高潮叫床怎么办| 我的老师免费观看完整版| 嫩草影院新地址| 亚洲天堂国产精品一区在线| 日日啪夜夜撸| 亚洲最大成人手机在线| 日本-黄色视频高清免费观看| 国产精品久久视频播放| 国产高清不卡午夜福利| 午夜福利在线观看吧| 日韩,欧美,国产一区二区三区 | 国产欧美日韩精品亚洲av| 国产综合懂色| 日韩欧美三级三区| 亚洲在线观看片| 亚洲av.av天堂| 亚洲精品乱码久久久v下载方式| 欧美日韩精品成人综合77777| 夜夜夜夜夜久久久久| 中文字幕av成人在线电影| 麻豆国产av国片精品| 亚洲国产高清在线一区二区三| 国产精品久久久久久av不卡| 国产又黄又爽又无遮挡在线| 国产乱人视频| 中文字幕精品亚洲无线码一区| 日韩欧美 国产精品| 老司机影院成人| 亚洲第一区二区三区不卡| 我要看日韩黄色一级片| av黄色大香蕉| 久久久久久伊人网av| 直男gayav资源| 一区福利在线观看| 亚洲精品乱码久久久v下载方式| 黄片wwwwww| 极品教师在线视频| 男人舔女人下体高潮全视频| 亚洲久久久久久中文字幕| 搡老岳熟女国产| 久久精品国产99精品国产亚洲性色| 99久久久亚洲精品蜜臀av| 在现免费观看毛片| av在线蜜桃| 午夜精品在线福利| 国产精品综合久久久久久久免费| 狂野欧美白嫩少妇大欣赏| 亚洲成人久久性| 亚洲欧美日韩高清专用| 国产高清视频在线播放一区| 午夜激情福利司机影院| 91久久精品电影网| 非洲黑人性xxxx精品又粗又长| 国产视频一区二区在线看| 最近中文字幕高清免费大全6| 精品少妇黑人巨大在线播放 | 麻豆国产97在线/欧美| 亚洲欧美成人综合另类久久久 | 三级毛片av免费| 日韩人妻高清精品专区| 午夜激情欧美在线| 久久人人精品亚洲av| 女人被狂操c到高潮| 亚洲自拍偷在线| 国产精品伦人一区二区| 少妇被粗大猛烈的视频| 欧美色欧美亚洲另类二区| 亚洲成人久久性| 成人性生交大片免费视频hd| av天堂在线播放| 99热全是精品| 免费看a级黄色片| 97超级碰碰碰精品色视频在线观看| 神马国产精品三级电影在线观看| 色视频www国产| 国产三级在线视频| 亚洲av不卡在线观看| 午夜a级毛片| 在线播放国产精品三级| 国产一区亚洲一区在线观看| 99久久精品一区二区三区| 国产精品三级大全| 尾随美女入室| 日本爱情动作片www.在线观看 | 日本精品一区二区三区蜜桃| 精品99又大又爽又粗少妇毛片| eeuss影院久久| 欧美极品一区二区三区四区| 久久精品夜夜夜夜夜久久蜜豆| 99热这里只有精品一区| 国产精品国产三级国产av玫瑰| 亚洲av成人精品一区久久| 99热只有精品国产| 免费一级毛片在线播放高清视频| 久久久精品欧美日韩精品| 亚洲欧美成人精品一区二区| 俄罗斯特黄特色一大片| 少妇高潮的动态图| 一本久久中文字幕| 日本免费一区二区三区高清不卡| 午夜激情欧美在线| 亚洲成人精品中文字幕电影| 亚洲精品乱码久久久v下载方式| 香蕉av资源在线| 国产在线精品亚洲第一网站| 久久久精品94久久精品| 特级一级黄色大片| 欧美xxxx黑人xx丫x性爽| 国产熟女欧美一区二区| 国产精品一区www在线观看| 欧美又色又爽又黄视频| 亚洲18禁久久av| 亚洲人成网站高清观看| 国产精品三级大全| 国产精品99久久久久久久久| 天天躁日日操中文字幕| 免费观看在线日韩| 亚洲成人久久爱视频| 欧美又色又爽又黄视频| 亚洲av电影不卡..在线观看| 日韩欧美三级三区| 亚洲av中文字字幕乱码综合| 男女之事视频高清在线观看| 日本成人三级电影网站| 最近在线观看免费完整版| 欧美+亚洲+日韩+国产| 欧美中文日本在线观看视频| 国产蜜桃级精品一区二区三区| 99热这里只有是精品50| 丰满乱子伦码专区| 日本熟妇午夜| 男插女下体视频免费在线播放| 两个人的视频大全免费| 最新在线观看一区二区三区| 久久99热这里只有精品18| 精品久久久久久久久久免费视频| 色综合站精品国产| 免费看光身美女| 小蜜桃在线观看免费完整版高清| 成人漫画全彩无遮挡| 久久鲁丝午夜福利片| 一边摸一边抽搐一进一小说| 三级国产精品欧美在线观看| 一级黄片播放器| 精品人妻视频免费看| 九九爱精品视频在线观看| 深夜精品福利| 久久久久免费精品人妻一区二区| 亚洲精品影视一区二区三区av| 亚洲欧美精品自产自拍| 国产日本99.免费观看| 看非洲黑人一级黄片| 丝袜喷水一区| 久久久精品欧美日韩精品| 欧美激情久久久久久爽电影| 99久国产av精品国产电影| 久久欧美精品欧美久久欧美| 午夜福利成人在线免费观看| 国产高清视频在线观看网站| 国产伦在线观看视频一区| 嫩草影院入口| 国产高潮美女av| 国产69精品久久久久777片| 男女做爰动态图高潮gif福利片| 久久久久久久亚洲中文字幕| av在线天堂中文字幕| h日本视频在线播放| 好男人在线观看高清免费视频| 亚洲精品亚洲一区二区| 99热这里只有精品一区| 国产国拍精品亚洲av在线观看| 色5月婷婷丁香| av卡一久久| 身体一侧抽搐| 国产69精品久久久久777片| 色综合亚洲欧美另类图片| av专区在线播放| 午夜亚洲福利在线播放| 在线观看免费视频日本深夜| 少妇熟女欧美另类| 久久久成人免费电影| 日韩一区二区视频免费看| 亚洲七黄色美女视频| 久久精品91蜜桃| aaaaa片日本免费| 最新中文字幕久久久久| 人妻丰满熟妇av一区二区三区| 久久综合国产亚洲精品| 你懂的网址亚洲精品在线观看 | 春色校园在线视频观看| 免费看a级黄色片| av黄色大香蕉| 听说在线观看完整版免费高清| 三级经典国产精品| 午夜爱爱视频在线播放| 男女下面进入的视频免费午夜| 国产成人a区在线观看| or卡值多少钱| 精品不卡国产一区二区三区| 韩国av在线不卡| 久久久成人免费电影| 国产在线男女| 国产精品久久视频播放| 国产精品久久久久久亚洲av鲁大| 久久精品国产鲁丝片午夜精品| 女生性感内裤真人,穿戴方法视频| 天天一区二区日本电影三级| 1024手机看黄色片| 久久国产乱子免费精品| 人妻丰满熟妇av一区二区三区| 嫩草影院新地址| 日本黄大片高清| 性欧美人与动物交配| 精品乱码久久久久久99久播| 国产一区二区三区在线臀色熟女| av福利片在线观看| 久久精品国产亚洲网站| 国产精品av视频在线免费观看| 真人做人爱边吃奶动态| 亚洲精品影视一区二区三区av| 日韩一本色道免费dvd| 一卡2卡三卡四卡精品乱码亚洲| 床上黄色一级片| 免费无遮挡裸体视频| 国产成人aa在线观看| 亚洲美女视频黄频| 久久久久久久久中文| 中国美女看黄片| av天堂在线播放| 日韩欧美国产在线观看| 最近视频中文字幕2019在线8| 国产视频一区二区在线看| 亚洲成a人片在线一区二区| 99久久中文字幕三级久久日本| av在线亚洲专区| 国产亚洲av嫩草精品影院| 亚洲国产精品sss在线观看| 国产色爽女视频免费观看| 日韩成人av中文字幕在线观看 | 成人鲁丝片一二三区免费| 久久精品国产清高在天天线| 伦精品一区二区三区| 久久午夜亚洲精品久久| 亚洲成人中文字幕在线播放| 联通29元200g的流量卡| 欧美激情在线99| 男女啪啪激烈高潮av片| av黄色大香蕉| 在线观看美女被高潮喷水网站| 不卡一级毛片| 成人亚洲精品av一区二区| 床上黄色一级片| 黄色日韩在线| 夜夜夜夜夜久久久久| 天堂动漫精品| 国产精品国产高清国产av| 51国产日韩欧美| 久久精品国产清高在天天线| 国产免费男女视频| 成人av在线播放网站| 久久精品国产亚洲av香蕉五月| 国产精品不卡视频一区二区| 国产亚洲av嫩草精品影院| 一级a爱片免费观看的视频| 99热精品在线国产| 男女视频在线观看网站免费| 日本五十路高清| 天堂网av新在线| 69av精品久久久久久| 91午夜精品亚洲一区二区三区| 日韩欧美国产在线观看| 精品免费久久久久久久清纯| 22中文网久久字幕| 天堂影院成人在线观看| 日本a在线网址| 男女啪啪激烈高潮av片| 久久久a久久爽久久v久久| 两个人视频免费观看高清| 91午夜精品亚洲一区二区三区| 国产综合懂色| 成人午夜高清在线视频| 亚洲精品日韩在线中文字幕 | 精品一区二区三区人妻视频| 国产精品久久电影中文字幕| 最近视频中文字幕2019在线8| 永久网站在线| 亚洲不卡免费看| 国产精品女同一区二区软件| 日韩欧美 国产精品| av黄色大香蕉| 69人妻影院| 亚洲精品亚洲一区二区| 国产精品一区二区免费欧美| 亚洲国产精品成人综合色| 亚洲av中文av极速乱| 淫妇啪啪啪对白视频| 久久韩国三级中文字幕| 色在线成人网| 亚洲av中文字字幕乱码综合| 日本一二三区视频观看|