CGCS2000坐标转换问题的思考

背景：按照自然资源部要求，2018年7月1日后，自然资源系统将全面使用2000国家大地坐标系；自2019年1月1日起，全面停止向社会提供1954年北京坐标系和1980西安坐标系基础测绘成果。

CGCS2000定义示意图

作者通过两年多的培训班授课，与一线从业人员做了大量交流和研讨，总结了CGCS2000坐标转换的若干重点问题与解决方法，基本涵盖了生产实践中的具体问题。

现将一些思考和心得记叙如下：

1、不是一个软件就能解决问题

测绘相关行业的从业人员涉及到的专业主要包括：大地测量、工程测量、导航、遥感、GIS，甚至计算机专业。当前，遥感和GIS专业很实用、很热门。但是，无论如何，大地测量专业都是测绘行业的基石。了解和掌握一些大地测量基础知识，可以在生产实践中做到心中有数。

与遥感和GIS专业偏重于软件应用和开发不同，大地测量是一门科学。而大地测量最重要、最基础、最难的内容就是基准，以及基准转换问题。主要包括：坐标基准、高程基准、重力基准。遥感和GIS专业的一线从业人员习惯于用软件解决问题。学习CGCS2000坐标转换的目的也很直接，就是获得一个软件，并学会使用。

软件我们当然会拷贝给学员，而且是永久性的免费使用。然而，与以往不同。由空间测量技术实现的高精度地心坐标必须考虑板块运动影响，同时也大幅度增加了坐标转换的复杂性。在实践中会遇到各种各样的具体情况，并不是一个软件就能解决问题的。必须深入理解，才能做到举一反三。

在测绘行业中，大地测量理论难度最大，但无法回避。对于不同专业背景的从业人员，学习大地测量知识，可以侧重于知识体系的构建，避免纠缠于具体公式。采用总结、归纳、分类、比较的方式，尽量做到打消神秘感、理解原理、理顺逻辑、解决困惑、掌握应用。

2、行业内的两极分化现象

当前，行业内存在明显的两极分化现象。一方面，大量专家做着高精尖的科研项目；另一方面，一线从业人员满足于采用固定的工作流程完成项目。这两极存在较严重的隔离现象。

例如坐标转换问题，我们有很多研究成果，实际上并不实用。这样说可能有失偏颇，但是，最起码这些程序没有给大家推广使用，只是发表了一些文章。而很多一线从业人员仅仅就知道个七参数转换、四参数转换，能完成任务就行。《技术规范》给出了8个经典模型，每个都有自己的特点和适用条件。我们也应该有所取舍，不能一个七参数模型走遍天下。

控制点转换到2000国家大地坐标系的坐标转换模型选取

市面上有很多坐标系方面的书籍和教材，很难找到一本能站在一线从业人员的角度来详细描述的。一是，大家的很多困惑和问题，专家们并不了解。二是，解决这些问题不出成果，不能显示水平。三是，在专家们的意识里，这根本就不是问题，不值得探讨。你怎么会连这些都不懂呢？

如果不是在授课中与大量一线从业人员的广泛交流和探讨，作者也不可能想象出这么多的实际存在的、具有普遍性的问题。事实上，经过两年多的坐标系主题授课，作者发现，要把CGCS2000坐标转换讲透彻，并不是一件容易的事。其中牵扯到各个专业、各种情况、方方面面。具体的问题都可以单独解决，但要系统化、条理化，让人听得懂，能灵活运用，还真的是剪不断、理还乱。为了能清晰表达，还不得不定义了一些新的名词。这些名词虽然还没有得到官方认可，但是随着作者的推广，以及大家的普遍使用，也将会成为约定俗成的说法。

当然，也有一些培训班会把《技术指南》或《技术要求》照本宣科一遍，这是极不负责任的态度。无论机构或者个人，口碑是最重要的。

3、软件和硬件傻瓜化，人不能傻瓜化

现在的软件和硬件越来越智能化，也可以称为傻瓜化。软件和硬件的傻瓜化，也会直接导致人的傻瓜化。

