無人機傾斜攝影在土方測量中的應用研究

林金標

（福建省地質測繪院，福建 福州 350011）

0 引言

目前，采用無人機傾斜攝影技術對土方量進行測量，克服了傳統土方測量易受復雜地形條件影響，造成工作緩慢、成本提高的問題。無人機傾斜攝影具有高分辨率、生產效率高等優點，能保證整個項目數據的一體性，大大提高土方測量的速度與精度。

1 土方計算原理

一般采用南方CASS10.1軟件對土方進行計算，常用的土方計算方法有DTM法、斷面法、方格網法、等高線法等。

1.1 DTM法土方計算

DTM法適用于地形復雜、不規則形狀區域的計算，但不適合地勢變化較大區域的計算。因軟件生成的三角網不能與高低起伏的復雜地形地貌相吻合，要由人工對一小部分三角網進行修改，使之與高低起伏的地形地貌貼合。通過對開工前和填挖后兩期地面點坐標和高程進行實測，利用軟件分別生成兩期三角網，將指定范圍分割為若干個三棱錐，并進行編號，計算每個三棱錐的體積，最后通過累加每個空間立體的體積，計算得到指定范圍內的填（挖）方量。

1.2 斷面法土方計算

斷面法適用于地形狹長的帶狀區域，計算精度取決于各條橫斷面的里程間距，橫斷面間距越小，計算精度就越高。對高低起伏變化較大的地形，應縮小橫斷面間距，增加橫斷面條數。通過獲取縱斷面線上每個里程處的實測橫斷面線與設計橫斷面線所組成的橫斷面圖，從而求出相鄰里程兩條橫斷面面積的平均值，再乘以這兩條橫斷面的間距，可計算出每兩條相鄰斷面間的土方量，匯總所有兩條相鄰斷面間的土方量，即可計算出該地塊的總填（挖）方量。

1.3 方格網法土方計算

方格網法適用于地勢平坦、起伏不大的地形或測區范圍大的區域，但不適合地形起伏變化較大的區域。方格網的寬度越小，計算精度就越高。通常情況下，方格網四個角點的高程是周邊點高程通過內插的方式得來的，不考慮地形變化。通過預先輸入方格網邊長，將計算范圍劃分成一個個方格，把每個方格四個角的高程相加后取其平均值，并與設計標高相減得出挖（填）高度，再把這個挖（填）高度乘以對應方格的面積，從而計算出各個長方體的體積，最后對所有長方體的體積進行累加，即可計算出該地塊的總填（挖）方量。

1.4 等高線法土方計算

等高線法適用于土方概算、地形坡度大的區域，該方法可計算任意兩條等高線之間的土方量，但所選等高線必須閉合才有效。通過軟件自動求出相鄰兩條等高線所圍的面積，即可得到兩條等高線之間的高差，根據相鄰等高線所圍的面積乘以其高差，從而計算出每相鄰兩條等高線間的土方量，最后累加所有相鄰兩條等高線間的土方量，即可計算出該地塊的總填（挖）方量。

2 無人機傾斜攝影的應用

本研究以福建省莆田市某地塊的土方測量為例，該地塊東西長約290 m，南北長約170 m（見圖1），該地塊現為旱地，植被稀疏低矮，地勢相對平坦，最低處高程約為21.0 m，最高處高程約為26.0 m，屬地勢平坦的地形，故采用方格網法來計算土方量。

2.1 外業數據采集

2.1.1 外業像控點布設。在測區四個角點向外側約25 m處均布設像控點，其他區域按兩點間距為150～180 m均勻布設8個像控點。像控點一般布設在清晰的線狀地物交點處、明顯地物的折角頂點；在無明顯地物標志的區域布點時，要在地面上埋設不小于30 cm×30 cm的水泥樁，且在水泥樁面上用紅油漆刷出“十”字形，十字中心要加噴直徑為4 cm的圓點作為像控點標識。

2.1.2傾斜多視影像獲取。本研究采取無人機傾斜攝影方法進行外業影像數據的采集，設備采用大疆無人飛行器DJI Phantom 4RTK，相對航高值設為50 m，測區平均高程位置影像地面分辨率優于1.5 cm，最低點影像分辨率優于2.0 cm。按照五向航線進行飛行，航向重疊度設為85%，旁向重疊度設為70%。為確保所攝影像能全面覆蓋作業區域，航向及旁向覆蓋均超出測區邊界線一個相對航高值的范圍，共計獲取傾斜多視影像800多張。

2.1.3 像控點測量。像控點測量分為平面控制測量和高程控制測量，坐標系采用2000國家大地坐標系，高程基準采用1985國家高程基準。在無人機進行傾斜攝影時，大疆無人飛行器可同步獲取像控點的平面坐標，在完成對像控點平面控制測量后，采用水準儀按四等水準精度進行高程控制測量。

2.2 內業數據處理

