现有测绘成果转换到2000国家大地坐标系技术指南
一、2000国家大地坐标系的定义
国家大地坐标系的定义包括坐标系的原点、三个坐标轴的指向、尺度以及地球椭球的4个基本参数的定义。2000国家大地坐标系的原点为包括海洋和大气的整个地球的质量中心;2000国家大地坐标系的Z轴由原点指向历元2000.0的地球参考极的方向,该历元的指向由国际时间局给定的历元为1984.0的初始指向推算,定向的时间演化保证相对于地壳不产生残余的全球旋转,X轴由原点指向格林尼治参考子午线与地球赤道面(历元2000.0)的交点,Y轴与Z轴、X轴构成右手正交坐标系。采用广义相对论意义下的尺度。2000国家大地坐标系采用的地球椭球参数的数值为:
长半轴a=6378137m
地心引力常数GM=3.986004418×1014m3s-2
自转角速度ω=7.292l15×10-5rads-1
其它参数见下表:
采用2000国家大地坐标系后仍采用无潮汐系统。
二、点位坐标转换方法
(一)模型选择
全国及省级范围的坐标转换选择二维七参数转换模型;省级以下的坐标转换可选择三维四参数模型或平面四参数模型。对于相对独立的平面坐标系统与2000国家大地坐标系的联系可采用平面四参数模型或多项式回归模型。坐标转换模型详见本指南第六部分。
(二)重合点选取
坐标重合点可采用在两个坐标系下均有坐标成果的点。但最终重合点还需根据所确定的转换参数,计算重合点坐标残差,根据其残差值的大小来确定,若残差大于3倍中误差则剔除,重新计算坐标转换参数,直到满足精度要求为止;用于计算转换参数的重合点数量与转换区域的大小有关,但不得少于5个。
(三)模型参数计算
用所确定的重合点坐标,根据坐标转换模型利用最小二乘法计算模型参数。
(四)精度评估与检核
用上述模型进行坐标转换时必须满足相应的精度指标,具体精度评估指标及评估方法见附件中相关内容。选择部分重合点作为外部检核点,不参与转换参数计算,用转换参数计算这些点的转换坐标与已知坐标进行比较进行外部检核。应选定至少6个均匀分布的重合点对坐标转换精度进行检核。
(五)数据库中点位坐标转换模型参数计算的区域选取
对于1980西安坐标系下的数据库,采用全国数据计算的一套模型参数可满足1:5万及1:25万比例尺数据库转换的精度要求;采用全国数据计算的六个分区的模型参数可满足1:1万比例尺数据库转换的精度要求。对于1954年北京坐标系下的数据库的转换,采用全国数据计算的六个分区的模型参数可满足1:5万及1:25万比例尺数据库转换的精度要求;按(2°×3°)进行分区计算模型参数可满足1:1万比例尺数据库转换的精度要求。
三、1:2.5-1:25万数据库的转换
(一)按国家基本比例尺地形图分幅组织的数据库
按国家基本比例尺地形图分幅组织的图形数据(DLG、DEM、DRG),依据以下方案进行转换。
1、1:2.5-1:10万DLG数据库转换
(1)1954年北京坐标系下1:2.5-1:10万DLG数据库转换
a、依据相应比例尺分幅进行区域划分,分两步完成坐标转换。首先进行椭球体变换,再利用对应的比例尺图幅区域的X、Y坐标平移量进行坐标平移;
b、依据2000国家大地坐标系下对应的比例尺标准分幅图廓进行数据裁切,区域边缘图幅中的数据空白区利用相邻图幅数据进行补充;
c、添加2000国家大地坐标系下的方里格网层,删除原方里格网数据层;
d、完成图廓更改、数据编辑、数据接边、拓扑重建、数据入库等数据后处理及建库工作;
e、图幅换带接边:采用右图(1954年北京坐标系)接左图(2000国家大地坐标系)时,先进行右图的椭球体与换带转换,在左带中利用左图的平移量进行右图的坐标平移,完成接边后保存在左带中的右图(备份)成果。返回右图取消先前换带接边加入的平移量,并进行投影变换,最后利用右带自身的平移量完成平移后,方可与其相邻的右图接边;
f、对基础地理信息数据库元数据相关条目进行更改。
(2)1980西安坐标系下1:2.5-1:10万DLG数据库转换
依据相应比例尺分幅进行区域划分,不考虑椭球体变换,直接利用对应的比例尺图幅区域的X、Y坐标平移量进行坐标平移;然后按照1954年北京坐标系下DLG数据库转换的b~f对应步骤进行。
2、1:2.5-1:10万DRG数据库转换
原数据为300~500dpi的原版印刷地图经扫描纠正生成的RGB栅格数据,无图幅间要素的接边处理。
(1)1954年北京坐标系下1:2.5-1:10万DRG数据库转换
a、考虑椭球变换及对应图廓角点的X、Y坐标平移量,计算1954年北京坐标系分幅图廓角点在2000国家大地坐标系下的坐标,并修改数据头文件中相应的定位坐标;
b、在DRG数据上叠加2000国家大地坐标系下新的大地控制基础层(图廓及方里格网等),新图廓中数据空白或数据出图区域不做图纹补充和裁减;
c、在图例中添加2000国家大地坐标系下新的控制基准说明条款;
d、完成数据合层,并保持DRG数据的原有分辨率;
e、更改元数据中相关内容,增加1954年北京坐标系标准分幅的图廓四角点在2000国家大地坐标系下坐标,计算2000国家大地坐标系标准分幅的图廓四角点的坐标。
转换后数据为2000国家大地坐标系坐标、1954年北京坐标系分幅。
(2)1980西安坐标系下1:2.5-1:10万DRG数据库转换
获取图幅对应比例尺图幅图廓角点的X、Y坐标平移量,根据平移量计算图幅定位坐标,修改数据头文件;然后按照1954年北京坐标系到2000国家大地坐标系的1:2.5-1:10万DRG数据库转换的b~e步骤进行。
转换后数据为2000国家大地坐标系坐标、1980年西安坐标系分幅。
3、1:2.5-1:10万DEM数据库转换
原数据为25米分辨率的灰阶(256个)栅格数据,建库数据图幅间接边处理完好。此数据转换可有两种方式:一种是依据2000国家大地坐标系下DLG相关图层数据(等高线、高程点)重新生成DEM(见DEM数据生产规范),一种是进行DEM数据的转换。以下给出DEM数据转换方法。
(1)1954年北京坐标系下1:2.5-1:10万DEM数据库转换
a、按照比例尺对应图幅分块,在需补充内容的邻接边各增加一个相应比例尺图幅;
b、考虑椭球变换及相应的比例尺图幅的X、Y坐标平移量,求得X、Y坐标改正值;
c、根据坐标改正值进行图幅坐标平移,同时,参考像素分辨率确定起算坐标进行数据重采样;
d、按2000国家大地坐标系新的图廓及重叠像素进行图幅裁切,更改数据头文件中定位坐标;
e、修改元数据相关条目。
(2)1980西安坐标系下1:2.5-1:10万DEM数据库转换
a、按照相应比例尺对应图幅分块,在需补充内容的邻接边各增加一个相应比例尺图幅;
b、依据相应的比例尺图幅的X、Y坐标平移量,进行图幅坐标平移,并参考像素分辨率确定起算坐标完成数据重采样;
c、d按1954年北京坐标系1:2.5-1:10万DEM数据库转换的d、e步骤进行。
4、1:25万DLG数据库转换
(1)将1:25万分幅的平面坐标平移量转换为对应的经、纬度平移量或直接获取对应图幅的经、纬度平移量;
(2)根据1:25万分幅的经、纬度平移量,完成1:25万经纬度数据到2000国家大地坐标系经纬度数据的转换(1954年北京坐标系需同时考虑椭球体变化和平移量);
(3)依据2000国家大地坐标系下对应的1:25万标准分幅图廓进行数据裁切,区域边缘图幅中的数据空白区利用相邻图幅数据进行补充;
(4)数据后处理,包括:图廓更改、新格网层添加、数据编辑、数据接边、拓扑重建、数据入库等;
(5)更改元数据文件。
5、1:25万DEM数据库转换
(1)利用2000国家大地坐标系对应的DLG数据层,重新内插生成DEM;
(2)依据新的DEM更改元数据文件。
(二)按其它方式建立的数据库
1、按区域建立的图形数据库
按区域(省、地区、流域等)建立的图形数据库(DLG、DEM、DRG),可先分带分块分层完成转换,参照以上相应比例尺基础地理信息数据库的转换方案转换后拼接合成。
1:10万-1:25万数据库,依1:25万数据库转换方案逐块进行转换,再整体拼接合成;按非高斯投影方式组织的,将原数据经纬网30′×30′或15′×15′交点作为坐标转换参考点,计算这些参考点在2000国家大地坐标系下的坐标,利用地理信息软件进行图形纠正,完成数据转换。
2、按线性条带建立的图形数据库
按线性条带(境界、河流、交通线、管道线等)建立的图形数据库,可依据条带的方向、长短等分段进行,再拼接合成;也可通过条带中一定密度地物点的两套坐标,通过软件逐点进行纠正。具体方法:
分块纠正:对于1:1万分块,按1:1万数据转换方案逐块纠正后接边合成;对于1:5万分块,按1:2.5-1:10万数据转换方案逐块纠正后接边合成;
逐点纠正:依据数据精度,建立一定密度(1:1万100米格网点、1:5万2000米格网点)的坐标转换参考点,计算这些参考点在新坐标系下的坐标,利用地理信息软件完成数据转换。
3、按无固定分幅分区建立的图形数据库
按无固定分幅分区建立的图形数据,根据坐标系、比例尺及数据主体所在的图幅、数据的组织方式、产品类型(DLG、DEM、DRG)等,参照相应比例尺的转换方案,实施数据转换。
4、DOM数据库转换
原数据为航空或航天遥感获取的黑白或彩色影像数据,是连续的灰度(全色)或RGB(彩色)栅格数据,分辨率有多种方式(主要包括用于1:5地形图测绘的各种分辨率航空影像,以及用于专题调查的10米、15米、30米等卫星影像)。影像数据转换可参照下列方式进行。
对于已按数据库组织方式加工与处理的DOM数据,可采用1:2.5-1:10万DEM的数据转换方法,也可采用计算各景影像有效图边的4点在2000国家大地坐标系下的坐标来重新定位的方式。
对于尚未按数据库组织方式加工与处理的DOM数据,可采用1:2.5-1:10万DRG的数据转换方法,不再添加新的控制基础信息。
分辨率5米-30米的数据,需依据其数据主体所在的1:25万图幅区域来选用1:25万对应图幅的综合坐标改正值;对于分辨率在2米到5米间的数据,需依据其数据主体所在的1:5万图幅区域来选用1:5万对应图幅的综合坐标改正值;由此确定各自的X、Y方向平移像素数对应的坐标值(直接取1:25万或1:5万综合坐标改正值,或由像素数×像素分辨率求得)。按高斯投影、分像对(分景)组织的高分辨率影像数据,参照1:1万DOM转换技术方案进行转换。
四、1:1万及1:5千基础地理信息数据库的转换
(一)1:1万及1:5千格网点坐标转换改正量计算
1、1980西安坐标系坐标转换改正量计算
1:1万以上大比例尺一般按(2°×3°)进行分区,并对每个分区向外扩充约20′,分别解算出各分区的转换参数后,利用确定的转换方法与转换模型分别计算全国1:1万及1:5千格网点的2000国家大地坐标系坐标B2000,L2000,进而求出各点的1980西安坐标系与2000国家大地坐标系的差值DB802000,DL802000(B2000-B80,L2000-L80),形成全国1:1万及1:5千格网点的1980西安坐标系与2000国家大地坐标系的转换改正量DB802000,DL802000。
2、1954年北京坐标系坐标转换改正量计算
全国1954年北京坐标系向2000国家大地坐标系转换改正量计算采用两步法:首先计算1954年北京坐标系转换向1980西安坐标系转换改正量,其次计算1980西安坐标系向2000国家大地坐标系转换改正量,最后将两改正量叠加形成1954年北京坐标系向2000国家大地坐标系转换坐标转换改正量。
①1954年北京坐标系向1980西安坐标系转换坐标改正量计算
新旧坐标系统(1954年北京坐标系与1980年西安坐标系)的转换
大地坐标改正量计算公式:
式中: △a,△e2分别为IAG-75椭球与克拉索夫斯基椭球长半径,第一偏心率平方之差。即。
则各个点在1980西安坐标系中的大地坐标为:
● 根据转换的B80,L80,采用高斯投影正算公式计算相应的高斯平面坐标 X80,Y80。
● 求取全国1:1万以大比例尺格网点的转换改正量
● 平差改正量的计算
1954年北京坐标系所提供的大地点成果没有经过整体平差,1980西安坐标系提供的大地点成果是经过整体平差的数据,所以新旧系统转换还要考虑平差改正量的问题。计算平差改正量比较麻烦,没有一定的数学模式,不同地区,平差改正量差别很大,在我国中部某些地区,平差改正量在1米以下,而在东北地区的某些图幅则在10米以上。在实际计算中,在全国均匀地选择一定数量的一、二等大地点,利用它们新(1980西安坐标系)旧(1954年北京坐标系)坐标系的坐标进行多种分析试算并剔除粗差点,然后分别计算它们的坐标差值,根据这些差值和它们的大地坐标分别绘制两张平差改正量分布图(即dX,dY分布图),这样在分布图上可以直接内插出全国1:1万以大比例尺格网点的平差改正量DX2,DY2。
● 根据全国1:1万以大比例尺格网点的转换改正量DX1,DY1和平差改正量DX2,DY2按下列公式计算1954年北京坐标系向1980西安坐标系转换坐标转换改正量DX,DY。
● 将DX,DY换算成1:1万以大比例尺格网点大地坐标转换改正量DB5480,DL5480。
②1954年北京坐标系向2000国家大地坐标系转换坐标转换改正量计算
将全国1:1万以大比例尺格网点的1954年北京坐标系向1980西安坐标系的转换改正量DB5480与1980西安坐标系向2000国家大地坐标系的转换改正量DB802000叠加,得到全国1:1万以大比例尺格网点1954年北京坐标系向2000国家大地坐标系转换的坐标转换改正量DB542000。
即:DB542000=DB5480+DB802000
DL542000=DL5480+DL802000
(二)1:1万及1:5千DLG数据库转换
转换流程如图1所示:
图1 1:1万及1:5千基础地理信息数据库转换技术流程
1、1954年北京坐标系下1:1万、1:5千DLG数据库转换
(1)每个图幅的四个图廓点坐标改正量选用1954年北京坐标系向2000国家大地坐标系转换方法计算;
(2)图幅内各要素点的坐标改正量根据选用的本图幅的四个图廓点坐标改正量,按双线性内插等方法计算;
(3)根据图幅四个图廓点坐标改正量和图幅内各要素点的坐标改正量,计算2000国家大地坐标系下的图幅四个图廓点坐标和图幅内各要素点的坐标;
(4)与周边图幅拼接;
(5)按照2000国家大地坐标系下对应1:1万、1:5千标准分幅计算新的公里格网数据,即添加2000国家大地坐标系下新的公里格网层;
(6)完成图廓更改、数据编辑、换带接边、拓扑重建;
(7)对空间数据库元数据相关条目进行更改;
(8)数据入库等数据后处理工作。
2、1980西安坐标系下1:1万、1:5千DLG数据库转换
(1) 每个图幅的四个图廓点坐标改正量选用1980西安坐标系向2000国家大地坐标系转换方法计算;
(2)~(8)参照1954年北京坐标系到2000国家大地坐标系1:1万、1:5千DLG数据库转换的对应步骤进行。
(三)1:1万及1:5千DRG数据库转换
在保持原分辨率不变的情况下,利用逐格网纠正的方法进行数据转换。
1、1954年北京坐标系下1:1万、1:5千DRG数据库转换
(1)在2000国家大地坐标系下生成图廓坐标及公里格网,逐公里格网点纠正1:1万、1:5千DRG数据;
(2)修改元数据相关条目;
(3)修改相关的图外整饰。
2、1980西安坐标系下1:1万、1:5千DRG数据库转换
参照1954年北京坐标系到2000国家大地坐标系1:1万、1:5千DRG数据库转换的对应步骤进行。
(四)1:1万及1:5千DEM数据库转换
利用DEM生产过程中形成的矢量数据与DEM离散点数据完成数据转换。
1、1954年北京坐标系下1:1万、1:5千DEM数据库转换
(1)矢量数据与DEM离散点数据转换方法;
a、每个图幅的四个图廓点坐标改正量选用1954年北京坐标系向2000国家大地坐标系转换方法计算;
b、图幅内各要素点的坐标改正量根据选用的本图幅的四个图廓点坐标改正量,按双线性内插等方法计算;
c、根据图幅四个图廓点坐标改正量和图幅内各要素点的坐标改正量,计算2000国家大地坐标系下的图幅四个图廓点坐标和图幅内各要素点的坐标;
d、与周边图幅拼接。
(2)构造TIN;
(3)按相关规范或规定内插DEM;
(4)对检查点坐标文件进行点对点坐标转换;
(5)修改元数据条目。
2、1980西安坐标系下1:1万、1:5千 DEM数据库转换
(1)每个图幅的四个图廓点坐标改正量选用1980西安坐标系向2000国家大地坐标系转换方法计算;
(2)~(5)参照1954年北京坐标系1:1万、1:5千 DEM数据库转换的对应步骤进行。
(五)1:1万及1:5千DOM数据库转换
原数据为航空或航天遥感获取的黑白或彩色影像数据,是连续的灰度或RGB栅格数据,在转换中应保持原影像分辨率。
1、在原DOM上叠加相应坐标系统的内图廓及公里格网,在2000国家大地坐标系下生成图廓坐标及公里格网,逐公里格网点纠正1:1万、1:5千DOM数据;
2、转换后,删除内图廓及公里格网;
3、修改元数据相关条目。
五、相对独立的平面坐标系统与2000国家大地坐标系建立联系的方法
(一)相对独立的平面坐标系统控制点建立联系的方法
可通过现行国家大地坐标系的平面坐标过渡,利用坐标转换方法将相对独立的平面坐标系统下控制点成果转换到2000国家大地坐标系下。
选取相对独立的平面坐标系统与2000国家大地坐标系的重合点的原则如下:择优选取地方控制网的起算点及高精度控制点、周围国家高精度的控制点,大中城市至少选取5个重合点(城外4个,市内中心1个);小城市在城市外围至少选取4个重合点,重合点要分布均匀,包围城市区域,并在城市内部选定至少6个均匀分布的重合点对坐标转换精度进行检核。
建立相对独立的平面坐标系统与2000国家大地坐标系联系时,坐标转换模型要同时适用于地方控制点转换和城市数字地图的转换。一般采用平面四参数转换模型,重合点较多时可采用多元逐步回归模型。当相对独立的平面坐标系统控制点和数字地图均为三维地心坐标时,采用Bursa七参数转换模型。坐标转换中误差应小于0.05米。
(二)相对独立的平面坐标系统下数字地形图转换
采用点对点转换法完成相对独立的平面坐标系统下数字地形图到2000国家大地坐标系的转换,转换后相邻图幅不存在接边问题。具体步骤如下:
利用控制点的转换模型和参数,对相对独立的平面坐标系统下数字地形图进行转换,形成2000国家大地坐标系地形图。
根据转换后的图幅四个图廓点在2000国家大地坐标系下的坐标,重新划分公里格网线,原公里格网线删除。
根据2000国家大地坐标系下的图廓坐标,对每幅图进行裁剪和补充。
六、坐标转换方法
(一)坐标转换模型
1、二维七参数转换模型
其中:
△B,△L---同一点位在两个坐标系下的纬度差、经度差,单位为弧度,
△a,△f ---椭球长半轴差(单位米)、扁率差(无量纲),
△X,△Y,△Z ---平移参数,单位为米,
εx,εy,εz---旋转参数,单位为弧度,
m---尺度参数(无量纲)。
2、平面四参数转换模型
属于两维坐标转换,对于三维坐标,需将坐标通过高斯投影变换得到平面坐标再计算转换参数。
平面直角坐标转换模型:
其中,x,y为平移参数,α为旋转参数,m为尺度参数。x2,y2为2000国家大地坐标系下的平面直角坐标,x1,y1为原坐标系下平面直角坐标。坐标单位为米。
3、综合法坐标转换
所谓综合法即就是在相似变换(Bursa七参数转换)的基础上,再对空间直角坐标残差进行多项式拟合,系统误差通过多项式系数得到消弱,使统一后的坐标系框架点坐标具有较好的一致性,从而提高坐标转换精度。
综合法转换模型及转换方法:
● 利用重合点先用相似变换转换Bursa七参数坐标转换模型
式中,3个平移参数 [△X△Y△Z]T,3个旋转参数[εxεyεz]T和1个尺度参数 m。
● 对相似变换后的重合点残差Vx,Vy,Vz采用多项式拟合
式中:B,L单位:弧度;K为拟合阶数;aij为系数,通过最小二乘求解。
4、三维七参数坐标转换模型
(二)高斯正反算公式
1、高斯投影正算公式
子午线弧长X计算见附录。
2、高斯投影反算公式
(三)坐标转换精度评定和评估方法
对于1954年北京坐标系、1980西安坐标系与2000国家大地坐标系转换分区转换及数据库转换点位的平均精度应小于图上的0.1mm。具体:
对于1:5千坐标转换,1980西安坐标系与2000国家大地坐标系转换分区转换平均精度≤0.5m;1954年北京坐标系与2000国家大地坐标系转换分区转换平均精度≤1.0m;
1:5万基础地理信息数据库坐标转换精度≤5.0m;
1:1万基础地理信息数据库坐标转换精度≤1.0m;
1:5千基础地理信息数据库坐标转换精度≤0.5m。
依据计算坐标转换模型参数的重合点的残差中误差评估坐标转换精度。对于n个点,坐标转换精度估计公式如下:
①V(残差)=重合点转换坐标-重合点已知坐标
②空间直角坐标X残差中误差
③空间直角坐标Y残差中误差
④空间直角坐标Z 残差中误差
点位中误差
⑤平面坐标x残差中误差
⑥平面坐标y残差中误差
⑦大地高H残差中误差
平面点位中误差为
2000 国家大地坐标系转换常见问题分析
摘要:在现有成果向 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 坐标部分省建立基础控制网采用的 GCS2000 起算 点坐标非国家测绘行政主管部门权威发布,或不能 量值溯源到国家测绘行政主管部门权威数据。
表 1 为某省 C 级控制网建立时采用的 IGS 站或 国家级 GNSS 基准站的 CGCS2000 坐标与国家权威 数据的较差统计值,这一差异导致转换数据整体出 现系统性偏差。
②起算点精度等级达不到相应等级控制点要求 部分省市坐标基准框架或基础控制网的建立未与 IGS 站或国家级 GNSS 基准站进行联测,仅与本 省市及周边省市 2000 国家 GPS 大地控制网( 三网 平差) 控制点进行联测,经约束平差获得 CGCS2000 坐标。这种方法在对 C、D 级 GPS 控制点进行转换 时确保了成果与本省区域内的其他成果的一致性。 省级 GNSS 基准站作为省级大地基准的骨干和主要 支撑,采用这种方法在当时历史条件下和过渡期内 实现了 CGCS2000 在省级的快速推广使用,但在现 阶段不利于维持省级三维、动态地心坐标系统,不利 于保证大地控制网点位三维地心坐标的精度、现势 性及全国的统一。
表 2 为某省采用坐标归算方法以 A、B 级点为 基准建立的 56 个省级 GNSS 基准站 CGCS2000 坐标 的外符合精度情况统计,可以明显看出在空间三维 方向上均存在一定的误差,并且均具有系统性偏差。 这些误差已经对 GNSS 基准站的服务造成了一定 影响。
③未考虑框架不同历元间由于板块运动引起的坐标变化值
部分省在进行历元归算后,未考虑计算框架所 对应历元下坐标从观测历元到需转换历元之间,由 于板块运动引起的坐标变化值,把坐标变化值带入 到转换成果中,引起转换成果误差。
ITRF 2005 转换到 ITRF 2000 框架时站的速度 场起主要作用,因此若所确定的速度场不准确对转 换结果有很大的影 响。 而 从 ITRF 2000 转 换 到 ITRF 97 框架下起主要作用的是框架之间的转换关 系,对所需转换的站的速度场要求不是很高[13]。
表 3 为部分点不同年代观测数据联合平差时,因板块运动而引起的坐标变化值无法消除,最终转换成果产生的系统误差。
( 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 转换精度损失,进而转换 精度超限。
表 4 为某省不同地形类别区域的 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 推广应用中的生 产和质量检查工作提供参考,确保成果转换的质量。
来源:兰冶工业设计研究院
------ End ------
领取专属 10元无门槛券
私享最新 技术干货