1、12附件：现有测绘成果转换到 2000 国家大地坐标系技术指南一、2000 国家大地坐标系的定义国家大地坐标系的定义包括坐标系的原点、三个坐标轴的指向、尺度以及地球椭球的4个基本参数的定义。2000国家大地坐标系的原点为包括海洋和大气的整个地球的质量中心；2000国家大地坐标系的Z轴由原点指向历元2000.0的地球参考极的方向，该历元的指向由国际时间局给定的历元为1984.0的初始指向推算，定向的时间演化保证相对于地壳不产生残余的全球旋转，X轴由原点指向格林尼治参考子午线与地球赤道面（历元2000.0）的交点，Y轴与Z 轴、X轴构成右手正交坐标系。采用广义相对论意义下的尺度。2000国家大地坐
2、标系采用的地球椭球参数的数值为：长半轴 a6378137m扁率 f=1/298.257222101 地心引力常数 GM3.98600441810 14m3s-2 自转角速度 7.292l1510 -5rad s-113其它参数见下表：短半径 b(m) 6356752.31414极曲率半径 c (m) 6399593.62586第一偏心率 e 0.0818191910428第一偏心率平方 e2 0.00669438002290第二偏心率 0.0820944381519第二偏心率平方 2e 0.006739496775481/4 子午圈的长度 Q(m) 10001965.7293椭球平均半径 R1
3、(m) 6371008.77138相同表面积的球半径 R2(m) 6371007.18092相同体积的球半径 R3(m) 6371000.78997椭球的正常位 U0(m2s-2) 62636851.7149动力形状因子 J2 0.001082629832258球谐系数 J4 -0.00000237091126球谐系数 J6 0.00000000608347球谐系数 J8 -0.000000000014272/mabGM0.00344978650678赤道正常重力值 e（伽） 9.7803253361两极正常重力值 p（伽） 9.8321849379正常重力平均值 （伽） 9.79764322
4、24纬度 45 度的正常重力值 45（伽） 9.806197769514采用 2000 国家大地坐标系后仍采用无潮汐系统。二、点位坐标转换方法（一）模型选择全国及省级范围的坐标转换选择二维七参数转换模型；省级以下的坐标转换可选择三维四参数模型或平面四参数模型。对于相对独立的平面坐标系统与 2000 国家大地坐标系的联系可采用平面四参数模型或多项式回归模型。坐标转换模型详见本指南第六部分。（二）重合点选取坐标重合点可采用在两个坐标系下均有坐标成果的点。但最终重合点还需根据所确定的转换参数，计算重合点坐标残差，根据其残差值的大小来确定，若残差大于 3 倍中误差则剔除，重新计算坐标转换参数，直到满足
5、精度要求为止；用于计算转换参数的重合点数量与转换区域的大小有关，但不得少于 5 个。（三）模型参数计算用所确定的重合点坐标，根据坐标转换模型利用最小二乘法计算模型参数。 （四）精度评估与检核用上述模型进行坐标转换时必须满足相应的精度指标，具体精度评估指标及评估方法见附件中相关内容。选择部分重合15点作为外部检核点，不参与转换参数计算，用转换参数计算这些点的转换坐标与已知坐标进行比较进行外部检核。应选定至少 6 个均匀分布的重合点对坐标转换精度进行检核。（五）数据库中点位坐标转换模型参数计算的区域选取对于 1980 西安坐标系下的数据库，采用全国数据计算的一套模型参数可满足 1:5 万及 1:2
6、5 万比例尺数据库转换的精度要求；采用全国数据计算的六个分区的模型参数可满足 1:1 万比例尺数据库转换的精度要求。对于 1954 年北京坐标系下的数据库的转换，采用全国数据计算的六个分区的模型参数可满足 1:5 万及 1:25 万比例尺数据库转换的精度要求；按（23 ）进行分区计算模型参数可满足 1:1 万比例尺数据库转换的精度要求。三、1:2.5-1:25 万数据库的转换（一）按国家基本比例尺地形图分幅组织的数据库按国家基本比例尺地形图分幅组织的图形数据（DLG、 DEM、DRG ） ，依据以下方案进行转换。1、1:2.5-1:10 万 DLG 数据库转换（1）1954 年北京坐标系下 1
7、:2.5-1:10 万 DLG 数据库转换a、依据相应比例尺分幅进行区域划分，分两步完成坐标转换。首先进行椭球体变换，再利用对应的比例尺图幅区域的X、Y 坐标平移量进行坐标平移；16b、依据 2000 国家大地坐标系下对应的比例尺标准分幅图廓进行数据裁切，区域边缘图幅中的数据空白区利用相邻图幅数据进行补充；c、添加 2000 国家大地坐标系下的方里格网层，删除原方里格网数据层；d、完成图廓更改、数据编辑、数据接边、拓扑重建、数据入库等数据后处理及建库工作；e、图幅换带接边：采用右图（1954 年北京坐标系）接左图（2000 国家大地坐标系）时，先进行右图的椭球体与换带转换，在左带中利用左图的平
8、移量进行右图的坐标平移，完成接边后保存在左带中的右图（备份）成果。返回右图取消先前换带接边加入的平移量，并进行投影变换，最后利用右带自身的平移量完成平移后，方可与其相邻的右图接边；f、对基础地理信息数据库元数据相关条目进行更改。（2）1980 西安坐标系下 1:2.5-1:10 万 DLG 数据库转换依据相应比例尺分幅进行区域划分，不考虑椭球体变换，直接利用对应的比例尺图幅区域的 X、Y 坐标平移量进行坐标平移；然后按照 1954 年北京坐标系下 DLG 数据库转换的bf 对应步骤进行。2、1:2.5-1:10 万 DRG 数据库转换原数据为 300500dpi 的原版印刷地图经扫描纠正生成的
9、17RGB 栅格数据，无图幅间要素的接边处理。（1）1954 年北京坐标系下 1:2.5-1:10 万 DRG 数据库转换a、考虑椭球变换及对应图廓角点的 X、Y 坐标平移量，计算 1954 年北京坐标系分幅图廓角点在 2000 国家大地坐标系下的坐标，并修改数据头文件中相应的定位坐标；b、在 DRG 数据上叠加 2000 国家大地坐标系下新的大地控制基础层（图廓及方里格网等） ，新图廓中数据空白或数据出图区域不做图纹补充和裁减；c、在图例中添加 2000 国家大地坐标系下新的控制基准说明条款；d、完成数据合层，并保持 DRG 数据的原有分辨率；e、更改元数据中相关内容，增加 1954 年北京
10、坐标系标准分幅的图廓四角点在 2000 国家大地坐标系下坐标，计算 2000国家大地坐标系标准分幅的图廓四角点的坐标。转换后数据为 2000 国家大地坐标系坐标、1954 年北京坐标系分幅。（2）1980 西安坐标系下 1:2.5-1:10 万 DRG 数据库转换获取图幅对应比例尺图幅图廓角点的 X、Y 坐标平移量，根据平移量计算图幅定位坐标，修改数据头文件；然后按照1954 年北京坐标系到 2000 国家大地坐标系的 1:2.5-1:10 万DRG 数据库转换的 be 步骤进行。18转换后数据为 2000 国家大地坐标系坐标、1980 年西安坐标系分幅。3、1:2.5-1:10 万 DEM
11、数据库转换原数据为 25 米分辨率的灰阶（256 个）栅格数据，建库数据图幅间接边处理完好。此数据转换可有两种方式：一种是依据 2000 国家大地坐标系下 DLG 相关图层数据（等高线、高程点）重新生成 DEM（见 DEM 数据生产规范）,一种是进行DEM 数据的转换。以下给出 DEM 数据转换方法。（1）1954 年北京坐标系下 1:2.5-1:10 万 DEM 数据库转换a、按照比例尺对应图幅分块，在需补充内容的邻接边各增加一个相应比例尺图幅；b、考虑椭球变换及相应的比例尺图幅的 X、Y 坐标平移量，求得 X、Y 坐标改正值；c、根据坐标改正值进行图幅坐标平移，同时，参考像素分辨率确定起算
12、坐标进行数据重采样；d、按 2000 国家大地坐标系新的图廓及重叠像素进行图幅裁切，更改数据头文件中定位坐标；e、修改元数据相关条目。（2）1980 西安坐标系下 1:2.5-1:10 万 DEM 数据库转换a、按照相应比例尺对应图幅分块，在需补充内容的邻接19边各增加一个相应比例尺图幅；b、依据相应的比例尺图幅的 X、Y 坐标平移量，进行图幅坐标平移，并参考像素分辨率确定起算坐标完成数据重采样；c、 d 按 1954 年北京坐标系 1:2.5-1:10 万 DEM 数据库转换的 d、e 步骤进行。4、1:25 万 DLG 数据库转换（1）将 1:25 万分幅的平面坐标平移量转换为对应的经、纬
13、度平移量或直接获取对应图幅的经、纬度平移量；（2）根据 1:25 万分幅的经、纬度平移量，完成 1:25 万经纬度数据到 2000 国家大地坐标系经纬度数据的转换（1954 年北京坐标系需同时考虑椭球体变化和平移量） ；（3）依据 2000 国家大地坐标系下对应的 1:25 万标准分幅图廓进行数据裁切，区域边缘图幅中的数据空白区利用相邻图幅数据进行补充；（4）数据后处理，包括：图廓更改、新格网层添加、数据编辑、数据接边、拓扑重建、数据入库等；（5）更改元数据文件。5、1:25 万 DEM 数据库转换（1）利用 2000 国家大地坐标系对应的 DLG 数据层，重新内插生成 DEM;20（2）依据
14、新的 DEM 更改元数据文件。（二）按其它方式建立的数据库1、按区域建立的图形数据库按区域（省、地区、流域等）建立的图形数据库（DLG、 DEM、DRG ） ，可先分带分块分层完成转换，参照以上相应比例尺基础地理信息数据库的转换方案转换后拼接合成。1:10 万-1:25 万数据库，依 1:25 万数据库转换方案逐块进行转换，再整体拼接合成；按非高斯投影方式组织的，将原数据经纬网 3030或 1515交点作为坐标转换参考点，计算这些参考点在 2000 国家大地坐标系下的坐标，利用地理信息软件进行图形纠正，完成数据转换。2、按线性条带建立的图形数据库按线性条带（境界、河流、交通线、管道线等）建立的
15、图形数据库，可依据条带的方向、长短等分段进行，再拼接合成；也可通过条带中一定密度地物点的两套坐标，通过软件逐点进行纠正。具体方法：分块纠正：对于 1:1 万分块，按 1:1 万数据转换方案逐块纠正后接边合成；对于 1:5 万分块，按 1:2.5-1:10 万数据转换方案逐块纠正后接边合成； 逐点纠正：依据数据精度，建立一定密度（1:1 万 100 米21格网点、1:5 万 2000 米格网点）的坐标转换参考点，计算这些参考点在新坐标系下的坐标，利用地理信息软件完成数据转换。3、按无固定分幅分区建立的图形数据库按无固定分幅分区建立的图形数据，根据坐标系、比例尺及数据主体所在的图幅、数据的组织方式
16、、产品类型（DLG、 DEM、DRG ）等，参照相应比例尺的转换方案，实施数据转换。4、DOM 数据库转换原数据为航空或航天遥感获取的黑白或彩色影像数据，是连续的灰度（全色）或 RGB（彩色）栅格数据，分辨率有多种方式（主要包括用于 1:5 地形图测绘的各种分辨率航空影像，以及用于专题调查的 10 米、15 米、30 米等卫星影像） 。影像数据转换可参照下列方式进行。对于已按数据库组织方式加工与处理的 DOM 数据，可采用 1:2.5-1:10 万 DEM 的数据转换方法， 也可采用计算各景影像有效图边的 4 点在 2000 国家大地坐标系下的坐标来重新定位的方式。对于尚未按数据库组织方式加工
17、与处理的 DOM 数据，可采用 1:2.5-1:10 万 DRG 的数据转换方法，不再添加新的控制基础信息。22分辨率 5 米-30 米的数据，需依据其数据主体所在的 1:25万图幅区域来选用 1:25 万对应图幅的综合坐标改正值；对于分辨率在 2 米到 5 米间的数据，需依据其数据主体所在的 1:5 万图幅区域来选用 1:5 万对应图幅的综合坐标改正值；由此确定各自的 X、Y 方向平移像素数对应的坐标值（直接取 1:25 万或1:5 万综合坐标改正值，或由像素数 像素分辨率求得） 。按高斯投影、分像对（分景）组织的高分辨率影像数据，参照 1:1万 DOM 转换技术方案进行转换。四、1:1 万
18、及 1:5 千基础地理信息数据库的转换（一）1:1 万及 1:5 千格网点坐标转换改正量计算1、1980 西安坐标系坐标转换改正量计算1:1 万以上大比例尺一般按(23) 进行分区，并对每个分区向外扩充约 20，分别解算出各分区的转换参数后，利用确定的转换方法与转换模型分别计算全国 1:1 万及 1:5 千格网点的 2000 国家大地坐标系坐标 B2000，L 2000，进而求出各点的1980 西安坐标系与 2000 国家大地坐标系的差值DB802000， DL802000（ B2000-B80，L 2000-L80） ，形成全国 1:1 万及1:5 千格网点的 1980 西安坐标系与 200
19、0 国家大地坐标系的转换改正量 DB802000， DL802000。2、1954 年北京坐标系坐标转换改正量计算全国 1954 年北京坐标系向 2000 国家大地坐标系转换改正23量计算采用两步法：首先计算 1954 年北京坐标系转换向 1980西安坐标系转换改正量，其次计算 1980 西安坐标系向 2000 国家大地坐标系转换改正量，最后将两改正量叠加形成 1954 年北京坐标系向 2000 国家大地坐标系转换坐标转换改正量。1954 年北京坐标系向 1980 西安坐标系转换坐标改正量计算新旧坐标系统（1954 年北京坐标系与 1980 年西安坐标系）的转换大地坐标改正量计算公式： BeW
20、NaeMBzLyBMxd cosin)si2(1cosinscosin 22)i(1xNL式中: 分别为 IAG-75 椭球与克拉索夫斯基椭球长半径，2,ea第一偏心率平方之差。即 2280548054,ae则各个点在 1980 西安坐标系中的大地坐标为：dB5480 dL5480 根据转换的 ，采用高斯投影正算公式计算相应的高80,L斯平面坐标 。80,YX 求取全国 1:1 万以大比例尺格网点的转换改正量54801D 平差改正量的计算241954 年北京坐标系所提供的大地点成果没有经过整体平差，1980 西安坐标系提供的大地点成果是经过整体平差的数据，所以新旧系统转换还要考虑平差改正量的问
21、题。计算平差改正量比较麻烦，没有一定的数学模式，不同地区，平差改正量差别很大，在我国中部某些地区，平差改正量在 1 米以下，而在东北地区的某些图幅则在 10 米以上。在实际计算中，在全国均匀地选择一定数量的一、二等大地点，利用它们新（1980 西安坐标系）旧（1954 年北京坐标系）坐标系的坐标进行多种分析试算并剔除粗差点，然后分别计算它们的坐标差值，根据这些差值和它们的大地坐标分别绘制两张平差改正量分布图（即dX，dY 分布图） ，这样在分布图上可以直接内插出全国 1:1 万以大比例尺格网点的平差改正量 DX2，DY 2。 根据全国 1:1 万以大比例尺格网点的转换改正量DX1， DY1和平
22、差改正量 DX2，DY 2按下列公式计算 1954 年北京坐标系向 1980 西安坐标系转换坐标转换改正量 DX，DY 。21DYX 将 DX，DY 换算成 1:1 万以大比例尺格网点大地坐标转换改正量 DB5480，DL 5480。1954 年北京坐标系向 2000 国家大地坐标系转换坐标转换改正量计算25将全国 1:1 万以大比例尺格网点的 1954 年北京坐标系向1980 西安坐标系的转换改正量 DB5480与 1980 西安坐标系向2000 国家大地坐标系的转换改正量 DB802000叠加，得到全国1:1 万以大比例尺格网点 1954 年北京坐标系向 2000 国家大地坐标系转换的坐标
23、转换改正量 DB542000。即：DB 542000=DB5480+DB802000DL542000=DL5480+DL802000（二）1:1 万及 1:5 千 DLG 数据库转换转换流程如图 1 所示：26图 1 1:1 万及 1:5 千基础地理信息数据库转换技术流程1、1954 年北京坐标系下 1:1 万、1:5 千 DLG 数据库转换（1）每个图幅的四个图廓点坐标改正量选用 1954 年北京坐标系向 2000 国家大地坐标系转换方法计算；（2）图幅内各要素点的坐标改正量根据选用的本图幅的四个图廓点坐标改正量，按双线性内插等方法计算；27（3）根据图幅四个图廓点坐标改正量和图幅内各要素点
24、的坐标改正量，计算 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 年北京坐标系
25、到 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 数据；28（2）修改元数据相关条目；（3）修改相关的图外整饰。2、1980 西安坐标系下 1:1 万、1:5 千 DRG 数据库转换参照 1954 年北京坐标系到 2000 国家大地坐标系 1:1 万、1:5千 DRG
26、 数据库转换的对应步骤进行。（四）1:1 万及 1:5 千 DEM 数据库转换利用 DEM 生产过程中形成的矢量数据与 DEM 离散点数据完成数据转换。1、1954 年北京坐标系下 1:1 万、1:5 千 DEM 数据库转换（1）矢量数据与 DEM 离散点数据转换方法；a、每个图幅的四个图廓点坐标改正量选用 1954 年北京坐标系向 2000 国家大地坐标系转换方法计算；b、图幅内各要素点的坐标改正量根据选用的本图幅的四个图廓点坐标改正量，按双线性内插等方法计算；c、根据图幅四个图廓点坐标改正量和图幅内各要素点的坐标改正量，计算 2000 国家大地坐标系下的图幅四个图廓点坐标和图幅内各要素点的
27、坐标；d、与周边图幅拼接。（2）构造 TIN；（3）按相关规范或规定内插 DEM；（4）对检查点坐标文件进行点对点坐标转换；29（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
28、 上叠加相应坐标系统的内图廓及公里格网，在 2000 国家大地坐标系下生成图廓坐标及公里格网，逐公里格网点纠正 1:1 万、 1:5 千 DOM 数据；2、转换后，删除内图廓及公里格网；3、修改元数据相关条目。五、相对独立的平面坐标系统与 2000 国家大地坐标系建立联系的方法（一）相对独立的平面坐标系统控制点建立联系的方法可通过现行国家大地坐标系的平面坐标过渡，利用坐标转换方法将相对独立的平面坐标系统下控制点成果转换到 2000国家大地坐标系下。30选取相对独立的平面坐标系统与 2000 国家大地坐标系的重合点的原则如下：择优选取地方控制网的起算点及高精度控制点、周围国家高精度的控制点，大中
29、城市至少选取 5 个重合点（城外 4 个，市内中心 1 个） ；小城市在城市外围至少选取4 个重合点，重合点要分布均匀，包围城市区域，并在城市内部选定至少 6 个均匀分布的重合点对坐标转换精度进行检核。建立相对独立的平面坐标系统与 2000 国家大地坐标系联系时，坐标转换模型要同时适用于地方控制点转换和城市数字地图的转换。一般采用平面四参数转换模型，重合点较多时可采用多元逐步回归模型。当相对独立的平面坐标系统控制点和数字地图均为三维地心坐标时，采用 Bursa 七参数转换模型。坐标转换中误差应小于 0.05 米。（二）相对独立的平面坐标系统下数字地形图转换采用点对点转换法完成相对独立的平面坐标
30、系统下数字地形图到 2000 国家大地坐标系的转换，转换后相邻图幅不存在接边问题。具体步骤如下： 利用控制点的转换模型和参数，对相对独立的平面坐标系统下数字地形图进行转换，形成 2000 国家大地坐标系地形图。 根据转换后的图幅四个图廓点在 2000 国家大地坐标系下的坐标，重新划分公里格网线，原公里格网线删除。31 根据 2000 国家大地坐标系下的图廓坐标，对每幅图进行裁剪和补充。六、坐标转换方法（一）坐标转换模型1、二维七参数转换模型222sincos“0coincos“csi1ino0sinc“0(si)sic“1xyzLLXNBBYZMMtgLt mNeBNeBMa ios“aff
31、其中：同 一 点 位 在 两 个 坐 标 系 下 的 纬 度 差 、 经 度 差 ， 单 位 为,BL弧 度 ，椭球长半轴差（单位米） 、扁率差（无量纲） ，,af平移参数，单位为米，,XYZ旋转参数，单位为弧度，,xyz尺度参数（无量纲） 。m2、平面四参数转换模型属于两维坐标转换，对于三维坐标，需将坐标通过高斯投32影变换得到平面坐标再计算转换参数。平面直角坐标转换模型：102 cossini)1(yxmyx其中，x 0，y 0为平移参数， 为旋转参数，m 为尺度参数。x2，y 2为 2000 国家大地坐标系下的平面直角坐标，x 1，y 1为原坐标系下平面直角坐标。坐标单位为米。3、综合法
32、坐标转换所谓综合法即就是在相似变换（Bursa 七参数转换）的基础上，再对空间直角坐标残差进行多项式拟合，系统误差通过多项式系数得到消弱，使统一后的坐标系框架点坐标具有较好的一致性，从而提高坐标转换精度。综合法转换模型及转换方法： 利用重合点先用相似变换转换Bursa 七参数坐标转换模型 SSZYXSST ZYmYZXZY00式中，3 个平移参数 ，3 个旋转参数 和TTZX1 个尺度参数 。m 对相似变换后的重合点残差 采用多项式拟合ZYXV,KiijjSiZYXLBaV0或或式中：B，L 单位：弧度；K 为拟合阶数； 为系数，通过最ija33小二乘求解。4、三维七参数坐标转换模型 2 22
33、sincos“0()co()incos“()ssin(1)1)cosi1sinLLXLNHBNBYBMHMZeetgBLtgBLNNH 22 222222( 0cosicsins0inos“()00(sin)cos“sinco“1(1) xyzBeemMNHeeBfaeBa 22iaf 同一点位在两个坐标系下的纬度差、经度差、大地,BLH高差，经纬度差单位为弧度，大地高差单位为米，弧度秒18036/椭球长半轴差，单位为米， a扁率差，无量纲，f平移参数，单位为米，,XYZ旋转参数，单位为弧度，,xyz34尺度参数，无量纲。m（二）高斯正反算公式1、高斯投影正算公式 4422422 cos)58
34、61(70cos)95(41.0cos lBtlBtlBNtx 422422185(0cos)(61 lNttlBtly子午线弧长 X 计算见附录。2、高斯投影反算公式 42222 )(59061(3)(935(1)( ffffffffff NytNyttNyMtB 4224222 88510)6cos ffffffffff ttytl 式中 、 分别为按 值计算的相应量， 的计算见附录。tfBfB（三）坐标转换精度评定和评估方法 对于 1954 年北京坐标系、1980 西安坐标系与 2000 国家大地坐标系转换分区转换及数据库转换点位的平均精度应小于图上的 0.1mm。具体：对于 1:5 千
35、坐标转换，1980 西安坐标系与 2000 国家大地坐标系转换分区转换平均精度0.5m； 1954 年北京坐标系与2000 国家大地坐标系转换分区转换平均精度1.0m；1:5 万基础地理信息数据库坐标转换精度5.0m；351:1 万基础地理信息数据库坐标转换精度1.0m；1:5 千基础地理信息数据库坐标转换精度0.5m。依据计算坐标转换模型参数的重合点的残差中误差评估坐标转换精度。对于 n 个点，坐标转换精度估计公式如下： V（残差）=重合点转换坐标-重合点已知坐标空间直角坐标 X 残差中误差1XvMn空间直角坐标 Y 残差中误差YY空间直角坐标 Z 残差中误差1ZZvn点位中误差 22ZYX
36、pM平面坐标 x 残差中误差1xxvn平面坐标 y 残差中误差yy大地高 H 残差中误差1HvMn平面点位中误差为 2yxp36附录：1、常用量定义a 为椭球长半轴， 1954 年北京坐标系为 6378245m1980 西安坐标系为 6378140mb 为椭球短半轴 为椭球扁率 , 1954 年北京坐标系为 1/298.3f1980 西安坐标系为 1/298.257abf21ee第一偏心率ab222efe 第二偏心率b2Be22costgeV2cs1221VBWinB 为纬度，单位弧度 bac237M子午圈曲率半径 332)1(VcWeaMN卯酉圈曲率半径 2、子午线弧长 X 设有子午线上两点
37、 p1 和 p2，p1 在赤道上 p2 的纬度为B，p1、p2 间的子午线弧长 X 计算公式 )12sin0si8sin 6sin4)(12 BGFBEDCarcAeX式中1210128106 1284 1062 1284 83654937268361054572937283661340551 549727238 863643105765 eGeF eeE eD eeCeeB eeA 3、底点纬度 迭代公式fB)(,)1(120 iii BFXBAeaX直到 小于某一个指定数值，即可停止迭代。iiB138式中 12sin0si8sin 6sin4)1()2 BGFBEDCarcAeF cos12cos10cos864)()2eaB
