关于GPS坐标转换（二）

 

    上一章已经简单说明了各坐标系的基本概念和相互的转换关系（其实也不叫学习笔记了，所有转换都做完了，只是一直在忙，没有时间一次写完，只能有时间就写一章，一点一点来，实在愧对关注的朋友）。

    这一章就说一下84坐标系下的地心大地坐标如何使用84坐标系的椭球参数转换为84坐标系下的地心直角坐标（好拗口……）。

    先来看看84坐标系的椭球参数：

        椭圆长半径A:6378137.000000

        椭圆短半径B: 6356752.314245179497

    根据椭圆长短半径，我们很容易就能算出椭圆第一偏心率e1，计算公式如下：

       e1 ＝ sqrt(A×A－B×B)/ A；

    注意这里sqrt是开方的意思。

    OK ，上面我们通过一个简单的公式就计算出了椭圆的第一 偏心率e1，下面，我们就利用这个参数把 地心大地坐标转换为地心直角坐标。

    已知：B、L、H地心大地坐标

    求：X、Y、Z地心直角坐标

    计算步骤：

1． 参数N：N＝A/sqrt（1－e1×e1×sin（B）×sin（B））；

2． X ：X＝（N＋H）×cos（B）×cos（L）；

3． Y ：Y＝（N＋H）×cos（B）×sin（L）；

4． Z ：Z＝（N×（1－e1×e1）＋H）×sin（B）；

    好了，就是这么简单，BLH方式表示的大地坐标一下就转换为XYZ表示的直角坐标了。不过大家要注意的是，就是BL这两个数值是角度值，在C编程中，cos和sin函数要用弧度值，这个角度转换问题要注意。