纬度/经度距离通过转换到UTM计算给出不同的结果比使用近似方法

我有两种方法用于计算在python地理参考坐标之间的距离：纬度/经度距离通过转换到UTM计算给出不同的结果比使用近似方法

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

我与以下值拨打：

lat1 = 34.866527 
lon1 = 69.674606 
lat2 = 34.864990 
lon2 = 69.657655 
print "approximation method: ", calc_distance(lat1, lon1, lat2, lon2) 
print "converting to utm method: ", calc_distance_convert_utm(lat1, lon1, lat2, lon2)

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

approximation method: 1.55593476881 
converting to utm method: 1928.21537269

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

我做错了什么？任何帮助表示赞赏！谢谢

来源

2017-07-07 M F

我发现错误... 在utm转换方法中，我已经在转换过程中切换了纬度/经度值。它应该是：

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

来源

2017-07-07 09:35:32

纬度/经度距离通过转换到UTM计算给出不同的结果比使用近似方法

回答

相关问题