基于ArcGIS的土方量計算與可視化

沈名威，徐良，袁偉俊

(1.常州市測繪院，江蘇常州 213002； 2.金壇市城鄉(xiāng)測繪院，江蘇金壇 213200)

基于ArcGIS的土方量計算與可視化

沈名威1?，徐良1，袁偉俊2

(1.常州市測繪院，江蘇常州 213002； 2.金壇市城鄉(xiāng)測繪院，江蘇金壇 213200)

利用ArcGIS三維構(gòu)面的特性，形象的構(gòu)建地表、河道現(xiàn)狀立體模型，同時根據(jù)場地、河道的設(shè)計標(biāo)高，構(gòu)建立體模型，然后選定統(tǒng)一參考面，計算兩個模型到參考面的體積之差，可以快捷、準(zhǔn)確地計算出挖填土方量。實踐證明，這一算法較傳統(tǒng)的根據(jù)等高線、高程點計算土方量的方法更加嚴(yán)密、準(zhǔn)確、便捷。同時這一方法能根據(jù)要求靈活調(diào)整、分割土方計算區(qū)域，不需要對調(diào)整后的邊界進行外業(yè)高程點數(shù)據(jù)處理就能快速實現(xiàn)分區(qū)域土方量計算。

ArcGIS；TIN；DTM；土方量算；分割量算

1 引 言

工程土方量(挖方與填方)的精確性直接影響到工程費用概算和工程甲方、乙方的利益，甚至產(chǎn)生矛盾糾紛。筆者參與過多個工程的土方量計算工作，對斷面法、方格網(wǎng)法、等高線[1]、DTM(不規(guī)則三角網(wǎng))法[2]等幾種常用的計算方法在不同的計算機軟件(AutoCAD類軟件與Arc-GIS類軟件)中做了嘗試和對比，實踐證明，在這幾種方法在AutoCAD類軟件中實現(xiàn)土方量計算都有一個共同的缺陷:當(dāng)需要對一個施工范圍進行分區(qū)計算土方量時，分區(qū)計算的土方量之和跟直接計算的總施工范圍土方量總是相差很大[3]。原因是在AutoCAD類軟件中這幾種算法都是依據(jù)外業(yè)測量數(shù)據(jù)構(gòu)DTM，要實地在分區(qū)邊界上測量采集高程信息，當(dāng)邊界發(fā)生變化時，需要對新的邊界重新實測數(shù)據(jù)，從而增加工作量；而且會導(dǎo)致分區(qū)構(gòu)建的計算模型與整個施工范圍的總計算模型不一致?；贏rcGIS類軟件從整個施工范圍一次構(gòu)建土方量計算模型，分區(qū)范圍可以任意調(diào)整，可以有效克服AutoCAD類軟件中計算的不足，還能為工程相關(guān)方提供形象、直觀的圖像表達。本文以ArcGIS為例，詳細描述了整個計算過程，從中我們可以看出基于ArcGIS的土方量計算有其準(zhǔn)確、便于調(diào)整分割等優(yōu)勢。

2 基于ArcGIS的工程土方計算原理

ArcGIS下土方計算的原理如下:①確定挖填方總邊界；②外業(yè)采集高程數(shù)據(jù)并導(dǎo)入；③構(gòu)TIN；④模型編輯與高程檢核；⑤選擇投影面，計算總的挖填方；⑥根據(jù)實際分割要求在不改變TIN的前提下任意計算分割部分的土方量。

ArcGIS軟件的3D分析通過對離散高程點構(gòu)建不規(guī)則三角網(wǎng)(TIN)，計算挖、填方時采用三棱柱法。在構(gòu)建TIN的過程中，添加多種條件對即將生成的TIN進行約束和細分，根據(jù)TIN的形狀，能很好地還原現(xiàn)場原貌，根據(jù)模型形狀能發(fā)現(xiàn)外業(yè)測繪時高程數(shù)據(jù)采集中產(chǎn)生的粗差和遺漏，從而剔除外業(yè)粗差。以達到和實際情況盡可能接近的目的。

從整個過程可以看出ArcGIS在保持構(gòu)TIN不變的情況下，根據(jù)標(biāo)段分割計算各標(biāo)段土方量相當(dāng)快捷簡便。最重要的是該方法能保持同一模型下的分割計算，而不會出現(xiàn)AutoCAD下分割計算需要調(diào)整構(gòu)TIN邊界、重新構(gòu)TIN等步驟，最終導(dǎo)致總土方量跟分割后計算的土方量之和不一致的現(xiàn)象。

