潘 俊,胡 芬,王 密,金淑英,李國元

1.武漢大學(xué)測繪遙感信息工程國家重點實驗室,湖北 武漢 430079;2.國家測繪地理信息局衛(wèi)星測繪應(yīng)用中心,北京 101300

一種非共線TDI CCD成像數(shù)據(jù)內(nèi)視場拼接方法

潘 俊1,胡 芬2,王 密1,金淑英1,李國元2

1.武漢大學(xué)測繪遙感信息工程國家重點實驗室,湖北 武漢 430079;2.國家測繪地理信息局衛(wèi)星測繪應(yīng)用中心,北京 101300

星載非共線TDI CCD成像數(shù)據(jù)的高質(zhì)量內(nèi)視場拼接是保證影像后續(xù)處理和應(yīng)用的基礎(chǔ)。本文提出一種非共線TDI CCD成像數(shù)據(jù)內(nèi)視場拼接方法,該方法基于物方投影面和相機的傳感器幾何模型,建立了從拼接影像像點到原始影像像點的坐標(biāo)轉(zhuǎn)換關(guān)系,進而采用間接法影像糾正的方式對原始影像進行重采樣生成拼接影像。本文在描述非共線TDI CCD嚴(yán)格成像幾何模型之后,詳細(xì)敘述了方法的基本原理、潛在的誤差來源以及方法的實現(xiàn)流程,并通過ZY-1 02C衛(wèi)星的高分辨率(HR)相機成像數(shù)據(jù)進行了試驗驗證,對拼接精度的分析與評價證明了方法的可行性。

非共線TDI CCD;內(nèi)視場拼接;傳感器幾何模型;物方投影面;精度評價

1 引 言

時間延時積分電荷耦合元件(TDI CCD)能在不降低成像空間分辨率的同時,有效地提高空間相機的成像靈敏度和信噪比,因而已逐漸取代常規(guī)CCD器件成為高分辨率光學(xué)推掃式相機的主流傳感器[1]。目前,IKONOS、Quick Bird、World View、Pleiades等高分辨率商業(yè)衛(wèi)星以及我國資源一號02B(CBERS-02B)、資源一號02C星(ZY-1 02C)、資源三號(ZY-3)等衛(wèi)星上的光學(xué)相機均采用TDI CCD器件作為成像傳感器[2-10]。由于單片TDI CCD像元個數(shù)有限,在空間相機的設(shè)計上,為了保證獲取較高空間分辨率和一定幅寬的影像,需要將多片TDI CCD連接排列以獲得較大的成像視場。然而受器件等物理因素的限制,多片TDI CCD很難直接在焦平面上按照共線方式進行物理排列,因而很多傳感器采用了非共線的設(shè)計方式進行焦平面視場拼接,即將多片TDI CCD在焦平面上裝配成上下交錯的雙列[11],例如,ZY-1 02C所搭載的高分辨率(HR)相機采用了典型的三片非共線設(shè)計,其焦平面排列方式如圖1所示。采用這種成像設(shè)計的TDI CCD相機獲取的原始影像,受傳感器焦平面擺放位置及其誤差、相機光學(xué)系統(tǒng)畸變、偏流角控制偏差、衛(wèi)星姿態(tài)變化、地面起伏、行積分時間變化等因素的影響,難以直接形成一幅覆蓋相機全視場的連續(xù)無縫影像。為了滿足后續(xù)影像產(chǎn)品的生產(chǎn)要求,必須對原始數(shù)據(jù)進行高精度內(nèi)視場拼接,而多片TDI CCD的非共線設(shè)計使得分片影像間的錯位和偏移特性較為顯著和復(fù)雜,大大增加了高精度拼接處理的難度[12-14],因此針對該類相機成像數(shù)據(jù)進行內(nèi)視場拼接的理論與方法研究,具有重要的理論意義和應(yīng)用價值。

圖1 三片非共線TDI CCD焦平面排列示意圖Fig.1 Focal plane assembly of three non-collinear TDI CCD chips

