山區(qū)高速公路帶狀測區(qū)應用無人機航測技術(shù)研究

張寧　楊潤書　甘淑

摘要：

以云南省臨滄市某高速公路選線勘測工程為例，試驗利用低空無人機遙感技術(shù)對山區(qū)高速公路帶狀測區(qū)進行航空攝影，內(nèi)業(yè)使用軟件Inpho進行空三加密，制作出DOM、DLG等產(chǎn)品。對無人機航攝外業(yè)像控測量和內(nèi)業(yè)數(shù)據(jù)處理主要流程做了詳細介紹。通過高質(zhì)量完成測區(qū)1：1000數(shù)字地形圖產(chǎn)品生產(chǎn)并進行精度分析，證明利用低空無人機航測技術(shù)能夠為山區(qū)高速公路的勘測設(shè)計提供準確的地形基礎(chǔ)數(shù)據(jù)，可滿足山區(qū)帶狀測區(qū)1∶1000地形圖精度要求。相比于傳統(tǒng)地形測量，展現(xiàn)出航測在山區(qū)高速公路勘測方面的優(yōu)勢，提高了生產(chǎn)效率。

關(guān)鍵詞：

無人機航測技術(shù)；帶狀測區(qū)；Inpho；空三加密；數(shù)字正射影像圖；數(shù)字地形圖

DOIDOI：10.11907/rjdk.172944

中圖分類號：TP319

文獻標識碼：A文章編號文章編號：16727800（2018）003015504

英文摘要Abstract：Taking a highway survey project in Lincang of Yunnan as an example， it summaries the application of aerial photography using low altitude UAV remote sensing technology for mountainous highway belt area ， through MatchAT to make the result like DOM and DLG. In this paper， the main process of UAV aerial field control and internal data processing are introduced in detail. Finally， the production of 1： 1000 digital topographic map is completed by high quality and the precision analysis is carried out . It is proved that the low aerial survey technology can provide accurate terrain data for the survey and design of mountainous expressway， which can meet the requirement of 1∶1000 topographic map in mountainous belt area. Compared with traditional terrain measurement， aerial survey shows the advantages in the survey of mountain expressway， it improve the production efficiency.

英文關(guān)鍵詞Key Words：UAV aerial survey； belt area； Inpho； MatchAT； DOM； DLG

0引言

無人機航測技術(shù)是近年來發(fā)展起來的一門新興技術(shù)，以其靈活機動、操作維護簡單，可以便捷地獲取高分辨率影像的特點，成為當下測繪技術(shù)的熱門[1]。伴隨著計算機技術(shù)、攝影測量技術(shù)的發(fā)展，無人機航測技術(shù)開始運用到越來越多的專業(yè)領(lǐng)域，特別是對于地形條件復雜、高差起伏很大的地域，基本代替了傳統(tǒng)的人工測量[2]。無人機航測內(nèi)業(yè)理論不斷完善，出現(xiàn)很多無人機航片處理軟件，如svs、inpho等，使得航測精度大大提升。

云南省地處祖國西南邊陲，主要地形為山區(qū)，地形復雜，高差起伏相對較大。在高速公路勘測方面，由于測區(qū)有很長的工程區(qū)間，范圍十分狹窄且地形復雜，精度要求苛刻，運用傳統(tǒng)的工程測量方法雖能較好地滿足要求，但費時費力，效率低下。無人機低空航攝系統(tǒng)以無人機為飛行平臺，利用高分辨率相機系統(tǒng)獲取遙感影像，利用空中和地面控制系統(tǒng)實現(xiàn)影像的自動拍攝和獲取，同時實現(xiàn)航跡規(guī)劃和監(jiān)控、信息數(shù)據(jù)壓縮和自動傳輸、影像預處理等功能，是具有高智能化程度、穩(wěn)定可靠及較強作業(yè)能力的低空遙感系統(tǒng)[3]。隨著無人機航測技術(shù)的發(fā)展，其已在許多方面展現(xiàn)出優(yōu)勢，比如在農(nóng)村土地承包經(jīng)營權(quán)、土地整治等方面[4]。本文通過云南省某高速公路勘測的無人機航測生產(chǎn)實例，探討了無人機航測在山區(qū)高速公路勘測中的作業(yè)流程，為高速公路勘測選線提供了正射影像圖和1∶1 000的地形圖。

