2000国家大地坐标系转换常见问题分析

张训虎，刘晋虎，何 川，邱 博，杨绪峰，李 淼

(国家测绘产品质量检验测试中心，北京 100830)

2000国家大地坐标系转换常见问题分析

张训虎，刘晋虎，何 川，邱 博，杨绪峰，李 淼

(国家测绘产品质量检验测试中心，北京 100830)

在现有成果向2000国家大地坐标系转换工作中，各省市都做了很多理论研究和实际转换工作。本文阐述了现有成果向2000国家大地坐标系转换的方法，对不同数据、不同方法转换过程中常见的问题进行汇总、整理，并分析问题产生的原因及对成果的影响，为今后现有成果向2000国家大地坐标系转换工作提供参考和建议，以确保成果转换的质量。

CGCS2000；坐标转换；大地控制点；基础地理信息数据；GNSS基准站；三角点；4D产品

自2008年7月1日起，我国启用2000国家大地坐标系(CGCS2000)，各地有关部门开展了现有各类测绘信息成果到CGCS2000的转换工作，积极推进CGCS2000的推广使用。为做好启用CGCS2000的实施工作，国家测绘地理信息局于2008年7月组织编制了《启用2000国家大地坐标系实施方案》和《现有测绘成果转换到2000国家大地坐标系技术指南》。为加快CGCS2000推广使用工作，进一步指导各部门、各单位顺利完成原有坐标系成果向CGCS2000转换，确保转换技术方法正确，国家测绘地理信息局于2013年6月组织编制了《2000国家大地坐标系推广使用技术指南》和《大地测量控制点坐标转换技术规程》。CGCS2000转换涉及的测绘地理信息成果主要包括大地控制点类成果(GNSS基准站、GNSS控制点、三角点)和基础地理信息数据成果(DLG、DOM、DEM、DRG)。文献[1—2]从总体上介绍了CGCS2000的构建、精化、维持和更新，以及可用于转换工作的国家级成果。文献[3—11]研究了GNSS基准站、大地控制点、4D产品的转换方法。文献[12] 探讨了转换的检查方法。本文对省级坐标转换中存在的常见问题进行梳理和分析。

一、坐标转换方法

1. 大地控制点类成果

(1) 坐标归算

本方法适用于非CGCS2000的省级GNSS基准站和卫星大地控制点坐标向CGCS2000的坐标转换。即对拟转换点采用与周边稳定的IGS站及国家级的GNSS连续运行基准站进行联测的方法，获得各站点在现ITRF框架、观测历元下的位置，经过历元归算、板块运动改正、框架转换[13]3个步骤进行坐标计算。用这种方法进行转换必须要知道网中各站的速度场。

(2) 参数转换

本方法适用于未联测已知点的卫星大地控制点和三角点坐标向CGCS2000的坐标转换。即按照转换区域选取适当的转换模型(布尔莎模型、三维七参数模型、二维七参数模型、三维四参数模型、二维四参数模型、多项式拟合模型等)，选择重合点，经粗差剔除后计算转换参数，进行坐标转换。

2. 基础地理信息数据成果

基础地理信息数据成果转换针对分幅数据或数据库实体数据，采用不同的坐标转换方法。比例尺大于1∶1万的DLG、DEM数据一般采用高分辨率的格网坐标改正量进行逐要素点转换的方法完成转换；DOM、DRG数据一般采用平移或纠正的方法完成转换，转换参数一般采用高分辨率的格网坐标改正量进行表达。

二、常见问题及分析

1. 大地控制点类成果

(1) 坐标归算

采用坐标归算方法进行坐标转换的关键工序主要有：基准控制点(IGS站、国家级GNSS基准站、国家GNSS大地控制点)的选取、高精度数据处理、板块运动改正、框架转换等。坐标归算常见问题主要集中在基准控制点的选取和板块运动改正两个方面。

①起算点坐标非国家测绘行政主管部门权威发布的CGCS2000坐标

