ZY-3衛(wèi)星異軌立體影像幾何定位精度分析

胡芬楊博唐新明高小明

（1 國(guó)家測(cè)繪地理信息局衛(wèi)星測(cè)繪應(yīng)用中心，北京 100048）

（2 衛(wèi)星測(cè)繪技術(shù)與應(yīng)用國(guó)家測(cè)繪地理信息局重點(diǎn)實(shí)驗(yàn)室，北京 100048）

（3 武漢大學(xué)計(jì)算機(jī)學(xué)院，武漢 430079）

（4 武漢大學(xué)地球空間信息技術(shù)協(xié)同創(chuàng)新中心，武漢430079）

ZY-3衛(wèi)星異軌立體影像幾何定位精度分析

胡芬1,2楊博3,4唐新明1,2高小明1,2

（1 國(guó)家測(cè)繪地理信息局衛(wèi)星測(cè)繪應(yīng)用中心，北京 100048）

（2 衛(wèi)星測(cè)繪技術(shù)與應(yīng)用國(guó)家測(cè)繪地理信息局重點(diǎn)實(shí)驗(yàn)室，北京 100048）

（3 武漢大學(xué)計(jì)算機(jī)學(xué)院，武漢 430079）

（4 武漢大學(xué)地球空間信息技術(shù)協(xié)同創(chuàng)新中心，武漢430079）

文章選取“資源三號(hào)”（ZY-3）衛(wèi)星在兩個(gè)軌道上分別側(cè)擺10°和–10.6°獲取的山東同一地區(qū)全色正視相機(jī)影像，構(gòu)成連續(xù)覆蓋 4景區(qū)域的異軌立體影像。以區(qū)域內(nèi)均勻分布的 35個(gè)高精度外業(yè)GPS點(diǎn)作為控制數(shù)據(jù)來(lái)源，進(jìn)行影像平差定向和幾何定位精度統(tǒng)計(jì)。實(shí)驗(yàn)結(jié)果表明，在基高比為0.36左右的非理想條件下，無(wú)控制點(diǎn)的平面誤差和高程誤差均值分別為 1.5m和 85m，中誤差分別為 2.3m和4.0m，高程精度明顯低于平面精度；基于少量控制點(diǎn)對(duì)系統(tǒng)誤差改正效果明顯，在單控制點(diǎn)條件下的平面誤差和高程誤差均值分別為2.4m和1.7m，中誤差分別為1.4m和4.0m，完全滿足1∶5萬(wàn)比例尺精度要求。文章的研究?jī)?nèi)容對(duì)于國(guó)產(chǎn)光學(xué)衛(wèi)星異軌立體影像獲取效能的發(fā)揮能起到一定的參考作用。

異軌立體影像 區(qū)域網(wǎng)平差 幾何定位精度分析 “資源三號(hào)”衛(wèi)星

0 引言

同軌觀測(cè)和異軌觀測(cè)是衛(wèi)星獲取立體影像的兩種典型方式[1-2]。早在1986年，法國(guó)發(fā)射的SPOT衛(wèi)星所搭載的高分辨率可見(jiàn)光傳感器HRV的地面分辨率達(dá)到10m，并能通過(guò)側(cè)視觀測(cè)在相鄰軌道間構(gòu)成異軌立體，其良好的基高比很適合立體測(cè)圖，對(duì)衛(wèi)星影像的測(cè)繪應(yīng)用產(chǎn)生了重大影響[3-4]。隨后的SPOT-5衛(wèi)星進(jìn)一步提高了立體成像能力，可以獲取同軌或異軌立體影像[5-7]。近年來(lái)，IKONOS，Geoeye、Pleiades-1A/1B、SPOT-6/7、Worldview-1/2/3等高分辨率商業(yè)遙感衛(wèi)星平臺(tái)均具備了大角度、快速姿態(tài)機(jī)動(dòng)能力，能靈活地實(shí)現(xiàn)同軌或異軌立體觀測(cè)[8]。國(guó)產(chǎn)光學(xué)衛(wèi)星立體影像獲取能力也逐步提高[9-10]，“資源三號(hào)”（ZY-3）、“天繪一號(hào)”等國(guó)產(chǎn)高分辨率光學(xué)測(cè)繪遙感衛(wèi)星沿用SPOT-5、IRS-P5、ALOS等衛(wèi)星的多線陣同軌立體攝影測(cè)量體制，即安裝不同指向的2臺(tái)或3臺(tái)線陣相機(jī)來(lái)完成立體觀測(cè)任務(wù)，具有立體觀測(cè)條件理想、數(shù)據(jù)獲取高效、利于大規(guī)模測(cè)繪產(chǎn)品生產(chǎn)等特點(diǎn)，盡管此類衛(wèi)星還具備一定程度的側(cè)擺能力，能夠獲取異軌立體影像，提升成像方式的靈活性，但是異軌立體影像的幾何定位精度和應(yīng)用潛力一直未被發(fā)掘。

