趙一方

（中鐵第一勘察設(shè)計(jì)院集團(tuán)有限公司，西安 710043）

無人機(jī)航空攝影測量具有機(jī)動(dòng)靈活、影像清晰和成本低廉等優(yōu)勢，可以便捷地獲取地物地貌信息，彌補(bǔ)了既有手段效率低、成本高、周期長的短板，已成為傳統(tǒng)攝影測量技術(shù)的有效補(bǔ)充[1-3]。但無人機(jī)攝影測量使用的立體模型仍是單個(gè)像對，觀測范圍有限，需要頻繁切換像對來進(jìn)行大范圍場景瀏覽，嚴(yán)重影響生產(chǎn)效率[4]。

為充分發(fā)揮無人機(jī)航測技術(shù)優(yōu)勢，進(jìn)一步提高鐵路無人機(jī)航測制圖生產(chǎn)效率，基于李德仁等[4-5]提出的航空影像大場景立體模型理論，本院自主研發(fā)了基于無人機(jī)影像大場景立體模型的數(shù)字化測圖軟件[6-8]。無人機(jī)大場景能夠大范圍地再現(xiàn)地形地物三維信息，彌補(bǔ)了單像對立體模型觀測范圍受限的缺點(diǎn)，在大場景模型上進(jìn)行測圖能夠極大地提高無人機(jī)航測制圖效率。

1 大場景立體模型

1.1 基本原理

大場景立體模型是由數(shù)字微分糾正鑲嵌后的大范圍無縫正射影像（左片）與引入左右視差后生成的立體輔助影像（右片）構(gòu)成的，可以看作是一個(gè)大范圍的單像對立體模型，其核心技術(shù)包括數(shù)字微分糾正鑲嵌與人為視差引入。數(shù)字微分糾正鑲嵌的目的是把大量離散的航攝像片制作成為大范圍無縫正射影像，但由于消除了地形投影差，大場景立體模型的地形理論上是一個(gè)平面，沒有起伏，僅地物保留了立體效果。因此需要人為在大場景右片上引入左右視差，來模擬消除了的地形投影差，達(dá)到恢復(fù)地形起伏的目的。

1.2 制作流程

無人機(jī)影像大場景立體模型的制作流程主要包括影像勻光勻色、DEM 生成、拼接線編輯和大場景立體模型生成4 個(gè)步驟，如圖1 所示。

圖1 無人機(jī)影像大場景立體模型制作流程圖

影像勻光勻色是指將輸入影像調(diào)整成為亮度適中、反差均勻、色彩一致的影像。由于航空攝影過程中時(shí)相、環(huán)境、姿態(tài)和對象等客觀因素的差異，不同測段、不同航線之間的影像容易出現(xiàn)過亮過暗、色彩不一致、噪聲過大和反差不均勻等現(xiàn)象。使用這樣的影像進(jìn)行數(shù)字鑲嵌會(huì)導(dǎo)致大場景立體模型接邊痕跡明顯，視覺效果差，給后續(xù)工作如影像判讀、數(shù)字化測圖等帶來不便[9]。為了保證無人機(jī)影像大場景立體模型的色彩一致性，對于色彩不均衡、反差不均勻的原始無人機(jī)影像，可以選擇進(jìn)行勻光勻色處理。

測區(qū)DEM 生成首先依靠嚴(yán)密的空三加密成果匹配和解算出海量加密點(diǎn)，其次通過點(diǎn)云分類技術(shù)篩選出位于地面的加密點(diǎn)，最后對地面點(diǎn)進(jìn)行重采樣生成測區(qū)DEM。由于點(diǎn)云自動(dòng)分類出的地面點(diǎn)含有錯(cuò)點(diǎn)，測區(qū)DEM 不可避免地含有誤差，但實(shí)際上對大場景立體模型量測精度影響較小[8]。DEM 的分辨率宜為無人機(jī)大場景模型的50 倍，如制作0.1 m 分辨率的大場景模型，DEM 格網(wǎng)間距為5 m 較為合適[7]。

由于正射影像無法消除地物投影差，如果在投影差大的地方（如高層建筑）進(jìn)行影像鑲嵌，那地物的一部分就會(huì)被裁切掉。通過對鑲嵌線進(jìn)行人工編輯，使之盡量避開建筑物，讓影像在沒有地物的區(qū)域進(jìn)行拼接，可以減少因鑲嵌造成的影像重要信息丟失。

2 項(xiàng)目實(shí)例

2.1 測區(qū)概況

項(xiàng)目位于山地丘陵地區(qū)，屬構(gòu)造侵蝕溶蝕中低山地貌，本文研究區(qū)域面積約12 km2，測區(qū)內(nèi)相對高差最大約400 m。

2.2 數(shù)據(jù)獲取

2.2.1 航空攝影

采用大白II 無人機(jī)搭載飛思IXU-RS1000 數(shù)碼航攝儀進(jìn)行測區(qū)范圍的航攝數(shù)據(jù)獲取，主要相機(jī)參數(shù)見表1。該系統(tǒng)配備GNSS RTK、減震云臺(tái)及高精度IMU，實(shí)現(xiàn)直接傳感器定向或輔助定向測量，可準(zhǔn)確獲取航攝儀曝光時(shí)刻的外方位元素。

表1 飛思IXU-RS1000 數(shù)碼航攝儀主要參數(shù)

2.2.2 像控選刺

像控點(diǎn)布設(shè)采用區(qū)域網(wǎng)布點(diǎn)，像控點(diǎn)的密度和位置依據(jù)CH/T 3003—2021《低空數(shù)字航空攝影測量外業(yè)規(guī)范》要求確定。以正射影像和像控點(diǎn)布設(shè)結(jié)合圖為參考，采用網(wǎng)絡(luò)RTK 定位技術(shù)完成像控點(diǎn)選刺與測量。

2.3 數(shù)據(jù)處理

2.3.1 空三加密

空三加密作業(yè)采用GodWork 軟件，依次進(jìn)行連接點(diǎn)自動(dòng)提取、自由網(wǎng)平差與漸次刪點(diǎn)、像控點(diǎn)量測和光束法區(qū)域網(wǎng)平差等工作。通過漸次刪點(diǎn)，將內(nèi)業(yè)加密點(diǎn)重投影誤差控制在1.5 個(gè)像素以內(nèi)，保證模型連接質(zhì)量。本項(xiàng)目像控點(diǎn)多為房角等具有一定高度的特征點(diǎn)，因此采用雙像立體量測，以減小控制點(diǎn)量測誤差。不斷優(yōu)化像控點(diǎn)量測位置，直至區(qū)域網(wǎng)平差后像控點(diǎn)平面和高程殘差滿足相應(yīng)規(guī)范為止[10]，最終像控點(diǎn)平面和高程殘差分別為0.084、0.070 m?？杖用芡瓿珊筝敵鯥npho 和Smart3D 格式成果，分別供DEM 和大場景制作使用。

2.3.2 影像勻光勻色

影像勻光勻色采用GeoDodging 軟件的批處理功能，處理模式選擇“整體自適應(yīng)”，處理方法為“濾波法”，對航攝像片進(jìn)行批量處理。利用原始影像與勻光勻色后的影像分別制作大場景立體模型，對比處理前后大場景視覺效果如圖2 所示。

圖2 影像勻光勻色前后無人機(jī)大場景（左片）視覺效果對比

2.3.3 生成測區(qū)DEM

測區(qū)DEM 采用Inpho 軟件的Match-T DSM 模塊制作，主要考慮到該模塊具有DEM 平滑功能，使用平滑后的DEM 進(jìn)行數(shù)字微分糾正可以避免正射影像出現(xiàn)扭曲拉花的現(xiàn)象，優(yōu)化大場景模型的視覺效果。在Match-T DSM 中定義格網(wǎng)大小為10 m，設(shè)置平滑選項(xiàng)為中度，運(yùn)行后生成成果為GRD 格式，需要轉(zhuǎn)換為TIF 格式，供大場景制作使用。

