(完整)克里金插值法

1、3)克里金插值法克里金插值法又稱空間局部插值法,是以變異函數(shù)理論和結(jié)構(gòu)分析為基礎(chǔ),在有限區(qū)域內(nèi)對區(qū)域化變量進(jìn)行無偏最優(yōu)估計(jì)的一種方法，是地統(tǒng)計(jì)學(xué)的主要內(nèi)容之一，由南非礦產(chǎn)工程師D.Matheron于1951年在尋找金礦時(shí)首次提出，法國著名統(tǒng)計(jì)學(xué)家G.Matheron隨后將該方法理論化、系統(tǒng)化，并命名為Kriging,即克里金插值法。1克里金插值法原理克里金插值法的適用范圍為區(qū)域化變量存在空間相關(guān)性，即如果變異函數(shù)和結(jié)構(gòu)分析的結(jié)果表明區(qū)域化變量存在空間相關(guān)性，則可以利用克里金插值法進(jìn)行內(nèi)插或外推。其實(shí)質(zhì)是利用區(qū)域化變量的原始數(shù)據(jù)和變異函數(shù)的結(jié)構(gòu)特點(diǎn)，對未知樣點(diǎn)進(jìn)行線性無偏、最優(yōu)估計(jì),無偏是指偏差

2、的數(shù)學(xué)期望為0，最優(yōu)是指估計(jì)值與實(shí)際值之差的平方和最小1.因此，克里金插值法是根據(jù)未知樣點(diǎn)有限領(lǐng)域內(nèi)的若干已知樣本點(diǎn)數(shù)據(jù),在考慮了樣本點(diǎn)的形狀、大小和空間方位，與未知樣點(diǎn)的相互空間關(guān)系,以及變異函數(shù)提供的結(jié)構(gòu)信息之后，對未知樣點(diǎn)進(jìn)行的一種線性無偏最優(yōu)估計(jì)。假設(shè)研究區(qū)域a上研究變量Z(x),在點(diǎn)xgA(i=1,2，n)處屬性值為Z(x),則待插點(diǎn)xeAii0處的屬性值Z(x)的克里金插值結(jié)果Z(x)是已知采樣點(diǎn)屬性值Z(x)(i=1,2，n)的加權(quán)和，00i即:Z*(x)二工九Z(x)(1)0iii=1式中是待定權(quán)重系數(shù).其中Z(x)之間存在一定的相關(guān)關(guān)系，這種相關(guān)性除與距離有關(guān)外，還與其相對方