為促進(jìn)衛(wèi)星數(shù)據(jù)獲取效能的發(fā)揮，本文選取ZY-3衛(wèi)星在兩個(gè)軌道上分別側(cè)擺10°和–10.6°獲取的山東同一地區(qū)全色正視相機(jī)影像，構(gòu)成連續(xù)覆蓋4景區(qū)域的異軌立體影像；以區(qū)域內(nèi)分布均勻的35個(gè)高精度外業(yè)GPS點(diǎn)作為控制數(shù)據(jù)來(lái)源，在控制點(diǎn)像點(diǎn)坐標(biāo)高精度量測(cè)和影像同名點(diǎn)密集匹配的基礎(chǔ)上，進(jìn)行影像平差定向和幾何定位精度統(tǒng)計(jì)。結(jié)果表明，基于少量控制點(diǎn)補(bǔ)償系統(tǒng)誤差之后，平面誤差和高程誤差均值有較大改善，分別為2.4m和1.7m，中誤差分別為1.4m和4.0m，滿足1∶5萬(wàn)比例尺精度要求，由此可見(jiàn)正視相機(jī)獲取的異軌立體影像能夠作為同軌立體影像的補(bǔ)充，提升數(shù)據(jù)來(lái)源的多樣性。

1 定位精度理論分析

影響衛(wèi)星立體影像幾何定位精度的因素很多[2]，概括起來(lái)主要包括以下幾個(gè)方面：相機(jī)分辨率、相機(jī)主點(diǎn)和畸變、CCD拼接精度、姿態(tài)角誤差、像點(diǎn)量測(cè)誤差和軌道誤差、時(shí)統(tǒng)誤差、攝影基高比等。在無(wú)控制點(diǎn)情況下，姿態(tài)測(cè)量精度和事后姿態(tài)處理精度是最重要的技術(shù)指標(biāo)，姿態(tài)測(cè)量精度的影響因素包括星敏感器和陀螺以及相應(yīng)的控制系統(tǒng)；在有控制點(diǎn)情況下，衛(wèi)星的姿態(tài)穩(wěn)定度是最重要的技術(shù)指標(biāo)，影響姿態(tài)穩(wěn)定度的則是衛(wèi)星上的活動(dòng)部件，以及相應(yīng)的姿態(tài)控制系統(tǒng)。通常，滿足高程精度相比平面精度更為困難。

在傳感器幾何關(guān)系經(jīng)過(guò)精確檢校、并且立體像對(duì)之間已轉(zhuǎn)換為核線幾何關(guān)系的理想情況下，高程誤差的理論計(jì)算公式為

式中 z為高程誤差；ε為像點(diǎn)匹配沿核線方向的誤差；B為攝影基線；Ht為軌道高度。可見(jiàn)，立體觀測(cè)由像對(duì)不同的視差產(chǎn)生，而基高比 B/Ht對(duì)高程測(cè)量精度有較大的影響。基高比越大，交會(huì)條件越好，高程量測(cè)的理論精度更高。傳統(tǒng)攝影測(cè)量在分辨率一定、量測(cè)精度受限的前提下，為了獲得較高的高程精度，基本上都采用大基高比（0.6～1）的方案，以實(shí)現(xiàn)地形三維信息的快速提取。小基高比情況下，像點(diǎn)量測(cè)誤差和匹配誤差對(duì)高程方向定位精度的影響將被放大，要達(dá)到同等精度的高程測(cè)量值，保證區(qū)域網(wǎng)平差精度的可靠性，對(duì)像點(diǎn)匹配算法及精度會(huì)有更高的要求[11]。

2 高精度幾何定位方法與流程

2.1 區(qū)域網(wǎng)平差定向的數(shù)學(xué)模型