出于對衛(wèi)星設(shè)計參數(shù)以及衛(wèi)星應(yīng)用技術(shù)的保密,商業(yè)高分辨率衛(wèi)星影像的前期視場拼接工作通常由數(shù)據(jù)供應(yīng)方完成,可供參考的文獻(xiàn)較少。目前,國內(nèi)學(xué)者針對以CBERS-02B為代表的非共線TDI CCD成像數(shù)據(jù)的內(nèi)視場拼接問題開展了一些研究工作。從理論基礎(chǔ)上看,現(xiàn)有的方法可以分為像方拼接和物方拼接兩大類。像方拼接本質(zhì)上是基于像面變換的片間配準(zhǔn)過程,如傳統(tǒng)的基于片間平移的方法[14-15]、基于分段仿射變換的方法[12,16]和基于行積分時間歸一化的方法[12]。像方拼接方法具有模型簡單、計算量小的優(yōu)點,通常能夠保證拼接后的影像目視無縫,但缺點在于拼接后的影像失去了攝影測量的投影幾何特性,無法對拼接影像進行高精度幾何建模,這在較大程度上制約了影像的后續(xù)應(yīng)用。相比而言,物方拼接具有嚴(yán)密的理論基礎(chǔ),其本質(zhì)上是從相機的成像幾何模型出發(fā),利用相鄰TDI CCD的物方攝影幾何約束條件建立嚴(yán)格的相機內(nèi)視場拼接模型,文獻(xiàn)[17-18]發(fā)展了基于虛擬線陣的方法[12],將其應(yīng)用于ALOS/PRISM掃描儀、ZY-3全色多光譜相機等共線TDI CCD相機成像數(shù)據(jù)的內(nèi)視場拼接處理中。文獻(xiàn)[12,19]提出了基于投影基準(zhǔn)面的方法,并基于模擬的非共線TDI CCD數(shù)據(jù)對拼接模型的正確性進行了驗證,但該方法采用了片間相對定向的策略,幾何模型并不嚴(yán)謹(jǐn)。本文在此基礎(chǔ)上,針對非共線TDI CCD成像數(shù)據(jù)的特點,提出了一種非共線TDI CCD成像數(shù)據(jù)內(nèi)視場拼接方法。該方法直接建立在原始影像的傳感器幾何模型之上,并考慮了各種潛在的誤差來源,對拼接處理的流程也進行了優(yōu)化,無需進行片間相對定向處理,適合于業(yè)務(wù)化的運行系統(tǒng)。

2 非共線TDI CCD的嚴(yán)格成像幾何模型

如圖2所示,以ZY-1 02C所搭載的HR相機為例,當(dāng)相機沿軌推掃成像時,在焦平面上按非共線排列的三片TDI CCD共有一套軌道姿態(tài)、相機安裝角、焦距和主點參數(shù),基于傳感器線陣與地面掃描投影線的映射關(guān)系,相機獲得的地面覆蓋范圍由品字形分布的3條掃描帶構(gòu)成。根據(jù)HR相機焦平面設(shè)計特點[7],在垂軌方向相鄰掃描帶之間有大約30個像元的重疊,在沿軌方向上兩側(cè)掃描帶相對于中間掃描帶有大約2600行的成像時延。三片TDI CCD影像按成像時間或行計數(shù)對齊后即構(gòu)成HR相機的原始影像。

在推掃成像時,每片TDI CCD都遵循線中心投影透視幾何原理,所有分片TDI CCD影像的嚴(yán)格成像幾何模型共同構(gòu)成原始影像的傳感器幾何模型。每片TDI CCD影像的嚴(yán)格成像幾何模型描述如下

圖2 HR相機成像模式和原始影像坐標(biāo)系示意圖Fig.2 Illustration of the imaging mode and the original image scene of HR camera

3 一種非共線TDI CCD成像數(shù)據(jù)內(nèi)視場拼接方法

3.1 基本原理

