| 查看: 2667 | 回复: 1 | |||
lb860木虫 (正式写手)
|
[交流]
半变异函数的求解 已有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个数据点对）？ 其实，我常常想自己编程计算一下来验证这两个软件的结果，但这也是极为费力的事情。不知道在此启开这个话题，大家一起讨论后，能否让大家真正的使用这两个软件进行正确的半变异函数计算和选择？ 突然又想起一个问题，根据半变异函数理论，只有在一定范围内的数据才具有空间相关性，但我们在实际区域数据分析中，比如相关分析、因子分析等等，却并不考虑数据之间的距离，这不是矛盾了么？？？ |
回复此楼» 猜你喜欢
地球化学论文润色/翻译怎么收费?
已经有86人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 物理模型中的复数解
已经有12人回复
- 求解高斯函数和超高斯函数的功率谱密度函数
已经有3人回复
- GS+获得的变异函数参数为何在Arcgis中表现那么差，如何选择Arcgis的克里格插值差数呢
已经有7人回复
- 一个高一函数题求解
已经有7人回复
- 求助高斯拟合的问题
已经有7人回复
- 求解导函数的连续性问题
已经有24人回复
- 请教函数求解问题
已经有11人回复
- matlab 隐函数求解
已经有6人回复
- 如何选择数据拟合的模型和样本点？
已经有3人回复
- 急救！…… 能量传递中的，（H-I）模型的拟合问题
已经有9人回复
- 求解二元原函数
已经有3人回复
- 半方差函数是否就是半变异函数
已经有5人回复
- GS+和Arcgis变异函数拟合
已经有9人回复
- 指数函数增长率求出来有偏差，找不到原因。。。。
已经有15人回复
- 关于GIS的差值问题，克里金差值和坡度模拟，请学长们帮忙支支招！
已经有8人回复
- 二元幂函数参数求解
已经有10人回复
- 拟合模型的好坏
已经有4人回复
- 求助matlab隐函数参数拟合
已经有16人回复
- 【求助】求教:如何用SURFER 求半方差和空间趋势分析
已经有6人回复
- 【求助】求解隐函数
已经有6人回复
- 【求助】Arcgis地统计学半变异函数模型拟合
已经有7人回复
- 【求助】贝塞尔函数求解
已经有4人回复
- 【求助】推求“稳健变异函数”可通过什么软件实现？
已经有4人回复
- 【求助】二元函数拟合
已经有12人回复
- 【求助】求在arcgis中做插值
已经有9人回复
- 【求助】递推公式的参数拟合
已经有5人回复
2楼2020-06-10 09:51:26
