我国现有测绘成果转换到2000国家大地坐标系参考

－12－现有测绘成果转换到2000国家大地坐标系技术指南一、2000国家大地坐标系的定义国家大地坐标系的定义包括坐标系的原点、三个坐标轴的指向、尺度以及地球椭球的4个基本参数的定义。2000国家大地坐标系的原点为包括海洋和大气的整个地球的质量中心；2000国家大地坐标系的Z轴由原点指向历元2000.0的地球参考极的方向，该历元的指向由国际时间局给定的历元为1984.0的初始指向推算，定向的时间演化保证相对于地壳不产生残余的全球旋转，X轴由原点指向格林尼治参考子午线与地球赤道面（历元2000.0）的交点，Y轴与Z轴、X轴构成右手正交坐标系。采用广义相对论意义下的尺度。2000国家大地坐标系采用的地球椭球参数的数值为：长半轴a＝6378137m扁率f=1/298.257222101地心引力常数GM＝3.986004418×1014m3s-2自转角速度ω＝7.292l15×10-5rads-1其它参数见下表：短半径b(m)6356752.31414极曲率半径c(m)6399593.62586第一偏心率e0.0818191910428第一偏心率平方e20.00669438002290第二偏心率e0.0820944381519第二偏心率平方e20.006739496775481/4子午圈的长度Q(m)10001965.7293椭球平均半径R1(m)6371008.77138相同表面积的球半径R2(m)6371007.18092相同体积的球半径R3(m)6371000.78997椭球的正常位U0(m2s-2)62636851.7149动力形状因子J20.001082629832258球谐系数J4-0.00000237091126球谐系数J60.00000000608347球谐系数J8-0.0000000000142722/mabGM0.00344978650678－13－赤道正常重力值γe（伽）9.7803253361两极正常重力值γp（伽）9.8321849379正常重力平均值γ（伽）9.7976432224纬度45度的正常重力值γ45°（伽）9.8061977695采用2000国家大地坐标系后仍采用无潮汐系统。二、点位坐标转换方法（一）模型选择全国及省级范围的坐标转换选择二维七参数转换模型；省级以下的坐标转换可选择三维四参数模型或平面四参数模型。对于相对独立的平面坐标系统与2000国家大地坐标系的联系可采用平面四参数模型或多项式回归模型。坐标转换模型详见本指南第六部分。（二）重合点选取坐标重合点可采用在两个坐标系下均有坐标成果的点。但最终重合点还需根据所确定的转换参数，计算重合点坐标残差，根据其残差值的大小来确定，若残差大于3倍中误差则剔除，重新计算坐标转换参数，直到满足精度要求为止；用于计算转换参数的重合点数量与转换区域的大小有关，但不得少于5个。（三）模型参数计算用所确定的重合点坐标，根据坐标转换模型利用最小二乘法计算模型参数。（四）精度评估与检核用上述模型进行坐标转换时必须满足相应的精度指标，具体精度评估指标及评估方法见附件中相关内容。选择部分重合点作为外部检核点，不参与转换参数计算，用转换参数计算这些点的转换坐标与已知坐标进行比较进行外部检核。应选定至少6个均匀分布的重合点对坐标转换精度进行检核。（五）数据库中点位坐标转换模型参数计算的区域选取－14－对于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、对基础地理信息数据库元数据相关条目进行更改。－15－（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步骤进行。－16－转换后数据为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万分幅的平面坐标平移量转换为对应的经、纬度平移量或直接获取对应图幅的经、纬度平移量；－17－（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、按线性条带建立的图形数据库按线性条带（境界、河流、交通线、管道线等）建立的图形数据库，可依据条带的方向、长短等分段进行，再拼接合成；也可通过条带中一定密度地物点的两套坐标，通过软件逐点进行纠－18－正。具体方法：分块纠正：对于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转换技术方案进行转换。－19－四、1:1万及1:5千基础地理信息数据库的转换（一）1:1万及1:5千格网点坐标转换改正量计算1、1980西安坐标系坐标转换改正量计算1:1万以上大比例尺一般按(2°×3°)进行分区，并对每个分区向外扩充约20′，分别解算出各分区的转换参数后，利用确定的转换方法与转换模型分别计算全国1:1万及1:5千格网点的2000国家大地坐标系坐标B2000，L2000，进而