Terra包中栅格按不规则多边形高效裁剪方法问询

你当前的代码逻辑可行，但可以通过调整步骤顺序和简化操作来大幅提升效率——先裁剪到多边形的外接范围，再执行遮罩，这样能减少mask函数需要处理的像素数量，尤其针对大尺寸栅格时效果更明显。

优化后的代码示例 #

# 读取数据
kml <- vect("pathtokml.kml")
r <- rast("pathtoraster.tif")

# 统一CRS：将矢量投影到栅格的坐标系
kml_proj <- project(kml, r)

# 先裁剪栅格到矢量外接矩形，再用矢量遮罩保留内部像素
crp_final <- mask(crop(r, kml_proj), kml_proj)

# 可视化结果
plot(crp_final)

关键优化点说明 #

• 步骤顺序调整：先执行crop把栅格缩小到矢量的外接矩形范围，再用mask保留多边形内部像素。相比先对全幅栅格做mask再裁剪，这种方式能避免大量不必要的像素计算，直接降低运算负载。
• 避免重复计算：将投影后的矢量对象存为独立变量kml_proj，避免在crop和mask中重复执行project操作，节省系统资源。

额外提速技巧 #

• 若KML包含多个多边形要素，可先用aggregate(kml_proj)合并为单个要素，减少遮罩时的计算开销。
• 针对多波段栅格，可通过subset(r, c(1,3))这类语法先筛选需要处理的波段，进一步降低计算量。

内容的提问来源于stack exchange，提问作者Andrew