为何月度均值栅格输出值被莫名放大10000倍？

问题根源 #

输出结果被放大10000倍的核心原因是**writeRaster函数默认的输出格式设置**：原TIFF文件是0-1的浮点型数据，但写入新TIFF时，函数自动选用了整型数据类型（比如16位无符号整型），为适配整型的取值范围，默认把0-1的数值放大10000倍存储。

解决步骤 #

1. 验证计算结果的正确性 #

在stackApply之后，先检查月度均值的数值范围，排除计算环节的问题：

# 查看均值结果的数值范围
cellStats(alb_mean, range)

如果输出显示范围在0-1左右，说明问题完全出在写入TIFF的步骤。

2. 修改writeRaster参数 #

写入文件时明确指定输出为浮点型，关闭自动缩放逻辑，替换原有循环代码：

for(n in names(alb_mean)){
  writeRaster(alb_mean[[n]], 
              paste0(n, ".tif"),
              datatype = "FLT4S",  # 指定32位浮点型，匹配原数据类型
              scale = 1,           # 强制缩放系数为1，不放大数值
              offset = 0,          # 强制偏移量为0
              overwrite = TRUE,    # 允许覆盖已存在文件
              options = c("COMPRESS=LZW")  # 可选，压缩文件节省空间
  )
}

• datatype="FLT4S"：确保输出为浮点型，避免整型转换导致的数值放大；
• scale=1和offset=0：明确禁止函数对数值做任何缩放调整，保留原始计算结果。

3. 可选：检查原数据元数据 #

若想确认原TIFF的属性，可查看单个文件的信息：

r <- raster(rast[1])
print(r)  # 查看数据类型、scale/offset等元数据
cellStats(r, range)  # 确认原数据的数值范围

如果原数据本身带有scale元数据，读取时raster包会自动转换为真实值，但写入时若不指定参数，可能反向应用缩放导致数值放大。

优化建议 #

尽量避免使用setwd，直接在函数中使用完整路径更稳妥，比如list.files(path="D:/Data", ...)和writeRaster(..., path="D:/Data/Output")，防止路径混乱。

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