问题标签 [sp]

如何按时间对spacetime::SDFDF（具有完整时空网格的时空数据）进行子集化？

到目前为止，我尝试过：

但这给了我error: ncol(i) == 2 is not TRUE.

有任何想法吗？

我有一个英国道路网络的形状文件。我还有一个连接到车辆的 GPS 跟踪设备的经纬度。我可以为每次旅行将 GPS 经纬度转换为 SpatialLine。然而，GPS 跟踪设备中存在一些错误的数据点。每个连续的经纬对彼此相距约 2 分钟。在地图上绘制经纬度时，很容易在视觉上注意到“错误”点。有没有办法通过评估创建的空间线与道路网络形状文件的交集来以编程方式执行此操作？即，鉴于我们知道这些点应该大致遵循道路的形状，是否有任何 R 包中的函数能够评估这一点？

我知道gInterserct，但我不只是在寻找真假......结果......我试图找出一系列点中的哪些点是“错误的”......

伙计们，

我需要一些专家意见。我一天中的大部分时间都在与这个问题作斗争，但我别无选择。

我正在尝试创建一张地图，描绘自上次火灾以来整个景观的时间

很明显有些地方是不会发生火灾的，这里用0表示，但是有点随意

我的代码如下所示：

数据可以在这里找到

如果你不想运行代码，这就是它在谷歌地球上的样子

这是主要问题，其中：

1）它几乎看起来像多边形边界信息被保留，因为棕褐色线是确切的边界，

2）这里并不总是报告“最大值”，尤其是在它出现的边缘附近？

3）我的目标是在物种分布建模中将此地图与其他栅格结合使用（因此使用栅格堆栈、模型和预测功能来制作地图）。当我尝试使用该图层制作地图时，由于 NA 的缘故，我最终得到了一堆洞......应该没有 NA 的地方。

如果我在没有这个时间的情况下运行模型，那么一切都很好。当我包含它时，生成的地图中的 NA 通常是您在第一张图片中看到棕褐色“多边形轮廓”的位置。

我怎样才能让这个光栅层有清晰的边缘，而不是像现在这样一团糟？？？

这就是我的最终地图物种地图的样子，而不是丢失的数据反映了火灾地图中的问题

我想找出我的对象的精确地图比例尺（以 m 或 km 为单位）SpatialPolygons并将其绘制在我的对象上，spplot以便读者了解距离。这意味着我想找出???在情节中输入的数字。

也许有什么sp::coordinates或Imap::gdist？任何先前的重新投影都SpatialPolygons需要找出这一点？

我有一些想要合并的相邻宗地集群，但我希望保留每个集群内的宗地边界。如果我使用类似的东西，Reduce(geom_list, gUnion)那么相邻的边界就会被溶解。

这是我想要的一个简单且可重复的示例，或多或少在纯 R 中实现：https ://gist.github.com/holstius/37c0ff01109e28ff817a 。GEOS如果还没有本机实现它，我会感到惊讶，但我无法确定rgeos可以做我想做的绑定。

是否有不会溶解相邻多边形的变rgeos::gUnion体？

（编辑：粗略的代码示例%>%已删除以支持https://gist.github.com/holstius/37c0ff01109e28ff817a）

是否可以像这样将 spplot（spplot 多边形）图例放置在地图的左下角？

我能得到的最接近的是这个（我没有发布我的数据，我只是使用示例数据，所以在这种情况下，尝试将图例放在地图的左上角）：

但是图例在页面中间，在地图之外。不幸的是，colorkeyargument 不支持“bottomleft”、x、y 或角参数（请参阅 参考资料?levelplot）。我也尝试使用key.space参数，但它似乎只在绘图时有效，SpatialPoints*但它似乎被忽略了 SpatialPolygons*（或上面示例中的 SpatialPixelsDataFrame）。

如何在 spplot 颜色键中有对数色标？即在这张地图中：

比例尺上的数字不是 0.0、0.2、0.4、0.6、0.8、1.0，而是像 0.0、0.01 0.1 0.2 0.5 1.0 或类似的东西？

我尝试了这个do.log论点，但它只适用于 SpatialPointsDataFrame （根据文档）。我还尝试了 scales 参数：

但它似乎只适用于 x 和 y 坐标，而不是显示的变量和色标。

我希望能够创建一个覆盖另一个 SpatialPolygon 的 SpatialPolygons 对象（它是一个六边形网格）。

我希望所有六边形的直径为 1 公里（理想情况下我可以改变它）并且所有六边形一起覆盖整个物体。下面的方法似乎只涵盖了其中的一小部分......

sp以下是我使用该软件包的尝试：

非常感谢任何帮助...

注意：在 Edzer Pebesma 的建议下，这个问题已被交叉发布到 R-sig-geo，在这里，它得到了一些很好的回应。

我遇到了以下意外结果checkPolygonsHoles：

这里一个恼人的副作用是原始对象wrld_simpl也被改变了。这个结果在我看来就像一个错误，还是我错过了什么？

PS：之前shape_with_hole编辑的对象checkPolygonsHoles，继续表现奇怪：

我用 utm32-shapefile 制作了一个等值线图，现在，我想添加第二层，其中包含来自另一个数据帧的数据。第二个数据帧只有 Lon/Lat 信息。我想将 Lon/Lat 转换为 utm32 坐标，因此我可以将它们与 ggplot2 一起用作另一个收敛层。

数据框如下所示：

首先，我使用gsub将“，”重新编码为“.”。然后我将变量转换为数字。

它没有用。我只是得到一份我的纬度/经度变量的副本。

我没有找到有关如何执行此操作的良好文档。也许其他人可以帮助我？