區(qū)域網(wǎng)平差定向是實(shí)現(xiàn)影像高精度幾何定位的有效途徑[12-15]。本文通過(guò)構(gòu)建基于有理多項(xiàng)式模型（RFM）的影像區(qū)域網(wǎng)平差數(shù)學(xué)模型，利用影像連接點(diǎn)以及外業(yè)測(cè)量獲取的控制點(diǎn)，進(jìn)行區(qū)域網(wǎng)平差解算得到精化后的 RFM 參數(shù)和連接點(diǎn)地面坐標(biāo)，進(jìn)而利用檢查點(diǎn)統(tǒng)計(jì)無(wú)控制或有控制條件下的立體影像幾何定位精度。

RFM是將像方影像坐標(biāo)表示為以相應(yīng)物方地面三維空間坐標(biāo)為自變量的多項(xiàng)式的比值，其基本公式如式（2）所示，其中地面坐標(biāo)和影像坐標(biāo)均為正則化之后的坐標(biāo)，以保證計(jì)算的穩(wěn)定性。

這里， (φn,λn,Hn)為正則化的地面點(diǎn)大地坐標(biāo)，φn為緯度，λn為經(jīng)度，Hn為高程坐標(biāo)；(ln,sn)為正則化的影像像點(diǎn)坐標(biāo)，其中l(wèi)n代表行方向坐標(biāo)，而sn則代表列方向坐標(biāo)； p1i jk,p2 ijk,p3 ijk,p4ijk（i=1,2,3; j=1,2,3;k =1,2,3）代表 RFM 的有理多項(xiàng)式系數(shù)(rational polynomial coefficients, RPC)。(φn,λn,Hn)、(ln,sn)與非正則化的地面點(diǎn)大地坐標(biāo)(φ,λ,H)和像點(diǎn)坐標(biāo)(l,s)之間的關(guān)系為

式中 (l0,s0)為像點(diǎn)坐標(biāo)的平移參數(shù)；(φ0,λ0,H0)為地面點(diǎn)大地坐標(biāo)的平移參數(shù)；(ls, ss)為像點(diǎn)坐標(biāo)的縮放參數(shù)；（φs, λs,Hs） 為地面點(diǎn)大地坐標(biāo)的縮放參數(shù)；l代表像點(diǎn)行方向坐標(biāo)；s代表列方向坐標(biāo)。

為了便于后續(xù)闡述，本文將上述的RFM用函數(shù)形式表示為

基于RFM的影像區(qū)域網(wǎng)平差數(shù)學(xué)模型為

這里，Δl和Δs代表RFM系統(tǒng)誤差的像方補(bǔ)償模型，可根據(jù)具體情況采用平移變換模型、線性變換模型或仿射變換模型，分別如式（6）～（8）所示，模型參數(shù) a0, b0, a1, b1, a2, b2即為像方附加參數(shù)。

基于影像間自動(dòng)匹配獲取的連接點(diǎn)信息和人工半自動(dòng)量測(cè)獲取的控制點(diǎn)信息構(gòu)建區(qū)域網(wǎng)平差模型時(shí)，其原始平差觀測(cè)值包括連接點(diǎn)像點(diǎn)坐標(biāo)和控制點(diǎn)像點(diǎn)坐標(biāo)兩類。對(duì)于控制點(diǎn)像點(diǎn)而言，由于其對(duì)應(yīng)的物方點(diǎn)坐標(biāo)精確已知，因此，所構(gòu)建的誤差方程中未知參數(shù)僅包括該像點(diǎn)對(duì)應(yīng)影像RFM的像方附加參數(shù)，顯然，對(duì)于像方附加參數(shù)而言，式（5）為線性方程而無(wú)需進(jìn)行線性化處理。因此，可構(gòu)建誤差方程為

式中 vl，vs表示像點(diǎn)行向和列向坐標(biāo)觀測(cè)值的殘差。

然而，對(duì)于連接點(diǎn)像點(diǎn)而言，由于其對(duì)應(yīng)的物方點(diǎn)坐標(biāo)未知，因此，其所構(gòu)建的誤差方程式中未知參數(shù)除了包括影像RFM的像方附加參數(shù)外，還包括其對(duì)應(yīng)物方點(diǎn)坐標(biāo)(φ,λ,H)。對(duì)于物方點(diǎn)坐標(biāo)(φ,λ,H)而言，式（5）為非線性方程，需要對(duì)其賦予合適的初值(φ,λ,H)0并進(jìn)行線性化處理。因此，所構(gòu)建的誤差方程為：

2.2 高精度幾何定位流程

區(qū)域網(wǎng)平差在無(wú)控制的條件下又被稱為自由網(wǎng)平差，本文分別在無(wú)控制及不同控制的條件下對(duì)ZY-3衛(wèi)星異軌立體影像采用區(qū)域網(wǎng)平差的方式進(jìn)行高精度的幾何定位，其流程如圖1所示。

