下視影像輔助大傾斜影像外方位角元素的解算

柯　玲，范冬娟，徐　彬

（1.武昌首義學(xué)院基礎(chǔ)科學(xué)部，武漢 430064；2.武漢大學(xué)遙感信息工程學(xué)院，武漢 430079；3.73608部隊，南京 210028）

?

下視影像輔助大傾斜影像外方位角元素的解算

柯玲1，2，范冬娟3，徐彬1

（1.武昌首義學(xué)院基礎(chǔ)科學(xué)部，武漢 430064；2.武漢大學(xué)遙感信息工程學(xué)院，武漢 430079；3.73608部隊，南京 210028）

摘要：本文在已有研究成果基礎(chǔ)上，采用下視和傾斜同步攝影方式獲取多視影像，改進(jìn)了傳統(tǒng)攝影方式，并在下視影像輔助下實現(xiàn)了無控制條件下大傾斜影像的外方位角元素解算。

關(guān)鍵詞：大傾斜影像；下視影像；POS；外方位元素；無控制

范冬娟，女，1976年生，湖北天門人，碩士，高級工程師，主要研究方向是數(shù)字?jǐn)z影測量、GIS。

徐 彬，女，1981年生，山東昌邑人，碩士，講師，主要研究方向是應(yīng)用數(shù)學(xué)方面的研究。

1 引言

相比傳統(tǒng)的豎直攝影方式，傾斜攝影獲取的影像更加接近一般的地面投影，其結(jié)果也更符合人們的日常感受，正受到越來越多的關(guān)注。傾斜攝影目前主要分為兩種形式：一是采用較低傾斜角度，以獲取飛行器側(cè)下方影像為主的傾斜攝影方式；另一種則是采用較大傾斜角度，以獲取飛行器側(cè)面至地平線遠(yuǎn)處影像為主要目的傾斜攝影方式。相對于多鏡頭的低傾斜角度攝影，利用單個相機(jī)進(jìn)行大傾角攝影的主要目的是用于軍事偵察或危機(jī)地域?qū)崟r信息獲取，因其傾角較大且沒有類似多鏡頭相機(jī)那樣的幾何約束關(guān)系，盡管其影像獲取簡單，但地理參考更加復(fù)雜，難以進(jìn)行自動匹配或常規(guī)攝影測量處理［3］。

2 外方位元素解算流程與方法

當(dāng)前，集DGPS（Differential GPS）技術(shù)和慣性導(dǎo)航（Inertial Navigation System）技術(shù)于一體的POS系統(tǒng)已廣泛應(yīng)用于無人機(jī)等低空遙感攝影系統(tǒng)，在近似水平攝影且攝影比例尺較大的情況下（大于1︰10萬比例尺），采用機(jī)載POS進(jìn)行直接的地理定位已可以滿足1︰5萬正射影像圖制作精度要求［9］，獲取的像片姿態(tài)角經(jīng)過轉(zhuǎn)換也可以作為外方位元素的初值使用，在這一基礎(chǔ)上，研究下視影像輔助下大傾角影像外方位角元素的解算方法，實現(xiàn)了無控制條件下大傾斜影像外方位元素求解和地面坐標(biāo)量測。其處理流程如圖1所示。

圖1　無控制條件下傾斜影像處理流程

當(dāng)前在航攝作業(yè)過程中，一般都會隨影像提供高精度的定位定向數(shù)據(jù)，也就是通常所說的POS （Position & Orientation System）數(shù)據(jù)，利用該數(shù)據(jù)可以實現(xiàn)對地面目標(biāo)的直接地理定位，即直接對POS系統(tǒng)獲取的GPS天線相位中心的空間坐標(biāo)及IMU系統(tǒng)獲取的側(cè)滾角、俯仰角、航跡角進(jìn)行數(shù)據(jù)處理，獲取航空影像曝光瞬間的攝站坐標(biāo)Xs, Ys, Zs及3個姿態(tài)角φ,ω,κ共6個外方位元素［10］，進(jìn)而解算地面坐標(biāo)。在6個外方位元素中，利用載波相位GPS動態(tài)定位技術(shù)，攝站位置的平面和高程精度可分別達(dá)到2cm和5cm，完全可以滿足攝影測量的精度要求，而3個角元素的精度則受限于IMU精度的影響，難以直接作為影像的外方位元素使用，因此對外方位元素的求解問題就可以簡化為求解3個角元素φ,ω,κ的問題。在攝影測量處理過程中，遙感影像的外方位元素獲取依賴于地面控制點，通過建立地面控制點和其像點之間的共線關(guān)系來解算外方位元素，因此在無控制條件下如何建立這樣的關(guān)系條件成為解算外方位元素的關(guān)鍵所在。