3 ArcGIS用于工程土方計算的過程

下面，以某一湖泊清淤工程為例，說明ArcGIS下土方計算過程。

某河道清淤工程，面積為1.51 km2。根據(jù)工程實施需要，該項工程分5個標(biāo)段進行施工，需要對每個標(biāo)段分別預(yù)算工程土方量。具體分割如圖1所示，設(shè)計清淤后標(biāo)準(zhǔn)斷面如圖2所示。

圖1 工程清淤前現(xiàn)狀與標(biāo)段劃分

圖2 工程設(shè)計清淤后斷面

(1)確定挖填方邊界

確定挖、填方邊界很重要，邊界確定后，就明確了外業(yè)采集離散高程點的大致范圍，外業(yè)高程采集時，邊界線上的每個頂點要有高程信息，確保在ArcGIS中建立的模型邊界準(zhǔn)確，插值后的計算模型在邊界上準(zhǔn)確。

(2)外業(yè)采集高程數(shù)據(jù)并導(dǎo)入

數(shù)據(jù)測量外業(yè)離散點高程精度是影響構(gòu)建TIN的重要因素，所以在實地采集高程信息時，儀器精度要高，可以使用全站儀或RTK，同時注意以下幾點:

①工程范圍邊界線上的高程點要能翔實反映高程變化，適當(dāng)在范圍線外采集高程點；

②工程范圍內(nèi)的高程點分布要均勻；

③對工程范圍內(nèi)的特征地物要采集范圍線及地形變化處的高程等特征信息。

測量成果可以采用AutoCAD的dwg格式或Excel的sheet導(dǎo)入到ArcGIS中，形成高程離散點Shp點文件和特征范圍線Shp線文件。

(3)構(gòu)TIN

構(gòu)TIN的過程在ArcGIS中操作起來相對簡單，在這里不做具體的闡述，圖3是對開挖前的測量成果構(gòu)建的TIN的分色表達，藍色部分的高程最低。

圖3 開挖前高程離散點構(gòu)建的TIN

(4)模型編輯與高程檢核

該部分工作是在TIN建立后有了一個圖像的界面，從分色圖上可以明顯的瀏覽出哪些高程點有異常的情況，對高程明顯異常的點，要檢查原始數(shù)據(jù)，確定高程值有誤的話，要刪除該高程點，再重新構(gòu)TIN[4]。在該土方量計算項目中，我們在施工范圍內(nèi)選擇了一些高程點不參與構(gòu)TIN，然后在構(gòu)建好的TIN的對應(yīng)位置采集插值后的高程信息，和測量的高程離散點進行了對比后發(fā)現(xiàn)，除TIN邊緣上的個別點插值后的高程與實測點的高程值有明顯差異(最大高程差為1.8 cm)，其他區(qū)域的點與檢查點的高程差不明顯。

(5)挖填方計算

有了上面4步的成果，在ArcGIS的面積與體積統(tǒng)計界面中，輸入計算參考面高程，即可快速計算出每個分區(qū)的挖方、填方體積，這里計算出來的挖、填方是該參考面高程中的挖、填方平衡后的一個凈挖方、填方值。圖4是根據(jù)設(shè)計標(biāo)準(zhǔn)斷面構(gòu)建的TIN模型，將現(xiàn)狀TIN和設(shè)計標(biāo)準(zhǔn)斷面TIN計算的高程基準(zhǔn)面選擇一致(這個值要小于設(shè)計的最低高程值)，用現(xiàn)狀TIN計算的結(jié)果減去設(shè)計斷面TIN計算結(jié)果即為該范圍內(nèi)的挖方量。給工程甲方提供的可視化模型圖片，如圖5所示。

圖4 設(shè)計斷面TIN

圖5 開挖清淤后立體高程模型

(6)預(yù)算、實際測算的土方量如表1、表2所示。

錢資湖清淤工程一期土方預(yù)算量成果 表1

錢資湖清淤工程一期土方量算量成果 表2

