PostgreSQL Breite Länge Abfrage

Question

ich habe latitude und longitude Spalten in location Tabelle in PostgreSQL-Datenbank, und ich versuche, Distanzabfrage mit einer PostgreSQL-Funktion auszuführen.

ich lesen Sie dieses Kapitel des Handbuchs:

https://www.postgresql.org/docs/current/static/earthdistance.html

aber ich denke, ich dort etwas fehlt bin.

Wie soll ich das tun? Sind weitere Beispiele verfügbar?

Answer 1

Dieses Modul ist optional und wird nicht in der standardmäßigen PostgreSQL Instalatlion installiert. Sie müssen es aus dem Verzeichnis contrib installieren.

Sie können die folgende Funktion verwenden, um die ungefähre Entfernung zwischen den Koordinaten (in Meilen) zu berechnen:

CREATE OR REPLACE FUNCTION distance(lat1 FLOAT, lon1 FLOAT, lat2 FLOAT, lon2 FLOAT) RETURNS FLOAT AS $$ 
DECLARE             
    x float = 69.1 * (lat2 - lat1);       
    y float = 69.1 * (lon2 - lon1) * cos(lat1/57.3);   
BEGIN              
    RETURN sqrt(x * x + y * y);        
END 
$$ LANGUAGE plpgsql; 

Answer 2

Vorausgesetzt, dass Sie das earthdistance Modul korrekt installiert haben, dies wird Ihnen die Entfernung in Meilen zwischen zwei Städten . Diese Methode verwendet die einfacheren punktbasierten Erdabstände. Beachten Sie, dass die Argumente für Punkt() zuerst Längengrad und dann Breitengrad sind.

create table lat_lon (
    city varchar(50) primary key, 
    lat float8 not null, 
    lon float8 not null 
); 

insert into lat_lon values 
('London, GB', 51.67234320, 0.14787970), 
('New York, NY', 40.91524130, -73.7002720); 

select 
    (
    (select point(lon,lat) from lat_lon where city = 'London, GB') <@> 
    (select point(lon,lat) from lat_lon where city = 'New York, NY') 
) as distance_miles 

distance_miles 
-- 
3447.58672105301 

Answer 3

Hier ist ein weiteres Beispiel der Punkt-Operator:

create extension cube; 
create extension earthdistance; 
select (point(-0.1277,51.5073) <@> point(-74.006,40.7144)) as distance; 

    distance  
------------------ 
3461.10547602474 
(1 row) 

Beachten Sie, dass points mit LONGITUDE FIRST erstellt werden. Per der documentation:

Punkte werden als (Länge, Breite) aufgenommen und nicht umgekehrt, weil Länge näher an die intuitive Idee der x-Achse und Breitengrad y-Achse ist.

Welches ist schreckliches Design ... aber so ist es.

Answer 4

Eine genauere Version von @ strkol Antwort, wobei den Eingang Haversine formula

CREATE OR REPLACE FUNCTION distance(
    lat1 double precision, 
    lon1 double precision, 
    lat2 double precision, 
    lon2 double precision) 
    RETURNS double precision AS 
$BODY$ 
DECLARE 
    R integer = 6371e3; -- Meters 
    rad double precision = 0.01745329252; 

    φ1 double precision = lat1 * rad; 
    φ2 double precision = lat2 * rad; 
    Δφ double precision = (lat2-lat1) * rad; 
    Δλ double precision = (lon2-lon1) * rad; 

    a double precision = sin(Δφ/2) * sin(Δφ/2) + cos(φ1) * cos(φ2) * sin(Δλ/2) * sin(Δλ/2); 
    c double precision = 2 * atan2(sqrt(a), sqrt(1-a));  
BEGIN              
    RETURN R * c;   
END 
$BODY$ 
    LANGUAGE plpgsql VOLATILE 
    COST 100; 

verwendet, ist in Graden (z.B. 52,34273489, 6,23847) und Ausgang ist in Metern.