我有两个shapefile - 让我们称之为shp1.shp和shp2.shp-我想在两者之间创建一个重叠的新shapefile。
基本上，我正在尝试确定shp1中的哪些地理区域属于shp2。大多数shp1多边形将包含shp2多边形，并且一些shp2多边形将落在多个shp1多边形内。
如果我从
开始 library(sp)
large_list <- over(shp1,shp2, returnList = TRUE)
这让我获得了大量共享地理位置列表。但是如何获取该列表并使用它来进行空间连接，并创建一个新的shapefile？
我是R的新手（特别是对于GIS），任何帮助都会受到赞赏。
答案 0 :(得分：0)
您可以在函数intersect中找到答案，来自包raster
library(raster)
newshape <- intersect(shape1, shape2)