例如，以前测地形图，需要布设控制网、测量图根点，还需要测量每个细部，然后还要内业数据处理、成图。现在只需要用RTK和无人机就可以很方便的完成任务。一个没有学过测量的人，也完全可以胜任大部分工作。遥感和GIS专业也存在这种问题，经过短时间培训就能熟练使用软件，懂不懂测绘好像也无关紧要。这些从很多学员提的外行问题中可以明确感觉到。

这就导致很多公司不愿意花更多的成本聘请专业人员。他们把CGCS2000坐标转换也看作是一个软件的问题。花钱报个培训班，派两个熟练工，学两天就要回去解决问题。这是不可能的，也不利于长远发展，对整个行业也是一种伤害。

再例如，很多软件可以直接用重合点的xyH坐标（平面坐标和高程）来求布尔莎七参数。但实际上用xyH坐标是不能直接求出布尔莎七参数的。在后台，软件用中央子午线将xyH坐标变成了BLH或XYZ坐标，然后才求出布尔莎七参数。但是这样就会给大家带来误解，在求七参数时产生很多困惑。这是一个非常普遍的问题。

因此，软件和硬件可以傻瓜化，人不能傻瓜化。应该保持对专业的尊重，重视对原理的理解。

4、一些容易混淆的概念

在我们的教材和课堂中，还在给大家讲着CIO和格林尼治天文台这些古典概念。事实上，我国从来没有使用过CIO。现代的地心坐标系使用的是IERS参考极IRP和IERS参考子午线IRM。

还有坐标系的尺度问题。为什么ITRF使用的是地球时TT，而CGCS2000使用的是地心坐标时TCG？尺度为什么由时间来定义？

CGCS2000的参考历元是2000.0，这个2000.0怎么来的？如何计算？和J2000.0有关系吗？

历元归算与框架转换

这些概念在很多书中只是简单罗列了一下概念，并没有给出明确和易于理解解释。看起来都很熟悉，其实不知道是怎么回事。

就连动态坐标系中最基础、最重要的历元和速度场的概念，大部分从业人员甚至都没听说过。

还有一个，我们现在关心的和迫切要解决的是生产中的坐标转换问题。但是很多培训班一上来就给大家讲天球参考系，又是岁差，又是章动的。我们做的是地球上的坐标转换，天球参考系是研究卫星轨道时才用的。要说拓展知识面吧，关键是很多人对地球参考系还稀里糊涂的，他会把天球参考系和地球参考系的概念混淆。就不止一个人问过我，地心坐标的历元归算是不是岁差引起的？

况且，天球参考系也是一个过时的、古典的概念。现在使用的是与地球本体无关的四维时空参考系，其定义与天极、赤道、黄道都无关。

5、西部地区的特殊性

很多理论和方法在东部地区可以很好的实现，但是到了西部地区就会出问题。除了自然条件以外，往往不是人的水平问题，高程因素是一个重要原因。

例如独立坐标系投影,用不同软件投影出来的平面坐标不同，差异随抵偿高程面的增加而变大，在西藏甚至差10米。

另外，在西藏，在中央子午线一致的前提下，加了抵偿高程面和没加抵偿高程面的平面坐标x值相差接近3公里。

再例如，在ArcGIS的布尔莎七参数坐标转换过程中，转换精度和测区平均大地高相关。西藏地区仅因大地高而引起的平面坐标转换误差就可能超过10cm。

这些都是客观存在的问题。我们首先应该承认它，直面这些问题，给出理论解释，提出解决方案。虽不能改变现状，但能打消大家的疑惑，做到心中有数，进而灵活处理问题。

6、常常被忽视的高程转换

遥感和GIS专业更多的处理的是平面地图，对于高程经常是一知半解，甚至认为XYZ坐标中的Z值就是高程！

高程具有区域性和复杂性特点。有不同的高程系统，还有不同的基准。85基准高程还有'99水准网'和'15水准网'两种不同的实现版本。不了解这些概念的话，就会出现诸如：为什么千寻测出的高程差好几米？这样的问题。

高程具有区域性和复杂性特点

高程转换的复杂程度并不比坐标系少。如果精度要求不高，就不需要测水准，直接用程序计算高程异常。

7、莫名其妙的54坐标系转换

按照《技术要求》，北京54坐标转换到CGCS2000坐标分两步走：北京54坐标先转换为西安80坐标，再把西安80坐标转换为CGCS2000坐标。

