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]])
## 使用矢量文件裁剪栅格数据
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]])
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】博客园社区专享云产品让利特惠,阿里云新客6.5折上折
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· DeepSeek “源神”启动!「GitHub 热点速览」
· 我与微信审核的“相爱相杀”看个人小程序副业
· 微软正式发布.NET 10 Preview 1:开启下一代开发框架新篇章
· 如何使用 Uni-app 实现视频聊天(源码,支持安卓、iOS)
· C# 集成 DeepSeek 模型实现 AI 私有化(本地部署与 API 调用教程)