根據(jù)相機的嚴(yán)格成像幾何模型,本文提出一種非共線TDI CCD成像數(shù)據(jù)內(nèi)視場拼接方法。如圖3所示,假設(shè)物方投影面是定義在成像覆蓋區(qū)域的平均高程面,沿衛(wèi)星推掃成像的目標(biāo)軌跡方向,對物方投影面上的成像覆蓋區(qū)域劃分等間隔格網(wǎng),格網(wǎng)大小與原始影像的地面分辨率相一致。通過二維平移、旋轉(zhuǎn)和縮放變換,將拼接影像像點與所定義的物方投影面格網(wǎng)點一一對應(yīng),然后基于原始影像的傳感器幾何模型,對物方投影面格網(wǎng)點進行反投影計算,得到其對應(yīng)的原始影像像點坐標(biāo),由此便建立了從拼接影像像點到原始影像像點的坐標(biāo)轉(zhuǎn)換關(guān)系,進而采用間接法影像糾正的方式對原始影像進行重采樣生成拼接影像。根據(jù)拼接影像像點在原始影像上的位置以及原始影像的傳感器幾何模型,得到拼接影像像點的等效攝影光線,由此構(gòu)建拼接影像的等效攝影幾何模型。該方法生成拼接影像的過程等效于沿衛(wèi)星推掃成像的目標(biāo)軌跡方向?qū)υ加跋襁M行了一次準(zhǔn)系統(tǒng)幾何校正。

圖3 拼接方法原理示意圖Fig.3 Principle of the stitching method

事實上,物方拼接就是利用拼接影像對地覆蓋區(qū)域的連續(xù)性和原始影像的成像幾何模型,建立拼接影像像點與原始影像像點的嚴(yán)格坐標(biāo)映射關(guān)系。基于上述方法的原理,不難發(fā)現(xiàn),良好的片間攝影幾何約束條件是實現(xiàn)非共線TDI CCD成像數(shù)據(jù)嚴(yán)格無縫內(nèi)視場拼接的幾何基礎(chǔ),即要求片間幾何定位精度具有一致性。如圖2所示,具體表現(xiàn)為,對于落在相鄰兩片TDI CCD成像地面覆蓋重疊區(qū)內(nèi)的任一物點P,基于原始影像的傳感器幾何模型分別反投影計算其在相鄰兩片影像上的像點,得到的一對像點p1和p2應(yīng)滿足共軛關(guān)系,即是一對同名像點。然而,在實際情況下,非共線TDI CCD在成像時存在的各類觀測誤差,均會不同程度地影響片間幾何定位精度的一致性,應(yīng)該在拼接處理之前對相機幾何畸變及其安裝角誤差、衛(wèi)星平臺穩(wěn)定性等參數(shù)進行檢校,對各類系統(tǒng)誤差進行補償[7,12]。

3.2 地物高差對拼接精度的影響

除相機幾何畸變及其安裝角誤差、衛(wèi)星平臺穩(wěn)定性等參數(shù)外,物點與平均高程面之間的高差對非共線TDI CCD成像數(shù)據(jù)沿軌方向上的拼接精度也有一定影響[12]。如圖4所示,對于片間重疊區(qū)域內(nèi)的任一物點P,假設(shè)其高程值與平均高程值的偏差為ΔH;φ1和φ2分別為相機兩列TDI CCD偏場角的設(shè)計值;Δq表示由ΔH引起的影像拼接誤差在物方的投影距離,則有

由式(2)可知,地物高差對物方拼接精度的影響程度與相鄰兩片TDI CCD的偏場角正切的差值成正比。對于HR相機,φ2≈1.15°,φ1≈0.7°,地面分辨率為2.36 m,根據(jù)式(2)計算,當(dāng)?shù)匚稂c相對于平均高程面的高差達(dá)到300 m時,將引起2.36 m即1個像元的拼接誤差。大多數(shù)情況下地物高差小于300 m,因此地物高差對拼接精度的影響均在一個像素以內(nèi),在保證精度的情況下將地物點的高程統(tǒng)一用平均高程替代,能降低地表模型的復(fù)雜度,提高物方拼接的效率。若相鄰兩片TDI CCD的偏場角非常相近乃至相等時, Δq取值接近于0,則在任何地形條件下的影響均可忽略不計。因此在非共線TDI CCD相機設(shè)計時,應(yīng)盡量縮小兩列TDI CCD沿軌方向的間距,以減弱地形起伏等因素對拼接精度的影響。

圖4 地物高差對拼接精度的影響Fig.4 Influence of the elevation error on stitching accuracy

3.3 方法流程