2.3.4 編輯鑲嵌線

在ArcGIS 中對鑲嵌線進(jìn)行編輯，以原始鑲嵌線制作的大場景左片為底圖，移動(dòng)頂點(diǎn)位置使鑲嵌線繞避房屋等拼接有問題的區(qū)域，讓影像在草木區(qū)進(jìn)行拼接，利用編輯過后的鑲嵌線重新制作大場景，恢復(fù)圖面建筑完整信息。鑲嵌線編輯前后無人機(jī)大場景模型拼接效果如圖3 所示。

圖3 鑲嵌線編輯前后無人機(jī)大場景（左片）局部拼接效果對比

2.3.5 無人機(jī)大場景模型制作

基于前述的Smart3D 格式空三加密成果、勻光勻色影像、測區(qū)平滑DEM 和編輯后的鑲嵌線，分別在中鐵第一勘察設(shè)計(jì)院（以下簡稱“我院”）自主研發(fā)軟件SAT_UAV_APP 和ImageStereo 中進(jìn)行無人機(jī)影像大場景立體模型制作與地貌地物數(shù)據(jù)采集，生成滿足項(xiàng)目技術(shù)質(zhì)量要求的1∶2000 比例尺地形圖，如圖4、圖5 所示。

圖4 我院自主研發(fā)的大場景生成軟件SAT_UAV_APP

圖5 我院自主研發(fā)的大場景立體測圖軟件ImageStereo

2.4 精度分析

2.4.1 統(tǒng)計(jì)方法

完成本項(xiàng)目地形圖測繪工作后，對地形圖成圖精度進(jìn)行驗(yàn)證，以測區(qū)內(nèi)像控點(diǎn)RTK 測量坐標(biāo)值為參考，統(tǒng)計(jì)地形圖成果中對應(yīng)點(diǎn)位相對于參考值的偏差，并計(jì)算平面位置中誤差和高程中誤差，驗(yàn)證基于無人機(jī)大場景技術(shù)測繪鐵路地形圖的數(shù)學(xué)精度，相應(yīng)計(jì)算公式如下

式中：x、y、h 分別為東、北、高坐標(biāo)值；mxy和mh分別為像控點(diǎn)平面位置中誤差和高程中誤差；Δ 為像控點(diǎn)各方向上大場景測量值與RTK 測量值之間的偏差；n 為測區(qū)像控點(diǎn)總數(shù)。

2.4.2 統(tǒng)計(jì)結(jié)果

在無人機(jī)影像大場景立體模型上測量所有像控點(diǎn)坐標(biāo)，與RTK 測量成果進(jìn)行比較，驗(yàn)證無人機(jī)大場景模型測繪鐵路地形圖的平面和高程數(shù)學(xué)精度，結(jié)果如圖6 所示。

圖6 無人機(jī)大場景模型測繪鐵路地形圖平面（左）和高程（右）殘差分布圖

統(tǒng)計(jì)可得，以RTK 測量坐標(biāo)值為參考，研究區(qū)域內(nèi)像控點(diǎn)的平面位置中誤差為0.109 m，高程中誤差為0.255 m，滿足國家標(biāo)準(zhǔn)中對于二級(jí)地形、1∶2000 比例尺地形圖地物平面中誤差0.6 m、高程中誤差0.5 m 的要求[11]。

3 結(jié)論

無人機(jī)大場景模型利用數(shù)字微分糾正鑲嵌技術(shù)彌補(bǔ)了單個(gè)無人機(jī)立體像對可視范圍窄小的缺點(diǎn)，極大地提高了三維場景的瀏覽和量測效率。本文借助某山地丘陵鐵路地形圖測繪項(xiàng)目，對無人機(jī)大場景模型的數(shù)學(xué)精度進(jìn)行驗(yàn)證，得出其滿足二級(jí)地形、1∶2000 比例尺地形圖國家標(biāo)準(zhǔn)精度要求的結(jié)論，為鐵路行業(yè)大比例尺地形圖測繪提供了新的解決方案。