我希望通过属性对 SpatialPolygonsDataFrame 进行子集化，但我想允许它返回一个空的 SpatialPolygonsDataFrame。

如果我们要将 SpatialPolygonsDataFrame 类型的对象像 data.frames 一样处理，如这里所讨论的，我们应该能够获取并使用空对象。

我很感兴趣，因为我想将它合并到一个函数中，该函数可能会尝试通过一个本质上不会选择任何特征的属性来进行子集化。

有没有办法使这项工作？如果没有，有没有办法初始化一个空的 SpatialPolygonsDataFrame 对象？我想在这样的对象上执行的未来操作涉及在现有地图上过度绘制，所以我希望无论如何都能生成图像，即使是空白的。

如何使用spplotor生成以下图ggplot。我尝试了两种方法，但不幸的是，它们中的任何一种都没有占上风。我面临的问题与地图的适当添加点有关。.shp,.prj,.dbf,.shx重现问题所需的所有必要数据 ( ) 均在此处提供

我的代码如下：

我通过执行上面的代码行得到的结果是：

来自的点pts1必须绘制在第一张地图上，而来自的点pts2则应分别描绘在第二张地图上。在我的情况下，您可以看到两个点都存在于两个地图中。知道如何解决问题吗？

最近，我找到了一种正确描绘点的方法，但是这样标题就消失了：

在这部分代码之后，我得到了想要的外观，这不是很漂亮但可以接受。 任何帮助都将受到高度评价！

我在 R 中绘制真实的地理范围像素时遇到问题。这些文件带有每日单个坐标和像素大小（区域）的列表。还有一个与此分开的 Z 元素。数据结构看起来是这样的：

X 和 Y 以经度和纬度为单位，而面积以平方公里为单位。我使用以下方法轻松创建点特征：

我使用面积值轻松绘制这些图以确定绘图的“cex”。Z 列是强度，我使用这些值来确定颜色。如何使用 R 中每个点的区域创建空间多边形特征？我将使用这些点来创建网格栅格。

