问题标签 [geostatistics]

我正在对我的数据进行空间分析，并根据经度和纬度对其进行定位。我用 geoR 和 gstat 计算了 R 中的样本变异函数，但我得到的是具有等级距离（经度和纬度）的样本变异函数。我想获得相同的信息，但距离以公里为单位，就像我使用 variogramST (gstat) 计算时空变异函数时发生的情况一样。

非常感谢那些将要回答的人，这对我有很大帮助！

对于克里金法，我需要计算mesh长度为的大型数组20000。下面的代码工作正常，特别是对于小mesh长度（< 100），但是，这么大的计算时间mesh非常长（大约 45 分钟）。范围的长度data在 300 到 800 之间。下面是一个工作代码：

有没有办法优化代码，从而减少计算时间？

我正在尝试使用 pykrige 以 3D 形式计算土壤体积。我正在使用普通克里金法，但是一旦我对数据进行了克里金法，我不确定如何将其写入 ascii 文件，这是我的代码：

我需要将 kriged 数据导出到 ascii 文件中，以便可以通过其他地方的切片对其进行可视化。

谢谢！

我试图读入我的形状文件然后绘制它，但似乎 RStudio 卡在了某些东西上并且不会完成运行绘图功能。现在我有：

但它不会产生情节，而且它似乎一直被困在某些东西上，因为它只说

在 R 控制台中。我是不是做错了什么导致这种情况，有没有更好的方法来读取和绘制 shapefile？

I'm trying to plot a variogram of a Spatial data frame I got but the axis in the plotted image are wrong, they don't allow all points to be plotted. The R-studio version I'm using is 1.3 in Ubuntu 20.

Here is the code:

Here is an extract of the data frame (of a total of 150 elements).

However, the plot is this. I'm learning R and I don't know what do the "cutoff" and "width" mean. How can I get the entire plot without those cuts in the distance axis?

Besides, the dput command shows this info: img1 img2 img3

我试图在一个图中绘制两个半变异函数，但不幸的是它不起作用。我想解决方案很简单，但我的拉丁文已经结束了。

这是我想放在一个图中的变异函数的代码：

其次，我还想在一个图中绘制两个半变异函数，但由于值不同，右侧有第二个 y 轴。

我只能显示每个变异函数的两个点，但不能在里面安装模型。这是我试过但不起作用的代码！

谢谢你的帮助！！！

我正在尝试绘制一张带有特定颜色“#2D3E50”的巴西地图。但是，在地图上使用这种颜色从最小区域（浅色阴影）到最大区域（深色阴影）进行着色是行不通的。请参阅下面我尝试为不同比例插入的所选颜色。

看到它自动填充绿色，而且颜色是相反的，因为较暗的颜色显示较低的区域。

我正在使用包 gstat 将指数模型拟合到经验半变异函数。我的数据集（数据）具有十进制度的纬度和经度值。使用的估计fit.variogram()范围与我通过查看模型图所期望的范围不匹配。请参阅下面的代码和数字。

这是来自的输出fit.variogram()：

从这里开始，我预计与 ~53 (26.66719+26.26643) 的窗台相对应的范围为 25.76976。然而，当查看该图时，对应于 ~53 的窗台的范围大于 60 公里。

我需要帮助。

我正在使用克里金函数来插入温度和流量数据。我有一个不规则间隔的大约 240 个数据的数据集。我正在使用 skgstat 的变异函数导入函数来计算我的变异函数。然后我使用 OrdinaryKriging 函数（也从 skgstat 导入）来计算我的克里金点，然后创建一个网格。这是我的代码：

然而，渲染并不令人满意。我的观点差距很大。我想知道如何增加克里金函数的插值半径。我找不到任何执行此操作的函数输入参数。

预先感谢您的回答

我在使用 RasterStack 循环时遇到了一点问题。我正在尝试提取过去 10 年中北卡罗来纳州每个县每小时的平均温度。我能够运行一些简单的代码并相对轻松地为每一层获取此代码，但是当我尝试循环此过程以获取光栅堆栈/砖中每个时间戳的平均值时，我运气不佳！我将不胜感激任何帮助或建议。

示例和更多信息：

我可以很容易地获得单个图层所需的内容：

但是当我尝试运行上面提到的循环时，它要么失败并且没有输出，要么给出以下错误：