余岸竹，姜　挺，龔　輝，郭文月，江剛武

1. 信息工程大學(xué)地理空間信息學(xué)院，河南 鄭州 450052； 2. 西安測繪研究所，陜西 西安 710054

?

線陣衛(wèi)星遙感影像外方位元素對偶四元數(shù)求解法

余岸竹1，姜挺1，龔輝2，郭文月1，江剛武1

1. 信息工程大學(xué)地理空間信息學(xué)院，河南 鄭州 450052； 2. 西安測繪研究所，陜西 西安 710054

Foundation support： The National Natural Science Foundation of China (Nos.41201477;41471387;41301526)

摘要：基于對偶四元數(shù)可統(tǒng)一描述位置與姿態(tài)的特點，提出了利用對偶四元數(shù)求解線陣衛(wèi)星遙感影像外方位元素的方法。該方法使用對偶四元數(shù)的實部描述傳感器的姿態(tài)，并利用對偶部和實部共同描述成像傳感器的位置。通過對位置和姿態(tài)的內(nèi)插建立了基于對偶四元數(shù)的外方位元素模型。為減少運算，將球面線性插值進(jìn)行化簡，進(jìn)而建立了基于線性插值的成像幾何模型。為求解外方位元素，首先對共線條件方程進(jìn)行了線性化，然后通過矩陣微分運算推導(dǎo)了線元素的虛擬觀測方程，并根據(jù)誤差傳播定律確定其權(quán)值，最終采用具有約束條件的參數(shù)平差法求解外方位元素。試驗結(jié)果表明本文方法正確可靠，與采用歐拉角和單位四元數(shù)的外方位元素求解方法相比，有更高的參數(shù)解算精度，同時也表明了準(zhǔn)確的初值和虛擬觀測方程對外方位元素求解的必要性。

關(guān)鍵詞：對偶四元數(shù)；線陣遙感影像；虛擬觀測方程；外方位元素；單位四元數(shù)

空間后方交會是衛(wèi)星攝影測量幾何處理中的重要環(huán)節(jié)，是衛(wèi)星遙感影像幾何定位、區(qū)域網(wǎng)平差和自檢校光束法平差等計算的關(guān)鍵步驟，其基本思想是利用合適的外方位元素模型描述成像傳感器的位置與姿態(tài)，建立線陣影像的成像幾何模型，對共線條件方程進(jìn)行線性化后迭代求解外方位元素。當(dāng)前，大多數(shù)衛(wèi)星采用線陣CCD傳感器成像，傳感器在一定時間內(nèi)的外方位元素通常被描述為某一成像行關(guān)于時間的多項式函數(shù)形式[1-3]，或者若干離散成像時刻外方位元素的插值形式[4-6]。對于運行比較平穩(wěn)的衛(wèi)星，這兩類模型均可較好地描述成像傳感器的成像位置與姿態(tài)，參數(shù)之間的相關(guān)性成為影響后方交會精度的主要因素，其成因是線陣傳感器的焦距較大、CCD線陣較短，使得成像時刻傳感器的瞬時視場角較小，共線條件方程經(jīng)線性化后產(chǎn)生了線元素和角元素的相關(guān)。

為了解決此類問題，學(xué)者們試圖從方程解算方法和外方位元素建模這兩個方面來提高外方位元素的求解精度與穩(wěn)定性。在方程解算方面，文獻(xiàn)[7]使用廣義嶺估計法求解外方位元素；文獻(xiàn)[8]采用嶺-壓縮組合估計的方法進(jìn)行了參數(shù)解算；文獻(xiàn)[9—10]對虛擬觀測方程中權(quán)值的確定進(jìn)行了改進(jìn)。以上文獻(xiàn)中均利用歐拉角描述角元素，因而這些算法僅能提高參數(shù)解算的精度，卻無法消除或者降低參數(shù)之間的相關(guān)性。在外方位元素建模方面，文獻(xiàn)[11]使用姿態(tài)四元數(shù)代替了歐拉角，建立了基于四元數(shù)球面線性(spherical linear interpolation,SLERP)的線陣影像外方位元素模型，實現(xiàn)了外方位元素的高精度解算；文獻(xiàn)[12]利用四元數(shù)微分理論，提出了一種基于四元數(shù)微分方程(quaternion differential equation,QDE)的外方位元素模型，使用Tikhonov正則化法進(jìn)一步提高了外方位元素求解的精度；文獻(xiàn)[13]利用信噪比定量比較了各種四元數(shù)外方位元素模型在參數(shù)求解中的差異。雖然利用四元數(shù)可以有效降低參數(shù)之間的相關(guān)性，但仍是將平移與旋轉(zhuǎn)分開考慮，位置與姿態(tài)之間沒有嚴(yán)格的約束條件。

