我有一个点数据集,我想使用Krige,最好使用开源软件包。如果可能的话,我也想在过程中选择半变异函数模型以改善估计。

#1 楼

根据您要应用的克里金类型,有不同的软件包可供选择:

GRASS-v.krige
SAGA-模块:普通克里格法

HPGL-normal-kriging()(PDF页面12)

简单Kriging

简单Kriging使用整个数据集的平均值,而普通Kriging使用本地平均值。因此,简单克里金法可能不太准确,但是通常会产生“更平滑”的结果。它的实现方式为:



HPGL-simple_kriging()(PDF第12页)
gstat krige

通用克里金法

通用克里金法可考虑数据漂移。包括在以下实现中:


SAGA-模块:通用克里金法.krige还支持块Kriging。

HPGL实现了许多鲜为人知的Kriging方法(有关详细信息,请查看手册): (IK)
局部变量均值Kriging(LVM克里格)
简单CoKriging(Markov模型1和2)
顺序指示器模拟(SIS)
心电图局部变量均值SIS(CLVM SIS) )
局部变量均值SIS(LVM SIS)
顺序高斯模拟(SGS)
截断高斯模拟(GTSIM)[在Python脚本集合中]

SAGA提供了不同普通和通用Kriging版本。

Gstat krige还支持块和点Kriging。

评论


伟大的克里格答案!

–拉吉·亚瑟(Ragi Yaser Burhum)
2011-11-22 4:37

您的动物区系链接无效...

– Alex Leith
13年13月13日在22:23

#2 楼

看起来GRASS GIS有几个选项。请查看GRASS Kriging Wiki页面:http://grass.osgeo.org/wiki/Kriging

2009年的Google Summer of Code项目制作了V.krige:http://grass.osgeo。 org / wiki / V.krige_GSoC_2009

GPL gstat软件包应单独运行或与GRASS GIS交互。 http://www.gstat.org/

Dylan Beaudette有一个很好的使用GRASS进行克里金法的例子。 http://casoilresource.lawr.ucdavis.edu/drupal/node/438(他的博客中充满了许多使用OpenSource GIS和统计工具的有趣示例!)

#3 楼

R项目拥有大量的空间统计软件包,但是R的学习曲线相当陡峭。

评论


人们总是这么说,但是我想知道:相对于什么陡峭?

–马特·帕克(Matt Parker)
2010年8月11日19:20

我已经多次在R上看到“深度学习曲线”注释-这对我来说毫无意义。当我发现R时,我与MATLAB的关系已经一年了。我发现R非常容易学习,我向MATLAB致以敬意,并立即放弃使用它。

–夏皮
2010年8月12日上午8:33

我认为这是因为人们很少尝试理解统计信息,并且因为它具有陡峭的学习曲线语法,所以很少会遇到问题

–多多巴
2010年8月12日在8:51

我认为从语法角度来看,这是一种较容易学习的语言。可以从命令行轻松学习的统计语言示例是什么。我认为人们抱怨,因为它不是Excel。

–TheSteve0
2010年8月14日下午4:39

与基于GUI的程序相比,它是步骤。如果您一生都使用Windows和基于GUI的程序,则在看到命令行外观时将运行。 Excel比较确实使它们运行。但是,如果有人可以向您展示基本技巧,R的使用非常简单。您必须准备好学习新概念,例如向量,矩阵,函数,循环,这在Excel / Windows世界中是不存在的。如果您以前使用过Linux,那将不是一个阶梯。

– Niculita Mihai
13年2月1日在12:27

#4 楼

如果您愿意将栅格读入一个numpy数组(gdal可以做到这一点),则可以使用Python或C / C ++的High Performance Geostatistics Library实现。

HPGL实现了以下算法:



简单Kriging(SK)
普通Kriging(OK)
指示器Kriging(IK)
局部均值Kriging(LVM Kriging) )
简单协同克里金法(Markov模型1和2)
顺序指示器模拟(SIS)
Corellogram局部均值SIS(CLVM SIS)
局部均值SIS(LVM SIS)
顺序高斯模拟(SGS)
截断的高斯模拟(GTSIM)[在Python脚本集合中]



我自己并没有使用它听说过有关它的好事,尤其是在速度方面。

#5 楼

查看这本免费书籍,它是有关在R中进行地统计的,并且还包含有关在SAGA和GRASS中进行地统计学的信息。
http://spatial-analyst.net/book/
http:// space-analyst.net/book/sites/default/files/Hengl_2009_GEOSTATe2c1w.pdf

#6 楼

我记得几年前使用SAGA进行此操作以获取一些洪水模型输出。
开源,非常值得一看。

#7 楼

gvSIG(另一个免费的GIS)确实允许使用Sextante进行克里金法。这与使用SAGA基本相同,但是gvSIG提供了更“典型”(即类似ESRI的)gis体验。

#8 楼

您可以尝试使用Surfpack 1.1版中的Kriging模型(我还是在DAKOTA团队中时写的),也可以尝试DAKOTA的“稳定”版本随附的最新最好的版本(Surfpack是DAKOTA的子包) ,它是从相关函数而不是半变异函数的角度进行通用Kriging的。

最近,用户Joel Guerrero将其与众多其他实现方式进行了对比,并指出
“总是与surfpack有关,我们正在将其与其他实现(包括商业实现)进行比较,到目前为止,它的表现都优于所有其他实现,以至于有时似乎在做黑魔法”。

#9 楼

GSLIB(地统计软件库)是由斯坦福大学开发并于1990年代发布的一流文件/命令驱动软件,在过去的十年中进行了一些维护。可以使用Fortran编译器在Linux / Windows上免费下载和编译源代码。有在线资源和一本书。

克里格软件是该软件的优势之一: jackknifing
SK,OK,英国,带有外部漂移的克里格格
cokriging
指示器克里格格