賈秋英
(中煤航測遙感集團有限公司,陜西 西安 710000)
傳統的地籍測繪采用全站儀、GPS-RTK 等設備進行數據采集,效率低、費時費力、成本高,已不能滿足當今社會對其現勢性的需求[1]。無人機航空攝影的出現,改變了大比例尺地形圖測繪的方式,但是其精度達不到地籍測繪精度需求。傾斜攝影可以獲取建構筑物多個角度的紋理信息,在空三加密環節,增加了更多的約束條件,使其精度有了很大的提升[2]。本文從傾斜攝影技術入手,對傾斜攝影測量技術、傾斜攝影系統的組成、傾斜攝影建模流程進行簡單說明,并以實際生產項目為例,對外業數據的獲取和內業數據處理流程進行了詳細介紹,并對生產的地籍圖精度進行了檢測,結果表明,本文的方法生產的地籍圖滿足地籍二級精度需求。
無人機傾斜攝影測量是指在低空無人機上,掛載多鏡頭航攝儀,從多個角度對地面進行拍攝,獲取盲區少、重疊度高的多角度影像[3]。一般掛載的都是非量測數碼相機,獲取影像存在畸變大的問題,一般在數據處理時,都需要對影像進行畸變糾正。在數據處理時,依據攝影關系和像控點,恢復影像地物上彼此之間的關系,用于測繪產品的生產以及測量。
無人機傾斜攝影系統組成主要有飛行平臺、傾斜航攝儀、電腦硬件和數據解算軟件。飛行平臺主要是指能夠低空飛行,能夠獲得超高分辨率的原始影像數據。傾斜航攝儀是指可以掛在多鏡頭航攝儀,完成對地面影像全方位、多角度、少盲區的影像采集。電腦硬件主要是指建模對電腦硬件要求高,需要大內存、高主頻、圖形渲染超強的顯卡。數據解算軟件是指傾斜攝影建模軟件,常見的有Photoscan、PIX4D、Smart3D2019、Context Capture 等,本次實驗使用Context Capture 進行。
利用傾斜攝影建模軟件,將航拍影像轉為實景三維模型的過程。在軟件運行過程中,主要涉及空中三角測量、多視影像密集匹配、不規則三角網優化與構建等。在整個作業過程中,具體包括飛前準備、外業數據獲取、內業數據處理、精度檢測等,具體流程見圖1。