2.2.1 實景三維模型生產。將像控點坐標數據、傾斜多視影像數據、POS數據導入到Context?Capture三維建模軟件平臺中，完成實景三維模型數據的生產。首先，在ContextCapture三維建模軟件中加載所有像控制點的坐標數據和多視影像數據。其次，區域網整體平差以光束法為基礎，以中心投影的共線方程為一個平差基礎方程，一張影像組成的一束光線作為一個平差單元，再將各光線束在空間中進行平移和旋轉，使模型間的公共光線實現最佳交會，把全部區域納入到控制坐標系里，進而利用影像密集匹配構建三角網（TIN），將紋理映射到對應的三角網格中，即可生成實景三維模型。

2.2.2 土方量測算。①高程點數據處理。基于生成的實景三維模型數據，通過軟件將其轉換為高密度的高程點數據，從而提高土方量計算的準確性。判斷實景三維模型上的高程點數據是否滿足地形地貌的變化情況，如不滿足，應補充采集明顯地貌變化處的高程。②土方量測算結果。該案例地勢平坦、起伏不大，適合采用南方CASS10.1軟件方格網法進行土方計算。經計算得出，挖方量為180 493.6 m3，填方量為0.0 m3。

3 精度評估

為評估無人機傾斜攝影開展土方測量的可行性、可靠性，采用像控點平面坐標、碎部點高程、土方量3個指標對案例測量的精度進行比較驗證。采用華測RTK儀器測量的像控點平面坐標值與傾斜攝影實景三維模型采集像控點的平面坐標值進行比較（見表1）；采用圖根水準測量的碎部點高程值與傾斜攝影實景三維模型采集碎部點的同點位的高程值進行比較（見表2）；采用方格網法土方計算方法、全野外數字化測量高程計算的土方量與傾斜攝影實景三維模型采集高程數據計算的土方量進行比較（見表3）。

表3 全野外數字化測量與實景三維測量的土方量比較

由表1可知，利用華測RTK測量像控點的平面坐標值與傾斜攝影實景三維模型采集像控點的平面坐標值之間較差的中誤差為±2.6 cm，最大為±3.6 cm，滿足規范規定點位中誤差不應大于±5.0 cm的精度要求。

表1 華測RTK測量與實景三維采集像控點的坐標值比較

由表2可知，利用圖根水準測量的碎部點高程值與同點位實景三維采集的內插高程值之間較差的中誤差為±4.1 cm，最大為±6.0 cm，滿足《城市測量規范》中土石方測量規定的檢測高程與原測內插高程的平均較差應在±10.0 cm內的精度要求。

表2 圖根水準測量與實景三維采集碎部點的高程值比較

由表3可知，采用同樣的土方計算方法，采用全野外數字化法測量的土方量與實景三維模型測量的土方量之間較差如下，填方量較差為0.0 m3、差值百分比為0%，挖方量較差為?4 210.2 m3、百分比較差為?2.33%，均滿足《城市測量規范》中規定的檢核計算成果與原計算成果的較差應不大于原計算成果3%的精度要求。

綜上所示，采用無人機傾斜攝影實景三維模型進行土方測量的計算結果是可靠的，滿足實際需求，且具有較高的實用性。

4 技術優點

與傳統的土方測量方式相比，無人機傾斜攝影在土方測量中具有靈活機動的優點，且不受地形的限制，在平緩或陡峭地區均適用，除植被茂密、密集建筑群等區域外，一般不用人工去現場補充測量。外業測量轉為內業解析，工作更加高效快捷，減少人員的投入，減少外業的工作量，提高生產效率，節約生產成本。在地形地貌特征明顯變化處無標注高程值時，可通過實景三維模型在地形地貌明顯特征變化處直接點取其高程值，有針對性地補充采集遺漏處的高程值，避免傳統土方測量方法容易出現高程點采集密度不夠的問題，從而對土方量計算的精度造成影響。無人機傾斜攝影模型數據客觀、真實，通過配套的軟件可進行高度、長度、面積、角度的量測，能提供高精度、高真實感的高密度三維點云、實景三維模型等多元可視化數據。可獲取DEM格式的高程點數據，可導入Civil 3D三維設計軟件中進行分析處理，減少中間作業環節，真正實現內、外業一體化的作業流程。測量精度更可靠，處理環節更加簡便，軟件能對數據處理過程自動檢查，減輕作業人員的工作強度，避免人為造成的錯誤。

5 結語

本研究詳細介紹無人機傾斜攝影在土方測量中應用的全過程，通過與全野外數字化測量在像控點坐標、碎部點高程、土方量等3個方面的差值進行比較，驗證無人機傾斜攝影土方測量的精度滿足規范要求，計算結果是可靠的，且具有較高的實用性和可行性，在一定程度上能優化土方測量及計算流程，提高土方量的計算精度，但對被茂密植被、高大建筑物遮擋區域對無人機傾斜攝影土方測量計算精度的影響程度還要進一步研究。