Python 两地之间距离计算

两地之间距离计算

两地之间距离计算（模块：numEx，所在文件名 num_hw.py，Level：★）
利用 Python 实现地球上两点之间的距离计算，地球上点的位置以经纬度坐标形式提供。
距离计算采用 Haversine 公式：

这里 r 是地球半径 6371Km， 代表点的（纬度，经度）坐标。
参考网站：https://en.wikipedia.org/wiki/Haversine_formula
完成距离计算函数：
函数原型：def sphere_distance(p1, p2)
参数 p1：tuple 元组类型，二元组，（纬度，经度），坐标精确到小数点后 7 位
参数 p2：tuple 元组类型，二元组，（纬度，经度），坐标精确到小数点后 7 位
纬度取值范围：[0-90]，经度取值范围：[0-180]，单位均为角度；而 Haversine 公式计算
采用的是弧度，注意转换。
返回值：如果输入的坐标数据合规，则返回两点之间的距离，单位为 Km，保留两位小
数；如果输入的坐标不合规，返回错误“Parameter Error.”

代码如下（函数引用）

# 题目：两地之间距离计算

# 计算两点p1, p2之间的距离
# p1：（纬度、经度）
# p2: （纬度、经度）
def sphere_distance(p1, p2):
    import math
    if p1[0] < 0 or p2[0] < 0 or p1[0] > 90 or p2[0] > 90 or p1[1] < 0 or p2[1] < 0 or p1[1] > 180 or p2[1] > 180:
        return 'Parameter Error.'
    r = 6371
    i1 = (p1[0] / 180) * math.pi
    j1 = (p1[1] / 180) * math.pi
    i2 = (p2[0] / 180) * math.pi
    j2 = (p2[1] / 180) * math.pi
    d = 2 * r * math.asin(math.sqrt((math.sin((i2 - i1) / 2)) ** 2 + math.cos(i1) * math.cos(i2) * (math.sin((j2 - j1) / 2)) ** 2))
    return d