亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        多視場(chǎng)遙感圖像虛擬焦面拼接理論誤差分析

        2015-10-11 02:22:54岳春宇何紅艷鮑云飛邢坤周楠
        航天返回與遙感 2015年2期
        關(guān)鍵詞:投影坐標(biāo)系高程

        岳春宇 何紅艷 鮑云飛 邢坤 周楠

        (北京空間機(jī)電研究所,北京100094)

        0 引言

        隨著遙感對(duì)地觀測(cè)技術(shù)的發(fā)展,國(guó)土調(diào)查、地圖測(cè)繪等應(yīng)用領(lǐng)域?qū)υ龃筮b感圖像覆蓋范圍的要求不斷提高,因此,擴(kuò)大視場(chǎng)、獲得大寬幅數(shù)據(jù)是空間光學(xué)遙感發(fā)展的一個(gè)重要方向。由于空間相機(jī)受衛(wèi)星平臺(tái)、空間環(huán)境等因素制約,以及CCD器件工藝的限制[1],一般采用多相機(jī)多視場(chǎng)圖像拼接的方式來(lái)獲取滿足較高空間分辨率且大幅寬的圖像[2]。在遙感應(yīng)用領(lǐng)域,目前遙感圖像拼接主要是通過(guò)圖像匹配,尋找同名信息建立圖像間的空間對(duì)應(yīng)關(guān)系,建立空間轉(zhuǎn)換模型,再將遙感圖像轉(zhuǎn)換到同一坐標(biāo)系下[3-5]。這種拼接方式需要在待拼接圖像中找到同名信息,當(dāng)待拼接圖像重疊區(qū)域提取同名特征困難時(shí),則無(wú)法準(zhǔn)確拼接[6]。而且這種方法計(jì)算量大,不同的待拼接圖像處理結(jié)果不同,不具有普遍性,不能夠?qū)崿F(xiàn)標(biāo)準(zhǔn)化生產(chǎn)。數(shù)據(jù)生產(chǎn)單位更希望通過(guò)成像幾何關(guān)系直接對(duì)不同相機(jī)圖像進(jìn)行拼接,這樣避免了復(fù)雜費(fèi)時(shí)的圖像匹配,同時(shí)可以按照統(tǒng)一流程高效生產(chǎn)圖像數(shù)據(jù)產(chǎn)品。

        本文從遙感圖像虛擬焦面拼接原理出發(fā),推導(dǎo)圖像拼接理論誤差公式,研究影響多視場(chǎng)遙感圖像拼接精度的因素,并通過(guò)仿真實(shí)驗(yàn)對(duì)圖像拼接誤差進(jìn)行分析。

        1 多視場(chǎng)遙感圖像拼接

        多視場(chǎng)遙感圖像拼接需要建立數(shù)學(xué)轉(zhuǎn)換模型,將不同相機(jī)圖像通過(guò)坐標(biāo)轉(zhuǎn)換以及重采樣,變換到拼接圖像坐標(biāo)系下。在常規(guī)遙感圖像產(chǎn)品生產(chǎn)過(guò)程中,地面點(diǎn)坐標(biāo)轉(zhuǎn)換到圖像坐標(biāo)系下,需要經(jīng)過(guò)多次坐標(biāo)轉(zhuǎn)換。

        1.1 多視場(chǎng)遙感圖像拼接空間坐標(biāo)系確定

        地面點(diǎn)P由物方坐標(biāo)計(jì)算其圖像坐標(biāo)的流程為[7]:

        1)根據(jù)成像時(shí)刻,從星歷文件中獲得的衛(wèi)星位置及速度矢量參數(shù),計(jì)算地心直角坐標(biāo)系到衛(wèi)星軌道坐標(biāo)系的變換矩陣,將P點(diǎn)物方坐標(biāo)轉(zhuǎn)換到衛(wèi)星軌道坐標(biāo)系下,得到P點(diǎn)在成像時(shí)所對(duì)應(yīng)的衛(wèi)星軌道坐標(biāo)系坐標(biāo);

        2)根據(jù)成像時(shí)刻星歷記載的成像時(shí)衛(wèi)星姿態(tài)參數(shù),計(jì)算衛(wèi)星軌道坐標(biāo)系與衛(wèi)星本體坐標(biāo)系間的正交變換矩陣,將P點(diǎn)的衛(wèi)星軌道坐標(biāo)系坐標(biāo)轉(zhuǎn)換到衛(wèi)星本體坐標(biāo)系,得到P點(diǎn)在衛(wèi)星本體坐標(biāo)下坐標(biāo);

        3)根據(jù)相機(jī)安置矩陣,得到從衛(wèi)星本體坐標(biāo)系到像空間坐標(biāo)系的正交變換矩陣,將P點(diǎn)衛(wèi)星本體坐標(biāo)轉(zhuǎn)換到相機(jī)坐標(biāo)系下,得到P點(diǎn)的像空間坐標(biāo);

        4)根據(jù)攝影比例尺,將P點(diǎn)的像空間坐標(biāo)投影到像平面坐標(biāo),即得到地面點(diǎn)P的圖像坐標(biāo)。

        在上述轉(zhuǎn)換過(guò)程中,第1、2步只與衛(wèi)星平臺(tái)有關(guān),坐標(biāo)轉(zhuǎn)換關(guān)系適用于待拼接的各臺(tái)相機(jī);第3、4步是相機(jī)與衛(wèi)星本體和相機(jī)內(nèi)部之間的坐標(biāo)轉(zhuǎn)換。在建立成像幾何模型時(shí),由于待拼接各臺(tái)相機(jī)第 3、4步中的坐標(biāo)轉(zhuǎn)換矩陣不同,所以本文在討論多視場(chǎng)遙感圖像拼接誤差時(shí),將成像幾何模型簡(jiǎn)化為由共線條件方程表示的衛(wèi)星本體坐標(biāo)系到像空間坐標(biāo)的轉(zhuǎn)換過(guò)程。

        為了方便后續(xù)討論,不考慮軌道坐標(biāo)系與地心坐標(biāo)系之間的轉(zhuǎn)換。對(duì)于兩臺(tái)相機(jī)拼接的情況(如圖1所示),將各相機(jī)攝影中心、地面點(diǎn)統(tǒng)一到一個(gè)物方坐標(biāo)系(O - XYZ)中,將各圖像統(tǒng)一到一個(gè)像空間坐標(biāo)系(O - x yz)中,則物方坐標(biāo)和圖像坐標(biāo)就在兩個(gè)坐標(biāo)系下轉(zhuǎn)換,滿足經(jīng)典攝影測(cè)量共線方程的條件。(O - XYZ)與(O - xyz)的原點(diǎn)重合,為左右攝影中心 SL及 SR的中點(diǎn),Y軸與y軸指向飛行方向,Z軸與z軸指向地心相反方向,X軸與x軸分別垂直于YOZ與yOz平面,并構(gòu)成右手坐標(biāo)系。

        圖1 圖像拼接坐標(biāo)系Fig.1 Imagemosaiccoordinatesystem

        1.2 多視場(chǎng)遙感圖像拼接原理

        為方便討論,本文以兩幅圖像拼接為例進(jìn)行研究分析。對(duì)于同一地面點(diǎn)P(X,Y,Z),在兩幅圖像上的圖像坐標(biāo)由共線方程計(jì)算得到[8]:

        式中 (xL, yL)與(xR,yR)分別為地面點(diǎn)P在左右圖像像空間坐標(biāo)系下坐標(biāo); fL、 fR為左右圖像攝影主距。在以虛擬攝影中心為原點(diǎn)的像空間輔助坐標(biāo)系中,他們的空間坐標(biāo)分別為和分別為左右圖像的攝影中心物方坐標(biāo);分別為左右圖像外方位角元素的方向余弦參數(shù),構(gòu)成的旋轉(zhuǎn)矩陣為:此時(shí)地面點(diǎn)P和左右圖像的攝影中心在同一物方坐標(biāo)系下。

        所有左右圖像上的點(diǎn)由共線方程投影到像空間坐標(biāo)系中的過(guò) ( 0,0,- fM)點(diǎn)且平行于平面xOy的平面上,即得到了虛擬拼接圖像。與作為地面點(diǎn)P在左右圖像中的同名點(diǎn),投影到像空間坐標(biāo)系中的高程為 - fM的平面上,圖像坐標(biāo)的差異即是拼接誤差。

        2 多視場(chǎng)遙感圖像拼接誤差

        2.1 基于像方的圖像拼接誤差

        基于像方的拼接方案就是將待拼接不同圖像所在的像空間坐標(biāo)系坐標(biāo)投影到虛擬拼接圖像焦面上。首先建立虛擬拼接圖像坐標(biāo)系,一般原點(diǎn)為攝影基線中點(diǎn),然后找到待拼接圖像像空間坐標(biāo)系與虛擬拼接圖像坐標(biāo)系之間的轉(zhuǎn)換關(guān)系,最后將待拼接圖像所有像素坐標(biāo)投影到虛擬拼接圖像焦面上。

        式中Lλ與Rλ分別為左右圖像的攝影比例尺,消除攝影比例尺,整理得:

        將式(1)、(2)帶入式(5)、(6)得,

        將左右圖像的坐標(biāo)歸化到虛擬拼接圖像坐標(biāo),才可以進(jìn)行比較。由于左右圖像的像主點(diǎn)在虛擬拼接圖像中的坐標(biāo)分別為和,所以與在虛擬拼接圖像中滿足多視場(chǎng)遙感圖像像方拼接在x方向上誤差,在y方向上誤差

        從式(11)可以看出,像方拼接圖像的誤差源為地面點(diǎn)坐標(biāo)、攝影基線及虛擬拼接圖像主距,在投影到虛擬像平面時(shí),當(dāng)待拼接圖像主距與虛擬圖像主距相等的情況下,相機(jī)的姿態(tài)角及主距不產(chǎn)生拼接誤差??紤]到實(shí)際成像時(shí)一般垂直攝影,兩臺(tái)相機(jī)攝影中心沒(méi)有高差,可以認(rèn)為攝影基線B在Z方向上的分量此時(shí)基于像方的圖像拼接誤差為,

        2.2 基于物方的圖像拼接誤差

        生產(chǎn)單位一般采用基于物方的圖像拼接?;谖锓降钠唇臃桨甘紫葘⒋唇訄D像投影到物方的某一高程面上,一般為覆蓋地區(qū)的平均高程,即給定物方投影面的坐標(biāo) ZM,再由物方坐標(biāo)(X ,Y,ZM)根據(jù)待拼接圖像的嚴(yán)格成像模型計(jì)算其在圖像中的圖像坐標(biāo)與,重采樣原始圖像,即得到系統(tǒng)幾何校正圖像,最后根據(jù)系統(tǒng)幾何校正后圖像像素的地理坐標(biāo)進(jìn)行圖像拼接。

        考慮到圖像坐標(biāo)系與物方坐標(biāo)系之間的轉(zhuǎn)換關(guān)系,此時(shí)

        考慮到實(shí)際成像時(shí)一般垂直攝影,衛(wèi)星沒(méi)有俯仰、偏航和側(cè)擺。令為左右攝影中心連線的中點(diǎn),攝影基線為,為了計(jì)算方便,令即待拼接的兩臺(tái)相機(jī)攝影時(shí)攝影中心高程值相同。此時(shí)式(12)、(13)簡(jiǎn)化為

        基于物方的拼接只進(jìn)行了一次投影變換,多視場(chǎng)遙感圖像像方拼接在x方向上誤差在y方向上誤差

        從上式可以看出,基于物方拼接圖像的誤差源為攝影基線、拼接圖像主距及物方投影高程。當(dāng)待拼接相機(jī)垂直成像且無(wú)側(cè)擺時(shí),且待拼接圖像主距與虛擬圖像主距相等的情況下,基于物方的圖像拼接誤差大小與基于像方的拼接方法相同。

        3 仿真實(shí)驗(yàn)及分析

        根據(jù)上述理論推導(dǎo)及分析,采用仿真數(shù)據(jù)實(shí)驗(yàn),對(duì)圖像拼接誤差進(jìn)行仿真實(shí)驗(yàn)分析。設(shè)定兩臺(tái)框幅式遙感相機(jī),攝影基線垂直于飛行方向的分量為0.2m,即攝影基線X方向上的分量 BX= 0 .2m,攝影基線平行于飛行方向的分量為0.01m,即攝影基線Y方向上的分量 BY= 0 .01m,BX、BY是相機(jī)安裝前后不對(duì)齊造成的。仿真實(shí)驗(yàn)中模仿待拼接的兩臺(tái)相機(jī)的攝影中心連線垂直于飛行方向,安裝時(shí)有0.01m的誤差,物方空間坐標(biāo)系下左右攝影中心坐標(biāo)分別為,,地面點(diǎn)P物方坐標(biāo),基于像方拼接的虛擬拼接圖像主距Mf=2000mm,垂直攝影情況下航高 500H= km,即攝影比例尺為1∶250000,焦面CCD像元數(shù)為60000×60000個(gè),像元分辨率為2.5m,圖像重疊區(qū)為4000個(gè)像元,基于物方拼接的物方投影高程為50m,則投影面MZ =–499950m,相機(jī)安裝側(cè)擺角為0°,CCD尺寸為10μm×10μm。選擇平面坐標(biāo)分布均勻地面點(diǎn),在高程不同時(shí),同名點(diǎn)在基于像方拼接圖像上的坐標(biāo)差異由式(11)計(jì)算,計(jì)算結(jié)果見(jiàn)表1和表2。同名點(diǎn)在基于物方拼接圖像上的坐標(biāo)差異由式(16)計(jì)算,計(jì)算結(jié)果為dx=8.0×10–7m=0.08像元,dy=0.4×10–7m=0.004像元。

        表1 高程為100m時(shí)地面點(diǎn)位置對(duì)基于像方拼接誤差的影響Tab.1 Influences of ground point position on virtual image geometric mosaic error with elevation value 100m

        實(shí)驗(yàn)中選擇16個(gè)地面點(diǎn)(見(jiàn)表1和表2),在高程分別為100m和1000m時(shí)的32種情況下來(lái)分析基于像方拼接的理論誤差。由于基于像方拼接直接將待拼接圖像坐標(biāo)投影到虛擬拼接圖像焦面上,地面點(diǎn)高程的影響較小。隨著地面點(diǎn)在拼接區(qū)內(nèi)遠(yuǎn)離虛擬拼接圖像主點(diǎn),拼接誤差增大;x方向與y方向的拼接誤差變化規(guī)律一致。當(dāng)?shù)孛孀鴺?biāo)中X或Y偏離拼接物方坐標(biāo)系原點(diǎn)100m(即40個(gè)像元)時(shí),拼接誤差為0.2個(gè)像元左右,可以直接根據(jù)兩幅圖像成像幾何模型計(jì)算同一地面點(diǎn)的像點(diǎn)坐標(biāo)進(jìn)行拼接。當(dāng)?shù)孛孀鴺?biāo)X或Y偏離物方坐標(biāo)系原點(diǎn)1000m(即400個(gè)像元)時(shí),拼接誤差為3個(gè)像元左右,不能對(duì)待拼接圖像進(jìn)行直接拼接。當(dāng)?shù)孛孀鴺?biāo)X或Y偏離物方坐標(biāo)系原點(diǎn)10000m(即4000個(gè)像元)時(shí),拼接誤差已達(dá)到20像元左右。隨著地面坐標(biāo)的變化,其在待拼接圖像中對(duì)應(yīng)的同名點(diǎn)投影到虛擬拼接圖像焦面上的坐標(biāo)也變化,且主要與地面點(diǎn)平面位置變化有關(guān)。一般圖像間重疊區(qū)域?yàn)?00個(gè)像元,即百米級(jí),此時(shí)拼接誤差不超過(guò)1個(gè)像元,所以當(dāng)待拼接圖像成像幾何模型穩(wěn)定,經(jīng)過(guò)系統(tǒng)幾何校正消除系統(tǒng)誤差后,可以將同名點(diǎn)根據(jù)成像幾何模型投影到虛擬圖像中直接拼接。當(dāng)嚴(yán)格檢校出攝影基線B后,可以根據(jù)式(11)計(jì)算待拼接圖像中同名點(diǎn)在虛擬拼接圖像中的拼接誤差值,進(jìn)而對(duì)拼接區(qū)重疊像元數(shù)設(shè)計(jì)提供參考。

        表2 高程為1000m時(shí)地面點(diǎn)位置對(duì)基于像方拼接誤差的影響Tab.2 Influences of ground point position on virtual image geometric mosaic error with elevation value 1 000m

        通過(guò)仿真實(shí)驗(yàn)分析得出,當(dāng)垂直攝影且待拼接圖像主距等于虛擬拼接圖像主距時(shí),在重疊區(qū)為100個(gè)像元左右的情況下,如果成像幾何模型穩(wěn)定,基于像方和物方的拼接方法拼接誤差都不到1個(gè)像元,且拼接誤差具有系統(tǒng)性,通過(guò)拼接理論誤差推導(dǎo),可以計(jì)算基于像方及物方的拼接方法拼接理論誤差,經(jīng)過(guò)系統(tǒng)誤差校正,可以大幅降低,不影響目視判讀。

        4 結(jié)束語(yǔ)

        本文由嚴(yán)格共線方程幾何定位模型推導(dǎo)了不同視場(chǎng)圖像基于虛擬焦面的拼接誤差。拼接圖像可以通過(guò)將待拼接圖像投影到虛擬像平面或物方高程面獲得。將待拼接圖像投影到虛擬像平面的基于像方的拼接方法,需要將待拼接圖像投影后坐標(biāo),轉(zhuǎn)換到虛擬像平面所在的圖像坐標(biāo)系下。拼接誤差主要受地面點(diǎn)坐標(biāo)、攝影基線及虛擬拼接圖像主距影響,當(dāng)待拼接圖像主距與虛擬圖像主距相等的情況下,待拼接圖像的姿態(tài)角及主距不產(chǎn)生拼接誤差。將待拼接圖像投影到物方高程面的基于物方的拼接方法直接給定地面點(diǎn)高程,計(jì)算地面點(diǎn)在虛擬拼接圖像中的圖像坐標(biāo)進(jìn)行拼接,拼接誤差主要受攝影基線、待拼接圖像主距及物方投影高程面影響。

        基于物方的圖像拼接方法,待拼接圖像中的同名點(diǎn)投影到固定高程面上,在不考慮俯仰、翻滾及側(cè)擺的影響時(shí),拼接誤差是一個(gè)系統(tǒng)值,僅與攝影基線、虛擬拼接圖像主距和虛擬投影高程面有關(guān),這三個(gè)因素都可以進(jìn)行地面標(biāo)定。仿真結(jié)果顯示,基于物方的圖像拼接方法誤差較小,在模型穩(wěn)定的情況下,可以直接根據(jù)成像幾何模型計(jì)算待拼接圖像的像點(diǎn)坐標(biāo)在虛擬拼接圖像中的坐標(biāo)位置,進(jìn)行圖像拼接。下一步的研究應(yīng)重點(diǎn)關(guān)注相機(jī)幾何定位模型誤差及相機(jī)狀態(tài)變化時(shí)對(duì)拼接誤差的影響,以及對(duì)不同地區(qū)真實(shí)遙感圖像進(jìn)行實(shí)驗(yàn)分析對(duì)理論進(jìn)行驗(yàn)證。

        References)

        [1]杜延峰.一種新型大面陣CCD航測(cè)相機(jī)等效中心投影影像生成技術(shù)研究[D].西安:長(zhǎng)安大學(xué),2010.DU Yanfeng. Study on Equivalent Center Projection Image Generation Technology on a New Large Array CCD Aerial Camera[D]. Xi’an: Chang’an University, 2010. (in Chinese)

        [2]白照廣.高分一號(hào)衛(wèi)星的技術(shù)特點(diǎn)[J].中國(guó)航天,2013,8:5-9.BAI Zhaoguang. GF-1 Satellite Technique Characteristics[J]. Aerospace China, 2013, 8: 5-9. (in Chinese)

        [3]HUA Shungang, ZENG Lingyi, OU Zongying. Fast Algorithm for Cylindrical Panoramic Image Mosaic[J]. Journal of Data Acquisition&Processing, 2006, 21(4): 434-438.

        [4]ZHANG Jing, HU Zhiping, LIU Zhitai. Image Automatic Mosaics Based on Contour Phase Correlation[J]. Frontiers of Mechanical Engineering in China. 2007, 2(2): 228-234.

        [5]周清華,潘俊,李德仁.遙感圖像鑲嵌接縫線自動(dòng)生成方法綜述[J].國(guó)土資源與遙感,2013,25(2):1-7.ZHOU Qinghua, PAN Jun, LI Deren. Overview of Automatic Generation of Mosaicking Seamlines for Remote Sensing Images[J]. Remote Sensing For Land Resources, 2013, 25(2): 1-7. (in Chinese)

        [6]岳春宇,鄭永超,陶宇亮.星載激光測(cè)高儀輔助衛(wèi)星攝影測(cè)量淺析[J].航天返回與遙感,2013,34(4):71-76.YUE Chunyu, ZHENG Yongchao, TAO Yuliang. Study on Space-borne Laser Altimeter Supported Satellite Photogrammetry[J]. Spacecraft Recovery Remote Sensing, 2013, 34(4): 71-76. (in Chinese)

        [7]張過(guò).缺少控制點(diǎn)的高分辨率衛(wèi)星遙感影像幾何糾正[D].武漢:武漢大學(xué),2005.ZHANG Guo. Rectification for High Resolution Remote Sensing Image Under Lack of Ground Control Points[D]. Wuhan: Wuhan University, 2005. (in Chinese)

        [8]余俊鵬.高分辨率衛(wèi)星遙感影像的精確幾何定位[D].武漢:武漢大學(xué),2009.YU Junpeng. Accurate Geometric Positioning of High Resolution Satellite Remote Sensing Imagery[D]. Wuhan: Wuhan University, 2009. (in Chinese)

        [9]高穎慧,沈振康.航空搖擺圖像組的配準(zhǔn)拼接算法[J].紅外與激光工程,2009,38(1):170-174.GAO Yinghui, SHEN Zhenkang. Registration and Stitching Algorithm for Aviation Swaying Image Set[J]. Infrared and Laser Engineering, 2009, 38(1): 170-174. (in Chinese)

        [10]王慧,吳云東,張永生.面陣CCD數(shù)字航測(cè)相機(jī)影像幾何拼接誤差模型與分析[J].測(cè)繪學(xué)院學(xué)報(bào),2003,20(4):257-262.WANG Hui, WU Yundong, ZHANG Yongsheng. Modeling and Analyzing of Geometric Joint Error for CCD Matrix Images of Digital Aerial Camera[J]. Journal of Institute of Surveying and Mapping, 2003, 20(4): 257-262. (in Chinese)

        [11]王慧.面陣CCD航測(cè)相機(jī)成像模型與處理技術(shù)[D].鄭州:解放軍信息工程大學(xué),2006.WANG Hui. Imaging Model and Processing of Large Format CCD Aerial Camera[D]. Zhengzhou: PLA University of Information Engineering, 2006. (in Chinese)

        [12]徐彭梅,楊樺,伏瑞敏,等.CBERS-1衛(wèi)星CCD相機(jī)的光學(xué)拼接、配準(zhǔn)和定焦[J].航天返回與遙感,2001,22(3):12-15.XU Pengmei, YANG Hua, FU Ruimin, etal. Optical Butting, Registration and Focus-fixing of CCD Camera for CBERS-1 Satellite[J]. Spacecraft Recovery Remote Sensing, 2001, 22(3): 12-15. (in Chinese)

        [13]溫紅艷.遙感圖像拼接算法研究[D].武漢:華中科技大學(xué),2009.WEN Hongyan. Study on Remote Sensing Images Mosaic[D] Wuhan: Huazhong University of Science and Technology, 2009. (in Chinese)

        [14]阮鵬.多圖像拼接算法研究[D].長(zhǎng)沙:中南大學(xué),2009.RUANPeng.Studyon Images Mosaic[D].Changsha:Zhongnan University,2009.(inChinese)

        [15]嚴(yán)磊.基于特征匹配的全自動(dòng)圖像拼接算法研究[D].合肥:中國(guó)科學(xué)技術(shù)大學(xué),2009.YAN Lei. Automatic Image Mosaic Method Based on Feature Matching[D]. Hefei: University of Science and Technology of China, 2009. (in Chinese)

        猜你喜歡
        投影坐標(biāo)系高程
        解變分不等式的一種二次投影算法
        8848.86m珠峰新高程
        基于最大相關(guān)熵的簇稀疏仿射投影算法
        找投影
        找投影
        解密坐標(biāo)系中的平移變換
        坐標(biāo)系背后的故事
        基于重心坐標(biāo)系的平面幾何證明的探討
        GPS控制網(wǎng)的高程異常擬合與應(yīng)用
        SDCORS高程代替等級(jí)水準(zhǔn)測(cè)量的研究
        亚洲国产日韩精品一区二区三区| 91日韩东京热中文字幕| 国产色av一区二区三区| 青青草小视频在线播放| 亚洲av福利天堂一区二区三| 无码爆乳护士让我爽| 国产精品无码一区二区在线看| 亚洲精品无码久久久久| 亚洲av鲁丝一区二区三区| 九九九影院| 天天做天天爱天天爽综合网| 老妇肥熟凸凹丰满刺激| 国产爽爽视频在线| 亚洲无码视频一区:| 日本一区二区三区在线观看免费| 亚洲性感毛片在线视频| 中国老熟女露脸老女人| 国产精品亚洲色婷婷99久久精品| 亚洲日本一区二区一本一道| 免费黄色电影在线观看| 无码中文av有码中文av| 中文字幕人成乱码中文乱码| 日韩一级137片内射视频播放| 99久久精品免费看国产一区二区三区| 亚洲精品久久久久久久不卡四虎 | 蜜臀人妻精品一区二区免费| 人妻少妇中文字幕在线观看| 高潮内射双龙视频| 激情综合一区二区三区| 91精品久久久久含羞草| 中文字幕久久熟女人妻av免费| 日本熟妇中出高潮视频| 妺妺窝人体色www在线| 久久精品国产亚洲av蜜臀| 在线精品国内视频秒播| 漂亮的小少妇诱惑内射系列| 日韩一区二区av伦理| 包皮上有一点一点白色的| 中国丰满人妻videoshd| 国产欧美日韩a片免费软件| 在线观看国产精品91|