geom.py 文件源码

python
阅读 31 收藏 0 点赞 0 评论 0

项目:pypgroutingloader 作者: danieluct 项目源码 文件源码
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
评论列表
文章目录


问题


面经


文章

微信
公众号

扫码关注公众号