基于R解决栅格影像行列数不一致的问题

当我们使用Matlab对长时间序列的影响做偏相关分析时，如果各影像的像元大小不匹配，会出现下面的这种报错。

解决方法：

根据列数的大小使用arcgis迭代工具或批量掩膜对影响进行提取，此时需要注意更改环境变量的参数，基于一个影像来统一大小。

使用代码

R语言，对于文件路径，三四行是参照的路径，七八是需要更改的影响的路径，10是输出位置。注意，文件夹内的影响像元大小必须相同。

library(tidyverse)
library(terra)
library(projects)
ras1=list.files("D:\\xiangyuan\\E\\") %>% 
  paste0("D:\\xiangyuan\\E\\",.) %>% rast()

ras4=list.files("E:\\tmp_Peng\\resample\\mask\\new\\",pattern=".tif$") %>% 
  paste0("E:\\tmp_Peng\\resample\\mask\\new\\",.) %>%
  rast() %>% project(.,ras1[[1]])
writeRaster(ras4,filename = paste0("D:\\xiangyuan\\tmp_end\\",
                                   names(ras4),".tif"))

如果您遇到问题或者下载处理遥感数据等方面的帮助，您可以评论或者私信。

最后，谢谢大家点点赞和关注~