Python坐标计算（Python坐标计算三角形周长）

from math import radians, cos, sin, asin, sqrt,pi

class Get_new_gps():

def __init__(self):

# 地球半径

self.R = 6371 * 1000

pass

def geodistance(self, lng1, lat1, lng2, lat2):

"""计算两点间距离"""

"""

:param lng1: 120.12802999999997

:param lat1: 30.28708,115

:param lng2: 115.86572000000001

:param lat2: 28.7427

:return: 两点间距离

"""

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)) * self.R # 地球平均半径，6371km

distance = round(distance, 3)

return distance

def get_new_lat(self, lng1, lat1, dist=500):

"""计算点经纬度南500米的点的经纬度"""

"""

:param lng1: 120.12803

:param lat1: 20.28708

:param dist:

:return: (120.12803, 20.291576608029594)

"""

lat2 = lat1 - 180 * dist / (self.R * pi)

return (lng1, lat2)

"""计算点经纬度正东1000米点的经纬度"""

def get_new_lng(self, lng1, lat1, dist=1000):

"""计算点经纬度正东1000米点的经纬度"""

"""

:param lng1: 116.55272514141352

:param lat1: 30.28708

:param dist:

:return: (116.55272514141352,30.28708)

"""

lng2 = 180 * dist / (self.R * pi * cos(radians(lat1))) + lng1

return (lng2, lat1)

def get_new_lng_angle(self, lng1, lat1, dist=500, angle=30):

"""计算点东北方向与正东某夹角某距离的经纬度"""

"""

:param lng1:116.55272514141352

:param lat1:30.28708

:param dist:指定距离

:param angle:指定角度

:return:（0.0091871843081617/pi + 116.498079 0.0122339171779312/pi + 39.752304）

"""

lat2 = 180 * dist*sin(radians(angle)) / (self.R * pi) + lat1

lng2 = 180 * dist*cos(radians(angle)) / (self.R * pi * cos(radians(lat1))) + lng1

return (lng2, lat2)

functions = Get_new_gps()

if __name__ == '__main__':

functions = Get_new_gps()

lng1, lat1 = [117.317732, 36.565515]

# 计算正北的点

# lng2, lat2 = functions.get_new_lat(lng1, lat1)

# 计算正东的点

lng3, lat3 = functions.get_new_lng(lng1, lat1)

# # 计算该点与正东逆时针夹角45度，距离500米点的经纬度值

# lng4, lat4 = functions.get_new_lng_angle(lng1, lat1, dist=1000, angle=45)

print("原始点的经纬度坐标", lng1, lat1)

# print("正北500米坐标点为%f,%f,距离计算为%f米" % (lng2, lat2, functions.geodistance(lng1, lat1, lng2, lat2)))

print("正东1000米坐标点为%f,%f,距离计算为%f米" % (lng3, lat3, functions.geodistance(lng1, lat1, lng3, lat3)))

# print("东北方夹角，距离500米坐标点为%f,%f,距离计算为%f米" % (float(lng3), float(lat3), functions.geodistance(lng1, lat1, lng3, lat3)))

# print(functions.geodistance(lng3, lat3, lng4, lat4))

lng, lat = [117.774649, 36.005881]

lng2, lat2 = functions.get_new_lng(lng, lat)

print("原始点的经纬度坐标", lng, lat)

print("正东1000米坐标点为%f,%f,距离计算为%f米" % (lng2, lat2, functions.geodistance(lng, lat, lng2, lat2)))

更多学习资料请关注公众号“读书同友会”