無人機攝影測量爲您介紹無人機航空攝影測量在土石方量計算中的應用

0　引　言

近年來，伴随國家經濟高速發展，各種基礎設施建設項目越來越多。土石方量計算是工程項目核心環節之一。爲瞭(le)能合理安排項目進度，準確(què)計算工程量大小與費用，通常需要高效、準確(què)地計算土石方量。因此選擇合适的測繪方法十分重要。

傳統的土石方測量方法有水準儀測量法、全站儀測量法和GPS測量法。水準儀測量法是通過使用水準儀測量事先在測區布設方格網的每個角點高程來計算土石方量的。該方法适用性單一，若測區不适合布設方格網，該方法就不适用瞭，且費時費力。全站儀測量法具有操作簡單，儀器要求低等優點，适合測量面積較小和通視良好的區域，反之，則會非常繁瑣，且效率低下
。GPS測量法是目前土石方測量中應用較多的一種方法，它不受距離和通視限制，且測量速度和精度較全站儀測量有所提高，但當測區有一些建築、樹木、電磁場等影響GPS信号時，該方法就不太适用瞭。因此傳統方法受場地影響大、效率低下、人工成本高，亟待尋求一種高效、安全且經濟的測量方法。

新興無人機航測技術爲解決上述難題開辟瞭(le)一條嶄新途徑。無人機航測作爲測繪發展的新技術，以其機動靈活、數據現勢性強、影像分辨率高、減輕勞動強度、提高生産效率等優點，已在工程勘測、設計、施工、竣工驗收及運行等多個環節中發揮瞭(le)重要作用。國内外衆多學者也已開始嘗試用該技術進行土石方量測量。該方法不受場地障礙影響，費用相對低廉，在對場地土石方量追蹤管理方面成本較低，同時由於(yú)避免瞭(le)大量人工現場作業，大幅提高瞭(le)測量人員的安全保障。

本文採(cǎi)用無人機航測技術，對測區改造前後的地形分别進行航空攝影測量，獲取場地變(biàn)化前後三維地形及影像數據，再将數據導入Autodesk 的Civil 3D軟件中，快速生成前後兩次三維地形模型，以對改造前後土石方量的計算進行分析。

1　實施技術路線

土石方量計算的目标是求取地表面體的體積差，關鍵在於(yú)對現狀地形和改造後地形進行表述
。因此，利用無人機航測技術進行土石方量計算主要是通過無人機獲取地形改造前後同一區域數字高程模型DEM 與數字正射影像DOM，然後将地形改造前後的DEM 與DOM 導入Autodesk的三維可視化分析軟件Civil  3D中，分别進行統計分析，再通過DOM 選取需要計算土石方量同一區域的DEM，對其進行統計分析，得出改造過程中土石方的填挖方量。總之，基於(yú)無人機航空攝影測量的土石方量計算作業流程主要包括：測區踏勘、航線規劃、地形改造前後航測數據採(cǎi)集、航測數據處理、DEM 制作、DOM 制作、地形改造前後數據分析、土石方量計算分析8項。

2　應用實例

本實例爲某一房建項目的土石方量計算任務。本任務航飛工作採(cǎi)用六旋翼微型無人機對改造前後的地形分别進行航測數據採(cǎi)集，内業航測數據處理使用Pix4D 進行，再使用Autodesk的Civil  3D軟件做前後兩次地形分析，無人機攝影測量計算改造前後地形變(biàn)化，得出場地平整過程中的土石方量。

2.1　外業數據採集

實施航測(cè)外業數據採(cǎi)集時，基本流程如下。

(1)航帶設計。對測(cè)區進行現場踏勘，確(què)定測(cè)區範圍與地形情況[6] ，根據《低空數字航空攝影測(cè)量外業規範》(CH/Z　3004—2010)中1∶500地形圖航測(cè)規範要求，對本項目外業進行航帶設計。根據要求，該項目確(què)定航飛高度爲200m，航向重疊度和旁向重疊度分别爲75%和45%。

(2)像控點布設。根據規範要求，本項目在測區四周以及中間均勻布設瞭(le)12個平高控制點，保證100張相片能有6個控制點，並(bìng)保證控制點至少能在兩張影像上同時找到。

(3)航測參數設定。根據航飛現場情況，確(què)定各項飛行參數
，完成航測數據採(cǎi)集工作。

2.2　航測數據處理

