如何用terra包将SpatRaster替换后的值导出为ASCII格式.txt文件？

问题1：导出SpatRaster的实际替换值 #

你的核心问题是分类栅格（categorical SpatRaster）在导出时默认存储类别ID，而非显示的分类值。使用subst()替换分类名称后，栅格仍保持分类属性，新的替换值会被作为新类别，底层存储的是自动分配的新ID（比如你输出里的4、5、6）。解决方法是将分类栅格转换为普通数值型栅格，让单元格直接存储替换后的数值：

修正后的代码 #

library(terra)

# 创建测试数据
r <- rast(nrows=3, ncols=3)
values(r) <- sample(3, ncell(r), replace=TRUE)
cls <- data.frame(id=1:3, cover=c("forest", "water", "urban"))
levels(r) <- cls

# 定义替换规则
replacement <- c(10,20,30)
# 创建类别名称到替换值的映射
cover_map <- setNames(replacement, cls$cover)

# 提取每个单元格的类别名称，替换为目标数值
cell_values <- cover_map[as.character(as.factor(values(r)))]
# 创建新的数值型栅格（去掉分类属性）
r_replaced <- r
values(r_replaced) <- as.integer(cell_values)

# 导出为AAIGrid格式的txt文件
writeRaster(r_replaced, "test.txt", filetype = "AAIGrid")

原理说明 #

• 通过映射表将原分类名称直接转为目标整数
• 赋值给新栅格后，栅格不再保持分类属性，底层存储的就是实际替换值
• 导出时writeRaster会直接写入这些数值，而非类别ID

问题2：导出为.txt文件的注意事项 #

• 必须指定filetype = "AAIGrid"：这是ASCII栅格的标准格式，writeRaster会自动生成正确的文件头（ncols、nrows等）和数值布局，仅改后缀为.txt不影响格式正确性，但必须确保参数正确。
• 栅格数据类型：目标模型要求整数的话，需将栅格转为整数型（as.integer()），避免导出时出现小数或科学计数法。
• NODATA值：默认NODATA值是255，若目标模型要求其他值，可通过NAflag参数指定（比如writeRaster(..., NAflag = -9999)）。
• 坐标信息：文件头的xllcorner、yllcorner等继承自原栅格，需确认这些坐标符合目标模型的坐标系要求，必要时先转换栅格坐标系。
• 文件名兼容性：虽然标准AAIGrid后缀是.asc，但目标模型要求.txt的话完全可行，只要格式结构正确即可。

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