问题标签 [sp]

从上面的以下代码和输出中，我可以发现log(zinc)~1并 log(zinc)~x+y返回 v 和 v1 的不同伽马值。这两种操作的基本区别是什么？

help(variogram)告诉我“定义响应向量和（可能的）回归量的公式，在没有回归量的情况下，使用例如 z~1;”。但是，我没有清楚地理解这句话！谁能详细告诉我应该何时使用z~1或何时使用z~LON+LAT？

我有一些空间数据，我使用下面的代码创建热图并将创建的关卡提取为多边形。我现在的问题是如何限制创建的关卡数量？我的目标是拥有例如 5 个不同的密度级别？

在第二步中，我想提取属于一个特定级别的多边形。由于我没有来自可重现示例中使用的正态分布的数据，因此可能存在具有相同密度级别的不同多边形。

这是我到目前为止的代码：

编辑：感谢@pHroc 的回答，我能够控制级别的数量，并且我还发现了如何提取具有相同级别的多边形。但是现在我遇到了一些创建区域非常小的问题。有没有办法控制一个区域的最小大小或每个区域至少应该包含的点数？

我正在使用 R使用 poly_center （pracma 包）计算在此文件中读取的图形定义的多边形的质心。

奇怪的是，我得到的结果超出了多边形边界，中心点低于 x（lon）和 y（lat）轴/坐标的最小值（参见下面的代码）。谁能提示我我做错了什么或这里发生了什么？

我有一个列表（从 Stata dta 文件导入），其中包含为地图定义空间多边形的条目，以提供一个想法：

我想把它转换成一个数据结构，我可以用它在 R 中创建一个地图（没有更大的难度，也不需要广泛的专业知识，以前从未做过）。我推断 sp 包的 SpacialPolygons 类型是最简单的选择。此外，从定义来看，似乎 SpacialPolygons 方法（在此包中定义，请参阅包的文档，第 79 页）似乎是从列表转换为此数据类型的正确方法。

不幸的是，该方法并不那么简单，我需要一些帮助。我的（天真）尝试产生了一个我不理解的错误，并且在谷歌搜索中没有出现任何有趣的结果：

上面的列表可以转换为 SpacialPolygon 吗？如果是这样，怎么做？如果不是，我应该选择什么格式？谢谢。

我有一个多边形形状文件和一系列数据框。我想要做的是创建一个循环，在形状文件和所有这些数据帧之间进行“空间连接”，但我不断收到错误消息。对此有什么想法吗？

这是我的循环代码

这是一个分三个步骤的可重现示例。

(1) 从随机多边形创建规则网格多边形

(2) 创建数据框

(3) 尝试循环

如果我一个一个地进行空间连接，我会得到正确的结果，但这并不能解决我的问题，因为我有太多的数据集

我有一个从具有要素几何属性的 MySQL 数据库中提取的数据框。例如，这是 geom 属性的一条记录值：

表中的所有要素都是多边形。

我一直在搜索，但找不到从这些坐标字符串创建空间数据框的干净方法。

我需要快速确定空间多边形和空间线是否相交。我目前正在将 polgon 转换为空间线并使用gIntersection(). 任何人都可以提出一种可能更快的方法吗？也许使用栅格而不是空间线或其他东西。我需要这样做数千次。

我正在尝试gDifference在两个多边形之间进行处理，但出现错误：

所以我这样做了：

但现在我遇到了另一个错误：

有没有办法来解决这个问题？在 ArcGIS 中使用 Erase 功能做同样的事情效果很好。

我想合并一个SpatialPolygonsDataFrame：

使用普通数据框：

目标是使用结果制作地图leaflet

我试过 sp::merge

但我收到一个错误：

使用 R中的sp包。我想知道何时使用 SpatialPoints 以及何时使用 SpatialPointsDataframe。在我看来没有太大区别？！

唯一的区别是在 SpatialPointsDataframe 中我可以存储更多属性吗？！如果是这样，我可以从现有数据帧（如果坐标存在于该数据帧中）创建一个 SpatialPointDataframe，而无需绕道创建空间点吗？