劉青豪+高浩原+李練武

【摘要】Arc GIS軟件中設(shè)立的地統(tǒng)計(jì)分析擴(kuò)展模塊應(yīng)用領(lǐng)域十分廣泛，在施工建設(shè)中發(fā)揮著顯著的作用。本文首先介紹了ArcGIS、地質(zhì)統(tǒng)計(jì)學(xué)的概念以及克里格插值（kriging）的各種方法，然后基于地統(tǒng)計(jì)的理論，運(yùn)用Arc GIS 軟件中地統(tǒng)計(jì)分析模塊，對(duì)Kriging插值法在土方計(jì)算中的應(yīng)用進(jìn)行了淺要的分析。

【關(guān)鍵詞】ArcGIS；地統(tǒng)計(jì)；土方量

1.引言

隨著地理信息技術(shù)在當(dāng)今社會(huì)的飛快發(fā)展，以及人類對(duì)空間數(shù)據(jù)質(zhì)量要求的提高，空間數(shù)據(jù)插值方法及其在實(shí)際中的應(yīng)用越來(lái)越受到人類的重視。地統(tǒng)計(jì)學(xué)里面的克里格插值方法，由于其具有插值和預(yù)測(cè)兩方面的優(yōu)勢(shì)，正在越來(lái)越多的領(lǐng)域中得到廣泛應(yīng)用，成為許多專業(yè)軟件、商業(yè)軟件的重要組成部分。社會(huì)飛速發(fā)展的今天，地統(tǒng)計(jì)的應(yīng)用也體現(xiàn)在生活的大小方面，本文以工程施工中的土方計(jì)算為例子做出了詳細(xì)的地統(tǒng)計(jì)分析。

2. 地質(zhì)統(tǒng)計(jì)學(xué)基礎(chǔ)

地質(zhì)統(tǒng)計(jì)學(xué)，又名地統(tǒng)計(jì)，是統(tǒng)計(jì)學(xué)與地質(zhì)學(xué)兩學(xué)科領(lǐng)域的交叉部分。經(jīng)過(guò)不斷完善和改革，目前已成為具有扎實(shí)理論根基和實(shí)用價(jià)值的數(shù)學(xué)工具。地統(tǒng)計(jì)運(yùn)用變異函數(shù)的理論，對(duì)自然現(xiàn)象進(jìn)行深入研究。研究對(duì)象既具有隨機(jī)性又具有結(jié)構(gòu)性，任何科學(xué)研究，凡是與空間數(shù)據(jù)的結(jié)構(gòu)性、隨機(jī)性等有關(guān)，皆可應(yīng)用地統(tǒng)計(jì)學(xué)的理論與方法。

2.1 ArcGIS 軟件介紹。ArcGIS誕生于Esri公司，是該公司集近50年地理信息系統(tǒng)咨詢和研發(fā)經(jīng)驗(yàn)的產(chǎn)物，下分多個(gè)模塊，具有強(qiáng)大的地圖制作、空間分析能力，并且在空間信息整合、發(fā)布與共享方面有著十足的優(yōu)勢(shì)。是一套完整的地理信息平臺(tái)產(chǎn)品。目前地理信息系統(tǒng)行業(yè)中，GIS是應(yīng)用最為廣泛、性能最為豐富的軟件，并且支持對(duì)地理信息進(jìn)行處理、歸類，從中提取生產(chǎn)所需的信息。

2.2 克里格插值法。克里格插值（Kriging），一般可理解為包括普通克里格方法、協(xié)同克里格法等在內(nèi)的眾多插值方法。其中，普通克里格應(yīng)用最為普遍，它以未知的常數(shù)作為假定均值。隨著Kriging與其它研究領(lǐng)域的相互交叉發(fā)展，在自身應(yīng)用發(fā)展的同時(shí)，也衍生了一些新的克里金方法。在克里格法的應(yīng)用過(guò)程中，Kriging考慮了空間屬性的各向異性，插值精度較高，分析結(jié)果可靠，對(duì)數(shù)據(jù)處理后的成果偏離度微乎其微。

3. ArcGIS地統(tǒng)計(jì)

隨著地理信息學(xué)科的發(fā)展，近年來(lái)，GIS 與地統(tǒng)計(jì)學(xué)之間的生產(chǎn)發(fā)展迅速，地統(tǒng)計(jì)分析模塊中產(chǎn)生的新表面，可以繼而使用于GIS建模以及諸如空間分析和3D分析等ArcGIS擴(kuò)展模塊的可視化當(dāng)中。地統(tǒng)計(jì)分析的地表擬合包括三個(gè)關(guān)鍵步驟：探索性空間數(shù)據(jù)分析、結(jié)構(gòu)分析（臨近點(diǎn)位空間屬性的計(jì)算與建模）、趨勢(shì)面預(yù)測(cè)和結(jié)果評(píng)價(jià)。運(yùn)用這些基本的分析處理模塊，不僅可以高效地完成多種地統(tǒng)計(jì)分析問(wèn)題，此外還可以創(chuàng)建形象的的專題地圖（表面預(yù)測(cè)）。

