24小时热门版块排行榜

返回列表

lb860

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 2519.6
散金: 4
红花: 1
帖子: 382
在线: 77.6小时
虫号: 845011
注册: 2009-09-10
性别: GG
专业: 实验地球化学和计算地球化

[交流] 半变异函数的求解已有1人参与

克里金差值首先需要求取半变异函数，它是矢量距离h的函数，但这个问题似乎一直是大家纠结的问题，我也很纠结。
实际工作中，采样点位并未位于正规网格节点上，甚至较为离散，所以在计算半变异函数值时，要考虑角度容差和距离容差；也就是说，在理论上，x+h数据是足够的，但实际上，x+h数据极少，因此必须考虑容差。在矢量h的角度容差和距离容差范围内，都可以看做是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，surfer选择的是对角线距离最大值的1/3。但这个数值我个人认为影响不大（只要不是太离谱），它影响的仅仅是点对数的多少（因为在实际工作中，各自距离的1/2和1/3都应该超出了样品的相关性范围）。不过对于搜索半径，我也看到一些资料说选择采样间隔的2.5倍到3倍。
（3）各向异性的整体考虑。GS+中，在半变异函数计算中并未整体考虑各向异性（我个人认为，不知道是否对），而surfer考虑了，但是surfer中的自动拟合参数似乎有些问题；而且，模型得自己选择并进行比较得出最优结果，而GS+默认选择的已经是最优的。
不知道上述观点大家是否同意？大家一起讨论讨论。此外，我个人还有四点不懂：
（1）两者的距离容差是怎么规定的？
（2）surfer中的Angular Divisions和Radial Divisions具体代表什么意义？
（3）surfer中模型选择中似乎块金效应和常用的直线、球状、指数三种模型是并列的？
（4）GS+中，h的选择是计算间隔的倍数，surfer中h是如何选择的？是最大滞后距离的25份（surfer默认好像是25个数据点对）？
其实，我常常想自己编程计算一下来验证这两个软件的结果，但这也是极为费力的事情。不知道在此启开这个话题，大家一起讨论后，能否让大家真正的使用这两个软件进行正确的半变异函数计算和选择？

突然又想起一个问题，根据半变异函数理论，只有在一定范围内的数据才具有空间相关性，但我们在实际区域数据分析中，比如相关分析、因子分析等等，却并不考虑数据之间的距离，这不是矛盾了么？？？

回复此楼