兰勃特等角圆锥(Lambert Conformal Conic)投影正反变换

1. 引言

Johann Heinrich Lambert(译为兰勃特,或兰伯特),瑞士裔德国科学家、哲学家,他首次给出了π为无理数的严格证明。1772年提出两种球面向投影面投影方式:等角圆锥投影和等积方位投影。通常所说的兰勃特投影指的是等角圆锥投影。

兰勃特等角圆锥投影是设想用一个正圆锥切于或割于球面,应用等角条件将地球面投影到圆锥面上,然后沿一母线展开成平面。投影后纬线为同心圆圆弧,经线为同心圆半径。兰勃特等角圆锥投影没有角度变形,经线长度比和纬线长度比相等。适于制作沿纬线分布的中纬度地区中、小比例尺地图。我国的分省图的兰伯特等角圆锥投影采用的两条标准纬度线为Q1=25度,Q2=45度。

兰勃特等角圆锥(Lambert Conformal Conic)投影正反变换_第1张图片

 兰勃特等角圆锥(Lambert Conformal Conic)投影正反变换_第2张图片

 在北半球(左)和南半球(右)上以标准纬线显示了兰勃特等角圆锥投影。

2. 投影变换中所用符号

北方基准纬线:\phi_N

南方基准纬线:\phi_S

原点处纬度:\phi_0

原点处经度:\lambda_0

变换点纬度:\phi

变换点经度:\lambda

原点在投影坐标系中虚北向坐标:N_0

原点在投影坐标系中虚东向坐标:E_0

变换点假北向坐标:N

变换点假东向坐标:E

椭圆的长半轴:a=6378137米(WGS-84椭球体)

椭圆的短半轴:b=6356752.3142米(WGS-84椭球体)

椭圆的扁率:\alpha = (a-b)/a

椭圆的第一偏心率:

e=\sqrt{\frac{a^2-b^2}{a^2}}=\sqrt{2\alpha - \alpha^2}

其它中间参数:

        Q_S = \frac{1}{2}\left [ \ln\left (\frac{1+\sin\phi_S}{1-\sin\phi_S} \right ) - e \ln\left (\frac{1+e\sin\phi_S}{1-e\sin\phi_S} \right ) \right ]

        W_S = \sqrt{1-e^2\sin^2\phi_S}

        Q_N = \frac{1}{2}\left [ \ln\left (\frac{1+\sin\phi_N}{1-\sin\phi_N} \right ) - e \ln\left (\frac{1+e\sin\phi_N}{1-e\sin\phi_N} \right ) \right ]

        W_N = \sqrt{1-e^2\sin^2\phi_N}

        Q_0 = \frac{1}{2}\left [ \ln\left (\frac{1+\sin\phi_0}{1-\sin\phi_0} \right ) - e \ln\left (\frac{1+e\sin\phi_0}{1-e\sin\phi_0} \right ) \right ]

        \phi_t = \sin^{-1}\left [ \frac{\ln{W_N\cos\phi_S} - \ln{W_S\cos\phi_N}}{Q_N - Q_S} \right ]

        K = \frac{a\cos\phi_S\exp(Q_S\sin\phi_t)}{W_S\sin\phi_t}

        R_0 = \frac{K}{\exp(Q_0\sin\phi_t)}  

3. 正变换(\large (\lambda, \phi) \Rightarrow (N,E)

        Q = \frac{1}{2}\left [ \ln\left (\frac{1+\sin\phi}{1-\sin\phi} \right ) - e \ln\left (\frac{1+e\sin\phi}{1-e\sin\phi} \right ) \right ] 

        R = \frac{K}{\exp(Q\sin\phi_t)}

        \gamma = -(\lambda_0 - \lambda)\sin\phi_t

        N = N_0 + R_0 - R\cos\gamma

        E = E_0 + R\sin\gamma

 4. 逆变换(\large (N,E) \Rightarrow (\lambda, \phi)

        R_t = R_0 - (N - N_0)

        E_t = E - E_0

        \gamma = \tan^{-1}\left(\frac{E_t}{R_t} \right )

被求点的经度为:

        \lambda = \lambda_0 + \frac{\gamma}{\sin\phi_t}

被求点的纬度需要进行迭代运算获得:

        R = \sqrt{R_t^2 + E_t^2}

        Q = \frac{\ln(K/R)}{\sin\phi_0}

step 1: 估计纬度\large \phi初值

        \phi = \sin^{-1}\left(\frac{\exp(2Q) - 1}{\exp(2Q) + 1} \right )

setp 2: 计算修正量\large -f_1/f_2

        f_1 = \frac{1}{2}\left [ \ln\left (\frac{1+\sin\phi}{1-\sin\phi} \right ) - e \ln\left (\frac{1+e\sin\phi}{1-e\sin\phi} \right ) \right ] - Q

        f_2 = \frac{1}{1-\sin^2\phi} - \frac{e^2}{1-e^2\sin^2\phi}

step 3: 更新纬度\large \phi估值

         \phi = \sin^{-1}\left ( \sin\phi - f_1/f_2\right )

重复step 2 与 3 两次以上,直至收敛。

你可能感兴趣的:(电磁波覆盖范围分析,数据库)