基于以上原理,本文非共線TDI CCD成像數(shù)據(jù)內(nèi)視場拼接方法流程如圖5所示。

(1)定義物方投影面。將物方投影面定義在成像區(qū)域所對應(yīng)的通用橫軸墨卡托(UTM)投影坐標(biāo)系下,高程取WGS-84大地坐標(biāo)系下的攝影區(qū)域平均橢球高H。

(2)建立從拼接影像到原始影像的像點坐標(biāo)映射關(guān)系。

① 基于原始影像的傳感器幾何模型、高程值H及UTM投影正算,分別計算原始影像第一列的首、尾像素在物方投影面上投影點坐標(biāo),將兩個投影點的連線方向定義為推掃成像時的目標(biāo)軌跡方向。

② 基于原始影像的傳感器幾何模型、高程值H及UTM投影正算,確定原始影像在物方投影面上的投影覆蓋區(qū)域AOP(area of projection)。

③ 確定拼接影像在物方投影面上的投影覆蓋區(qū)域AOI(area of interest),沿推掃成像時的目標(biāo)軌跡方向,確定AOP的最小外廓矩形AOI,即拼接影像在物方投影面上的投影覆蓋區(qū)域。

④ 根據(jù)原始影像的地面分辨率確定物方投影面采樣間隔,對AOI劃分等間隔格網(wǎng),將拼接影像的像點與格網(wǎng)點一一對應(yīng),沿目標(biāo)軌跡和垂直于目標(biāo)軌跡方向的格網(wǎng)數(shù)目分別對應(yīng)拼接影像的行數(shù)和列數(shù)。拼接影像與物方投影面之間的坐標(biāo)轉(zhuǎn)換關(guān)系如式(3)所示

由此即建立了從拼接影像到原始影像的像點坐標(biāo)映射關(guān)系。

(3)對原始影像進行重采樣,生成拼接影像。根據(jù)拼接影像到原始影像的像點坐標(biāo)映射關(guān)系,采用間接法影像糾正的方式對原始影像進行重采樣,生成拼接影像。對于高差較大區(qū)域,可在片間重疊區(qū)進一步提取若干連接點,并在一定范圍進行影像局部糾正[16],以達(dá)到子像素級的目視拼接精度。

(4)構(gòu)建拼接影像的成像幾何模型。

圖5 拼接方法流程圖Fig.5 Workflow of the stitching method

4 試驗與分析

為驗證本文方法的效果和可行性,本文采用ZY-1 02C衛(wèi)星HR相機成像數(shù)據(jù)對方法進行試驗驗證。ZY-1 02C衛(wèi)星上所搭載的HR相機每片TDI CCD有4096像元,在沿軌方向上CCD1、CCD3與CCD2有大約2600行的成像時延,在垂軌方向上有大約30像素左右的重疊,相機地面分辨率為2.36 m,幅寬27 km。本文給出了該相機獲取的不同地形條件下的兩景數(shù)據(jù)的試驗情況(圖6、圖7),表1列出了試驗數(shù)據(jù)的相關(guān)信息。

圖6 數(shù)據(jù)1試驗結(jié)果Fig.6 Experimental results of data 1

圖7 數(shù)據(jù)2試驗結(jié)果Fig.7 Experimental results of data 2

表1 試驗數(shù)據(jù)基本信息Tab.1 Specification of test data

本文分別采用目視(定性)評價以及統(tǒng)計分析拼接影像的RFM擬合精度及其像點定位中誤差(定量),對視場拼接的效果以及精度進行評價。

4.1 視場拼接及目視評價

對表1中的兩景數(shù)據(jù)分別進行拼接試驗,檢查所生成的拼接影像是否滿足目視無縫的要求,即在拼接區(qū)域是否存在錯位現(xiàn)象。

數(shù)據(jù)1(圖6)所對應(yīng)的北京地區(qū)地勢較為平坦,地物高差所引起的拼接誤差很小,遠(yuǎn)在1個像素以內(nèi);數(shù)據(jù)2(圖7)所對應(yīng)的河南登封地區(qū)屬于山地和高山地地形,地物相對于平均高程的高差最大達(dá)到了1100 m以上,根據(jù)地物高差對拼接精度的影響可知,由此會引起最大3.67個像素的局部拼接誤差。為了克服該影響,在片間重疊區(qū)提取了若干連接點,并在拼接區(qū)域進行局部像面糾正[16]從而實現(xiàn)子像素級的目視拼接精度。數(shù)據(jù)1和數(shù)據(jù)2拼接前后的影像分別如圖6(b)和圖7(b)所示,由拼接區(qū)域的放大圖可以看出,拼接影像達(dá)到了目視無縫的精度要求。

