如何在terra R包中实现raster::deratify的等效功能？

terra包中没有和raster::deratify完全对应的函数，但可以通过分类栅格的属性表操作结合值替换/重分类实现相同需求，以下是完整实现方案：

1. 构建带RAT的分类栅格 #

先创建基础栅格并添加属性表，步骤和raster包类似：

library(terra)
# 创建10x10的基础栅格
r <- rast(nrow=10, ncol=10)
values(r) <- 1
r[51:100] <- 2
r[3:6, 1:5] <- 3

# 转为分类栅格并添加RAT
r <- categories(r)
rat <- levels(r)[[1]]
rat$landcover <- c("Pine", "Oak", "Meadow")
rat$code <- c(0.002, 0.0056, 0.0124)
levels(r) <- rat

2. 基于RAT指定列生成新栅格 #

这里提供两种简洁的实现方式：

方法一：使用subst()进行值替换（推荐） #

通过构建"原分类值→目标列值"的映射规则，直接替换栅格值：

# 构建替换映射：原分类值对应code列的数值
replace_map <- setNames(rat$code, rat$value)
# 执行替换生成新栅格
x <- subst(r, from = as.integer(names(replace_map)), to = replace_map)

# 验证结果：提取指定单元格的值
extract(x, c(1, 2, 25, 100))

方法二：使用classify()进行重分类 #

通过分类矩阵定义值的映射关系，适合更复杂的范围类重分类场景：

# 构建分类矩阵：[原值起始, 原值结束, 目标值]
class_matrix <- cbind(rat$value, rat$value, rat$code)
# 执行重分类生成新栅格
x <- classify(r, class_matrix)

# 验证结果
extract(x, c(1, 2, 25, 100))

核心注意点 #

• terra中categories()对应raster的ratify()，用于将普通栅格转为带RAT的分类栅格
• 查看或修改RAT仍使用levels()函数，用法与raster包保持一致
• subst()更适合一对一的精确值替换，是该需求下最直接的实现方式

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