CGCS2000坐标转换问题的思考(转)

本文从其他网站转来,对CGCS2000至一些问题可以有更深的认识,文末有原始连接

作者简介:中科地信首席专家刘光明研究员,研究方向为空间大地测量。曾获国家科技进步二等奖,排名第三。应测绘出版社约稿,作为独立作者撰写《CGCS2000坐标转换》等专著两部。应多个学术机构邀请,在各省会城市做“CGCS2000坐标转换”和“北斗导航卫星系统”等专题学术报告80多次。

专栏作者 刘光明 | 撰文

王臻 | 编辑

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

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

1

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

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

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

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

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

2

行业内的两极分化现象

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

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

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

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

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

3

软件和硬件傻瓜化,人不能傻瓜化

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

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

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

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

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

4

一些容易混淆的概念

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

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

CGCS2000的参考历元是2000.0[6],这个2000.0怎么来的?如何计算?和J2000.0有关系吗?

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

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

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

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

5

西部地区的特殊性

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

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

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

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

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

6

常常被忽视的高程转换

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

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

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

7

莫名其妙的54坐标系转换

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

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

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

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

8

令人纠结的WGS84坐标系

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

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

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

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

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

9

在学习中要有批判精神

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

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

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

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

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

结束语:

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

参考文献

[1]大地测量控制点坐标转换技术规范.测绘行业标准:CH/T 2014—2016[S].北京:国家测绘地理信息局,2016.

[2]现有测绘成果转换到2000国家大地坐标系技术指南[S].北京:国家测绘地理信息局,2008.

[3]国土资源数据2000国家大地坐标系转换技术要求[S].北京:国家测绘地理信息局,2017.

[4]吕志平,乔书波.大地测量学基础[M].北京:测绘出版社,2016.

[5]赵铭.天体测量学导论[M].北京:中国科学技术出版社,2012.

[6]2000中国大地测量系统.国家军用标准:GJB 6304—2008[S].北京:中国人民解放军总装备部,2008.

[7]王若璞,赵东明译.时空参考系[M].北京:科学出版社,2015.

[8]城市坐标系统建设规范.国家标准:GB/T 28584—2012[S].北京:国家质量监督检验检疫总局,2012.

[9]王文利等.全国一等水准点高程近20年变化分析[J].测绘学报,2019,48(1):1-8.

[10]刘光明等.山区正常高测量新方法[J].测绘通报,2007,(6):27-29.

[11]郑祖良.大地坐标系的建立与统一[M].北京:解放军出版社,1993.

[12]魏子卿.2000中国大地坐标系及其与WGS-84的比较[J].大地测量与地球动力学,2008,28(5):1-5.

[13]北斗卫星导航系统空间信号接口控制文件公开服务信号B1C、B2a(测试版)[S].北京:中国卫星导航系统管理办公室,2017.

[14]刘光明.四维坐标系速度场建模技术研究[M].兰州大学硕士论文,2012年6月.

本文来自:中科地信(北京)遥感信息技术研究院《自然资源科技学术精要总集》,
这个来源为其他网站标注的来源,可惜的是这个来源里并没有找到这个总集,甚至这篇文章。
有关这个总集,网上只有这个一篇文章

发表回复

您的电子邮箱地址不会被公开。 必填项已用 * 标注

鄂ICP备号2023007360