2.1 模擬控制點

在航攝過程中采用了豎直攝影和傾斜攝影相結(jié)合的方式，在獲取傾斜影像的同時，獲取豎直影像，如圖2所示。

圖2　豎直和傾斜同步攝影示意圖

在圖2中，O0為傾斜攝影相機(jī)的像主點，S0為其攝站點；O1為豎直攝影相機(jī)的像主點，S1為其攝站點，N點是豎直攝影相機(jī)攝站點和像主點的連線與地面的交點，也就是O1對應(yīng)的地面點；n為N在傾斜影像上的對應(yīng)像點，可知，n是O1的同名像點。對于豎直航攝影像來說，攝站點經(jīng)像主點投射線豎直向下，此時可以根據(jù)式（1）得到像主點對應(yīng)的地面點坐標(biāo)。

式中，（X,Y,Z）為N點坐標(biāo)；Xs，Ys，Zs為攝站坐標(biāo)；H為航高，考慮到航高存在的誤差以及地面起伏的影響，可采用DEM內(nèi)插的方式獲得N點Z坐標(biāo)。在求得N點坐標(biāo)后，可將N點作為控制點，再進(jìn)一步得到n點的像點坐標(biāo)，即可建立S0，n，N三點共線的數(shù)學(xué)方程。根據(jù)航攝過程中的航線重疊和旁向重疊情況，在一張傾斜影像上可以建立不少于四組的共線方程（甚至可以建立九組方程），從而滿足外方位角元素答解的要求。

2.2 匹配同名像點

由于傾斜影像和豎直影像之間攝影角度差異較大，導(dǎo)致兩幅影像的投影成像存在較大差別，一般的模板匹配誤差較大，因此采用ASIFT算法完成匹配過程。

圖3　ASIFT中的傾斜角

ASIFT算法是Jean-Michel Morel等人在SIFT算法基礎(chǔ)上，針對兩幅影像投影差別較大的問題而提出的改進(jìn)算法［11］，即以影像的成像模型為基礎(chǔ)，把待匹配的兩幅影像分別看成是相機(jī)在以被攝目標(biāo)為中心，攝影距離為半徑的球面上成的像，這樣就可以通過球面極坐標(biāo)的兩個角度（θ,φ）和相機(jī)自身的旋轉(zhuǎn)角度ω來確定相機(jī)成像時的傾斜矩陣A，如圖3所示。其中的旋轉(zhuǎn)矩陣則按照下列公式進(jìn)行計算：

文獻(xiàn)［12］中即以該式為基礎(chǔ)，把兩幅圖像分別按照一定的間隔對θ,φ,ω三個角度進(jìn)行循環(huán)，對于每一組角度可利用成像模型得到傾斜后的圖像，并對這一系列圖像兩兩組合進(jìn)行匹配，最后統(tǒng)計匹配成功的同名點。

根據(jù)ASIFT算法的計算流程可知，ASIFT算法的計算效率相較SIFT算法的計算效率，要低一個數(shù)量級，因此在實際應(yīng)用過程中還需充分利用已知條件，減少模擬影像數(shù)量和匹配的計算量，提高匹配效率。采用ASIFT算法的匹配步驟如下：

⊙ 基于豎直影像P1，根據(jù)已知的傾斜大角度，模擬得到其傾斜影像P1'。

