CGCS2000----2000中国大地坐标系

书 关于书的成语关于读书的排比句社区图书漂流公约怎么写关于读书的小报汉书pdf 书书 第２８卷第６期 ２００８年１２月 大 地 测 量 与 地 球 动 力 学 ＪＯＵＲＮＡＬＯＦＧＥＯＤＥＳＹＡＮＤＧＥＯＤＹＮＡＭＩＣＳ Ｖｏｌ．２８Ｎｏ．６ 　 Ｄｅｃ．，２００８ 　　文章编号：１６７１５９４２（２００８）０６０００１０５ ２０００中国大地坐标系 魏子卿 （西安测绘研究所，西安　７１００５４） 摘　要　２０００中国大地坐标系是我国新一代大地坐标系。在扼要叙述我国新一代大地坐标系建立的背景和基 本原则之后，详细介绍２０００中国大地坐标系的定义和实现，给出参考椭球的定义常数和导出常数，以及相关的正 常重力公式，最后对坐标系进行几点说明。 关键词　２０００中国大地坐标系（ＣＧＣＳ２０００）；地心大地坐标系；参考椭球；正常重力公式；常数 中图分类号：Ｐ２２７　　　　文献标识码：Ａ ＣＨＩＮＡＧＥＯＤＥＴＩＣＣＯＯＲＤＩＮＡＴＥＳＹＳＴＥＭ２０００ ＷｅｉＺｉｑｉｎｇ （Ｘｉ’ａｎＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅｏｆＳｕｒｖｅｙｉｎｇａｎｄＭａｐｐｉｎｇ，Ｘｉ’ａｎ　７１００５４） Ａｂｓｔｒａｃｔ　ＴｈｅＣｈｉｎａＧｅｏｄｅｔｉｃＣｏｏｒｄｉｎａｔｅＳｙｓｔｅｍ２０００ｉｓａｎｅｗｇｅｎｅｒａｔｉｏｎｏｆｔｈｅｎａｔｉｏｎａｌｇｅｏｄｅｔｉｃｃｏｏｒｄｉｎａｔｅ ｓｙｓｔｅｍ．Ｈａｖｉｎｇｅｘｐｌａｉｎｅｄｂｒｉｅｆｌｙｔｈｅｂａｃｋｇｒｏｕｎｄａｎｄｔｈｅｂａｓｉｃｐｒｉｎｃｉｐｌｅｓｏｆｓｅｔｔｉｎｇｕｐｔｈｅｎｅｗｇｅｎｅｒａｔｉｏｎｎａｔｉｏｎａｌ ｇｅｏｄｅｔｉｃｃｏｏｒｄｉｎａｔｅｓｙｓｔｅｍ，ｔｈｉｓｐａｐｅｒｅｘｐｏｕｎｄｓｉｔｓｄｅｆｉｎｉｔｉｏｎａｎｄｒｅａｌｉｚａｔｉｏｎ，ｇｉｖｅｓｔｈｅｄｅｆｉｎｉｎｇｃｏｎｓｔａｎｔｓａｎｄｄｅ ｒｉｖｅｄｃｏｎｓｔａｎｔｓｏｆｔｈｅｒｅｆｅｒｅｎｃｅｅｌｌｉｐｓｏｉｄｕｓｅｄａｎｄｔｈｅｒｅｌａｔｅｄｆｏｒｍｕｌａｓｆｏｒｔｈｅｎｏｒｍａｌｇｒａｖｉｔｙ，ａｎｄｆｉｎａｌｌｙａｄｄｒｅｓｓｅｓ ｓｅｖｅｒａｌｉｔｅｍｓｎｏｔｉｎｇｏｎｔｈｅｃｏｏｒｄｉｎａｔｅｓｙｓｔｅｍ． Ｋｅｙｗｏｒｄｓ：ＣｈｉｎａＧｅｏｄｅｔｉｃＣｏｏｒｄｉｎａｔｅＳｙｓｔｅｍ２０００；ｇｅｏｃｅｎｔｒｉｃｇｅｏｄｅｔｉｃｃｏｏｒｄｉｎａｔｅｓｙｓｔｅｍ；ｒｅｆｅｒｅｎｃｅｅｌｌｉｐｓｏｉｄ； ｎｏｒｍａｌｇｒａｖｉｔｙｆｏｒｍｕｌａｓ；ｃｏｎｓｔａｎｔ １　引言 如所周知，２０世纪５０年代，为满足测绘工作的 迫切需要，我国采用了１９５４年北京坐标系，后来随 着天文大地网布设任务的完成，通过天文大地网整 体平差，于８０年代初我国又建立了１９８０西安坐标 系。应当肯定，１９５４北京坐标系和１９８０西安坐标 系在我国的经济建设和国防建设中发挥了巨大作 用。同时也应当看到，随着情况的变化和时间的推 移，这两个以经典测量技术为基础的局部大地坐标 系，目前已经不能适应科学技术特别是空间技术发 展，不能适应我国经济建设和国防建设需要。我国 大地坐标系的更新换代，是经济建设、国防建设、社 会发展和科技发展的客观需要。 以地球质量中心为原点的地心大地坐标系，是 当今空间时代全球通用的基本大地坐标系。以空间 技术为基础的地心大地坐标系，是我国新一代大地 坐标系的适宜选择。地心大地坐标系可以满足大地 测量、地球物理、天文、导航和航天应用以及经济、社 会发展的广泛需求。多年以来，我国测绘、地震部门 和科学院有关单位为建立我国新一代大地坐标系作 了大量基础性工作，近年又先后建成全国 ＧＰＳ一、 二级网，国家ＧＰＳＡ、Ｂ级网，中国地壳运动观测网 络和许多地壳形变网，为地心大地坐标系的实现奠  收稿日期：２００８０９３０ 作者简介：魏子卿，男，１９３７年生，中国工程院院士，近期主要研究大地坐标系和大地边值问题．Ｅ－ｍａｉｌ：ｚｉｑｉｎｇｗ＠ｐｕｂｌｉｃ．ｘａ．ｓｎ．ｃｎ 大地测量与地球动力学 ２８卷 定了较好的基础。我国大地坐标系更新换代的条件 业已具备。 我国新一代大地坐标系建立的基本原则是： １）坐标系应尽可能对准 ＩＴＲＦ（国际地球参考 架）［１］； ２）坐标系应由空间大地网在某参考历元的坐 标和速度体现； ３）参考椭球的定义参数选用长半轴、扁率、地 球地心引力常数和地球角速度，其参数值采用ＩＵＧＧ （国际大地测量与地球物理联合会）或 ＩＥＲＳ（国际 地球旋转与参考系服务局）的采用值或推荐值。 ２０００中国大地坐标系（ＣｈｉｎａＧｅｏｄｅｔｉｃＣｏｏｒｄｉ ｎａｔｅＳｙｓｔｅｍ２０００，ＣＧＣＳ２０００），国人又称之为 ２０００ 国家大地坐标系，是我国新一代大地坐标系，现已在 全国正式实施［２］。本文意在介绍２０００中国大地坐 标系的定义和实现，给出参考椭球的定义常数和导 出常数以及相关的正常重力公式，并对坐标系作几 点说明。 ２　２０００中国大地坐标系 ２．１　定义 ２０００中国大地坐标系符合 ＩＴＲＳ（国际地球参 考系）的如下定义［１］： １）原点在包括海洋和大气的整个地球的质量 中心； ２）长度单位为米（ＳＩ）。这一尺度同地心局部 框架的ＴＣＧ（地心坐标时）时间坐标一致； ３）定向在１９８４．０时与 ＢＩＨ（国际时间局）的定 向一致； ４）定向随时间的演变由整个地球的水平构造 运动无净旋转条件保证。 以上定义对应一个直角坐标系，它的原点和轴 定义如下（图１）： １）原点：地球的质量中心； ２）Ｚ轴：指向ＩＥＲＳ参考极方向； 　　３）Ｘ轴：ＩＥＲＳ参考子午面与通过原点且同 Ｚ 轴正交的赤道面的交线； ４）Ｙ轴：完成右手地心地固直角坐标系。 图１　ＣＧＣＳ２０００定义的示意图 Ｆｉｇ．１　ＳｋｅｔｃｈｆｏｒｔｈｅｄｅｆｉｎｉｔｉｏｎｏｆＣＧＣＳ２０００ ＣＧＣＳ２０００的参考椭球为一等位旋转椭球。等 位椭球（或水准椭球）定义为其椭球面是一等位面 的椭球［３］。ＣＧＣＳ２０００的参考椭球的几何中心与坐 标系的原点重合，旋转轴与坐标系的 Ｚ轴一致。参 考椭球既是几何应用的参考面，又是地球 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 面上及 空间正常重力场的参考面。等位旋转椭球由４个独 立常数定义。ＣＧＣＳ２０００参考椭球的定义常数 是［２］： 长半轴ａ＝６３７８１３７．０ｍ； 扁率ｆ＝１／２９８．２５７２２２１０１； 地球的地心引力常数（包含大气层）ＧＭ＝ ３９８６００４．４１８×１０８ｍ３ｓ－２； 地球角速度ω＝７２９２１１５．０×１０－１１ｒａｄｓ－１。 根据以上４个定义常数，利用文献［３，４］给出 的公式，可以算出一系列导出常数。常用的导出几 何常数列于表１，导出物理常数列于表２。 ２．２　实现 ＣＧＣＳ２０００由２０００国家 ＧＰＳ大地网［５］在历元 ２０００．０的点位坐标和速度具体实现。２０００国家 ＧＰＳ大地网是下列 ＧＰＳ网经联合平差得到的一个 全国规模的ＧＰＳ大地网（图２）： 表１　ＣＧＣＳ２０００参考椭球的导出几何常数值 Ｔａｂ．１　ＤｅｒｉｖｅｄｇｅｏｍｅｔｒｉｃａｌｃｏｎｓｔａｎｔｓｆｏｒｔｈｅＣＧＣＳ２０００ｅｌｌｉｐｓｏｉｄ 常数名 值 常数名 值 短半轴 ｂ ６３５６７５２．３１４１ｍ 轴比 ｂ／ａ ０．９９６６４７１８９３１９ 线偏心率Ｅ ５２１８５４．００９７００２５ｍ 子午圈一象限弧长 Ｑ １０００１９６５．７２９３ｍ 极曲率半径 ｃ ６３９９５９３．６２５９ｍ 椭球体积 Ｖ １０８３２０７３１９７８３．５４６ｋｍ３ 第一偏心率平方 ｅ２ ０．００６６９４３８００２２９０ 椭球表面积 Ｓ ５１００６５６２１．７１８ｋｍ２ 第一偏心率 ｅ ０．０８１８１９１９１０４２８１６ 平均半径 Ｒ１ ６３７１００８．７７１４ｍ 第二偏心率平方 ｅ′２ ０．００６７３９４９６７７５４８ 同面积之球的半径Ｒ２ ６３７１００７．１８０９ｍ 第二偏心率 ｅ′ ０．０８２０９４４３８１５１９１７ 同体积之球的半径Ｒ３ ６３７１０００．７９００ｍ ２ 　第６期 魏子卿：２０００中国大地坐标系 表２　ＣＧＣＳ２０００参考椭球的导出物理常数值 Ｔａｂ．２　ＤｅｒｉｖｅｄｐｈｙｓｉｃａｌｃｏｎｓｔａｎｔｓｆｏｒｔｈｅＣＧＣＳ２０００ｅｌｌｉｐｓｏｉｄ 常数名 值 常数名 值 椭球面正常位 Ｕ０ ６２６３６８５１．７１４９ｍ２ｓ－２ 赤道正常重力 γｅ ９．７８０３２５３３６１ｍｓ－２ ２阶带谐系数 Ｊ２ ０．１０８２６２９８３２２５８×１０－２ 极正常重力 γｐ ９．８３２１８４９３７９ｍｓ－２ ４阶带谐系数 Ｊ４ －０．２３７０９１１２５６１４１×１０－５ 平均正常重力珔γ ９．７９７６４３２２２４ｍｓ－２ ６阶带谐系数 Ｊ６ ０．６０８３４６５２５８８９２×１０－８ ｍ ＝ω２ａ２ｂ／（ＧＭ） ０．００３４４９７８６５０６７８ ８阶带谐系数 Ｊ８ －０．１４２６８１１００９７９８×１０－１０ ｋ＝ｂγｐ／ａγｅ－１ ０．００１９３１８５２６１９３１ １０阶带谐系数 Ｊ１０ ０．１２１４３９３３８３３４３×１０－１３ 地球质量Ｍ（包括大气） ５．９７３３３１９６×１０２４ｋｇ 图２　２０００国家ＧＰＳ大地网的点位分布图 Ｆｉｇ．２　ＤｉｓｔｒｉｂｕｔｉｏｎｏｆｐｏｉｎｔｓｏｆｔｈｅｎａｔｉｏｎａｌＧＰＳｇｅｏｄｅｔｉｃｎｅｔｗｏｒｋ２０００ 　　１）全国 ＧＰＳ一、二级网（５３４点），１９９１—１９９７ 年布设； ２）国家 ＧＰＳＡ、Ｂ级网（８１８点），１９９１—１９９６ 年布设； ３）地壳运动监测网（３３６点），包括全国ＧＰＳ地 壳运动监测网（２１点），１９９４、１９９６和１９９９年进行３ 次观测；还包括９个区域性地壳形变监测网，１９８８— １９９８年施测； ４）中国地壳运动观测网络（１０８１点），包括基 准网（２５点）、基本网（５６点）和区域网（１０００点）。 基准网自１９９９年１月连续观测；基本网和区域网于 １９９９、２０００和２００１年进行３次施测。基本网每次 连续观测１０昼夜，区域网连续观测４昼夜。平差数 据截止于２００１年底。 联合平差从中国地壳运动观测网络平差开始。 它的平差分两步进行。第一步进行单时段（一般为 ２４ｈ）基线解。第二步将得到的所有单时段基线解 与全球１００多个 ＩＧＳ站的 Ｈ文件数据进行整体平 差，此时将全球分布的４７个 ＩＧＳ核心站的 ＩＴＲＦ９７ 坐标和速度施以１σ的约束，得到该网络点位三维 坐标和速度的最小二乘解及协方差矩阵。坐标分量 的平均中误差为 σｘ＝０．０７０ｃｍ，σｙ＝０．１１０ｃｍ，σｚ ＝０．０９０ｃｍ；速度分量误差为１～４ｍｍ／ａ。最后用 得到的速度将平差坐标归算到历元２０００．０，作为随 ３ 大地测量与地球动力学 ２８卷 后诸网联合平差的基准数据。对于其余 ３个 ＧＰＳ 网，分别同样首先进行单时段基线解，然后将它们的 单时段基线解与先前得到的中国地壳运动观测网络 的站坐标解一起进行联合平差，平差中固定ＩＧＳ站、 基准网和基本网的站坐标，对区域网站的南北、东西 和高度坐标分量分别给予５ｍｍ、５ｍｍ和１０ｍｍ的 约束。通过平差得到２０００国家 ＧＰＳ大地网的最后 坐标。主要精度统计数据是： 坐标平均中误差：σｘ＝０．８４ｃｍ，σｙ＝１．８２ｃｍ， σｚ＝１．３０ｃｍ，σＢ＝０．４０ｃｍ，σＬ＝０．５２ｃｍ，σｈ＝２．３１ ｃｍ； 位置平均中误差：σｐ＝２．４２ｃｍ； 基线长度（不计２０ｋｍ以内的基线，平均长度 为１０６ｋｍ）平均误差：０．０３０１×１０－６。 以上平差数据显示，ＣＧＣＳ２０００的实现精度，在 历元２０００．０时，位置精度好于３ｃｍ；水平速度约３ ｍｍ／ａ。 ３　ＣＧＣＳ２０００正常重力公式 ３．１　椭球面上正常重力公式 椭球面上正常重力由Ｓｏｍｉｇｌｉａｎａ公式［３，６］给出： γ０＝γｅ １＋ｋｓｉｎ２Ｂ １－ｅ２ｓｉｎ２槡 Ｂ （１） 式中；ｋ＝ｂγｐ／ａγｅ－１；γｅ、γｐ为赤道重力和极重力； ａ、ｂ为椭球长半轴和短半轴；ｅ２为第一偏心率（平 方），Ｂ为大地纬度。γｅ、γｐ和ｋ之值见表２；ｂ、ｅ ２之 值见表１。 实践中，通常将封闭式（１）展成级数［３，４］。对于 ＣＧＣＳ２０００，级数展开式是（单位：ｍｓ－２）： γ０＝９．７８０３２５３３６１（１＋０．００５２７９０４２６３１ｓｉｎ ２Ｂ＋ ０．００００２３２７１７９９ｓｉｎ４Ｂ＋ ０．００００００１２６２１８ｓｉｎ６Ｂ＋ ０．０００００００００７３０ｓｉｎ８Ｂ＋ ０．０００００００００００４ｓｉｎ１０Ｂ） （２） 该式的误差为０．００１×１０－８ｍｓ－２。传统的简化 级数式，对于ＣＧＣＳ２０００是（单位：ｍｓ－２）： γ０＝９．７８０３２５３３６１（１＋０．００５３０２４４ｓｉｎ ２Ｂ－ ０．０００００５８２ｓｉｎ２２Ｂ） （３） 该式的误差小于０．１×１０－５ｍｓ－２。 ３．２　椭球外部正常重力公式 椭球外部任意点的正常重力，由以下公式给 出［４，７］： γ＝ γ２ｕ＋γ ２槡 β （４） 式中： γｕ＝－ １ ｗ ＧＭ ｕ２＋Ｅ２ ＋ω ２ａ２Ｅ ｕ２＋Ｅ２ ｑ′ ｑ０ １ ２ｓｉｎ ２β－( )[ ]１６ ＋ １ ｗω ２ｕｃｏｓ２β （４ａ） γβ＝－ １ ｗ － ω ２ａ２ ｕ２＋Ｅ槡 ２ ｑ ｑ０ ＋ω２ ｕ２＋Ｅ槡( )２ ｓｉｎβｃｏｓβ （４ｂ） ｗ＝ ｕ ２＋Ｅ２ｓｉｎ２β ｕ２＋Ｅ槡 ２ （４ｃ） ｕ＝ １２ ｘ ２＋ｙ２＋ｚ２－Ｅ２[{ ＋ （ｘ２＋ｙ２＋ｚ２－Ｅ２）２＋４Ｅ２ｚ槡 ] }２ １／２ （４ｄ） β＝ｔａｎ－１ ｚ ｕ ２＋Ｅ槡 ２ ｕ ｘ２＋ｙ槡( )２ （４ｅ） ｑ′＝３１＋ｕ ２ Ｅ( )２ １－ｕＥｔａｎ－１ Ｅ( )[ ]ｕ －１ （４ｆ） ｑ＝１２ １＋３ｕ ２ Ｅ( )２ ｔａｎ－１ Ｅ( )ｕ －３ｕ[ ]Ｅ （４ｇ） ｑ０＝ １ ２ １＋３ｂ ２ Ｅ( )２ ｔａｎ－１ Ｅ( )ｂ －３ｂ[ ]Ｅ （４ｈ） 式中ｘ、ｙ、ｚ和 ｕ、β、λ分别代表计算点的地固直角坐 标和椭球坐标系坐标；ｂ为参考椭球短半轴，Ｅ为线 偏心率，它们的值见表１。 封闭公式（４）可按高度 ｈ展成级数［６，７］。对于 ＣＧＣＳ２０００，至ｈ的４次项的级数式是： γ＝γ０－（３．０８３３８７８８８７１×１０ －６＋４．４２９７４３９６３ ×１０－９ｃｏｓ２Ｂ－１．９９６４６１４×１０－１１ｃｏｓ４Ｂ）ｈ＋ （７．２４４２７７７９９９×１０－１３＋ ２．１１６０６２×１０－１５ｃｏｓ２Ｂ－３．３４３０６×１０－１７ｃｏｓ４Ｂ－ １．９０８×１０－１９ｃｏｓ６Ｂ－４．８６×１０－２２ｃｏｓ８Ｂ）ｈ２－ （１．５１１２４９２２×１０－１９＋１．１４８６２４×１０－２１ｃｏｓ２Ｂ＋ １．４９７５×１０－２３ｃｏｓ４Ｂ＋１．６６×１０－２５ｃｏｓ６Ｂ）ｈ３＋ （２．９５２３９×１０－２６＋４．１６７×１０－２８ｃｏｓ２Ｂ）ｈ４ （５） 式中γ０为椭球面正常重力（单位为ｍｓ －２），按式（２） 计算；Ｂ为计算点的大地纬度，ｈ为它的大地高（单 位为ｍ）。 用式（５）计算正常重力的误差，当ｈ达到２０ｋｍ 时，小于０．１×１０－８ｍｓ－２；当ｈ达到７０ｋｍ时，小于１ ×１０－８ｍｓ－２。 ４　关于ＣＧＣＳ２０００的几点说明 关于ＣＧＣＳ２０００，扼要作以下几点说明： １）ＣＧＣＳ２０００的定义和 ＩＴＲＳ的定义一致。 ＣＧＣＳ２０００实现的实质是使 ＣＧＣＳ２０００框架对准 ＩＴＲＦ９７。相对ＩＴＲＦ９７，ＣＧＣＳ２０００的实现精度，对于 水平坐标达到１ｃｍ量级（见第２节）。因此可以认 为ＣＧＣＳ２０００与ＩＴＲＦ９７（ＩＴＲＦ２０００和ＩＴＲＦ２００５）在 ｃｍ级 水 平 上 是 一 致 的。但 是，倘 若 一 点 的 ４ 　第６期 魏子卿：２０００中国大地坐标系 ＣＧＣＳ２０００坐标精度达不到ｃｍ级水平，这两个坐标 系的坐标不可认为是一致的。 ２）ＣＧＣＳ２０００的定义和 ＷＧＳ８４［８，９］实质上是一 致的。它们采用的参考椭球非常接近，仅扁率存在 微小差异（１／２９８．２５７２２２１０１对 １／２９８．２５７２２３ ５６３［８，９］）。扁率差异引起椭球面上的纬度和高度坐 标变化最大达 ０．１０５ｍｍ，正常重力变化最大达 ０．０１６×１０－８ｍｓ－２［８，９］。在当前的测量精度范围内， 忽略这样小的坐标和重力变化是容许的。ＷＧＳ８４ （Ｇ１１１５０）［１０］是 ＷＧＳ８４的最近（２００２年）一次实 现。实现精度（相对ＩＴＲＦ２０００），对于每一坐标分量 为１ｃｍ量级［１０］，与 ＣＧＣＳ２０００的实现精度（相对 ＩＴＲＦ９７）大体相当。因此可以说，ＣＧＣＳ２０００与 ＷＧＳ８４（Ｇ１１５０）相容至ｃｍ级水平。但是，倘若一点 的坐标精度达不到ｃｍ水平，则不可认为 ＣＧＣＳ２０００ 和ＷＧＳ８４的坐标是相容的。 ３）ＣＧＣＳ２０００和１９５４年北京坐标系或１９８０西 安坐标系，在定义与实现上有根本区别。局部坐标 和地心坐标之间的变换是不可避免的。坐标变换可 以通过联合平差实现，而一般通过一定变换模型实 现。当采用模型变换时，变换模型的选择应依精度 要求而定。对于高精度（好于０．５ｍ）要求，可采用 基于最小曲率法［１１］或其他方法的格网模型；对于中 等精度（０．５～５ｍ）要求，可采用七参数（３个平移、１ 个尺度和３个旋转）模型；对于低精度（５～１０ｍ）要 求，可采用四参数（３个平移和１个尺度）或三参数 （３个平移）模型。 ４）ＣＧＣＳ２０００的实现为２０００国家 ＧＰＳ大地网 在历元 ２０００．０的点位坐标和速度，即是说， ＣＧＣＳ２０００框架由大约２５００个 ＧＰＳ点在历元 ｔ０＝ ２０００．０的坐标Ｘ → （ｔ０）和速度Ｘ →· （ｔ０）构成。这些点在 历元ｔ（ｔ≠ｔ０）的坐标Ｘ → （ｔ）按下式计算： Ｘ → （ｔ）＝Ｘ → （ｔ０）＋（ｔ－ｔ０）Ｘ →· （ｔ０） （６） 如果速度Ｘ →· （ｔ０）未知，还应根据周围已知点的速 度内插得到，或按现成的板块运动模型计算。这方 面的细节已超出本文的范围。 地面上一点的坐标，受地壳运动和潮汐影响，在 时间上和空间上都是变化的。数据表明，中国大陆 不同点在ＩＴＲＦ框架内的运动速度为３～５ｃｍ［９］，运 动方向大致从西北向东南。在 ｃｍ级精度要求的应 用（如ＧＰＳ网平差）中，记住要将已知点的坐标从历 元ｔ０归算至观测历元ｔ。 致谢　ＣＧＣＳ２０００椭球的导出常数和椭球面正常重 力公式的系数由李迎春副研究员进行了认真检核。 参 考 文 献 １　ＤｅｎｎｉｓＤ．ＭｃＣａｒｔｈｙａｎｄＧéｒａｒａＰｅｔｉｔ（ｅｄｓ．）．ＩＥＲＳＣｏｎ ｖｅｎｔｉｏｎｓ（２００３）［Ｒ］．ＩＥＲＳＴｅｃｈｎｉｃａｌＮｏｔｅＮｏ．３２，Ｖｅｒｌａｇ ｄｅｓＢｕｎｄｅｓａｍｔｓｆüｒＫａｒｔｏｇｒａｐｈｉｅｕｎｄＧｅｏｄｓｉｅＦｒａｎｋｆｕｒｔａｍ Ｍａｉｎ２００４． ２　国家测绘局．国家测绘局公告（２００８年第２号）．２００８年 ６月２７日中国测绘报． ２　ＳｔａｔｅＢｕｒｅａｕｏｆＳｕｒｖｅｙｉｎｇａｎｄＭａｐｐｉｎｇ．Ａｎｎｏｕｎｃｅｍｅｎｔｏｆ ｔｈｅＳｔａｔｅＢｕｒｅａｕｏｆＳｕｒｖｅｙｉｎｇａｎｄＭａｐｐｉｎｇ（Ｎｏ．２ｏｆ２００８）． ＣｈｉｎａＳｕｒｖｅｙｉｎｇａｎｄＭａｐｐｉｎｇＮｅｗｓ，２７Ｊｕｎｅ２００８． ３　ＭｏｒｉｔｚＨ．Ｇｅｏｄｅｔｉｃｒｅｆｅｒｅｎｃｅｓｙｓｔｅｍ１９８０［Ｊ］．Ｊｏｕｒｎａｌｏｆ Ｇｅｏｄｅｓｙ，２０００，７４：１２８－１３３． ４　魏子卿．关于２０００中国大地坐标系的建议［Ｊ］．大地测 量与地球动力学，２００６，２６（２）：１－４． ４　ＷｅｉＺｉｑｉｎｇ．ＰｒｏｐｏｓａｌｃｏｎｃｅｒｎｉｎｇＣｈｉｎａｇｅｏｄｅｔｉｃｃｏｏｒｄｉｎａｔｅ ｓｙｓｔｅｍ２０００［Ｊ］．ＪｏｕｒｎａｌｏｆＧｅｏｄｅｓｙａｎｄＧｅｏｄｙｎａｍｉｃｓ， ２００６，２６（２）：１－４．（ｉｎＣｈｉｎｅｓｅ） ５　陈俊，等．２０００国家大地控制网的构建和它的技术进步 ［Ｊ］．测绘学报，２００７，３６（１）：１－８． ５　ＣｈｅｎＪｕｎｙｏｎｇ，ｅｔａｌ．Ｅｓｔａｂｌｉｓｈｍｅｎｔｏｆ２０００ｎａｔｉｏｎａｌｇｅｏｄｅｔｉｃ ｃｏｎｔｒｏｌｎｅｔｗｏｒｋｏｆＣｈｉｎａａｎｄｉｔ’ｓｔｅｃｈｎｏｌｏｇｉｃａｌｐｒｏｇｒｅｓｓ［Ｊ］． ＡｃｔａＧｅｏｄａｅｔｉｃａｅｔＣａｒｔｏｇｒａｐｈｉｃａＳｉｎｉｃａ，２００７，３６（１）：１－ ８．（ｉｎＣｈｉｎｅｓｅ） ６　ＨｏｆｍａｎｎＷｅｌｌｅｎｈｏｆＢａｎｄＭｏｒｉｔｚＨ．Ｐｈｙｓｉｃａｌｇｅｏｄｅｓｙ［Ｍ］． ｓｅｃｏｎｄｅｄｉｔｉｏｎ，ＳｐｒｉｎｇｅｒＷｉｅｎＮｅｗＹｏｒｋ，２００６． ７　魏子卿．正常重力公式［Ｊ］．测绘学报，２００３，３２（２）：９５ －１０１． ７　ＷｅｉＺｉｑｉｎｇ．Ｎｏｒｍａｌｇｒａｖｉｔｙｆｏｒｍｕｌａｅ［Ｊ］．ＡｃｔａＧｅｏｄａｅｔｉｃａ ｅｔＣａｒｔｏｇｒａｐｈｉｃａＳｉｎｉｃａ，２００３，３２（２）：９５－１０１．（ｉｎＣｈｉ ｎｅｓｅ） ８　ＮａｔｉｏｎａｌＩｍａｇｅｒｙａｎｄＭａｐｐｉｎｇＡｇｅｎｃｙ．Ｄｅｐａｒｔｍｅｎｔｏｆｄｅ ｆｅｎｓｅｗｏｒｌｄｇｅｏｄｅｔｉｃｓｙｓｔｅｍ１９８４：Ｉｔｓｄｅｆｉｎｉｔｉｏｎａｎｄｒｅｌａｔｉｏｎ ｓｈｉｐｓｗｉｔｈｌｏｃａｌｇｅｏｄｅｔｉｃｓｙｓｔｅｍｓ［Ｒ］．ＮａｔｉｏｎａｌＩｍａｇｅｒｙａｎｄ ＭａｐｐｉｎｇＡｇｅｎｃｙｔｅｃｈｎｉｃａｌｒｅｐｏｒｔ８３５０．２ＴｈｉｒｄＥｄｉｔｉｏｎ． ９　魏子卿．２０００中国大地坐标系及其与 ＷＧＳ８４的比较 ［Ｊ］．大地测量与地球动力学，２００８，２８（５）：１－５． ９　ＷｅｉＺｉｑｉｎｇ．Ｃｈｉｎａｇｅｏｄｅｔｉｃｃｏｏｒｄｉｎａｔｅｓｙｓｔｅｍ２０００ａｎｄｉｔｓ ｃｏｍｐａｒｉｓｏｎｗｉｔｈＷＧＳ８４［Ｊ］．ＪｏｕｒｎａｌｏｆＧｅｏｄｅｓｙａｎｄＧｅｏｄｙ ｎａｍｉｃｓ，２００８，２８（５）：１－５．（ｉｎＣｈｉｎｅｓｅ） １０　ＡｄｄｅｎｄｕｍｔｏＮＩＭＡＴＲ８３５０．２：Ｉｍｐｌｅｍｅｎｔａｔｉｏｎｏｆｔｈｅ ＷｏｒｌｄＧｅｏｄｅｔｉｃＳｙｓｔｅｍ１９８４（ＷＧＳ８４）ＲｅｆｅｒｅｎｃｅＦｒａｍｅ Ｇ１１５０［ＥＢ／ＯＬ］．ｈｔｔｐ：／／ｅａｒｔｈ－ｉｎｆｏ．ｎｇａ．ｍｉｌ／ＧａｎｄＧ／ ｐｕｂｌｉｃａｔｉｏｎｓ／ｔｒ８３５０．２／Ａｄｄｅｎｄｕｍ％２０ＮＩＭＡ％２０ＴＲ８３５０． ２．ｐｄｆ／ １１　ＢｒｉｇｇｓＩＣ．Ｍａｃｈｉｎｅｃｏｕｔｏｕｒｉｎｇｕｓｉｎｇｍｉｎｉｍｕｍｃｕｒｖａｔｕｒｅ ［Ｊ］．Ｇｅｏｐｈｙｓｉｃｓ，１９７４，３９（１）：３９－４８． ５