子午线收敛角python代码

最新推荐文章于 2024-04-12 12:30:42 发布

泽几98

最新推荐文章于 2024-04-12 12:30:42 发布

阅读量1k

点赞数

分类专栏：自记文章标签：几何学算法计算机视觉

本文链接：https://blog.csdn.net/baidu_41639643/article/details/122371126

版权

自记专栏收录该内容

1 篇文章 0 订阅

订阅专栏

子午线收敛角是指通过该点的子午线投影与过该点的纵坐标线之间的夹角，规定西偏为正，东偏为负。在投影带的中央经线以东的图幅均为东偏，以西的图幅均为西偏。
子午线收敛角（convergence of meridian，Gauss grid convergence）是地球椭球体面上一点的真子午线与位于此点所在的投影带的中央子午线之间的夹角。即在高斯平面上的真子午线与坐标纵线的夹角。
此角有正、负之分：以真子午线北方向为准，当坐标纵轴线北端位于以东时称东偏，其角值为正；位于以西时称西偏，其角值为负。
子午线收敛角示意图
大地坐标计算：
在这里插入图片描述

def Magnetic_A_Norm(B,L,L0,e2_2):    
    #  大地坐标计算子午线收敛角
    # B：纬度
    # L：经度 
    # L0：中央子午线 
    # e2_2：第二离心率的平方 
    # 1980年国家大地坐标系
    t=tan(radians(B))
    l=L-L0
    l=radians(l)
    eta_2=e2_2*(cos(radians(B)**2))
    gamma=sin(radians(B))*l*((1+float(1/3)*(1+3*(eta_2)+2*(eta_2**2))*(cos(radians(B))**2)*(l**2))+float(1/15)*(2-(t**2))*(cos(radians(B))**4)*(l**4))

    return math.degrees(gamma)