Encontrar puntos dentro de un área elíptica utilizando PostGIS?

Question

Usando PostGIS 2.0.0 necesito filtrar los resultados de una función, que devuelve muchos PUNTOS, para encontrar los que se encuentran dentro de un área elíptica (como una línea de gran círculo con un buffer alrededor) definida por dos puntos y una distancia (todos los parámetros pasados en la consulta).

Un filtro usando ST_Distance hace el truco, pero me pregunto si hay alguna manera de acelerarlo.

Consulta simplificada:

SELECT candidate_point
FROM my_function(some, params)
WHERE
    ST_Distance(candidate_point, 'SRID=4326;POINT(1 2)'::geography) +
    ST_Distance(candidate_point, 'SRID=4326;POINT(3 4)'::geography) <= 60000

Intenté usar ST_DWithin como filtro adicional, pero eso lo hace diez veces más lento. Dado que candidate_point proviene de una función, no de una tabla, no hay ningún índice que utilizar.

Answer 1

Si sus datos están en una proyección de área equidistante, podría hacer esto:

SELECT gid,the_geog FROM mytable WHERE ST_DWithIn(the_geog,ST_MakeLine(point1,point2),corridor_width_m);

ST_DWithIn también funciona con los tipos de geografía, pero incluso si usted lanza a la geometría todavía será lento porque usted no tiene índices de geometría.

Si quieres velocidad, tienes que usar tipos de geometría.

Answer 2

No he probado esto pero podría ayudar algo así:

SELECT candidate_point
FROM my_function(some, params)
WHERE
    ST_Distance(candidate_point, 'SRID=4326;POINT(1 2)'::geography) <= 60000 &&
    ST_Distance(candidate_point, 'SRID=4326;POINT(3 4)'::geography) <= 60000 &&
    ST_Distance(candidate_point, 'SRID=4326;POINT(1 2)'::geography) +
    ST_Distance(candidate_point, 'SRID=4326;POINT(3 4)'::geography) <= 60000

Alternativamente, si pudieras calcular el cuadro delimitador de la elipse antes y usar eso para filtrar los puntos antes de hacer las llamadas a ST_Distance, esperaría un aumento de velocidad.

Answer 3

Prueba a eliminar ::geografía. Eso podría ayudar