1測區(qū)概況

本項目測區(qū)位于云南省臨滄市轄區(qū)內(nèi)。測區(qū)呈帶狀，總長約90km，寬約500m，大部分測區(qū)位于山區(qū)，最低海拔約1 050m，最高海拔約1 850m，高差達800m，任務(wù)要求生產(chǎn)1∶1 000的正射影像圖和數(shù)字線畫圖，見圖1。

2工程實例

2.1實施流程

工程數(shù)據(jù)的獲取是內(nèi)業(yè)處理的前提，只有成像清晰、質(zhì)量良好的影像才可進行空中三角測量，滿足進度要求。無人機影像獲取流程以及后續(xù)處理流程見圖2。

2.2航線設(shè)計及航飛

本項目采用DB-2型無人機，航程約150km，最大抗風能力為6級，環(huán)境適應能力較強，無人機搭載尼康D810相機進行航拍，相機焦距為35mm，像素3630萬，像元大小4.8微米。航線設(shè)計將Google地球影像圖和SRTM World Elevation Data[5]疊加后進行分區(qū)設(shè)計，根據(jù)規(guī)范要求，像片航向重疊度一般為60%～80%，最小不應低于53%。旁向重疊度一般為15%～60%，最小不應低于8%。本次航線設(shè)計航向重疊度大于70%，旁向重疊度大于50%。航高根據(jù)公式（1）、（2）計算：

H=f×GSDa（1）

H0=H+h（2）

尼康D810焦距為35mm，分辨率為0.1m，像元大小為0.004 8mm，經(jīng)計算得出攝影航高H為729m。

部分航線分布如圖3所示。由于測區(qū)是一條狹長的帶狀，所以設(shè)計航線時必須結(jié)合飛行軌跡和測區(qū)位置綜合考慮，保證航測范圍略大于帶狀測區(qū)。

2.3外業(yè)像控點測量

像片控制點測量，為航測內(nèi)業(yè)成圖提供必要的像片平面坐標和高程控制，是保證空中三角測量數(shù)學精度的關(guān)鍵環(huán)節(jié)。外業(yè)控制點的布點方案采用區(qū)域網(wǎng)布點方法，參照規(guī)范GB/T 7391-2008 《1∶500 1∶1 000 1∶2 000地形圖航空攝影測量外業(yè)規(guī)范》[6]，像片控制點選點條件：

（1）在帶狀測區(qū)的條件下，像控點布設(shè)應滿足以下要求：①像控點位目標要求清晰，易于判別，選點以地面、較堅固、不易變化的目標為主；②像控點應均勻分布，且在旁向重疊范圍內(nèi)，布設(shè)的控制點盡量共用。一般在帶狀測區(qū)兩側(cè)200～250m布設(shè)，順著測區(qū)每400～650m布設(shè)一個點；③線狀地物交叉點（正交最佳），明顯地物拐角點；④選點位置應避開各類高大建筑物、高壓線、較大水面等，弧形地物及陰影不應選作點位目標。

（2）像片條件應滿足以下要求：點位距像片邊緣不應小于150像素。當目標條件與像片條件矛盾時應著重考慮目標條件。對于無明顯特征地物的區(qū)域，需要在航飛前進行人工布置十字標記，然后進行采集，可以把提前標記的像控點標在Google地圖上，航飛之后便可在影像上找到該點，然后刺點。

本項目還采用了GNSS輔助空中三角測量的方法，就是利用由差分GNSS相位觀測值進行相對動態(tài)定位所獲得的攝站坐標，作為區(qū)域網(wǎng)平差中的附加非攝影測量觀測值，以空中控制來減少地面控制的方法進行區(qū)域網(wǎng)平差[7]。此次像控點多選擇特征明顯色彩差異較大的地物，在相鄰測區(qū)也布設(shè)有公共控制點，以確保測區(qū)接邊的精度。帶狀測區(qū)部分像控點分布見圖4。

2.4內(nèi)業(yè)數(shù)據(jù)處理

在內(nèi)業(yè)完成空三加密前需對影像進行預處理，包括導入pos進行航帶編輯、畸變糾正、勻光勻色、影像旋轉(zhuǎn)、金字塔生成及圖像增強，然后進行空三加密、編輯DEM、正射糾正制作DOM、DLG。內(nèi)業(yè)采用專業(yè)的航空攝影測量軟件inpho進行空三加密，流程見圖5。

