80、54、84坐标系七参数转换算法

2019独角兽企业重金招聘Python工程师标准>>>

一、为什么要进行坐标转换

   我们所在地球是一个不规则的椭球，地表凹凸不平，地底密度不均，因此很难用一个简单模型来概括。国际上根据建模坐标系的原点不同分为参心坐标系和地心坐标系，其中参心坐标系是指参考椭球的几何中心坐标基准，地心，通常分为：参心空间直角坐标系（以x，y，z为其坐标元素）和参心大地坐标系（以B，L，H为其坐标元素），比如我国的北京54坐标系和西安80坐标系；坐标系是用地球质心作为原点，通常分为地心空间直角坐标系(以x，y，z为其坐标元素)和地心大地坐标系(以B，L，H为其坐标元素)，比如GPS采用的WGS84坐标系。

   除了原点以外，任何一个椭球的模型还会有长轴、短轴和扁率三个参数，我整理了下54,80,84三个坐标的椭球模型如下表:

坐标系	长轴（单位：米）	短轴（单位：米）	扁率
北京54坐标系	6378245	6356863	1/298.3
西安80坐标系	6378140	6356755	1/298.25722101
WGS84坐标系	6378137	6356752.314	1/298.25722356

 

        所以，我们可见这三个坐标系是原点不一，椭球模型不一的坐标系，就如同三个达芬奇的鸡蛋，虽然都是椭球状，但是形状还是不一样的。所以如果直接做两个坐标系之间的转换是明显不行的。

   那么有没有办法进行转换呢？答案是肯定的。在数学上，任何两个空间直角坐标系都可以通过七个参数，即XYZ三个轴的平移，绕XYZ三个轴的旋转角以及两个坐标系的大小比例因子的转换后重叠在一起。因此，我们在做着三个坐标系的转换基本思路就是在空间直角坐标系里，通过七参数的偏移校正，实现坐标系之间参数的转换。

 二、怎样坐标转换

        好了，那我们如何转换呢。首先，我再区分一下地理上的空间直角坐标系、大地坐标系和平面坐标系的区别。

   空间直角坐标系和数学的空间直角坐标系基本相同，是以一个点为原点，三个互相垂直的直线作为三个轴延伸出来一个坐标系，通常表示为（x,y,