你发现的这个问题确实是典型的二进制浮点表示精度限制导致的——像5.3、7.1这类十进制小数没法被默认的FLT4S（32位单精度浮点）格式精确存储，写入文件再读回时就会显现出微小的偏差。你的缩放转整数方法是有效的，但还有更简便、更优雅的解决方案，也可以利用Terra的内置属性来处理这类问题：

方案1：直接使用更高精度的浮点格式写入 #

Terra默认用32位单精度浮点存储栅格，这是偏差的根源。你可以在writeRaster时指定datatype="FLT8S"（64位双精度浮点），这样就能极大降低甚至完全避免这类小数的精度损失：

library(terra)
r <- rast(ncol=2, nrow=2, vals=c(5.3, 7.1, 3, 1.2))
sprintf("%.8f", values(r)) # [1] "5.30000000" "7.10000000" "3.00000000" "1.20000000"

# 用64位双精度格式写入
t <- writeRaster(r, "test_double.tif", overwrite=TRUE, datatype="FLT8S")
sprintf("%.8f", values(t)) # [1] "5.30000000" "7.10000000" "3.00000000" "1.20000000"

这个方法不需要额外的缩放计算，操作最简洁，适合需要保留原始浮点格式的场景。唯一的小缺点是文件体积会比32位浮点大一倍，但在现代存储条件下这个代价基本可以忽略。

方案2：利用Terra的scale和offset元数据属性 #

Terra支持为栅格关联scale（缩放因子）和offset（偏移量）元数据，本质和你手动缩放的思路一致，但不需要手动修改栅格值，Terra会自动处理数值的存储与读取：

r <- rast(ncol=2, nrow=2, vals=c(5.3, 7.1, 3, 1.2))
# 设置缩放因子为10，偏移量为0（先确保元数据正确计算）
r <- setMinMax(r)
r@ptr$setScale(10)
r@ptr$setOffset(0)

# 写入时指定整数格式，Terra会自动应用缩放规则
t <- writeRaster(r, "test_scaled.tif", overwrite=TRUE, datatype="INT4S")
sprintf("%.8f", values(t)) # [1] "5.30000000" "7.10000000" "3.00000000" "1.20000000"

这种方法的优势很明显：

• 栅格值始终保持原始的浮点数值，不用手动做乘除缩放的操作
• 缩放和偏移的元数据会被写入文件，其他支持GDAL的工具（比如QGIS）读取时也能自动还原正确数值
• 同样能获得整数存储带来的文件体积优势

关于Terra是否支持精度属性的问题 #

Terra底层依赖GDAL库，GDAL确实支持为栅格添加类似STATISTICS_PRECISION的描述性元数据，但这类属性只是用来标注精度信息，没法直接解决浮点存储的精度损失问题。上面提到的「选用更高精度的浮点格式」和「利用scale/offset结合整数存储」，才是处理这类问题的实际有效手段。

总结一下：

• 如果追求最简便的操作，优先选datatype="FLT8S"写入双精度浮点
• 如果需要严格控制文件体积，推荐用scale/offset结合整数格式的方法，比手动缩放更优雅

内容的提问来源于stack exchange，提问作者Jean-Luc Dupouey