24小时热门版块排行榜    

查看: 3309  |  回复: 13

孔子无敌老子

金虫 (小有名气)

[交流] 关于学习地质统计学的一些建议 已有11人参与

最近收到很多虫友关于地质统计学方面的咨询,也都逐个回复了,但是还是觉得非常有必要单开一个帖子说点什么。其实我学习地质统计学的时间也不是很长,到目前为止才两年不到,所以说不上懂行,如果下面有说的不对的地方还请各位批评指正。
相信很多虫友接触到地质统计学都是因为要对研究对象进行空间变异分析,或者克里格插值,其中最为关键的应该就是变异函数(又称半方差函数、半变异函数)拟合了。相信很多朋友也发现了,介绍变异函数或者地质统计学的著作很多,但是几乎没有对应的商业型拟合软件。这一点主要是因为变异函数拟合存在很大的不确定性:由于原始的变异函数模型及参数未知,实际过程中我们只能通过有限的观测数据去推断,因此结果不可避免的会存在不确定性。变异函数拟合的方法有很多种,具体的参见文献《Comparison of inference methods for estimating semivariogram model parameters and their uncertainty: The case of small data sets》以及其中涉及到的一些文献。在回复一些虫友的咨询中我也提到,现在国外比较认可的拟合方法主要有加权最小二乘法、极大似然估计等,对应的程序有VARFIT、MLREML,其中VARFIT已经有虫友上传到小木虫了。这两个程序都是同一个作者编写的,关于程序的介绍都有相应的文献,对这方面感兴趣的朋友可以多阅读一下这个作者的其他文章,大多是关于变异函数拟合、不确定性分析、克里格插值方面的,相信对大家会很有帮助。另外,希望大家在进行地质统计分析之前多看看专业相关的著作,中文的推荐《地统计学概论》、《空间变异理论及应用》,英文的有《GSLIB》、《Geostatistics for environmental scientists》、《Introduction to Geostatistics: Applications in Hydrogeology》、《Geostatistics: modeling spatial uncertainty》,英文的可以重点看第一和最后一本。其中GSLIB是一套完整的体系,包括几乎地质统计的各个模块,里面的很多小程序都非常实用,学习起来也比较简单。最后,跟大家提一下GS+和ArcGIS。GS+由于傻瓜式的操作,使得用户不需要太多地了解地质统计学的相关知识就可以进行空间变异分析,尤其是融合了变异函数拟合功能,使得它广受一些学者的喜欢。但是,GS+在国外的认可度是比较低的,因为GS+在拟合变异函数时采用的是最小二乘法,其拟合出来的结果精度往往会很不理想,甚至会得出明显错误的结论,而如果没有地质统计学知识背景的话,很有可能就会将错误的结果运用到自己的论文当中(有些已经刊出的文章就发现过这个问题)。ArcGIS相比于GS+,结果的可信度就更差了(至于为什么会出现这个问题,需要查看ArcGIS中地统计学模块的相关代码,奈何自己编程水平太菜,所以就留给其他人吧),但ArcGIS胜在出图美观,因此,可以外部(比如用VARFIT、MLREML)拟合好变异函数后,将相应的模型及参数代入到ArcGIS中。
回复此楼
幸福总是隐藏于平凡之中!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

孔子无敌老子

金虫 (小有名气)

引用回帖:
2楼: Originally posted by qmsdlily at 2016-03-08 06:58:32
感谢楼主!如果算出来变程是四万多米,就是明显错的吧。VARFIT软件是DOS版的,界面不够友好,而且大家看专业的英文文献多少是费劲的,它对数据格式也有较多规定。请楼主辛苦一下给大家讲讲。或者把您平常用的软件分 ...

变程的大小跟你研究区大小有关,不能说4万多米就一定错。另外英文文献看起来真的一点都不费劲,别自己吓唬自己。我用的软件主要就是GSLIB里的源程序,以及我跟你提到的那个作者编写的一些小程序。
幸福总是隐藏于平凡之中!
3楼2016-03-08 11:45:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

qmsdlily

金虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
感谢楼主!如果算出来变程是四万多米,就是明显错的吧。VARFIT软件是DOS版的,界面不够友好,而且大家看专业的英文文献多少是费劲的,它对数据格式也有较多规定。请楼主辛苦一下给大家讲讲。或者把您平常用的软件分享一下。相信大家都会感激您的!

发自小木虫Android客户端
感谢所有帮助我的人,我会努力奋斗的!
2楼2016-03-08 06:58:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qmsdlily

金虫 (著名写手)

感谢楼主!!
感谢所有帮助我的人,我会努力奋斗的!
4楼2016-03-08 17:54:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
感谢分享!
5楼2016-03-10 10:27:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sophiawater

新虫 (初入文坛)



小木虫: 金币+0.5, 给个红包,谢谢回帖
美国南达科达矿业学院 招收地质统计学方面的公派及联合培养博士生[联系:liangpingli@utexas.edu]
6楼2016-04-10 11:08:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

没有如果2015

铜虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
想请教一下,GS+做半方差分析 Active Lag  Distance   Lag Distance  interal   应该怎么确定?
7楼2016-06-23 16:57:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bsstzw

新虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
7楼: Originally posted by 没有如果2015 at 2016-06-23 16:57:40
想请教一下,GS+做半方差分析 Active Lag  Distance   Lag Distance  interal   应该怎么确定?

请问你知道怎么确定了吗,我也一直很困惑这个问题
8楼2016-11-23 10:58:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hearterode

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
8楼: Originally posted by bsstzw at 2016-11-23 10:58:53
请问你知道怎么确定了吗,我也一直很困惑这个问题...

怎么确定定了吗?求解啊
9楼2017-03-11 23:08:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

丹尼mm

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
内容已删除
10楼2017-12-25 13:00:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 孔子无敌老子 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 环境学硕288求调剂 +8 皮皮皮123456 2026-03-22 8/400 2026-03-23 23:47 by 热情沙漠
[考研] 306求0703调剂一志愿华中师范 +7 纸鱼ly 2026-03-21 8/400 2026-03-23 23:31 by chixmc
[考研] 材料292调剂 +4 橘颂思美人 2026-03-23 4/200 2026-03-23 23:16 by peike
[考研] 335分 | 材料与化工专硕 | GPA 4.07 | 有科研经历 +4 cccchenso 2026-03-23 4/200 2026-03-23 23:00 by 徐ckkk
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 070300,一志愿北航320求调剂 +3 Jerry0216 2026-03-22 5/250 2026-03-23 09:16 by 。。堂堂
[考研] 260求调剂 +3 朱芷琳 2026-03-20 4/200 2026-03-22 15:12 by 朱芷琳
[考研] 286求调剂 +10 Faune 2026-03-21 10/500 2026-03-21 23:34 by 314126402
[考研] 一志愿南大,0703化学,分数336,求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 302求调剂 +12 呼呼呼。。。。 2026-03-17 12/600 2026-03-21 17:29 by ColorlessPI
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 二本跨考郑大材料306英一数二 +3 z1z2z3879 2026-03-17 3/150 2026-03-21 02:29 by JourneyLucky
[考研] 求调剂 +3 Ma_xt 2026-03-17 3/150 2026-03-21 02:05 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 288求调剂 +16 于海海海海 2026-03-19 16/800 2026-03-20 22:28 by JourneyLucky
[考研] 材料与化工 322求调剂 +4 然11 2026-03-19 4/200 2026-03-20 22:12 by luoyongfeng
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
信息提示
请填处理意见