24小时热门版块排行榜    

查看: 2832  |  回复: 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个数据点对)?
其实,我常常想自己编程计算一下来验证这两个软件的结果,但这也是极为费力的事情。不知道在此启开这个话题,大家一起讨论后,能否让大家真正的使用这两个软件进行正确的半变异函数计算和选择?

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

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

糯糯鸽子

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
请问解决了吗?也想了解这个
2楼2020-06-10 09:51:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lb860 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +14 爱看书的可乐 2026-08-26 16/800 2026-08-26 18:28 by skyjudy
[基金申请] 我不理解! +9 Edward_pc 2026-08-26 14/700 2026-08-26 18:03 by hsding2026
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +6 process2012 2026-08-23 9/450 2026-08-26 15:33 by process2012
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 只有每年这种时候来逛逛小木虫 +26 yaoyewhu2008 2026-08-20 28/1400 2026-08-26 13:49 by zzuzxg
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见