摘要：土方量算在新興技術的發展下數據獲取與量算方式多種多樣，土方量算精度參差不一。文章闡述了通過傳統測繪、無人機傾斜攝影測量與無人機載激光雷達掃描對同一土方進行數據采集，不同測繪方式數據采集方法。采用南方Cass、Cass3D、Context Capture、DP-Modeler和CoProcess多種數據處理軟件進行土方量算，分析不同數據不同計算方法土方量算精度。分析結果表明：利用無人機傾斜攝影測量技術與無人機機載激光雷達掃描技術有效簡化了外業數據采集流程，數據處理土方量算以2.5m邊界采樣間距DTM法土方量算精度準確；植被剔除傾斜三維模型與地面點分類后點云模型直接量算法能快速高效地獲取土方量并滿足土方量算規范。
關鍵詞：土方；RTK；傾斜攝影；無人機載激光雷達；精度分析
1 引 言
計算機技術的發展為測繪帶來了更多的數據收集、提取和處理手段也給我們帶來了更多關于土方計算方法的思考。隨著測繪技術的不斷更新，擁有了多種獲取測繪數據的方式方法，在工程建設過程中土方量算的方式變得多種多樣。傳統的土方測量方法工作量大，不易在計算機上實現，不能有效利用現有的數據。在工程建設過程中土方量算精度是否準確直接關系到工程建設中各方面經濟利益，故要確保土方量算的準確性。土方量算中常用的量算方法有DTM(Digital Terrain Model; 數字地面模型)法、斷面法、方格網法和等高線法，不同的量算方法其量算結果也有差異。DTM法是利用所測土方特征點與高程點構建不規則三角網進行土方量算；數據所建立的不規則三角網模型能更好反映實際地形。斷面法主要是將實際地形以某分割路線分割若干橫斷面，通過計算橫斷面面積根據菱形體積計算方法量算土方量；斷面法更適用于沿縱向延伸及不規則地形土方量算。方格網法是利用一定邊長正方形格網對土方量算區域分割，內插計算各格網點高程及差值，通過高差擬合計算每個格網土方填挖方量，總和得到總土方填挖量；其計算相對簡單，適用性廣泛。等高線法是通過地形特征點形成一組地面高程等值線顯示地形，利用相鄰兩條閉合等高線內插高程所圍成的面積與高差計算土方量。DSM(Digital Surface Model)即數字表面模型，它不同于DTM數字地形模型，其包含了地表建筑樹木等除地面信息以外的其他特征點信息；對于無人機傾斜攝影測量所建立的三維模型與無人機載激光雷達點云模型實則為數字表面模型，通過對數字表面模型直接范圍量算可獲取土方量。對土方進行量算時，每種量算方法都有各自的優點和局限性，考慮各方面因素選擇合適的量算方法，達到最優的目的。土方量算以外業實測特征離散點坐標為基礎，通過建立相關土方計算模型來獲取土方量，因此，土方量的計算精度取決于所測量特征點位的精度和密度。本文將對研究區域采用DTM法、方格網法與DSM直接量算法對土方量算及準確性進行分析。對于土方數據采集，傳統采集方式主要是利用RTK(real time kinematic)技術或全站儀人工采集，測繪設備技術不斷更新后，基于無人機傾斜攝影測量技術與無人機載激光雷達掃描技術的數據采集迅速興起。無人機傾斜攝影測量技術以范圍廣、精度高等方式感知空間地理信息，建立的三維可視化數據更加真實地反映空間地理信息的地理位置、外觀、高度等屬性特征，從多方位、多角度地對地形地物進行數據采集獲取更多三維信息。擁有“全要素、真紋理”的全息表達能力、高精度、豐富的可測性、低成本、高效率等優勢。無人機載激光雷達掃描技術通過對不規則、復雜地表進行連續、快速、大面積、非接觸掃描，將單站掃描擴展為整體掃描，無人機機載激光雷達掃描技術能將地表上的點、線、面體、地表三維模型等幾何數據快速重構出來，可對調查對象的整體變化、細節特征進行詳細了解。
2 數據采集
本文以昆明市某高校的場地平整土方量算工程為例進行土方數據外業采集，研究區域平場面積約 7856 m2平場設計標高為 1883 m, 大部分區域為裸露土地帶有零星植被，研究區域現狀如圖1所示，研究區域分別通過RTK人工測量、無人機傾斜攝影測量與無人機機載激光雷達測量進行數據采集。
圖1 研究區域概況及范圍線
2.1 人工RTK數據獲取
傳統土方外業數據獲取主要以RTK技術人工測繪的方法進行地形點測量，需測量人員抵達每一個地形特征點，逐點采集獲取數據，其工作效率不高，地形特征點密度較低。本研究區域采用中海達RTK在已知3個地面控制點上進行參數校正，對研究區域以 5 m一個地形特征點進行數據采集。
2.2 無人機傾斜攝影數據獲取
無人機傾斜攝影測量技術是非接觸式的，外業數據采集具有速度快、效率高、精度高、范圍大的特點，采集得到的圖像分辨率高，經過處理可以生成三維表面模型。研究區域采用DJI PHANTOM 4 PRO消費級無人機對研究區域進行數據采集，外業數據采集流程如圖2所示。研究區域共布設了8個像控點，利用RTK以對中桿的形式對像控點進行坐標采集，采用昆明87坐標系。根據研究區域現場情況為滿足土方量算精度需求，航線規劃步驟中無人機飛行高度設置為 60 m，采用單鏡頭飛五鏡頭飛行模式，航向重疊率為90%，旁向重疊率為80%，相機云臺俯仰角度正射為90°傾斜為45°，經過對原始影像數據檢查篩選出1055張影像進行研究區三維模型建立。
圖2 無人機外業數據采集流程圖
2.3 無人機機載激光雷達數據獲取
激光雷達是新型激光傳感器。在對地觀測領域，其最初目的是為獲取高精度數字表面模型。經點云去噪、濾波，獲得剔除植被、人工建筑等以后的數字地面模型。與傳統的測繪方法相比，激光雷達顯示出許多有前途的優點，例如收集高密度和大規模點云的效率高，有利于地形的詳細表示。研究區采用華測無人機載激光雷達進行數據采集，并在現場進行采集數據下載、POS數據解算、點云數據解算和坐標轉換的預處理，外業數據采集流程如圖3所示。
圖3 無人機載激光雷達外業數據采集流程圖
3 數據處理
3.1 傳統數據處理方法
RTK技術人工測繪的土方外業數據采用南方CASS 10.1 軟件進行土方量算。為了對比不同數據的土方量精度的準確性，采用DTM法與方格網法，邊界采樣間距以 2.5 m、5 m、10 m 和 20 m 進行土方量算，量算結果如圖4所示。
圖4 RTK技術人工測繪土方內業數據處理圖
(1)三維傾斜模型建立
處理無人機航片數據時采用Context Capture軟件進行三維傾斜模型建立。Context Capture實景建模軟件是一套不需太多人工干預的通過影像自動生成高分辨率的三維模型軟件，模型建立可為后續模型修模提供支持。研究區域航片數據導入軟件中加載云臺相機參數并導入像控點數據進行影像刺點。刺點完成進行空中三角測量如圖5所示，其目的是通過測量相應的像控點等像點坐標以解析或數字的形式建立三維立體模型并進行數值解算。解算完成后進行三維模型構建，根據電腦運行內存研究區分為16個瓦片，設置其余參數后生成OSGB與OBJ格式的三維模型如圖5所示。
圖5 空中三角測量結果(左)與三維模型(右)圖
(2)模型土方量算
研究區域三維模型首先采用南方Cass 3D軟件進行土方量算。對模型進行地物特征點采集，同人工測繪的方法一致，對研究區域進行DTM法與方格網法的4種邊界采樣間距土方量算。建立的研究區域三維模型還可采用DSM直接量算法在Context Capture Viewer軟件土方量算如圖6所示。
圖6 DSM直接量算法傾斜模型土方量算
基于傾斜攝影測量技術構建的三維模型中含有懸浮物碎片、植被等地物對DSM直接量算法所得土方量的準確性產生較大干擾，因此利用DP-Modeler軟件對研究區域三維模型進行干擾因素如植被的剔除。DP-Modeler模型修飾在OBJ格式三維模型上進行，通過踏平、刪除、補洞及貼圖等操作，改善維懸浮物、道路凹凸、紋理拉花及水面缺失等問題。在DP-Modeler工程中，首先分別導入OSGB和OBJ格式瓦片數據，然后在OBJ數據上進行模型修飾，完成修飾的OBJ瓦片后再轉換成OSGB格式。修改后的三維模型再進行土方量算，模型修改前后如圖7所示。
圖7 模型植被剔除前(上圖)后(下圖)對比圖
3.3 無人機載激光雷達數據處理方法
對激光雷達數據后處理采用CoProcess華測點云后處理軟件，對研究區域數據進行分割刪除、噪點處理、壓縮采樣等。地面點分類是將點云濾除地面以上的地物，如植被建筑物等僅保留地面的點云，是土方量算不可或缺的步驟，處理后點云模型如圖8所示。處理后的點云模型可直接進行軟件土方量算。
圖8 處理后激光雷達點云數據圖
處理后的研究區域點云模型還可以輸出LAS點云格式，將點云數據導入Context Capture Viewer軟件進行OSGB格式的三維模型的建立，建立的研究區域三維模型可采用DSM直接量算法土方量算，量算結果如圖9所示。輸出后三維模型再在Cass3D軟件上，同人工測繪的方法一致，對研究區域進行DTM法與方格網法的4種邊界采樣間距土方量算。
圖9 DSM直接量算法模型土方量算圖
4 精度分析
4.1 模型精度檢驗
在研究區域選取的52個地物地貌特征點對人工測繪、無人機傾斜模型和機載點云生成的三維模型高程數據進行對比分析。根據《工程測量規范》(GB50026-2007)每幅圖檢測點數不少于20個且均勻分布，其高程中誤差公式為：
式中，Mh為模型的高程中誤差(m)；n為檢查點個數；△hi為檢測高程與模型高程較差(m)。其高程殘差精度統計圖如圖10所示。
圖10 高程殘差精度統計圖
人工測繪檢查數據與無人機傾斜模型數據高程中誤差計算得 0.050016，人工測繪檢查數據與機載點云生成的三維模型數據高程中誤差計算得 0.065691；根據GB50026-2007《工程測量規范》數字高程模型格網點的高程中誤差限差為1/3hd(hd為基本等高距；一般地區)即研究區域高程中誤差限差為0.3，模型精度滿足規范要求。
4.2 土方量算精度分析
為分析不同數據采集的數據土方量精度的準確性，通過Cass.10.1 Cass3D軟件對多種模型進行DTM法與方格網法的4種邊界采樣間距土方量算，其DTM法量算結果如表1所示，方格網法如表2所示。
由量算表格可知DTM法傾斜攝影測量模型最大挖方誤差率為1.4%，最大填方誤差率為3.5%，機載激光雷達模型最大挖方誤差率為0.7%，最大填方誤差率為2.2%；方格網法傾斜攝影測量模型最大挖方誤差率為1.5%，最大填方誤差率為3.3%，機載激光雷達模型最大挖方誤差率為1.9%，最大填方誤差率為3.9%滿足土方量算要求；通過兩種土方量算方法看出邊界采樣間距越小，模型量算越可靠。
而通過DSM模型直接量算法對無人機傾斜模型量算結果挖方量為 23571.925 m3填方量為 1 923.186 m3，無人機載激光雷達模型量算結果挖方量為 23 537.354 m3，填方量為 1938.500 m3；對模型進行修飾、植被剔除后，土方直接量算精度可靠。
5 結 論
土方量算在工程建設過程中是必不可少的環節，土方量算精度是否準確直接關系到工程建設中各方面經濟利益。本文研究了基于多種測繪技術對同一土方量算方法，探討了不同數據模型在不同邊界采樣間距與模型直接量算的土方量算精度的準確性，其中以DTM法所計算的 2.5m 邊界采樣間距的各類模型土方量算精度最高。通過研究區域土方量算結果看出，基于無人機傾斜攝影測量技術與無人機載激光雷達技術的數據采集比傳統數據采集方式其作業迅速，勞動強度低，成本可控。在工程建設過程中以傳統技術手段與新興技術相結合為土方量算提供了高效的數據支撐。
來源：城市勘測 2022(02) 作者：太浩宇 夏永華
掃一掃關注我微信
© 2016-2020 河南中圖空間信息技術有限公司 技術支持：木薪工作室 Designed By 太子召
聯系地址：鄭州市金水科教園區楊金路199號新科技產業園10號樓606室 | 營業執照 |