近些年，部分學(xué)者引入對偶四元數(shù)對外方位元素進(jìn)行描述，利用其同時描述平移與旋轉(zhuǎn)的特性，已在遙感影像幾何處理中取得了初步的成果：文獻(xiàn)[14]利用對偶四元數(shù)實現(xiàn)了面陣影像的空間后方交會；文獻(xiàn)[15]驗證了對偶四元數(shù)在區(qū)域網(wǎng)平差求解中的優(yōu)勢；文獻(xiàn)[16]利用對偶四元數(shù)進(jìn)行了面陣影像的直接定向；文獻(xiàn)[17]利用線性蒙皮混合算法進(jìn)行外方位元素建模，提高了GeoEye-1衛(wèi)星遙感影像的立體定位結(jié)果，但是該衛(wèi)星隱藏了定軌和定姿數(shù)據(jù)，因而該文獻(xiàn)未研究定軌、定姿數(shù)據(jù)在后方交會中的應(yīng)用，難以保證迭代在真值附近進(jìn)行。

實際上，線陣衛(wèi)星遙感影像的后方交會計算依賴于未知數(shù)的初值[18]，應(yīng)當(dāng)在外方位元素準(zhǔn)確建模的基礎(chǔ)上，充分利用衛(wèi)星定軌、定姿設(shè)備的觀測結(jié)果，準(zhǔn)確給定各類觀測方程的權(quán)值[19]，才能確保外方位元素的正確求解。本文試圖利用對偶四元數(shù)同時描述成像傳感器的位置和姿態(tài)，引入虛擬觀測方程，進(jìn)一步提高外方位元素的求解精度。

1對偶四元數(shù)

對偶四元數(shù)由單位四元數(shù)和對偶數(shù)理論發(fā)展起來，其定義為

q=r+εs

(1)

(2)

因此對偶四元數(shù)中僅有6個自由度，與平移和旋轉(zhuǎn)未知數(shù)的數(shù)量相同。利用對偶四元數(shù)的性質(zhì)[20]，旋轉(zhuǎn)矩陣R與平移矢量t可利用對偶四元數(shù)的實部與對偶部分別表示為[16]

(3)

(4)

(5)

式(5)與單位四元數(shù)[11]所構(gòu)成的旋轉(zhuǎn)矩陣是一樣的。同理可知，平移矢量t的具體表達(dá)式為

2Trs=t

(6)

式中，(XS,YS,ZS)T為平移矢量的3個分量，在傳統(tǒng)攝影測量中表示線元素。當(dāng)線元素和姿態(tài)歐拉角已知時，對偶四元數(shù)實部r可按歐拉角計算單位四元數(shù)的方法確定[21]，對偶部s可由利用矢量t和r計算為

(7)

由式(5)和式(6)可知，兩個坐標(biāo)系之間的旋轉(zhuǎn)矩陣R可以用由對偶四元數(shù)的實部r計算得到，坐標(biāo)系原點之間的平移矢量t可由實部r與對偶s共同計算得到，因而利用對偶四元數(shù)建立線陣影像的外方位元素模型是可行的。

2基于對偶四元數(shù)的外方位元素求解

2.1基于對偶四元數(shù)的成像幾何模型

設(shè)q1=r1+εs1和q2=r2+εs2分別為影像首行和末行外方位元素所對應(yīng)的對偶四元數(shù)，rk=[r0krxkrykrzk]T與sk=[s0ksxksykszk]T(k=1,2)分別表示影像首行和末行外方位元素的實部和對偶部。在衛(wèi)星成像過程中，衛(wèi)星的位置和姿態(tài)變化比較平穩(wěn)，可利用低階插值計算第i個成像行的外方位元素為

(8)

式中，pt=i/n，n為首末兩行之間成像行的數(shù)量；t1與t2分別表示首末兩行的外方位線元素；C1和C2為四元數(shù)的插值系數(shù)。為了實現(xiàn)四元數(shù)在插值時角度勻速變化[6,11]，通常采用四元數(shù)球面線性插值，即

