基于WebGL的土方在線自動(dòng)算量方法研究

文｜山東同圓數(shù)字科技有限公司 孫梅強(qiáng) 侯曉穎

0 引言

土方量的測(cè)量和計(jì)算作為工程建設(shè)前期階段一項(xiàng)不可或缺的工作，土方工程量在整個(gè)的項(xiàng)目的工程量中往往占有相當(dāng)?shù)谋壤渲型恋仄秸こ掏顿Y占總投資的比例甚至可達(dá)40%～80%[1]，其計(jì)算結(jié)果的精確性和準(zhǔn)確性直接關(guān)系到工程量預(yù)算、施工組織設(shè)計(jì)和施工現(xiàn)場(chǎng)的工作安排，對(duì)工程建設(shè)具有重要影響[2]。因此，如何在保證計(jì)算精度、效率基礎(chǔ)上，簡(jiǎn)化土方算量流程，降低成本，擺脫平臺(tái)限制成為一個(gè)值得研究的問題。

WebGL 技術(shù)的出現(xiàn)和硬件性能不斷提升，使得Web 端加載三維場(chǎng)景變得容易，且不需要安裝任何插件或軟件，被越來越多的應(yīng)用到各個(gè)領(lǐng)域。例如，曾應(yīng)用WebGL 技術(shù)對(duì)機(jī)載激光點(diǎn)云數(shù)據(jù)三維可視化進(jìn)行了研究，結(jié)果表明在網(wǎng)絡(luò)瀏覽器中能夠提供很好的三維可視化體驗(yàn)[3]；其后周達(dá)、邢峰、朱海林等學(xué)者分別在WebGIS、激光點(diǎn)云、BIM 可視化等領(lǐng)域應(yīng)用WebGL 技術(shù)做了不同程度的研究[6-8]。WebGL 技術(shù)研究的不斷深入為線上實(shí)現(xiàn)土方量計(jì)算提供了底層支持，傳統(tǒng)的土方量計(jì)算方法主要有方格網(wǎng)法、斷面法和等高線法，輔以相應(yīng)的南方CASS、MapGIS、ArcGIS 等桌面應(yīng)用程序進(jìn)行計(jì)算，方?。?012）等使用ArcGIS 軟件計(jì)算了實(shí)際項(xiàng)目的土方工程量，結(jié)果表明計(jì)算精度與DEM 模型精度密切相關(guān)[10]；昂龍等則使用Autodesk Civil 3D 軟件對(duì)施工場(chǎng)區(qū)內(nèi)的土方進(jìn)行數(shù)字信息建模和挖填土方量模擬計(jì)算，可快速、準(zhǔn)確地計(jì)算土石方量，以達(dá)到挖填平衡效果[11]。

鑒于以上意義和背景，本文在對(duì)比分析傳統(tǒng)土方量計(jì)算方法和方式的基礎(chǔ)上，通過研究各計(jì)算方法的適用性和精度，提出一種基于WebGL 技術(shù)的免插件土方在線自動(dòng)算量方法，可同時(shí)多項(xiàng)目并行管理，在此基礎(chǔ)上可進(jìn)行相關(guān)應(yīng)用拓展。最后以實(shí)際工程項(xiàng)目為例，進(jìn)行實(shí)踐驗(yàn)證，結(jié)果表明：該方法實(shí)現(xiàn)了預(yù)期目標(biāo)，簡(jiǎn)化了土方量計(jì)算分批、分步的繁瑣操作流程，計(jì)算效率大幅提高，成本明顯降低，同時(shí)提供直觀的三維可視化展示。

1 土方量計(jì)算方法與對(duì)比分析

土方工程量計(jì)算的本質(zhì)就是計(jì)算同一地塊開挖前后的填挖方量，即是計(jì)算土體的體積。根據(jù)不同的項(xiàng)目現(xiàn)場(chǎng)情況，目前土方量計(jì)算常用的方法有橫斷面法、方格網(wǎng)法、DEM 法等。

1.1 橫斷面法

橫斷面法是根據(jù)地形走勢(shì)或者根據(jù)規(guī)劃設(shè)計(jì)的重要軸線，沿著橫斷面線每間隔一定距離垂直向下剖切所形成的斷面。

上面公式中：i≥2,Li 為斷面間距； Ai-1和Ai 分別為第i 段起斷面和終斷面的填挖方面積；n 為斷面?zhèn)€數(shù)。

本方法通常適用于線性土方工程，對(duì)地形較為復(fù)雜區(qū)域適用性較高，

1.2 方格網(wǎng)法

