昆明市CGCS2000坐标计算方法研究.docVIP

昆明市CGCS2000坐标计算方法研究 　　摘要： 高精度定位和定规软件GAMIT/GLOBK，在处理长基线和连续时段的静态定位GPS数据方面有很多优势。本文阐述了GAMIT/GLOBK的GLOBK模块获取全球框架下坐标的基本步骤。然后结合云南省部分CORS数据，阐述获得区域站ITRF坐标的数据处理过程，说明可以通过GAMIT/GLOBK得到区域站的高精度ITRF坐标。最后用MATLAB实现将区域站的ITRF05解算成果转换为高精度CGCS2000坐标。 　　Abstract： As a high precision positioning and orbit determination software， GAMIT/GLOBK has many advantages in handling long baseline and continuous-time static positioning GPS data. In this paper， the basic steps of data processing in solving global framework coordinates are described. CORS data of Yunnan Province is processed to demonstrate the data processing of getting ITRF coordinates of regional stations. It is concluded that high-precision ITRF coordinates of the regional stations can be obtained by GAMIT/GLOBK. At last， GCCS2000 coordinates with high-precision are obtained on the basis of the regional station on ITRF05 by MATLAB. 　　关键词： GAMIT/GLOBK；CGCS2000；CORS；数据处理；坐标转换 　　Key words： GAMIT/GLOBK；CGCS2000；CORS；data processing；coordinate transformation 　　中图分类号：P228.4 文献标识码：A 文章编号：1006-4311（2013）25-0168-03 　　0 引言 　　GAMIT由许多程序模块组成，在解算长基线方面应用广泛，GLOBK主要利用卡尔曼滤波方法来综合处理多元数据。为了不对ITRF实际点位坐标造成较大影响，应该进行ITRF坐标向CGCS2000的转换。国内外很多学者对求解ITRF坐标及ITRF坐标向CGCS2000的转换进行了研究。Herring T A， King R W等[1]对GAMIT/GLOBK的安装和使用给予详细的说明，对获得ITRF坐标也进行了一定的研究；占伟等[2]对GLOBK数据处理结果的误差分布进行了分析，得出了误差分布比较均匀，N、E分量方向的精度明显优于U分量方向的精度的结论；魏子卿[3]详细介绍2000中国大地坐标系的定义和实现，最后对CGCS2000坐标系的实质进行了说明。 　　ITRF地球参考框架是国际公认的精度最高、稳定性最好的参考框架，所以有必要将GPS区域网统一到ITRF中来[4]。本文主要具体实践运用GLOBK模块求区域站的ITRF坐标，并从中总结精度的规律，最后用MATLAB实现ITRF坐标向CGCS2000的转换。 　　1 GLOBK与转换模型 　　1.1 GLOBK介绍 GLOBK主要由htoglb、glred、globk和glorg四个模块组成。 htoglb是将GAMIT的H-文件转为GLOBK认可的二进制的h文件；glred用于计算单天解的重复性；globk可以估计多天解的测站坐标和速度，根据区域网和全球网中的公共测站和卫星的协方差阵进行绑定，得到WGS84框架下每天的坐标.prt文件；glorg根据指定IGS站的两套坐标得到ITRF坐标.org文件。本文主要利用GLOBK模块获得区域站的ITRF05坐标。 　　1.2 转换模型 由于缓慢的地壳运动，GPS测站在参考框架内的位置随时间变化。当两个框架处在不同的历元上时，需要根据地壳运动速度先将框架的坐标归算到统一的历元上来，目标历元下的坐标计算公式如下： 　　X（t）=X（EPOCH）+X*（t-EPOCH）（1） 　　其中X（t）和X（EPOCH）分别是GPS站在归算历元t和参考历元EPOCH的坐标，X是测站速率。 　　历元统一后进行框架转换。坐标转换模型采用布尔