将栅格裁剪到 R 中矢量 (.shp) 的精确轮廓？

【问题标题】：Clip a raster to the exact outline of a vector (.shp) in R?将栅格裁剪到 R 中矢量 (.shp) 的精确轮廓？
【发布时间】：2021-06-01 20:11:43
【问题描述】：

我想在 R 中将光栅层精确地裁剪为 shpfile 的轮廓。

  crop_rast      = crop (raster , extent(vector))                 # Crop
  mask_rast      = mask (crop_rast, vector)                       # Mask

这不起作用。

【问题讨论】：

这能回答你的问题吗？ Clipping raster using shapefile in R, but keeping the geometry of the shapefile
@Val 否，此解决方案是保持矢量的几何形状（即 WGS84）。它对我不起作用。
您需要发布一个包含一些数据的可重现示例 (stackoverflow.com/questions/5963269/…)，以便其他人可以重新创建问题并修复它。

标签： r shapefile r-raster

【解决方案1】：

在不提供证据的情况下说“这行不通”并不是很有用。请提供一个独立的最小可重现示例。如果您不能这样做，请至少显示一些信息，如show(raster) 和show(vector)，并显示您使用的软件包library( )。最重要的是，显示错误消息。这可能对您不起作用的一个原因可能是矢量和栅格数据具有不同的坐标参考系。

如果没有，我只能说它确实工作，无论是 raster 或 terra 包。这里用terra说明：

library(terra)
v <- vect(system.file("ex/lux.shp", package="terra"))
r <- rast(system.file("ex/elev.tif", package="terra"))
v <- v[3,]
plot(r)
lines(v, lwd=2)

x <- crop(r, ext(v) + .01)
y <- mask(x, v)
plot(y)
lines(v)

【讨论】：