方格網(wǎng)法是根據(jù)測(cè)繪數(shù)據(jù)獲取高程點(diǎn)的（X，Y，Z）坐標(biāo)生成規(guī)則格網(wǎng)，然后計(jì)算每一個(gè)規(guī)則格網(wǎng)的填挖方量，然后將每一個(gè)規(guī)則格網(wǎng)計(jì)算的填挖方量求和來得到特定地塊范圍的挖填方量，最后生成挖填方分界線。每一個(gè)方格內(nèi)的挖填方量是通過格網(wǎng)面積乘以高差來計(jì)算得到的，每個(gè)規(guī)則格網(wǎng)的面積是利用預(yù)設(shè)的規(guī)則格網(wǎng)邊長求得，而高差則是利用規(guī)則格網(wǎng)四個(gè)角點(diǎn)的高程與設(shè)計(jì)高程來求取，其公式如式[13]：

上面公式中：ij≥1,Hij 為第i 行j 列的規(guī)則格網(wǎng)的高程差；a、b 分別為規(guī)則格網(wǎng)邊長；n 為區(qū)域內(nèi)格網(wǎng)數(shù)量。

本方法一般適用于地形平坦、坡度變化較緩、起伏不大的大范圍土方量計(jì)算。

1.3 DEM 法

基于DEM 的土方工程量計(jì)算方法是通過建立數(shù)字模擬地形來實(shí)現(xiàn)。該方法幾乎適用于任何地形，且計(jì)算精度較高，但同時(shí)它的計(jì)算精度又與測(cè)繪成果的精確程度有關(guān)，并要求認(rèn)真分析和編輯依據(jù)高程數(shù)據(jù)建立的三角網(wǎng)，以保證三角網(wǎng)的準(zhǔn)確性。所以該DEM 法的缺點(diǎn)就在于模型建立較為復(fù)雜，需實(shí)地采集大量的數(shù)據(jù)點(diǎn)三維坐標(biāo)，占用大量的存儲(chǔ)空間，計(jì)算時(shí)間較長。

通過對(duì)上述土方量計(jì)算方法的對(duì)比研究分析，為達(dá)到土方算量方法普適性、計(jì)算結(jié)果準(zhǔn)確性的目標(biāo)，我們采用適用性最強(qiáng)的DEM 法作為本文土方量計(jì)算方法的基本原理，進(jìn)行基于WebGL 的土方量在線自動(dòng)算法研究。

2 實(shí)施路徑

圖1 技術(shù)路線圖

圖2 土方計(jì)算原理

2.1 技術(shù)路線

本文的技術(shù)路線主要分為數(shù)據(jù)預(yù)處理、數(shù)據(jù)格式在線轉(zhuǎn)換、建立數(shù)學(xué)模型、可視化展示（含報(bào)告）四個(gè)階段。

2.2 方法研究

2.2.1 土方量計(jì)算原理和算法模型

基本原理：通過地形數(shù)字高程模型（Digital Elevation Model，DEM）來計(jì)算土石方工程量的方法主要是依托各種測(cè)繪技術(shù)采集區(qū)域內(nèi)現(xiàn)狀地面點(diǎn)坐標(biāo)（X，Y，Z）和規(guī)劃設(shè)計(jì)地形高程，根據(jù)算法構(gòu)建區(qū)域內(nèi)不規(guī)則三角網(wǎng)或規(guī)則格網(wǎng)來計(jì)算每一個(gè)小區(qū)域內(nèi)的填方工程量與挖方工程量，最后累計(jì)求和即可得到規(guī)定范圍內(nèi)總體填方與挖方的工程量，并生成區(qū)域內(nèi)填方與挖方區(qū)域分布圖。

算法模型：計(jì)算每個(gè)填/挖區(qū)域的土體體積。其中針對(duì)單個(gè)計(jì)算單元，體積公式為：

其中，

Voli表示單個(gè)計(jì)算單元面積；

Areaj表示單個(gè)計(jì)算單元面積；

ΔZ 表示同一計(jì)算單元現(xiàn)狀與設(shè)計(jì)高程的差值；

ZBefore表示現(xiàn)狀高程值；

ZAfter表示設(shè)計(jì)高程值。

2.2.2 數(shù)據(jù)輕量化處理

現(xiàn)狀數(shù)據(jù)、設(shè)計(jì)數(shù)據(jù)以及用地地塊作為土方算量計(jì)算的數(shù)據(jù)基礎(chǔ)，數(shù)據(jù)本身的體量、精度直接影響著土方量計(jì)算的精度和效率?，F(xiàn)狀數(shù)據(jù)往往需要轉(zhuǎn)換成點(diǎn)云數(shù)據(jù)進(jìn)行計(jì)算，而點(diǎn)云數(shù)據(jù)最為突出的特點(diǎn)是點(diǎn)數(shù)量多，精度高，但數(shù)據(jù)量極其龐大，對(duì)軟件、硬件都較高的要求。