由表1、表2結(jié)合清淤前后的模型可以看出，實際開挖清淤時，標(biāo)段1和標(biāo)段2由于保留了幾個湖中島嶼，而設(shè)計時這幾個島嶼是計劃清淤一起挖掉的，所以在這個區(qū)域的實際清淤土方量明顯少于預(yù)估土方量。標(biāo)段3在湖的南邊清淤不到位，所以實際工程土方量少于預(yù)估土方量。標(biāo)段4由于實際淤泥厚度超過了設(shè)計標(biāo)高，同時清淤比較到位，最終出現(xiàn)了超挖的情況，因此土方量也比預(yù)計的有所超出。標(biāo)段5的實際清淤情況跟設(shè)計沒有變化，因此實際清淤土方量也跟設(shè)計預(yù)估土方量比較接近。這些工程實際情況也反過來證明了ArcGIS下土方計算的可靠性。

4 ArcGIS模型糾錯的特點

實踐證明，這一方法能根據(jù)外業(yè)高程數(shù)據(jù)在構(gòu)TIN的時候很好的控制土方邊界范圍，同時能很好還原現(xiàn)場原貌，根據(jù)模型形狀能發(fā)現(xiàn)外業(yè)測繪時高程數(shù)據(jù)采集中產(chǎn)生的粗差和遺漏[5]。實際計算過程中，我們發(fā)現(xiàn)根據(jù)現(xiàn)場采集的測量數(shù)據(jù)在生成DTM時，在幾個湖心島處DTM與現(xiàn)實不一致，根據(jù)現(xiàn)場測繪的地形圖形狀，完整的復(fù)原了這幾處湖心島的DTM，從而使參與計算的三維模型與實際保持一致，再加上ArcGIS下構(gòu)建三維立體的嚴(yán)密性，最終保證土方計算的正確性。

5 結(jié) 語

ArcGIS在工程土方量算上的應(yīng)用，解決了工程土方量分區(qū)計算的一大難題。為工程量算提供了一種靈活的計算方法。該方法的優(yōu)勢有兩點:一是能生成直觀的三維模型，通過模型對比，能發(fā)現(xiàn)外業(yè)測繪的粗差，可以根據(jù)地形圖數(shù)據(jù)調(diào)整DTM，從而使參與計算的三維模型與實際保持一致。二是ArcGIS提供了便捷的邊界調(diào)整功能，能夠在不調(diào)整數(shù)據(jù)的情況下，對總的工程范圍進行任意分割，單獨計算各個分區(qū)的工程土方量，真正實現(xiàn)在一個三維模型下的立體體積分割。保證了各個分區(qū)土方量與總土方量的嚴(yán)格一致，保證了計算成果的準(zhǔn)確性。

[1] 龔鍵雅.地理信息系統(tǒng)基礎(chǔ)[M].北京:科學(xué)出版社，2001.

[2] 慕永峰，朱昌永，李建.三角網(wǎng)結(jié)構(gòu)DTM的土方計算及應(yīng)用[J].測繪工程，2000，9(1):52～56.

[3] 王斌.DEM土方量計算的誤差探討[J].鐵路航測，1999 (3):21～23.

[4] 王永明.地形可視化[J].中國圖象圖形學(xué)報，2000(6).

[5] 魯成樹，吳次芳，汪峰等.農(nóng)地整理土方工程量計算及調(diào)配中地理信息系統(tǒng)的應(yīng)用[J].農(nóng)業(yè)工程學(xué)報，2003(6).

Earthwork Calculation and Visualization Base on ArcGIS

Shen Mingwei1，Xu Liang1，Yuan Weijun2
(1.Changzhou Surveying and Mapping Institute，Changzhou 213002，China；2.Jintan Surveying and Mapping Institute of Urban and Rural，Jintan 213200，China)

Taking advantages of ArcGIS in 3D modeling，it is easy to construct 3D models of earth surface and the current situation of rivers according to the designed space and elevation of rivers.And with a selected reference plane，the volume differences between twomodels can be calculated，which can help figure out the amountof cut and fill earthworks quickly and accurately.It has been proved by practice that the algorithm in this article ismuch more accurate and convenient than the traditionalmethods using contours and elevation points in calculating the amounts of earthworks.At the same time，thismethod can be adjusted flexibly and spilt the areas for earthwork calculations as required.It can help achieve the calculations of subregional earthworks quickly without processing the elevation data points outside the industry for the adjusted boundary.

ArcGIS；TIN；DTM；earthwork calculation；segmentation and calculating

1672-8262(2013)05-140-03

P209

B

2013—01—06

沈名威(1958—)，男，工程師，主要從事工程測量與測量管理工作。