克里金差值首先需要求取半变异函数，它是矢量距离h的函数，但这个问题似乎一直是大家纠结的问题，我也很纠结。实际工作中，采样点位并未位于正规网格节点上，甚至较为离散，所以在计算半变异函数值时，要考虑角度容差和距离容差；也就是说，在理论上，x+h数据是足够的，但实际上，x+h数据极少，因此必须考虑容差。在矢量h的角度容差和距离容差范围内，都可以看做是2009上非诚勿扰
x+h，这样才能计算半变异函数值。在半变异函数的求解中，最方便又常用的软件就是GS+和Surfer（不要提ArgGIS），两者区别在哪？个人认为主要在以下三个方面：（1）容差。我们知道，在看各向异性时，一般都是以0度（即x轴正向）为始，45度为间隔，看8斯蒂文斯皮尔伯格个方向上的各向异性。在GS+中，默认角度容差为22.5度，这个数字化刚刚好（这个容易理解），而Surfer中默认为90度，那也就是说surfer中考虑各向异性仅仅考虑x轴正向和x轴负向两个方向，当然这个似乎可以改变。（2）距离选择。GS+中有两个距离，一个是最大滞后距离，一个是计算间隔，其中计算间隔才是决定半变异函数模型的主要参数；surfer中只有一个，是最大滞后距离。最大滞后距离（是否也就是搜索半径呢？我个人认为是），GS+选择的是x、y轴两者最大距离的1/2，su rfer选择的是对角线距离最大值的1/3。但这个数值我个人认为影响不大（只要不是太离谱），它影响的仅仅是点对数的多少（因为在实际工作中，各自距离的1/2和1/3都应该超出了样品的相关性范围）。不过对于搜索半径，我也看到一些资料说选择采样间隔的2.5倍到3倍。
（3）各向异性的整体考虑。GS+中，在半变异函数计算中并未整体考虑各向异性（我个人认为，不知道是否对），而surfer考虑了，但是surfer中的自动拟合参数似乎有些问题；而且，模型得自己选择并进行比较得出最优结果，而GS+默认选择的已经是最优的。
不知道上述观点大家是否同意？大家一起讨论讨论。此外，我个人还有四点不懂：
（1）两者的距离容差是怎么规定的？
（2）msn shellsurfer中的Angular Divisions和Radial Divisions具体代表什么意义？
（3）surfer中模型选择中似乎块金效应和常用的直线、球状、指数三种模型是并列的？
（4张博树）GS+中，h的选择是计算间隔的倍数，surfer中h是如何选择的？是最大滞后距离的25份（surfer默认好像是25个数据点对）？
其实，我常常想自己编程计算一下来验证这两个软件的结果，但这也是极为费力的事情。不知道在此启开这个话题，大家一起讨论后，能否让大家真正的使用这两个软件进行正确
的半变异函数计算和选择？
突然又想起一个问题，根据半变异函数理论，只有在一定范围内的数据才具有空间相关性，但我们在实际区域数据分析中，比如相关分析、因子分析等等，却并不考虑数据之间的距离，这不是矛盾了么？？？
在GS+里对变异函数进行拟合，选择出最优的变异函数。然后在Arcgis里进行克里格插值，比较了不同的变异函数，进行交叉验证。验证结果发现，GS+里得到的最优的变异函数，交叉验证的结果并不是最好的。这个时候该要怎么取舍和解决呢？
拿来主义课堂实录
3、如何将GS+软件作出的半变异函数图表通过输出数据用Excel做出来？？？急！ 已有1人参与 | |
| |
初始悬赏金币 5 个
由于要使用ArcGIS的Ordinary Kriging插值，所以需要先用GS+获得相应的参数。
在GS+做成半变异函数模型后，因为软件给出的图实在是太丑了（请原谅我没事事）。。然后想用Excel做一下
初步探索后，可以点击list导出一个average distance，average semivariance和pair的一组数据点，然后在Excel中以average distance为X，average semivariance为Y 就得到experimental semivariogram（一些点）
然后关键就是怎么把模型的那条线弄出来了，以X=100,20000 然后根据得出的最佳模型（线性、指数、高斯、球面），按照它的公式计算出相应的y值，然后这些密密麻麻的点就形成了一条线
但是我用Excel作出的图和GS+生物学杂志本身给出的图不大一样，主要就是Excel作出的图点和线的贴合度不是很好，那些试验点离那条线比较远，求大神指导到底怎么做，这样做对不对，Cd和Cr的图和原图差不多，Pb和Cu差的比较多
回复：可以用ISATIS软件做变差函数图。
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议