但實(shí)際工作中點(diǎn)云數(shù)據(jù)本身并非所有數(shù)據(jù)都是對(duì)曲面構(gòu)建有用的，為了提高效率以及獲得更好的網(wǎng)格化效果，我們通過設(shè)定相關(guān)標(biāo)準(zhǔn)算法，對(duì)無效數(shù)據(jù)進(jìn)行清洗，地勢(shì)起伏較小的數(shù)據(jù)進(jìn)行數(shù)據(jù)抽稀，然后通過“分布式計(jì)算”的方式，實(shí)現(xiàn)土方計(jì)算所需數(shù)據(jù)的輕量化處理。付瑋在2015年提出了基于三維激光點(diǎn)云數(shù)據(jù)的輕量化計(jì)算方法，將劃定范圍內(nèi)局部區(qū)域和整體區(qū)域采集的數(shù)據(jù)進(jìn)行互相融合，該方法既考慮到了局部細(xì)節(jié)，又考慮到了整體特征，可以更加有效的對(duì)所采集到的三維激光點(diǎn)云數(shù)據(jù)進(jìn)行輕量化處理，我們可以引用到對(duì)于地勢(shì)起伏較大的地形點(diǎn)云數(shù)據(jù)的輕量化處理中[14]。

圖3 角度偏差法

圖4 非均勻網(wǎng)格法簡(jiǎn)化數(shù)據(jù)點(diǎn)

圖5 點(diǎn)云輕量化算法流程

通過“坡度濾波”算法進(jìn)行數(shù)據(jù)抽稀，該方法計(jì)算邏輯簡(jiǎn)單，但是需要預(yù)先知道地形坡度，而且對(duì)于每一個(gè)點(diǎn)進(jìn)行坡度計(jì)算，這樣勢(shì)必造成計(jì)算量大進(jìn)程緩慢的問題，同時(shí)高差閥值的選擇是整個(gè)算法的關(guān)鍵，這些過濾閥值的設(shè)置取決于測(cè)區(qū)的實(shí)際地形狀況設(shè)置不同的過濾參數(shù)值。進(jìn)行兩次濾波計(jì)算，可以對(duì)數(shù)據(jù)進(jìn)行更深度的輕量化處理。

2.2.3 多方法封裝達(dá)到自動(dòng)化計(jì)算

圖6 在線計(jì)算與線下處理流程對(duì)比

鑒于目前土方量計(jì)算流程的繁瑣及效率低下等問題，本研究通過將坐標(biāo)動(dòng)態(tài)投影、數(shù)據(jù)轉(zhuǎn)換模型、土方量計(jì)算算法的統(tǒng)一集成，用戶只需要通過瀏覽器登錄系統(tǒng)，根據(jù)提示上傳標(biāo)準(zhǔn)化后的現(xiàn)狀數(shù)據(jù)、設(shè)計(jì)數(shù)據(jù)和用地地塊（含紅線范圍），即可通過在線入口實(shí)現(xiàn)土方量的自動(dòng)化計(jì)算。

本研究的自動(dòng)化計(jì)算流程中“坐標(biāo)動(dòng)態(tài)投影模型”的實(shí)現(xiàn)采用以下公式[15]：

2.2.4 Web 端在線豎向設(shè)計(jì)調(diào)整

精確的土方計(jì)算不僅可以量化土方工程量，還可以為設(shè)計(jì)師豎向設(shè)計(jì)方案提供數(shù)據(jù)支撐，避免施工階段后期調(diào)整困難，達(dá)到節(jié)約成本的目的。本研究提出了一種有效的節(jié)約時(shí)間成本方法，在Web 端面實(shí)時(shí)、動(dòng)態(tài)調(diào)整局部豎向設(shè)計(jì)高程，該方法既考慮了設(shè)計(jì)師設(shè)計(jì)過程對(duì)圖紙的調(diào)整，又考慮了溝通過程中信息折損，可以有效提高工作效率，并可通過窗口對(duì)兩次或多次計(jì)算結(jié)果進(jìn)行可視化展示和數(shù)據(jù)對(duì)比，同時(shí)生成豎向設(shè)計(jì)方案調(diào)改記錄、土方算量統(tǒng)計(jì)分析報(bào)告，單位公頃點(diǎn)云面積的計(jì)算效率按地勢(shì)復(fù)雜情況可達(dá)0.08 min～1 min。