4.2 拼接影像RFM擬合精度評價

按照本文3.3節(jié)步驟4的描述構(gòu)建上述兩景拼接影像的RFM,評價RFM對等效攝影幾何模型的擬合精度。采用256像素×256像素的像面格網(wǎng),高程分層為5,這些格網(wǎng)點全部作為控制點;再平移半個格網(wǎng)距離、格網(wǎng)大小不變,再計算出另一批虛擬格網(wǎng)點作為檢查點?？刂泣c用于求解分母不相等的三階RFM模型參數(shù)[20],檢查點用于統(tǒng)計RFM擬合精度。如表2所示,對兩景拼接影像構(gòu)建RFM的擬合精度均在0.000 3個像素以內(nèi),完全可滿足影像生產(chǎn)應(yīng)用的需求。

表2 拼接影像RFM擬合精度評價Tab.2 Fitting error for RFM of stitched images像素

4.3 拼接影像幾何精度分析

本文采用統(tǒng)計像點定位中誤差的方式對拼接影像的幾何精度進行分析評價。其過程如下:

(1)利用提供的參考數(shù)據(jù),在拼接影像范圍內(nèi)提取均勻分布的控制點。對于數(shù)據(jù)1,在拼接影像范圍選取了268個控制點;對于數(shù)據(jù)2,在拼接影像范圍內(nèi)選取了205個控制點。

(2)基于拼接影像的RFM和控制點的地面經(jīng)緯度坐標(biāo)反投影計算其像點坐標(biāo),與其量測的像點坐標(biāo)進行比較求差,得到無控制點情況下像點定位誤差。如表3所示,經(jīng)統(tǒng)計,數(shù)據(jù)1行和列兩個方向上的像點定位中誤差分別為1.92個像素和2.46個像素,數(shù)據(jù)2行和列兩個方向上的像點定位中誤差分別為1.57個像素和2.47個像素。

(3)利用控制點數(shù)據(jù)求解拼接影像RFM的像方仿射變換補償參數(shù)[21],進而統(tǒng)計有控制點情況下像點定位誤差。如表3所示,經(jīng)統(tǒng)計,數(shù)據(jù)1行和列兩個方向上的像點定位中誤差分別為2.37個像素和1.69個像素,數(shù)據(jù)2行和列兩個方向上的像點定位中誤差分別為1.42個像素和2.44個像素。

根據(jù)方法的原理,拼接影像無控制點幾何精度與拼接過程中所采用的原始影像嚴(yán)格成像幾何模型的精度是一致的,通過引入幾何檢校對模型誤差進行補償,在實現(xiàn)子像素級目視無縫拼接精度的同時,還可以消除拼接影像中的幾何畸變。因此基于該方法生成的拼接影像,用戶可以無須控制點或只選取極少量控制點就能獲得具有較高幾何質(zhì)量的拼接影像,滿足后續(xù)應(yīng)用的需求。試驗也表明,對于數(shù)據(jù)1和數(shù)據(jù)2,拼接影像無控制點情況下的像點定位中誤差在2.5個像素以內(nèi),即使利用控制點對拼接影像的RFM誤差進行補償之后,反映拼接影像內(nèi)部畸變的中誤差指標(biāo)值也沒有明顯變化,這說明通過在原始影像的傳感器幾何模型中引入相機幾何檢校結(jié)果已經(jīng)較好地去除了系統(tǒng)性的幾何畸變,因此在有控制點情況下,控制點對內(nèi)部畸變的改善效果不明顯。

表3 拼接影像幾何精度評價Tab.3 Geometric accuracy evaluation of the stitched images 像素

4.4 與像方拼接方法的比較

