Início / dba / Perguntas / 156167
Accepted
user606521
Asked: 2016-11-24 06:18:26 +0800 CST

Por que o operador earthdistance <@> retorna um valor diferente da função earth_distance?

  • 772

Vendo um valor diferente com <@>eearth_distance

with points as (
  select 
    '(1,1)'::point as p1
    ,'(2,2)'::point as p2
)
select 
  ((p1 <@> p2) * 1609.34::double precision) as d1
  ,(earth_distance(ll_to_earth(p1[1], p1[0]), ll_to_earth(p2[1], p2[0]))) as d2
from points

Resultado é:

157224.717583288 157402.326073052

Quase 200m de diferença - por que isso?

postgresql spatial

1 respostas

  1. Best Answer

    Eles fazem coisas diferentes,

    • earth_distanceRetorna a distância do grande círculo entre dois pontos na superfície da Terra.
    • <@>Dá a distância em milhas terrestres entre dois pontos na superfície da Terra.

    Portanto, esta é a diferença entre a média dos dois (supondo que nenhum deles esteja correto) é dado como,

    SELECT ((d2-d1) / ((d2+d1)/2)) * 100 AS distance_diff FROM results;
   distance_diff   
-------------------
 0.112652419817465
(1 row)

    Que é 0,11%. Pela distância do Grande Círculo na wikipedia

    A Terra é quase esférica (veja o raio da Terra), então as fórmulas de distância do grande círculo dão a distância entre os pontos na superfície da Terra (como o vôo do corvo) correta dentro de 0,5% ou mais.

    Então, qual é mais preciso? <@>, claramente. Mas, podemos fazer melhor, como dizem os documentos,

    Neste módulo, a Terra é considerada perfeitamente esférica. (Se isso for muito impreciso para você, você pode querer olhar para o projeto PostGIS.)

    SELECT ST_Distance(
  ST_MakePoint(1,1)::geography,
  ST_MakePoint(2,2)::geography
);
   st_distance   
-----------------
 156876.14940189
(1 row)

    Mas se seus coordenadores são SRID 4326, a distância mais correta é 156876.14940189.

    • 4

relate perguntas