问题标签 [geostatistics]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
0 回答
241 浏览

r - geoRglm R包中二项式数据的回归克里金法

我正在使用binom.krige()R 包geoRglm的函数来确定具有多个连续和离散协变量的二元 (0, 1) 响应变量的空间预测。

使用glm()二项式 logit 链接函数,我发现响应变量显示出对几个协变量的显着依赖性。

我将趋势包含在binom.krige()使用krige.glm.control()中,我将两个趋势模型指定为

让我感到困惑的问题是,何时进入trend.d并最终进入,它实际上是否适合带有二项式 logit的glm或只是线性模型(因为上述方程似乎是线性模型)?trend.lkrige.glm.control()binom.krige()

0 投票
0 回答
39 浏览

r - 从 gstat 中的变异函数获取滞后向量

我想从 R 中的一组数据计算变异函数。我正在使用 gstat 包中的函数“变异函数”。

现在,我想从变异函数中获取滞后向量。问题是 myvariogram$dist 返回所有点之间距离的平均值。

我怎样才能得到滞后向量呢?

我的数据是二维的:坐标 x 和 y 与 z 值

0 投票
0 回答
61 浏览

statistics - 一些Neighbourless空间点的调整

请找到随附的 R 代码。

**主要问题是这些地区有一些无邻域的空间点。因此,当我尝试通过生成的空间权重矩阵计算 moran's I 索引时,它显示了无邻域的错误。

我没有得到,如果一些空间点没有邻居,如何调整空间权重矩阵。

请帮帮我。提前致谢

0 投票
1 回答
359 浏览

r - 时间变异函数

我想建立一个时间变异函数,而不指定空间坐标,因为我只有一个我正在研究的量规。这是我的数据示例:

日期采用 POSIXct 形式,我每 360 秒或 6 分钟测量一次 RR6(感兴趣的变量)。我怎样才能做到这一点 ?谢谢 !

0 投票
1 回答
36 浏览

r - 在 R 的 gstat 包中定义新的相关模型?

vgm是否可以在 R 中的 package 函数中定义和使用新的相关模型gstat,以模拟随机场?

我需要对对数转换的指数相关空间数据使用相关函数(如此处所建议)。可用的选项vgm有“Exp”、“Sph”、“Gau”、“Mat”

0 投票
1 回答
435 浏览

r - 带有“gp”平滑器的 GAM:如何检索变异函数参数?

我正在使用以下geoadditive模型

在这里,我使用范围 = 10 和幂 = 1 ( m=c(2,10,1)) 的指数协方差函数。如何从结果中检索变异函数参数(块金、窗台)?我在模型输出中找不到任何东西。

0 投票
1 回答
47 浏览

geolocation - 一旦用户停止移动一段合理的时间,确定 GPS 跟踪用户位置的最佳策略是什么?

一旦认为用户没有移动,从一系列 GPS 定位中确定用户位置的最佳策略是什么?

在跟踪用户时,如果他们应该停止移动,随后将在同一位置进行一系列修复。

如果可能的话,我不仅想使用最后的修复,还想考虑以前的修复,以便为它们计算更准确的位置。

我认为需要考虑的因素:

  • 确定用户静止不动的最佳方法(从 GPS 定位的经验速度来看不够可靠)
  • 每个修复都有一个准确性,如何将其考虑在内?
  • 是否有可以使用的完善的算法/库?

任何建议都非常感谢

0 投票
1 回答
278 浏览

r - 如何在R中计算变异函数时更改滞后距离

我正在尝试计算不同滞后距离的 Experimental-Variogram 值,所以我正在使用variogramm命令

但我不知道如何指定所需的滞后距离。例如,我想在h=(15, 30, 45, 60)

0 投票
0 回答
48 浏览

python - 如何在网格数据上拟合方向变异函数,并为给定的方位角返回估计范围?

我有一个地形特征的数字高程模型。我的总体目标是通过回归给定方位角的方向变异函数的范围来制作预测模型。是否有任何好的 python 库可以适应网格数据?

0 投票
0 回答
29 浏览

kriging - 使用 scikit-geostat 时,金块值始终为零

在此处输入图像描述

从图中可以清楚地看出,金块值不为零,但它给出的值为 0

输出: