标签：__ cos 经纬度 latitude1 距离 180 计算 pi math

根据经纬度计算距离是有固定公式的，按照公式写个函数很简单。

import math
def cal_dis(lat1, lon1,lat2, lon2):
    latitude1 = (math.pi/180)*lat1
    latitude2 = (math.pi/180)*lat2
    longitude1 = (math.pi/180)*lon1
    longitude2= (math.pi/180)*lon2
    #因此AB两点的球面距离为:{arccos[sinb*siny+cosb*cosy*cos(a-x)]}*R
    #地球半径
    R = 6378.137
    d = math.acos(math.sin(latitude1)*math.sin(latitude2)+ math.cos(latitude1)*math.cos(latitude2)*math.cos(longitude2-longitude1))*R
    return d


if __name__ == '__main__':
    print cal_dis(23.0,101.1,23.06,113.34)

注意，经纬度是角度，而三角函数的输入是弧度，角度与弧度的关系  180=pi

有一次我在一个项目中发现 matlab 的计算距离和上面这个函数有些出入，虽然相差不大，但是确实不一样，当时也折腾了好久，后来发现 python 有个自带的函数计算结果和 matlab 相同。

from geopy.distance import geodesic

geodesic((30.28708,120.12802999999997), (28.7427,115.86572000000001)).m

标签：__,cos,经纬度,latitude1,距离,180,计算,pi,math
来源： https://www.cnblogs.com/yanshw/p/11382501.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。