Quão caro é ST_GeomFromText

Question

No postgis, é o ST_GeomFromText ligar muito caro?Eu pergunto principalmente porque tenho uma consulta chamada com frequência que tenta encontrar o ponto mais próximo de outro ponto que corresponda a alguns critérios e que também esteja a uma certa distância desse outro ponto, e da maneira como escrevi atualmente, está fazendo o mesmo ST_GeomFromText duas vezes:

 $findNearIDMatchStmt = $postconn->prepare(
    "SELECT     internalid " .
    "FROM       waypoint " .
    "WHERE      id = ? AND " .
    "           category = ? AND ".
    "           (b.category in (1, 3) OR type like ?) AND ".
    "           ST_DWithin(point, ST_GeomFromText(?," . SRID .
    "           ),".  SMALL_EPSILON . ") " .
    "           ORDER BY ST_Distance(point, ST_GeomFromText(?,", SRID .
    "           )) " .
    "           LIMIT 1");

Existe uma maneira melhor de reescrever isso?

Um pouco AT:Na tela de visualização, todos os meus sublinhados estão sendo renderizados como & # 9 5 ; - Espero que isso não apareça assim no post.

Answer 1

eu não acredito ST_GeomFromText() é particularmente caro, embora no passado eu tenha otimizado PostGIS consultas criando uma função, declarando uma variável e então atribuindo o resultado de ST_GeomFromText para a variável.

Você já tentou verificar o plano de execução da sua consulta com uma variedade de parâmetros diferentes, porque isso deve lhe dar uma ideia definitiva de quais partes da consulta estão demorando?

Acredito que a maior parte do tempo de execução estará nas chamadas para ST_DWithin() e ST_Distance(), embora se as colunas id e categoria não estiverem indexadas, poderá estar fazendo uma varredura de tabela interessante.

Answer 2

@Ubiguch parece que ST_DWithin usa o índice espacial, o que parece reduzir rapidamente o número de pontos a serem consultados.

 navaid=> explain select internalid from waypoint where id != 'KROC' AND ST_DWithin(point,                                                                  ST_GeomFromText('POINT(-77.6723888888889 43.1188611111111)',4326), 0.05) order by st_distance(point, st_geomfromtext('POINT(-77.6723888888889 43.1188611111111)',4326)) limit 1;
                                                                                                                                                                                                                                                      QUERY PLAN                                                                                                                                                                                                                                                       
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 Limit  (cost=8.37..8.38 rows=1 width=104)
   ->  Sort  (cost=8.37..8.38 rows=1 width=104)
         Sort Key: (st_distance(point, '0101000020E61000002FFE676B086B53C0847E44D7368F4540'::geometry))
         ->  Index Scan using waypoint_point_idx on waypoint  (cost=0.00..8.36 rows=1 width=104)
               Index Cond: (point && '0103000020E61000000100000005000000000000C03B6E53C000000060D0884540000000C03B6E53C0000000409D95454000000020D56753C0000000409D95454000000020D56753C000000060D0884540000000C03B6E53C000000060D0884540'::geometry)
               Filter: (((id)::text <> 'KROC'::text) AND (point && '0103000020E61000000100000005000000000000C03B6E53C000000060D0884540000000C03B6E53C0000000409D95454000000020D56753C0000000409D95454000000020D56753C000000060D0884540000000C03B6E53C000000060D0884540'::geometry) AND ('0101000020E61000002FFE676B086B53C0847E44D7368F4540'::geometry && st_expand(point, 0.05::double precision)) AND (st_distance(point, '0101000020E61000002FFE676B086B53C0847E44D7368F4540'::geometry) < 0.05::double precision))
(6 rows)

Sem o order by e a limit, parece que uma consulta típica está retornando apenas de 5 a 10 pontos de referência no máximo.Portanto, provavelmente não devo me preocupar com o custo adicional do filtro aplicado aos pontos retornados.