非量測相機傾斜航空影像空三加密精度評價

馬晨，楊遼，池夢群，劉瑩，鄭宏偉

（1.中國科學院新疆生態與地理研究所，烏魯木齊830011；2.中國科學院大學，北京100049）

非量測相機傾斜航空影像空三加密精度評價

馬晨1，2，楊遼1，池夢群1，2，劉瑩1，2，鄭宏偉1

（1.中國科學院新疆生態與地理研究所，烏魯木齊830011；2.中國科學院大學，北京100049）

鑒于傾斜空三加密及其精度評價工作的重要性，使用非量測型相機獲取傾斜航空影像，利用計算機視覺的SIFT特征匹配算子對影像的空三加密進行改進，然后對加密結果進行精度評價。結果表明：對于人工構造物密集地區的立體影像，特征點的數量密集，空三加密精度傾斜影像檢查點的精度相較下視影像精度高；影像的分辨率為10cm，量測結果平面x、y的中誤差都在精度范圍中，高程z值中誤差略大；增加控制點數目后，平面x與y的精度明顯提高。分析影響試驗結果的誤差來源有：航攝質量、控制點布設與量測精度、內業人員操作、區域網平差方法等。

空三加密；傾斜航空影像；SIFT；計算機視覺；航空攝影測量；地面控制點

0 引 言

近年來，隨著對城市三維模型的需求，傳統的下視鏡頭可以獲得主體建筑的高度信息和頂部紋理信息，但在獲取建筑物的側面紋理上能力有限，難以滿足城市建筑物三維數據獲取的要求。而傾斜攝影可以實現建筑物的正直影像與其立面紋理的傾斜影像相結合，是城市景觀三維建模的主流趨勢［1］。傾斜航空攝影是在同一飛行平臺上搭載多臺或多種傳感器，同時從多個角度采集地面影像，從而克服了傳統航空攝影技術只能從垂直角度進行拍攝的局限性，能夠更加真實地反映地物的實際情況，彌補了正射影像的不足［2］。傾斜影像技術的應用，使目前高昂的三維城市建模成本大幅降低，大大提高三維城市建模的速度［3］。

空中三角測量又稱作空三加密，是三維重建的一個重要環節?？罩腥菧y量的數據質量決定著最終產品的數學精度，因此對空中三角測量數據的質量評定顯得尤為重要。近年來全數字自動化空中三角測量的技術日益成熟，作業平臺日益實用完善，現基本上進入到規?；纳a階段，但目前的文獻中對傾斜空三加密的精度還沒有一個確切評價標準，對傾斜的精度評價工作較少，然而空三加密的精度對最后模型的建成有很大的影響。所以對其數據質量的檢查驗收及質量評定已提到日程上來［4］。

1 基于計算機視覺的空三加密原理和精度評價方法

1.1 基于計算機視覺的空中三角測量

基于計算機視覺的影像空三加密，是通過計算機視覺技術的相關方法在影像群之間建立一種較為嚴格的相對關系，通過影像信息反推出拍攝影像時的真實場景，一個完整的立體視覺系統通?？煞譃閳D像獲取、攝像機定標、特征提取、立體匹配、深度確定及內插等6個部分［5］。通過影像信息反推出拍攝影像時的真實場景，也就是對于一組圖像數據中的兩幅影像，通過一種空間變換把一幅圖像映射到另一幅圖像，使兩幅圖像中在空間中同一位置的點一一對應。

（1）特征點檢測與匹配

首先檢測每張影像的特征點。SIFT（Scale Invariant Feature Transform）算子是計算機視覺領域非常著名的特征算子，可以很好地避免由于影像變形等原因造成特征點信息的丟失，可以用于模式識別和影像匹配，是一種基于尺度空間的、對圖像縮放、旋轉甚至仿射變換保持不變性的圖像局部特征描述算子。

（2）計算相對定向參數

使用運動模型算法進行相對定向，運動模型算法能夠從一組只有二維信息圖像中提取出三維信息，可以把它轉化為如下問題：有m幅圖像，n個空間中的點，有如下方程：