3.1 探索性空間數(shù)據(jù)分析。探索性空間數(shù)據(jù)分析模塊是目前發(fā)展應(yīng)用迅猛的一個(gè)領(lǐng)域，興起于20世紀(jì)后半葉是一種新的統(tǒng)計(jì)分析方法。在得到樣本點(diǎn)數(shù)據(jù)之后，首先我們需要做的是對(duì)樣本數(shù)據(jù)進(jìn)行分析檢驗(yàn)，從而分析數(shù)據(jù)的趨勢(shì)性。在ArcGIS中設(shè)立的探索性空間數(shù)據(jù)分析模塊，為操作者提供了一系列的工具來(lái)檢查數(shù)據(jù)，提高了解決問(wèn)題的高效性。

筆者采用了若干樣本高程點(diǎn)，并對(duì)其進(jìn)行了樣本分析與擬合表面（詳見下圖），對(duì)數(shù)據(jù)的分布狀況進(jìn)行詳細(xì)了解，不僅反應(yīng)著工作的嚴(yán)謹(jǐn)程度，更關(guān)系著輸出結(jié)果的準(zhǔn)確性，是運(yùn)用地統(tǒng)計(jì)分析建立表面之前的重要工作。Arc GIS 中的數(shù)據(jù)分析模塊，為人們提供了分布直方圖法和正態(tài)QQPlot 圖法兩種檢驗(yàn)數(shù)據(jù)的方法。兩種方法在高程點(diǎn)擬合過(guò)程中均有著很好的表現(xiàn)。下圖為對(duì)研究區(qū)內(nèi)的樣本高程點(diǎn)進(jìn)行分析及擬合表面的情況。

3.2 使用普通克里格法創(chuàng)建表面

普通克里格方法進(jìn)行表面內(nèi)插的方法簡(jiǎn)便，可分為兩步進(jìn)行，分別是：

（1）生成變異函數(shù)和協(xié)方差函數(shù)

（2）預(yù)測(cè)未知點(diǎn)的值

為了表達(dá)研究區(qū)域的原始地形，通過(guò)生成一個(gè)關(guān)于高程的克里格插值圖是最直觀的方法。如上兩步操作，便可以生成一幅基于Kriging里格插值法的地理信息預(yù)測(cè)圖。本文重點(diǎn)在于討論地統(tǒng)計(jì)在實(shí)際中的應(yīng)用，在此不過(guò)多闡述。

4. 土方量計(jì)算

隨著工程建設(shè)以及社會(huì)各個(gè)領(lǐng)域的發(fā)展，地統(tǒng)計(jì)在施工建設(shè)、環(huán)境調(diào)查、土地管理、農(nóng)業(yè)開發(fā)等等方面越來(lái)越凸顯出其影響。筆者以土方量計(jì)算為典型，介紹地統(tǒng)計(jì)在實(shí)際中的應(yīng)用。土方量的計(jì)算是建筑工程施工的一個(gè)重要步驟，工程施工前的設(shè)計(jì)階段必須對(duì)土石方量進(jìn)行估算， 方格網(wǎng)法、等高線法、DTM法和平均高程法等均為比較常用到的幾種計(jì)算土方量。ArcGIS計(jì)算土方量的原理簡(jiǎn)單卻很嚴(yán)謹(jǐn)，整個(gè)流程是以 DEM 為基礎(chǔ)，先后分別通過(guò)建模分析、區(qū)域求差、統(tǒng)計(jì)單區(qū)域土方量、計(jì)算工程土方量來(lái)完成最終工程需求，此種流程方法已在工程施工中得到了廣泛應(yīng)用。

5. 結(jié)語(yǔ)

計(jì)算土方量的方法非常之多，隨著工程施工技術(shù)的發(fā)展，無(wú)論在土方計(jì)算的精度還是效率方面都呈現(xiàn)逐步提高的形勢(shì)，然而，大量的實(shí)踐證明，基于 Arc GIS地統(tǒng)計(jì)分析模塊的計(jì)算方法，其計(jì)算與地形狀況無(wú)關(guān)，并且計(jì)算簡(jiǎn)單、高效，可以適用于各種地形項(xiàng)目。在今后的施工建設(shè)或其他領(lǐng)域的發(fā)展中，筆者堅(jiān)信地統(tǒng)計(jì)的應(yīng)用范圍會(huì)越來(lái)越廣。

參考文獻(xiàn)：

[1]王艷妮，謝金梅，郭祥. ArcGIS中的地統(tǒng)計(jì)克里格插值法及其應(yīng)用[J]. 軟件導(dǎo)刊，2008，12：36-38.

[2]陳勇. 利用ArcGIS地統(tǒng)計(jì)分析進(jìn)行土地平整土方量計(jì)算的研究[J]. 安徽農(nóng)業(yè)科學(xué)，2007，01：70-71+103.

[3]徐嬋，徐漢超，包德高. 基于ArcGIS的土方量計(jì)算[J]. 廣西水利水電，2013，01：81-83.

[4]王先鵬，曹榮林. 土方量計(jì)算的原理與方法及ArcGIS的應(yīng)用前景[J]. 地理空間信息，2009，04：139-141.