3、向變化有關(guān)，克里金i插值方法將研究的對象稱“區(qū)域化變量針對克里金方法無偏、最小方差條件可得到無偏條件可得待定權(quán)系數(shù)(i=1,2,，n)滿足關(guān)系式：(2)工九=1ii=12，n)以無偏為前提，kriging方差為最小可得到求解待定權(quán)系數(shù)i的方程組:9X)+P=C(X9X)(j=1,2,j0j工九=1ii=1式中，C(x,x)是Z(x)和乙(x)的協(xié)方差函數(shù).ijij2方法步驟(完整)克里金插值法克里金插值法的應(yīng)用步驟如下:1、輸入原始數(shù)據(jù)，即采樣點(diǎn)，下面以輸入三個(gè)采樣點(diǎn)求待估插值為例來進(jìn)行說明。如圖1所示:2、網(wǎng)格化，選擇區(qū)域的范圍和網(wǎng)格的大小，對區(qū)域進(jìn)行網(wǎng)格化處理。3、數(shù)據(jù)檢驗(yàn)與分析，根據(jù)采樣

4、值是否合乎實(shí)際情況，剔除明顯差異點(diǎn)。4、直方圖的計(jì)算，直方圖有助于掌握區(qū)域變化的分布規(guī)律，以便決定是否對原始數(shù)據(jù)進(jìn)行轉(zhuǎn)換。5、利用變異函數(shù)進(jìn)行變異函數(shù)計(jì)算，了解變量的空間結(jié)構(gòu)。6、克里金插值估計(jì)待估點(diǎn)權(quán)重系數(shù)估計(jì)利用多邊形估計(jì)的方法，首先確定離待估點(diǎn)最近的采樣點(diǎn)的權(quán)重，根據(jù)公式(4)進(jìn)行采樣點(diǎn)權(quán)重估計(jì):(4)c+dwi=1根據(jù)搜索策略選擇合適的參估點(diǎn)，如圖2：圖2參估點(diǎn)圖示根據(jù)已經(jīng)求出的變異函數(shù)以及采樣點(diǎn)數(shù)量，三個(gè)采樣點(diǎn)列出三個(gè)等式，求出方程組的系數(shù)，公式為:C(1,1)C(1,2)C(1,3)1-C(0,1)C(2,1)C(2,2)C(2,3)1九2=C(0,2)C(3,1)C(3,2)C(

5、3,3)九L3C(0,3)_(4)分析在各向同性條件下改變塊金值與在塊金值相同條件下改變各向異性對權(quán)重值的影響。各向同性條件下改變塊金值時(shí)對權(quán)重值的影響效果如圖3(a)，在塊金值相同條件下改變各向異性對權(quán)重值帶來的影響如圖3(b):(a)(b)圖3各向同性條件下改變塊金值與在塊金值相同條件下改變各向異性對權(quán)重值的影響(5)根據(jù)求出的權(quán)重值，代入公式(1)，即可求得評估領(lǐng)域內(nèi)n個(gè)采樣值的線性組合2克里金插值法的方法路線圖如下：圖4方法路線圖3克里金插值法分類及適用類型克里金插值法主要有以下幾種類型：普通克里金(OrdinaryKriging)、簡單克里金(SimpleKriging)、泛克里金(

6、UniversalKriging)、協(xié)同克里金(Co-Kriging)、對數(shù)正態(tài)克里金(LogisticNormalKriging)、指示克里金(IndicatorKriging)、概率克里金(ProbabilityKriging)和析取克里金(DisjunctiveKriging)等1.克里金插值法可以簡單地表達(dá)為：Z(s)=p(s)+e(s)(6)式中，s為不同位置的點(diǎn)，可以人為是用經(jīng)緯度表示的空間坐標(biāo)；Z(s)為$處的變量值，它可以分解為確定趨勢值卩(s)和自相關(guān)隨機(jī)誤差e(s)。通過對這個(gè)公式進(jìn)行變化，可以生成克里金插值法的不同類型。首先，對于趨勢值卩(s)，可以簡單地賦予一個(gè)常量，即

7、在任何位置s處卩(s)=卩，如果卩是未知的，這便是普通克里金基本模型；卩(s)也可表示為空間坐標(biāo)的線性函數(shù)，如：卩(s)=P+Px+Py+Px2+Py2+(3xy(7)012345如果趨勢面方程中的回歸系數(shù)是未知的，則形成泛克里金模型；如果在任何時(shí)候趨勢已知的(如所有系數(shù)和協(xié)方差均已知)，無論趨勢常量與否，都會形成簡單克里金模型。其次，無論趨勢如何復(fù)雜，卩(s)仍無法獲得很好的預(yù)測，在這種情況下需要對誤差項(xiàng)e(s)進(jìn)行一些假設(shè),即假設(shè)誤差項(xiàng)e(s)的期望均值為0,且e(s)和e(s+h)之間的自相關(guān)不取決于s點(diǎn)的位置，而取決于位移量h。為了確保自相關(guān)方程有解，必須允許某兩點(diǎn)間自相關(guān)可以相等。然

8、后，可以對方程式左邊Z(s)進(jìn)行變換.例如，可以將其轉(zhuǎn)換成指示變量，即如果Z(s)低于一定的閾值，則將其值轉(zhuǎn)換為0，將高于閾值的部分轉(zhuǎn)換為1，然后對高于閾值部分作出預(yù)測，基于此模型作出預(yù)測便形成了指示克里金模型。如果將指示值轉(zhuǎn)變成含有變量的函數(shù)f(Z(s)，即形成析取克里金的指示函數(shù).最后，如果有多個(gè)變量的情況，則模型為：Z(s)=卩(s)+e(s)，其中j表示第j個(gè)變量。除了為每個(gè)jjj變量考慮不同的趨勢卩(s)外，隨機(jī)誤差e(s)之間還存在交叉相關(guān)性。這種基于多個(gè)變量的克里金模型即為jj協(xié)同克里金模型。不同的方法有其適用的條件，當(dāng)數(shù)據(jù)不服從正態(tài)分布時(shí)，若服從對數(shù)正態(tài)分布，則選用對數(shù)正態(tài)克里

9、金；若不服從簡單分布時(shí)，選用析取克里金；當(dāng)數(shù)據(jù)存在主導(dǎo)趨勢時(shí)，選用泛克里金；當(dāng)只需要了解屬性值是否超過某一閾值時(shí)，選用指示克里金；當(dāng)同一事物的兩種屬性存在相關(guān)關(guān)系時(shí)，且一種屬性不易獲取時(shí)，選用協(xié)同克里金，借助另一屬性實(shí)現(xiàn)該屬性的空間內(nèi)插；當(dāng)假設(shè)屬性值的期望值為某一已知常數(shù)時(shí)，選用簡單克里金；當(dāng)假設(shè)屬性值的期望值是未知的，選用普通克里金。（完整）克里金插值法4國內(nèi)外研究進(jìn)展從克里金方法被提出到現(xiàn)在已有完善的理論，并在很多領(lǐng)域得到了實(shí)際的應(yīng)用，在某些領(lǐng)域的應(yīng)用又推動了克里金理論的發(fā)展3。它的發(fā)展可歸納為四個(gè)時(shí)期，每個(gè)時(shí)期都是以每一屆地質(zhì)統(tǒng)計(jì)學(xué)大會的召開為標(biāo)志。第一時(shí)期，初次提出了地質(zhì)統(tǒng)計(jì)學(xué)理論，將

10、地質(zhì)統(tǒng)計(jì)學(xué)與傳統(tǒng)的統(tǒng)計(jì)學(xué)分開，且提出了區(qū)域化變量、簡單克里金、普通克里金、泛克里金的概念。第二時(shí)期，地質(zhì)統(tǒng)計(jì)學(xué)的理論逐步的幵始改進(jìn)和完善.第三時(shí)期，地質(zhì)統(tǒng)計(jì)學(xué)克里金在實(shí)踐應(yīng)用的發(fā)展相對理論發(fā)展更快，形成了兩種類型的理論體系：一類是有參數(shù)的克里金方法，另一類是沒有參數(shù)的克里金方法,有參數(shù)的克里金方法是指所研究的數(shù)據(jù)必須符合正態(tài)分布，如析取克里金；而沒有參數(shù)的克里金方法對所研究的變量的分布沒有特殊要求，如指示克里金和概率克里金.第四時(shí)期,克里金方法的應(yīng)用領(lǐng)域不斷擴(kuò)展壯大，在研究中有很多新的課題產(chǎn)生，克里金所研究對象已經(jīng)不再局限于空間領(lǐng)域的變量，隨著某些領(lǐng)域的需求，正在向時(shí)間-空間領(lǐng)域擴(kuò)展4。從目前

11、來看，克里金技術(shù)的發(fā)展可以概括如下：（1）形成了一套完整的理論體系。線性平穩(wěn)地質(zhì)統(tǒng)計(jì)學(xué)是地質(zhì)統(tǒng)計(jì)學(xué)的基礎(chǔ)部分，包含基本概念：區(qū)域化變量理論;基本工具：變差函數(shù);基本假設(shè):二階平穩(wěn)假設(shè)和本征假設(shè)；基本公式：估計(jì)反差和普通克里金法；線性非平穩(wěn)地質(zhì)統(tǒng)計(jì)學(xué)包括了泛克里金和K階本征函數(shù)法等。平穩(wěn)非線性地質(zhì)統(tǒng)計(jì)學(xué)包含析取克里金等。（2）編制了一些實(shí)際有效的程序以及軟件。例如斯坦福大學(xué)的GeostatisticalEarthModelingSoftware.（3）地質(zhì)統(tǒng)計(jì)學(xué)的提出原本是為了解決礦產(chǎn)儲量的估計(jì)，但是隨著地質(zhì)統(tǒng)計(jì)學(xué)的發(fā)展，人們發(fā)現(xiàn)其研究對象存在于很多種自然現(xiàn)象中.于是，地質(zhì)統(tǒng)計(jì)學(xué)不再是研究地質(zhì)領(lǐng)

12、域的特有方法，而成為研究某類自然現(xiàn)象通用的方法，例如降水量的分布、水文層的滲透率和孔隙度等屬性值、在醫(yī)學(xué)上對骨豁的三維重建5等等。目前國內(nèi)外學(xué)者利用克里金插值法做了大量研究。翟進(jìn)乾應(yīng)用克里金插值方法對煤層分布監(jiān)測進(jìn)行了系統(tǒng)分析研究6；張蕾、陳曉宏將克里金插值方法用于珠江三角洲網(wǎng)河區(qū)水位空間插值7;尚慶生、郭建文等將克里金插值方法用于計(jì)算青藏鐵路鉆孔地溫?cái)?shù)據(jù),實(shí)現(xiàn)了數(shù)據(jù)的體視化8；顏輝武，祝國瑞等采用克里金插值方法建立水文地質(zhì)層三維模型9，并利用體繪制技術(shù)進(jìn)行可視化表達(dá)，取得了良好的效果;劉承香、阮雙深、伍小芹提出基于克里金插值方法進(jìn)行水深數(shù)據(jù)插值形成規(guī)則網(wǎng)格數(shù)字高程模型的算法，對海底數(shù)字地圖的

13、模擬具有重要參考價(jià)值，數(shù)字仿真結(jié)果證明該算法可行10。參考文獻(xiàn)：1湯國安，楊昕.ArcGIS地理信息系統(tǒng)空間分析實(shí)驗(yàn)教程M。北京：科學(xué)出版社，2011。2孟俊貞?？死锝鸩逯到凭W(wǎng)格算法在柵格數(shù)據(jù)投影變換中的應(yīng)用D。長沙：中南大學(xué)，2009。3曲壽利，王鑫國內(nèi)外物探技術(shù)現(xiàn)狀與展望M.石油工業(yè)出版社，2003.4姚興苗快速三維克里金插值方法研究及實(shí)現(xiàn)D。成都：電子科技大學(xué)，2013。5胡巖，王田苗，王君臣?；贙riging算法的手術(shù)導(dǎo)航三維形變技術(shù)J北京航空航天大學(xué)學(xué)報(bào)，2010,5:12.6翟進(jìn)乾?？死锝穑╧riging）插值方法在煤層分布檢測中的應(yīng)用研究D。太原：太原理工大學(xué)，2008.(完整)克里金插值法7張蕾，陳曉宏珠江三角洲網(wǎng)河區(qū)水位空間插值的kriging方法J。中山大學(xué)學(xué)報(bào)(自然科學(xué)版