なぜ私のPython haversine距離計算は、オンラインツールやGoogleマップに比べて間違っていますか？

Question

小さな自律RCカープロジェクトの一環として、haversine distanceとangle calculatorをPythonで書いています。私の2つの試験場所は38.63594444444444,-90.2315と38.63594444444444,-90.23211111111111です。

ほとんどのオンライン電卓（および私の個人的なTI-89）の距離は約0.05308kmです。しかし、私のPython関数は0.06795 kmの距離を得ています。それは〜15メートルオフであり、それは小さなRCカーで巨大です。

私の軸受け計算機能points2angleは、toDegrees機能で浮動鋳造を行うまで失敗していました。整数部が私を苦しめていた。

私のpoints2angleとpoints2distance関数は（度、分、秒）のタプルが必要であることに注意してください。 2つのテスト場所は、その形式の(38, 38, 9.4), (-90, 13, 53.4)と(38, 38, 9.4), (-90, 13, 55.6)です。

編集：ありがとうございますMSeifert。私はちょうど私の緯度と経度が混ざった。私はpoints2angleコードを以下のように修正しましたが、私のpoints2distanceコードにエラーを残しました。私の間違ったコードとその答えの違いはまだ明らかです。

マイ距離計算（間違った距離が返さ）：

def points2distance(start, end): 
     start_long = math.radians(toDegrees(start[0])) 
     start_latt = math.radians(toDegrees(start[1])) 
     end_long = math.radians(toDegrees(end[0])) 
     end_latt = math.radians(toDegrees(end[1])) 
     d_latt = float(end_latt - start_latt) 
     d_long = float(end_long - start_long) 
     a = (math.sin(d_latt/2)**2) + math.cos(start_latt) * math.cos(end_latt)* (math.sin(d_long/2)**2) 
     c = 2 * math.atan2(math.sqrt(a), math.sqrt(1-a)) 
     return 6371 * c 

マイ小数度変換機能（作動）：

def toDegrees(coord): 
    degrees = float(math.fabs(coord[0])) + float(coord[1]/60) + float(coord[2]/3600) 
    if coord[0] < 0: 
     degrees = degrees*-1 
    return degrees 

マイ方位角計算（作業）：

def points2angle(start, end): 
    start_long = math.radians(toDegrees(start[1])) 
    start_latt = math.radians(toDegrees(start[0])) 
    end_long = math.radians(toDegrees(end[1])) 
    end_latt = math.radians(toDegrees(end[0])) 
    d_latt = end_latt - start_latt 
    d_long = end_long - start_long 
    y = math.sin(d_long)*math.sin(end_latt) 
    x = (math.cos(start_latt)*math.sin(end_latt)) - (math.sin(start_latt)*math.cos(end_latt)*math.cos(d_long)) 
    brng = math.degrees(math.atan2(y,x)) 
    compBear = (brng+360) % 360; 
    return compBear 

Answer 1

あなたの経度と緯度が混在しているだけです（あるいは、それらがスワップされた場合はswそれらをそこにap）それは働く：

def points2distance(start, end): 
    start_long = math.radians(toDegrees(start[1])) 
    start_latt = math.radians(toDegrees(start[0])) 
    end_long = math.radians(toDegrees(end[1])) 
    end_latt = math.radians(toDegrees(end[0])) 
    d_latt = float(end_latt - start_latt) 
    d_long = float(end_long - start_long) 
    a = (math.sin(d_latt/2)**2) + math.cos(start_latt) * math.cos(end_latt)* (math.sin(d_long/2)**2) 
    c = 2 * math.atan2(math.sqrt(a), math.sqrt(1-a)) 
    return 6371 * c 

points2distance([(38, 38, 9.4), (-90, 13, 53.4)], [(38, 38, 9.4), (-90, 13, 55.6)]) 
# returns: 0.053079628495340196