淺析ArcGIS中的地統(tǒng)計克里格插值法及其應(yīng)用

劉青豪+高浩原+李練武

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

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

1.引言

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

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

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

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

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

3. ArcGIS地統(tǒng)計

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

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

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

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

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

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

（2）預(yù)測未知點的值

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

4. 土方量計算

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

5. 結(jié)語

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

參考文獻(xiàn)：

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

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

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

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