8

ユーザーから2つの地理座標を取得し、それらの間の距離を計算する小さなプログラムを設定しています（地球の曲率を考慮に入れて）。だから私は式がここにあるものについてウィキペディアを調べました。

私は基本的にそれに基づいてPython関数を設定し、これが私が思いついたものです：

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

問題は、結果が本当に不正確になることです。私はPythonを初めて使用するので、ヘルプやアドバイスをいただければ幸いです。

4

4 に答える 4

12

あなたは4つまたは5つまたは6つの問題を抱えています：

（1 end_lat = math.radians(end_long)）end_lat = math.radians(end_lat)

（2）誰かがすでに述べたように、あなたはいくつかのものを見逃しています。

（3）コードが判読できない（行が長すぎる、括弧が冗長である、「数学」の17の無意味なインスタンス）。

（4）ウィキペディアの記事での使用に関するコメントに気づかなかったatan2()

（5）座標を入力するときに緯度と経度を入れ替えていた可能性があります

（6）delta(latitude)不必要に計算されます。数式には表示されません

すべてを一緒に入れて：

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
于 2012-01-14T02:25:06.543 に答える
4

距離計算用の組み込み関数を備えたgeopyモジュールを使用できます。下のリンクの「CalculatingDistances」までスクロールダウンしてください： https ：//pypi.python.org/pypi/geopy

于 2012-01-13T23:59:55.633 に答える
4

これは機能します（printfは実際の例@http: //en.wikipedia.org/wiki/Great-circle_distanceに従って2887.26kmを返します）：

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

提出物の次の行に注意してください。end_lat = math.radians(end_long)

于 2012-01-14T01:04:31.623 に答える
3

私はあなたが最初にmath.sin（d_long）を逃したと思います、多分これであるはずです：

 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))) )
于 2012-01-14T00:01:26.693 に答える