#######################################################
## 使用矢量文件裁剪栅格数据
raster_data # 待裁剪的栅格数据
shapefile_data # 读入的shapefile文件
# 检查栅格数据和shapefile数据的投影是否一致
projection(raster_data)
projection(shapefile_data)
# 如果投影不一致,则需要将其中一个数据集投影到另一个的投影坐标系
shapefile_data_proj <- spTransform(shapefile_data, crs(raster_data))
# 获取shapefile数据的Extent
ext <- extent(shapefile_data_proj)
# 使用shapefile数据裁剪栅格数据
cropped_raster <- crop(raster_data, ext)
# 最后,使用mask()函数根据shapefile数据裁剪栅格数据
clipped_raster <- mask(cropped_raster, shapefile_data_proj)
clipped_raster
# 查看结果
plot(clipped_raster[[1]])