圖1 ZY-3衛(wèi)星異軌立體影像高精度幾何定位流程Fig.1 Workflow of high-accuracy geo-positioning of ZY-3 cross-track stereo-images

首先，通過(guò)利用SIFT特征匹配方法提取大量均勻分布的連接點(diǎn)，充分利用ZY-3衛(wèi)星軌道、姿態(tài)等參數(shù)，建立相鄰軌道影像間的相對(duì)幾何關(guān)系，并利用該相對(duì)幾何關(guān)系和基于縮放/旋轉(zhuǎn)幾何不變特征的高精度定位及多影像匹配算法，實(shí)現(xiàn)立體影像密集匹配，提高匹配的效率和可靠性，降低誤匹配率，對(duì)多種性能互補(bǔ)的匹配方法進(jìn)行融合，兼顧影像的局部信息和全局信息；通過(guò)自由網(wǎng)平差定向，對(duì)觀測(cè)像點(diǎn)的殘差進(jìn)行統(tǒng)計(jì)分析，進(jìn)而采用選權(quán)迭代法實(shí)現(xiàn)誤匹配連接點(diǎn)的自動(dòng)識(shí)別與剔除。

然后，進(jìn)行控制點(diǎn)半自動(dòng)量測(cè)，通過(guò)區(qū)域網(wǎng)平差定向和定向殘差分析，迭代計(jì)算得到區(qū)域網(wǎng)平差定向結(jié)果。

最后，在自由網(wǎng)平差或區(qū)域網(wǎng)平差的基礎(chǔ)上，利用檢查點(diǎn)對(duì)立體影像無(wú)控制或有控制條件下的幾何定位誤差進(jìn)行統(tǒng)計(jì)分析。

3 實(shí)驗(yàn)結(jié)果分析

3.1 實(shí)驗(yàn)區(qū)及數(shù)據(jù)選擇

ZY-3衛(wèi)星裝載的三線陣全色相機(jī)由具有22°夾角的前視、正視和后視3臺(tái)線陣CCD相機(jī)構(gòu)成，在軌狀態(tài)下從前、正和后三個(gè)方向獲取同一地物 100%地面重疊率的立體影像，前后視相機(jī)攝影基高比為0.89，具備了較為理想的同軌立體觀測(cè)條件[2]，為基于光學(xué)立體遙感影像的應(yīng)用服務(wù)提供了保障。盡管衛(wèi)星具備側(cè)擺±32°的能力，但出于對(duì)衛(wèi)星壽命的保障及降低損耗的考慮，現(xiàn)階段 ZY-3衛(wèi)星的工作機(jī)制對(duì)整星側(cè)擺的幅度有所限制，一般不超過(guò)15°，并且側(cè)擺大于5°時(shí)前后視全色相機(jī)不開(kāi)機(jī)成像。為此，本文暫且針對(duì)正視全色相機(jī)獲取的異軌立體影像展開(kāi)定位精度分析，并從實(shí)際情況出發(fā)，盡可能地選擇側(cè)擺幅度較大的異軌立體影像為研究對(duì)象。通過(guò)定位精度分析，得出對(duì)后續(xù)研究有參考價(jià)值的結(jié)論。

本文所選的實(shí)驗(yàn)區(qū)域位于山東省境內(nèi)，地形以丘陵為主。2014年10月24日的15502軌有覆蓋該地區(qū)的4景正視相機(jī)（NAD）影像，衛(wèi)星側(cè)擺約10°；為了構(gòu)成滿足要求的異軌立體影像，綜合考慮實(shí)際情況，2015年4月5日的17979軌，衛(wèi)星通過(guò)側(cè)擺–10.6°，再次獲取了覆蓋該地區(qū)的4景正視相機(jī)影像。以這8幅影像的傳感器校正（SC）級(jí)產(chǎn)品作為數(shù)據(jù)來(lái)源，構(gòu)成4個(gè)立體像對(duì)。影像基本信息如表1所示。SC級(jí)產(chǎn)品是對(duì)原始衛(wèi)星下傳影像經(jīng)過(guò)輻射校正和傳感器幾何校正（包括CCD條帶拼接、積分時(shí)間規(guī)劃、內(nèi)方位元素規(guī)劃等），但沒(méi)有經(jīng)過(guò)幾何處理的一級(jí)影像產(chǎn)品，產(chǎn)品帶有高精度的 RFM 參數(shù)，是進(jìn)行實(shí)際應(yīng)用的一種最基礎(chǔ)的分發(fā)產(chǎn)品[16]。影像質(zhì)量良好，云量小于 5%，異軌立體重疊度優(yōu)于 80%，基高比約為0.36。

