R中替代rgdal包导出带颜色表的单波段栅格TIFF方法

terra是rgdal和raster包的官方继任者，对栅格颜色表的支持更完善，导出的TIFF能被QGIS直接识别颜色映射。

实现步骤 #

• 用rast()读取栅格文件
• 设置缺失值（对应原mvFlag参数）
• 定义颜色表与类别名称并绑定到栅格
• 用writeRaster()导出TIFF

示例代码：

# 加载terra包
library(terra)

# 读取栅格
r <- rast("my_raster.tif")

# 设置缺失值为0（对应原代码的mvFlag=0）
NAflag(r) <- 0

# 定义颜色表和类别名称（顺序对应栅格值0到5）
color_table <- c("#00000000", "#FF0000", "#FFA500", "#7FFF00", "#006300", "#0000ff")
cat_names <- c("class_a", "class_b", "class_c", "class_d", "class_e", "class_e")

# 为栅格绑定颜色表与类别信息
coltab(r) <- data.frame(value = 0:5, color = color_table, name = cat_names)

# 导出TIFF，指定数据类型为Byte（对应原type="Byte"）
writeRaster(r, "my_new_raster.tif", datatype = "INT1U", format = "GTiff", overwrite = TRUE)

如果仍依赖raster包，需直接在writeRaster()中指定colorTable参数（修改@legend插槽仅用于R内可视化，不会写入文件）。

示例代码：

# 加载raster包
library(raster)

# 读取栅格
r <- raster("my_raster.tif")

# 定义颜色表（顺序对应栅格值0到5）
color_table <- c("#00000000", "#FF0000", "#FFA500", "#7FFF00", "#006300", "#0000ff")

# 导出时直接传入colorTable参数
writeRaster(r, "my_new_raster.tif", datatype = "INT1U", format = "GTiff", 
            colorTable = color_table, overwrite = TRUE)

# 若要保留类别名称，可生成QGIS识别的辅助文件
cat_names <- c("class_a", "class_b", "class_c", "class_d", "class_e", "class_e")
writeLines(paste0(0:5, ":", cat_names), "my_new_raster.tif.aux.xml")

注意事项 #

• terra包的方法更可靠，颜色表直接嵌入TIFF元数据，无需额外辅助文件
• 确保栅格为整数类型（Byte/INT1U），与颜色表的离散类别匹配
• QGIS打开导出的TIFF时会自动加载颜色映射，无需手动配置

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