顯然球面插值會導(dǎo)致C1和C2中包含三角函數(shù)與反三角函數(shù)的復(fù)合運算，使得線性化求解外方位元素的表達(dá)式過于復(fù)雜。由于衛(wèi)星平臺運行較為平穩(wěn)，首末兩行的姿態(tài)四元數(shù)差異很小，因而θ→0+，sinθ與θ是等價無窮小，此時插值系數(shù)可化簡為

(9)

式(9)表明，當(dāng)衛(wèi)星運行比較平穩(wěn)，兩成像行的成像間隔不太長時，可利用線性插值代替球面線性插值計算任意成像行的姿態(tài)。本文將利用式(9)所示的線性表達(dá)式構(gòu)建的外方位元素模型稱為線性插值的對偶四元數(shù)模型(linear interpolation based dual quaternion model,LDQ)，將利用球面線性插值構(gòu)建的外方位元素模型稱為球面線性插值的對偶四元數(shù)模型(SLERP based dual quaternion model,SDQ)。

本文以LDQ模型為例推導(dǎo)求解外方位元素的誤差方程式。不考慮主點偏移與焦距變化時，第個成像行的瞬時構(gòu)像方程可以寫成

(10)

式中，f為相機(jī)焦距；(X,Y)為像點坐標(biāo)，可利用傳感器沿軌道方向的指向角ψX和垂直于軌道方向的指向角ψY表示為(-ftanψY,ftanψX)；系數(shù)aji、bij、cji(j=1,2,3)為該行對應(yīng)的旋轉(zhuǎn)矩陣中的元素，可內(nèi)插計算ri后，由式(5)計算旋轉(zhuǎn)矩陣得到；ti=[XsiYsiZsi]T是成像行攝影中心在物方坐標(biāo)系內(nèi)的坐標(biāo)，可由式(6)和式(8)聯(lián)立計算得到。式(10)即為利用LDQ模型構(gòu)建的線陣影像成像幾何方程。

為求解外方位元素，將式(10)在初值處展開至一次項，可得誤差方程為

(11)

(12)

誤差方程(11)的矩陣形式為

V=CX-lE

(13)

此外，影像首行和末行對應(yīng)的對偶四元數(shù)均滿足式(2)，對該式進(jìn)行線性化并寫成矩陣形式可得

BX+W=0

(14)

式中

誤差方程(13)與條件方程(14)構(gòu)成了外方位元素計算的基本平差模型。

2.2對偶四元數(shù)虛擬觀測方程的建立

虛擬觀測方程在攝影測量中有很高的實際應(yīng)用價值，在線陣影像的后方交會和光束法平差中被廣泛采用。由于該方法可以將外方位元素的改正數(shù)限制在一個合理的范圍內(nèi)，因而于有偏估計這種方法更為嚴(yán)密[18]。當(dāng)前大多數(shù)衛(wèi)星可以提供高精度的定軌數(shù)據(jù)，而定姿數(shù)據(jù)通常含有系統(tǒng)誤差，且誤差的構(gòu)成較為復(fù)雜[22]。本文僅對線元素設(shè)置虛擬觀測方程，適當(dāng)放寬姿態(tài)改正數(shù)的限制。

對式(6)求微分可得

dt=2×[MrMs]X

(15)

式中，Mr和Ms的表達(dá)式分別為

對線元素增設(shè)的虛擬觀測方程可表示為

(16)

式中，PM為虛擬觀測方程的權(quán)矩陣。本文在迭代過程中動態(tài)調(diào)整權(quán)值[18]，即先利用驗前權(quán)值進(jìn)行預(yù)平差，并進(jìn)行方差分量估計值，然后計算第次迭代中參數(shù)di在虛擬觀測方程中的權(quán)值為

(17)

(18)

利用nG個地面控制點，按照式(18)總共可以列出2nG+6個誤差方程和4個條件方程，應(yīng)有nG≥3。當(dāng)nG<6時，式(18)實質(zhì)上是將定軌數(shù)據(jù)作為控制條件參與方程的求解，僅在精確給定初值時正確收斂，通常情況下應(yīng)有nG≥6。利用約束條件的參數(shù)平差方法可求得改正數(shù)X，通過迭代即可求得利用對偶四元數(shù)描述的外方位元素值。

3試驗與分析

