Uso de postgis para aplanar todos los polígonos de la tabla

Question

Tengo una tabla con muchos polígonos. Todos los polígonos tienen un nombre y un valor.

Quiero aplanar estos polígonos para que cada borde de los polígonos permanezca como borde de los polígonos resultantes. En lugar de tener una pila de polígonos me quedaré con un mosaico de polígonos basados en las mismas líneas.

Resultando en esto:

Los nombres son insignificantes. Quiero sumar los valores de cada polígono original dentro de los diferentes mosaicos del mosaico. También sería relevante calcular el número de solapamientos.

Estoy buscando una solución utilizando PostGIS. He resuelto esto usando FME, pero el tamaño de mis datos (más de 20.000 objetos con millones de solapamientos) hace que esto consuma demasiado tiempo. El función de unión en QGIS lo hace de forma bastante sencilla, tal y como ilustran sus páginas web de documentos.

Para los que necesiten un contexto, los datos son en realidad topes de varios anillos alrededor de los sitios de anidación de aves marinas en Noruega. Los conjuntos de datos son por especie y estamos viendo varios miles de sitios de anidación por especie.

Puede descargar un conjunto de datos de prueba aquí .

Answer 1

Descargue los complementos de PostGIS desde este enlace: https://github.com/pedrogit/postgisaddons

Instale ejecutando el archivo postgis_addons.sql.

Pruebe ejecutando el archivo postgis_addons_test.sql.

Aquí hay un ejemplo de un problema similar al suyo:

WITH geomtable AS (
  SELECT 1 id, ST_GeomFromText('POLYGON((0 0, 0 2, 2 2, 2 0, 0 0), (0.2 0.5, 0.2 1.5, 0.8 1.5, 0.8 0.5, 0.2 0.5))') geom
  UNION ALL
  SELECT 2 id, ST_GeomFromText('POLYGON((1 0.2, 1 1, 3 1, 3 0.2, 1 0.2))') geom
  UNION ALL
  SELECT 3 id, ST_GeomFromText('POLYGON((1.5 0.8, 1.5 1.2, 2.5 1.2, 2.5 0.8, 1.5 0.8))') geom
  UNION ALL
  SELECT 4 id, ST_GeomFromText('MULTIPOLYGON(((3 0, 3 2, 5 2, 5 0, 3 0)), ((4 3, 4 4, 5 4, 5 3, 4 3)))') geom
)
SELECT a.id, unnest(ST_SplitAgg(a.geom, b.geom, 0.00001)) geom
FROM geomtable a,
     geomtable b
WHERE ST_Equals(a.geom, b.geom) OR
      ST_Contains(a.geom, b.geom) OR
      ST_Contains(b.geom, a.geom) OR
      ST_Overlaps(a.geom, b.geom)
GROUP BY a.id;

En su consulta, sólo tiene que sustituir geomtable por su tabla:

SELECT a.id, unnest(ST_SplitAgg(a.geom, b.geom, 0.00001)) geom
FROM polytest a,
     polytest b
WHERE ST_Equals(a.geom, b.geom) OR
      ST_Contains(a.geom, b.geom) OR
      ST_Contains(b.geom, a.geom) OR
      ST_Overlaps(a.geom, b.geom)
GROUP BY a.id;

Debería funcionar para miles de polígonos y cuando hay más de dos solapamientos.

Answer 2

Puedes ir por partes... Utilice el enfoque mencionado anteriormente:

CREATE TABLE tbl_intersection AS
SELECT
    ST_INTERSECTION(t1.geom, t2.geom) as geom
FROM original_dataset t1, original_dataset t2 
WHERE t1.ID < t2.ID AND ST_INTERSECTS(t1.geom, t2.geom);

Y luego se utiliza la tabla anterior y el conjunto de datos original para ejecutar el ST_Diferencia (El orden es importante en este caso) entre entonces. Debería ponerte en la pista

SELECT
    ST_Difference(t2.geom, t1.geom) AS geom
FROM
  original_dataset t1, tbl_intersection t2

Answer 3

Si has añadido todos los índices correctos (en el ID así como el índice GIST en geom), lo siguiente debería hacer el trabajo - siempre y cuando tu servidor esté a la altura;

   select st_intersection(a.geom, b.geom) as geom
   from table a, table b
   where st_intersects(a.geom, b.geom)
   group by st_intersection(a.geom, b.geom)