Raster_PostGIS_II

 

 

 

 

 

 

 

 

 

 

 

Vamos a continuar con la segunda parte del Tutorial para el manejo de datos raster en PostGIS que empezamos ya hace unas semanas – Tutorial Raster PostGIS I. Mientras que en la anterior entrada nos centramos en la carga de datos Raster en PostGIS hoy nos vamos a centra en el uso combinado de datos raster y vector mediante Geoprocesos en PostGIS.  Aunque dejaré para la descarga el tutorial con varios ejemplos de uso (vector+raster) me voy a centrar en explicar una sentencia en concreto que hace:

Calcula los valores de un raster sobre una capa de puntos pero con el filtro de un municipio concreto:

Para ello lo que haremos será devolver cada pixel que pertenezca al municipio de “Campeche”con la siguiente sentencia:

Explicado con más detalle…

1 – Seleccionamos los campos de la capa vector que queremos en el resultado

 

2 – Calculamos el valor del raster para cada punto concreto definido por la geometría del vector (geom) y lo guardaremos en un campo llamado i_deforestacion

 

3 – Establecemos la condición de que estén en el municipio de Campeche

 

Listo!!

PostGIS_Vector_Raster

En el tutorial que os dejamos en este enlace podréis ver cómo se construyen las siguientes consultas:

  • —  Metadata: extrae todos los valores de metadatos para el ráster “indice_deforestacion_p”
  • —  Histograma: estadísticas para todas las teselas (campo rid) del ráster en cuestión.
  • —  Loc_2010_rast_stats: agregación de nuevos campos a la capa loc_2010 con valores de deforestación y además las estadísticas de la tesela en la que cae cada punto.
  • —  Estados_buffer: se extrae la geometría ráster de un buffer de 200 metros del centroide de cada tesela del ráster, intersectando este con la capa estados.
  • —  Stats_clip_estado: obtiene las estadísticas de las teselas que intersectan con el estado de “Campeche“

Archivos a utilizar: Para el desarrollo del módulo vamos a utilizar diferente tipo de cartografía, tanto ráster como vectorial. —Los datos vectoriales de origen están en formato shapefile y los datos ráster en Gtiff. Toda la cartografía se encuentra en el sistema de referencia EPSG: 26916. Enlace a los archivos

Presentación: Presentación

Herramientas: Contamos con el programa PostgreSQL y su extensión espacial PostGIS para el desarrollo del módulo.