GS地統(tǒng)計(jì)和ARCGIS克里格插值過程

1、由于是初學(xué)地統(tǒng)計(jì)和克里格插值，現(xiàn)將自己處理數(shù)據(jù)的過程和步驟列出，中間有幾個(gè)問題很是迷惑，還請(qǐng)相關(guān)的專家們給點(diǎn)指導(dǎo)，或者同行們討論一下，對(duì)我處理的過程有什么不合理的地方，還請(qǐng)指出，謝謝！1、在GS+7.0中進(jìn)行地統(tǒng)計(jì)分析，將經(jīng)緯度坐標(biāo)轉(zhuǎn)換成平面坐標(biāo)，Z值為土壤鹽分?jǐn)?shù)據(jù)，導(dǎo)入到軟件中，重計(jì)算后如下圖1：2、查看數(shù)據(jù)的分布，發(fā)現(xiàn)進(jìn)行l(wèi)og變換后數(shù)據(jù)的分布狀態(tài)還不如元數(shù)據(jù)，所以未進(jìn)行變換（圖2）：3、接下來進(jìn)行半方差分析，初始界面如圖3：4、進(jìn)行計(jì)算（圖4）：5、查看模型信息，顯示最優(yōu)模型為高斯模型，以及各種參數(shù)，這里有點(diǎn)不明白的是那個(gè)RSS值怎么會(huì)那么大？圖5：6、然后再ARCGIS中進(jìn)行克里格插值

2、，初始界面如圖6：7、選擇普通克里格，數(shù)據(jù)不進(jìn)行變換，圖7：8、下一步，圖8顯示的是ARCGIS自動(dòng)給出的各個(gè)參數(shù)和模型9、根據(jù)GS+7.0中的參數(shù)對(duì)圖8進(jìn)行修改，修改后的界面為圖9，主要修改了塊金值、變程、偏基臺(tái)值、模型類型以及l(fā)ag size：10、下一步，沒做改動(dòng) 圖10：11、下一步顯示交互驗(yàn)證結(jié)果，圖11：12、最終的出圖顯示，圖12：我最后將交叉驗(yàn)證的屬性表導(dǎo)出來之后，計(jì)算各點(diǎn)的真實(shí)值和預(yù)測(cè)值的相關(guān)系數(shù)，僅為0.2多，這算是好嗎？還有就是，我分別按照GS+中給出的其他模型的參數(shù)輸入到ARCGIS插值過程中，最后得到的交叉驗(yàn)證結(jié)果為下圖13和14，圖13為指數(shù)模型，圖14為球狀模型，比較三者，發(fā)現(xiàn)指數(shù)模型的交叉驗(yàn)證結(jié)果最好，但是指數(shù)模型中真實(shí)值和預(yù)測(cè)值的相關(guān)系數(shù)僅為0.19啊，這都怎么回事啊？最后我用ARCGIS默認(rèn)的各個(gè)參數(shù)進(jìn)行插值，得出的交叉驗(yàn)證結(jié)果為圖15，比指數(shù)模型的效果差，而且相關(guān)系數(shù)為0.14，都是那么低啊。指數(shù)模型、球狀模型和ARCGIS默認(rèn)參數(shù)的最終效果為圖16、17、18就是這樣了，請(qǐng)大家積極討論啊，相互學(xué)習(xí)！1.jpg2.jpg3.jpg4.jpg5.jpg6.jpg7.jpg8.jpg9.jpg10.jpg11.jp