关于使用R语言堆叠Landsat影像特定TIFF文件的技术问询

Hey there! 你已经成功筛选出了Landsat 5场景中需要的B1-B5、B7波段TIFF文件，接下来我会一步步带你完成栅格堆叠的实现，同时提供两种主流工具包的方案供你选择：

1. 准备必要的R包 #

处理栅格数据我们常用terra（raster包的升级版本，性能更优）或者经典的raster包，先把它们准备好：

# 首次使用先安装包
install.packages(c("terra", "raster"))

# 加载你需要的包（二选一即可，推荐terra）
library(terra)
# library(raster)

2. 执行堆叠操作 #

方案一：使用terra包（推荐） #

terra的rast()函数可以直接读取你筛选好的文件列表，一键生成栅格栈：

# 基于已筛选的文件列表创建栅格栈
may2006_stack <- rast(may2006)

# 查看堆叠后的栅格信息，确认波段数量、投影等是否正确
may2006_stack

方案二：使用raster包 #

如果你更习惯用经典的raster包，代码如下：

library(raster)

# 先逐个读取每个波段的TIFF文件
may2006_rasters <- lapply(may2006, raster)

# 将单个栅格对象堆叠成一个栅格栈
may2006_stack <- stack(may2006_rasters)

# 查看堆叠结果信息
may2006_stack

3. 验证堆叠结果 #

完成堆叠后，可以通过以下方式确认结果是否符合预期：

• 检查波段数量：使用nlyr(may2006_stack)（terra）或nlayers(may2006_stack)（raster），正常应该返回6（对应B1-B5、B7）
• 查看波段名称：输入names(may2006_stack)，确认每个波段的标识和文件名对应
• 快速可视化：执行plot(may2006_stack)，查看各波段的影像是否正常显示

4. 保存堆叠后的文件 #

如果需要把堆叠好的栅格保存为新的TIFF文件，执行以下代码即可：

# terra包写法
writeRaster(may2006_stack, "may2006_stacked.tif", overwrite = TRUE)

# raster包写法
writeRaster(may2006_stack, "may2006_stacked.tif", overwrite = TRUE)

小提示 #

• Landsat同一场景的各波段文件默认是同投影、同分辨率、同范围的，所以堆叠一般不会出错，但如果遇到报错，优先检查文件名是否完整（比如你列表里最后一个B7文件显示不全，确认文件名没有截断）
• 如果后续需要做波段运算（比如计算NDVI），栅格栈可以直接作为输入使用

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