如何在terra栅格的指定位置嵌入指定文本？

可以利用terra包的绘图系统结合R基础绘图的text()函数实现，步骤简单直接：

• 1. 准备栅格数据
  先加载terra包，创建或读取目标SpatRaster对象：

  library(terra)
  # 示例：创建一个随机栅格
  r <- rast(ncol=10, nrow=10)
  values(r) <- runif(ncell(r))
  

• 2. 在指定位置嵌入文本
  先绘制栅格，再用text()函数指定坐标位置添加文本。坐标可以是栅格的地理坐标，也可以通过行列号转换得到：

  # 绘制栅格
  plot(r)
  # 示例：在地理坐标(x=0, y=0)位置添加文本 "Sample Text"
  text(x=0, y=0, labels="Sample Text", col="red", cex=1.2, font=2)
  
  # 如果知道目标单元格的行列号，先转换为地理坐标
  cell_coord <- xyFromCell(r, cell=55) # 获取第55个单元格的坐标
  text(x=cell_coord[,1], y=cell_coord[,2], labels="Cell 55", col="blue", cex=1)
  

  参数说明：

  • x/y: 文本放置的坐标位置
  • labels: 要嵌入的指定文本
  • col: 文本颜色
  • cex: 文本大小
  • font: 文本字体样式（1=常规，2=粗体，3=斜体，4=粗斜体）

• 3. 保存带文本的栅格图像
  如果需要保存嵌入文本后的可视化结果，可使用R的图形设备函数：

  png("raster_with_text.png", width=800, height=600)
  plot(r)
  text(x=0, y=0, labels="Sample Text", col="red", cex=1.2)
  dev.off()
  

若要在栅格文件的元数据中嵌入文本（非可视化层），可使用terra::writeRaster()的options参数添加GDAL元数据：

writeRaster(r, "raster_with_metadata.tif", options=c("TIFFTAG_IMAGEDESCRIPTION=Embedded Text Here"))

这种方式的文本会存储在栅格文件的属性信息中，而非图像显示层。

内容的提问来源于stack exchange，提问作者Jean-Luc Dupouey