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

最新推荐文章于 2024-06-09 23:44:15 发布

云兮001

最新推荐文章于 2024-06-09 23:44:15 发布

阅读量2.6k

点赞数 3

文章标签： r语言

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

版权

本文介绍了如何利用R语言的rgdal和raster库来裁剪和处理大量栅格数据，通过读取shp矢量图文件和tif栅格文件，进行投影匹配和裁剪操作，然后对结果进行可视化和NA值处理，以便进一步分析如平均值、最大值和最小值等。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

由于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文件都可以进行处理。无论大的小的。谢谢您的阅读。