在R中以Shapefile矢量边界裁剪栅格的技术求助

嘿，我明白你的困扰！crop()函数确实只能按矩形范围裁剪，要实现矢量边界的精准裁剪，咱们得搭配mask()函数来完成——不管是用经典的raster包还是现在更高效的terra包都能搞定，下面给你一步步拆解：

方法一：使用raster包（经典方案） #

这是之前常用的工具链，兼容性很好：

• 第一步：加载依赖包并读取数据

# 加载所需包
library(raster)
library(rgdal) # 用于读取Shapefile

# 读取栅格数据（替换成你的文件路径）
raster_data <- raster("your_raster.tif")

# 读取Shapefile（dsn是文件夹路径，layer是shp文件名，不要加.shp后缀）
shape_data <- readOGR(dsn = "path/to/shapefile/folder", layer = "shapefile_name")

• 第二步：统一坐标系（关键！）
  如果你的栅格和矢量坐标系不一致，必须先转换：

# 将矢量数据转换为栅格的坐标系
shape_data <- spTransform(shape_data, crs(raster_data))

• 第三步：先缩小范围再精准裁剪
  先通过crop()缩小栅格范围（减少后续计算量，可选但推荐），再用mask()按矢量边界裁剪，超出边界的部分会被设为NA：

# 按矢量范围裁剪出矩形区域
cropped_raster <- crop(raster_data, extent(shape_data))

# 用矢量边界精准裁剪
masked_raster <- mask(cropped_raster, shape_data)

方法二：使用terra包（更高效的新方案） #

terra是raster包的替代工具，处理速度更快，语法更简洁：

• 第一步：加载包并读取数据

library(terra)

# 读取栅格
raster_terra <- rast("your_raster.tif")

# 读取Shapefile（直接传入shp文件路径即可）
shape_terra <- vect("path/to/your_shapefile.shp")

• 第二步：统一坐标系
  同样要确保两者坐标系一致：

# 转换矢量坐标系匹配栅格
shape_terra <- project(shape_terra, crs(raster_terra))

• 第三步：裁剪+精准 masking

# 先裁剪到矢量的矩形范围
cropped_terra <- crop(raster_terra, shape_terra)

# 按矢量边界精准裁剪
masked_terra <- mask(cropped_terra, shape_terra)

关键注意事项 #

• 坐标系必须一致：这是最容易踩坑的点，如果CRS不匹配，裁剪结果会完全错误
• 筛选单个要素：如果你的Shapefile里有多个矢量要素，只想裁剪其中一个，可以先筛选：

  # 比如筛选属性列"name"等于"Beijing"的要素（raster包示例）
  subset_shape <- shape_data[shape_data$name == "Beijing", ]
  

• 保存结果：裁剪完成后可以导出为tif文件：

  # raster包保存
  writeRaster(masked_raster, "cropped_masked_raster.tif", format = "GTiff", overwrite = TRUE)
  
  # terra包保存
  writeRaster(masked_terra, "cropped_masked_raster.tif", overwrite = TRUE)
  

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