為驗證提出方法的正確性與有效性，本文將使用兩組高分辨率衛(wèi)星遙感影像進(jìn)行對比試驗驗證。第1組試驗數(shù)據(jù)(影像Ⅰ)為河南某地SPOT-5 HRS立體影像，共有13 376條掃描行，每個掃描行有12 000個像元，影像內(nèi)包含25個平面高程控制點(圖1(a))。第2組試驗數(shù)據(jù)(影像Ⅱ)為中國西部某地丘陵地區(qū)的SPOT-5 HRS立體影像，影像內(nèi)共有12 000條掃描行，共包含10個平高控制點(圖1(b))。

3.1試驗1——影像I的外方位元素求解

本試驗用于比較對偶四元數(shù)外方位元素模型、歐拉角外方位模型和單位四元數(shù)外方位元素模型在外方位元素求解中的差異。為了驗證式(9)的合理性，本試驗還將比較SDQ模型和LDQ

模型在后方交會計算中的區(qū)別。試驗評價指標(biāo)有兩種：一是利用求解出的外方位元素反算地面點的像點坐標(biāo)，并計算像點重投影誤差的中誤差；二是利用兩張影像的外方位元素進(jìn)行空間前方交會，計算三維坐標(biāo)殘差的中誤差。試驗中分別采用6、7、14和20個控制點求解外方位元素。表1和表2為不同控制點數(shù)量時檢查點的像點坐標(biāo)重投影的中誤差及單位權(quán)方差統(tǒng)計結(jié)果；表3和表4為求解外方位元素值按空間前方交會計算地面點的三維坐標(biāo)，計算檢查點坐標(biāo)殘差的誤差統(tǒng)計結(jié)果。

圖1　SPOT-5 HRS影像控制點分布示意圖　Fig.1　Distribution of ground control points for SPOT-5 HRS images

影像控制點數(shù)檢查點數(shù)x/像素y/像素EulerSLERPQDESDQLDQEulerSLERPQDESDQLDQ前視影像6192.7121.7991.7661.7351.7221.5291.5261.4471.5301.5317181.7341.8561.8471.7931.7761.4181.3171.4811.3221.32614112.2021.6671.5981.5181.4691.2561.2571.4521.2611.2622051.1961.3031.2761.2071.1950.7890.7540.7670.7530.752后視影像619不收斂1.5821.4141.5541.538不收斂1.2611.4391.2611.2617181.5761.5171.5361.4951.4811.7681.2381.3351.2381.23814111.6761.6561.4851.6141.5851.5541.4441.4871.4441.4442050.9291.2100.9971.1421.0921.4361.1931.1831.1911.191

注：Euler為利用歐拉角求解外方位元素的算法，SLERP為文獻(xiàn)[11]中利用四元數(shù)球面線性插值求解外方位元素算法，QDE為文獻(xiàn)[12]中的四元數(shù)微分方程算法，下同。

從試驗1的結(jié)果可知：

(1) 引入了虛擬觀測方程后，除6個點時Euler算法在后視影像上未能收斂，其余情況下所有算法均可正確收斂。當(dāng)7個和14個控制點參與平差時，Euler算法的立體定位精度較低。主要因為傳統(tǒng)外方位元素描述中，角元素、線元素之間的相關(guān)性很強，誤差方程系數(shù)矩陣呈現(xiàn)嚴(yán)重病態(tài)。當(dāng)控制點數(shù)量較少時，方程不能穩(wěn)定求解，使得第1次迭代的改正數(shù)比較大，求解過程未能在初值附近收斂，反而得到臨近區(qū)間的局部最優(yōu)值。這也從側(cè)面表明了合理的初值和虛擬觀測方程的必要性。

(2) SDQ算法和LDQ算法的立體定位平面精度略高于SLERP算法，三者的高程精度相當(dāng)，僅有厘米級的差別；在6個、7個和14個控制點參與運算時，三者的立體定位精度均優(yōu)于QDE算法，像點重投影的中誤差相當(dāng)，這表明利用對偶四元數(shù)描述外方位元素模型正確可靠，且在外方位元素求解中具有一定優(yōu)勢，主要是因為對偶四元數(shù)的實部與對偶部滿足正交關(guān)系，使得r與s之間的相關(guān)性很低，對參數(shù)求解較為有利。值得注意的是，試驗中QDE的計算精度略低于SLERP算法，這與文獻(xiàn)[13]中結(jié)果有出入，原因是為確保試驗變量唯一，本文試驗中未使用Tikhonov正則化求解參數(shù)，因而平差結(jié)果有所差異。

