R使用矢量文件裁剪栅格数据

#######################################################
## 使用矢量文件裁剪栅格数据
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]])


posted @   junjunang  阅读(97)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· DeepSeek “源神”启动!「GitHub 热点速览」
· 我与微信审核的“相爱相杀”看个人小程序副业
· 微软正式发布.NET 10 Preview 1:开启下一代开发框架新篇章
· 如何使用 Uni-app 实现视频聊天(源码,支持安卓、iOS)
· C# 集成 DeepSeek 模型实现 AI 私有化(本地部署与 API 调用教程)
点击右上角即可分享
微信分享提示