¿Cómo etiquetar celdas ráster en QGIS?

Intento mostrar un archivo ráster de resolución pequeña en la herramienta QGIS Print Manager. Me gustaría mostrar una pequeña etiqueta con cada valor de celda (posición, alrededor del centro de cada celda).

¿Cómo puedo hacerlo? ¿Hay algún complemento que pueda ayudarme a realizar esta tarea?

Mi idea es usar gdal2xyz, cargar dicho archivo usando el complemento CSV y luego guardarlo como un archivo vectorial. Desafortunadamente es mucho trabajo, porque tengo muchos archivos raster.

Otra idea es usar gdal2xyz.py, luego ogr2ogr para guardarlo como archivo vectorial.

¿Hay alguna forma de agrupar este proceso? ¿Hay algún programa de código abierto donde pueda hacerlo mucho más fácilmente?

No parece haber un complemento o funcionalidad en QGIS para etiquetar celdas.

El uso de un archivo vectorial convertido para etiquetar los centros de las celdas es probablemente su mejor opción. El paquete de estadísticas de código abierto R tiene muchas herramientas espaciales excelentes y puede procesar por lotes o procesar rápidamente archivos ráster para crear archivos de forma (vector / puntos).

library(maptools)
library(raster)

# Load the raster from a file
r <- raster("/workspace/TEMP/raster.asc")

# Convert to spatial points
p <- as(r, "SpatialPointsDataFrame")

# Save as a shapefile
writeSpatialShape(p, "/workspace/TEMP/raster_points")

El archivo de forma tendrá una columna con los valores ráster para cada punto. Los puntos se pueden cargar en QGIS, dado un tamaño de 0 y etiquetados adecuadamente; Aparecen en el centro de las células.

Para revisar, por ejemplo, todos los archivos TIF en un directorio:

for (file in dir("/workspace/TEMP/", pattern="*.tif")) { # list all .tif files
  r <- raster(file)
  p <- as(r, "SpatialPointsDataFrame")
  writeSpatialShape(p, substr(file, start = 1, stop = nchar(file) -4)) # substr() removes extension.
}

En caso de que R no esté disponible (o desee hacerlo dentro del entorno QGIS), ahora hay una herramienta SAGA en la caja de herramientas de procesamiento, valores ráster a puntos para extraer valores de celda.

Podemos etiquetar la Shapescapa de salida y superponerla con el ráster original.