R语言之批量重分类—1

R语言作为统计学中非常重要的语言，随着计量地理学的兴起……，算了不多说这么多没用的，总之随着地理空间数据分析需求量的增大，R语言对GISers已经越来越重要。

关于R语言的介绍和下载就不多说。这篇文章主要介绍一下如何读取矢量数据（shp数据格式等……），栅格数据（以tiff数据为例）。

当然最重要的是，如何批量处理数据，减少我们的工作量，节约我们宝贵的时间，解放双手……

1、读取矢量和栅格数据  (R语言还是很简洁的)

# 读取矢量数据
shp<- vect("G:/MaoW/Watersheds/Watersheds.shp")
# 读取栅格数据
rast <- rast("G:/MaoW/jsdata/rast.tif")

2、进行可视化，R语言中进行简单的可视化非常便捷，只需要一行代码,这里进行可视化的是全国的NDVI数据。

plot（rast）

3、批量重分类，重分类是ArcGIS中非常常用的工具，如何快捷方便的进行批量处理呢？R语言可以很轻松的完成这项任务。正好最近又重新制作了赣州市近20年的产水深度数据，就拿它来实验。这是他一开始的样子（局部）。

基本流程：先读取文件夹下的需要重分类的多个栅格数据，存储在列表中；再进行重分类，重分类包括建立矩阵，作为重分类的参数。代码如下:

# 引入相应的包
library(terra)
library(trend)
# 存入列表中
WYield = rast(fs::dir_ls("G:/MaoW/r1/", regexp = ".tif$")[1:21])
#构造一个矩阵
m = c(0,300,1, 300,600,2, 600,1200,3, 1200,1500,4, 1500,2200,5)
m = matrix(m,ncol = 3,byrow = T)
# 开始执行重分类
rcWYield = classify(WYield,m,right = T)
#进行可视化
plot(rcWYield)

下面是重分类后的结果：

4、其实，R语言在可视化方面也是很好用的，熟练掌握的话可以和轻松的画出好看的地图。

比如换一种方式画上面的图

技术学习需要日积月累，长夜漫漫，愿你我同行，我是加拿大一枝黄花，感谢关注梧桐凉月！

梧桐凉月

GIS与诗词，分享地理信息知识、诗词，每天学习一点小技巧

97篇原创内容

公众号