(3) SDQ算法和LDQ算法的計算結(jié)果基本一樣，立體定位的結(jié)果僅有厘米級的差別，說明利用式(9)簡化球面線性插值是正確可靠的。同時，LDQ算法中共線條件方程經(jīng)線性化后的系數(shù)矩陣中不包括三角函數(shù)運算，在保證精度的同時提高了運算效率和參數(shù)求解精度，這對實際應(yīng)用有重要意義。但是這種近似是有條件的，應(yīng)針對不同的衛(wèi)星影像和不同成像條件區(qū)分對待。

(4) 當(dāng)6、7和14個控制點參與平差解算時，5種算法立體定位的精度并無明顯提高，直到20個點參與平差后精度才達(dá)到1個像元左右。這主要因為試驗所用控制點均在立體環(huán)境下人工量測得到，存在一定的辨識誤差和量測誤差，并且試驗中未針對系統(tǒng)誤差進(jìn)行處理，欲進(jìn)一步提高立體定位的精度，需要顧及控制點的誤差，引入自檢校參數(shù)進(jìn)行自檢校光束法平差計算。

表2影像Ⅰ中不同控制點數(shù)量下平面重投影誤差與單位權(quán)方差統(tǒng)計結(jié)果

Tab.2Statisticsresultofplanerootmeansquareerrorforre-projectionandunitweightvariancewithdifferentnumbersofGCPsforimageⅠ

影像控制點數(shù)檢查點數(shù)平面/像素^σ0/像素EulerSLERPQDESDQLDQEulerSLERPQDESDQLDQ前視影像6193.1142.3592.2832.3132.3040.3912.5732.0070.9940.9657182.2402.2762.3702.2272.2160.4761.8891.5661.1551.13114112.5362.0872.1591.9731.9370.6771.6351.4801.3041.2992051.4331.5051.4891.4231.4120.8771.7621.5651.3741.358后視影像619不收斂2.0232.0182.0021.989不收斂0.9160.8390.5170.5097182.3681.9582.0351.9411.9310.2240.9980.6680.7390.73614112.2852.1972.1012.1652.1440.4411.0510.9200.8710.8602051.7101.6991.5481.6511.6160.4361.6721.2791.1821.163

表3　影像Ⅰ中不同控制點數(shù)量下求解地面三維坐標(biāo)水平方向誤差統(tǒng)計

表4　影像Ⅰ中不同控制點數(shù)量下求解地面三維坐標(biāo)平面和高程誤差統(tǒng)計

3.2試驗2——影像Ⅱ外方位元素求解

本試驗用于比較本文外方位元素模型和文獻(xiàn)[17]中基于線性蒙皮混合的對偶四元數(shù)外方位元素模型(linear blending dual quaternion model,BDQ)在后方交會中的區(qū)別。由于文獻(xiàn)[17]在求解外方位元素時沒有準(zhǔn)確給定初值，且未增設(shè)虛擬觀測方程，試驗中將對本文方法與原始BDQ算法、賦初值的BDQ算法和賦初值并增設(shè)虛擬觀測方程的BDQ算法在求解外方位元素時的差異。分別使用6、7和8個控制點求解影像Ⅱ的外方位元素并統(tǒng)計重投影誤差和立體定位精度，結(jié)果見表5。

表5　影像Ⅱ求解外方位元素的重投影誤差和立體定位精度統(tǒng)計

由試驗2的結(jié)果可知：

(1) 當(dāng)沒有準(zhǔn)確給定初值時，BDQ算法在3次試驗中均不能正確收斂。這是因為線陣衛(wèi)星遙感影像定向參數(shù)的求解過程中通常使用地心地固坐標(biāo)系作為物方坐標(biāo)系[9]，像方坐標(biāo)系和物方坐標(biāo)系之間的旋轉(zhuǎn)矩陣不再近似為單位陣。同時，傳感器的成像姿態(tài)會隨著星下點經(jīng)度和緯度的不同而變化，應(yīng)由定軌、定姿數(shù)據(jù)給定外方位元素初值。

(2) 加入初值后，BDQ算法在8個控制點參與運算時可以收斂，但是重投影誤差較大，立體定位的精度較低。在此基礎(chǔ)上引入虛擬觀測方程后，BDQ算法在各次試驗中均正確收斂，重投影誤差和立體定位結(jié)果也有一定的改善，這表明了求解線陣影像外方位元素時增設(shè)虛擬觀測方程的必要性。