部分省建立基础控制网采用的CGCS2000起算点坐标非国家测绘行政主管部门权威发布，或不能量值溯源到国家测绘行政主管部门权威数据。

表1为某省C级控制网建立时采用的IGS站或国家级GNSS基准站的CGCS2000坐标与国家权威数据的较差统计值，这一差异导致转换数据整体出现系统性偏差。

表1 起算点坐标较差统计 m

②起算点精度等级达不到相应等级控制点要求

部分省市坐标基准框架或基础控制网的建立未与IGS站或国家级GNSS基准站进行联测，仅与本省市及周边省市2000国家GPS大地控制网(三网平差)控制点进行联测，经约束平差获得CGCS2000坐标。这种方法在对C、D级GPS控制点进行转换时确保了成果与本省区域内的其他成果的一致性。省级GNSS基准站作为省级大地基准的骨干和主要支撑，采用这种方法在当时历史条件下和过渡期内实现了CGCS2000在省级的快速推广使用，但在现阶段不利于维持省级三维、动态地心坐标系统，不利于保证大地控制网点位三维地心坐标的精度、现势性及全国的统一。

表2为某省采用坐标归算方法以A、B级点为基准建立的56个省级GNSS基准站CGCS2000坐标的外符合精度情况统计，可以明显看出在空间三维方向上均存在一定的误差，并且均具有系统性偏差。这些误差已经对GNSS基准站的服务造成了一定影响。

表2 某省级GNSS基准站网外符合精度统计 m

③未考虑框架不同历元间由于板块运动引起的坐标变化值

部分省在进行历元归算后，未考虑计算框架所对应历元下坐标从观测历元到需转换历元之间，由于板块运动引起的坐标变化值，把坐标变化值带入到转换成果中，引起转换成果误差。

ITRF 2005转换到ITRF 2000框架时站的速度场起主要作用，因此若所确定的速度场不准确对转换结果有很大的影响。而从ITRF 2000转换到ITRF 97框架下起主要作用的是框架之间的转换关系，对所需转换的站的速度场要求不是很高[13]。

表3为部分点不同年代观测数据联合平差时，因板块运动而引起的坐标变化值无法消除，最终转换成果产生的系统误差。

表3 不同年代同名点坐标差异 m

(2) 参数转换

采用参数转换方法进行坐标转换的关键工序主要有：转换分区、转换模型的选取、重合点的选择和剔除、转换参数计算、外部检核等。参数转换常见问题主要集中在转换分区、重合点覆盖范围、重合点的剔除、检核点分布4个方面。

1) 由于1954北京坐标系的坐标是采用逐级控制分区平差的方法推算的，存在明显的平差变形，甚至个别地区在分区或锁网接合部点出现了成果不一致或产生了裂缝[14]，因此在这类地区不宜采用一个分区和一套转换参数，以避免产生较大误差。

2) 重合点选取的基本原则为等级高、精度高、局部变形小、分布均匀及覆盖整个转换区域。当重合点不能覆盖整个转换区域时，不能覆盖的区域转换参数只能通过外推得到，但转换精度可能随外推距离放大而急剧损失，导致转换后的成果与邻省成果间存在不接边的情况。对于从国家申领的具有CGCS2000坐标的一、二、三、四等天文大地网点，不加区别全部用于转换模型的计算，造成了重合点利用的等级和精度不统一，转换精度不高，局部变形较大。

3) 粗差点剔除不严密，不严格按照大于3倍点位中误差进行，易导致局部转换参数的变形。粗差点的剔除还应包括造成重合点分布不均匀的点，如应最大限度避免模型中狭长三角形的出现，这种点可作为外部检核点使用。

4) 利用未参与计算转换参数的重合点作为外部检核点，其点数应不少于6个且分布均匀。外部检核点不足时应进行野外实测检核，尤其应注意对转换区域边缘的检核。

