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

    多波束測(cè)深數(shù)據(jù)處理的抗差最小二乘配置迭代解法

    2017-06-07 10:32:40王樂(lè)洋陳漢清
    測(cè)繪學(xué)報(bào) 2017年5期
    關(guān)鍵詞:抗差協(xié)方差波束

    王樂(lè)洋,陳漢清

    1. 東華理工大學(xué)測(cè)繪工程學(xué)院,江西 南昌 330013; 2. 流域生態(tài)與地理環(huán)境監(jiān)測(cè)國(guó)家測(cè)繪地理信息局重點(diǎn)實(shí)驗(yàn)室,江西 南昌 330013; 3. 江西省數(shù)字國(guó)土重點(diǎn)實(shí)驗(yàn)室,江西 南昌330013

    ?

    多波束測(cè)深數(shù)據(jù)處理的抗差最小二乘配置迭代解法

    王樂(lè)洋1,2,3,陳漢清1,2

    1. 東華理工大學(xué)測(cè)繪工程學(xué)院,江西 南昌 330013; 2. 流域生態(tài)與地理環(huán)境監(jiān)測(cè)國(guó)家測(cè)繪地理信息局重點(diǎn)實(shí)驗(yàn)室,江西 南昌 330013; 3. 江西省數(shù)字國(guó)土重點(diǎn)實(shí)驗(yàn)室,江西 南昌330013

    針對(duì)利用最小二乘配置處理多波束測(cè)深數(shù)據(jù),存在二次曲面數(shù)學(xué)模型通常無(wú)法精確表征海底地形的整體變化趨勢(shì)以及觀測(cè)數(shù)據(jù)存在粗差或異常點(diǎn)時(shí),常規(guī)方法給出的協(xié)方差函數(shù)不能精確表征其統(tǒng)計(jì)特性的問(wèn)題,本文提出了一種抗差最小二乘配置迭代解法。該方法首先進(jìn)行協(xié)方差函數(shù)和觀測(cè)值方差陣初始化,以多面函數(shù)擬合趨勢(shì)項(xiàng),然后應(yīng)用等價(jià)權(quán)抗差估計(jì)并通過(guò)迭代計(jì)算,最終給出穩(wěn)健的協(xié)方差函數(shù)參數(shù)解及最小二乘配置解。利用本文提出的方法及傳統(tǒng)的方法處理實(shí)測(cè)的多波束測(cè)深數(shù)據(jù),試驗(yàn)結(jié)果表明,相比于傳統(tǒng)的方法,本文提出的方法能夠較好地表征海底地形的整體變化趨勢(shì),一定程度上克服了多波束測(cè)深數(shù)據(jù)中粗差或異常點(diǎn)的影響。相比于傳統(tǒng)的抗差方法,本文方法更為有效地識(shí)別出測(cè)深數(shù)據(jù)中異常點(diǎn),推估效果較好,具有穩(wěn)健性。

    最小二乘配置;海底地形生成;協(xié)方差函數(shù);抗差;趨勢(shì)項(xiàng)

    近年來(lái),在海底地形測(cè)量技術(shù)中多波束測(cè)深系統(tǒng)成為海洋測(cè)量的主要設(shè)備,給研究區(qū)域的海底地形構(gòu)建帶來(lái)了豐富的數(shù)據(jù)資料[1]。高精度的海底地形地貌信息是保障海上船只航行及艦艇作戰(zhàn)訓(xùn)練安全的必備資料,也是海洋工程、資源開(kāi)發(fā)和海洋科學(xué)研究等多個(gè)應(yīng)用領(lǐng)域的重要基礎(chǔ)信息[2]。為了提高海底地形構(gòu)建的精度,必須采用合適的算法對(duì)多波束測(cè)深數(shù)據(jù)進(jìn)行處理,眾多學(xué)者作了大量的研究工作[3]。文獻(xiàn)[3]提出了抗差Kriging擬合方法,該方法在精確估計(jì)的同時(shí)能夠檢測(cè)測(cè)深數(shù)據(jù)的異常點(diǎn),但是該方法需要測(cè)深數(shù)據(jù)滿足平穩(wěn)假設(shè)或內(nèi)蘊(yùn)假設(shè),只能對(duì)平坦的海底數(shù)據(jù)進(jìn)行處理。

    最小二乘配置最初是由Trarup提出,是用來(lái)研究地球形狀和重力場(chǎng)的一種數(shù)學(xué)方法[4],隨后Moritz和Wolf對(duì)最小二乘配置作了進(jìn)一步研究[5-6]。如今,最小二乘配置已被廣泛應(yīng)用于地殼形變分析[7-9]、高程異常[10]、重力異常[11]、坐標(biāo)轉(zhuǎn)換[12]等領(lǐng)域。文獻(xiàn)[13]將最小二乘配置應(yīng)用于海道測(cè)量非定位點(diǎn)的估計(jì)中,利用區(qū)域深度場(chǎng)信息解算出非定位點(diǎn)的水深值,該方法方便精度評(píng)估,取得了不錯(cuò)的效果。最小二乘配置能夠顧及多波束測(cè)深數(shù)據(jù)具有趨勢(shì)性和隨機(jī)性的影響,但無(wú)法探測(cè)出多波束測(cè)深數(shù)據(jù)中的粗差或異常點(diǎn),針對(duì)于此問(wèn)題,文獻(xiàn)[14]提出了一種抗差最小二乘配置方法進(jìn)行處理,取得了良好的效果。文獻(xiàn)[15]對(duì)多波束測(cè)深數(shù)據(jù)處理的傳統(tǒng)抗差最小二乘配置方法和抗差Kriging擬合方法進(jìn)行了比較分析,通過(guò)試驗(yàn)驗(yàn)證了在多波束測(cè)深數(shù)據(jù)處理中,最小二乘配置在模型上優(yōu)于Kriging擬合方法。但現(xiàn)有的多波束測(cè)深數(shù)據(jù)處理的配置法均基于二次多項(xiàng)式擬合趨勢(shì)項(xiàng),一般情況下,特別是在海底地形起伏較大時(shí),二次多項(xiàng)式難以表征海底地形的整體變化趨勢(shì),甚至擬合的地形可能會(huì)出現(xiàn)“龍格現(xiàn)象”[16],從而導(dǎo)致后續(xù)的經(jīng)驗(yàn)協(xié)方差函數(shù)估計(jì)存在偏差;另一方面,由于海況環(huán)境、測(cè)量?jī)x器以及人為等因素的影響,多波束測(cè)深數(shù)據(jù)不可避免的存在粗差點(diǎn)或異常點(diǎn)的情況,迫切需要一種有效的抗差方法進(jìn)行處理。

    為了解決以上兩個(gè)問(wèn)題,提高多波束測(cè)深數(shù)據(jù)處理的精度,本文提出一種抗差最小二乘配置迭代解法,該方法以多面函數(shù)代替?zhèn)鹘y(tǒng)的二次多項(xiàng)式擬合海底地形的趨勢(shì)項(xiàng),應(yīng)用Huber權(quán)函數(shù)并通過(guò)迭代的方式求出經(jīng)驗(yàn)協(xié)方差函數(shù)參數(shù)估計(jì)和最小二乘配置解。最后利用本文方法、傳統(tǒng)的抗差最小二乘配置方法及最小二乘配置方法進(jìn)行實(shí)際的多波束測(cè)深數(shù)據(jù)處理,驗(yàn)證了本文方法的有效性與可行性。

    1 多波束測(cè)深數(shù)據(jù)處理的最小二乘配置方法

    最小二乘配置方法能夠同時(shí)顧及到海底地形的趨勢(shì)項(xiàng)及隨機(jī)項(xiàng),根據(jù)海底地形生成的特點(diǎn),建立多波束測(cè)深觀測(cè)模型[17-18]

    Z=GY+BX+Δ

    (1)

    (2)

    式中,DΔ為測(cè)深值方差陣,PΔ為測(cè)深值權(quán)陣。將式(1)改寫(xiě)成誤差方程形式

    (3)

    根據(jù)廣義最小二乘原則,構(gòu)造目標(biāo)函數(shù)

    (4)

    式中,PX為多波束測(cè)深已測(cè)點(diǎn)與未測(cè)點(diǎn)的權(quán)陣,由式(4)求自由極值,并通過(guò)整理,得

    (5)

    將式(3)代入式(5),通過(guò)法方程解算,得到多波束測(cè)深數(shù)據(jù)處理的最小二乘配置解[17]

    (6)

    式中,P=(DΔ+DS)-1;DS為測(cè)深值已測(cè)點(diǎn)信號(hào)協(xié)方差;DS′S為測(cè)深值未測(cè)點(diǎn)信號(hào)與已測(cè)點(diǎn)信號(hào)的協(xié)方差。其中,DS和DS′S由經(jīng)驗(yàn)協(xié)方差函數(shù)給出,常選定高斯函數(shù)進(jìn)行擬合[8-9,14]

    (7)

    (8)

    (9)

    式中,(xi,yi)為i點(diǎn)測(cè)深值的坐標(biāo),Y=[abc

    def]T為二次多項(xiàng)式模型待估參數(shù)。但一般情況下,特別是在海底地形起伏較大的區(qū)域,利用二次多項(xiàng)式擬合海底的整體連續(xù)變化趨勢(shì)可能會(huì)出現(xiàn)“龍格現(xiàn)象”[16],不能夠精確地表征海底地形的整體連續(xù)變化情況。本文針對(duì)此問(wèn)題,提出利用多面函數(shù)代替?zhèn)鹘y(tǒng)的二次多項(xiàng)式進(jìn)行海底地形整體變化趨勢(shì)的擬合。

    文獻(xiàn)[19—20]提出了多面函數(shù)模型來(lái)逼近不規(guī)則表面的思想,即任何的一個(gè)數(shù)學(xué)表面均可由一系列的數(shù)學(xué)表面的總和,以任意的精度進(jìn)行逼近。文獻(xiàn)[21]提出了基于多面函數(shù)的最小二乘配置改進(jìn)方法,該方法利用多面函數(shù)擬合地殼形變的趨勢(shì)項(xiàng),并得到了較好的效果。假設(shè)海底地形的整體趨勢(shì)是一個(gè)連續(xù)不規(guī)則變化的表面,基于多面函數(shù)的優(yōu)良特性,利用多面函數(shù)對(duì)海底地形的整體趨勢(shì)項(xiàng)進(jìn)行逼近[20]

    (10)

    式中,γi為待定系數(shù),R為核函數(shù),(xi,yi)為核函數(shù)的選定點(diǎn)坐標(biāo)。其中,核函數(shù)R的表達(dá)形式為

    (11)

    h=Rγ

    (12)

    以多面函數(shù)擬合海底地形的整體連續(xù)變化趨勢(shì),在最小二乘配置的框架下,把多面函數(shù)融入最小二乘配置中,此時(shí)基于多面函數(shù)的最小二乘配置函數(shù)模型為

    Z=Rγ+BX+Δ

    (13)

    式中,Rγ表示海底地形的整體連續(xù)變化部分;R為n×m維的核函數(shù)矩陣;γ為m×1維的待定系數(shù)向量;BX為海底地形不規(guī)則變化的隨機(jī)部分。同理,根據(jù)廣義最小二乘原理,解得基于多面函數(shù)的最小二乘配置解[21]

    (14)

    2 多波束測(cè)深數(shù)據(jù)處理的抗差最小二乘配置迭代解法

    2.1 觀測(cè)值初始權(quán)陣的確定

    多波束測(cè)深系統(tǒng)是現(xiàn)今海洋測(cè)深的主要數(shù)據(jù)采集設(shè)備,它獲取的數(shù)據(jù)覆蓋范圍廣,采樣密度高等特點(diǎn)而被廣泛應(yīng)用。但由于海洋測(cè)量相比于陸地測(cè)量存在較多的不穩(wěn)定因素,導(dǎo)致多波束測(cè)深數(shù)據(jù)存在粗差或異常點(diǎn)的情況較多,當(dāng)多波束測(cè)深數(shù)據(jù)含有粗差或異常點(diǎn)時(shí),給出的最小二乘配置解并非是一個(gè)最優(yōu)解。針對(duì)多波束測(cè)深數(shù)據(jù)含有粗差或異常點(diǎn)的問(wèn)題,在文獻(xiàn)[21—24]的基礎(chǔ)上,本文提出了抗差最小二乘配置方法處理含有粗差或異常點(diǎn)的多波束測(cè)深數(shù)據(jù)。該方法首先進(jìn)行觀測(cè)值權(quán)陣的初始化,具體權(quán)陣初始化步驟如下:

    (1) 對(duì)多波束測(cè)深數(shù)據(jù)消除整體連續(xù)變化趨勢(shì),求取初始?xì)埐钪?/p>

    (15)

    (2) 根據(jù)給出的殘差初始值,計(jì)算初始方差因子,利用該式去計(jì)算參數(shù)的抗差解具有50%的高崩潰污染率

    (16)

    (3) 選取Huber權(quán)函數(shù)計(jì)算權(quán)因子[14]

    (17)

    2.2 協(xié)方差函數(shù)的抗差擬合方法

    經(jīng)驗(yàn)協(xié)方差函數(shù)擬合的準(zhǔn)確性是求解最小二乘配置問(wèn)題的關(guān)鍵。經(jīng)驗(yàn)協(xié)方差函數(shù)的擬合主要為求解協(xié)方差函數(shù)(式(7))的兩個(gè)參數(shù)解D(0)和k。由上述得到的更新后的觀測(cè)值權(quán)陣,代入式(8),得到抗差后的D(0)解

    (18)

    (19)

    (20)

    (1) 首先對(duì)式(20)左邊的k2開(kāi)根號(hào)處理,并取其中位數(shù),計(jì)算其殘差值,即

    (21)

    (2) 根據(jù)給出的殘差初始值,計(jì)算初始方差因子

    (22)

    (3) 選取Huber權(quán)函數(shù)計(jì)算權(quán)因子

    (23)

    最后通過(guò)迭代計(jì)算,給出穩(wěn)健的協(xié)方差函數(shù)參數(shù)k的值。抗差最小二乘配置迭代解法的計(jì)算步驟如下:

    (1) 初始化協(xié)方差函數(shù),給出協(xié)方差函數(shù)的初始值D(0)(0)=1,k(0)=0。利用觀測(cè)值中位數(shù),初始觀測(cè)值方差陣,以多面函數(shù)擬合趨勢(shì)項(xiàng),計(jì)算得到最小二乘配置的初始解;

    (2) 由步驟(1)得到的初始最小二乘配置解,去除趨勢(shì)項(xiàng),根據(jù)Huber權(quán)函數(shù)更新觀測(cè)值方差陣,隨后進(jìn)行信號(hào)協(xié)方差統(tǒng)計(jì),根據(jù)3.2節(jié)的方法得到協(xié)方差函數(shù)的參數(shù)解D(0)和k;

    (3) 由步驟(2)得到的新的協(xié)方差函數(shù)以及更新后的觀測(cè)值方差陣,計(jì)算信號(hào)協(xié)方差Ds,把Ds和更新后的觀測(cè)值方差陣代入式(14),計(jì)算得到新的最小二乘配置解。

    3 算例分析

    本文數(shù)據(jù)來(lái)自于中國(guó)某海域,水深范圍為1100~1400 m,海底地形起伏變化較大。本文選取其中一個(gè)條帶作為試驗(yàn)數(shù)據(jù),由于數(shù)據(jù)量較大,冗余點(diǎn)較多,首先進(jìn)行了數(shù)據(jù)抽稀處理,共選取886個(gè)多波束測(cè)深數(shù)據(jù)點(diǎn),其中隨機(jī)抽取不含粗差的150個(gè)點(diǎn)作為檢核點(diǎn),所有數(shù)據(jù)如圖1所示,星點(diǎn)為檢核點(diǎn),其余為已測(cè)點(diǎn)。測(cè)深值等高線圖如圖2所示。

    圖1 多波束測(cè)深數(shù)據(jù)分布圖Fig.1 Multi-beam data distribution map

    圖2 原始數(shù)據(jù)測(cè)深等高線圖Fig.2 Contour map of raw data

    為了驗(yàn)證本文算法的有效性,分別利用本文算法、傳統(tǒng)抗差最小二乘配置方法、基于多面函數(shù)的最小二乘配置方法和基于二次多項(xiàng)式的最小二乘配置方法進(jìn)行多波束測(cè)深數(shù)據(jù)處理,設(shè)計(jì)的4種方案如下所示:

    方案1:基于二次多項(xiàng)式的最小二乘配置法,該法以二次多項(xiàng)式擬合趨勢(shì)項(xiàng),計(jì)算檢核點(diǎn)的水深值和外符精度;

    方案2:基于二次多項(xiàng)式的傳統(tǒng)抗差最小二乘配置法,計(jì)算檢核點(diǎn)的水深值和外符精度;

    方案3:基于多面函數(shù)的最小二乘配置法,該法以多面函數(shù)擬合趨勢(shì)項(xiàng),計(jì)算檢核點(diǎn)的水深值和外符精度;

    方案4:基于多面函數(shù)的抗差最小二乘配置迭代解法,該法以多面的函數(shù)擬合趨勢(shì)面,結(jié)合抗差的方法進(jìn)行協(xié)方差函數(shù)擬合,計(jì)算檢核點(diǎn)的水深值和外符精度。

    通過(guò)以上4種方案對(duì)實(shí)際多波束測(cè)深數(shù)據(jù)進(jìn)行處理,4種方案的外符精度如表1所示,殘差統(tǒng)計(jì)如圖3所示。

    表1 4種方案檢核點(diǎn)的殘差統(tǒng)計(jì)情況

    Tab.1 Residual statistics of external check points using four schemes m

    通過(guò)表1可以看出,基于多面函數(shù)的最小二乘配置方法的各項(xiàng)指標(biāo)均比傳統(tǒng)的基于二次多項(xiàng)式的最小二乘配置要優(yōu)。在外符精度上,基于多面函數(shù)的最小二乘配置(方案3)和基于二次多項(xiàng)式的最小二乘配置方法(方案1)的外符精度分別為0.935 86 m和1.205 97 m,相比于傳統(tǒng)的方法,基于多面函數(shù)的最小二乘配置方法在外符精度上提高了0.270 11 m,說(shuō)明了以多面函數(shù)擬合海底地形的趨勢(shì)項(xiàng)更能反映研究區(qū)域海底地形的整體連續(xù)變化的趨勢(shì),去除趨勢(shì)項(xiàng)后得到的隨機(jī)變化部分也更為平穩(wěn),因此可以利用該方法得到更為正確的信號(hào)值。通過(guò)比較兩種方法的殘差統(tǒng)計(jì)圖(圖3),可以看出基于多面函數(shù)的最小二乘配置方法給出的殘差小于傳統(tǒng)的基于二次多項(xiàng)式的最小二乘配置方法的殘差,進(jìn)一步驗(yàn)證多面函數(shù)擬合趨勢(shì)項(xiàng)的有效性。另一方面,通過(guò)比較分析傳統(tǒng)的抗差最小二乘配置方法(方案2)與基于多面函數(shù)的最小二乘配置方法(方案3)可以發(fā)現(xiàn),雖然傳統(tǒng)的抗差方法進(jìn)行了抗差處理,但是由于二次多項(xiàng)式不能較好地表征海底地形的整個(gè)連續(xù)變化趨勢(shì),使得后續(xù)由扣除趨勢(shì)項(xiàng)后得到的殘差進(jìn)行擬合得到的協(xié)方差函數(shù)發(fā)生了偏差,導(dǎo)致推估的效果比沒(méi)有考慮到抗差的基于多面函數(shù)的最小二乘配置方法要差。最后,本文抗差最小二乘配置方法的推估效果最好,在外符精度上,比基于二次多項(xiàng)式的最小二乘配置方法提高了0.308 13 m,比傳統(tǒng)的抗差最小二乘配置方法提高了0.317 94 m,比基于多面函數(shù)的最小二乘配置方法提高了0.038 02 m,說(shuō)明了本文抗差方法能在一定程度上克服多波束測(cè)深數(shù)據(jù)中粗差或異常點(diǎn)對(duì)最小二乘配置結(jié)果的影響;同時(shí),相比傳統(tǒng)的抗差方法,本文方法推估精度更高。同樣的,為了進(jìn)一步驗(yàn)證本文抗差方法的有效性的,給出四種方案的經(jīng)驗(yàn)協(xié)方差函數(shù)的參數(shù)解,具體如表2所示。本文抗差最小二乘配置能夠在一定程度上抵抗多波束測(cè)深數(shù)據(jù)中粗差或異常點(diǎn)的影響,給出協(xié)方差函數(shù)更符合實(shí)際情況,更為客觀的反映點(diǎn)位之間的相關(guān)性。

    表2 4種方案的協(xié)方差函數(shù)參數(shù)解

    Tab.2 Fitting results of covariance function parameters using four schemes

    參數(shù)方案1方案2方案3方案4D(0)13.9217813.006262.440060.93110k26.95240×10-86.87077×10-85.97079×10-83.35709×10-7

    圖3 4種方案檢核點(diǎn)殘差統(tǒng)計(jì)圖Fig.3 Residual statistics chart of external check points using four schemes

    圖4 原始多波束測(cè)深數(shù)據(jù)TIN圖Fig.4 TIN chart of original multi-beam data

    圖5 傳統(tǒng)抗差方法處理后的多波束測(cè)深數(shù)據(jù)TIN圖Fig.5 TIN chart of multi-beam data using traditional robust least squares collocation

    圖6 經(jīng)本文抗差方法處理后的TIN圖Fig.6 TIN chart of multi-beam data using robust least squares collocation was proposed in this paper

    4 結(jié) 論

    為了進(jìn)一步提高多波束測(cè)深數(shù)據(jù)處理的精度,本文提出了一種抗差最小二乘配置迭代解法,給出了具體的求解公式以及迭代算法,通過(guò)算例驗(yàn)證了本文方法的有效性,得出了以下的結(jié)論:

    (1) 針對(duì)現(xiàn)有配置法的趨勢(shì)項(xiàng)的二次曲面數(shù)學(xué)模型通常無(wú)法精確表征海底地形的整體連續(xù)變化趨勢(shì)的問(wèn)題,本文提出以多面函數(shù)代替?zhèn)鹘y(tǒng)的方法擬合海底地形的趨勢(shì)項(xiàng),通過(guò)算例驗(yàn)證了該方法的有效性,進(jìn)一步提高了配置法在多波束測(cè)深數(shù)據(jù)處理中的推估精度。

    圖7 兩種抗差方法測(cè)深點(diǎn)等價(jià)權(quán)值概況Fig.7 Equivalent weight of sounding points for two robust methods

    (2) 針對(duì)多波束測(cè)深數(shù)據(jù)存在粗差或異常點(diǎn)時(shí),常規(guī)最小二乘配置方法給出的協(xié)方差函數(shù)不能精確表征其統(tǒng)計(jì)特性的問(wèn)題,提出了一種抗差最小二乘配置迭代解法,以迭代計(jì)算的方式給出了經(jīng)驗(yàn)協(xié)方差函數(shù)的參數(shù)估計(jì)和最小二乘配置解,通過(guò)算例驗(yàn)證了本文方法在一定程度上能夠克服多波束測(cè)深數(shù)據(jù)中粗差或異常點(diǎn)的影響,提高了推估的精度,相比于傳統(tǒng)的抗差方法,本文抗差方法能夠更為有效地識(shí)別出測(cè)深數(shù)據(jù)中的異常點(diǎn)。

    最小二乘配置問(wèn)題的關(guān)鍵是確定合理的協(xié)方差函數(shù),接下來(lái)還需進(jìn)一步研究協(xié)方差函數(shù)的確定方法,從而進(jìn)一步提高最小二乘配置推估的精度。

    致謝:感謝東華理工大學(xué)王勝平博士提供多波束測(cè)深數(shù)據(jù)。

    [1] 朱慶, 李德仁. 多波束測(cè)深數(shù)據(jù)的誤差分析與處理[J]. 武漢測(cè)繪科技大學(xué)學(xué)報(bào), 1998, 23(1): 1-4, 46. ZHU Qing, LI Deren. Error Analysis and Processing of Multibeam Soundings[J]. Journal of Wuhan Technical University of Surveying and Mapping, 1998, 23(1): 1-4, 46.

    [2] 陸秀平, 黃謨濤, 翟國(guó)君, 等. 多波束測(cè)深數(shù)據(jù)處理關(guān)鍵技術(shù)研究進(jìn)展與展望[J]. 海洋測(cè)繪, 2016, 36(4): 1-6, 11.LU Xiuping, HUANG Motao, ZHAI Guojun, et al. Development and Prospect of Key Technologies for Multibeam Echosounding Data Processing[J]. Hydrographic Surveying and Charting, 2016, 36(4): 1-6, 11.

    [3] 王海棟, 柴洪洲, 王敏. 多波束測(cè)深數(shù)據(jù)的抗差Kriging擬合[J]. 測(cè)繪學(xué)報(bào), 2011, 40(2): 238-242, 248. WANG Haidong, CHAI Hongzhou, WANG Min. Multibeam Bathymetry Fitting Based on Robust Kriging[J]. Acta Geodaetica et Cartographica Sinica, 2011, 40(2): 238-242, 248.

    [4] BORRE K. Mathematical Foundation of Geodesy: Selected Papers of Torben Krarup[M]. Berlin Heidelberg: Springer, 2006.

    [5] MORITZ H. Advanced Physical Geodesy[M]. Karlsruhe: Wichmann, Abacus Press, 1980.

    [6] WOLF H.über Verallgemeinerte Kollokation[J]. Zeitschrift für Vermessungswesen, 1974, 99(11): 475-478.

    [7] EL-FIKY G S, KATO T. Continuous Distribution of the Horizontal Strain in the Tohoku District, Japan, Predicted by Least-squares Collocation[J]. Journal of Geodynamics, 1998, 27(2): 213-236.

    [8] 柴洪洲, 崔岳, 明鋒. 最小二乘配置方法確定中國(guó)大陸主要塊體運(yùn)動(dòng)模型[J]. 測(cè)繪學(xué)報(bào), 2009, 38(1): 61-65. DOI: 10.3321/j.issn:1001-1595.2009.01.011. CHAI Hongzhou,CUI Yue,MING Feng.The Determination of Chinese Mainland Crustal Movement Model Using Least-squares Collocation[J]. Acta Geodaetica et Cartographica Sinica, 2009, 38(1): 61-65. DOI: 10.3321/j.issn:1001-1595.2009.01.011.

    [9] 江在森, 劉經(jīng)南. 應(yīng)用最小二乘配置建立地殼運(yùn)動(dòng)速度場(chǎng)與應(yīng)變場(chǎng)的方法[J]. 地球物理學(xué)報(bào), 2010, 53(5): 1109-1117. JIANG Zaisen, LIU Jingnan. The Method in Establishing Strain Field and Velocity Field of Crustal Movement Using Least Squares Collocation[J]. Chinese Journal of Geophysics, 2010, 53(5): 1109-1117.

    [10] YANG Y X,ZENG A M,ZHANG J.Adaptive Collocation with Application in Height System Transformation[J]. Journal of Geodesy, 2009, 83(5): 403-410.

    [11] JARMOLOWSKI W. Least Squares Collocation with Uncorrelated Heterogeneous Noise Estimated by Restricted Maximum Likelihood[J]. Journal of Geodesy, 2015, 89(6): 577-589.

    [12] YOU R J, HWANG H W. Coordinate Transformation between Two Geodetic Datums of Taiwan by Least-Squares Collocation[J]. Journal of Surveying Engineering, 2006, 132(2): 64-70.

    [13] 徐衛(wèi)明, 劉雁春. 海道測(cè)量中非定位點(diǎn)的精密確定方法[J]. 系統(tǒng)仿真學(xué)報(bào), 2007, 19(20): 4612-4615. XU Weiming, LIU Yanchun. Methods of Determining Non-Position Points in Hydrographic Survey[J]. Journal of System Simulation, 2007, 19(20): 4612-4615.

    [14] 王海棟, 柴洪洲, 趙東明. 基于抗差最小二乘配置的海底地形生成研究[J]. 系統(tǒng)仿真學(xué)報(bào), 2010, 22(9): 2091-2094, 2099. WANG Haidong, CHAI Hongzhou, ZHAO Dongming. Research on Seabed Terrain Generation Based on Robust Least-squares Collocation[J]. Journal of System Simulation, 2010, 22(9): 2091-2094, 2099.

    [15] 王海棟. 多波束系統(tǒng)測(cè)深異常處理理論與方法研究[D]. 鄭州: 信息工程大學(xué), 2010. WANG Haidong.Research on the Theories and Algorithms of Processing the Bathymetry Outlier in Multibeam Echosounder System[D]. Zhengzhou: Information Engineering University, 2010.

    [16] 薛樹(shù)強(qiáng), 楊元喜. 廣義反距離加權(quán)空間推估法[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2013, 38(12): 1435-1439. XUE Shuqiang, YANG Yuanxi. Generalized Inverse Distance Weighting Method for Spatial Interpolation[J]. Geomatics and Information Science of Wuhan University, 2013, 38(12): 1435-1439.

    [17] 黃維彬. 近代平差理論及其應(yīng)用[M]. 北京: 解放軍出版社, 1992: 92-118. HUANG Weibin. Theory of Adjustment and Its Applications in Geodesy[M]. Beijing: The PLA Press, 1992: 92-118.

    [18] YANG Yuanxi. Robustfying Collocation[J]. Manuscripta Geodaetica, 1992, 17(1): 21-28.

    [19] HARDY R L. Multiquadric Equations of Topography and Other Irregular Surfaces[J]. Journal of Geophysical Research, 1971, 76(8): 1905-1915.

    [20] HARDY R L. The Application of Multiquadric Equations and Point Mass Anomaly Models to Crustal Movement Studies[R]. NOAA Technical Report NOS 76 NGS 11. Rockville, Maryland: US Department of Commerce, National Oceanic and Atmospheric Administration, National Ocean Survey, National Geodetic Survey, 1978.

    [21] 謝曦霖, 許才軍, 溫?fù)P茂, 等. 一種基于多面函數(shù)的改進(jìn)最小二乘配置方法[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版). DOI: 10.13203/j.whugis20150664. XIE Xilin, XU Caijun, WEN Yangmao, et al. A Refined Least Squares Collocation Method Based on Multiquadric Function[J]. Geomatics and Information Science of Wuhan University. DOI: 10.13203/j.whugis20150664.

    [22] 楊元喜. 抗差估計(jì)理論及其應(yīng)用[M]. 北京: 八一出版社, 1993. YANG Yuanxi. Robust Estimation Theory and Its Applications[M]. Beijing: Bayi Publishing House, 1993.

    [23] SCHAFFRIN B, BOCK Y. Geodetic Deformation Analysis Based on Robust Inverse Theory[J]. Manuscripta Geodaetica, 1994, 19(1): 31-44.

    [24] 劉念. 協(xié)方差函數(shù)的抗差擬合[J]. 測(cè)繪科學(xué), 2001, 26(3): 25-28. LIU Nian. Robust Fitting of Covariance Function[J]. Science of Surveying and Mapping, 2001, 26(3): 25-28.

    [25] BROUNS G, DE WULF A, CONSTALES D. Multibeam Data Processing: Adding and Deleting Vertices in a Delaunay Triangulation[J]. Hydrographic Journal, 2001(101): 3-9.

    (責(zé)任編輯:陳品馨)

    WANG Leyang (1983—), male, PhD, associate professor, majors in geodetic inversion and geodetic data processing.

    Multi-beam Bathymetry Data Processing Using Iterative Algorithm of Robust Least Squares Collocation

    WANG Leyang1,2,3,CHEN Hanqing1,2

    1. Faculty of Geomatics, East China University of Technology, Nanchang 330013, China; 2. Key Laboratory of Watershed Ecology and Geographical Environment Monitoring, NASG, Nanchang 330013, China; 3. Key Laboratory for Digital Land and Resources of Jiangxi Province, Nanchang 330013, China

    In the process of dealing with multi-beam bathymetry data by least squares collocation, the quadric curved mathematical model of trend term can not express accurately the whole variation trend of seafloor topography in general. Moreover, the covariance function estimated by general method is incapable of accurately expressing statistical characteristics with the multi-beam bathymetry data contains gross errors or outliers. So the iterative algorithm of robust least squares collocation is proposed in this paper. Firstly, the initial weight matrix of observations and the initial parameters of covariance function are both given in this method, then the trend term is fitted by polyhedral function and equivalent weights scheme is applied into robust estimation in this method. Finally, the robust parameters of covariance function and solutions of least squares collocation are iteratively calculated. The experimental results show that the method proposed in this paper can express well the whole variation trend of seafloor topography and overcome the effect of gross error or outlier in multi-beam bathymetry data to a certain extent. Compared with the conventional robust method, the proposed method in this paper more effectively probes the outliers in bathymetry data with the robust and better predicted results.

    least squares collocation; seabed terrain generation; covariance function; robust; trend term

    National Natural Science Foundation of China (Nos.41664001; 41204003); Support Program for Outstanding Youth Talents in Jiangxi Province (No.20162BCB23050); National Key Research and Development Program(No.2016YFB0501405); Science and Technology Project of the Education Department of Jiangxi Province (No.GJJ150595); the Project of Key Laboratory for Digital Land and Resources of Jiangxi Province(No.DLLJ201705).

    王樂(lè)洋,陳漢清.多波束測(cè)深數(shù)據(jù)處理的抗差最小二乘配置迭代解法[J].測(cè)繪學(xué)報(bào),2017,46(5):658-665.

    10.11947/j.AGCS.2017.20160491. WANG Leyang, CHEN Hanqing.Multi-beam Bathymetry Data Processing Using Iterative Algorithm of Robust Least Squares Collocation[J]. Acta Geodaetica et Cartographica Sinica,2017,46(5):658-665. DOI:10.11947/j.AGCS.2017.20160491.

    2016-10-17

    王樂(lè)洋(1983—),男,博士,副教授,主要研究方向?yàn)榇蟮販y(cè)量反演及大地測(cè)量數(shù)據(jù)處理。

    E-mail: wleyang@163.com

    P229

    A

    1001-1595(2017)05-0658-08

    國(guó)家自然科學(xué)基金 (41664001; 41204003);江西省杰出青年人才資助計(jì)劃項(xiàng)目(20162BCB23050);國(guó)家重點(diǎn)研發(fā)計(jì)劃(2016YFB0501405);江西省教育廳科技項(xiàng)目(GJJ150595) ;江西省數(shù)字國(guó)土重點(diǎn)實(shí)驗(yàn)室開(kāi)放研究基金(DLLJ201705)

    修回日期: 2017-04-13

    猜你喜歡
    抗差協(xié)方差波束
    毫米波大規(guī)模陣列天線波束掃描研究*
    圓陣多波束測(cè)角探究
    Helix陣匹配場(chǎng)三維波束形成
    改善單頻PPP參數(shù)收斂速度的抗差估計(jì)方法
    不確定系統(tǒng)改進(jìn)的魯棒協(xié)方差交叉融合穩(wěn)態(tài)Kalman預(yù)報(bào)器
    一種基于廣義協(xié)方差矩陣的欠定盲辨識(shí)方法
    地形簡(jiǎn)化對(duì)DEM不確定性的抗差性研究
    基于抗差最小均方估計(jì)的輸電線路參數(shù)辨識(shí)
    基于非正交變換的局域波束空時(shí)自適應(yīng)處理
    縱向數(shù)據(jù)分析中使用滑動(dòng)平均Cholesky分解對(duì)回歸均值和協(xié)方差矩陣進(jìn)行同時(shí)半?yún)?shù)建模
    欧美激情高清一区二区三区| 中文字幕最新亚洲高清| 搡老熟女国产l中国老女人| 亚洲色图av天堂| 韩国精品一区二区三区| 免费在线观看完整版高清| 国产精品久久久久久人妻精品电影| bbb黄色大片| 黄色a级毛片大全视频| 成人国产综合亚洲| 在线视频色国产色| 女性生殖器流出的白浆| 真人一进一出gif抽搐免费| 国产av一区在线观看免费| 亚洲一区二区三区不卡视频| 亚洲自拍偷在线| 欧美性猛交黑人性爽| 真人一进一出gif抽搐免费| svipshipincom国产片| 99国产精品一区二区三区| 午夜福利在线在线| av有码第一页| 亚洲欧美激情综合另类| 真人做人爱边吃奶动态| 亚洲国产欧美日韩在线播放| 国产精品日韩av在线免费观看| 后天国语完整版免费观看| 中文在线观看免费www的网站 | 成人特级黄色片久久久久久久| 亚洲真实伦在线观看| 美女免费视频网站| 午夜免费成人在线视频| 午夜影院日韩av| 国产亚洲欧美在线一区二区| 国内精品久久久久久久电影| 韩国av一区二区三区四区| 999久久久国产精品视频| 亚洲精品色激情综合| 久久精品成人免费网站| 国产久久久一区二区三区| 色婷婷久久久亚洲欧美| 天天躁狠狠躁夜夜躁狠狠躁| 69av精品久久久久久| 老司机午夜福利在线观看视频| 亚洲国产日韩欧美精品在线观看 | 成年女人毛片免费观看观看9| 在线观看舔阴道视频| 国产精品综合久久久久久久免费| 国产亚洲欧美在线一区二区| 性欧美人与动物交配| 欧美乱码精品一区二区三区| 国产精品久久久久久亚洲av鲁大| 久久久久久久午夜电影| 国产极品粉嫩免费观看在线| 国产熟女xx| 久久久久亚洲av毛片大全| 久久久久久久午夜电影| 久久这里只有精品19| 一级毛片高清免费大全| 91字幕亚洲| 国产三级黄色录像| 国产极品粉嫩免费观看在线| 老司机深夜福利视频在线观看| 2021天堂中文幕一二区在线观 | av免费在线观看网站| 午夜a级毛片| 成人午夜高清在线视频 | 亚洲久久久国产精品| 亚洲五月婷婷丁香| 黄片播放在线免费| 亚洲精品国产区一区二| 天堂动漫精品| 老司机午夜十八禁免费视频| 女性被躁到高潮视频| 97人妻精品一区二区三区麻豆 | 亚洲av中文字字幕乱码综合 | 一卡2卡三卡四卡精品乱码亚洲| 欧美绝顶高潮抽搐喷水| 美女免费视频网站| 99久久无色码亚洲精品果冻| 啦啦啦观看免费观看视频高清| 老司机靠b影院| 久久精品影院6| 麻豆成人午夜福利视频| 日日摸夜夜添夜夜添小说| 日韩免费av在线播放| 淫秽高清视频在线观看| 97超级碰碰碰精品色视频在线观看| 欧美日本视频| 国产免费av片在线观看野外av| 中国美女看黄片| 日本一本二区三区精品| 999精品在线视频| 天堂√8在线中文| 99精品久久久久人妻精品| 丁香六月欧美| 97人妻精品一区二区三区麻豆 | 韩国av一区二区三区四区| 亚洲国产毛片av蜜桃av| 欧美在线一区亚洲| 正在播放国产对白刺激| 亚洲 欧美 日韩 在线 免费| a级毛片a级免费在线| 在线av久久热| 精品不卡国产一区二区三区| 国产一区二区激情短视频| 久久香蕉激情| 中文字幕精品免费在线观看视频| 黄片小视频在线播放| 亚洲成av人片免费观看| 男男h啪啪无遮挡| 午夜免费成人在线视频| 大型av网站在线播放| 91麻豆av在线| 色av中文字幕| 老汉色∧v一级毛片| 亚洲精品国产精品久久久不卡| 美女 人体艺术 gogo| 精品卡一卡二卡四卡免费| 国产野战对白在线观看| 欧美zozozo另类| 亚洲第一电影网av| 亚洲国产高清在线一区二区三 | 国产成人啪精品午夜网站| 午夜激情av网站| 国产91精品成人一区二区三区| 久久性视频一级片| 可以在线观看毛片的网站| 嫩草影视91久久| 欧美日韩亚洲综合一区二区三区_| 久久精品国产99精品国产亚洲性色| 久久久久国产精品人妻aⅴ院| 久久精品aⅴ一区二区三区四区| www.www免费av| 99精品在免费线老司机午夜| 亚洲国产看品久久| 视频在线观看一区二区三区| 麻豆一二三区av精品| 亚洲精品一区av在线观看| 少妇裸体淫交视频免费看高清 | 久久国产精品人妻蜜桃| 午夜视频精品福利| 99国产综合亚洲精品| 亚洲欧美精品综合久久99| 婷婷精品国产亚洲av| 在线免费观看的www视频| 日本一区二区免费在线视频| 免费在线观看视频国产中文字幕亚洲| 一区福利在线观看| 欧美色视频一区免费| 色综合婷婷激情| 97人妻精品一区二区三区麻豆 | 欧美激情久久久久久爽电影| 一二三四社区在线视频社区8| 三级毛片av免费| 老司机深夜福利视频在线观看| 他把我摸到了高潮在线观看| x7x7x7水蜜桃| 亚洲五月色婷婷综合| 国产一区二区在线av高清观看| 在线av久久热| 老司机深夜福利视频在线观看| 18禁观看日本| 午夜福利欧美成人| 天堂动漫精品| 黄频高清免费视频| 亚洲人成电影免费在线| 国产又黄又爽又无遮挡在线| 女性被躁到高潮视频| 国产亚洲精品第一综合不卡| 久久香蕉国产精品| 两性夫妻黄色片| 在线播放国产精品三级| 757午夜福利合集在线观看| 亚洲国产毛片av蜜桃av| 久久久久久久午夜电影| 无人区码免费观看不卡| 亚洲 欧美一区二区三区| 国产一级毛片七仙女欲春2 | 99热这里只有精品一区 | 99久久国产精品久久久| 亚洲av电影不卡..在线观看| 成人精品一区二区免费| av中文乱码字幕在线| 久久久精品欧美日韩精品| 人人妻人人澡人人看| 黄色丝袜av网址大全| 日韩精品中文字幕看吧| 精品人妻1区二区| 国产黄a三级三级三级人| 黄片大片在线免费观看| 法律面前人人平等表现在哪些方面| 美国免费a级毛片| 久久精品人妻少妇| 十八禁网站免费在线| 一本久久中文字幕| 天堂影院成人在线观看| 老熟妇乱子伦视频在线观看| 中文字幕人成人乱码亚洲影| 搡老妇女老女人老熟妇| av片东京热男人的天堂| 麻豆一二三区av精品| 国产精品综合久久久久久久免费| 国产黄a三级三级三级人| 国产真实乱freesex| 波多野结衣巨乳人妻| 99re在线观看精品视频| 99在线视频只有这里精品首页| 搞女人的毛片| 国产一区二区激情短视频| 亚洲精品国产一区二区精华液| 久久久久久免费高清国产稀缺| 中文字幕另类日韩欧美亚洲嫩草| 99精品在免费线老司机午夜| 亚洲第一青青草原| 亚洲第一欧美日韩一区二区三区| 精品电影一区二区在线| 一个人观看的视频www高清免费观看 | 国产蜜桃级精品一区二区三区| 黄色片一级片一级黄色片| 成人国产综合亚洲| 视频区欧美日本亚洲| 成人18禁在线播放| 波多野结衣高清无吗| 久久精品成人免费网站| 成人永久免费在线观看视频| 国产午夜精品久久久久久| 精品午夜福利视频在线观看一区| 好男人电影高清在线观看| 免费人成视频x8x8入口观看| 99国产精品99久久久久| 国产av一区二区精品久久| 人人妻人人澡人人看| 又紧又爽又黄一区二区| 久久久精品欧美日韩精品| 成人18禁高潮啪啪吃奶动态图| 午夜免费成人在线视频| 久久九九热精品免费| 熟妇人妻久久中文字幕3abv| 久久久久久久午夜电影| 最新美女视频免费是黄的| 最好的美女福利视频网| 给我免费播放毛片高清在线观看| 白带黄色成豆腐渣| videosex国产| 高潮久久久久久久久久久不卡| 亚洲国产精品999在线| 制服诱惑二区| 变态另类成人亚洲欧美熟女| 欧美日韩亚洲综合一区二区三区_| 国产高清激情床上av| 色播亚洲综合网| 亚洲 欧美一区二区三区| 婷婷六月久久综合丁香| e午夜精品久久久久久久| 亚洲精品中文字幕在线视频| 桃红色精品国产亚洲av| 少妇熟女aⅴ在线视频| 天堂√8在线中文| 成人18禁在线播放| 国产成人系列免费观看| 一级作爱视频免费观看| 久热爱精品视频在线9| 久久久国产精品麻豆| 精品第一国产精品| 国产主播在线观看一区二区| 午夜福利成人在线免费观看| 国产成人系列免费观看| 国产成人一区二区三区免费视频网站| 色哟哟哟哟哟哟| 精品久久久久久久末码| 欧美国产精品va在线观看不卡| 日韩三级视频一区二区三区| 成人18禁高潮啪啪吃奶动态图| 国产一区在线观看成人免费| 可以免费在线观看a视频的电影网站| 国产精品综合久久久久久久免费| 亚洲色图 男人天堂 中文字幕| 老司机在亚洲福利影院| 欧美日韩乱码在线| 在线国产一区二区在线| 国产精品二区激情视频| 亚洲国产精品久久男人天堂| 国产精品香港三级国产av潘金莲| 精品午夜福利视频在线观看一区| 国产精品久久久av美女十八| 99久久国产精品久久久| 久久婷婷成人综合色麻豆| 18美女黄网站色大片免费观看| 女人爽到高潮嗷嗷叫在线视频| 十八禁网站免费在线| 亚洲第一电影网av| 男人舔女人下体高潮全视频| 成年女人毛片免费观看观看9| 亚洲av成人不卡在线观看播放网| 无人区码免费观看不卡| 99久久久亚洲精品蜜臀av| 中国美女看黄片| 成人国产综合亚洲| 日本a在线网址| 精品国产国语对白av| 亚洲久久久国产精品| 黄色丝袜av网址大全| 中文字幕人成人乱码亚洲影| 成人三级做爰电影| 午夜影院日韩av| 午夜老司机福利片| 99久久久亚洲精品蜜臀av| 久久精品成人免费网站| 亚洲中文日韩欧美视频| 久久人人精品亚洲av| av中文乱码字幕在线| 2021天堂中文幕一二区在线观 | 色综合欧美亚洲国产小说| 欧美日韩精品网址| 国产成年人精品一区二区| 免费在线观看亚洲国产| 亚洲av五月六月丁香网| 亚洲男人天堂网一区| 国内精品久久久久久久电影| netflix在线观看网站| 男女做爰动态图高潮gif福利片| 欧美成狂野欧美在线观看| 人人妻人人看人人澡| 午夜福利欧美成人| 麻豆国产av国片精品| 欧美绝顶高潮抽搐喷水| 精品乱码久久久久久99久播| 亚洲免费av在线视频| 国语自产精品视频在线第100页| 哪里可以看免费的av片| 色av中文字幕| 国产午夜福利久久久久久| 国产高清有码在线观看视频 | 777久久人妻少妇嫩草av网站| 日韩欧美免费精品| 真人做人爱边吃奶动态| 色综合站精品国产| 国产欧美日韩一区二区精品| 少妇粗大呻吟视频| 欧美日韩福利视频一区二区| 欧美+亚洲+日韩+国产| 一级作爱视频免费观看| 99国产精品一区二区三区| 免费看美女性在线毛片视频| 亚洲自偷自拍图片 自拍| 欧美性猛交黑人性爽| 国产精品野战在线观看| 视频区欧美日本亚洲| 中文字幕久久专区| 亚洲久久久国产精品| 亚洲天堂国产精品一区在线| 岛国视频午夜一区免费看| 熟女少妇亚洲综合色aaa.| 午夜福利欧美成人| 国产激情久久老熟女| 国产在线观看jvid| 国产区一区二久久| 精品少妇一区二区三区视频日本电影| 久久久精品欧美日韩精品| 国产在线精品亚洲第一网站| 国产精品,欧美在线| 国产精品永久免费网站| 91麻豆精品激情在线观看国产| 成人三级做爰电影| 久久亚洲精品不卡| 亚洲欧洲精品一区二区精品久久久| 久久亚洲精品不卡| 正在播放国产对白刺激| 在线观看一区二区三区| av片东京热男人的天堂| 在线视频色国产色| 午夜福利成人在线免费观看| 一区福利在线观看| 国产高清videossex| 香蕉av资源在线| 成人18禁在线播放| 亚洲第一青青草原| 免费观看精品视频网站| 女性生殖器流出的白浆| 韩国av一区二区三区四区| 丝袜美腿诱惑在线| 一本大道久久a久久精品| 欧美日韩亚洲综合一区二区三区_| av电影中文网址| 叶爱在线成人免费视频播放| 老司机深夜福利视频在线观看| 一夜夜www| 欧美精品啪啪一区二区三区| 久久精品91无色码中文字幕| 亚洲五月色婷婷综合| 午夜免费激情av| 亚洲三区欧美一区| 午夜精品在线福利| 天天躁夜夜躁狠狠躁躁| 三级毛片av免费| 一区二区三区高清视频在线| 三级毛片av免费| 午夜福利在线在线| 精品久久蜜臀av无| 在线观看一区二区三区| 中文在线观看免费www的网站 | a在线观看视频网站| 动漫黄色视频在线观看| 给我免费播放毛片高清在线观看| 日本免费一区二区三区高清不卡| 国内精品久久久久久久电影| 97人妻精品一区二区三区麻豆 | 久久久久亚洲av毛片大全| 两人在一起打扑克的视频| 国产亚洲精品第一综合不卡| 99热这里只有精品一区 | 中文亚洲av片在线观看爽| 高清毛片免费观看视频网站| 色哟哟哟哟哟哟| 女警被强在线播放| 人人妻人人澡欧美一区二区| 久久精品国产99精品国产亚洲性色| 日韩 欧美 亚洲 中文字幕| 亚洲中文日韩欧美视频| 久久精品aⅴ一区二区三区四区| 久久香蕉激情| 亚洲五月色婷婷综合| 两个人视频免费观看高清| 国产激情偷乱视频一区二区| 久久精品国产清高在天天线| 首页视频小说图片口味搜索| 免费在线观看影片大全网站| 一个人观看的视频www高清免费观看 | 又紧又爽又黄一区二区| 大型av网站在线播放| 9191精品国产免费久久| 国产欧美日韩精品亚洲av| 欧洲精品卡2卡3卡4卡5卡区| 亚洲真实伦在线观看| 精品熟女少妇八av免费久了| 久久人人精品亚洲av| 国产亚洲欧美精品永久| 丰满的人妻完整版| 高潮久久久久久久久久久不卡| 国产精品日韩av在线免费观看| 国产熟女xx| 51午夜福利影视在线观看| 曰老女人黄片| 十八禁人妻一区二区| 国产高清视频在线播放一区| 国产日本99.免费观看| 精品久久久久久,| www.精华液| 18禁国产床啪视频网站| 国产激情欧美一区二区| av在线天堂中文字幕| 曰老女人黄片| 男人操女人黄网站| 国产精品影院久久| 国产精品久久电影中文字幕| 精品乱码久久久久久99久播| 午夜日韩欧美国产| 最近最新中文字幕大全免费视频| 欧美性长视频在线观看| 久久精品成人免费网站| 18禁观看日本| 成人三级做爰电影| 日本精品一区二区三区蜜桃| 国产激情偷乱视频一区二区| 国产高清有码在线观看视频 | 99久久久亚洲精品蜜臀av| 欧美色欧美亚洲另类二区| 亚洲成人精品中文字幕电影| 热99re8久久精品国产| 老熟妇仑乱视频hdxx| 久久久久九九精品影院| 少妇熟女aⅴ在线视频| 国产一区二区三区在线臀色熟女| 99精品在免费线老司机午夜| 免费电影在线观看免费观看| 熟女少妇亚洲综合色aaa.| 美女大奶头视频| 女人被狂操c到高潮| 亚洲真实伦在线观看| 99热6这里只有精品| 少妇粗大呻吟视频| a在线观看视频网站| 在线观看免费视频日本深夜| 日本 欧美在线| 婷婷精品国产亚洲av| 免费搜索国产男女视频| 曰老女人黄片| 亚洲av熟女| 日本熟妇午夜| 日日干狠狠操夜夜爽| 美国免费a级毛片| 欧美日韩亚洲综合一区二区三区_| 国产熟女xx| 日韩欧美免费精品| 中文字幕精品免费在线观看视频| 精品国产乱码久久久久久男人| 香蕉久久夜色| 一个人免费在线观看的高清视频| 一级毛片高清免费大全| av欧美777| 亚洲国产中文字幕在线视频| 男女下面进入的视频免费午夜 | 熟妇人妻久久中文字幕3abv| 制服丝袜大香蕉在线| 欧美精品啪啪一区二区三区| 日韩av在线大香蕉| 久久久久国产精品人妻aⅴ院| 51午夜福利影视在线观看| 在线十欧美十亚洲十日本专区| 黄频高清免费视频| 国产精品国产高清国产av| 亚洲aⅴ乱码一区二区在线播放 | 亚洲成人免费电影在线观看| 午夜精品在线福利| 大型黄色视频在线免费观看| 别揉我奶头~嗯~啊~动态视频| 国产99白浆流出| 可以在线观看的亚洲视频| 日韩精品青青久久久久久| 一进一出抽搐动态| 12—13女人毛片做爰片一| 日本一区二区免费在线视频| 91字幕亚洲| 亚洲成a人片在线一区二区| 免费看a级黄色片| 亚洲免费av在线视频| www日本黄色视频网| √禁漫天堂资源中文www| 久久亚洲精品不卡| xxx96com| 国产三级黄色录像| 美国免费a级毛片| 日韩av在线大香蕉| 黄色成人免费大全| 这个男人来自地球电影免费观看| 岛国在线观看网站| 夜夜爽天天搞| 免费无遮挡裸体视频| 男女之事视频高清在线观看| 国产精品亚洲美女久久久| 久久精品国产亚洲av高清一级| 色老头精品视频在线观看| 免费搜索国产男女视频| 又紧又爽又黄一区二区| 成人手机av| 欧美三级亚洲精品| 18禁裸乳无遮挡免费网站照片 | 美女午夜性视频免费| 国产v大片淫在线免费观看| 亚洲片人在线观看| 午夜成年电影在线免费观看| 亚洲激情在线av| 国内精品久久久久精免费| 国产午夜精品久久久久久| 精品国产国语对白av| 亚洲av第一区精品v没综合| 国产黄片美女视频| 男人的好看免费观看在线视频 | 日韩精品青青久久久久久| 中文字幕精品亚洲无线码一区 | av视频在线观看入口| 国产精品99久久99久久久不卡| 桃红色精品国产亚洲av| 成人av一区二区三区在线看| 亚洲欧美激情综合另类| 巨乳人妻的诱惑在线观看| 99热这里只有精品一区 | 亚洲中文日韩欧美视频| cao死你这个sao货| 欧美成人午夜精品| 欧美日韩福利视频一区二区| av中文乱码字幕在线| 久久精品国产99精品国产亚洲性色| 午夜福利欧美成人| www.精华液| 国产免费av片在线观看野外av| 国产又黄又爽又无遮挡在线| 天堂动漫精品| 一级毛片高清免费大全| 亚洲无线在线观看| 久久精品国产综合久久久| av在线天堂中文字幕| 国产免费男女视频| 人成视频在线观看免费观看| 久久久久国产精品人妻aⅴ院| 色精品久久人妻99蜜桃| 亚洲国产高清在线一区二区三 | 国产久久久一区二区三区| 久久国产精品男人的天堂亚洲| 成人三级做爰电影| 久久精品国产亚洲av香蕉五月| 国产成人av激情在线播放| 国产精品久久视频播放| 国产亚洲精品久久久久久毛片| 日日干狠狠操夜夜爽| 熟妇人妻久久中文字幕3abv| 精品久久久久久久末码| 国产成人欧美| 在线av久久热| 他把我摸到了高潮在线观看| 窝窝影院91人妻| 黄频高清免费视频| 一边摸一边做爽爽视频免费| 99精品久久久久人妻精品| 日韩成人在线观看一区二区三区| 黄色丝袜av网址大全| 国产一区二区激情短视频| 国产精品乱码一区二三区的特点| 国产成人影院久久av| 老熟妇仑乱视频hdxx| 可以免费在线观看a视频的电影网站| 亚洲激情在线av| 久久久久久久久中文| 免费人成视频x8x8入口观看|