(3) 準(zhǔn)確給定初值并加入虛擬觀測方程之后，BDQ算法與本文方法的精度相近，但是本文方法的精度更高，這種精度的差異是由外方位元素模型的區(qū)別引起的。BDQ算法中的使用插值算法是一種近似算法[23]，對外方位元素的對偶部進(jìn)行插值，物理意義并不明確；本文方法依據(jù)衛(wèi)星成像時運動平穩(wěn)的特性，對線元素進(jìn)行插值來構(gòu)建外方位元素模型，具有明確的物理意義，因而在3組試驗中重投影和立體定位的精度更高。

4結(jié)論

本文利用對偶四元數(shù)構(gòu)建的SDQ模型和LDQ模型可實現(xiàn)線陣衛(wèi)星遙感影像的外方位元素求解。試驗結(jié)果表明，相比于利用歐拉角和利用單位四元數(shù)的外方位元素求解方法，本文方法更具優(yōu)勢；與已有的對偶四元數(shù)外方位元素求解方法相比，本文的外方位元素模型與平差手段更為合理，是線陣衛(wèi)星遙感影像外方位元素求解的一條新途徑。然而在試驗過程中也發(fā)現(xiàn)了本文方法的局限性，控制點的精度和權(quán)初值會影響平差計算的結(jié)果，應(yīng)引入抗差估計等手段進(jìn)一步提高算法的性能。當(dāng)然，欲將對偶四元數(shù)實際應(yīng)用到衛(wèi)星遙感影像的幾何處理當(dāng)中，需要進(jìn)行更為深入和系統(tǒng)的研究，如何利用對偶四元數(shù)進(jìn)行自檢校區(qū)域網(wǎng)平差計算和構(gòu)建定向片模型等問題，都是今后研究的重點內(nèi)容。

參考文獻(xiàn)：

[1]張永生, 鞏丹超, 劉軍, 等. 高分辨率遙感衛(wèi)星應(yīng)用[M]. 北京: 科學(xué)出版社, 2004.

ZHANG Yongsheng, GONG Danchao, LIU Jun, et al. Application of High-resolution Remote Sensing Satellites[M]. Beijing: Science Press, 2004.

[2]袁修孝, 余翔. 高分辨率衛(wèi)星遙感影像姿態(tài)角系統(tǒng)誤差檢校[J]. 測繪學(xué)報, 2012, 41(3): 385-392.

