spherical distance

sherlockjjobs

已于 2023-05-05 11:47:10 修改

阅读量88

点赞数 1

分类专栏： GIS Algorithms 文章标签： python 开发语言

于 2023-05-05 11:44:34 首次发布

本文链接：https://blog.csdn.net/weixin_43456810/article/details/130503804

版权

GIS Algorithms 专栏收录该内容

8 篇文章 0 订阅

订阅专栏

spherical distance

请添加图片描述

import math


def spherical_distance(lat1, lon1, lat2, lon2):
    """
    Calculates the great circle distance given
    the latitudes and longitudes of two points.
    Input
     lat1, lon1: lat and long in degrees for the first point
     lat2, lon2: lat and long in degrees for the second point
    Output
     d: great circle distance
    """
    D = 3959   # earth radius in miles
    phi1 = math.radians(lat1)
    lambda1 = math.radians(lon1)
    phi2 = math.radians(lat2)
    lambda2 = math.radians(lon2)
    dlambda = lambda2 - lambda1

    dphi = phi2 - phi1
    sinlat = math.sin(dphi/2.0)
    sinlong = math.sin(dlambda/2.0)

    alpha = (sinlat*sinlat) + math.cos(phi1) * \
        math.cos(phi2) * (sinlong*sinlong)
    c = 2 * math.asin(min(1, math.sqrt(alpha)))
    d = D*c
    return d


if __name__ == "__main__":
    # Columbus, OH lat2, lon2 = 39.91, 116.56 # Beijing print spdist(lat1, lon1, lat2, lon2)
    lat1, lon1 = 40, -83
    lat2, lon2 = 39.91, 116.56
    print(spherical_distance(lat1, lon1, lat2, lon2))