1、大地测量控制点坐标转换技术规程国家测绘地理信息局2013 年 6 月目 录1 范围.752 规范性引用文件.753 术语和定义.754 控制点类型.794.1 控制点分类.794.2 控制点精度.804.3 控制点用途.805 控制点坐标转换模型及适用范围.815.1 总则.815.2 转换模型.815.3 适用范围.825.4 模型选取.836 控制点坐标转换.836.1 省级GNSS连续运行基准站点坐标归算.836.2 省市级卫星大地控制网点坐标转换.856.3 其他控制点坐标转换.857 精度评定.867.1 省级GNSS连续运行基准站点坐标归算精度评定.867.2 其他控制点坐标转换精
2、度评定.86附录A87附录B89附录C95附录D98附录E100大地测量控制点坐标转换技术规程1范围本规程规定了各种坐标系控制点坐标转换到2000国家大地坐标系时控制点选取、坐标转换模型、转换方法、精度评价等方面应遵循的原则、适用范围和精度要求。本规程适用于大地测量控制点坐标转换过程中的重合点选取、模型选择和转换方法。2规范性引用文件下列文件对于本文件的应用是必不可少的。凡是注日期的引用文件，仅注日期的版本适用于本文件。凡是不注日期的引用文件，其最新版本（包括所有的修改单）适用于本文件。GB/T171592009 大地测量术语GB/T183142009 全球定位系统（GPS）测量规范GB/T1
3、93912003 全球定位系统（GPS）术语及定义GB/T220212008 国家大地测量基本技术规定GB/T 285882012 全球导航卫星系统连续运行基准站网技术规范CH/T10042005 测绘技术设计规定2010（第三版） 测绘学名词3术语和定义下列术语和定义适用于本文件。3.1 控制点control point以一定精度测定其位置为其他测绘工作提供依据的固定点。3.2 地球椭球earth ellipsoid近似表示地球的形状和大小，并且其表面为等位面的旋转椭球。3.3 参考椭球reference ellipsoid最符合一定区域的大地水准面，具有一定大小和定位参数的旋转地球椭球。3
4、.4 国际地球参考系统 International TerrestrialReference System，ITRS由国际地球自转服务（IERS）给出的地球坐标系统的定义和大地测量常数。3.5 国际地球参考框架International Terrestrial ReferenceFrame,ITRF国际地球参考系统（ITRS）的实现。由国际地球自转服务局（IERS）根据空间大地测量技术，包括甚长基线干涉测量（VLBI）、卫星激光测距（SLR）、多里斯系统（DORIS）、全球定位系统（GPS）等，所确定的地面点的坐标所构成的集合。3.6 大地坐标系geodetic coordinate syst
5、em以地球椭球中心为原点、起始子午面和赤道面为基准面的地球坐标系。3.7 1954 年北京坐标系 Beijing Geodetic CoordinateSystem 1954将我国大地控制网与苏联1942普尔科沃大地坐标系联结后建立的我国过渡性大地坐标系。3.8 1980西安坐标系Xian Geodetic Coordinate System1980采用IAG 1975国际椭球，以JYD1968.0系统为椭球定向基准，大地原点设在陕西泾阳县永乐镇，采用多点定位所建立的大地坐标系。3.9 地心坐标系geocentric coordinate system原点O与地球质心重合，Z轴指向地球北极，X
6、轴指向格林尼治平子午面与地球赤道的交点E，Y轴垂直于XOZ的平面构成右手坐标系。3.10 1984 世界大地坐标系 World Geodetic System1984,WGS84美国军用大地坐标系统，坐标系定义和国际地球参考系统（ITRS）一致。3.11 2000国家大地坐标系China Geodetic CoordinateSystem 2000,CGCS2000原点为包括海洋和大气的整个地球的质量中心。Z轴为国际地球自转局（IERS）定义的参考极方向，X轴为参考子午面与垂直于Z轴的赤道面的交线，Y轴与Z轴和X轴构成右手正交坐标系。3.12 坐标转换coordinate transforma
7、tion包含坐标系变换和椭球基准变换两层含义。在测量数据处理过程中，采用适用的转换模型和转换方法，空间点从某一参考椭球基准下的坐标转换到另一坐标系统下的坐标。坐标转换过程就是转换参数的求解过程。3.13 坐标系变换coordinate conversion同一椭球基准下，空间点的不同坐标表示形式间进行变换。包括大地坐标系和空间直角坐标系的互相转换、空间直角坐标系与站心坐标系间的转换和高斯投影坐标正反算。3.14 椭球变换ellipsoid conversion空间点在不同的参考椭球间的坐标变换。3.15 重合点superposition point同时拥有不同坐标系坐标的控制点。3.16 平移
8、参数translation parameters两坐标系转换时，新坐标系原点在原坐标系中的坐标分量。3.17 旋转参数rotation parameters两坐标系转换时，把原坐标系中的各坐标轴左旋转到与新坐标系相应的坐标轴重合或平行时坐标系各轴依次转过的角度。3.18 尺度参数scale parameter两坐标系转换时引入的两坐标系中长度变化参数。3.19 连续运行基准站Continuously Operating ReferenceStation，CORS连续接收和发送本站坐标及其变化，GNSS星历，星钟差等信息的地面固定站。3.20 国际GNSS服务International GNSS
9、 Service，IGS提供全球导航卫星系统，包括 GPS、GLONASS、GALILEO等卫星星历，卫星钟差以及相应卫星系统的地面基准站坐标等方面信息的国际组织。3.21 参考历元reference epoch在天文学上，历元是为指定天球坐标或轨道参数而规定的某一特定时刻。在天文学和卫星定位中，所获数据对应的时刻也称为历元。3.22 观测历元observation epoch为比较不同时刻的观测结果，需要注明观测资料所对应的观测时刻，称之为观测历元。3.23 板块tectonic plate板块是指岩石圈板块，即地壳和软流圈以上的地幔顶部。3.24 2000国家GPS大地控制网Nationa
10、l GPS GeodeticControl Network 20002000国家GPS大地控制网由A、B级网、GPS一、二级网以及中国地壳运动观测网联合平差后，得到的以CGCS2000三维地心坐标为特征的高精度全国性大地控制网。3.25 最小二乘法least square method在残差满足VTPV为最小的条件下解算测量估值或参数估值并进行精度估算的方法，其中V为残差向量，P为其权阵。4控制点类型4.1控制点分类a）国家级GNSS连续运行基准站点。b）2000国家GPS大地控制网点。c）国家一、二、三、四等天文大地点。d）省级GNSS连续运行基准站点。e）省市卫星大地控制网C级点、D级点。
11、f）其他1954年北京坐标系、1980西安坐标系及相对独立的平面坐标系下的控制点。上述控制点中，a）、b）、c）三类控制点已有2000国家大地坐标系坐标；d）、e）两类控制点需要归算到2000国家大地坐标系；f）类控制点需转换到2000国家大地坐标系。4.2控制点精度表1 控制点的2000国家大地坐标系坐标精度要求点位类型 点位坐标精度GNSS连续运行基准站点 国家级GNSS连续运行基准站点 0.005m省级GNSS连续运行基准站点 0.01m2000国家GPS大地控制网点 0.03m国家天文大地点 一、二等点 0.11m三、四等点 0.20m省市级卫星大地控制网点 C级点内符精度 0.01m
12、D级点内符精度 0.02m4.3控制点用途高等级控制点可用于低等级控制网的外部控制；可用于1954年北京坐标系、1980西安坐标系坐标成果转换为2000国家大地坐标系坐标成果时计算坐标转换参数。a）国家级GNSS连续运行基准站点：可作为省级连续运行基准站网建设的控制点。b）2000国家GPS大地控制网点：可作为天文大地点控制点及相对独立的平面坐标系建立控制点。c）国家一、二等天文大地点：可作为三、四等天文大地点的控制点使用。d）国家三、四等天文大地点：可作为测图控制点使用；三等及以上天文大地点坐标成果可作为像控点的起算点。e）省级GNSS连续运行基准站点：点位坐标归算到2000国家大地坐标系后
13、，可作为省级、市、县城市基础建设控制网点。f）省市卫星大地控制网C级点、D级点：点位坐标归算到2000国家大地坐标系后，可作为建立相对独立的平面坐标系的控制点。g）其他大地测量控制点：点位坐标转换到2000国家大地坐标系后，可作为建立相对独立的平面坐标系的控制点。5控制点坐标转换模型及适用范围5.1总则在控制点坐标转换过程中，控制点所在坐标系涉及到1954年北京坐标系、1980西安坐标系、2000国家大地坐标系，上述坐标系的参考椭球基准及参数参见附录A。5.2转换模型a）布尔沙模型：用于不同地球椭球基准下的空间直角坐标系间的点位坐标转换。涉及七个参数，即三个平移参数，三个旋转参数和一个尺度变化
14、参数。b）三维七参数转换模型：用于不同地球椭球基准下的大地坐标系间的点位坐标转换，涉及三个平移参数，三个旋转参数和一个尺度变化参数，同时需顾及两种大地坐标系所对应的两个地球椭球长半轴和扁率差。c）二维七参数转换模型：用于不同地球椭球基准下的地心坐标系向大地坐标系的点位坐标转换，涉及三个平移参数，三个旋转参数和一个尺度变化参数。d）三维四参数转换模型：用于局部区域内、不同地球椭球基准下的地心坐标系向大地坐标系间的坐标转换，涉及三个平移参数和一个旋转参数。e）二维四参数转换模型：用于局部区域内、不同高斯投影平面坐标转换，涉及两个平移参数，一个旋转参数和一个尺度参数。对于三维坐标，需将坐标通过高斯投
15、影变换得到平面坐标，再计算转换参数。f）多项式拟合模型：不同范围的坐标转换均可用多项式拟合。有椭球面和平面两种形式。椭球面上多项式拟合模型适用于全国或大范围的拟合；平面拟合多用于相对独立的平面坐标系统转换。上述转换模型参见附录B。5.3适用范围a）布尔沙模型：适用于省级及全国范围的控制点空间直角坐标转换。b）三维七参数转换模型：适用于椭球面3及以上的省级及全国范围的控制点坐标转换。c）二维七参数转换模型：适用于椭球面3及以上的省级及全国范围控制点坐标转换。d）三维四参数转换模型：适用于2以内局部范围的控制点坐标转换。e）二维四参数转换模型：适用于小范围的控制点平面坐标转换、相对独立的平面坐标系
16、统与2000国家大地坐标系的联系。f）多项式拟合模型：椭球面多项式拟合适用于全国或大范围的大地坐标拟合；平面拟合多用于相对独立的平面坐标系统转换。5.4模型选取表2 控制点转换到2000国家大地坐标系的坐标转换模型选取控制点 转换模型 适用区域范围所属坐标系 坐标类型1980西安坐标系 大地坐标 三维七参数 椭球面3及以上的省级及全国范围二维七参数椭球面多项式拟合空间直角坐标 布尔沙 全国及省级范围三维四参数 2以内局部区域平面坐标 二维四参数 局部区域1954年北京坐标系 大地坐标 三维七参数 椭球面3及以上的省级及全国范围二维七参数椭球面多项式拟合空间直角坐标 布尔沙 全国及省级范围三维四
17、参数 2以内局部区域平面坐标 二维四参数 局部区域相对独立的平面坐标系 平面坐标 二维四参数 局部区域平面多项式拟合 局部区域6控制点坐标转换6.1省级GNSS连续运行基准站点坐标归算本方法适用于非2000国家大地坐标系下的省级GNSS连续运行基准站点坐标向2000国家大地坐标系的归算。6.1.1基准控制点选取选取我国周边稳定的IGS站、国内IGS站及国家级GNSS连续运行基准站至少10个作为省级GNSS连续运行基准站的最高级基准控制站。选取原则有：a）连续性：测站连续观测近3年(或以上)；b）稳定性：站点坐标时序稳定性好，具有稳定、可知的点位速度；c）高精度：速度场精度优于3mm/y；d）多
18、种解：至少3个不同分析中心的速度场残差优于3mm/y；e）平衡性：站点尽量全球分布；f）精度一致性：站点位置和速度的精度应当一致。6.1.2数据处理采用国内外成熟的高精度数据处理软件对GNSS连续运行基准站观测数据进行处理与平差，获得各站点在现ITRF框架、观测历元下的位置和基线向量。省级GNSS连续运行基准站作为省级2000国家大地坐标系框架基准，需将相邻省份的相邻GNSS连续运行基准站纳入本省GNSS连续运行基准站一同处理。6.1.3框架转换GNSS连续运行基准站点获得的观测历元的坐标转换为2000国家大地坐标系的坐标成果，需经历元归算、板块运动改正、框架转换三个步骤。a）历元归算。不同I
19、TRF框架对应的历元不同，需将不同ITRF框架下各参数（参见附录C）归算到同一历元下。框架间如无直接转换关系，可按间接方法转换；b）板块运动改正。计算框架点坐标从观测历元到需转换历元期间，由于板块运动引起的坐标变化值；c）框架转换。利用框架转换公式进行基准站坐标计算。框架转换公式参见附录D。6.2省市级卫星大地控制网点坐标转换省市级卫星大地控制网C级点坐标、D级点坐标转换到2000国家大地坐标系的方法参见6.1。6.3其他控制点坐标转换本方法适用于除上述控制点类型外的、其他1954年北京坐标系、1980西安坐标系及相对独立的平面坐标系下的大地测量控制点坐标向2000国家大地坐标系的转换。6.3
20、.1重合点选取原则选用具有原坐标系坐标和2000国家大地坐标系坐标的控制点作为重合点。重合点选取的基本原则为：等级高、精度高、分布均匀、覆盖整个转换区域、局部变形小。6.3.2转换参数计算转换参数计算按以下步骤进行：a）按转换区域选取适当的转换模型；b）按重合点选取原则选择重合点；c）将重合点代入转换模型，利用最小二乘法计算转换参数；d）用得到的转换参数计算重合点坐标残差；e）剔除残差大于3倍点位中误差的重合点；f）重复上述b）到e）的计算过程，直至重合点坐标残差均小于3倍点位中误差；g）最终用于计算转换参数的重合点数量与转换区域大小有关，但不得少于6个；h）根据最终确定的重合点，利用最小二乘
21、法计算转换参数。6.3.3坐标转换利用计算得到的转换参数，进行坐标转换，求得各控制点在2000国家大地坐标系下的坐标。6.3.4外部检核点数量要求利用未参与计算转换参数的重合点作为外部检核点，其点数不少于6个且均匀分布；7精度评定7.1省级GNSS连续运行基准站点坐标归算精度评定省级GNSS连续运行基准站点坐标归算后坐标精度在平面3cm以内。7.2其他控制点坐标转换精度评定采用内符合精度和外符合精度评定，依据计算转换参数的重合点残差中误差评估坐标转换精度，残差小于3倍点位中误差的点位精度满足要求。内符合精度计算公式参见附录E。附录A （资料性附录）常用坐标系椭球基准及参数A1 WGS 84大地
22、坐标系的地球椭球基本参数A1.1 地球椭球基本参数长半径a=6 378 137m扁率f=1/298.257223563地球引力常数（含大气层）GM=3 986 0051014m3s2正常化二阶带谐系数C2.0=484.166 85106地球自转角速度=7 292115105rad.s1A1.2 WGS84（G730）大地坐标系，其他地球椭球基本参数A1.1同规定。A2 2000国家大地坐标系的地球椭球基本参数长半径a=6378137m扁率f=1/298.257222101地球引力常数（含大气层）GM=3.986 0044181014m3s-2正常化二阶带谐系数J2=1.082629832258
23、10-3地球自转角速度=7.292 115105rad/sA3 1980西安坐标系的参考椭球基本参数长半径a=6378140m扁率f=1/298.257地球引力常数（含大气层）GM=3.9860051014m3s-2二阶带谐系数J2=1.0826310-3地球自转角速度=7.292115105rad/sA4 1954年北京坐标系参考椭球基本参数长半径a=6378245m扁率f=1/298.3附录B （规范性附录）坐标转换模型B1 布尔沙模型2 1 12 1 12 1 1x Z Yy Z Xz Y XX X T D R R XY Y T R D R YZ Z T R R D Z （B.1）式中：
24、1X ， 1Y ， 1Z 原坐标系坐标2X ， 2Y ， 2Z 目标坐标系坐标xT ， yT ， zT ，D， xR ， yR ， zR 七参数B2 三维七参数转换模型 faBBeaM BBf BeBeaN BBeMaN DBNeHN BBeMN RRRLBBNe LHM BNeHNLBBNe LHM BNeHN LtgBHN HeNLtgBHN HeN TTTB HM BLB HM LBLB HM LB BHN LBHN LHBL zyxzyx22222222 222 2 222 22 22 sin)sin1(1 “cossin1 )sin2()sin1( “cossin 00 sin)(
25、“cossin0 00“/coscossin cossin)(“/sincossin sinsin)( 1sin)1(cos)1( sin “)( cossinsin “)( sinsincoscos “)( cossin 0“cos)( cos“cos)( sin （B.2）式中：2e 第一偏心率的平方, 无量纲 2 22e f f M 地球椭球子午圈曲率半径，单位为米， 2/3222 sin1)1( BeeaM N 地球椭球卯酉圈曲率半径，单位为米 2/122 sin1 BeaN ，B，L，H 点位纬度、经度、大地高，经纬度单位为弧度，大地高单位为米B ， L ， H 点位在两个坐标系下纬
26、度差、经度差、大地高差。经纬度差值单位为弧度秒，大地高差值单位为米 角度与弧度间转换量，单位弧度秒， 180 3600/ a， a 椭球长半轴和长半轴差，单位为米f ， f 椭球扁率和扁率差，无量纲xT ， yT ， zT 平移参数，单位为米xR ， yR ， zR 旋转参数，单位为弧度秒D尺度参数，无量纲B3 二维七参数转换模型 faBBf BeBBeMaN DBBeMNRRRLL LtgBLtgB TTTMBM LBM LB BN LBN LBL zyx zyx “cossin1 )sin2(“cossin 00 “cossin0 0cossin 1sincos “cos“sinsin“c
27、ossin 0“coscos“cossin222 2 （B.3）式中：2e 第一偏心率的平方,无量纲M ，N 子午圈和卯酉圈的曲率半径，单位为米B，L， B ， L 点位纬度、经度，及其在两个坐标系下的纬度差、经度差。经纬度单位为弧度，其差值单位为弧度秒 角度与弧度间转换量，单位弧度秒， 180 3600/ a， a 椭球长半轴和长半轴差，单位为米f ， f 椭球扁率和扁率差，无量纲xT ， yT ， zT 平移参数，单位为米xR ， yR ， zR 旋转参数，单位为弧度秒D尺度参数，无量纲B4 三维四参数转换模型坐标转换模型实质上是采用 xT， yT， zT3个坐标平移量和1个控制网水平定向
28、旋转量作为参数。是以区域中心P0点法线为旋转轴的控制网水平定向旋转量，顾及1954年北京坐标系或1980西安坐标系平面坐标由于起始定向与2000国家大地坐标系的差异引起的坐标变化。 0000 000 000 sincoscoscos sincoscos sinsincos LBXLBY BXLBZ BYLBZTTTZYXZYX CC CC CCzyxCCCGGG （B.4）式中： GGG ZYX , 2000国家大地坐标系下的坐标，单位为米0B ， 0L 区域中心P0点的大地经、纬度，单位为弧度CCC ZYX , 大地坐标系（1954年北京坐标系或1980西安坐标系）坐标，单位为米xT ， y
29、T ， zT 坐标平移量，单位为米 为旋转参数，单位为弧度B5 二维四参数转换模型 1122 cossin sincos)1( yxmyxyx （B.5）式中：1x , 1y 原坐标系下平面直角坐标，单位为米2x , 2y 2000国家大地坐标系下的平面直角坐标，单位为米x , y 为平移参数，单位为米 为旋转参数，单位为弧度m尺度参数，无量纲B6 多项式拟合模型椭球面上拟合公式： 25243210 25243210 BbLbBLbBbLbbdL LaBaBLaLaBaadB （B.6）式中：B，L纬度、经度，单位为弧度ia ， ib 多项式拟合系数，通过最小二乘求解。平面拟合公式 yyy x
30、xx 12 12 （B.7）式中：1x 、 1y 原平面直角坐标2x 、 2y 目标平面直角坐标x 、 y 坐标转换改正量，用下式计算； 3828 273625423210 yaxya yxaxayaxyaxayaxaayx或其中 ia 为系数，通过最小二乘求解。B7 高斯投影正算公式 4442224222222 44442222422222 cos)5814185(1201cos)1(611cos cos)5861(7201cos)495(2415.0cos lBNtttlBtlBNy lBttlBtlBNtx （B.8）式中符号含义及计算方法参见B8高斯投影反算公式。B8 高斯投影反算公式
31、 422242222 44222222 )(8624285(1201)(21(611)(cos )(459061(3601)(935(1211)(2 fffffffffff ffffffffffff NytttNytNyBl NyttNyttNyyMtBB （B.9）式中：2e 第一偏心率的平方, 2 22e f f ， a bf a ， 21 eab ， bac 2e 第二偏心率， b bae 22 ， Be 222 cos ， tgBt W 计算中间过程变量， BeW 22 sin1 ，V 计算中间过程变量， BeV 22 cos1M 子午圈曲率半径， 33 2)1( VcW eaM N卯
32、酉圈曲率半径， VcWa X 子午线弧长。设有子午线上两点p1和p2，p1在赤道上，p2纬度为B，p1、p2间的子午线弧长X计算公式 )12sin10sin8sin6sin4sin2sin)(1( 2 BGBFBEBDBCBBarcBAeaX 式中：121210 12108 121086 1210864 12108642 12108642 83886081001524288090091310720693 8388608990995242883465131072315 104857655055262144103954096105307235 838860814864856553610395163
33、842205102410525615 52428829729713107272765409622051024525321583 1048576693693655364365916384110252561756445431 eG eeF eeeE eeeeD eeeeeC eeeeeeB eeeeeeA 底点纬度 fB 迭代公式：)( )(,)1( 120 i iii BF BFXBBAea XB 直到 ii BB 1 小于某一指定数值，即可停止迭代。式中12sin10sin8sin 6sin4sin2sin)1()( 2 BGBFBE BDBCBBarcBAeaBF 12cos1210cos1
34、08cos8 6cos64cos42cos2)1()( 2 BGBFBE BDBCBBAeaBF 附录C （资料性附录）ITRF 框架转换参数及其速率C1 转换参数符号表示及单位表C.1转换参数符号表示及单位转换参数 平移参数 尺度参数 旋转参数 参考历元符号表示 xT yT zT D xR yR zR t0计量单位 毫米（mm） 10-9(ppb) 毫角秒（mas，0.001“）速率符号 xT yT zT D xR yR zR计量单位 毫米/年（mm/y） 10-9 /年(ppb/y) 毫角秒/年(0.001“/y)C2从ITRF2000转换到ITRFyy的转换参数与速率表C.2从ITRF2
35、000转换到ITRFyy的转换参数与速率ITRFyy xT yT zT D xR yR zR t0速率(rates) xT yT zT D xR yR zRITRF97 6.7 6.1 -18.5 1.55 0.00 0.00 0.00 1997.0rates 0.0 -0.6 -1.4 0.01 0.00 0.00 0.02ITRF96 6.7 6.1 -18.5 1.55 0.00 0.00 0.00 1997.0rates 0.0 -0.6 -1.4 0.01 0.00 0.00 0.02ITRF94 6.7 6.1 -18.5 1.55 0.00 0.00 0.00 1997.0rat
36、es 0.0 -0.6 -1.4 0.01 0.00 0.00 0.02ITRF93 12.7 6.5 -20.9 1.95 -0.39 0.80 -1.14 1988.0rates -2.9 -0.2 -0.6 0.01 -0.11 -0.19 0.07ITRF92 14.7 13.5 -13.9 0.75 0.00 0.00 -0.18 1988.0rates 0.0 -0.6 -1.4 0.01 0.00 0.00 0.02ITRF91 26.7 27.5 -19.9 2.15 0.00 0.00 -0.18 1988.0rates 0.0 -0.6 -1.4 0.01 0.00 0.0
37、0 0.02ITRF90 24.7 23.5 -35.9 2.45 0.00 0.00 -0.18 1988.0rates 0.0 -0.6 -1.4 0.01 0.00 0.00 0.02ITRF89 29.7 47.5 -73.9 5.85 0.00 0.00 -0.18 1988.0rates 0.0 -0.6 -1.4 0.01 0.00 0.00 0.02ITRF88 24.7 11.5 -97.9 8.95 0.10 0.00 -0.18 1988.0rates 0.0 -0.6 -1.4 0.01 0.00 0.00 0.02C3ITRF2005到 ITRF2000的转换参数及其
38、速率表C.3ITRF2005到ITRF2000的转换参数及其速率转换参数 xT yT zT D xR yR zR t0速率(rates) xT yT zT D xR yR zR0.1 -0.8 -5.8 0.40 0.000 0.000 0.000 2000.0+/- 0.3 0.3 0.3 0.05 0.012 0.012 0.012rates -0.2 0.1 -1.8 0.08 0.000 0.000 0.000+/- 0.3 0.3 0.3 0.05 0.012 0.012 0.012C4ITRF2008到ITRFyy的转换参数及其速率表C.4ITRF2008到ITRFyy的转换参数及
39、其速率ITRFyy xT yT zT D xR yR zR t0速率(rates) xT yT zT D xR yR zRITRF2005 -2.0 -0.9 -4.7 0.94 0.00 0.00 0.00 2000.0rates 0.3 0.0 0.0 0.00 0.00 0.00 0.00ITRF2000 -1.9 -1.7 -10.5 1.34 0.00 0.00 0.00 2000.0rates 0.1 0.1 -1.8 0.08 0.00 0.00 0.00ITRF97 4.8 2.6 -33.2 2.92 0.00 0.00 0.06 2000.0rates 0.1 -0.5 -
40、3.2 0.09 0.00 0.00 0.02ITRF96 4.8 2.6 -33.2 2.92 0.00 0.00 0.06 2000.0rates 0.1 -0.5 -3.2 0.09 0.00 0.00 0.02ITRF94 4.8 2.6 -33.2 2.92 0.00 0.00 0.06 2000.0rates 0.1 -0.5 -3.2 0.09 0.00 0.00 0.02ITRF93 -24.0 2.4 -38.6 3.41 -1.71 -1.48 -0.30 2000.0rates -2.8 -0.1 -2.4 0.09 -0.11 -0.19 0.07ITRF92 12.8
41、 4.6 -41.2 2.21 0.00 0.00 0.06 2000.0rates 0.1 -0.5 -3.2 0.09 0.00 0.00 0.02ITRF91 24.8 18.6 -47.2 3.61 0.00 0.00 0.06 2000.0rates 0.1 -0.5 -3.2 0.09 0.00 0.00 0.02ITRF90 22.8 14.6 -63.2 3.91 0.00 0.00 0.06 2000.0rates 0.1 -0.5 -3.2 0.09 0.00 0.00 0.02ITRF89 27.8 38.6 -101.2 7.31 0.00 0.00 0.06 2000
42、.0rates 0.1 -0.5 -3.2 0.09 0.00 0.00 0.02ITRF88 22.8 2.6 -125.2 10.41 0.10 0.00 0.06 2000.0rates 0.1 -0.5 -3.2 0.09 0.00 0.00 0.02附录D （规范性附录）ITRF 框架与 2000 国家大地坐标系框架间转换D1 历元归算ITRF框架对应的转换关系是定义在具体历元基础上，参见附录C。ITRF2000与之前框架的转换关系为历元1997.0或1988.0，将附录C中参数转换归算到同一历元（2000.0）。首先将框架的参考历元t0（1997.0或1988.0）历元转换为t（2
43、000.0）历元。转换公式：0 00 00 00 00 0( 2000.0) ( ) (2000.0 )( 2000.0) ( ) (2000.0 )( 2000.0) ( ) (2000.0 )( 2000.0) ( ) (2000.0 )( 2000.0) ( ) (2000.0 )( 2000.0) x x xy y yz z zx x x ry y y rzT t T t T tT t T t T tT t T t T tR t R t R t mR t R t R t mR t R 0 00 0( ) (2000.0 )( 2000.0) ( ) (2000.0 )z z rt R
44、t mD t D t D t （D.1）式中：-94.84813681 10 rm 由毫弧秒(mas)到弧度的转换因子。D2框架间转换的间接方法如果框架间没有直接关系，可采用间接方法进行转换。例如将基于ITRF2005框架、当前历元的坐标成果转换到2000国家大地坐标系（ITRF97框架、2000.0历元）下坐标时，框架间的转换可采用如下方法：ITRF2005-ITRF97=ITRF2005-ITRF2000+ITRF2000-ITRF97即：分两步进行转换，首先从ITRF2005框架、当前历元转换到ITRF2000框架、2000.0历元；然后再从ITRF2000框架、2000.0历元转换到I
45、TRF97框架、2000.0历元。D3 板块运动改正不同框架、不同历元下的成果转换到2000国家大地坐标系所在的ITRF97框架和2000.0历元下，方法如下：根据IGS站速度值和观测历元 ct与需转换历元t的历元差，求出框架所对应历元下的坐标由于板块运动引起的坐标变化值，将观测历元转换为转换公式为( ) ccctt st t c st stXX XY Y t t YZ ZZ （D.2）式中：ct 观测历元t 需转换历元tX ， tY ， tZ 历元t时的坐标ctX ， ctY ， ctZ 观测历元 ct 时的坐标sX ， sY ， sZ IGS站速度值D4 框架点坐标计算根据框架间转换关系进行转换x z yy z xz y xITRF97 ITRFyy ITRFyyX X T D -R R XY = Y + T + R D -R YZ Z T -R R D Z （D.3）式中： TITRFyyX,Y,Z ITRFyy 框架下坐标 TITRF97X,Y,Z 2000国家大地坐标系下坐标xT ， yT ， zT ，D， xR ， yR ， zR 七参数附录E （规范性附录）精度评定公式E1 重合点残差VV = 重合点转换坐标值-重合点已知坐标值