YUAN Xiuxiao, YU Xiang. Calibration of Angular Systematic Errors for High Resolution Satellite Imagery[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(3): 385-392.

[3]POLI D. Modelling of Spaceborne Linear Array Sensors[D]. Zurich: Swiss Federal Institute of Technology, 2005.

[4]HOFMANN O, NAVE P, EBNER H. DPS-A Digital Photogrammetric System for Producing Digital Elevation Models and Orthophotos by Means of Linear Array Scanner Imagery[J]. Photogrammetric Engineering and Remote Sensing, 1984, 50(8): 1135-1142.

[5]劉軍, 王冬紅, 劉敬賢, 等. IMU/DGPS 系統(tǒng)輔助ADS40三線陣影像的區(qū)域網(wǎng)平差[J]. 測繪學(xué)報, 2009, 38(1): 55-60.

LIU Jun, WANG Donghong, LIU Jingxian, et al. IMU/DGPS Supported Block Adjustment of ADS40 Three Line Imagery Data[J]. Acta Geodaetica et Cartographica Sinica, 2009, 38(1): 55-60.

[6]劉軍, 王冬紅, 張永生, 等. 基于單位四元數(shù)的機(jī)載三線陣影像光束法平差[J]. 測繪學(xué)報, 2008, 37(4): 451-457.

LIU Jun, WANG Donghong, ZHANG Yongsheng, et al. Bundle Adjustment of Airborne Three Line Array Imagery Based on Unit Quaternion[J]. Acta Geodaetica et Cartographica Sinica, 2008, 37(4): 451-457.

[7]郭海濤, 張保明, 歸慶明. 廣義嶺估計在解算單線陣CCD衛(wèi)星影像外方位元素中的應(yīng)用[J]. 武漢大學(xué)學(xué)報(信息科學(xué)版), 2003, 28(4): 444-447.

GUO Haitao, ZHANG Baoming, GUI Qingming. Application of Generalized Ridge Estimate to Computing the Exterior Orientation Elements of Satellite Linear Array Scanner Imagery[J]. Geomatics and Information Science of Wuhan University, 2003, 28(4): 444-447.

[8]張艷, 王濤, 朱述龍, 等. 嶺-壓縮組合估計在線陣推掃式影像外定向中的應(yīng)用[J]. 武漢大學(xué)學(xué)報(信息科學(xué)版), 2004, 29(10): 893-896.

ZHANG Yan, WANG Tao, ZHU Shulong, et al. Application of Combined Ridge-stein Estimator to Linear Pushbroom Imagery Exterior Orientation[J]. Geomatics and Information Science of Wuhan University, 2004, 29(10): 893-896.

[9]余岸竹. 遙感衛(wèi)星攝影測量參數(shù)動態(tài)檢測關(guān)鍵技術(shù)研究[D]. 鄭州: 信息工程大學(xué), 2014.

YU Anzhu. Research on the Key Technology of Dynamic Calibration of Photogrammetric Parameters for Spaceborne Sensors[D]. Zhengzhou: Information Engineering University, 2014.

[10]余俊鵬. 基于先驗求權(quán)的虛擬觀測值法及其在線陣影像定向中的應(yīng)用[J]. 測繪科學(xué), 2009, 34(3): 124-126.

YU Junpeng. The Virtual Observation Method with Weight Prior-determined and Its Application in Orientation of the Linear Array Imagery[J]. Science of Surveying and Mapping, 2009, 34(3): 124-126.

[11]JIANG Gangwu, JIANG Ting, GONG Hui, et al. Exterior Orientation of Line -array CCD Images Based on Quaternion Spherical Linear Interpolation[C]∥ISPRS TC VII Symposium-100 Years ISPRS. Vienna: [s.n.], 2010, 38: 71-76.

[12]龔輝, 姜挺, 江剛武, 等. 四元數(shù)微分方程的高分辨率衛(wèi)星遙感影像外方位元素求解[J]. 測繪學(xué)報, 2012, 41(3): 409-416.

GONG Hui,JIANG Ting,JIANG Gangwu, et al. Solution of Exterior Orientation Parameters for High-resolution Satellite Imagery Based on Quaternion Differential Equation[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(3): 409-416.

[13]龔輝. 基于四元數(shù)的高分辨率衛(wèi)星遙感影像定位理論與方法研究[D]. 鄭州: 信息工程大學(xué), 2011.

GONG Hui. Study on Theory and Method of Geopositioning for High-resolution Satellite Imagery Based on Quaternion[D]. Zhengzhou: Information Engineering University, 2011.

[14]姬亭, 盛慶紅, 王惠南, 等. 對偶四元數(shù)單片空間后方交會算法[J]. 中國圖象圖形學(xué)報, 2012, 17(4): 494-503.

JI Ting, SHENG Qinghong, WANG Huinan, et al. Dual Quaternion of Space Resection with Single-image[J]. Journal of Image and Graphics, 2012, 17(4): 494-503.

[15]龔輝, 姜挺, 江剛武, 等. 利用單位對偶四元數(shù)進(jìn)行航空影像區(qū)域網(wǎng)平差解算[J]. 武漢大學(xué)學(xué)報(信息科學(xué)版), 2012, 37(2): 154-159.

GONG Hui, JIANG Ting, JIANG Gangwu, et al. Bundle Block Adjustment of Aerial Imagery Based on Unit Dual Quaternion[J]. Geomatics and Information Science of Wuhan University, 2012, 37(2): 154-159.

[16]龔輝, 江剛武, 姜挺, 等. 基于對偶四元數(shù)的絕對定向直接解法[J]. 測繪科學(xué)技術(shù)學(xué)報, 2009, 26(6): 434-438.

GONG Hui, JIANG Gangwu, JIANG Ting, et al. Close-form Solution of Absolute Orientation Based on Dual Quaternion[J]. Journal of Geomatics Science and Technology, 2009, 26(6): 434-438.

[17]盛慶紅, 姬亭, 劉微微, 等. 對偶四元數(shù)線陣遙感影像幾何定位[J]. 中國圖象圖形學(xué)報, 2012, 17(10): 1319-1326.

SHENG Qinghong, JI Ting, LIU Weiwei, et al. Geo-positioning Line-array CCD Images with Dual Quaternion[J]. Journal of Image and Graphics, 2012, 17(10): 1319-1326.

[18]袁修孝, 曹金山. 高分辨率衛(wèi)星遙感精確對地目標(biāo)定位理論與方法[M]. 北京: 科學(xué)出版社, 2012.

YUAN Xiuxiao,CAO Jinshan. Theory and Method of High Resolution Satellite Remote Sensing for Precise Orientation[M]. Science Press, 2012.

[19]隋立芬, 宋力杰, 柴洪洲. 誤差理論與測量平差基礎(chǔ)[M]. 北京: 測繪出版社, 2010.

SUI Lifen, SONG Lijie, CHAI Hongzhou. Error Theory and Foundation of Surveying Adjustment[M]. Beijing: Surveying and Mapping Press, 2010.

[20]WALKER M W, SHAO Lejun, VOLZ R A. Estimating 3-D Location Parameters Using Dual Number Quaternions[J]. CVGIP: Image Understanding, 1991, 54(3): 358-367.

[21]江剛武. 空間目標(biāo)相對位置和姿態(tài)的抗差四元數(shù)估計[D]. 鄭州: 信息工程大學(xué), 2009.

JIANG Gangwu. A Robust Estimation Using Quaternions for Relative Position of Space Object[D]. Zhengzhou: Information Engineering University, 2009.

[22]王任享. 三線陣CCD影像衛(wèi)星攝影測量原理[M]. 北京: 測繪出版社, 2006.

WANG Renxiang. Satellite Photogrammetric Principle for Three-line-array CCD Imagery[M]. Beijing: Surveying and Mapping Press, 2006.

(責(zé)任編輯：張艷玲)

修回日期： 2015-08-29

Dual Quaternion Method to Solve Exterior Orientation Parameters for Satellite Linear Array Images

YU Anzhu1,JIANG Ting1,GONG Hui2,GUO Wenyue1,JIANG Gangwu1

1. Institute of Surveying and Mapping, Information Engineering University, Zhengzhou 450052,China； 2. Xi’an Research Institute of Surveying and Mapping, Xi’an 710054, China

Abstract：According to the characteristics that dual quaternion can describe the position and orientation simultaneously, a method using dual quaternion to calculate the exterior orientation parameters for satellite linear array images is proposed. In this method the real part of dual quaternion is used to represent attitude of the sensor, and both the real part and the dual part are used to represent the position of the sensor.The model describing exterior parameters is established by interpolating the position and attitude. To reduce calculation, the simplified spherical linear interpolation algorithm is proposed, and the geometric model is established.In order to calculate the exterior parameters, the collinear equation is linearized at first. Then the virtual observation equations for position parameters are deduced according to the theory of matrix derivative, and the weight of the equation is calculated using error propagation principle. The exterior orientation parameters are calculated iteratively using parameter adjustment with constraints at last. Experimental results indicate that the proposed method is reliable and correct, and the precision of parameter estimation is higher than traditional methods using Euler angles or unit quaternion. The necessity of accuracy initial values and virtual observation equations for calculating exterior parameters is also revealed.

Key words：dual quaternion； linear array images;virtual observation equation; exterior orientation parameters; unit quaternion

基金項目：國家自然科學(xué)基金(41201477；41471387；41301526)

中圖分類號：P236

文獻(xiàn)標(biāo)識碼：A

文章編號：1001-1595(2016)02-0186-08

作者簡介：第一 陳敏(1986—)，男，博士，講師，研究方向為多源遙感影像處理與分析。

收稿日期：2015-01-20

First author： CHEN Min(1986—)，male，PhD，lecturer，majors in multi-source remote sensing images processing and analysis.

E-mail： minchen@home.swjtu.edu.cn

引文格式：余岸竹，姜挺，龔輝,等.線陣衛(wèi)星遙感影像外方位元素對偶四元數(shù)求解法[J].測繪學(xué)報，2016,45(2)：186-193. DOI:10.11947/j.AGCS.2016.20140643.

YU Anzhu,JIANG Ting,GONG Hui,et al.Dual Quaternion Method to Solve Exterior Orientation Parameters for Satellite Linear Array Images[J]. Acta Geodaetica et Cartographica Sinica,2016,45(2):186-193. DOI:10.11947/j.AGCS.2016.20140643.