R语言学习笔记（三）R语言处理遥感影像

R语言批量加载遥感影像

tif_file_path <- list.files(r"(D:\degree\path\paper-jinxing\paper1\invest)", pattern = ".tif$", full.names = TRUE, ignore.case = TRUE)#设置路径
tif_file_path #输出路径
tif_file_all <- stack(tif_file_path) #将全部栅格数据放入同一变量
plot(tif_file_all)#绘制图像 但是所有的图像必须是相同边界大小

如图：

tif_all_mean <- mean(tif_file_all, na.rm = TRUE)
#计算像元的平均值，在tif_file_all后面不加[]，的话，相当于在计算每个影像自己的mean,
tif_all_mean <- mean(tif_file_all[], na.rm = TRUE)
#加了[]之后，相当于计算所有影像的mean,效果完全不一样

效果如下图：

tif_file_all <- stack(tif_file_path[c(1)])
#表示画出数据集的第一个数据

效果如图