在R中用SpatialPolygons覆盖SpatialLines

Question

我有一个代表道路的SpatialLines对象和一个SpatialPolygons对象 包含城市。

我想知道如何通过SpatialLines覆盖SpatialPolygons对象 R.中的对象。

我想知道经过这两个城市和哪些城市的线路？

有可能吗？

Answer 1

您需要包装GEOS几何操作库的rgeos包。

然后gIntersection或gIntersects可能会做你想要的。

Answer 2

通常，空间对象的叠加层由over包中的sp函数处理。在over（?over）文档的“方法”部分中，列出了over的方法，其中不包含Lines和{{1}的组合}。幸运的是，Polygons的文档说通过安装over包，可以使用这些方法。