其中，Xij為第j個點在第i幅圖像中的二維信息，Pi為第i幅圖像的投影矩陣，Xj為第j個點的三維坐標，已知m×n個二維信息，要估算m個投影矩陣和n個點的三維空間坐標，根據特征點的匹配結果已知圖像特征點間的對應關系，所以可以利用運動模型算法求得這些特征點的空間坐標。

相對定向用7個參數來確定，其中6個外參，有

為像點的相空間輔助坐標［6］。

然后對運動模型相對定向的結果進行優化，加入附加參數包括相機畸變參數，曝光點pos數據，控制點數據。在影像匹配過程中，從金字塔各個層面，分別進行匹配，進行多次迭代。

（3）區域網平差

區域網平差采用光束法平差方法。由于光束法區域網平差是從原始的影像坐標觀測值出發建立平差數學模型的，所以只有在光束法平差中才能最佳的顧及和改正影像系統誤差的影響［6］。

空三加密的流程如圖1所示。

圖1 空三加密流程

1.2 空中三角測量的精度評價方法

空三加密的實際精度與理論精度是有一定差異的。測量中偶然誤差的影響與點位的分布有關，而實際情況是復雜的，往往要受到偶然誤差和殘余系統誤差的綜合影響。下面是用來衡量空中三角測量實際精度的估算公式：

其中，X控，Y控，Z控為控制點的外業測量坐標，X加，Y加，Z加為控制點的內業加密坐標，n為控制點的個數［7］。

2 試驗方案及試驗結果分析

2.1 試驗區概況

本次試驗位于新疆維吾爾自治區吐魯番市市區，試驗區地勢平坦，主要以平地為主，測區面積約為15km2，測區范圍：89°09′29.26″E～89°12′10.11″E，42°56′08.82″N～42°57′35.40″N試驗區概況圖如圖2所示。

圖2 試驗區概況圖

2.2 影像獲取

本文采用的是非量測型相機，為自制的傾斜航攝儀，是由5個鏡頭組成的傾斜航攝儀器（圖3）。傾斜航攝儀獲取兩種影像，一種是傳統的垂直向下的影像，另一種是四張傾角為40°到60°之間的傾斜影像，同時拍攝東南西北4個方向（圖4），并配合定位系統獲取高精度的地理位置和姿態信息。

量測型相機的內方位元素是已知的，其影像具有明確的幾何位置關系。但量測型相機價格昂貴，儀器笨重，進行外業拍攝時很不方便。非量測型數碼相機采用存儲卡存儲影像，可直接與計算機連接，與量測型相機相比，信息處理周期短，在地形復雜地區進行外業拍攝具有很強的機動靈活性［8］。但非量測型相機不能提供內方位元素，并且相機安裝、調試過程中難免會產生一些殘余像差，造成物鏡構象畸變，所以在飛行前需要對相機進行標定，計算出相機的鏡頭畸變參數和內方位元素值。采用DLT（直接線性變換方法）對非量測型相機進行檢校，計算出內外方位元素。5個相機的相機畸變參數（依次為相機焦距，徑向畸變參數，切向畸變參數，像主點偏移）如表1所示。

圖3 傾斜相機鏡頭示意圖

圖4 傾斜相機5個鏡頭影像示意圖

2.3 試驗方案

試驗采用的航空攝影相機焦距為35mm，影像幅面6000像素×4000像素，像元分辨率為5.966μm，航向重疊度為75%，旁向重疊度為70%，航線方向為東西走向，共22條航線，相鄰航線間距離為100m，每條航線73個曝光點。航攝相對航高250m，設計GSD（航攝像機CCD一個像元對應的地面尺寸大?。┫乱暈?cm，其余鏡頭為10cm。測區內布設4行，7列共計28個地面控制點。試驗區航線與控制點分布示意圖如圖5所示。

表1 相機鏡頭物理參數

圖5 試驗區航線與控制點分布示意圖

試驗一選取的區域范圍為圖中紅框標記區域范圍，使用左視鏡頭、下視鏡頭、右視鏡頭，共3個鏡頭的數據，總共1892張航空影像，12個控制點，其中6個控制點，6個檢查點。

