Python根据经纬度坐标点计算该点正东正西正南正北一定距离的坐标点
业务需求需要计算经纬度点正东正西正南正北的一定距离的经纬度点,但没有百度到实现的函数,自己又温习了弧长公式,球的弧长公式计算,推导过程很简单,百度也能找到,最终的公式结果为:
假设地球是个完美的球体,半径为R,地球周长2πR,需要计算经纬度的距离为dist。
正北和正南距离dist 的经度不变,纬度分别加减 (dist* 360)/2πR 度
正东和正西距离dist 的纬度不变,经度分别加减 (dist* 360)/2π(cos(纬度)* R) 度
其他代码只要套用这个公式就可以实现
Python的代码实现如下
计算点经纬度正北500米的点的经纬度
from math import radians, cos, sin, asin, sqrt
def get_new_lat(lng1, lat1, dist=500):
"""
:param lng1: 116.498079
:param lat1: 39.752304
:param dist:需要计算的距离,默认500米
:return: (116.498079,39.756801)
"""
#6371 * 1000 为地球半径,单位米
lat2 = 180 * dist / (6371 * 1000 * pi) + lat1
return (lng1, lat2)
计算点经纬度正东500米的点的经纬度
def get_new_lng(lng1, lat1, dist=500):
"""
:param lng1: 116.498079
:param lat1: 39.752304
:param dist:需要计算的距离,默认500米
:return: (116.503928,39.752304)
"""
lng2 = 180 * dist / (6371 * 1000 * pi * cos(radians(lat1))) + lng1
return (lng2, lat1)
```