2.2.5 基于WebGL 的三維可視化

圖7 調(diào)整前

WebGL 技術(shù)的應(yīng)用解決了Web 3D 應(yīng)用面臨的插件和圖形渲染問題，使得大場(chǎng)景、大規(guī)模、多類型三維數(shù)據(jù)在Web 端的快速呈現(xiàn)成為現(xiàn)實(shí)，這也是本文最重要研究?jī)?nèi)容之一?；贐/S 架構(gòu)，采用WebGL技術(shù)構(gòu)建了一種無插件的渲染方法。通過測(cè)試，本文將土方算量結(jié)果數(shù)據(jù)轉(zhuǎn)化為Geojson 的數(shù)據(jù)結(jié)構(gòu)，加載到三維場(chǎng)景中，其結(jié)果表現(xiàn)出良好效果，能夠直觀展示區(qū)域內(nèi)挖填分布及方量情況，生成紅線內(nèi)給地塊土方填挖量結(jié)果，并且生成項(xiàng)目土方計(jì)算報(bào)告。

圖8 標(biāo)高調(diào)整

圖9 方案對(duì)比

圖10 三維可視化

3 應(yīng)用實(shí)例

3.1 項(xiàng)目介紹

濟(jì)南萬達(dá)文化體育旅游城項(xiàng)目位于濟(jì)南市高新區(qū)，面積為357.56hm2。開發(fā)商提供了項(xiàng)目首開區(qū)施工現(xiàn)場(chǎng)實(shí)際土方開挖數(shù)值，挖方量126628 m3，填方量156588 m3，填挖總量為283216 m3。

3.2 土方量計(jì)算及結(jié)果

通過激光點(diǎn)云掃描獲取了項(xiàng)目范圍現(xiàn)狀0.5m 精度的現(xiàn)狀點(diǎn)云數(shù)據(jù)，根據(jù)設(shè)計(jì)豎向方案圖紙?zhí)崛∝Q向設(shè)計(jì)數(shù)據(jù)以及用地范圍，分別通過線上線下的方式計(jì)算了該項(xiàng)目區(qū)域的土方量。計(jì)算結(jié)果如表1所示，分別得到了線上線下項(xiàng)目場(chǎng)地的填挖方總量。其中PC 端為基于 線下土方計(jì)算軟件方法進(jìn)行計(jì)算，Web 端為基于WebGL 的土方在線自動(dòng)算量方法進(jìn)行計(jì)算。

表1 土方算量計(jì)算結(jié)果精確度對(duì)比

通過對(duì)比分析兩種方法的土方量計(jì)算結(jié)果，在像元大小為1 m×1 m 時(shí)，填挖總量相對(duì)誤差為0.07%和0.09%，差異值為0.02%，結(jié)算結(jié)果相近，誤差在可接受范圍內(nèi)，計(jì)算結(jié)果可靠，能夠滿足施工要求。

4 結(jié)語

通過實(shí)際工程應(yīng)用可以得到：基于WebGL 的在線自動(dòng)算量方法的研究在精準(zhǔn)度、計(jì)算效率、可視化、易獲取等關(guān)鍵特征取得了預(yù)期的結(jié)果，能夠?yàn)楹罄m(xù)土方程中提供精細(xì)化算量與方案優(yōu)化服務(wù)。多線程計(jì)算、數(shù)據(jù)輕量化、數(shù)據(jù)可視等技術(shù)的應(yīng)用，明顯提高了土方算量的計(jì)算效率，降低了土方工程成本，在提供可視化動(dòng)態(tài)展示的同時(shí)，能夠保證方法的易用性，對(duì)改變目前業(yè)主與施工單位信息不對(duì)稱產(chǎn)生諸多問題的現(xiàn)狀，推動(dòng)數(shù)字化應(yīng)用與管理，使其由粗放型管理與施工向精細(xì)化管理與施工邁進(jìn)具有一定建設(shè)意義。當(dāng)然，該計(jì)算方法也有一定的局限性，一是其計(jì)算精度主要取決于點(diǎn)云數(shù)據(jù)精度；二是當(dāng)前高精度點(diǎn)云數(shù)據(jù)采集成本相對(duì)較高，但我們相信隨著新一代測(cè)繪技術(shù)與信息技術(shù)的進(jìn)步，該方法一定能夠在工程建設(shè)項(xiàng)目中得到推廣應(yīng)用，發(fā)揮價(jià)值。

表2 土方算量計(jì)算效能對(duì)比