⊙ 以O(shè)1在模擬傾斜影像上的對應(yīng)像點O1'為中心得到小區(qū)域影像P1"。

⊙ 同樣根據(jù)O1'像點坐標(biāo)在傾斜影像P2上得到以其為中心的小區(qū)域影像P2'。

⊙ 把P1"和P2'采用SIFT算法進(jìn)行匹配，得到O1'匹配結(jié)果；如果經(jīng)過匹配無法得到O1'匹配結(jié)果，則擴(kuò)大P2影像上的區(qū)域范圍，再進(jìn)行匹配。

⊙ 重復(fù)上述步驟，直至完成整個匹配過程。

2.3 外方位元素答解

通過匹配得到控制點對應(yīng)的像點坐標(biāo)后，因為通過POS數(shù)據(jù)可得到精確的攝站坐標(biāo)，因此可認(rèn)為外方位元素中的線元素已知，則可以采用角錐體方法，通過線元素和角元素分別答解的方式，按照（3）式建立方程組，答解外方位角元素［12］［8］。

式中，（XA，YA，ZA）為模擬的控制點坐標(biāo)；（xa，ya）為通過匹配得到的地面控制點對應(yīng)像點坐標(biāo)；Sa和SA分別為攝站點到像點和地面點的距離；a1…a3，b1…b3，c1…c3為旋轉(zhuǎn)矩陣。

3 實驗情況

3.1 實驗區(qū)概況

實驗區(qū)面積約為8平方千米，為高原平地，起伏較小，平均海拔約為1,100米。

3.2 相機(jī)參數(shù)

采用兩個自制相機(jī)（一個相機(jī)豎直攝影，另一個相機(jī)傾斜攝影，傾斜角度45度），具體相機(jī)參數(shù)見表1。

表1　相機(jī)參數(shù)

3.3 飛行參數(shù)

獲取的影像分布如圖4（a）所示，黑色點表示相鄰豎直影像像主點所對應(yīng)地面控制點在傾斜影像上（紅色線框表示）的對應(yīng)像點。一張傾斜像片上，根據(jù)航向多張像片和旁向多條航帶像片攝站點模擬得到的控制點中選取9個作為像片控制點，并從中選擇分布均勻的4個控制點作為解算控制點（如圖4（b）所示），其攝站坐標(biāo)和像坐標(biāo)數(shù)據(jù)羅列在表3中。

表2　飛行參數(shù)

表3　四個模擬控制的坐標(biāo)數(shù)據(jù)

解算得到的外方位元素如表4所列。

表4　外方位元素解算值

為了驗證外方位元素的有效性，本文采用解算得到的外方位元素，對剩余的5個點分別根據(jù)像點坐標(biāo)，結(jié)合5米間隔的DEM數(shù)據(jù)，解算其對應(yīng)的平面坐標(biāo)（X,Y），并和通過攝站坐標(biāo)模擬的地面控制點坐標(biāo)（X’,Y’）進(jìn)行比較（高程值通過平面坐標(biāo)迭代或內(nèi)插得到），計算和對比情況如表5所列。

圖4　航攝影像分布示意圖

在表5的實驗驗證計算結(jié)果中，通過外方位元素反求的地面點坐標(biāo)和模擬的控制點坐標(biāo)之間的最大誤差為-9.266米，接近1︰10萬地形圖的量測精度，能夠滿足一般情況下采用傾斜影像對地面目標(biāo)進(jìn)行定位的精度需求。

表5　五個驗算點的計算數(shù)

4 結(jié)束語

當(dāng)前國內(nèi)POS數(shù)據(jù)中IMU精度不高，不能直接作為影像外方位角元素使用的現(xiàn)狀，與在應(yīng)急處置和跨境偵察等實際應(yīng)用中采用大傾斜影像進(jìn)行無控制定位的需求之間存在矛盾。本文針對此問題，吸收借鑒現(xiàn)有研究成果，對傳統(tǒng)的單鏡頭傾斜攝影方式進(jìn)行了改進(jìn)，采用下視和傾斜兩個相機(jī)組合攝影模式同時獲取影像，研究了在下視影像輔助下大傾斜影像的自動化處理技術(shù)與流程，實現(xiàn)了無控制條件下的大傾斜影像外方位角元素解算和地面坐標(biāo)量測。實驗情況表明該方法計算簡單，精度較好，能夠滿足無控制條件下傾斜影像對地觀測的需求?！?/p>