2. 基础地理信息数据成果

(1) DEM转换

由于生产DEM成果的过程数据(等高线、特征线、高程点等)一般不存在，DEM转换不能按照相关生产技术规程构造TIN并内插重新生成DEM，一般选用高分辨率格网坐标改正量并采用平移或双线性内插的方法对图幅进行坐标转换，同时参考像素分辨率确定起算坐标进行数据重采样，按CGCS2000新的图廓及重叠像素进行图幅裁切，更改数据头文件中定位坐标，修改元数据相关条目。DEM转换常见问题主要有以下几个方面：

1) 采用平移方法进行DEM转换，以图幅4个角点平移量的平均值作为图幅左下角点改正量，不进行数据重采样，DEM数据仍以原坐标系图廓范围进行单幅存储。

这种方法的图幅起始点坐标为非格网间距的整数倍，因相邻图幅坐标平移量不一致导致图幅不接边。在后期DEM数据应用时，接边区域内高程仍需处理，并造成重采样精度损失。

2) 采用平移方法进行DEM转换，坐标平移量归整化为DEM格网间距的整数倍，不进行数据重采样。

这种方法会产生DEM局部相邻图幅间相差一排(一列)DEM格网点，导致局部图幅接边处格网数值不唯一，出现少一排(一列)或重合一排(一列)的情况(如图1所示 )。因坐标平移量规整化为格网点间距整数倍，导致DEM转换精度损失，进而转换精度超限。

图1 DEM转换相邻图幅间差一列

表4为某省不同地形类别区域的DEM转换精度统计，可见这种方法在山区容易导致部分图幅转换精度超限。

表4 1∶1万DEM转换精度统计

(2) DOM转换

DOM转换一般选用高分辨率格网坐标改正量采用平移或纠正的方法对图幅进行坐标转换，按CGCS2000规定的新的图廓及重叠像素进行图幅裁切，按像素关系计算移动量(像素数)，更改数据头文件中定位坐标，修改元数据相关条目。DOM转换常见问题主要有以下几个方面。

DOM转换过程中将平移量规整化为DOM地面采样间距的整数倍后对整图进行坐标平移，以及局部相邻图幅间相差一排(一列)DOM栅格点，导致局部相邻图幅接边区域数值不唯一。这种转换方法虽不会对DOM转换精度造成重大影响，但转换工作并未全面完成，宜对接边成果进行重采样处理，完善转换工作。

三、建议及措施

1. 大地控制点类成果

1) 平差计算过程中的起算控制点CGCS2000成果不能仅利用向国家申领的2000国家GPS大地控制网成果(三网平差成果，地心坐标精度平均优于3 cm)中的大地点成果，需要更加充分利用精度更高的2000国家GPS大地控制网中的GNSS连续运行基准站坐标(坐标精度为毫米级)。

2) 在坐标归算过程中顾及板块运动的特性和不同历元间框架的严格转换关系，充分利用可用于转换工作的国家级最新速度场成果CGCS2000板块运动模型(China Plate Model-CGCS2000，CPM-CGCS2000)和CGCS2000格网速度场模型。CPM-CGCS2000是目前国内最精确的相关模型，解决了CGCS2000动态维持及我国已有基础测绘成果转换的难题，适用于基于ITRF框架非2000历元下各类GNSS定位成果到CGCS2000的转换。

3) 各省与邻省进行重合点数据交换，建立参数转换模型的重合点尽可能覆盖全部转换区域；不能覆盖的转换区域转换参数可通过平滑外推得到，但要加强外部检核工作和邻省的接边工作。

4) 采用不同转换模型进行比较分析，绘制点位残差分布图和点位残差等值线图，选择最优模型进行坐标转换。采用多项式拟合模型进行坐标转换，还应参考布尔莎模型、二维七参数转换模型等适合于全国及省级范围的转换模型进行精度分析，剔除残差较大点。

2. 基础地理信息数据成果