表1 影像信息列表Tab.1 Information list of experimental images

3.2 實(shí)驗(yàn)過(guò)程與結(jié)果分析

在平差解算前采用密集的方式對(duì)立體影像進(jìn)行連接點(diǎn)匹配，得到匹配點(diǎn)的分布如圖2所示。由于第3個(gè)立體像對(duì)上存在少量的朵云，局部區(qū)域無(wú)法提取影像連接點(diǎn)。此外，在實(shí)驗(yàn)區(qū)域內(nèi)通過(guò)GPS外業(yè)實(shí)測(cè)手段按一定格網(wǎng)間距測(cè)量了35個(gè)均勻分布的高精度控制點(diǎn)（如圖2）。由于基高比相對(duì)較小，控制點(diǎn)量測(cè)精度對(duì)高程精度驗(yàn)證的影響較大，針對(duì)此問(wèn)題，在進(jìn)行控制點(diǎn)外業(yè)測(cè)量時(shí)直接以ZY-3影像為參考數(shù)據(jù)，控制點(diǎn)選取影像上易于辨識(shí)的、位于地形平坦區(qū)域的明顯地物特征點(diǎn)，如道路交叉口等，像點(diǎn)量測(cè)精度能夠達(dá)到子像素級(jí)。

圖2 實(shí)驗(yàn)區(qū)范圍及控制點(diǎn)、連接點(diǎn)分布Fig.2 Test site area and distribution of GCPs and tie points

平差定向的結(jié)果與所構(gòu)建的誤差方程緊密有關(guān)，而誤差方程構(gòu)建又主要取決于參與平差的影像數(shù)據(jù)，對(duì)此，本文實(shí)驗(yàn)采取兩種平差定向方式：1）對(duì)4個(gè)立體像對(duì)逐一平差定向，然后分別統(tǒng)計(jì)每個(gè)立體像對(duì)的幾何定位精度；2）對(duì)所有影像進(jìn)行整體區(qū)域網(wǎng)平差，然后從總體上統(tǒng)計(jì)幾何定位精度。相對(duì)于前者，后者能充分利用立體像對(duì)之間的空間幾何約束關(guān)系。

對(duì)單個(gè)立體像對(duì)逐一定向，得到的幾何定位精度統(tǒng)計(jì)結(jié)果如表2所示。

表2 單個(gè)立體像對(duì)平差幾何定位精度Tab.2 Geo-positioning accuracy of each single stereo-image adjustment

對(duì)8幅影像進(jìn)行整體區(qū)域網(wǎng)平差，得到的幾何定位精度統(tǒng)計(jì)結(jié)果如表3和圖3所示。

表3 整體區(qū)域網(wǎng)平差幾何定位精度Tab.3 Geo-positioning accuracy of overall block adjustment

圖3 整體區(qū)域網(wǎng)平差幾何定位誤差分布圖Fig.3 Adjustment residuals maps of overall block adjustment

從實(shí)驗(yàn)結(jié)果可以看出，在無(wú)控制點(diǎn)的情況下，ZY-3衛(wèi)星幾何定位的平面誤差均值為 1.5m，高程誤差均值為85m，高程方向上存在明顯的系統(tǒng)誤差。當(dāng)利用一個(gè)控制點(diǎn)參與平差時(shí)，高程方向系統(tǒng)誤差得到有效補(bǔ)償，高程精度得到明顯改善，誤差均值在1.7m左右，而平面精度則未有明顯改變。隨著控制點(diǎn)個(gè)數(shù)的增加，平面和高程的定位精度均有所提升。

無(wú)控制點(diǎn)的高程精度明顯低于平面精度，實(shí)驗(yàn)數(shù)據(jù)的基高比較小是一個(gè)主要原因。經(jīng)過(guò)整體平差后區(qū)域網(wǎng)內(nèi)部精度較高，僅利用單個(gè)控制點(diǎn)即可滿足要求。無(wú)論是單個(gè)立體像對(duì)平差還是整體區(qū)域網(wǎng)平差，基于單個(gè)控制點(diǎn)對(duì)系統(tǒng)誤差改正效果明顯。無(wú)控制點(diǎn)的平面中誤差和高程中誤差分別為2.3m和4.0m，在單控制點(diǎn)條件下，中誤差分別為1.4m和4.0m。增加控制點(diǎn)對(duì)中誤差改善不明顯，表明影像產(chǎn)品本身的內(nèi)部幾何畸變較小，達(dá)到了較高的精度水平。這說(shuō)明，當(dāng)影像內(nèi)部幾何畸變較小時(shí)，即使異軌立體影像的基高比不在理想范圍內(nèi)，基于少量控制點(diǎn)也可能實(shí)現(xiàn)較高的幾何定位精度。

