在 R 中裁剪空间文件

Question

我目前正在尝试将一个空间文件剪辑到另一个空间文件，但是我得到的空间文件不是正确的形状。我知道使用光栅必须使用蒙版，是否有类似的命令用于 2 个形状文件？

c1 <- crop(spatial1, spatial2)

Answer 1

使用 SpatialPolgyonDataFrame 对象

library(raster)
p <- shapefile(system.file("external/lux.shp", package="raster"))
b <- as(extent(6, 6.4, 49.75, 50), 'SpatialPolygons')
crs(b) <- crs(r)

你可以的

pc <- crop(p, b)

您需要安装包rgdal 和rgeos。