¿Cómo etiquetar las celdas raster en QGIS?

Question

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

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

Mi idea es usar gdal2xyz, cargar dicho archivo usando el plugin CSV y luego guardarlo como archivo vectorial. Por desgracia, es mucho trabajo, porque tengo muchos archivos de trama.

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

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

Answer 1

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

La mejor opción es utilizar un archivo vectorial convertido para etiquetar los centros de las celdas. El paquete estadístico de código abierto R tiene muchas herramientas espaciales excelentes y puede procesar por lotes o rápidamente archivos raster para crear shapefiles (vectores/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 shapefile tendrá una columna con los valores ráster de cada punto. Los puntos se pueden cargar en QGIS, dándoles un tamaño de 0 y etiquetándolos adecuadamente; aparecen en el centro de las celdas.

Para recorrer, por ejemplo, todos los archivos TIF de 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.
}

Answer 2

En caso de que R no esté disponible (o, quiera hacerlo dentro del entorno de QGIS), ahora existe la herramienta SAGA en la caja de herramientas de procesamiento, Valores rasterizados a puntos para extraer los valores de las celdas.

Podemos etiquetar la salida Shapes y superponerla a la trama original.