參考文獻(xiàn)

［1］ 李德仁，劉立坤，邵振峰.集成傾斜航空攝影測量和地面移動測量技術(shù)的城市環(huán)境監(jiān)測［J］。武漢大學(xué)學(xué)報·信息科學(xué)版，2015,40（4）, 427-435

［2］ Ad a m R.Mosbr ucker.Applicat ion of Dig it al Oblique Photogrammetry for Monitoring Erosion and Sedimentation at Mount ST.HELNS,Washington,USA.2nd Joint Federal Interagency Conference,Las Vegas,NV. June 27-July1,2010.

［3］ G.T.Grenzdorffer,M.Guretzk,I.Friedlander. Photogrammetric Image Acquisition and Image Aralysis of Oblique Imagery -Anew Challenge for the Digital Airborne System PFIFF.

［4］ K. Jacobsen.Geomet r y of ver t ical and oblique i mage combinations［OL］,2012.03,http：//spiedigitallibrary.org/proceedings/ resource/2/psisdg/5286/1/191_1.

［5］ MassIMilianoMolinari，Stefano Medda and samir Villani. Vertical Measurements in Oblique Aerial Imagery［J］.ISPRS Int.J.Geo-Inf.2014,3：914-928.

［6］ 江剛武，姜挺，王勇等.基于單位四元素的無初值依賴空間后方交會［J］.測繪學(xué)報，2007,31（2）：169-175

［7］ 官云蘭，程曉輝，周世健等.基于單位四元素的空間后方交會解算［J］.測繪學(xué)報，2008,37（1）：30-35

［8］ 張一平，龔志輝，王勃等.基于大傾角的角錐體空間后方交會研究［J］.測繪工程，2012,21（2）：17-20

［9］ 王建超，郭大海，鄭雄偉等.機(jī)載POS直接地理定位的精度分析［J］.國土資源遙感，2007,（4）：33-37

［10］ 郭大海，吳立新，王建超等.機(jī)載POS系統(tǒng)對地定位方法初探［J］.國土資源遙，2004,（2）： 26-30

［11］ Jean-Michel More, GuoShen Yu. ASIFT： A New Framework For Fully Affne Invariant Image Comparison. SIAM J.Imaging Sciences. 2009,2（2）：438-469.

［12］ 江延川.解析攝影測量學(xué)［M］.中國人民解放軍測繪學(xué)院，1991.5：37-43

Calculated Exterior Orientation Angle Elements of the Large Oblique Image with Vertical Images

Ke Ling1，2，F(xiàn)an Dongjuan2，Xu Bin3
（1.Basic science department of Wuchang Shouyi University, Wuhan, 430064； 2. School of Remote Sensing and information Engineering, Wuhan University, Wuhan, 430079； 3. 73608 PLA Troop, Nanjing, 210028）

Abstract:In this paper, based on the existing research results, using vertical and oblique photography to obtain multiple images, the traditional method of photography is improved. Under the visual image assistant, the calculation of the exterior orientation angle elements of the large oblique image is realized.

Keywords:large oblique image； visual image； POS； exterior orientation elements； without control points

doi：10.3969/J.ISSN.1672-7274.2016.06.003

中圖分類號：TP391

文獻(xiàn)標(biāo)識碼：A 文章編碼：1672-7274（2016）06-0007-04

基金項目：湖北省教育廳2015年度科研計劃指導(dǎo)性項目（B2015200）。

作者簡介：柯 玲，女，1980年生，湖北鄂州人，博士研究生，講師，主要研究方向是數(shù)字?jǐn)z影測量與遙感應(yīng)用、圖形圖像處理、數(shù)學(xué)應(yīng)用。