在R中使用相交的shapefile裁剪shapefile

Question

我正在尝试裁剪(不是根据范围裁剪)美国州和县的shapefiles，以匹配这个DEM (数字高程模型)裁剪的栅格的边界。

​

​

我的裁剪(根据范围)状态shapefile如下所示：

​

如何裁剪状态shapefile以匹配DEM的复杂轮廓？我试着将两者相交，这很有效，但这不能区分州和分水岭的边界，使得地图非常难以阅读(特别是与县！)。

​

huc6_cty_intersect <- raster::intersect(huc6s, counties)  

我需要能够覆盖的州和县分别，所以我可以改变他们的边界宽度，颜色等，所以地图是可读的。

我试过对州和县的shapefile进行栅格化，然后用dem进行掩蔽，但这似乎需要很长时间(过夜的时间还不够)。

counties_rast <- rasterize(counties, dem)
r4 <- raster::mask(counties_rast, huc6_cty_intersect)

我试图得到下面的图像，但所有分水岭以外的州边界都被砍掉了。如何做到这一点的任何帮助都将是非常好的。谢谢。

​

注意:分水岭边界内的白色立面是故意遮挡的。

Answer 1

如果我没听错的话，你是在找raster::crop

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

pb <- crop(p, b)

如果我不能很好地理解你，那么请像我一样用简单的例子数据来阐明你的答案。引用其他人没有的对象没有什么价值。