本文還將基于行積分時間歸一化的方法[12]以及基于分段仿射變換的方法[12,16]等兩種像方方法與本文方法進行了對比。這兩種像方拼接方法對數(shù)據(jù)1和數(shù)據(jù)2所生成的拼接影像均能滿足目視無縫的要求。但由于像方拼接后影像實際上失去了攝影測量的投影幾何特性,因此只能對其構(gòu)建了一種近似等效攝影幾何模型[12],本文在此基礎(chǔ)上進行了精度統(tǒng)計,統(tǒng)計的像點定位中誤差見表3。從表3可以看出,無控制點情況下,基于行積分時間歸一化的方法像點定位中誤差最大接近21個像素,基于分段仿射變換的方法像點定位中誤差最大也接近12個像素,表明影像內(nèi)部幾何精度較低。而利用控制點采用像方仿射變換模型補償RFM系統(tǒng)誤差之后,像點定位中誤差有所改善,基于行積分時間歸一化的方法最大仍接近7個像素,基于分段仿射變換的方法中誤差則與本文方法接近。采用像方多項式變換模型補償RFM系統(tǒng)誤差之后,像點定位中誤差均進一步減小,與本文方法接近。但由于選取的控制點較多,在實際應(yīng)用中工作量太大,因此本文方法更具優(yōu)越性。

5 結(jié) 論

本文從相機的嚴(yán)格成像幾何模型出發(fā),提出了一種非共線TDI CCD成像數(shù)據(jù)內(nèi)視場拼接方法。該方法根據(jù)拼接影像像點在原始影像上的位置以及原始影像的傳感器幾何模型,得到拼接影像像點的等效攝影光線,由此構(gòu)建拼接影像的等效攝影幾何模型。為了驗證方法的效果及可行性,本文用國產(chǎn)ZY-1 02C衛(wèi)星高分辨率相機的成像數(shù)據(jù)進行了試驗驗證。試驗結(jié)果及定性、定量評價表明,本文方法切實可行,拼接影像的生成過程等效于沿衛(wèi)星推掃成像的目標(biāo)軌跡方向?qū)υ加跋襁M行了一次準(zhǔn)系統(tǒng)幾何校正,能滿足影像后續(xù)處理和應(yīng)用的需求。目前,該方法已成功應(yīng)用于國產(chǎn)資源衛(wèi)星地面處理系統(tǒng)。本文方法不僅適用于ZY-1 02C衛(wèi)星高分辨率相機的成像數(shù)據(jù)的視場拼接處理,也同樣適用于其他非共線設(shè)計的TDI CCD相機高分辨率成像數(shù)據(jù)的視場拼接處理。如果具有高精度的相機幾何定標(biāo)數(shù)據(jù),可進一步提高本文方法拼接處理的幾何精度。

[1] HUANG Qiaolin,JIANG Wei.A Study of the New Direction of Space-borne Hi-resolution Optical Remote Sensor[J].Spacecraft Recovery&Remote Sensing,2007,28(4): 48-50.(黃巧林,姜偉.高分辨率航天光學(xué)遙感器發(fā)展新思路研究[J].航天返回與遙感,2007,28(4):48-50.)

[2] YANG Bingxin.Characteristics and Main Specifications of IKONOS and QuickBird2 Satellite Camera:Some Points for Developing such Like Satellite Camera[J].Spacecraft Recovery&Remote Sensing,2002,23(4):14-16.(楊秉新.美國IKONOS和QuickBird-2衛(wèi)星相機的主要性能和特點分析及看法[J].航天返回與遙感,2002,23(4):14-16.)

[3] UPDIKE T,COMP C.Radiometric Use of World View-2 Imagery:Technical Note[DB/OL].Colorado:[s.n.],2010-11-01.http:∥www.digitalglobe.com/sites/default/files/.

[4] JACOBSEN K.Calibration of Optical Satellite Sensors[C]∥Proceedings of International Calibration and Orientation Workshop EuroCOW.Casteldefels:[s.n.],2006.

[5] MATERNE A,BARDOUX A,GEOFFRAY H,et al.Backthinned TDI CCD Image Sensor Design and Performance for the Pleiades High Resolution Earth Observation Satellites [C]∥Proceedings of Sixth International Conference on Space Optics.Noordwijk:European Space Agency,2006.

