使用多边形裁剪栅格不选择 terra 包中的所有单元格

Question

我有一个栅格，我正在尝试使用多边形对其进行裁剪

这是我的多边形：

我的光栅文件包含 0.25 度分辨率的整个地球，所以我使用多边形裁剪了光栅：

 library(terra)
 crop_raster <- terra::crop(temp_rast, county_ref)  
 plot(crop_raster)
 plot(county_ref, add = T)

这种裁剪看起来很奇怪，因为只选择了 2 个栅格单元格。我期待像下面这样的东西，我会使用

snap = "out"

参数。

crop_raster_snap <- terra::crop(temp_rast, county_ref, snap = "out") 
plot(crop_raster_snap)
plot(county_ref, add = T)

为什么作物论证一开始没有奏效？ “快照”论点的目的是什么？

Answer 1

crop

的文档，带有 SpatRaster 参数

x

和 SpatExtent（或另一个空间对象的范围）参数

y

；那个论点“啪”是

“近”、“内”或“外”之一。用于将 y 与 x

的几何对齐

y

需要与

x

的几何形状对齐，因为您不能裁剪部分行或列。您只能裁剪整行或整列。所以你有三个选择，进入内部（较小的区域），外部（较大的区域）或最近的那个。