def get_angle_between_points(point1, point2, point3):
m_point1 = pyproj.transform(PROJ_WGS_84, PROJ_MERCATOR,
point1[LON], point1[LAT])
m_point2 = pyproj.transform(PROJ_WGS_84, PROJ_MERCATOR,
point2[LON], point2[LAT])
m_point3 = pyproj.transform(PROJ_WGS_84, PROJ_MERCATOR,
point3[LON], point3[LAT])
v1x = (m_point1[LON] - m_point2[LON]) # / COORDINATE_PRECISION
v1y = m_point1[LAT] - m_point2[LAT]
v2x = (m_point3[LON] - m_point2[LON]) # / COORDINATE_PRECISION
v2y = m_point3[LAT] - m_point2[LAT]
angle = (atan2(v2y, v2x) - atan2(v1y, v1x)) * 180. / pi
while angle < 0:
angle += 360.
return angle
评论列表
文章目录