Interpolación de datos ráster con PostGIS

Question

Actualmente estoy utilizando CDO para extraer datos de archivos netCDF que contienen datos meteorológicos por satélite. En el archivo netCDF los puntos tienen una resolución de 0,5 grados, así que cuando se piden datos para una coordenada específica x,y, uso cdo -remapbil,lon=x/lat=y para generar un archivo netCDF temporal del que extraer los datos relativos a x,y. Ahora, a medida que se añaden nuevos datos, el proceso es cada vez más lento.

Buscando nuevas formas de realizar la tarea, veo que PostGIS puede guardar datos raster (de hecho he podido importar algunos datos, convirtiendo netCDF a GeoTIFF y luego usando raster2pgsql para importar los datos).

Ahora estoy atascado porque veo que el raster se importa como píxeles y no sé cómo formular una consulta para una coordenada x,y genérica que devuelva los mismos datos que con CDO. ¿Cómo lo afronto?

Answer 1

SELECT ST_Value(rast, ST_SetSRID(ST_MakePoint(X, Y), ST_SRID(rast))) val
FROM rasterTable;

Si ha importado su trama sin matizarla (con la opción -t) o:

SELECT ST_Value(rast, ST_SetSRID(ST_MakePoint(X, Y), ST_SRID(rast))) val
FROM rasterTable
WHERE ST_Intersects(rast, ST_SetSRID(ST_MakePoint(X, Y), ST_SRID(rast)));

Si al importar el ráster lo has hecho en forma de mosaico.