R语言shp矢量图裁剪栅格tif(同arcgis掩膜裁剪）

最新推荐文章于 2023-06-27 19:45:44 发布

云兮001

最新推荐文章于 2023-06-27 19:45:44 发布

阅读量1.9k

点赞数 3

文章标签： r语言 Powered by 金山文档

本文链接：https://blog.csdn.net/m0_53956136/article/details/129474738

版权

由于arcgis对于大量栅格数据的裁剪速度较慢而且数据导出也比较麻烦，在这里我运用R语言进行shp矢量图裁剪栅格文件。下面列出矢量图裁剪一个栅格文件的代码，对于多个裁剪，参照我前面写的几篇博客进行组合即可。

library("rgdal")
library("raster")
myshp <- readOGR("languanhe.shp")
myraster <- raster("X1979.01.01.10.30.tif")
projection(myshp) <- CRS("+proj=longlat +datum=WGS84 +no_defs")
masked <- mask( crop(myraster,myshp, snap="out"), mask =myshp)
plot(masked)
c <- masked@data@values

这样就可以提取出裁剪的数据了，对于里面有NA的数据可以直接用na.omit函数对na值进行去除，在这里na值就是没有数据的值。这样就可以对提取的数据进行，平均值，最大值，最小值等等的处理。在这里对于nc文件的处理方法已经完结。对于任何nc文件都可以进行处理。无论大的小的。谢谢您的阅读。