我有一个栅格，正尝试使用面裁剪它
这是我的多边形：
我的栅格文件包含分辨率为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)
为什么作物的论点在第一个例子中不起作用？“抢购”的论点的目的是什么？
1条答案
按热度按时间06odsfpq1#
crop的文档，具有SpatRaster参数x和SpatExtent（或另一空间对象的范围）参数y;这个论点“snap”是“near”、“in”或“out”之一。用于将y与x的几何对齐
y需要与x的几何图形对齐，因为您不能裁剪部分行或列。您只能裁剪整行或整列。因此您有三个选项，向内（较小区域）、向外（较大区域）或到最近的边界。