我想使用 [LeSage & Pace 提出的模型规范] ( http://www4.fe.uc.pt/spatial/doc/lecture7.pdf ) (p.8, p.8,等式[10]）：

其中o代表始发地，d代表目的地。

我认为可以将此模型称为空间自回归交互模型。

该模型的关键在于，我们有起点和终点的空间滞后响应，以及SAR 方程中起点和终点的空间滞后预测因子。

---

• 您对如何在 R 中实现这样的模型有任何想法吗？
• 有没有办法以spdep::lagsarlm特定的方式或至少以类似的方式指定 a ，或者我应该使用不同的包/功能？据我所知，那里spdep::lagsarlm只允许一个 listw参数。

我希望能够从 R 中的等高线创建一个高程图。我对使用形状文件非常陌生

目前我已经从这里下载了数据， 这些数据提供了适用于整个英国的 .shp 文件。

它还提供了等高线，总结了英国的拓扑结构。

对于高程图，我想要一个data.frame或data.table多个均匀间隔的点（彼此相距 100m）来产生一个数据输出，给出一个 x、y 和 z 值。其中 x 和 y 代表纬度和经度（或东向和北向），z 代表高度（以米为单位）。

我认为可能有一些工具会自动为您执行插值，但不确定它如何处理地理空间数据。

这是我的基本开始...

这可能是一个简单的问题，但是，我想查看地图上的每个输出区域，从这里下载了寄宿生。选择 2011 OA 并获得完整的分辨率，领域的范围。

我已经对其进行了压缩...并使用以下命令读取了 shp 文件

如何绘制或查看每个输出区域的坐标？

这是我在 StackOverflow 上发布的第一个问题，因此请更正我犯的任何协议错误。但是我已经从对网站上许多其他问题的回答中受益匪浅，所以在我继续之前，我要感谢为该网站做出贡献的社区。

概述

我在 Windows 7 机器上使用 R、RStudio、包和包（devtools问题结尾处完整）roxygen2session_info

我正在尝试编写自己的包，它SpatialLines从包中定义了 S4 类的子类sp，并且还使用了包rgeos。document当我使用 devtools和 load_all我的包时，我的问题与我不太了解且无法完全解决的警告消息有关。

警告信息

第一条消息是

当我解决这个问题时，通过在 NAMESPACE 文件中添加一个条目（不完全理解为什么它是必要的），然后对DESCRIPTION 文件进行进一步的更改，我收到了第二条消息，到目前为止我已经被打败了。

我所尝试的，以及我留下的问题

我已经用谷歌搜索并搜索了几个寻找提示的网站，但我找到的最接近的是 2011 年 12 月的一个简短线程，询问 SpatialLinesNULL 来自哪里 --- 答案是rgeos。 http://lists.r-forge.r-project.org/pipermail/rspatial-devel/2011-December/000033.html

我已经尝试自己尝试了解出了什么问题，并在调试器中逐步执行代码，直到我迷路为止，这比我能够很好地理解应该发生的事情要深入得多。在此过程中产生了一个小（非）问题的工作示例，如下所示。

我有 3 个问题。它们在此处进行了总结，但实际上仅在下面我的（非）工作玩具示例的上下文中才有意义：

• 问题 1（在 test1 之后）为什么 R 仍然需要（虚拟？）超类Spatial，为什么只有在我的代码中定义了第二代 MyClass2时？
  澄清补充说@importClassesFrom，我曾认为，一旦我编辑了我实际使用的类，R 打包命名空间机制就会满足查找和访问前项/超类的任何需要。
• 问题 2（在 test3 之后）我假设无论我的第一条警告消息的原因是什么，也是这个类似消息的原因。
  我在正确的轨道上吗？
• 问题 3（在 test4 之后）来自 test3 的警告和随后的错误是我可以自己解决的吗（如果可以的话）？还是我需要向rgeos软件包的维护者寻求帮助（例如要求他们导出SpatialLinesNULL）？

代码示例

我的玩具工作示例的 r 代码（仅去除 roxygen 注释以节省空间）是：

MyClass2 第二代子类似乎很重要（没有它，警告就会消失）。我不懂为什么。

devtools我使用and 对此代码运行了一系列测试roxygen2。测试因 NAMESPACE 文件中的内容以及DESCRIPTION 文件中的 Imports 指令中的内容而异。执行测试包括在 RStudio 中运行devtools::load_all()（或沿途devtools::document()调用）。load_all()发生的警告消息已在上面发布。我在每次测试之间在干净的全局环境中重新启动我的 R 会话（警告仅在我第一次运行时出现load_all或document）。

在跟进下面给出的答案时，我意识到我也应该尝试从 shell 命令行使用 R 构建和安装我的玩具包。

测试1

test1 有DESCRIPTION文件条目（完整的DESCRIPTION文件见这个问题的结尾）：

完整的 NAMESPACE 文件（从上面省略的 roxygen 注释生成）是：

尝试在重新启动的 R 会话中运行devtools::load_all()（或devtools::document()）会生成第一个警告（如上所示）。

问题 1为什么 R 需要 ~~(virtual?)~~ 超类Spatial，为什么只有MyClass2在我的代码中定义了第二代时？ 澄清补充说@importClassesFrom，我曾认为，一旦我编辑了我实际使用的类，R 打包命名空间机制就会满足查找和访问前项/超类的任何需要。

测试2

test2通过更新 NAMESPACE 文件来解决问题，如警告消息所示。没有进行其他更改。

当完整的 NAMESPACE 文件为：

测试3

但是，当我对DESCRIPTION 文件中的Imports 指令（即添加rgeos. 这是从 test2 创建 test3 场景的唯一更改。

现在我收到一个警告，看起来像是来自 test1 的问题，即：

据我所知，SpatialLinesNULL是一个新的超类Spatial，定义在rgeos它可以提供一些通用函数和方法。

问题 2我假设无论我的第一条警告消息的原因是什么，这也是此类似消息的原因。
我在正确的轨道上吗？

测试4

在 test4 中，我尝试应用与 test2 中相同的解决方案，即将缺失的条目添加到 NAMESPACE 文件中。这是 test4 的完整 NAMESPACE 文件。

在哪个阶段我收到最后的错误消息：

在这个阶段我被卡住了。

问题 3来自 test3 的警告和随后的错误是我可以自己解决的吗（如果可以的话）？还是我需要向rgeos软件包的维护者寻求帮助（例如要求他们导出SpatialLinesNULL）？

session_info（来自 test4）

test4 的完整文件

R/myclasses.r

描述文件

命名空间文件

结论

非常感谢您提供的任何帮助、建议或指导！即使这只是关于我如何提出更好的问题:-)

我想绘制一张我有两个变量的地图。一个是多边形，另一个是点。

多边形 - 省边界（name1、name2、name3、name4），

点 - 类（a、b、c、d）

因为我必须将它们绘制在一起，所以我需要两个图例。下面是我得到多边形特征图例的代码。请建议如何为点特征添加图例。

我得到的地图：

如您所见，只有PROVINCE（多边形）的图例。我怎样才能将它添加到点特征呢？

我想spdep::lagsarlm在 R 中估计一个模型（空间自回归）。我的观察结果（n=447）是多边形，每个代表柏林的一个行政区域。

然而，问题在于这些地区的居民数量差异很大（在 500 到 32000 之间）。因此，我想用它的居民数量来衡量每个观测值。这很容易，因为lm它接受可选参数weights=......

我怎样才能做类似的事情spdep::lagsarlm？有解决方法吗？

如果我有一组这样创建的 SpatialPolygons：

这样它看起来像这样：

以及由以下内容创建的点向量：

如果有的话，我如何检测哪些点位于哪个多边形内？

谢谢