對採(cǎi)集完成的兩期航測(cè)數據採(cǎi)用Pix4Dmapper軟件進行數據處理。

2.3　土石方量計(jì)算與(yǔ)分析

在計算機處理過程中
，地表面模型DEM 的主要數字表現形式是不規則三角網TIN，故土石方量計算主要是利用TIN進行的。因此，需要将改造前地表DEM 生成的TIN 和改造後生成的TIN 疊加，形成交線，即是場(chǎng)地開挖過程中開挖區與回填區的分界線。通過計算每條分界線所包圍成的封閉(bì)區域的體積大小，就可以計算每一個開挖區與回填區的體積
，然後對每個開挖區與回填區進行累加，即爲所求區域裏土石方的挖方與填方量。因此，改造前後兩期地形的土石方量計算主要流程如下。

(1)将改造前的點雲數據DEM 導出爲文本格式，並(bìng)導入Autodesk的Civil3D軟件中，生成TIN曲面模型。對改造前的地形進行分析，統計各處(chù)最大高程、最小高程以及各高程區間範圍。

(2)将改造後的點雲數據DEM 導出爲文本格式，並(bìng)導入Autodesk的Civil3D軟件中，生成TIN曲面模型。對改造後的地形進行分析，統計各處(chù)最大高程、最小高程以及各高程區間範圍。

(3)在正射影像上圈出需要計算土石方量的區域，並(bìng)用確(què)定的邊界對該區域地形的TIN曲面模型進行提取。無人機航測可同時得到測量區域影像，可直接在影像上確(què)定土石方量的計算範圍。

(4)在Civil3D中打開需要計算土石方量區域的兩期地形TIN 模型，利用軟件中的地形分析工具，以改造前的地形爲參(cān)考基準，計算改造後(hòu)地形的填方量和挖方量(圖5)。

(5)統計分析兩期地形變(biàn)化，得出填挖方量結果，輸出土石方量報(bào)告：投影面積14　814.51m2 ，挖土石方量18740.80m3 ，填土石方量47.28m3。

3　精度分析

爲瞭(le)分析無人機航測技術測算土石方量的精度，本文将該方法與傳統GPS測量法進行比較。採(cǎi)用傳統GPS-RTK方法獲取測區範圍内兩期地貌特征點各60個，記錄各點平面坐标、高程及每個點對應的實地位置作爲檢查點。從兩期地形模型中提取各檢查點對應位置的模型點坐标與檢查點坐标作比對
，進行精度評估。

先利用兩期測(cè)量的各60個檢查點(diǎn)分别對兩次航測(cè)的平面、高程進行精度分析
，結果爲地形改造前平面中誤差±0.39m，高程中誤差±0.32m，改造後平面中誤差±0.36m，高程中誤差±0.29m，滿足1∶500大比例尺低空數字航空攝影成圖要求。

再利用CASS9.0軟件中土方計算功能完成異常點剔除、計算邊界界定以及基準確(què)定等工作後，基於(yú)方格網法計算改造前後的土石方量。然後将傳統GPS測量法計算得到土石方量結果與無人機航測計算結果進行比較
，結果爲兩種方案的挖方量相差80.95m3 ，占總挖方量的0.4%，填方量相差1.59m3，占總填方量的3.3%。

比較分析結果顯示，採(cǎi)用無人機航測得到的土石方量精度可以滿足土石方量計算的一般工程測量規範要求，可用於(yú)工程項目土石方量計算。

4　結　語

研究瞭(le)採用無人機航測技術進行土石方量計算的流程以及數據處理法，結果表明，将無人機航測技術運用於(yú)工程土石方量計算有如下優勢：

① 相對於(yú)傳統土石方測(cè)量方法，無人機航測(cè)技術更加機動、靈活，不受地形限制，在平緩、陡峭地區均适用;

② 數據採(cǎi)集更加快速，傳統方法數據採(cǎi)集通常需要數周，該方法一般僅需1d就能完成，特别是當測(cè)量面積較大時其優勢更加明顯;

③ 在測得地形高程數據同時，該方法獲取瞭(le)影像數據，可更加精確(què)界定土石方的計算範圍，使計算結果更加精準;

④ 所得的DEM 爲數字形式，可直接導(dǎo)入商業軟件中進行計算分析，提高瞭(le)計算效率;

⑤ 無人機攝影測(cè)量減少瞭(le)人員投入，減輕瞭(le)外業工作量，節約瞭(le)生産成本。