栅格重分类结果异常求助：QGIS与R操作均存数值/可视化问题

QGIS端异常原因及解决 #

问题本质 #

栅格计算器直接用"wc2.1_30s_elev@1" >= 2300输出的布尔结果，会被QGIS默认存储为带符号整型（如Int16），此时True被编码为-1、False为0，这就是你看到的奇怪负数。虽然识别工具会把-1解析为布尔真，但后续栅格计算会因数值类型异常出错。

解决方案 #

1. 显式将布尔结果转为1/0：

   ("wc2.1_30s_elev@1" >= 2300) * 1
   

2. 输出时指定数据类型为无符号整型（Byte或UInt16）：在栅格计算器的输出设置里选择对应类型，避免负数值产生。

R端异常原因及解决 #

核心问题 #

不管是索引赋值还是raster::reclassify，操作后栅格的元数据（值域min/max、数据类型）未自动更新，导致可视化和后续terra::xyFromCell读取错误的数值范围信息；另外重分类区间边界设置遗漏了临界值2300。

分步解决 #

1. 修正重分类区间
   原代码中2300这个临界值被漏掉，应调整区间包含该值：

   reclass_df <- c(0, 2300, 0,
                   2300, 6651, 1)
   reclass <- matrix(reclass_df, ncol = 3, byrow = TRUE)
   r_reclassif <- reclassify(r, reclass)
   

2. 强制更新栅格元数据
   操作后手动更新值域信息，确保工具读取正确的数值范围：

   # raster包
   r_reclassif <- setMinMax(r_reclassif)
   # terra包（推荐替代raster）
   r_reclassif <- terra::setMinMax(r_reclassif)
   

3. 显式指定数据类型
   将重分类后的栅格转为整型，避免浮点型精度问题：

   # raster包
   r_reclassif <- as.integer(r_reclassif)
   # terra包更简洁的重分类方式
   library(terra)
   r <- rast("wc2.1_30s_elev.tif")
   r_reclass <- ifel(r >= 2300, 1, 0)
   r_reclass <- as.integer(r_reclass)
   r_reclass <- setMinMax(r_reclass)
   

4. 修复可视化问题
   直接指定色阶断点，确保0和1被正确区分：

   plot(r_reclassif, col = c("lightgray", "darkgreen"), breaks = c(-0.5, 0.5, 1.5))
   

5. 确保terra::xyFromCell正常工作
   先检查栅格的实际值是否正确：

   unique(values(r_reclassif))
   

   如果存在NA或异常值，需在重分类时保留NA：

   # terra包处理NA
   r_reclass <- ifel(r >= 2300, 1, ifel(r < 2300, 0, NA))
   

总结 #

所有异常的核心诱因是：

• 栅格数据类型不匹配导致布尔结果错误编码（QGIS）
• 重分类后元数据未更新，导致工具读取错误的数值范围（R）
• 重分类区间边界设置遗漏临界值

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