[6] ZHOU Shengli.Analysis of Three-linear TDI CCD Array Application in Cartography Satellite[J].Spacecraft Engineering,2007,16(4):19-22.(周勝利.三線陣TDI CCD在測圖衛(wèi)星中的應(yīng)用分析[J].航天器工程,2007,16 (4):19-22.)

[7] YANG Bo,WANG Mi.On-orbit Geometric Calibration Method of ZY-1 02C Panchromatic Camera[J].Journal of Remote Sensing,2013,17(5):1175-1190.(楊博,王密.資源一號02C衛(wèi)星全色相機在軌幾何定標(biāo)方法[J].遙感學(xué)報,2013,17(5):1175-1190.)

[8] LI Deren.China’s First Civilian Three-line-array Stereo Mapping Satellite:ZY-3[J].Acta Geodaetica et Cartographica Sinica,2012,41(3):317-322.(李德仁.我國第一顆民用三線陣立體測圖衛(wèi)星:資源三號衛(wèi)星[J].測繪學(xué)報,2012,41(3):317-322.)

[9] JIANG Yonghua,ZHANG Guo,TANG Xingming,et al.High Accuracy Geometric Calibration of ZY-3 Three Line Image[J].Acta Geodaetica et Cartographica Sinica,2013, 42(4):523-529.(蔣永華,張過,唐新明,等.資源三號測繪衛(wèi)星三線陣影像高精度幾何檢校[J].測繪學(xué)報,2013, 42(4):523-529.)

[10] TANG X M,ZHANG G,ZHU X Y,et al.Triple Linear-array Image Geometry Model of Ziyuan-3 Surveying Satellite and Its Validation[J].International Journal of Image and Data Fusion(IJIDF),2013,4(1):33-51.

[11] YANG Hua,GUO Yue,FU Ruimin.Study on Field Butting of TDI CCD[J].Optical Technique,2003,29(2):226-228.(楊樺,郭悅,伏瑞敏.TDI CCD的內(nèi)視場拼接[J].光學(xué)技術(shù),2003,29(2):226-228.)

[12] H U Fen.Research on Inner Fo V Stitching Theories and Algorithms for Sub-images of Three Non-collinear TDI CCD Chips[D].Wuhan:Wuhan University,2010.(胡芬.三片非共線TDI CCD成像數(shù)據(jù)內(nèi)視場拼接理論與算法研究[D].武漢:武漢大學(xué),2010.)

[13] LONG Xiaoxiang,WANG Xiaoyan,ZHONG Huimin.Analysis of Image Quality and Processing Method of a Space-borne Focal Plane View Splicing TDI CCD Camera[J].Scientia Sinica:Informationis,2011,41(sup):19-31.(龍小祥,王小燕,鐘慧敏.星載焦面視場拼接TDI CCD相機成像質(zhì)量及處理方法分析[J].中國科學(xué):信息科學(xué),2011,41(增刊):19-31.)

[14] LI Shiwei,LIU Tuanjie,WANG Hongqi.Image Mosaic for TDI CCD Push-broom Camera Image Based on Image Matching[J].Remote Sensing Technology and Application,2009,24(3):374-378.(李世威,劉團結(jié),王宏琦.基于圖像匹配的CBERS-02B衛(wèi)星HR相機圖像拼接方法[J].遙感技術(shù)與應(yīng)用,2009,24(3):374-378.)

[15] LU Jinbo.Automatic Mosaic Method of Large Field View and Multi-channel Remote Sensing Images of TDI CCD Cameras[D].Jilin:Changchun Institute of Optics,Fine Mechanics and Physics of CAS,2010.(祿金波.TDI CCD大視場多通道圖像自動拼接方法[D].吉林:中國科學(xué)院長春光學(xué)精密機械與物理研究所,2010.)

[16] HU Fen,WANG Mi,JIN Shuying.Satellite-borne Threenon-colinear TDI CCD(time delayed and integration chargecoupled device)Image Splicing Method Based on Segmented Affine Transformation:CN 101799293 B[P].2011-11-02.(胡芬,王密,金淑英.基于分段仿射變換的星載三片非共線TDI CCD影像拼接方法:中國,101799293 B[P].2011-11-02.)