單個(gè)立體像對(duì)平差和整體區(qū)域網(wǎng)模型平差后的幾何定位精度水平總體相當(dāng)，局部也無(wú)明顯差異。這表明利用 ZY-3衛(wèi)星影像產(chǎn)品進(jìn)行區(qū)域網(wǎng)平差時(shí)，幾何誤差在區(qū)域網(wǎng)內(nèi)部未發(fā)生明顯的傳遞與累積，區(qū)域網(wǎng)內(nèi)部的幾何精度具有高度的剛性結(jié)構(gòu)特征，不存在扭曲變形。

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

為促進(jìn)國(guó)產(chǎn)光學(xué)衛(wèi)星異軌立體影像獲取及應(yīng)用效能的發(fā)揮，提升數(shù)據(jù)來(lái)源的多樣性，本文對(duì) ZY-3衛(wèi)星異軌立體影像幾何定位精度展開(kāi)了初步分析。在基于少量控制點(diǎn)進(jìn)行系統(tǒng)誤差改正之后，實(shí)驗(yàn)數(shù)據(jù)的幾何定位平面誤差和高程誤差均值分別為2.4m和1.7m，中誤差分別為1.4m和4.0m，滿足1 ∶5萬(wàn)比例尺精度要求，表明 ZY-3衛(wèi)星具備高精度異軌立體觀測(cè)能力。從理論上講，假如增大立體成像的基高比，影像的定位精度有望更高。

由于時(shí)間倉(cāng)促且受實(shí)驗(yàn)條件等限制，研究仍有待深入，今后將主要從以下幾個(gè)方面進(jìn)行改進(jìn)：1）注重?cái)?shù)據(jù)選擇的典型性和廣泛性，例如用多個(gè)不同地形特點(diǎn)的實(shí)驗(yàn)區(qū)影像進(jìn)行更充分的驗(yàn)證，側(cè)擺角及基高比設(shè)置盡量靈活多變；2）研究多種立體像對(duì)構(gòu)建方式下的幾何定位精度，例如考慮利用前視和后視相機(jī)影像構(gòu)建異軌立體像對(duì)；3）探索基于多種控制數(shù)據(jù)來(lái)源的幾何定位精度改善方法，如考慮引入激光測(cè)高數(shù)據(jù)作為廣義高程控制點(diǎn)等；4）利用異軌立體影像生成DSM和DOM，并對(duì)產(chǎn)品精度進(jìn)行評(píng)價(jià)。

References)

[1]張永生, 鞏丹超, 劉軍, 等. 高分辨率遙感衛(wèi)星應(yīng)用: 成像模型、處理算法及應(yīng)用技術(shù)[M]. 北京: 科學(xué)出版社, 2004. ZHANG Yongsheng, GONG Danchao, LIU Jun, et al. The Application of High Resolution Satellite——the Imaging Model, Processing Algorithm & Application Technology[M]. Beijing: Science Press, 2004. (in Chinese)

[2]唐新明, 謝俊鋒, 張過(guò). 測(cè)繪衛(wèi)星技術(shù)總體發(fā)展和現(xiàn)狀[J]. 航天返回與遙感, 2012, 33(3): 17-24. TANG Xinming, XIE Junfeng, ZHANG Guo. Development and Status of Mapping Satellite Technology[J]. Spacecraft Recovery & Remote Sensing, 2012, 33(3): 17-24. (in Chinese)

[3]AI-ROUSAN N, CHENG P, PETRIEG, et al. Automated DEM Extraction and Orthoimage Generation from SPOT Level 1B Imagery [J]. Photogrammetry Engineering & Remote Sensing, 1997, 63(8): 965-974.

[4]燕琴, 張祖勛, 張劍清. 異軌遙感CCD 影像外方位元素的解求[J]. 武漢大學(xué)學(xué)報(bào)·信息科學(xué)版, 2001, 26(3): 270-274. YAN Qin, ZHANG Zuxun, ZHANG Jianqing. Orientation of Remote Sensing Images Taken by CCD from Different Orbits[J]. Geomatics and Information Science of Wuhan University, 2001, 26(3): 270-274. (in Chinese)

