Eliminación de geometrías no intersecantes en PostGIS

Answer

1. 

   davy.ai

   July 23, 2023 at 11:37 am

   Una posible solución a este problema es utilizar la función ST_Intersects en PostGIS. Esta función se puede utilizar para verificar si una geometría se cruza con otra geometría. En este caso, podemos usarla para verificar si cada mosaico se cruza con el polígono que contiene los datos.

   Aquí hay una forma de implementar esta solución:

   1. Primero, necesitamos crear un polígono que cubra el área donde tenemos los datos. Este polígono debe tener la misma proyección que los mosaicos. Llamemos a este polígono “data_polygon”.

   2. A continuación, podemos utilizar la función ST_Intersects para verificar si cada mosaico se cruza con data_polygon. Así es como podemos modificar la consulta original para incluir esta verificación:

   SELECT (ST_Dump(ST_Intersection(grid.geom, data_polygon.geom))).geom as geom
   FROM (
       SELECT ST_SquareGrid(0.2, geom) AS geom
       FROM some_table
       WHERE some_condition
   ) AS grid
   JOIN data_polygon ON ST_Intersects(grid.geom, data_polygon.geom);
   

   Lo que estamos haciendo aquí es unir la cuadrícula de mosaicos con data_polygon en la condición de que se crucen. Esto solo devolverá los mosaicos que cubren el área con datos.

   1. Por último, podemos agregar un índice en data_polygon para acelerar la verificación de intersección. Podemos crear un índice espacial utilizando el siguiente comando:

   CREATE INDEX data_polygon_idx ON data_polygon USING gist (geom);
   

   Esto creará un índice que debería hacer que la verificación de la intersección sea mucho más rápida.

   Con estas modificaciones, la consulta solo debería devolver mosaicos que cubren el área con datos y no crear mosaicos vacíos.

Comments are closed.