现有成果转换工作量大，且各省现有基础地理信息数据成果为过渡性成果，随着基础测绘工作的持续开展，新的基础地理信息数据成果宜直接采用CGCS2000生产，各省市对现有成果的转换工作，以满足实际应用为目的进行，转换过程中应保证转换数据的完整性、一致性、唯一性，确保转换到位、接边到位。

加强转换工作技术方案和技术路线的全面质量评估，避免数据转换出现重大质量问题。

四、结束语

现有成果向CGCS2000转换工作是一项系统工作，在转换过程中，尽管各省市结合自身情况，开展了一系列的理论研究和实际转换工作，但是在转换中仍然存在一些问题。本文对检查中发现的转换工作相关问题进行整理分析，旨在为今后的CGCS2000转换工作、为CGCS2000推广应用中的生产和质量检查工作提供参考，确保成果转换的质量。

[1] 陈俊勇,杨元喜,王敏，等. 2000国家大地控制网的构建和它的技术进步[J]. 测绘学报,2007,36(1):1-8.

[2] 宁津生,王华,程鹏飞,等.2000国家大地坐标系框架体系建设及其进展[J].武汉大学学报(信息科学版),2015,40(5):569-573.

[3] 王文利,程传录,李东,等.1954年北京坐标系和1980西安坐标系精度分析及其适用性研究[J].大地测量与地球动力学,2012,32(5):68-71，77.

[4] 成英燕,程鹏飞,秘金钟,等.基于现框架下的省市级CORS站到CGCS2000的转换[J].测绘通报,2011(7):1-3，14.

[5] 成英燕,程鹏飞,顾旦生,等.三维4参数模型实现地图到CGCS2000的转换[J].武汉大学学报(信息科学版),2010,35(6):747-751，755.

[6] 郭春喜,王文利,白贵霞,等.坐标系转换中全国高精度高分辨率格网改正量的确定[J].测绘科学,2013,38(2):5-7.

[7] 吕志平,魏子卿,李军,等.CGCS2000高精度坐标转换格网模型的建立[J].测绘学报,2013,42(6):791-797.

[8] 成英燕,程鹏飞,秘金钟,等.大尺度空间域下1980西安坐标系与WGS-84坐标系转换方法研究[J].测绘通报,2007(12):5-8.

[9] 韩买侠,郭春喜,王文利,等.小比例尺基础数据由1980坐标系向2000坐标系转换的方法[J].测绘科学,2014,39(1):32-34，58.

[10] 庞尚益,郭春喜,程传录.国家基本比例尺DLG数据坐标转换方法的研究[J].测绘科学,2006,31(S1):28-29.

[11] 郭春喜,韩买侠. 数字高程模型(DEM)和数字线划图(DLG)的坐标转换方法[J].测绘通报,2013(1):57-59.

[12] 张训虎.现有成果转换为2000国家大地坐标系成果检查方法探讨[J].工程勘察,2015(2):85-89.

[13] 程鹏飞,成英燕,文汉江,等.2000国家大地坐标系实用宝典[M].北京：测绘出版社,2008：130-132.

[14] 董鸿闻. 1954年北京坐标系的历史注释[J]. 东北测绘,2001,24(2):16-17，23.

AnalysisofCommonProblemsinCGCS2000Conversion

ZHANG Xunhu，LIU Jinhu，HE Chuan，QIU Bo，YANG Xufeng，LI Miao

张训虎，刘晋虎，何川，等.2000国家大地坐标系转换常见问题分析[J].测绘通报，2016(9)：52-55.

10.13474/j.cnki.11-2246.2016.0291.

P282

B

0494-0911(2016)09-0052-04

2015-11-04；

2016-05-24

测绘地理信息公益性行业科研专项经费(201512018)

张训虎(1977—)，男，高级工程师，主要从事测绘产品质量检验工作及大地测量、工程测量研究。E-mail：zhangxh@sbsm.gov.cn