[5]張永軍, 張勇. SPOT-5 HRS 立體影像無(wú)(稀少)控制絕對(duì)定位技術(shù)研究[J]. 武漢大學(xué)學(xué)報(bào)·信息科學(xué)版, 2006, 31(11): 941-944. ZHANG Yongjun, ZHANG Yong. Direct Georeferencing of SPOT-5 HRS Imagery Without (or with a Few) Ground Control Point[J]. Geomatics and Information Science of Wuhan University, 2006, 31(11): 941-944. (in Chinese)

[6]徐岳仁, 張軍龍, 陳長(zhǎng)云. DGPS和SPOT-5異軌立體像對(duì)在岷江源構(gòu)造地貌研究中的應(yīng)用[J]. 北京大學(xué)學(xué)報(bào)(自然科學(xué)版), 2012, 48(4): 574-582. XU Yueren, ZHANG Junlong, CHEN Changyun. Applications of DGPS Method and SPOT-5 Across-track Stereo Images on Tectonic Geomorphology Research in Origin of Minjiang River[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2012, 48(4): 574-582. (in Chinese)

[7]賈秀鵬, 焦偉利, 李丹. 基于SPOT-5異軌立體像對(duì)提取DEM試驗(yàn)與精度評(píng)估[J]. 測(cè)繪信息與工程, 2006, 31(2): 32-34. JIA Xiupeng, JIAO Weili, LI Dan. DEM Extraction and Accuracy Evaluation based on SPOT-5 Across-track Stereo Images[J]. Journal of Geomatics, 2006, 31(2):32-34. (in Chinese)

[8]朱紅, 劉維佳, 張愛(ài)兵. 光學(xué)遙感立體測(cè)繪技術(shù)綜述及發(fā)展趨勢(shì)[J]. 現(xiàn)代雷達(dá), 2014, 36(6): 6-12. ZHU Hong, LIU Weiwei, ZHANG Aibing. Overview and Development Trend of Optical remote Sensing Stereo Mapping Techniques [J]. Modern Radar, 2014, 36(6): 6-12.

[9]冉瓊, 遲耀斌, 王智勇, 等. 基于“北京一號(hào)”小衛(wèi)星異軌立體像對(duì)提取 DEM[C]. 第十二屆中國(guó)體視學(xué)與圖像分析學(xué)術(shù)年會(huì)論文集, 2008: 425-428. RAN Qiong, CHI Yaobin, WANG Zhiyong, et al. DEM Extraction Based on Across-track Beijing-1 Small Satellite Stereo Images[C]. The Proceedings of the 12th Annual Conference of the Chinese Academy of Stereo Vision and Image Analysis, 2008: 425-428. (in Chinese)

[10]胡莘, 曹喜濱. 三線陣立體測(cè)繪衛(wèi)星的測(cè)繪精度分析[J]. 哈爾濱工業(yè)大學(xué)學(xué)報(bào), 2008, 40(5): 695-699. HU Xin, CAO Xibin. Journal of Harbin Institute of Technology, 2008, 40(5): 695-699. (in Chinese)

[11]MORGAN G L K, LIU Jianguo, YAN Hongshi. Precise Subpixel Disparity Measurement from Very Narrow Baseline Stereo [J]. IEEE Transactions on Geoscience and Remote Sensing, 2010, 48(9): 3424-3433.

[12]DIAL G, GRODECKI J. Adjustment with Rational Polynomial Camera Models[C]. ACSM-ASPRS 2002 Annual Conference Procedings, 2002.

[13]元建勝. 稀少控制點(diǎn)異軌立體衛(wèi)星影像區(qū)域網(wǎng)平差應(yīng)用研究[J]. 海洋測(cè)繪, 2012, 32(1): 45-48. YUAN Jiansheng. Research on Block Adjustment to Space Borne Remote Sensing Stereo Images Obtained from Different Orbits with Limited GCPs[J]. Hydrographic Surveying and Charting, 2012, 32(1): 45-48. (in Chinese)

[14]張永軍, 張劍清. 異軌遙感立體像對(duì)外方位元素的求解算法[J]. 武漢大學(xué)學(xué)報(bào)（信息科學(xué)版）, 2003, 28(5): 521-524. ZHANG Yongjun, ZHANG Jianqing. Geomatics and Information Science of Wuhan University[J]. 2003, 28(5): 521-524. (in Chinese)

[15]呂爭(zhēng), 傅俏燕, 王小燕. 資源三號(hào)衛(wèi)星正視影像區(qū)域網(wǎng)平差[J]. 航天返回與遙感, 2014, 35(1): 72-80. LV Zheng, FU Qiaoyan, WANG Xiaoyan. Research on Block Adjustment of ZY-3 Satellite NAD Images[J]. Spacecraft Recovery & Remote Sensing, 2014, 35(1): 72-80. (in Chinese)

[16]劉斌, 孫喜亮, 邸凱昌, 等. 資源三號(hào)衛(wèi)星傳感器校正產(chǎn)品定位精度驗(yàn)證與分析[J]. 國(guó)土資源遙感, 2012(4): 36-40. LIU Bin, SUN Xiliang, DI Kaichang, et al. Accuracy Analysis and Validation of ZY-3’s Sensor Corrected Products[J]. Remote Sensing for Land & Resources, 2012(4): 36-40. (in Chinese)

[17]唐新明. 資源三號(hào)測(cè)繪衛(wèi)星三線陣成像幾何模型構(gòu)建與精度初步驗(yàn)證[J]. 測(cè)繪學(xué)報(bào), 2012, 41(2): 191-198. TANG Xinming. Triple Linear-array Imaging Geometry Model of ZY-3 Mapping Satellite and Its Validation[J]. Acta Geodaeticaet Cartographica Sinica, 2012, 41(2): 191-198.(in Chinese)

[18]李松明, 李巖, 李勁東. “天繪一號(hào)”傳輸型攝影測(cè)量與遙感衛(wèi)星[J]. 遙感學(xué)報(bào), 2012, 16(增刊): 10-16. LI Songming, LI Yan, LI Jindong. Mapping Satellite-1 Transmission Type Photogrammetric and Remote Sensing Satellite[J]. Journal of Remote Sensing, 2012, 16(S1): 10–16. (in Chinese)

Geo-positioning Accuracy Analysis of ZY-3 Cross-track Stereo-images

HU Fen1,2YANG Bo3,4TANG Xinming1,2GAO Xiaoming1,2

（1 Satellite Surveying and Mapping Application Center, National Administration of Surveying, Mapping and Geoinformation, Beijing 100048, China）
（2 Key Laboratory of Satellite Mapping Technology and Application, NASG, Beijing 100048, China）
（3 Computer School of Wuhan University, Wuhan 430079, China）
（4 Collaborative Innovation Center for Geospatial Technology, Wuhan 430079, China）

In this paper, four continuous scenes of nadir panchromatic images covering a certain district of Shandong Province in China are captured twice from two orbits and then are chosen to form four pairs of cross-track stereo-images, with 10 and –10.6 degrees side-looking angles separately. Thirty-five accurate and evenly-distributed GPS points collected in field are used as GCPs, the image points of which are measured accurately on the images, and the geo-positioning accuracy of these cross-track stereo-images is analyzed. It is demonstrated by the experimental results that when the B/Htratio is 0.36 which is not ideal, the horizontal and vertical systematic error without GCPs are 1.5m and 85m, and the RMS error are 2.3m and 4m, separately, in which the horizontal accuracy is much better than that in the vertical direction. With only a few GCPs, the systematic error can be obviously compensated, that is, with only one GCP the mean error values are reduced to 2.4m and 1.7m in horizontal and vertical directions, with the RMS error at the 1.4m and 4m level, respectively.Overall, the 1∶50000 scale mapping requirement can be thoroughly reached. This research to a certain degree gives reference to the effectiveness exertion of cross-track stereo-images acquisition of domestic optical satellites.

cross-track stereo-images; block adjustment; geo-positioning accuracy analysis ; ZY-3 satellite

P237

A

1009-8518(2016)01-0071-09

10.3969/j.issn.1009-8518.2016.01.009

胡芬，女，1984年生，博士，2010年獲武漢大學(xué)攝影測(cè)量與遙感專業(yè)博士學(xué)位，副研究員。研究方向?yàn)樾滦蛡鞲衅鞯臄?shù)字?jǐn)z影測(cè)量、高分辨率衛(wèi)星影像幾何處理算法等。E-mail：huf@sasmac.cn。

（編輯：夏淑密）

2015-11-06

測(cè)繪地理信息公益性行業(yè)科研專項(xiàng)（201412007）；2015年基礎(chǔ)測(cè)繪科技項(xiàng)目（A1501-3）