【求助】Geostaticstial 中kriging插值问题（一系列，很多人都想知道），高手都进来指 - 地学 - 小木虫 - 学术 科研 第一站
&& 查看话题
【求助】Geostaticstial 中kriging插值问题（一系列，很多人都想知道），高手都进来指
首先说明一点，我也有这样的问题，但这个贴不是本人原创，转自http://bbs./ESRI/viewthread.php?tid=48633&highlight=kriging
看了不少资料，还是很模糊，我把插值主要过程写下来，同时附上一些问题，麻烦高手都来指点一二，不胜感激。我想很多新手朋友也能从中学到一些东西。
1.先对数据检验，正态分布，很好。
2.做趋势分析，发现南北就方向有呈U形，不过不能完全确定是否真存在趋势，按教程上说投影点在投影面上要比较集中，我的图中投影点在南北方向并不是很集中。我也围换过location，无论在那个角度，投影点都比较分散。
3.空间自相关及方向变异查看，采用半变异和协方差云图。说实在的，没有看出多少东东来，感觉变异分布关于第一象限45度对称，是不是与趋势分析中不一至？
从教程中的例子，它能从半变异和协方差云图中看到变程值为16，的确在h为16左右时，rh比较稳定，并且h越小，rh也越小，相关性越强；但我的数据并不是这种趋势，看不到变程h，很分散，相反，与教程中所说没有相关性的那一个图更符合，是不是我的数据没有很好的相关性？注意此处的lag size为1549，与后面插值的不一样。
4.暂按有趋势的进行插值，按教程里说数据存在主导趋势时，采样泛克立格，不知道这里选择是否合理。趋势选择2次。
5.插值的第2步，不知道这里的参数要不要修改？我一般都是采用默认。
6.插值第3步，此处我的做法是：如果存在各向异性，则打开Anisotropy选项，至于参数，让系统自动生成。从这里的半变异图可以看到，系统生成的变程a在2041左右，从图上也可以确定在这个位置趋于稳定。间于这是插值最重要的一步，我有三个问题。
& & 1）.块金值Nugget、偏基台值Partial Sill、变程Major range是否要改变？如果要改变，如何设置？如果从前面第3步能看得到变程，在这里就应该得设置变程？
& & 2）.lag size和number of lags 是否需要修改，如何修改，原由？
& & 3）.模型的选择，一般是球面模型，但也有用指数模型及其它的。请问这里的选择是根据数据就可以确定呢还是采样不同模型进行插值比较最终取的最优的那一个模型？
7.插值第4步，这里的问题是图中圈的那部分，就是包含的点数是否采用默认？还是根据当前样本点数有一个经验性的数据？另外sector type的选择，我也没有弄明白。如果前面第6步选择了各向异性，在该页面椭圆的参数也是自动生成，没有问题吧？
至于后面的步骤就不再说了，看验证指标Mean,Root-Mean-Square等就可以了。选择插值后验证指标最好的一个。
&&发这个帖花了不少时间，小弟刚入行不久，很多东西不了解，希望高手朋友们指点指点，发表一些实质性的意见，讲讲你们插值的经验，非常感谢！
& & 谢谢大家能帮忙！:tiger28:本人BB不多了，请凉解！
这个是ArcGis的统计分析？？ Originally posted by crane_yi at
这个是ArcGis的统计分析？？ 嗯，是 我也是做了这个&&但其实也没看出什么东西啦&&
就选个把它一步一步点下去 生成了插值图&&
所以也很矛盾&&但现在也快来不及了&&就这样了 Originally posted by chenhb13 at
我也是做了这个&&但其实也没看出什么东西啦&&
就选个把它一步一步点下去 生成了插值图&&
所以也很矛盾&&但现在也快来不及了&&就这样了 你没有这样的感觉吗？：一个东西可以用很多方法做，也就是说参数不同，不同方法之间的比较，谁好就用谁的问题，还有就是，只靠倒数第二步那几个指标来比较吗，很茫然，，:shuai: Originally posted by duhb655 at
你没有这样的感觉吗？：一个东西可以用很多方法做，也就是说参数不同，不同方法之间的比较，谁好就用谁的问题，还有就是，只靠倒数第二步那几个指标来比较吗，很茫然，，:shuai: 对啊，是可以通过很多方法来比较，比如说你可以比较不同模型之间，不同阶趋势等等，但我做的就是找一个最佳的，不过是不是最佳呢，其实我真的不懂
现在也没多少时间了 ，就想混个毕业的了 也不想太去深究了
觉得主要的问题是在变程的设置上，你有没看过很多文章都没把变程说得很明白 Originally posted by chenhb13 at
对啊，是可以通过很多方法来比较，比如说你可以比较不同模型之间，不同阶趋势等等，但我做的就是找一个最佳的，不过是不是最佳呢，其实我真的不懂
现在也没多少时间了 ，就想混个毕业的了 也不想太去深究了
... 我看了个，是设置lag size,而这个通过spatial analysis tools实现，range通过lag size就生成了，我现在就是在趋势分析这块有点糊涂还，不知道怎么样，或者说定哪个方向才是最好的，我也看了个贴子，也用spatial analysis tools确定区域内数据（weight）的方向，但我试了一下，一个区域内不同的数，用这个做出来的方向居然都一样。。。。不可信，真是悲剧，难道geostatistics 这块没人会啊，一直等着神人出现，:sweat: help help F1 这一块我也遇到了同样的问题 有没有高人帮忙解决一下&&万分感激 看帮助文档吧，帮助文档上说的很清楚。 这个就是要根据经验和数据以及最后的结果来选模型。。。看看王政权的那本地统计在生态学中的应用 小弟是即将步入硕士的娃娃，本科毕设就是关于土壤微量元素的空间变异分析的。关于楼主的问题可以通过GS+9.0这个软件计算，步长、块金值、基台值等等，然后在ArcGIS中输入对应参数，即可。 Originally posted by 湘雨儿1 at
小弟是即将步入硕士的娃娃，本科毕设就是关于土壤微量元素的空间变异分析的。关于楼主的问题可以通过GS+9.0这个软件计算，步长、块金值、基台值等等，然后在ArcGIS中输入对应参数，即可。 没用过GS+9这个软件，它是什么呢？关键是步长和数组吧？？！！这两个确定了，块金值和基台值就确定了。而步长可以在ARCGIS里面算出来。而且，当具有各向异性时，方向的选择呢？就是趋势方向，这个在ARCGIS里也可以算出一个数据的方向，但不知道对不对。 Originally posted by duhb655 at
没用过GS+9这个软件，它是什么呢？关键是步长和数组吧？？！！这两个确定了，块金值和基台值就确定了。而步长可以在ARCGIS里面算出来。而且，当具有各向异性时，方向的选择呢？就是趋势方向，这个在ARCGIS里也可 ... GS+&&是一款挺好用的半方差变异分析的软件，你可以到网上去搜，或者给我发给你。
有些学者认为ArcGIS地统计模块受数据计算的限制，所以会采用GS软件进行变异函数的计算和理论模型的拟合。
所以，块金值和基台值两款软件计算出来的是不一样的。各向同性还是各向异性我个人觉得看看最后拟合精度选择吧~~~ Originally posted by 湘雨儿1 at
GS+&&是一款挺好用的半方差变异分析的软件，你可以到网上去搜，或者给我发给你。
有些学者认为ArcGIS地统计模块受数据计算的限制，所以会采用GS软件进行变异函数的计算和理论模型的拟合。
所以，块金值和 ... 谢谢你了，有的话而且方便的话就发我邮箱吧：
想再表达一下，那个各向异性应该是在数据检查或数据趋势处理时就能看出来吧，或根据数据在研究区的特性就能确定是，比如一个地区的臭氧浓度，如果当地存在主导风向，那肯定得选向异性了，接下来就是确定方向问题了，不知道我这种理解对不对，不管怎样，谢谢你了！ Originally posted by duhb655 at
谢谢你了，有的话而且方便的话就发我邮箱吧：
想再表达一下，那个各向异性应该是在数据检查或数据趋势处理时就能看出来吧，或根据数据在研究区的特性就能确定是，比如一个地区的 ... http://emuch.net/bbs/viewthread.php?tid=1082815&&这个帖子有链接&&呵呵 给你发给破解文件吧 就是复制到安装目录下运行即可&&
你这样说是对的，我对ArcGIS里的Trend Analyst怎么用还不太了解。 这个插值实际很有争议，边界对插值也会有很大的影响。 Originally posted by crane_yi at
这个是ArcGis的统计分析？？ ArcGis几:o？ Originally posted by muliptpt at
ArcGis几:o？ 9.2就有，不知道从哪个版开始的， ArcGIS的帮助中写的也很含糊不清，这一块关于参数的设定，大家的意见差异很大
var cpro_id = 'u1216994';
欢迎监督和反馈：本帖内容由
提供，小木虫为个人免费站点，仅提供交流平台，不对该内容负责。欢迎协助我们监督管理，共同维护互联网健康，如果您对该内容有异议，请立即发邮件到
联系通知管理员，也可以通过QQ周知，我们的QQ号为：8835100
我们保证在1个工作日内给予处理和答复，谢谢您的监督。
小木虫，学术科研第一站，为中国学术科研研究提供免费动力
欢迎监督，发现不妥请立即
E-mail： & QQ：8835100}