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

    ICESat激光高程點(diǎn)輔助的天繪一號(hào)衛(wèi)星影像立體區(qū)域網(wǎng)平差

    2018-03-31 01:43:21張祖勛陶鵬杰宋孟肖
    測(cè)繪學(xué)報(bào) 2018年3期
    關(guān)鍵詞:網(wǎng)平控制點(diǎn)立體

    王 晉,張 勇,2,張祖勛,李 曉,陶鵬杰,4,宋孟肖

    1. 武漢大學(xué)遙感信息工程學(xué)院,湖北 武漢 430079; 2. 地球空間信息技術(shù)協(xié)同創(chuàng)新中心,湖北 武漢 430079; 3. 中國(guó)天繪衛(wèi)星中心,北京 102102; 4. 武漢大學(xué)資源與環(huán)境科學(xué)學(xué)院,湖北 武漢 430079

    衛(wèi)星影像的精度不僅取決于影像測(cè)量精度,而且取決于事后處理控制點(diǎn)的幾何分布和精度。然而,全球測(cè)圖難免遇到無(wú)控制點(diǎn)或少控制點(diǎn)區(qū)域,如我國(guó)青藏高原、云貴高山峽谷、西北沙漠等外業(yè)測(cè)控困難區(qū)域以及境外大量無(wú)可用控制點(diǎn)區(qū)域。所以,近年來(lái),立體衛(wèi)星影像無(wú)控區(qū)域網(wǎng)平差逐漸引起了廣泛的關(guān)注。衛(wèi)星影像無(wú)控區(qū)域網(wǎng)平差是指在缺少外業(yè)控制的條件下,實(shí)現(xiàn)衛(wèi)星影像定位精度的相對(duì)一致性且盡可能地提高對(duì)地定位精度。

    無(wú)控區(qū)域網(wǎng)平差的對(duì)地定位精度受到衛(wèi)星影像質(zhì)量、初始定位精度、地形起伏等因素的影響。為了提高對(duì)地定位精度,國(guó)內(nèi)外研究人員開展了一系列的相關(guān)研究。文獻(xiàn)[1]研究了采用多次覆蓋的衛(wèi)星影像進(jìn)行無(wú)控區(qū)域網(wǎng)平差方法,可改善衛(wèi)星影像無(wú)控定位精度;文獻(xiàn)[2]研究并實(shí)現(xiàn)了覆蓋全國(guó)范圍的8802個(gè)資源三號(hào)三線陣立體像對(duì)在無(wú)控制條件下整體一張網(wǎng)的平差,達(dá)到無(wú)控制測(cè)圖精度優(yōu)于5 m的目標(biāo);文獻(xiàn)[3]利用更高定位精度的高分辨率衛(wèi)星影像(如GeoEye、WorldView等)或航空影像作為控制條件,也取得了較好的效果。除此之外,已有的地理信息數(shù)據(jù),如數(shù)字正射影像(digital orthophoto map,DOM)和數(shù)字高程模型(digital elevation model,DEM),也可以作為控制數(shù)據(jù),輔助衛(wèi)星影像區(qū)域網(wǎng)平差[4]。隨著公眾地理信息數(shù)據(jù)種類的多樣化和精度的提升,平差過(guò)程中引入公開且可穩(wěn)定獲取的公眾地理信息數(shù)據(jù)(如Google Earth影像、SRTM DEM)也逐漸成為衛(wèi)星影像無(wú)控平差的可選擇方式。其中ICESat衛(wèi)星獲取的激光高程點(diǎn)數(shù)據(jù)就是一種重要的公眾地理信息數(shù)據(jù)。ICESat衛(wèi)星是NASA于2003年發(fā)射的一顆地球觀測(cè)衛(wèi)星,星上攜帶有NASA設(shè)計(jì)的地面高程激光測(cè)量系統(tǒng)(geoscience laser altimeter system,GLAS),ICESat衛(wèi)星兼顧了地形及地面植被高度測(cè)量、云層高度及其垂直結(jié)構(gòu)測(cè)量、海洋高程測(cè)量等科研任務(wù)。雖然ICESat衛(wèi)星已于2009年10月停止了工作,但是所獲取的大量的激光觀測(cè)數(shù)據(jù)仍然具有重要的應(yīng)用價(jià)值,這些觀測(cè)數(shù)據(jù)已經(jīng)被廣泛應(yīng)用于冰川高度變化監(jiān)測(cè)[5-7]、湖泊河流水位高度變化監(jiān)測(cè)、植被高度統(tǒng)計(jì)[8-9]、數(shù)字高程模型的精度評(píng)價(jià)等方面[10-11]。

    雖然ICESat激光高程點(diǎn)具有較高的測(cè)量精度,但是利用其輔助衛(wèi)星影像區(qū)域網(wǎng)平差方面的研究才剛剛起步。文獻(xiàn)[12]開展了ICESat激光高程控制點(diǎn)自動(dòng)提取方法的研究,取得了初步的成果,能有效保證提取的高程控制點(diǎn)的絕對(duì)高程精度。文獻(xiàn)[13—14]討論了二線陣影像聯(lián)合激光測(cè)距數(shù)據(jù)進(jìn)行光束法平差的可行性,并進(jìn)行了模擬數(shù)據(jù)試驗(yàn),結(jié)果表明,利用激光測(cè)距數(shù)據(jù)參與二線陣影像光束法平差,能有效改善航線模型系統(tǒng)變形并保持較小的上下視差。文獻(xiàn)[15]利用嫦娥一號(hào)衛(wèi)星三線陣光學(xué)影像和激光測(cè)高數(shù)據(jù),采用局部表面約束的聯(lián)合平差思想進(jìn)行集成處理,最終生成了月球表面三維數(shù)字表面模型,為探測(cè)器落月選址奠定基礎(chǔ)。文獻(xiàn)[16—17]采用經(jīng)過(guò)篩選后的ICESat激光測(cè)高數(shù)據(jù)作為廣義高程控制點(diǎn),將平面和高程控制分離,輔助資源三號(hào)進(jìn)行三線陣立體區(qū)域網(wǎng)平差,驗(yàn)證了ICESat激光點(diǎn)輔助衛(wèi)星影像區(qū)域網(wǎng)平差的正確性和可行性;其試驗(yàn)結(jié)論為未來(lái)我國(guó)發(fā)射激光測(cè)高衛(wèi)星,開展衛(wèi)星激光測(cè)高與光學(xué)立體影像聯(lián)合處理提供了很好的參考。

    綜上所述,在衛(wèi)星影像攝影測(cè)量處理中引入ICESat激光測(cè)高數(shù)據(jù)是十分有意義的。然而,激光高程點(diǎn)輔助衛(wèi)星影像區(qū)域網(wǎng)平差,目前僅停留在平差原理和模型的研究以及小范圍數(shù)據(jù)試驗(yàn)階段。本文設(shè)計(jì)了一種ICESat激光高程點(diǎn)自動(dòng)篩選提取方案,并提出了一種模型法衛(wèi)星影像立體區(qū)域網(wǎng)平差方法。為驗(yàn)證該方法的有效性和可行性,本文采用山東全省的天繪一號(hào)衛(wèi)星三線陣影像和ICESat激光測(cè)高數(shù)據(jù)進(jìn)行了區(qū)域網(wǎng)平差試驗(yàn),利用外業(yè)檢查點(diǎn)對(duì)其進(jìn)行精度評(píng)定,并討論了不同激光高程控制點(diǎn)布設(shè)方案對(duì)區(qū)域網(wǎng)平差結(jié)果的影響。

    1 ICESat/GLAS激光高程點(diǎn)的自動(dòng)提取

    雖然GLAS標(biāo)稱的測(cè)距精度為15 cm[9](測(cè)距方向?yàn)榻拼怪狈较?,在平坦地區(qū)地形高度變化較小,測(cè)距精度近似于高程精度),但在實(shí)際測(cè)量過(guò)程中,由于激光脈沖往返經(jīng)過(guò)大氣層,大氣分子及氣溶膠會(huì)對(duì)激光脈沖產(chǎn)生散射、能量衰減、光斑漂移、波形展寬等多重影響,地表的高反射率和云層的前向散射也會(huì)引起回波波形的過(guò)飽和,使得測(cè)距精度下降。因此,有必要從龐大的點(diǎn)數(shù)據(jù)庫(kù)中提取出適用于高程控制的GLAS激光點(diǎn)。文獻(xiàn)[12]研究發(fā)現(xiàn),利用GLAS L14數(shù)據(jù)記錄的激光點(diǎn)質(zhì)量評(píng)價(jià)參數(shù)(如飽和度改正標(biāo)記i_satCorrFlg,反射率參數(shù)i_reflctUncor等)和其他地表測(cè)量數(shù)據(jù)(如SRTM DEM,ASTER DEM)組合對(duì)GLAS點(diǎn)進(jìn)行篩選,可以達(dá)到較好的篩選效果。但是該研究采用90 m分辨率的SRTM DEM(SRTM 3″)作為參考標(biāo)準(zhǔn),而GLAS光斑本身只有70 m直徑,在最好的情況下,單個(gè)GLAS光斑也只能覆蓋4個(gè)參考DEM格網(wǎng),導(dǎo)致地面的坡度與粗糙度估計(jì)困難,給GLAS點(diǎn)篩選造成不確定性。從2014年起美國(guó)陸續(xù)向全球用戶免費(fèi)開放30 m分辨率的SRTM DEM數(shù)據(jù)SRTM 1″,與SRTM 3″一樣,SRTM 1″覆蓋了全球56°S—60°N范圍,標(biāo)稱的絕對(duì)高程精度為16 m(LE90),相對(duì)高程精度為10 m(LE90),但是具有更高的分辨率和數(shù)據(jù)質(zhì)量[18]。針對(duì)以上問(wèn)題,本文采用30 m分辨率的SRTM數(shù)據(jù)計(jì)算得到的地表坡度作為限制條件,同時(shí)根據(jù)GLAS本身的質(zhì)量評(píng)價(jià)信息進(jìn)行控制點(diǎn)篩選,具體流程如圖1所示。

    圖1 GLAS激光高程點(diǎn)自動(dòng)提取流程Fig.1 The workflow of automatic extraction of GLAS laser points

    在實(shí)際使用中發(fā)現(xiàn),很多GLAS激光高程點(diǎn)與SRTM的高程之差達(dá)1000多米,部分原因是GLAS發(fā)射的脈沖遇到厚云即反射,導(dǎo)致測(cè)量高度異常;此外,地面坡度較大的地區(qū),如果星下點(diǎn)誤差較大,也導(dǎo)致高程點(diǎn)精度異常,于是,必須首先進(jìn)行坡度估計(jì),并剔除兩類高程異常點(diǎn)。最后,通過(guò)GLAS點(diǎn)本身記錄的點(diǎn)質(zhì)量評(píng)價(jià)信息對(duì)剩余的點(diǎn)做進(jìn)一步篩選,即可獲取較為可靠的高程控制點(diǎn)。本文采用上述方法對(duì)山東測(cè)區(qū)的ICESat原始數(shù)據(jù)進(jìn)行了可用點(diǎn)提取,并以1∶10 000比例尺的DEM產(chǎn)品為參考,對(duì)提取的可用點(diǎn)進(jìn)行高程精度評(píng)定。參考DEM覆蓋山東全省約15.8萬(wàn)km2,地形以平地為主(約占55%),丘陵和山地為輔(各占13%和16%),其余為其他地形,評(píng)價(jià)的GLAS激光高程點(diǎn)共11 395個(gè)。高程點(diǎn)篩選情況如表1所示,利用SRTM高差、SRTM坡度、大氣散射增益和飽和度改正參數(shù)4個(gè)限制條件各剔除了3.6%、43.2%、26.2%和10.6%,利用最終保留的1836個(gè)高程點(diǎn)進(jìn)行精度評(píng)定,與參考DEM的高差的均值為0.296 m,中誤差為1.223 m,評(píng)定結(jié)果充分證明了可用點(diǎn)提取方法的有效性。

    表1基于SRTM的GLAS激光高程點(diǎn)篩選依據(jù)與結(jié)果

    Tab.1ThestatisticalresultofGLASlaserpointsbasedonSRTM

    約束項(xiàng)數(shù)據(jù)保留個(gè)數(shù)數(shù)據(jù)剔除率/(%)SRTM高差109853.6SRTM坡度603343.4大氣散射增益304526.2飽和度改正參數(shù)183610.6與DEM的高差/m均 值0.296中誤差1.223

    從表1也可以看出,利用SRTM坡度條件剔除的比例最高,達(dá)到了43.2%,因此,保留的高程點(diǎn)都位于地形平緩的平地區(qū)域,將其作為高程控制點(diǎn)更加合理。

    2 激光高程點(diǎn)輔助的衛(wèi)星影像立體區(qū)域網(wǎng)平差模型

    2.1 衛(wèi)星影像有理函數(shù)模型

    由于具有傳感器參數(shù)保密性、成像模型簡(jiǎn)單性與通用性的優(yōu)點(diǎn)[19],目前高分辨率衛(wèi)星影像的成像模型一般采用有理函數(shù)模型(rational function model,RFM)。RFM是將像點(diǎn)坐標(biāo)(Sample,Line)表示為以相應(yīng)地面點(diǎn)空間坐標(biāo)(P,L,H)(其中P為緯度,L為經(jīng)度,H為高程)為自變量的多項(xiàng)式的比值

    (1)

    式中

    (2)

    式中,b0和d0通常為1;(Pn,Ln,Hn)為正則化的地面坐標(biāo);(cn,rn)為正則化的影像坐標(biāo)

    (3)

    式中,LINE_OFF、SAMP_OFF、LONG_OFF、LAT_OFF、HEIGHT_OFF為正則化平移參數(shù);LINE_SCALE、SAMP_SCALE、LONG_SCALE、LAT_SCALE、HEIGHT_SCALE為正則化比例參數(shù)。上述式(2)中多項(xiàng)式的系數(shù)ai、bi、ci、di(0≤i≤19)和式(3)中的正則化系數(shù)稱為有理多項(xiàng)式系數(shù)(rational polynomial coefficient,RPC)。

    綜合式(1)與式(3),有理函數(shù)模型建立的像素坐標(biāo)(Sample,Line)與空間坐標(biāo)(P,L,H)的關(guān)系式如式(4)所示

    (4)

    2.2 模型法區(qū)域網(wǎng)平差

    一般情況下,衛(wèi)星影像的RPC參數(shù)是由在軌幾何標(biāo)定后的軌道和姿態(tài)數(shù)據(jù)擬合計(jì)算獲得,雖然衛(wèi)星軌道和姿態(tài)數(shù)據(jù)經(jīng)過(guò)在軌幾何標(biāo)定后精度有了很大提高,但是依然會(huì)殘留一定的系統(tǒng)誤差,所以由此擬合得到的RPC參數(shù)往往也具有系統(tǒng)誤差。研究表明,該系統(tǒng)誤差可通過(guò)定義在像方的一個(gè)仿射變換來(lái)進(jìn)一步改正[20]

    (5)

    式中,Line和Sample是連接點(diǎn)或控制點(diǎn)的像方量測(cè)坐標(biāo);Line′和Sample′則是根據(jù)式(4)計(jì)算的連接點(diǎn)或控制點(diǎn)的投影坐標(biāo);a0、as、aL、b0、bs、bL為像方仿射變換參數(shù)。

    衛(wèi)星影像立體區(qū)域網(wǎng)平差的目的在于:①消除衛(wèi)星影像立體像對(duì)內(nèi)部各影像之間的上下視差;②使得不同立體像對(duì)之間的同名觀測(cè)具有相同的空間坐標(biāo);③在有控的條件下,平差解算的加密點(diǎn)空間坐標(biāo)與其控制點(diǎn)空間坐標(biāo)差異性最小。也就是通過(guò)平差調(diào)整衛(wèi)星影像的定向參數(shù),使得同名光線相交至一致位置,且與實(shí)際地面的差異性最小。

    與文獻(xiàn)[20]提出的衛(wèi)星影像光束法區(qū)域網(wǎng)平差方法不同,本文提出衛(wèi)星影像模型法立體區(qū)域網(wǎng)平差方法。其基本原理為:以衛(wèi)星立體像對(duì)(模型)為單元,在原始RPC的基礎(chǔ)上,通過(guò)空間前方交會(huì)計(jì)算連接點(diǎn)在所有立體像對(duì)中的空間坐標(biāo),并將空間坐標(biāo)的均值作為其空間加密坐標(biāo),然后使用式(5)逐一計(jì)算每張影像的像方仿射變換參數(shù);重復(fù)以上過(guò)程,直至前后兩次迭代的像點(diǎn)中誤差變化量小于一定的閾值為止。其原理示意如圖2所示。

    圖2 衛(wèi)星影像模型法區(qū)域網(wǎng)平差原理示意Fig.2 The principle of block adjustment based on stereo model of satellite images

    圖2中3條紅線和黃線分別表示兩個(gè)三視衛(wèi)星立體像對(duì)的3張影像,分別組成模型1和模型2。在平差之前,對(duì)于兩個(gè)立體像對(duì)上的一個(gè)連接點(diǎn)利用原始RPC參數(shù)進(jìn)行空間前方交會(huì),由于各影像的RPC參數(shù)具有一定系統(tǒng)誤差,導(dǎo)致模型1和模型2中像點(diǎn)的對(duì)應(yīng)光線不能交會(huì)至一致位置。如圖2(a)所示,連接點(diǎn)在模型1和模型2中,最優(yōu)交會(huì)空間點(diǎn)分別為P1和P2兩點(diǎn)。平差過(guò)程中,計(jì)算P1和P2兩點(diǎn)的平均位置P,并將其作為控制點(diǎn),根據(jù)式(5)調(diào)整每張影像的定向參數(shù)(即像方的仿射變換參數(shù),相當(dāng)于空間后方交會(huì)),從而使得同名光線相交至點(diǎn)P所在位置。

    衛(wèi)星影像模型法立體區(qū)域網(wǎng)平差的流程如圖3 所示。

    該流程可概括為“空間前方交會(huì)”—“均值化連接點(diǎn)地面坐標(biāo)”—“單片定向”3個(gè)關(guān)鍵性步驟,具體如下:

    (1) 空間前方交會(huì):以立體模型為單位,使用影像的原始RPC參數(shù)和仿射變換改正參數(shù)(每張影像的像方仿射變換改正參數(shù)初值為0),對(duì)連接點(diǎn)進(jìn)行空間前方交會(huì),計(jì)算連接點(diǎn)在每個(gè)模型中的空間坐標(biāo)(本文稱為模型坐標(biāo))。將式(3)代入式(1),然后將其按照泰勒公式展開至一次項(xiàng),得到式(6)

    圖3 衛(wèi)星影像模型法區(qū)域網(wǎng)平差流程Fig.3 The workflow of block adjustment based on stereo model of satellite images

    (6)

    于是,空間前方交會(huì)的誤差方程式為式(7),利用多個(gè)像點(diǎn)觀測(cè)聯(lián)立方程式(7),即可解算像點(diǎn)的物方空間坐標(biāo)

    (7)

    (2) 均值化連接點(diǎn)地面坐標(biāo):由于不同立體模型的直接對(duì)地定位精度不同,所以步驟(1)計(jì)算的同一連接點(diǎn)在不同模型中的前方交會(huì)坐標(biāo)各不相同。通過(guò)計(jì)算連接點(diǎn)在所有立體模型的模型坐標(biāo)均值,并將其作為連接點(diǎn)的地面坐標(biāo),從而逐步消除模型坐標(biāo)的不一致。在有控的情況下,通過(guò)空間相似變換,將連接點(diǎn)地面坐標(biāo)轉(zhuǎn)換為絕對(duì)地面坐標(biāo)。由于測(cè)區(qū)范圍較大,空間相似變換使用地心地固坐標(biāo)。

    (3) 單片定向:根據(jù)式(5)可得單片定向的誤差方程式(8),將步驟(2)得到的連接點(diǎn)作為“虛擬控制點(diǎn)”,按式(8)逐影像進(jìn)行單片定向,計(jì)算每張影像的仿射變換改正參數(shù)

    (8)

    式中,ΔL和ΔS為像點(diǎn)殘差。

    采用迭代計(jì)算,每次迭代結(jié)束后更新影像的仿射變換參數(shù)。計(jì)算相鄰兩次像點(diǎn)中誤差的變化值,若中誤差變化小于預(yù)定閾值,迭代結(jié)束。

    從上述過(guò)程可以看出,本文提出的衛(wèi)星影像模型法區(qū)域網(wǎng)平差的實(shí)質(zhì)為:通過(guò)一致化連接點(diǎn)在各立體模型中前方交會(huì)的地面坐標(biāo),實(shí)現(xiàn)自由網(wǎng)定向(衛(wèi)星影像區(qū)域網(wǎng)內(nèi)部一致性);通過(guò)一致化連接點(diǎn)的地面坐標(biāo)與對(duì)應(yīng)控制點(diǎn)的外業(yè)坐標(biāo),實(shí)現(xiàn)絕對(duì)定向(衛(wèi)星影像區(qū)域網(wǎng)與地面控制點(diǎn)間的外部一致性)。

    2.3 ICESat激光高程點(diǎn)輔助的平差精化

    ICESat激光高程點(diǎn)標(biāo)稱的高程精度為15 cm,該精度是現(xiàn)有國(guó)內(nèi)外高分辨率立體衛(wèi)星影像無(wú)控定位所難以達(dá)到的。天繪一號(hào)01星的無(wú)控定位的平面和高程中誤差分別為10.3 m和5.7 m[21],資源三號(hào)衛(wèi)星三線陣立體影像經(jīng)過(guò)在軌幾何標(biāo)定之后無(wú)控定位的平面和高程中誤差約為15 m和8 m[22-23],GeoEye-1衛(wèi)星GeoStereo級(jí)別立體像對(duì)無(wú)控定位的LE90精度為6 m[24]。因此,ICESat激光高程點(diǎn)可以用來(lái)作為高程控制進(jìn)行輔助衛(wèi)星影像的區(qū)域網(wǎng)平差。

    ICESat激光高程點(diǎn)輔助的平差精化的核心是自動(dòng)提取ICESat高程控制點(diǎn)??紤]到ICESat光斑直徑約70 m,平面精度為10.6±4.5 m[25],當(dāng)衛(wèi)星影像(自由網(wǎng)平差后)平面精度較高(優(yōu)于15 m),可認(rèn)為激光點(diǎn)與影像區(qū)域網(wǎng)平面位置套合較好,此時(shí),在平坦區(qū)域,由于高程變化小,可將激光點(diǎn)高程賦予其光斑范圍內(nèi)空三加密點(diǎn),作為高程控制條件。本文提取ICESat高程控制點(diǎn)的原理如圖4所示。

    對(duì)于ICESat激光高程點(diǎn)Pi,其地面坐標(biāo)為(Xi,Yi,Zi),首先,將其投影至基準(zhǔn)衛(wèi)星影像的pi點(diǎn),由于點(diǎn)pi不一定為顯著特征點(diǎn),因此以點(diǎn)pi為中心,在ICESat激光高程點(diǎn)光斑對(duì)應(yīng)的影像范圍內(nèi)提取影像特征點(diǎn)pf;其次,通過(guò)影像匹配自動(dòng)獲取特征點(diǎn)pf在其他影像上的同名像點(diǎn),通過(guò)多張影像的前方交會(huì)計(jì)算該同名像點(diǎn)對(duì)應(yīng)地面點(diǎn)Pf的坐標(biāo)(Xf,Yf,Zf);最后,使用ICESat激光高程點(diǎn)Pi的高程取代點(diǎn)Pf的高程,形成地面控制點(diǎn)Pi_control(Xf,Yf,Zi)。

    根據(jù)以上方法依次處理測(cè)區(qū)范圍內(nèi)的所有ICESat激光高程點(diǎn),得到多個(gè)ICESat高程控制點(diǎn),重新進(jìn)行有控條件下的區(qū)域網(wǎng)平差,計(jì)算ICESat激光高程點(diǎn)輔助條件下的影像定向參數(shù)。

    3 試驗(yàn)數(shù)據(jù)與結(jié)果分析

    3.1 試驗(yàn)數(shù)據(jù)介紹

    天繪一號(hào)是我國(guó)第一代傳輸型立體測(cè)繪衛(wèi)星,主要用于科學(xué)研究、國(guó)土資源普查、地圖測(cè)繪等領(lǐng)域的科學(xué)試驗(yàn)任務(wù)。天繪一號(hào)01星、02星、03星分別于2010年8月24日、2012年5月6日和2015年10月26日發(fā)射成功,現(xiàn)已組網(wǎng)運(yùn)行。天繪一號(hào)衛(wèi)星采用GPS進(jìn)行定軌,其中01星與02星采用單頻GPS,定軌精度為2~3 m,03星采用雙頻GPS,定軌精度優(yōu)于1 m。

    本文試驗(yàn)數(shù)據(jù)為天繪一號(hào)衛(wèi)星三線陣立體(正視與前后視之間的交會(huì)角為25°)影像數(shù)據(jù)和ICESat/GLAS激光測(cè)高數(shù)據(jù),其中影像地面分辨率為5 m,共2166張(包含01星、02星、03星3顆衛(wèi)星的數(shù)據(jù),初始定位精度不一致),組成了722個(gè)立體模型。影像攝自2010年至2016年間,涉及不同時(shí)相與不同季節(jié),平均同一個(gè)地方有3次覆蓋,地面范圍覆蓋整個(gè)山東省(約15.8萬(wàn)km2)。試驗(yàn)區(qū)域地形復(fù)雜,中部山地突起,東部主要為丘陵,西部主要為平原。為評(píng)定本文平差方法處理后影像的絕對(duì)定位精度,通過(guò)外業(yè)GPS測(cè)量了96個(gè)高精度的野外控制點(diǎn)(平面精度與高程精度均優(yōu)于0.1 m),作為檢查點(diǎn)使用。

    天繪一號(hào)衛(wèi)星影像和ICESat/GLAS激光高程控制點(diǎn)的地理位置如圖5所示,其中藍(lán)色方框表示影像中心,黑色實(shí)心點(diǎn)表示激光高程控制點(diǎn)。

    3.2 試驗(yàn)方案

    (1) 自動(dòng)匹配連接點(diǎn):采用多級(jí)金字塔物方連接點(diǎn)匹配方法[26]匹配影像間的連接點(diǎn),在匹配過(guò)程中使用SRTM數(shù)據(jù)進(jìn)行輔助。SRTM的作用為:首先,用于改正地形引起的影像幾何變形,提高匹配的精度和成功率;其次,用于預(yù)測(cè)匹配點(diǎn)的初始位置,對(duì)于兩張待匹配影像中基準(zhǔn)影像上的像點(diǎn),將其前方交會(huì)至SRTM獲得地面點(diǎn)坐標(biāo),然后再投影至匹配影像,投影位置即為匹配位置的初值。

    (2) 無(wú)控自由網(wǎng)平差:使用自動(dòng)匹配獲得連接點(diǎn),進(jìn)行無(wú)控制條件下的模型法區(qū)域網(wǎng)平差,計(jì)算影像像方仿射變換改正參數(shù)和連接點(diǎn)的地面坐標(biāo)(即空三加密坐標(biāo))。

    (3) ICESat高程控制點(diǎn)自動(dòng)選取與輔助平差:按前文介紹的ICESat激光高程點(diǎn)輔助平差精化方法,自動(dòng)選擇ICESat高程點(diǎn)并將其高程賦予對(duì)應(yīng)連接點(diǎn),作為控制點(diǎn),進(jìn)行有控條件下的區(qū)域網(wǎng)平差。

    (4) 精度驗(yàn)證:以外業(yè)實(shí)測(cè)控制點(diǎn)作為平差精度評(píng)定的檢查點(diǎn),分別檢查無(wú)控條件和ICESat激光高程點(diǎn)輔助條件下兩種平差模式的高程精度。

    3.3 試驗(yàn)結(jié)果及分析

    3.3.1 ICESat高程控制點(diǎn)輔助區(qū)域網(wǎng)平差試驗(yàn)

    試驗(yàn)1:完全無(wú)控條件下的區(qū)域網(wǎng)平差。

    試驗(yàn)2:ICESat高程控制點(diǎn)輔助的區(qū)域網(wǎng)平差。利用全測(cè)區(qū)內(nèi)自動(dòng)提取的1839個(gè)激光高程控制點(diǎn)輔助區(qū)域網(wǎng)平差。

    以上兩組試驗(yàn)均使用96個(gè)外業(yè)檢查點(diǎn)進(jìn)行精度評(píng)定,兩種情況下檢查點(diǎn)的高程殘差分布分別如圖6與圖7所示,其中X軸為經(jīng)度,Y軸為緯度,單位為度,紅色方向表示高程殘差方向,朝上為正(紅色),朝下為負(fù)(藍(lán)色),箭頭矢量長(zhǎng)度表示高程殘差的大小。

    兩種情況下檢查點(diǎn)的精度比較見(jiàn)表2。

    試驗(yàn)結(jié)果表明:

    (1) ICESat激光高程點(diǎn)輔助區(qū)域網(wǎng)平差可顯著提高高程定位精度。在完全無(wú)控條件下,檢查點(diǎn)的平面精度X、Y方向分別為8.45 m和6.86 m,高程精度為5.88 m;ICESat激光高程點(diǎn)輔助區(qū)域網(wǎng)平差后,檢查點(diǎn)平面精度基本保持不變,高程精度則有顯著改善,提高至2.51 m。

    表2完全無(wú)控與ICESat高程點(diǎn)輔助兩種條件下區(qū)域網(wǎng)平差的檢查點(diǎn)精度

    Tab.2 The accuracies of check points after adjustments under two conditions m

    (2) 可有效消除無(wú)控區(qū)域網(wǎng)內(nèi)部扭曲。如圖6所示,在完全無(wú)控情況下,測(cè)區(qū)左右兩端檢查點(diǎn)高程殘差為負(fù),而中部檢查點(diǎn)高程殘差為正,且殘差值較大,說(shuō)明在完全無(wú)控情況下,區(qū)域網(wǎng)內(nèi)部具有一定程度的扭曲;ICESat激光高程點(diǎn)輔助區(qū)域網(wǎng)平差后,整個(gè)測(cè)區(qū)檢查點(diǎn)的高程殘差大小基本一致,見(jiàn)圖7,且與所處的區(qū)域無(wú)關(guān),由此說(shuō)明區(qū)域網(wǎng)內(nèi)部扭曲問(wèn)題得到有效消除。

    3.3.2 不同ICESat高程控制點(diǎn)布設(shè)方案的平差試驗(yàn)

    采用不同的ICESat激光高程控制點(diǎn)布設(shè)方案,采用23個(gè)GLAS條帶的激光高程點(diǎn)進(jìn)行輔助區(qū)域網(wǎng)平差試驗(yàn)。提取了1839個(gè)激光高程點(diǎn),垂直條帶方向以100 km間隔,選取9個(gè)GLAS條帶,對(duì)于每個(gè)條帶沿軌道方向按不同間隔選取激光高程點(diǎn),共選取5個(gè)布設(shè)方案,其點(diǎn)位分布情況和平差結(jié)果統(tǒng)計(jì)見(jiàn)表3,高程殘差分布見(jiàn)圖8。

    圖4 ICESat高程控制點(diǎn)提取原理示意Fig.4 The principle of the extraction of ICESat laser control points

    圖5 天繪一號(hào)衛(wèi)星影像和ICESat/GLAS激光高程點(diǎn)位分布Fig.5 The distribution of Mapping Satellite-1 images and ICESat/GLAS laser points

    圖6 無(wú)控情況下檢查點(diǎn)高程殘差分布Fig.6 The distribution of vertical errors of check points without ground control points

    圖7 ICESat激光高程控制點(diǎn)輔助條件下檢查點(diǎn)高程殘差分布Fig.7 The distribution of vertical errors of check points with ICESat laser control points

    圖8 5種不同高程控制點(diǎn)布設(shè)方案的高程殘差分布Fig. 8 The distributions of vertical errors of five layout schemes for ICESat laser control points

    ICESat激光高程控制點(diǎn)布設(shè)方案高程控制點(diǎn)個(gè)數(shù)控制點(diǎn)高程中誤差/m檢查點(diǎn)高程中誤差/m1軌道兩端布點(diǎn)183.615.052沿軌方向均勻布點(diǎn)(80km間隔)463.064.143沿軌方向均勻布點(diǎn)(40km間隔)782.394.144沿軌方向均勻布點(diǎn)(20km間隔)1602.163.465沿軌方向均勻布點(diǎn)(20km間隔)+所有山區(qū)點(diǎn)2371.923.18

    從上述對(duì)比試驗(yàn)可以得出以下結(jié)論:

    (1) 在軌道兩端布點(diǎn)時(shí),全測(cè)區(qū)高程控制點(diǎn)僅分布在測(cè)區(qū)邊緣,雖然檢查點(diǎn)的高程精度從5.88 m提升至5.05 m,但是提升幅度不大,且由于測(cè)區(qū)內(nèi)部沒(méi)有高程控制點(diǎn),平差后仍然會(huì)出現(xiàn)局部扭曲(圖8(a)中部偏左處紅色箭頭區(qū)域)。

    (2) 隨著ICESat激光高程點(diǎn)數(shù)量增多,檢查點(diǎn)的高程精度將會(huì)提高(如圖8(b)、(c)、(d)所示),而且區(qū)域網(wǎng)的內(nèi)部扭曲也會(huì)得到抑制;當(dāng)采用沿軌方向以20 km間隔均勻布點(diǎn)時(shí),區(qū)域網(wǎng)內(nèi)部扭曲基本消除。然而,以上選點(diǎn)方案所能達(dá)到的高程精度都不及使用全測(cè)區(qū)所有的1839個(gè)激光高程點(diǎn)進(jìn)行整體平差的精度,因此,建議在實(shí)際生產(chǎn)中使用滿足質(zhì)量要求的所有ICESat激光高程點(diǎn)。

    (3) 增加山區(qū)激光高程控制點(diǎn)對(duì)高程精度有較大改善作用。方案5在方案4基礎(chǔ)上增加了位于山區(qū)的激光高程點(diǎn),檢查點(diǎn)的高程精度由3.46 m提升至3.18 m,且區(qū)域網(wǎng)的局部扭曲得到進(jìn)一步抑制,圖8(d)中部偏下區(qū)域的最大正誤差(最長(zhǎng)的紅色箭頭)在圖8(e)中明顯變小。

    4 結(jié) 論

    為了改善衛(wèi)星影像無(wú)控區(qū)域網(wǎng)平差的高程精度,本文使用ICESat/GLAS激光高程點(diǎn)作為輔助,提出了一種ICESat激光高程點(diǎn)輔助的衛(wèi)星影像模型法立體區(qū)域網(wǎng)平差方法,并利用山東全省的天繪一號(hào)衛(wèi)星三線陣影像進(jìn)行了多組區(qū)域網(wǎng)平差試驗(yàn),試驗(yàn)表明:①無(wú)控情況下,檢查點(diǎn)高程精度為5.88 m;利用測(cè)區(qū)內(nèi)自動(dòng)提取的1839個(gè)ICESat高程控制點(diǎn)輔助區(qū)域網(wǎng)平差時(shí),高程精度可提升至2.51 m;②在測(cè)區(qū)內(nèi)選擇布設(shè)ICESat高程控制點(diǎn)時(shí),平坦地區(qū)沿軌道方向間隔80 km即可,而在山區(qū)需要布設(shè)盡可能多的ICESat高程控制點(diǎn)。

    本文方法對(duì)于未來(lái)改善國(guó)產(chǎn)高分辨率立體測(cè)繪衛(wèi)星無(wú)控測(cè)圖精度和提高生產(chǎn)效率具有參考價(jià)值。然而,目前將激光測(cè)高數(shù)據(jù)用于高分辨率衛(wèi)星影像區(qū)域網(wǎng)平差仍然存在著一些問(wèn)題,例如,ICESat與天繪一號(hào)衛(wèi)星軌道不一致,不能保證參與平差的激光高程控制點(diǎn)均勻分布。隨著未來(lái)ICESat-2衛(wèi)星和高分七號(hào)(GF-7)衛(wèi)星的發(fā)射,我們將有機(jī)會(huì)使用更多的激光測(cè)高數(shù)據(jù)和國(guó)產(chǎn)高分辨率立體測(cè)繪衛(wèi)星數(shù)據(jù)進(jìn)行試驗(yàn),以進(jìn)一步驗(yàn)證和提高本文方法的穩(wěn)定性。

    [1] 劉昆波.提高天繪衛(wèi)星影像無(wú)控區(qū)域網(wǎng)平差定位精度的研究[D].武漢:武漢大學(xué),2015.

    LIU Kunbo.The Research of Improving the Positioning Accuracy of Block Adjustment without GCPs for Images from Mapping Satellite-I[D].Wuhan:Wuhan University,2015.

    [2] 王密,楊博,李德仁,等.資源三號(hào)全國(guó)無(wú)控制整體區(qū)域網(wǎng)平差關(guān)鍵技術(shù)及應(yīng)用[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2017,42(4):427-433.

    WANG Mi,YANG Bo,LI Deren,et al.Technologies and Applications of Block Adjustment Without Control for ZY-3 Images Covering China[J].Geomatics and Information Science of Wuhan University,2017,42(4):427-433.

    [3] LI Rongxing,DESHPANDE S,NIU Xutong,et al.Geometric Integration of Aerial and High-Resolution Satellite Imagery and Application in Shoreline Mapping[J].Marine Geodesy,2008,31(3):143-159.

    [4] ZHANG Zuxun,LU Luping,TAO Pengjie,et al.Registration of CBERS-02B Satellite Imagery in Quick GIS Updating[C]∥Proceedings of SPIE 8006,MIPPR 2011:Remote Sensing Image Processing,Geographic Information Systems,and Other Applications.Guilin:SPIE,2011,8006:80060C.

    [5] 黃海蘭,王正濤,金濤勇,等.利用ICESat激光測(cè)高數(shù)據(jù)確定極地冰蓋高程變化[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2012,37(10):1221-1223.

    HUANG Hailan,WANG Zhengtao,JIN Taoyong,et al.Determination of Polar Ice Sheet Height Change from ICESat Altimetry Data[J].Geomatics and Information Science of Wuhan University,2012,37(10):1221-1223.

    [6] 李建成,范春波,褚永海,等.ICESAT衛(wèi)星確定南極冰蓋高程模型研究[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2008,33(3):226-228.

    LI Jiancheng,F(xiàn)AN Chunbo,CHU Yonghai,et al.Using ICESAT Altimeter Data to Determine the Antarctic Ice Sheet Elevation Model[J].Geomatics and Information Science of Wuhan University,2008,33(3):226-228.

    [7] 鄂棟臣,沈強(qiáng),徐瑩,等.基于ASTER立體數(shù)據(jù)和ICESat/GLAS測(cè)高數(shù)據(jù)融合高精度提取南極地區(qū)地形信息[J].中國(guó)科學(xué) D輯:地球科學(xué),2009,39(3):351-359.

    E Dongchen,SHEN Qiang,XU Ying,et al.High-accuracy Topographical Information Extraction Based on Fusion of ASTER Stereo-data and ICESat/GLAS Data in Antarctica[J].Science in China Series D-Earth Sciences,2009,52(5):714-722.

    [8] WANG Xianwen,CHENG Xiao,GONG Peng,et al.Earth Science Applications of ICESat/GLAS:A Review[J].International Journal of Remote Sensing,2011,32(23):8837-8864.

    [9] ZWALLY H J,SCHUTZ B,ABDALATI W,et al.ICESat’s Laser Measurements of Polar Ice,Atmosphere,Ocean,and Land[J].Journal of Geodynamics,2002,34(3-4):405-445.

    [10] GONZALEZ J H,BACHMANN M,SCHEIBER R,et al.Definition of ICESat Selection Criteria for Their Use as Height References for TanDEM-X[J].IEEE Transactions on Geoscience and Remote Sensing,2010,48(6):2750-2757.

    [11] 萬(wàn)杰,廖靜娟,許濤,等.基于ICESat/GLAS高度計(jì)數(shù)據(jù)的SRTM數(shù)據(jù)精度評(píng)估——以青藏高原地區(qū)為例[J].國(guó)土資源遙感,2015,27(1):100-105.

    WAN Jie,LIAO Jingjuan,XU Tao,et al.Accuracy Evaluation of SRTM Data based on ICESat/GLAS Altimeter Data:A Case Study in the Tibetan Plateau[J].Remote Sensing for Land & Resources,2015,27(1):100-105.

    [12] 李國(guó)元,唐新明,張重陽(yáng),等.多準(zhǔn)則約束的ICESat/GLAS高程控制點(diǎn)篩選[J].遙感學(xué)報(bào),2017,21(1):96-104.

    LI Guoyuan,TANG Xinming,ZHANG Chongyang,et al.Multi-criteria Constraint Algorithm for Selecting ICESat/GLAS Data as Elevation Control Points[J].Journal of Remote Sensing,2017,21(1):96-104.

    [13] 王任享,王建榮.二線陣CCD衛(wèi)星影像聯(lián)合激光測(cè)距數(shù)據(jù)光束法平差技術(shù)[J].測(cè)繪科學(xué)技術(shù)學(xué)報(bào),2014,31(1):1-4.

    WANG Renxiang,WANG Jianrong.Technology of Bundle Adjustment Using Two-line-array CCD Satellite Image Combined Laser Ranging Data[J].Journal of Geomatics Science and Technology,2014,31(1):1-4.

    [14] 王建榮,王任享,胡莘,等.利用激光測(cè)距數(shù)據(jù)處理線陣衛(wèi)星攝影測(cè)量影像[J].測(cè)繪科學(xué),2013,38(2):15-16.

    WANG Jianrong,WANG Renxiang,HU Xin,et al.Application of Laser Distance Measurement Data in Linear Array Satellite Photogrammetry[J].Science of Surveying and Mapping,2013,38(2):15-16.

    [15] WU Bo,GUO Jian,HU Han,et al.Co-registration of Lunar Topographic Models Derived from Chang’E-1,SELENE,and LRO Laser Altimeter Data Based on a Novel Surface Matching Method[J].Earth and Planetary Science Letters,2013,364:68-84.

    [16] 唐新明,李國(guó)元,高小明,等.衛(wèi)星激光測(cè)高嚴(yán)密幾何模型構(gòu)建及精度初步驗(yàn)證[J].測(cè)繪學(xué)報(bào),2016,45(10):1182-1191.DOI:10.11947/j.AGCS.2016.20150357.

    TANG Xinming,LI Guoyuan,GAO Xiaoming,et al.The Rigorous Geometric Model of Satellite Laser Altimeter and Preliminarily Accuracy Validation[J].Acta Geodaetica et Cartographica Sinica,2016,45(10):1182-1191.DOI:10.11947/j.AGCS.2016.20150357.

    [17] 李國(guó)元,唐新明,王華斌,等.GLAS激光測(cè)高數(shù)據(jù)輔助的資源三號(hào)三線陣區(qū)域網(wǎng)平差研究[C]∥第三屆高分辨率對(duì)地觀測(cè)學(xué)術(shù)年會(huì).北京:高分辨率對(duì)地觀測(cè)系統(tǒng)重大專項(xiàng)管理辦公室,2014:586-600.

    LI Guoyuan,TANG Xinming,WANG Huabin,et al.Research on the ZY-3 Block Adjustment Supported by the GLAS Laser Altimetry Data[C]∥Proceedings of the 3rd China High Resolution Earth Observation Conference.Beijing:Major Special Management Office of High Resolution to Earth Observation System,2014:586-600.

    [18] HU Zhihua,PENG Jianwei,HOU Yaolin,et al.Evaluation of Recently Released Open Global Digital Elevation Models of Hubei,China[J].Remote Sensing,2017,9(3):262.

    [19] TAO C V,HU Y.A Comprehensive Study of the Rational Function Model for Photogrammetric Processing[J].Photogrammetric Engineering & Remote Sensing,2001,67(12):1347-1357.

    [20] GRODECKI J,DIAL G.Block Adjustment of High-resolution Satellite Images Described by Rational Polynomials[J].Photogrammetric Engineering & Remote Sensing,2003,69(1):59-68.

    [21] 付勇,鄒松柏,劉會(huì)安.“天繪一號(hào)”01星立體影像定位精度檢測(cè)[J].遙感學(xué)報(bào),2012(S1):94-97.

    FU Yong,ZOU Songbai,LIU Hui’an.Evaluation of the Location Accuracy of the Mapping Satellite-1 Stereo Image[J].Journal of Remote Sensing,2012(S1):94-97.

    [22] 潘紅播,張過(guò),唐新明,等.資源三號(hào)測(cè)繪衛(wèi)星影像產(chǎn)品精度分析與驗(yàn)證[J].測(cè)繪學(xué)報(bào),2013,42(5):738-744,751.

    PAN Hongbo,ZHANG Guo,TANG Xinming,et al.Accuracy analysis and verification of ZY-3 Products[J].Acta Geodaetica et Cartographica Sinica,2013,42(5):738-744,751.

    [23] 蔣永華,張過(guò),唐新明,等.資源三號(hào)測(cè)繪衛(wèi)星三線陣影像高精度幾何檢校[J].測(cè)繪學(xué)報(bào),2013,42(4):523-529,553.

    JIANG Yonghua,ZHANG Guo,TANG Xinming,et al.High Accuracy Geometric Calibration of ZY-3 Three-Line Image[J].Acta Geodaetica et Cartographica Sinica,2013,42(4):523-529,553.

    [24] GeoEye I.GeoEye Product Guide[EB/OL].[2011-12-21].http:∥www.geoeye.com/CorpSite/assets/docs/brochures/ GeoEye_Product_Guide.pdf.

    [25] MAGRUDER L A,WEBB C E,URBAN T J,et al.ICESat Altimetry Data Product Verification at White Sands Space Harbor[J].IEEE Transactions on Geoscience and Remote Sensing,2007,45(1):147-155.

    [26] TAO Pengjie,LU Luping,ZHANG Yong,et al.On-Orbit Geometric Calibration of the Panchromatic/Multispectral Camera of the ZY-1 02C Satellite Based on Public Geographic Data[J].Photogrammetric Engineering & Remote Sensing,2014,80(6):505-517.

    猜你喜歡
    網(wǎng)平控制點(diǎn)立體
    兩種大區(qū)域GNSS網(wǎng)平差方法對(duì)比分析
    念個(gè)立體咒
    立體登陸
    NFFD控制點(diǎn)分布對(duì)氣動(dòng)外形優(yōu)化的影響
    利用RFM的天繪一號(hào)衛(wèi)星影像區(qū)域網(wǎng)平差
    基于風(fēng)險(xiǎn)管理下的項(xiàng)目建設(shè)內(nèi)部控制點(diǎn)思考
    炫酷立體卡
    相似材料模型中控制點(diǎn)像點(diǎn)坐標(biāo)定位研究
    不同分辨率線陣衛(wèi)星影像聯(lián)合區(qū)域網(wǎng)平差
    遙感信息(2015年3期)2015-12-13 07:26:50
    SDCORS在基礎(chǔ)地理信息控制點(diǎn)補(bǔ)測(cè)中的應(yīng)用
    国产片特级美女逼逼视频| 99精国产麻豆久久婷婷| 亚洲欧美清纯卡通| 可以免费在线观看a视频的电影网站| 80岁老熟妇乱子伦牲交| videos熟女内射| 悠悠久久av| 精品少妇一区二区三区视频日本电影| 三上悠亚av全集在线观看| tube8黄色片| 一区福利在线观看| av视频免费观看在线观看| 欧美成人精品欧美一级黄| 在线 av 中文字幕| 青青草视频在线视频观看| av网站在线播放免费| 国产成人免费观看mmmm| 丝袜美腿诱惑在线| 欧美人与善性xxx| 日日夜夜操网爽| 免费女性裸体啪啪无遮挡网站| 日韩免费高清中文字幕av| 亚洲色图 男人天堂 中文字幕| 一边摸一边做爽爽视频免费| 国产男女内射视频| 一区二区三区精品91| 天堂中文最新版在线下载| 1024香蕉在线观看| 国产精品偷伦视频观看了| 国产一卡二卡三卡精品| 赤兔流量卡办理| 亚洲专区中文字幕在线| 亚洲中文字幕日韩| 免费在线观看视频国产中文字幕亚洲 | 久久精品aⅴ一区二区三区四区| 热99久久久久精品小说推荐| 国产精品九九99| 国产福利在线免费观看视频| 精品福利永久在线观看| 满18在线观看网站| 啦啦啦 在线观看视频| 午夜免费观看性视频| 国产成人系列免费观看| 国产xxxxx性猛交| 欧美av亚洲av综合av国产av| 99久久人妻综合| 深夜精品福利| 这个男人来自地球电影免费观看| 少妇的丰满在线观看| 免费av中文字幕在线| 高清欧美精品videossex| 婷婷色av中文字幕| 热99久久久久精品小说推荐| 日韩人妻精品一区2区三区| 国产一区亚洲一区在线观看| 亚洲色图综合在线观看| 一区二区三区激情视频| 国产精品久久久人人做人人爽| 肉色欧美久久久久久久蜜桃| 亚洲欧洲国产日韩| 欧美日韩亚洲国产一区二区在线观看 | 午夜av观看不卡| 国产欧美日韩综合在线一区二区| 国产视频首页在线观看| 久久精品国产综合久久久| 国产精品久久久久久精品电影小说| 午夜精品国产一区二区电影| 亚洲人成网站在线观看播放| 国产精品久久久久久人妻精品电影 | 欧美成人精品欧美一级黄| 日韩精品免费视频一区二区三区| 国产女主播在线喷水免费视频网站| 少妇人妻久久综合中文| 美女福利国产在线| 好男人电影高清在线观看| 国产精品二区激情视频| 十八禁人妻一区二区| 2018国产大陆天天弄谢| 不卡av一区二区三区| 97在线人人人人妻| 国产av精品麻豆| 久久这里只有精品19| av在线播放精品| 亚洲中文字幕日韩| 欧美老熟妇乱子伦牲交| 最新的欧美精品一区二区| 国产精品免费大片| 一边摸一边抽搐一进一出视频| 日本91视频免费播放| 日韩熟女老妇一区二区性免费视频| 久久精品亚洲av国产电影网| 亚洲av电影在线观看一区二区三区| 一二三四在线观看免费中文在| 午夜精品国产一区二区电影| 免费黄频网站在线观看国产| 91成人精品电影| 水蜜桃什么品种好| 不卡av一区二区三区| 99久久99久久久精品蜜桃| 欧美日韩国产mv在线观看视频| 午夜免费鲁丝| 国产成人av激情在线播放| www.精华液| 男女边吃奶边做爰视频| 精品一区在线观看国产| 人人澡人人妻人| 亚洲欧美一区二区三区久久| 中文字幕高清在线视频| 免费在线观看日本一区| 亚洲精品av麻豆狂野| 9191精品国产免费久久| 侵犯人妻中文字幕一二三四区| 人人妻,人人澡人人爽秒播 | 99re6热这里在线精品视频| 久久 成人 亚洲| 久久久久精品国产欧美久久久 | 欧美 日韩 精品 国产| 免费看av在线观看网站| 日韩制服丝袜自拍偷拍| 日韩电影二区| 国产又色又爽无遮挡免| 午夜影院在线不卡| 国产真人三级小视频在线观看| 国产主播在线观看一区二区 | 亚洲欧洲国产日韩| 国产精品亚洲av一区麻豆| 最新在线观看一区二区三区 | 免费一级毛片在线播放高清视频 | 亚洲,一卡二卡三卡| 免费久久久久久久精品成人欧美视频| 黄色a级毛片大全视频| 久久精品久久精品一区二区三区| 99国产精品一区二区蜜桃av | 高清不卡的av网站| 久9热在线精品视频| 日韩制服丝袜自拍偷拍| 老汉色∧v一级毛片| 国产一区有黄有色的免费视频| 国产91精品成人一区二区三区 | 黄色毛片三级朝国网站| 亚洲欧美激情在线| 亚洲精品美女久久av网站| 免费在线观看黄色视频的| 99久久综合免费| 精品亚洲乱码少妇综合久久| 国产伦人伦偷精品视频| 黑人猛操日本美女一级片| 黄网站色视频无遮挡免费观看| 深夜精品福利| 免费不卡黄色视频| 男人爽女人下面视频在线观看| 久久久精品免费免费高清| 久久精品成人免费网站| 国精品久久久久久国模美| 亚洲图色成人| 欧美xxⅹ黑人| 夫妻午夜视频| 久久久久久久大尺度免费视频| 国产亚洲欧美在线一区二区| 中文字幕色久视频| 亚洲精品一卡2卡三卡4卡5卡 | 精品人妻在线不人妻| 肉色欧美久久久久久久蜜桃| 国产成人91sexporn| 91精品国产国语对白视频| 晚上一个人看的免费电影| 国产在线观看jvid| 亚洲欧美精品自产自拍| 久久午夜综合久久蜜桃| 一级黄片播放器| 国产在线视频一区二区| 大陆偷拍与自拍| 精品国产乱码久久久久久小说| 午夜日韩欧美国产| 日本五十路高清| 免费女性裸体啪啪无遮挡网站| 日本vs欧美在线观看视频| 99国产精品一区二区三区| 午夜精品国产一区二区电影| 男女之事视频高清在线观看 | 亚洲一区中文字幕在线| 久久久精品免费免费高清| 国产免费福利视频在线观看| 精品一区在线观看国产| 一区二区三区激情视频| 黄色一级大片看看| 久久精品久久久久久噜噜老黄| 91九色精品人成在线观看| 亚洲 国产 在线| 免费女性裸体啪啪无遮挡网站| 老司机午夜十八禁免费视频| 欧美性长视频在线观看| 搡老岳熟女国产| 亚洲精品日韩在线中文字幕| 大香蕉久久成人网| 亚洲中文日韩欧美视频| 成年av动漫网址| 午夜福利免费观看在线| 大码成人一级视频| 王馨瑶露胸无遮挡在线观看| 制服人妻中文乱码| 免费看av在线观看网站| 国产有黄有色有爽视频| 在线av久久热| 亚洲中文日韩欧美视频| 亚洲精品中文字幕在线视频| 久久久国产欧美日韩av| 人体艺术视频欧美日本| 免费观看a级毛片全部| 亚洲精品久久午夜乱码| 日本欧美国产在线视频| 亚洲国产精品成人久久小说| 涩涩av久久男人的天堂| 欧美日本中文国产一区发布| 久久免费观看电影| 狂野欧美激情性xxxx| 国产精品99久久99久久久不卡| 日韩,欧美,国产一区二区三区| 久久久久视频综合| 亚洲精品中文字幕在线视频| 国产女主播在线喷水免费视频网站| 久久人人爽av亚洲精品天堂| 日本五十路高清| 色精品久久人妻99蜜桃| 亚洲成人手机| 久久毛片免费看一区二区三区| 一区在线观看完整版| 2021少妇久久久久久久久久久| 免费看不卡的av| 18禁黄网站禁片午夜丰满| 下体分泌物呈黄色| 性高湖久久久久久久久免费观看| 美女扒开内裤让男人捅视频| 91字幕亚洲| 国产成人精品久久二区二区91| 人人妻人人添人人爽欧美一区卜| videos熟女内射| 最新的欧美精品一区二区| 99国产精品99久久久久| 日韩欧美一区视频在线观看| 亚洲国产av新网站| 侵犯人妻中文字幕一二三四区| av有码第一页| 在线观看国产h片| 精品人妻一区二区三区麻豆| 青草久久国产| 男人爽女人下面视频在线观看| 在线天堂中文资源库| 成人国产av品久久久| 国产一区二区激情短视频 | 在线观看一区二区三区激情| 亚洲欧洲日产国产| 七月丁香在线播放| 久久久久精品国产欧美久久久 | 久久亚洲精品不卡| 亚洲国产欧美网| av在线播放精品| 妹子高潮喷水视频| 成人午夜精彩视频在线观看| 亚洲欧美精品自产自拍| 大片电影免费在线观看免费| 夜夜骑夜夜射夜夜干| 尾随美女入室| av又黄又爽大尺度在线免费看| 久久久久久久大尺度免费视频| 亚洲av欧美aⅴ国产| 午夜免费鲁丝| 亚洲精品久久成人aⅴ小说| bbb黄色大片| 国产精品免费视频内射| 国产亚洲精品久久久久5区| 欧美久久黑人一区二区| 国产日韩欧美在线精品| 免费观看a级毛片全部| 亚洲成人手机| 两个人免费观看高清视频| 欧美人与性动交α欧美精品济南到| 国产精品熟女久久久久浪| www.999成人在线观看| 一边亲一边摸免费视频| 妹子高潮喷水视频| 欧美激情高清一区二区三区| 啦啦啦中文免费视频观看日本| 观看av在线不卡| 99国产精品99久久久久| 男女床上黄色一级片免费看| 午夜久久久在线观看| 成人国产一区最新在线观看 | 色视频在线一区二区三区| 无遮挡黄片免费观看| 天天躁日日躁夜夜躁夜夜| 国产真人三级小视频在线观看| 久久精品久久久久久噜噜老黄| 亚洲专区国产一区二区| 免费观看人在逋| 亚洲欧洲日产国产| 亚洲第一av免费看| 一级毛片 在线播放| 五月天丁香电影| 一本一本久久a久久精品综合妖精| 亚洲欧美一区二区三区黑人| 亚洲av电影在线观看一区二区三区| 天天影视国产精品| 一区二区av电影网| 一边摸一边抽搐一进一出视频| 亚洲人成77777在线视频| 国产日韩一区二区三区精品不卡| 2021少妇久久久久久久久久久| 日韩视频在线欧美| 老司机影院毛片| bbb黄色大片| 欧美久久黑人一区二区| 日韩伦理黄色片| 99香蕉大伊视频| 日本猛色少妇xxxxx猛交久久| 亚洲一区中文字幕在线| 国产真人三级小视频在线观看| 满18在线观看网站| 91麻豆av在线| 中文乱码字字幕精品一区二区三区| 日韩欧美一区视频在线观看| 精品久久蜜臀av无| 午夜免费成人在线视频| 亚洲人成电影观看| 亚洲精品第二区| 中文字幕另类日韩欧美亚洲嫩草| 少妇猛男粗大的猛烈进出视频| 亚洲成人免费av在线播放| 欧美人与性动交α欧美精品济南到| 校园人妻丝袜中文字幕| 欧美97在线视频| 少妇人妻久久综合中文| 精品卡一卡二卡四卡免费| 欧美黑人欧美精品刺激| 国产成人精品在线电影| 久久久欧美国产精品| 999精品在线视频| 新久久久久国产一级毛片| 国产在线观看jvid| 久久女婷五月综合色啪小说| 亚洲熟女精品中文字幕| 大码成人一级视频| 久久久久网色| videos熟女内射| 中文字幕人妻丝袜制服| 99久久综合免费| 精品一品国产午夜福利视频| 成人三级做爰电影| 国产精品香港三级国产av潘金莲 | 人妻人人澡人人爽人人| 大香蕉久久网| 91麻豆av在线| 美女中出高潮动态图| 少妇人妻 视频| 午夜日韩欧美国产| 啦啦啦 在线观看视频| 久久热在线av| 国产又色又爽无遮挡免| 久久久久视频综合| 日韩视频在线欧美| 亚洲精品久久久久久婷婷小说| 咕卡用的链子| 国产亚洲精品久久久久5区| 咕卡用的链子| 人妻一区二区av| www.熟女人妻精品国产| 99热全是精品| 久久九九热精品免费| 日韩中文字幕欧美一区二区 | 久久 成人 亚洲| 老汉色av国产亚洲站长工具| 国产精品亚洲av一区麻豆| 国产深夜福利视频在线观看| 日韩一本色道免费dvd| 欧美日韩精品网址| 男女下面插进去视频免费观看| 中国国产av一级| 久久久欧美国产精品| 伊人亚洲综合成人网| 亚洲一区中文字幕在线| 黄频高清免费视频| 久久 成人 亚洲| 亚洲精品美女久久av网站| 青春草亚洲视频在线观看| www.自偷自拍.com| 国产免费福利视频在线观看| 下体分泌物呈黄色| 欧美黑人欧美精品刺激| 高潮久久久久久久久久久不卡| 国产在线观看jvid| 国产一区二区在线观看av| 国产高清国产精品国产三级| 国产欧美日韩一区二区三 | 国产一区二区三区综合在线观看| 麻豆国产av国片精品| 欧美成人精品欧美一级黄| 天天操日日干夜夜撸| 99国产精品一区二区三区| 1024香蕉在线观看| 美女扒开内裤让男人捅视频| 另类亚洲欧美激情| 国产精品秋霞免费鲁丝片| 午夜视频精品福利| 亚洲av日韩在线播放| 久久久久久久久免费视频了| 亚洲久久久国产精品| 爱豆传媒免费全集在线观看| 老司机亚洲免费影院| 高清黄色对白视频在线免费看| 50天的宝宝边吃奶边哭怎么回事| 久久久久久久大尺度免费视频| 中文欧美无线码| 精品久久久精品久久久| 美女脱内裤让男人舔精品视频| 激情五月婷婷亚洲| 国产亚洲一区二区精品| 久久精品成人免费网站| 99九九在线精品视频| 久久女婷五月综合色啪小说| 中国国产av一级| av在线老鸭窝| 女人精品久久久久毛片| 老司机靠b影院| 电影成人av| 青青草视频在线视频观看| 男人添女人高潮全过程视频| 爱豆传媒免费全集在线观看| 亚洲五月色婷婷综合| 欧美激情极品国产一区二区三区| 王馨瑶露胸无遮挡在线观看| 人妻 亚洲 视频| 亚洲欧美精品自产自拍| 久久 成人 亚洲| 波多野结衣av一区二区av| 两人在一起打扑克的视频| 亚洲中文字幕日韩| 久久女婷五月综合色啪小说| 欧美精品av麻豆av| 久热爱精品视频在线9| 一个人免费看片子| 亚洲成色77777| 两性夫妻黄色片| 99热国产这里只有精品6| 丝袜美足系列| 每晚都被弄得嗷嗷叫到高潮| 老司机靠b影院| 19禁男女啪啪无遮挡网站| 国产成人影院久久av| 多毛熟女@视频| 丝袜在线中文字幕| 51午夜福利影视在线观看| 精品一区二区三区四区五区乱码 | 51午夜福利影视在线观看| 久久综合国产亚洲精品| 日本黄色日本黄色录像| 国产一区亚洲一区在线观看| 午夜免费男女啪啪视频观看| 久久ye,这里只有精品| av天堂在线播放| av福利片在线| 久久久国产欧美日韩av| 爱豆传媒免费全集在线观看| 成人国产一区最新在线观看 | 自线自在国产av| 香蕉国产在线看| 日韩制服丝袜自拍偷拍| 国产精品一国产av| 欧美日韩黄片免| 国产一区二区三区av在线| 最近手机中文字幕大全| 母亲3免费完整高清在线观看| 丰满饥渴人妻一区二区三| 国产成人精品久久二区二区91| a级毛片黄视频| 国产精品一区二区在线观看99| 热99国产精品久久久久久7| 久久久久久久久久久久大奶| 中文乱码字字幕精品一区二区三区| 亚洲一码二码三码区别大吗| 午夜日韩欧美国产| 在线看a的网站| 亚洲三区欧美一区| 久久久亚洲精品成人影院| 亚洲精品中文字幕在线视频| 老司机午夜十八禁免费视频| 韩国精品一区二区三区| 首页视频小说图片口味搜索 | 男女高潮啪啪啪动态图| 赤兔流量卡办理| 看免费av毛片| 90打野战视频偷拍视频| 精品福利观看| 婷婷色麻豆天堂久久| 晚上一个人看的免费电影| 国产精品秋霞免费鲁丝片| 精品一区二区三卡| 天堂8中文在线网| 90打野战视频偷拍视频| 成年人午夜在线观看视频| 国产成人影院久久av| 91字幕亚洲| 啦啦啦在线观看免费高清www| 久久久精品免费免费高清| 精品少妇久久久久久888优播| 在线天堂中文资源库| 国产精品久久久久成人av| 最新的欧美精品一区二区| 在线亚洲精品国产二区图片欧美| 国产精品人妻久久久影院| 日韩 欧美 亚洲 中文字幕| 精品第一国产精品| 午夜福利在线免费观看网站| 9热在线视频观看99| 国产色视频综合| 亚洲专区中文字幕在线| 免费久久久久久久精品成人欧美视频| 欧美人与善性xxx| 日韩免费高清中文字幕av| 又大又爽又粗| 在线亚洲精品国产二区图片欧美| 成人三级做爰电影| 欧美日本中文国产一区发布| 婷婷色综合www| 超色免费av| 韩国精品一区二区三区| 中文精品一卡2卡3卡4更新| 午夜福利一区二区在线看| 亚洲国产av影院在线观看| 中文字幕高清在线视频| 国产福利在线免费观看视频| 欧美日韩黄片免| 久久久久久亚洲精品国产蜜桃av| 久久毛片免费看一区二区三区| 啦啦啦在线免费观看视频4| 国产精品二区激情视频| 一本综合久久免费| 视频在线观看一区二区三区| 水蜜桃什么品种好| 99香蕉大伊视频| 国产一区有黄有色的免费视频| 国产一区二区在线观看av| 美女午夜性视频免费| 国产国语露脸激情在线看| 国产成人系列免费观看| 在线av久久热| 国产伦理片在线播放av一区| 国产欧美日韩综合在线一区二区| 欧美成人精品欧美一级黄| 久久久久国产一级毛片高清牌| 亚洲图色成人| 欧美另类一区| 国产精品香港三级国产av潘金莲 | 女人被躁到高潮嗷嗷叫费观| 日韩大码丰满熟妇| 97精品久久久久久久久久精品| 国产精品久久久av美女十八| 9色porny在线观看| 一区二区日韩欧美中文字幕| 国产99久久九九免费精品| 黄片小视频在线播放| 丰满少妇做爰视频| 国产不卡av网站在线观看| 丝袜美足系列| 亚洲成色77777| 亚洲欧美一区二区三区国产| bbb黄色大片| 日本五十路高清| 男女国产视频网站| 18禁国产床啪视频网站| 免费在线观看日本一区| 交换朋友夫妻互换小说| 久久精品久久久久久久性| 日韩制服骚丝袜av| 精品国产一区二区三区四区第35| 午夜激情久久久久久久| 无遮挡黄片免费观看| 欧美 日韩 精品 国产| 亚洲精品在线美女| 2018国产大陆天天弄谢| 在线观看免费视频网站a站| 69精品国产乱码久久久| 五月天丁香电影| 精品人妻一区二区三区麻豆| 黄色一级大片看看| 国产成人啪精品午夜网站| 欧美日韩亚洲高清精品| 国产男女超爽视频在线观看| av天堂久久9| 亚洲精品久久成人aⅴ小说| 国产成人影院久久av| 亚洲午夜精品一区,二区,三区| 国产熟女欧美一区二区| 啦啦啦在线免费观看视频4| 国产午夜精品一二区理论片| 国产97色在线日韩免费| 黑人巨大精品欧美一区二区蜜桃| 国产一区二区三区av在线| 亚洲精品第二区| 亚洲精品自拍成人| 少妇 在线观看| 永久免费av网站大全| 一本久久精品| 人体艺术视频欧美日本| 黑人猛操日本美女一级片| 国产日韩欧美在线精品| 国产极品粉嫩免费观看在线| 欧美日韩国产mv在线观看视频| 人人妻人人爽人人添夜夜欢视频| 男女国产视频网站| 嫩草影视91久久| 国产欧美日韩综合在线一区二区| 午夜日韩欧美国产| 老汉色av国产亚洲站长工具| 免费看十八禁软件| 97精品久久久久久久久久精品|