PostGIS Intersección comportamiento extraño

Question

Observo un comportamiento extraño con mi consulta PostGIS. Estoy tratando de lograr lo siguiente:

• Tengo una capa poligonal
• Tengo 2 capas de línea diferentes
  • capa con líneas horizontales
  • capa con líneas de 45
• el objetivo es cortar las líneas en los polígonos para obtener las partes de línea de cada línea por polígono

He aquí un ejemplo:

Estoy utilizando una consulta que se parece a esto:

CREATE TABLE split AS
    ( SELECT row_number() over() AS id,
             a.id AS old_id,
             st_intersection(a.geom, b.geom) AS geom
     FROM lines a,
          polygons b
     WHERE st_intersects(a.geom, b.geom) ) ;

La consulta tarda unos 3 minutos para las líneas horizontales y verticales, pero más de 5 horas (canceladas ahí) para las líneas de 45°.

¿A alguien se le ocurre una posible causa? Todas las capas tienen índices espaciales, pk y demás. Estoy utilizando un Postgres 9.5 con PostGIS 2.2.

Realmente no tengo ni idea de ese comportamiento.

Answer 1

He reproducido tu problema y he hecho algunas pruebas. He creado una intersección entre una capa de polígonos y con líneas horizontales (cada km), luego con líneas a 45°, y luego con las mismas líneas a 45° divididas cada km verticalmente (lo que equivale a sqrt(2) = 1.414 km).

Estos son los resultados de la consulta de intersección:

• horizontal : 8,7 seg (sin índice)
• horizontal : 8,9 seg (índice)
• 45° : 60 seg (sin índice)
• 45° : 61 seg (índice)
• 45° split : 120 seg (sin índice)
• 45° split : 32,2 seg (índice)

Así que parece que el índice espacial no tiene mucho impacto con líneas largas. Resulta más útil con líneas más cortas.