Truy vấn kinh độ vĩ độ PostgreSQL

Question

tôi có latitude và longitude cột trong location bảng trong cơ sở dữ liệu PostgreSQL, và tôi đang cố thực hiện truy vấn khoảng cách bằng hàm PostgreSQL.

tôi đọc chương này của cuốn cẩm nang:

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

nhưng tôi nghĩ rằng tôi thiếu cái gì đó.

Tôi nên làm như thế nào? Có các ví dụ khác có sẵn

Answer 1

Mô-đun này là tùy chọn và không được cài đặt trong mặc định PostgreSQL instalatlion. Bạn phải cài đặt nó từ thư mục contrib.

Bạn có thể sử dụng chức năng sau đây để tính toán khoảng cách tương đối giữa các tọa độ (trong dặm):

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

Giả sử bạn đã cài đặt các module earthdistance một cách chính xác, điều này sẽ cung cấp cho bạn khoảng cách trong dặm giữa hai thành phố . Phương pháp này sử dụng khoảng cách trái đất dựa trên điểm đơn giản hơn. Lưu ý rằng đối số cho điểm() là kinh độ đầu tiên, sau đó là vĩ độ.

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

Dưới đây là một ví dụ khác sử dụng toán tử điểm:

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

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

Lưu ý rằng points được tạo ra với LONGITUDE FIRST. Theo documentation:

điểm được thực hiện như (kinh độ, vĩ độ) và không phải ngược lại vì kinh độ là gần gũi hơn với ý tưởng trực giác của trục x và vĩ độ để trục y.

Thiết kế khủng khiếp ... nhưng đó là cách thực hiện.

Answer 4

Một phiên bản chính xác hơn về câu trả lời @ strkol, sử dụng giao 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; 

Input là trong độ (ví dụ 52,34273489, 6,23847) và đầu ra là tính bằng mét.