PostGIS intersección y resumir los atributos

Question

Normalmente soy un usuario de ArcGIS desktop, pero con ganas de empezar a usar PostGIS más y tengo una gran bits de procesamiento de hacer. No estoy seguro de qué funciones a utilizar, espero que alguien pueda ayudar.

Tengo un conjunto de datos poligonales (de varios millones de características), basada en un tipo de uso de la tierra/ cubierta vegetal classfication (20 categorías). Tengo un número de regiones en otro conjunto de datos.

Para cada una de las regiones, me gustaría saber el área de cada landcover classfication.

En ArcGIS (si se trataba de un conjunto de datos más pequeño) me imagino que la primera adición de la región para cada uno de los polígonos en la tabla de atributos mediante una combinación. A continuación, el uso de "resumir" en la mesa por región y por la cubierta vegetal de la clasificación.

No está seguro de por dónde empezar a hacer esto en PostGIS / SQL.

Actualización:

Wow gracias por lo que ha sido una gran ayuda.

Ha estado corriendo un largo período de tiempo (44 horas!) y ahora me sale:

NOTICE:  TopologyException: found non-noded intersection between LINESTRING (coords
edited) and LINESTRING (coords edited) at (position edited)
ERROR:  GEOS Intersection() threw an error!

********** Error **********

ERROR: GEOS Intersection() threw an error!
SQL state: XX000

Supongo que esto es un problema en los datos originales - sólo un caso de la revisión de los datos originales o puedo comprobar primero la topología de alguna forma por el conjunto de los datos? Hay algo acerca de la aceptación de ciertos errores de procesamiento / tolerancias?

Answer 1

Suponiendo que tiene el siguiente diseño de tabla

landcover(id,type,the_geom)
region(id,name,the_geom)

Valores de área por tipo de ambientes por región pueden calcularse utilizando

SELECT r.id, r.name, r.the_geom, l.type, 
       Sum(ST_Area(ST_Intersection(l.the_geom,r.the_geom)))
FROM landcover l,
     region r
WHERE l.the_geom && r.the_geom
GROUP by r.id, r.name, r.the_geom, l.type

ST_Intersection se utiliza para tener en cuenta para polígonos de ambientes que son sólo parcialmente dentro de una región.