[17] ZH ANG Guo,LIU Bin,JIANG Wanshou.Inner Fo V Stitching Algorithm of Spaceborne Optical Sensor Based on the Virtual CCD Line[J].Journal of Image and Graphics, 2012,17(6):696-701.(張過,劉斌,江萬壽.虛擬CCD線陣星載光學(xué)傳感器內(nèi)視場拼接[J].中國圖象圖形學(xué)報, 2012,17(6):696-701.)

[18] PAN Hongbo,ZHANG Guo,TANG Xinming,et al.The Geometrical Model of Sensor Corrected Products for ZY-3 Satellite[J].Acta Geodaetica et Cartographica Sinica, 2013,42(4):516-522.(潘紅播,張過,唐新明,等.資源三號測繪衛(wèi)星傳感器校正產(chǎn)品幾何模型[J].測繪學(xué)報, 2013,42(4):516-522.)

[19] HU F,WANG M,JIN S Y.An Algorithm for Mosaicking Non-collinear TDI CCD Chip Images Based on Reference Plane in Object-space[C]∥Proceedings of the 6th International Symposium on Digital Earth(ISDE6).Beijing: SPIE,2009.

[20] ZH ANG Guo.Rectification for High Resolution Remote Sensing Image under Lack of Ground Control Points[D].Wuhan:Wuhan University,2005.(張過.缺少控制點的高分辨率衛(wèi)星遙感影像幾何糾正[D].武漢:武漢大學(xué),2005.)

[21] 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.

(責(zé)任編輯:叢樹平)

An Inner FOV Stitching Method for Non-collinear TDI CCD Images

PAN Jun1,HU Fen2,WANG Mi1,JIN Shuying1,LI Guoyuan2
1.State Key Laboratory of Information Engineering in Surveying,Mapping and Remote Sensing,Wuhan University, Wuhan 430079,China;2.Satellite Surveying and Mapping Application Center,National Administration of Surveying,Mapping and Geoinformation,Beijing 101300,China

High-quality inner FOV(field of view)stitching is a necessary step for post-process and application of images captured by spaceborne non-collinear TDI CCD camera.This paper presents an inner FoV stitching method for non-collinear TDICCDimages.The method utilizes the projection plane in object space and the sensor geometric model of the camera to establish the coordinates transformation relation between pixels of the stitched image and pixels of the orginal image,and then performs image resampling to generate the stitched image by indirect image geometric rectification.After describing the strict geometric models of non-collinear TDI CCDimages,the basic principle,the potential error sources and the workflow of the method are discussed in detail.Finally,experimental verification is carried out on the high-resolution (HR)camera images of ZY-1 02C satellite,and the accuracy evaluation of the generated stitched images indicates the feasibility of the presented method.

non-collinear TDI CCD;inner FOV stitching;sensor geometric model;projection plane in object-space;evaluation of accuracy

PAN Jun(1979—),male,PhD,associate professor,majors in image quality improvement,relative radiometric normalization and seamless mosaicking for high resolution remote sensed images.

HU Fen

TP751

A

1001-1595(2014)11-1163-09

國家973計劃(2014CB744201;2012CB719901);國家自然科學(xué)基金(41371430;40901210);全國博士學(xué)位論文作者專項資金(201249);長江學(xué)者和創(chuàng)新團隊發(fā)展計劃(IRT1278);測繪地理信息公益性行業(yè)科研專項(201412007)

2013-12-06

潘俊(1979—),男,博士,副教授,研究方向為高分辨率遙感圖像質(zhì)量改善、相對輻射校正以及無縫鑲嵌等。

E-mail:panjun1215＠whu.edu.cn

胡芬

E-mail:huf＠sasmac.cn

PAN Jun,HU Fen,WANG Mi,et al.An Inner FOV Stitching Method for Non-collinear TDI CCD Images[J].Acta Geodaetica et Cartographica Sinica,2014,43(11):1165-1173.(潘俊,胡芬,王密,等.一種非共線TDI CCD成像數(shù)據(jù)內(nèi)視場拼接方法[J].測繪學(xué)報, 2014,43(11):1165-1173.)

10.13485/j.cnki.11-2089.2014.0180

修回日期:2014-04-01