转换为utm的Lat/lon距离计算结果与近似方法不同

Question

我有两种方法来计算python中地理参考坐标之间的距离：

from pyproj import Proj
import math

def calc_distance(lat1, lon1, lat2, lon2):
    """
    Calculate the great circle distance between two points
    on the earth (specified in decimal degrees)
    """
    # convert decimal degrees to radians
    lon1, lat1, lon2, lat2 = map(math.radians, [lon1, lat1, lon2, lat2])

    # haversine formula
    dlon = lon2 - lon1
    dlat = lat2 - lat1
    a = math.sin(dlat / 2) ** 2 + math.cos(lat1) * math.cos(lat2) * math.sin(dlon / 2) ** 2
    c = 2 * math.asin(math.sqrt(a))
    km = 6371 * c

    return km

def calc_distance_convert_utm(lat1, lon1, lat2, lon2):
    myProj = Proj("+proj=utm +zone=42, +north +ellps=WGS84 +datum=WGS84 +units=m +no_defs")

    # convert to utm 
    utm_x1, utm_y1 = myProj(lat1, lon1)
    utm_x2, utm_y2 = myProj(lat2, lon2)

    diff_x = abs(utm_x1 - utm_x2)
    diff_y = abs(utm_y1 - utm_y2)

    distance = math.sqrt(diff_x**2 + diff_y**2)

    return distance

我称之为以下值：

^{pr2}$

但是，如果比较结果，会得到两个不同的值：

approximation method:  1.55593476881
converting to utm method:  1928.21537269

注意，第一种方法以公里为单位返回距离，而第二种方法以米为单位返回距离。 我把结果与你可以在网上找到的距离计算器进行了比较，似乎第一种方法（近似法）是“更正确”的答案，因为这是大多数在线计算器返回的值。我想知道，为什么第二个方法（首先转换为utm）没有返回更相似的结果（比如1555.9347…）。我有将近0.5公里的差距，这在我看来差不多。在

我做错什么了吗？ 感谢任何帮助！谢谢

Answer 1

我发现了错误。。。 在utm转换方法中，我切换了转换过程中的lat/lon值。它应该是：

utm_x1, utm_y1 = myProj(lon1, lat1)
utm_x2, utm_y2 = myProj(lon2, lat2)