terra包批量读取多栅格文件生成多层spatRaster对象的方法

最简便的零循环写法 #

完全不需要手动写循环逐文件加载再拼接，rast()函数本身就支持直接传入栅格文件路径向量，自动完成多文件读取、网格校验、拼接为多层spatRaster的操作，完全对应raster中stack(filelist)的功能：

temp_rasters <- rast(filelist_temp)

这是迁移过程中最推荐的实现。相比raster包的stack，这套实现全程基于磁盘引用读取，不会提前把所有栅格像元值加载到内存，处理大体积栅格的效率会高很多，输出的多层栅格层顺序和输入的文件列表顺序完全一致。

自定义写法的问题说明 #

你之前尝试运行的

temp_rasters <- c(lapply(filelist_temp, rast))

没有得到预期多层栅格的核心原因是：lapply返回的是存放了25个独立spatRaster对象的R列表，而terra的c()函数默认不会自动拆解列表中的元素，直接传入完整列表时，最终会生成一个SpatRasterCollection类型的容器对象——这个对象只是零散存放各个spatRaster的集合，并不是可以直接参与栅格计算的多层栅格栈。

基于已生成spatRaster列表的拼接方法 #

如果你的流程中已经产出了独立的spatRaster列表（比如不是从磁盘直接读取的文件，而是中间处理环节生成的栅格），不需要回退修改文件读取逻辑，可以用do.call()把列表中的spatRaster逐个拆解作为参数传给c()，即可得到正确的多层栅格：

temp_rasters <- do.call(c, lapply(filelist_temp, rast))

迁移过程的常见注意点 #

• 两种拼接方法得到的多层spatRaster，后续操作逻辑和raster包的stack基本一致，比如读取像素值、执行分层运算、赋值层名等操作都可以直接对应迁移
• terra在拼接栅格时，会严格校验所有输入栅格的坐标系、空间范围、分辨率、像元对齐情况，只要有任何不匹配就会直接抛出错误，不会像raster包那样自动重采样对齐到统一网格，这是迁移过程中最容易踩的逻辑差异点
• 拼好的多层栅格可以直接用names(temp_rasters) <- 自定义层名向量修改层名，操作逻辑和raster完全一致

内容的提问来源于stack exchange，提问作者canderson156