按照《技术指南》，北京54坐标转换到西安80坐标也需要分两步：首先计算大地坐标改正量，然后计算平差改正量。用这两个改正量就可以将北京54坐标转换到西安80坐标。

坐系参数

仅依据《技术要求》和《技术指南》根本无法理解这个莫名其妙的转换流程，原因在于这两个文件中都没有提及'1954北京坐标系（整体平差转换值）'这个坐标系。引入这个坐标系以后，整个转换过程就非常清晰了。

明白了原理，解决了困惑，这样还不够。用户往往不具备资源条件，例如现在很难找到诸如平差改正数分布图这些资料了。因此，我们还需要提出可行的解决方案。

8、令人纠结的WGS84坐标系

大量文献中都探讨了WGS84坐标系，并与CGCS2000坐标系做了详细比较。但是于事无补，我们想知道WGS84坐标与CGCS2000坐标到底差多少？关于这个问题，各种说法都有。有些网络RTK给用户既提供WGS84坐标，也提供CGCS2000坐标。有些WGS84坐标与CGCS2000坐标差几十公分，有些则很接近。WGS84坐标与CGCS2000坐标的差异不能一概而论，需要从历元、框架、精度和实现四个方面去比较，就可以得到明确的结论了。

WGS84坐标系与北斗坐标系并没有本质区别，为什么北斗坐标系不存在这些问题呢？这又是一个历史的误会。

过去把GPS测得的坐标都叫做WGS-84坐标，这种观念在早期码伪距单点定位时是正确的。而且由于其十米多的定位精度，也不需要考虑坐标的时变性。后来发展了高精度的相对定位和精密单点定位技术，这时解出的坐标已经不是WGS-84坐标了。但是这种观念和习惯一直延续至今。

WGS-84坐标如何转换到CGCS2000坐标？也不能一概而论，你必须先了解你的WGS-84坐标是如何得到的。或者通过与CGCS2000坐标比较，推断出WGS-84坐标的历元。最起码要保证一组WGS-84坐标是同一时期的坐标。这样才能合理选择采用动态或静态的方法转换。而《技术规范》里关于WGS-84坐标转换的叙述有明显错误。

CGS2000椭球与GRS80和WGS84椭球推导几何参数比较

WGS-84坐标与CGCS2000坐标到底是什么关系？很简单，CGCS2000坐标是约定了历元和框架的ITRS坐标，WGS-84坐标则是没有约定历元和框架的ITRS坐标。同北斗坐标系一样，WGS-84坐标系只是一个卫星导航坐标系。

9、在学习中要有批判精神

现在的一些期刊、教科书，甚至规范中也存在错漏和不合理之处。我们在学习中一定要有批判精神，不能全信。有时候可能并不是你自己的错，但是只有自己搞懂了，才能心中有数、灵活处理。

某些资料中提到，西安80坐标转到CGCS2000坐标，椭球面积每平方公里应该减小0.939平方米。这个结论是用椭球面积公式计算出来的，只考虑了椭球大小变化对面积的影响，没有考虑椭球的径向平移，也没有考虑尺度缩放。因此，没有参考价值。事实上有些地方是每平方公里增加9平方米。这就需要一个更合理的估算公式。

在分丘图、宗地图、房屋平面图的坐标转换中，要求坐标转换到CGCS2000后面积不能改变。但是无论是坐标转换，还是实测CGCS2000坐标，面积必然会改变。那么只能牺牲精度来实现等面积转换了。

按照《技术规范》要求，动态转换时，应该先做历元归算，后做框架转换。对于2020.0历元的坐标，如果先做框架转换，后做历元归算，就会发现两种不同转换顺序得到的坐标差约10cm。从逻辑上看，两个动态转换顺序都没有问题。那么，这10cm差异从哪里来？

速度是用坐标拟合出来的，因此速度也有对应的参考框架。《技术规范》里的速度场模型忽略了速度的参考框架。就会导致这10cm的差异。

结束语：

CGCS2000坐标转换千头万绪，会有各种情况。限于篇幅，这里只能列出一些普遍的共性问题。总之，只有理解原理，才能举一反三。