RStudio中为SpatRaster设置标签报错及栅格分类问题

一、tmap自定义标签报错修复 #

你遇到的invalid 'labels'; length 3 should be 1 or 2错误，核心原因是：

• 你的ski_suit栅格中，NA被设为缺失值，tmap默认不会将NA作为一个分类类别，实际仅识别到1和2两个有效类别
• 你传入了3个标签，长度与实际类别数不匹配

同时原代码存在语法错误：labels数据.frame的括号未闭合，需补充)。

修正后的代码 #

场景1：不显示NA类别

ski_suit <- dem_e*slope_e*distance_e
ski_suit[ski_suit<=5] <- NA
ski_suit[ski_suit>=6&ski_suit<=9] <- 1
ski_suit[ski_suit>9] <- 2

# 修正语法错误：闭合数据.frame括号
labels <- data.frame(
  values=c(NA,1,2),
  labels=c('NA','Not Suitable','Suitable')
)

# 仅传入2个标签，对应1和2
tm_shape(ski_suit)+
  tm_raster(style='cat', labels=c('Not Suitable','Suitable'), palette='-BrBG')

场景2：显示NA类别

如果需要显示NA的标签，需开启showNA=TRUE，同时传入包含NA的标签向量：

tm_shape(ski_suit)+
  tm_raster(style='cat', 
            labels=c('Not Suitable','Suitable','NA'), 
            palette='-BrBG',
            showNA=TRUE)

二、SpatRaster重分类错误修复 #

你遇到的两个重分类错误原因如下：

1. terra包的classify函数需要rcl参数定义重分类规则，而非单独的breaks和labels
2. reclassify是raster包的专属函数，不适用于terra包的SpatRaster对象，需使用terra的classify

修正后的重分类代码 #

library(terra)

# 定义重分类矩阵：每行格式为[起始值, 结束值, 新值]
rcl_matrix <- matrix(c(
  6, 9, 1,   # 6-9区间重分类为1
  9, Inf, 2  # >9区间重分类为2
), ncol=3, byrow=TRUE)

# 执行重分类
ski_rcl <- classify(ski_suit, rcl=rcl_matrix)

# 为栅格添加标签
levels(ski_rcl) <- data.frame(
  ID=c(1,2),
  labels=c('Not Suitable', 'Suitable')
)

# 若需包含NA标签，可追加
levels(ski_rcl) <- rbind(levels(ski_rcl), data.frame(ID=NA, labels='NA'))

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