GPS坐标间距离计算

# -*- coding:utf-8 -*-
# /usr/bin/python

import warnings
warnings.filterwarnings("ignore")
from math import radians, cos, sin, asin, sqrt

# 方法1
def geodistance(lat1,lng1,lat2,lng2):
    ''''数学转换'''
    lng1, lat1, lng2, lat2 = map(radians, [float(lng1), float(lat1), float(lng2), float(lat2)]) # 经纬度转换成弧度
    dlon=lng2-lng1
    dlat=lat2-lat1
    a=sin(dlat/2)**2 + cos(lat1) * cos(lat2) * sin(dlon/2)**2
    distance=2*asin(sqrt(a))*6371*1000 # 地球平均半径，6371km
    distance=round(distance/1000,3)
    return distance

# 方法2
from geopy.distance import geodesic

def fromGeodesic(lat1,lng1,lat2,lng2):
    '''调用包'''
    distance=geodesic((lat1,lng1), (lat2,lng2))
    return distance

# 方法3
def simpleDistance(lat1,lng1,lat2,lng2):
    dx = lng1 - lng2 # 经度差值
    dy = lat1 - lat2 # 纬度差值
    b = (lat1 + lat2) / 2.0 # 平均纬度
    Lx = radians(dx) * 6367.0 * cos(radians(b))# 东西距离
    Ly = 6367.0 * radians(dy) #南北距离
    distance = sqrt(Lx * Lx + Ly * Ly)
    return distance



print(geodistance(39.941,116.45,39.94,116.451))
print(fromGeodesic(39.941,116.45,39.94,116.451))
print(simpleDistance(39.941,116.45,39.94,116.451))

描述
radians() 方法将角度转换为弧度。

语法
以下是 radians() 方法的语法:

角度和弧度关系是：2π 弧度 = 360°。从而 1°≈0.0174533 弧度，1 弧度≈57.29578°。

1) 角度转换为弧度公式：弧度=角度÷180×π
2) 弧度转换为角度公式： 角度=弧度×180÷π
import math

math.radians(x)
注意：radians()是不能直接访问的，需要导入 math 模块，然后通过 math 静态对象调用该方法。

参数
x -- 一个角度数值，默认单位是角度 °。
返回值
返回一个角度的弧度值。

arduino 可以参考上面代码实现GPS坐标间距离测量