Cần trợ giúp tính khoảng cách địa lý

Question

Tôi đang thiết lập một chương trình nhỏ để lấy 2 toạ độ địa lý từ người dùng và sau đó tính khoảng cách giữa chúng (tính đến độ cong của trái đất). Vì vậy, tôi tra cứu wikipedia về công thức là here.

tôi về cơ bản thiết lập chức năng python của tôi dựa trên đó và đây là những gì tôi đã đưa ra:

def geocalc(start_lat, start_long, end_lat, end_long): 
    start_lat = math.radians(start_lat) 
    start_long = math.radians(start_long) 
    end_lat = math.radians(end_long) 
    end_long = math.radians(end_long) 

    d_lat = start_lat - end_lat 
    d_long = start_long - end_long 

    EARTH_R = 6372.8 

    c = math.atan((math.sqrt((math.cos(end_lat)*d_long)**2 +((math.cos(start_lat)*math.sin(end_lat)) - (math.sin(start_lat)*math.cos(end_lat)*math.cos(d_long)))**2))/((math.sin(start_lat)*math.sin(end_lat)) + (math.cos(start_lat)*math.cos(end_lat)*math.cos(d_long)))) 

    return EARTH_R*c 

Vấn đề là kết quả đi ra thực sự không chính xác. Tôi mới để python nên một số trợ giúp hoặc tư vấn sẽ được đánh giá rất nhiều!

Answer 1

Bạn đã có 4 hoặc 5 hoặc 6 vấn đề:

(1) end_lat = math.radians(end_long) nên end_lat = math.radians(end_lat)

(2) bạn đang thiếu một số nội dung như ai đó đã đề cập, có lẽ hầu hết vì

(3) mã của bạn là không đọc được (dòng quá dài, ngoặc thừa, 17 trường hợp vô nghĩa của "toán học").

(4), bạn không để ý nhận xét trong bài viết Wikipedia về việc sử dụng atan2()

(5) Bạn có thể đã được trao đổi lat và lon khi nhập tọa độ của bạn

(6) delta(latitude) được tính cách không cần thiết; nó không xuất hiện trong công thức

Đưa nó tất cả cùng nhau:

from math import radians, sqrt, sin, cos, atan2 

def geocalc(lat1, lon1, lat2, lon2): 
    lat1 = radians(lat1) 
    lon1 = radians(lon1) 
    lat2 = radians(lat2) 
    lon2 = radians(lon2) 

    dlon = lon1 - lon2 

    EARTH_R = 6372.8 

    y = sqrt(
     (cos(lat2) * sin(dlon)) ** 2 
     + (cos(lat1) * sin(lat2) - sin(lat1) * cos(lat2) * cos(dlon)) ** 2 
     ) 
    x = sin(lat1) * sin(lat2) + cos(lat1) * cos(lat2) * cos(dlon) 
    c = atan2(y, x) 
    return EARTH_R * c 



>>> geocalc(36.12, -86.67, 33.94, -118.40) 
2887.2599506071115 
>>> geocalc(-6.508, 55.071, -8.886, 51.622) 
463.09798886300376 
>>> geocalc(55.071, -6.508, 51.622, -8.886) 
414.7830891822618 

Answer 2

Bạn có thể sử dụng các mô-đun geopy trong đó có một chức năng built-in để tính toán khoảng cách, di chuyển xuống "Khoảng cách tính toán" trong liên kết dưới đây: https://pypi.python.org/pypi/geopy

Answer 3

Tôi nghĩ rằng bạn bỏ lỡ một Math.sin (d_long) về phía đầu, có lẽ phải được điều này:

c = math.atan((math.sqrt((math.cos(end_lat)*math.sin(d_long))**2 +((math.cos(start_lat)*math.sin(end_lat)) - (math.sin(start_lat)*math.cos(end_lat)*math.cos(d_long)))**2))/((math.sin(start_lat)*math.sin(end_lat)) + (math.cos(start_lat)*math.cos(end_lat)*math.cos(d_long)))) 

Answer 4

này hoạt động (in f trả 2.887,26 km theo ví dụ cách @http://en.wikipedia.org/wiki/Great-circle_distance):

import math 

def geocalc(start_lat, start_long, end_lat, end_long): 

    start_lat = math.radians(start_lat) 
    start_long = math.radians(start_long) 
    end_lat = math.radians(end_lat) 
    end_long = math.radians(end_long) 

    d_lat = math.fabs(start_lat - end_lat) 
    d_long = math.fabs(start_long - end_long) 

    EARTH_R = 6372.8 

    y = ((math.sin(start_lat)*math.sin(end_lat)) + (math.cos(start_lat)*math.cos(end_lat)*math.cos(d_long))) 

    x = math.sqrt((math.cos(end_lat)*math.sin(d_long))**2 + ((math.cos(start_lat)*math.sin(end_lat)) - (math.sin(start_lat)*math.cos(end_lat)*math.cos(d_long)))**2) 

    c = math.atan(x/y) 

    return EARTH_R*c 

f = geocalc(36.12, -86.67, 33.94, -118.40) 
print f 

Thông báo dòng này trong trình của bạn: end_lat = math.radians(end_long)