如何在R中重投影栅格时保持像元整数值不变？

你的问题出在projectRaster默认使用的插值方法上——默认的双线性插值（bilinear）会计算周边像元的加权平均值，自然会产生浮点数值。而你的土地利用数据属于类别型数据（像元值0-13是分类编码），应该使用最近邻插值，这种方法直接取距离最近的原像元值，能完整保留整数类别。

修改代码如下：

# 使用最近邻插值重投影，同时指定数据类型为整数
LU_AOI <- projectRaster(LU_AOI, crs = crs("EPSG:3005"), method = "ngb", datatype = "INT1U")

关键参数说明： #

• method = "ngb"：指定使用最近邻（nearest neighbor）插值，这是类别型栅格重投影的标准选择，不会改变原整数编码。
• datatype = "INT1U"：强制输出栅格为无符号1字节整数，匹配你原数据的整数范围（0-13完全在0-255范围内）。如果你的编码范围更大，可以选择对应类型，比如"INT2U"（无符号2字节）或"INT4S"（有符号4字节）。

额外注意： #

如果需要更精准控制输出栅格的分辨率、对齐方式，可以先创建一个目标栅格模板（比如基于已有同投影栅格），再用projectRaster的to参数匹配，示例：

# 假设已有一个3005投影的栅格target_raster作为模板
LU_AOI <- projectRaster(LU_AOI, to = target_raster, method = "ngb", datatype = "INT1U")

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