Как рассчитать расстояние от широты и долготы

Question

В базе данных Postgres 9.1+ содержит точки маршрута:

create table doksyndm (
id serial primary key, 
readingtime timestamp with time zone not null, 
latitude float, 
longitude float 
); 

Очки сохраняются в браузере с помощью navigator.geolocation вызовы из JavaScript.

Как рассчитать длину маршрута с этих данных? Маршрут составляет менее 300 км, поэтому можно суммировать прямые линии между точками. Он должен работать в Postgres 9.1 и новее.

Результат будет одним номером, общим расстоянием между точками в километрах. Точки отправления кода базы данных и javascript могут быть пересчитаны, если это разумно.

Ответы в Calculate distance between two points in google maps V3 как реализовать это на продуральных языках.

Они используют цикл, который недоступен в SQL. Maby SUM() с функцией окна может использоваться для реализации этого?

Answer 1

Я предполагаю, что вы хотите, расстояние от первой точки до следующего самого раннего, и так далее ...

По Postgresql 9.3, есть элегантное решение с помощью lateral join.

Во-первых, установите расширение PostGIS, так что вы можете использовать свои функции расстояния:

CREATE EXTENSION POSTGIS; 

Теперь запрос:

SELECT SUM(st_distance(st_point(a.longitude,a.latitude)::geography, 
         st_point(b.longitude,b.latitude)::geography) 
      )/1000 AS distance_km FROM doksyndm a 
LEFT JOIN LATERAL 
(SELECT * FROM doksyndm WHERE readingtime > a.readingtime limit 1) b ON TRUE; 

Этот запрос с использованием distinct on должны работать для всех версий:

SELECT SUM (st_distance(st_point(alon,alat)::geography, 
       st_point(blon,blat)::geography))/1000 as distance_km FROM 
    (select distinct ON (a.readingtime) a.readingtime, 
             a.latitude as alat, 
             a.longitude as alon, 
             b.latitude as blat , 
             b.longitude as blon 
    FROM doksyndm a inner join doksyndm b on a.readingtime < b.readingtime 
ORDER by a.readingtime,b.readingtime) x