試驗二用為圖中所有區域范圍，采用前視鏡頭、后視鏡頭、左視鏡頭、右視鏡頭、下視鏡頭，共5個鏡頭的數據，總共8030張航空影像，使用27個控制點，其中13控制點，14個檢查點。

2.4 精度分析

試驗采用SIFT特征點檢測對多視角影像進行檢測，在試驗中不同的地物場景所對應的各個視角影像SIFT特征點檢測的情況不同，對于人工構造物密集的地區的立體影像，特征點的數量密集，在地物比較稀疏的地區，特征點數量不足。在特征點檢測時，結構紋理特征不明顯的區域結果不好，不能保證特征點數量，需要額外的在這種區域上增加控制點，使得在區域網平差中，以權重數值的方式分配到各個匹配的同名點和控制點上，這樣就可以保證紋理特征不好的區域的匹配精度，優化整體的定向參數。城區的特征點匹配過程示意圖如圖6所示。

將得到的相機焦距，畸變參數，經過優化后的各角度影像和影像的6個外方位元素，進行立體量測。試驗一量測時分別在下視的影像對中和傾斜的影像對中進行量測。下視影像量測結果見表2，傾斜影像結果見表3。試驗二的量測結果見表4。

圖6 特征點匹配過程

表2 試驗一對下視影像對量測結果／m

表3 試驗一對傾斜影像對量測結果／m

從量測結果中可以看出在下視影像和傾斜影像中，平面x中誤差、平面y中誤差、高程z中誤差（m）都符合航空攝影測量的要求，在控制點的精度上，下視影像的精度比傾斜影像的精度高，而檢查點的精度則是傾斜影像的精度比較好。影像的分辨率為10cm，量測結果平面x，y的中誤差都在精度范圍中，高程z值中誤差略大。

表4 試驗二傾斜影像量測結果／m

在傾斜影像的量測結果中，由于控制點的數目增多，平面x與y的精度明顯比試驗一的精度高，控制點的精度比檢查點的精度好，高程z中誤差略大。

通過計算及精度統計，認為影響加密成果精度的因素有多種，其中主要因素為影像質量包括分辨率、清晰度、攝影比例。另外還包括飛行質量，像控點測量精度，外業像控的選點、刺點、描述準確情況；內業加密時對像控點點位的判定準確程度；內業加密人員利用光標對像控點的切準情況；內業加密點的選取、在不同像片上的判位；平差方法不同等［9］。

（1）加密精度受到航攝質量的影響

航攝影像是空三加密最原始的資料和基礎，飛行的質量影響著加密的成果，如：影像的質量（包括影像的灰度、色調、分辨率）、影像的重疊度、飛行航高的穩定度、影像的傾角、旋偏角、航線的質量，都是保證空三加密順利進行的重要因素。

（2）加密精度受到控制點布設與量測精度的影響

控制點在外業對點的選取，刺點以及對點的描述情況上也是影響加密成果的原因之一，控制點選取在影像清晰，特征明顯容易判讀的地物上，內業在刺點過程中的準確度會提高，反之會給加密結果帶來一定的誤差，高精度的控制點可以有效防止加密誤差的傳遞和積累。此外，控制點數目與分布也會造成誤差，區域網中每個位置上的像控點在平差中所起的作用是不同的，如：4個角點的像控點主要是歸劃比例尺和改正模型變形，航線中間的像控點主要是改正模型彎曲等，所以區域網的大小、形狀以及像控點的布設方案及其在像片上的位置等因素對空三成果精度有著舉足輕重的作用［10－11］。

（3）精度受到內業人員操作的影響

內業加密人員對控制點位置判讀的準確程度對加密精度也會造成影響，不同影像上對同一控制點判讀位置的誤差，內業加密人員的工作經驗、對影像的判讀能力，立體感的好壞，都會直接對加密結果有影響，這種誤差是不可避免的。

（4）精度受到區域網平差方法的影響

實際生產中常用平差軟件所遵循的理論依據相同且平差方法均采用光束法區域網平差，但各軟件在同名點匹配、平差計算過程中相關參數權值設置方法、區域網粗差剔除算法等方面也不盡相同。因此，針對初始條件完全相同的區域網平差工作，采用不同平差軟件，其計算結果會有所差異［10－11］。

