问题标签 [geor]

我有相关经度和纬度的数据。如何根据点之间的大圆距离获得该数据的变异函数？

这个简单的示例包含赤道上的所有数据：

第一个点和最后一个点实际上仅相隔 1 度，但我天真的尝试导致variog认为它们相隔 359 度。

摘要：我有一个包含 NA 值的栅格数据集，并且想要计算它的变异函数，忽略 NA。我怎样才能做到这一点？

我有一张使用该readGDAL函数加载到 R 中的图像，存储为im. 为了使这个可重现，dput图像上的结果可在https://gist.github.com/2780792获得。我正在尝试显示此数据的变异函数并且正在苦苦挣扎。我将介绍到目前为止我尝试过的内容：

我尝试了该gstat软件包，但似乎无法获得有效的函数调用。我已经收集到基本上我需要的是数据值本身（im@data$band1）和坐标（coordinates(im)）。我尝试了各种命令，例如：

和

我在这里做错了什么？

由于这似乎不起作用，我尝试了这个geoR包，我称之为使用：

该错误看起来与其中包含 NA 的数据有关，因此我尝试使用 删除它们na.omit，但这会将所有 NA 留在其中。这有点道理，因为光栅文件必须在每个网格正方形中有一些东西。有没有办法以某种方式删除 NA，或者至少让variog命令忽略它们？

任何帮助将非常感激。

我正在SyntaxError: keyword can't be an expression编写我正在编写的脚本。

我在 Python 中使用rpy2（和 R 包geoR）来处理使用 R 的数据库函数从 PostgreSQL 数据库中检索的数据表。数据是具有坐标和 2 列数值数据的空间数据，将在地统计模型中使用。

数据库查询调用后，数据框对象 x 如下所示：

我为 geoR 函数创建一个对象，如下所示：

从 rpy2.robjects.packages 导入 importr geo = importr('geoR')

调用 geoR 函数为

有效，但没有参数 data.col，它将位置属性分配为数据属性。（坐标属性后的第一列是默认值。）

试：

产生：

我似乎无法绕过它。我在这里查看了一些帖子并在网上四处查看，但我无法弄清楚这一点。

是否有可以生成具有截断分布的随机字段的 R 包或函数？

我正在尝试模拟对数正态空间随机场，但我需要一定范围内的模拟值。所以我需要一些易于使用的函数来生成截断的高斯场。具体来说，我需要一个GaussRF来自 RandomFields 包或 geoR 包grf的函数来生成一个具有截断边缘分布的随机场和一个具有指定范围直接的相关结构。

如果没有可用的read-to-use函数或包，我是否可以很容易地编写自己的函数或包？

嗨，我一直在尝试使用简单的数据集创建简单的轮廓。

使用的数据集如下：

我只是将上述数据合并到一个名为elevation.csv 的文件中

我使用 loess 和 expand.grid 函数进行插值。黄土模型中的度数和跨度如何选择？

我用来绘制的代码如下：

我不确定我得到的结果是否准确。任何人都可以用任何其他技术验证这个结果并分享这个观点吗？实际上，我必须处理一个大型数据集，并且想先从简单的事情开始。

我用上面的代码得到的输出如下：

这是处理大量数据集的好方法吗？

谢谢。

我想对 geoR 包中的 ca20-Dataset 进行交叉验证。例如，使用 meuse 数据集，这工作正常，但对于这个数据集，我遇到了 SpatialPointsDataFrame 尺寸的奇怪问题。也许您可以自己尝试一下并解释为什么 autoKrige.cv 函数不起作用（我尝试了几个 nfold-values 但这只会更改错误消息的位置值...）：

我希望有人能重现这个问题，我的 R 版本是 2.15，所有包都是最新的（至少不超过一个月左右......）。

谢谢你的帮助！！

我有一个带有离散点的 x、y 和 z 数据的情况。我使用 akima 包插入 x、y 和 z 点，如下所示。

当我查看变量 jd 的摘要时显示如下：

似乎我有一个由 x 和 y 点组成的网格，并且对于每个 x,y 数据集，我都有对应的 z。我想知道是否可以从此列表中创建一个数据框。

我尝试使用以下代码来创建数据框。

但是，当我尝试使用以下命令查看数据时，出现错误。

我得到的错误是

数据可以在https://www.dropbox.com/s/dbzulncdfz7pqwr/contour.csv上找到

编辑：关于如何解释数据的想法。

我认为数据格式如下所示：

所以，我相信我需要创建一个数据框，这样我应该有 1600 个点与 X 和 Y 数据的组合。因为我有 40 X 和 40 Y，所以应该可以创建 1600 个网格点，然后所有的 z 点将是这些点的对应 z。我可能错了。

谢谢。贾巴巴

我使用“geoR”包拟合了我的半变异函数，如下所示：

有没有办法知道指数模型的拟合优度（轻松）？

提前致谢 ...

我想知道是否有人可以向我解释 geoR 包如何计算协方差函数？我的意思是你会怎么用手做？

你得到：

然而，人们可能也想在 Matlab 中做这件事。

你们中的一些人可能已经看过超越“苏打水、汽水或可乐”。我面临着类似的问题，并想创建一个这样的情节。就我而言，我有大量的地理编码观测（超过 100 万）和一个二进制属性x。我想在p(x=1) 的色标范围为 0 到 1 的地图上显示x的分布。

我对其他方法持开放态度，但这里描述了 Katz 的超越“苏打水、汽水或可乐”的方法，并使用这些包：字段、地图、mapproj、plyr、RANN、RColorBrewer、比例和邮政编码。他的方法依赖于使用高斯核的 k 最近邻核平滑。他首先定义了地图上每个位置t到所有观测值的距离，然后对p(x=1|t)使用距离加权估计（x 为 1的概率取决于位置）。公式在这里。

当我正确理解这一点时，在 R 中创建这样的地图涉及以下步骤：

1. 构建覆盖 shapefile 整个区域的网格（我们将网格中的点称为t）。我尝试使用这种方法polygrid，但到目前为止失败了。代码如下。
2. 对于每个t，计算到所有观测值的距离（或者只找到 k 个最近点并计算该子集的距离）
3. 根据此处定义的公式计算p(x=1|t)
4. 使用范围从 0 到 1 的适当色标绘制所有t

这是一些示例数据和两个具体问题。首先，如何解决我的第 1 步问题？正如下面的第二张地图所示，我目前的方法失败了。这是一个明确的 R 实现问题，一旦解决，我应该能够完成其他步骤。其次，更广泛地说，这是正确的方法还是您会建议一种不同的方法来创建具有属性值分布的热图？

加载库并打开 shapefile 和包

绘制数据

1）构建覆盖shapefile整个区域的网格