內(nèi)業(yè)進行空三加密（見圖6），即自動匹配圖像同名像點（連接點），然后進行區(qū)域網(wǎng)平差，這是整個流程中最關(guān)鍵的環(huán)節(jié)，匹配過程自動化；然后進行二次加密，在空三轉(zhuǎn)點的基礎(chǔ)上進一步加密轉(zhuǎn)點以得到更多的連接點；再進行空三編輯即人工交互編輯，主要操作包括：調(diào)節(jié)像片的對應關(guān)系、量測控制點、刪除連接點、增加連接點等。

空中三角測量依據(jù)CH/Z 3003-2010《低空數(shù)字航空攝影測量內(nèi)業(yè)規(guī)范》要求，具體精度要求滿足其規(guī)定，見圖7。

經(jīng)過空三編輯后進行DEM制作，數(shù)字高程模型是用一組有序數(shù)值陣列形式表示地面高程的一種實體地面模型[8]。本次采用空三加密輸出的結(jié)果進行密集匹配，然后根據(jù)密集匹配結(jié)果設(shè)置數(shù)字高程模型格網(wǎng)尺寸，用軟件自動生成DEM。對生成的DEM采用軟件VirtuoZo進行檢查，對DEM出現(xiàn)異常區(qū)域的部分進行修正。

本項目采用inpho的光束法區(qū)域網(wǎng)平差法進行空三加密，建立工程后匹配連接點，匹配策略采用特征匹配和最小二乘匹配穿插進行。由于之前進行過像片的畸變矯正，所以未采用軟件的自檢校參數(shù)。完成匹配連接點后sigma值是0.5，符合規(guī)范要求。

光束法區(qū)域網(wǎng)平差的基本思想是：以每張像片所組成的一束光束作為平差的基本單元，以共線方程作為平差的基礎(chǔ)方程。通過各個光束在空中的旋轉(zhuǎn)和平移，使模型之間公共點的光線實現(xiàn)最佳交會，并使整個區(qū)域納入到已知的控制點地面坐標系中去[9]。根據(jù)每張像片上控制點、待定點的像點坐標，按照共線方程條件列出誤差方程式，整體求解全區(qū)域內(nèi)每張像片的外方位元素以及所有加密點的地面坐標，主要內(nèi)容包括：

（1）獲取每張影像的外方位元素及連接點點坐標的近似值。

（2）由每張像片上的控制點和加密點的像點坐標，根據(jù)共線方程式列誤差方程式。

共線方程式是光束法平差的基礎(chǔ)數(shù)學模型[10]，見式（3）、式（4）。

x-x0=-fa1 （XA-XS）+b1（YA-YS ）+c1 （ZA-ZS）a3 （XA-XS ）+b3 （YA-YS ）+c3 （ZA-ZS） （3）

y-y0=-fa2 （XA-XS ）+b2 （YA-YS ）+c2 （ZA-ZS ）（a3 （XA-XS ）+b3 （YA-YS ）+c3 （ZA-ZS ） （4）

式（3）、式（4）中：x、y為像點的像平面坐標；x0、y0、f為影像的內(nèi)方位元素；XS、YS、ZS為攝站點的物方空間坐標；XA、YA、ZA為物方點的物方空間坐標；ai、bi、ci（i=1，2，3）為像片的3個外方位角元素組成的9個方向余弦。

完成匹配連接點后，開始量刺像控點。首先測區(qū)4個角刺4個控制點，然后平差一次，余下的像控點軟件就會預測到比較準確的位置，然后一一進行微調(diào)刺準即可，如圖8、圖9所示。

2.5DEM編輯

由于測區(qū)有植被覆蓋，房屋較多，所以在制作DOM前先對DEM進行編輯，使DEM切準地面，以DOM影像完整、美觀為準。在inpho中可以對生成的DEM進行手動編輯，利用inpho的點云內(nèi)插功能對濾波效果較差的區(qū)域進行內(nèi)插平滑處理。

2.6DOM制作

DOM制作過程：①選用色彩亮度較為合適的影像作為模板，對畸變處理后的影像進行整體勻光勻色；②利用修正后的DEM對勻光勻色的影像進行數(shù)字微分糾正，完成單幅DOM的制作；③對單幅DOM進行影像鑲嵌；④對鑲嵌后影像進行編輯，包括拼接線的編輯和影像拉花變形編輯。

拼接線應根據(jù)地物的情況選取，如拼接區(qū)域通過居民地時，拼接線可以繞過居民地，最終保證拼接線兩邊的地物色彩平衡，亮度反差適中，見圖8。

2.7DLG制作

采用“先內(nèi)后外”的作業(yè)方式，指先在室內(nèi)利用已有的影像圖對地物實行判讀和采集，然后把內(nèi)業(yè)初步生成的DLG成果疊加到DOM上制作調(diào)繪片，供外業(yè)到實地進行調(diào)繪，對內(nèi)業(yè)的采集結(jié)果進行校對和改正，最后再對外業(yè)成果進行整理編輯并最終成圖[11]，見圖9。

3數(shù)據(jù)精度分析

為保證最終成果的可靠性，在外業(yè)采集一系列的高程點作為精度的檢驗條件，采用GPS-RTK在測區(qū)均勻采集302個點。高程精度檢查點主要是道路交叉口、田坎交叉處、水溝溝底高程等，檢查點分布位置均勻并具有代表性，檢查位置點302個誤差為0.36m。可以看出，生產(chǎn)的地形圖精度滿足測繪規(guī)范對1∶1 000地形圖的精度要求，見表3。

誤差精度分析公式：

MΔH=∑ΔH2n-1（5）

式（5）中，m為點中誤差值，單位為米（m）；Δ為內(nèi)外業(yè)采集高程點殘差，單位為米（m）；n 為參與精度評定的點數(shù)。

4結(jié)語

云南省屬于祖國西南高原地區(qū)，地形條件復雜，屬于山地地形，常規(guī)測量十分艱苦。本文基于無人機航攝系統(tǒng)獲得影像數(shù)據(jù)和POS數(shù)據(jù)，探討了航測山地帶狀測區(qū)的關(guān)鍵技術(shù)。通過臨滄市某段高速公路工程的研究實踐，運用inpho軟件進行空三加密，按照規(guī)范生成1∶1 000的正射影像圖，利用其制作了地形圖。將實地量測檢查點坐標和基于模型生成的地形圖中檢查點的高程坐標進行對比，證實無人機航攝在高速公路勘測方面可以滿足實際要求。與傳統(tǒng)測量技術(shù)相比，大大減少了外業(yè)工作量，降低了生產(chǎn)成本，提高了作業(yè)效率。

參考文獻參考文獻：

[1]楊潤書，吳亞鵬，李加明，等.無人機航攝系統(tǒng)的特點及應用前景探討[J].地礦測繪，2011（1）：812.

[2]李德仁，李明.無人機遙感系統(tǒng)的研究進展與應用前景[J].武漢大學學報，2014（5）：505513.

[3]李永樹.基于無人機技術(shù)的地形圖測繪研究[J].測繪，2011，34（4）：147151.

[4]萬保峰，楊潤書，王亞男，等.基于無人機技術(shù)的山區(qū)1∶2000地形圖的成圖精度控制[J].地礦測繪，2015（3）：3449.

[5]楊堯蘭.天狼星無人機攝影測量系統(tǒng)在大比例尺成圖中的應用研究[D].昆明：昆明理工大學，2016.

[6]曹明蘭，薄志毅，李亞東.無控制點數(shù)據(jù)的無人機影像DOM快速制作[J].測繪學報，2016（8）：3538.

[7]劉葛.基于Matrix軟件系統(tǒng)正射影像圖的研究[D].昆明：昆明理工大學，2014.

[8]陳天祎.基于CIPS的低空無人機遙感影像處理研究[D].撫州：東華理工大學，2013.

[9]徐紅巖.基于DGPS技術(shù)的低空攝影測量應用研究[D].贛州：江西理工大學，2015.

[10]肖亮明，陳建忠.無人機航測技術(shù)在國土資源管理中的應用[J].測繪地理信息，2017，42（5）：9699.

[11]石平，鹿荻，申勇智.Inpho與Photomod軟件在無人機數(shù)字正射影像圖制作中的對比[J].城市勘測，2016（3）：7988.

責任編輯（責任編輯：杜能鋼）