3 結束語

空中三角測量是航空攝影測量的重要步驟，它的精度結果直接影響著后期數字高程模型、數字地面模型、數字正射影像、三維建模的質量結果。以多次的試驗和生產經驗對傾斜航空影像的空三精度進行評價，對造成精度結果誤差的可能來源進行分析，主要有航攝質量、控制點的數目及布設、內業人員的操作、區域網平差的方法和軟件等。找到問題的關鍵點，避免操作過程中的誤差積累，才能保證空中三角測量的質量和精度，提高工作效率。

［1］桂德竹，林宗堅，張成成.傾斜航空影像的城市建筑物三維模型構建研究［J］.測繪科學，2012，37（4）：140－142.

［2］李祎峰，宮晉平，楊新海，等.機載傾斜攝影數據在三維建模及單斜片測量中的應用［J］.遙感信息，2013，28（3）：102－106.

［3］李鎮洲，張學之.基于傾斜攝影測量技術快速建立城市3維模型研究［J］.測繪與空間地理信息，2012（4）：117－119.

［4］許婭婭，王月莉.全數字空中三角測量數據檢查驗收及質量評價方法的探討［J］.測繪通報，2006（1）：33－35.

［5］張祖勛.數字攝影測量與計算機視覺［J］.武漢大學學報（信息科學版），2005，29（12）：1035－1039.

［6］張劍清，潘勵，王樹根.攝影測量學［M］.武漢：武漢大學出版社，2003.

［7］盧紅艷.淺談空中三角測量的精度［J］.中國地名，2010（7）：74.

［8］楊軍.非量測數碼相機檢校方法的研究［J］.測繪科學，2009（s1）：89－91.

［9］朱肇光，孫護，崔炳光.攝影測量學［M］.武漢：武漢大學出版社，1994.

［10］吳正鵬，張曉東，奚歌.淺談全數字空中三角測量精度影響因素及質量控制［J］.測繪與空間地理信息，2013，36（11）：184－186.

［11］高文革，范永杰，宋倩，等.空三加密在數字攝影測量中的精度分析［J］.山東國土資源，2010，26（12）：26－29.

Accuracy Assessment of Oblique Aerial Images and Aerial Triangulation Measurement Based on Non－metric Digital Camera

MA Chen1，2，YANG Liao1，CHI Meng－qun1，2，LIU Ying1，2，ZHENG Hong－wei1
（1.Xinjiang Institute of Ecology and Geography，Chinese Academy of Sciences，Urumqi 830011；2.University of Chinese Academy of Sciences，Beijing100049）

Aerial triangulation is a crucial step in three－dimensional（3D）modeling.The traditional photogrammetry cannot satisfy the requirement of 3Dmodeling in current commercial applications.However，recently，the rise of oblique photogrammetry has a good complement to the traditional aerial photogrammetry.In this paper，we proposed a novel approach based on the utilization of oblique aerial images taken by agroup of non－metric digital cameras.We improved the process of aerial triangulation via a SIFT matching computer vision approach in the first step.Secondly，we analyzed the accuracy of aerial triangulation as well as the sources of several errors which affect the result of the experiment.We found that feature points are dense when artificial structure areas on stereo images are few.The accuracy of the aerial triangulation and oblique aerial images is better than that of aerial images.Resolution of the image is 10cm，error in plane x and y are in the range of precision，the elevation error is slightly larger in the z value；the precision of plane x and y has improved significantly after increase the number of control points.The sources of errors which affect the result of the experiment are quality of aerophotography，control points setting and accuracy in measurement，personnel operating，method of block adjustment，etc.

aerial triangulation；oblique aerial image；SIFT；computer vision；aerial photogrammetry；ground control point

10.3969／j.issn.1000－3177.2015.06.014

P231

A

1000－3177（2015）142－0071－05

2014－10－20

2014－11－25

千人項目（Y474161）；國家863項目（2013AA122302）。

馬晨（1990—），女，碩士，研究方向為攝影測量與遙感。

E－mail：sensen810＠163.com