圖1 無人機傾斜攝影建模流程圖
本次為某農村地籍圖測繪項目,測區房屋較多且密集,傳統作業方式難度大且效率低。采用無人機傾斜攝影技術完成測區地籍圖項目的生產。結合地形和地籍圖精度要求,選用多旋翼無人機,搭載5 鏡頭相機,地面采樣分辨率設置為0.015 米,航向、旁向重疊度均設置為85%,航高80 米,測區非禁飛區,屬于空域開放區域。選擇某一晴天,上午十點到下午兩點之間進行航飛,航飛2 架次,獲得有效影像數12250 張。測區面積約1km2,按照200 米間距均勻布設并采集像控點30 個,采集檢測點15 個。
本次實驗內業用到的軟件主要有Bentley 公司的Context Capture 實景三維模型生產軟件,北京清華三維的EPS地籍測繪軟件。
2.2.1 Context Capture 軟件。是目前實景三維模型生產的主流軟件,具有高效率、可集群作業、生產的模型質量高、精度高等特點,是目前建模軟件中用戶群體最多的一款軟件。該軟件版本更新速度快,各模塊功能較完善。
2.2.2 EPS軟件。是一款綜合地理信息工作站軟件,具有模塊較多,基本上涵蓋地理信息各個行業的應用,可根據客戶需求定制開發,本次主要使用其三維測圖模塊完成地籍圖的采集與編輯。
外業數據獲取主要包括影像數據采集前的測區勘查、空域申請、已有資料收集、航線規劃、影像數據采集和像控點采集。
2.3.1 測區勘查。主要是對任務區的地形進行勘查,對當地的人情風俗進行了解,便于后期工作的展開。
2.3.2 空域申請。如果測區是禁飛區或者空域監管區,需要按照空域申請流程進行航飛攝影空域申請,在未獲得批準前,不得違規飛行,以免帶來某些意外。在獲得空域批準后,需根據批準文件里要求的時間、地點和范圍進行航飛作業。本次測區為空域開放區,因此在航飛時未進行空域申請。
2.3.3 已有資料收集。主要收集測區已有的測繪產品成果,本次收集到了測區已有的數字正射影像圖,可作為工作底圖;收集到周邊已有的高精度埋石點,可作為GPS-RTK精度檢核的已知點。
2.3.4 航線規劃。利用地面站軟件設置航向、旁向重疊度為85%,地面采樣分辨率為0.015 米,提交任務,無人機按照設置好的參數進行航攝作業。
2.3.5 影像數據采集。主要是按照規劃好的航線執行,但是需要在地面站上時刻觀察飛機飛行狀態和影像采集狀態,避免發生數據丟失、偏離航線等問題。在正式起飛采集影像前,需在地面對相機拍照進行檢測,確保相機正常拍照,內存卡正常存儲影像數據。
2.3.6 像控點采集。首先按照200 米間距,結合測區形狀均勻布設30 個像控點,然后利用紅色油漆對像控點點位進行噴涂,這里需要注意的是,影像數據采集前就需要完成像控點點位的噴涂。利用GPS-RTK對已知的埋石點進行采集,然后和已有準確的坐標進行比對,誤差在2 公分內即可。在檢核儀器后,對噴涂好的像控點點位進行坐標采集,每個點位至少采集三次,三次較差在2 公分之內,需要等到固定解狀態下進行采集,在采集的時候,通過多個角度對像控點點位進行拍照。采集完成后,對像控點進行平均值計算,對采集的影像按像控點進行整理。
2.4.1 數據預處理。只要是刪除2 架次的無效影像和POS,整理POS數據,刪除無用的參數,只保留相機曝光時的經緯度坐標值和高程值。無人機由于姿態較差,引入姿態可能對后續空三解算帶來不便,因此這里對姿態角也進行刪除,利用軟件空中三角測量解算,來計算每張影像的空間位置和姿態。利用EPT軟件,對原始影像進行勻色處理,提高影像的質量。利用重命名軟件對2 架次影像進行重命名,確保2 架次影像無重名。更改POS的點號名,使之與影像一一對應,確保可一次性引入POS數據。
2.4.2 空中三角測量。由于傾斜攝影依據的也是共線方程,但是共線方程是以垂直地面攝影理論來計算的,加上側視相機本來存在影像變形嚴重的問題,空三算法并不能對傾斜數據進行有效解算。為了提升傾斜影像空中三角測量解算的成功率,首先利用少量影像進行空中三角測量,得到5 個相機精確地相機參數,然后將得到的參數應用到所有影像解算中來,這樣可以提高空中三角測量的解算成功率。本次實驗,利用5 組鏡頭共240 張影像對5 個鏡頭相機參數進行解算,得到5 個鏡頭精確地相機參數,然后進行整體空中三角測量。空三計算完成后,通過人機交互的方式查看空三成果無分層,無彎曲等問題。空三報告中顯示加密點中誤差為0.39 個像素,加密點重投影中誤差為0.46 個像素,成果精度符合要求,自由網空三解算合格。
設置坐標系,導入30 個像控點,對30 個像控點全部進行轉刺。本次使用的Context Capture 軟件版本較新,支持在轉刺過程中實時計算像控點位置。通過對一個點轉刺3 張照片,可快速計算剩余照片的精準點位,有效提升了像控點轉刺的效率。在完成像控點的轉刺后,再次提交空三任務,連接點選擇保持,完成帶像控點的平差調整。通過平差后的空三報告,查看像控點的精度,平面位置中誤差為0.008 米,高程中誤差為0.006 米,像控點精度符合項目需求規范。
2.4.3 三維模型生產。實驗電腦內存為64G,為了提高建模效率且保證不出現內存不夠的情況,這里設置建模使用內存為電腦內存的一半,按照平面格網劃分瓦片。瓦片大小為100 米,預計使用內存31G,模型輸出分辨率設置為0.015 米,設置輸出模型格式為OSGB,提交建模任務,計算機自動完成實景三維模型的生產。在三維模型生產完成后,提交正射影像輸出任務,得到數字真正射影像產品。
2.4.4 地籍圖采集。在清華三維EPS軟件中,在三維測圖模塊下,對OSGB格式的模型進行轉換并加載,然后加載正射影像。這里主要以實景三維模型和數字真正射影像為參考,進行地籍圖的采集。在采集無屋檐的房屋時,利用真正射影像直接采集;在采集需做屋檐改正的房屋時,利用實景三維模型,在建構筑物的主體結構上進行采集,直接對屋檐進行了改正。在采集道路上時,可利用正射影像為底圖,直接在正射影像上進行道路采集。對門墩、圍墻等按照地形圖要求進行采集。在采集完成后,對由于模型問題內業未采集的,可通過外業進行補充采集,確保最后地籍圖成果完整可用。
對15 個檢測點,按照同精度中誤差計算公式進行計算,得到15 個點的平面位置中誤差為0.053 米,高程中誤差為0.046 米,地籍圖規范要求一級地籍精度5 公分,二級地籍精度10 公分,三級地籍精度15 公分,本次實驗精度符合地籍二級精度要求,可作為地籍測繪項目生產的一種作業方式。15 個檢測點的具體平面位置較差和高程位置較差見表1。

表1 檢測點位置較差統計表
本文首先對傾斜攝影測量做一簡單介紹,并以實際生產項目為例,通過對外業數據獲取,內業數據處理內容進行說明。并使用GPS-RTK采集的像控點對采集的地籍圖精度進行檢測,并按照同精度中誤差計算檢測點的中誤差。檢測結果表明,按照本文生產的地籍圖,精度可以滿足地籍精度需求,為地籍圖生產的作業人員提供了一種新的作業思路,提升了地籍圖測繪的效率,有助于確保地籍圖生產按時按點完成。