不同itrf参考框架neu坐标系统下的速度场转换

《不同itrf参考框架neu坐标系统下的速度场转换》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、第33卷第6期大地测量与地球动力学Vo1．33No．62013年12月JOURNALOFGEODESYANDGEODYNAMICSDec．，2013文章编号：1671—5942(2013)06-0130-03不同ITRF参考框架NEU坐标系统下的速度场转换郝明朱良玉(中国地震局第二监测中心，西安710054)摘要利用GLOBK输出的NEU坐标系与空间直角坐标系、大地坐标系之间的几何关系，顾及不同ITRF参考框架下三维速度值之间的转换参数，推导了不同ITRF参考框架下地壳水平运动速度之间的直接转换公式，过程清晰明了，易于编程实现。结果表明，ITRF2005与I

2、TRF2008下水平运动速度场存在明显的系统偏差。关键词ITRF2005；ITRF2008；NEU坐标系统；速度场转换；GPS中图分类号：P226．3文献标识码：ATRANSFoRMATIoNoFVELoCITYFIELDWITHINNEUCooRDINATESYSTEMBETWEENVAIUoUSITRFSHaoMingandZhuLiangyu(SecondMonitoringCenter，ChinaEarthquakeAdministration，Xi’an710054)AbstractUsingthegeometricrelationshipamon

3、gNEUcoordinatesystemobtainedwiththeGLOBKsoftware，Cartesiancoordinatesystem(XYZ)andellipsoidalcoordinatesystem(BLH)，takingintoaccountthethree-dimen—sionalvelocitytransformationparametersbetweenvariousITRFs，adirectconversionformulaofcrustalhorizontalvelocitywithvariousITRFsiSderived．

4、ThistransformationprocessiSclearandeasytobeprogrammed．TheresultdemonstratesthatthereexistsasignificantsystematicbiasinthehorizontalmovementvelocityfieldwithintheITRF2005andITRF2008．Keywords：ITRF2005；ITRF2008；NEUcoordinatesystem；velocityfieldtransformation；GPS了l2个版本。1引言每当IERS公布新的ITR

5、F版本，就会涉及将原随着空间定位技术的发展，使建立和维持基于有参考框架下的坐标和速度转换到新的框架下。目空间定位技术的、长期稳定和具有较高密度的、动态前对于不同参考框架下的坐标转换已有了较多的研的全球性坐标参考框架成为可能。ITRF就是基于究J，但对于速度场转换的研究却不是很多。VLBI、SLR、GPS和DORIS等空间技术所建立起来如果用GAMIT／GLOBK软件处理GPS资料获的全球地面参考框架，是协议地球参考框架的具体取运动速度场时，得到的速度场结果是在NEU坐标实现。它提供了全球统一的、地心的、三维的和系统下，这也是我们在研究断裂带滑动速率、地壳运动

6、态的高精度地面坐标参考基准。从ITRF1988开动、板块运动和地球动力学等方面时，经常采用的一始到ITRF2008，国际地球自转服务(IERS)已公布种三维地壳运动速度场表现形式。对于不同参考框收稿日期：2013-07-06基金项目：中国地震局地震行业科研重点专项(201208009，201208006)作者简介：郝明，男，1982年生，博士，研究方向：GPS数据处理及其在地形监测中的应用．E—mail：ha_rag@163．COB第6期郝明等：不同ITRF参考框架NEU坐标系统下的速度场转换131料而。。㈩Altamimi等给出的ITRF2005与ITRF

7、2008Vx+(7)2转换步骤和转换参数~v,N,j～I']l∞=。[[萋]∞+]=(8)，X=(No+H)cosBcosL?Y=(N。+H)cosBsinL(1)根据式(8)还可以得出ITRF08下的协方【Z=[(1一e)Ⅳ。+H]sinB差矩阵为DITR阳8DITRF05+1M2DT朋T1(9)其中，平移参数的协方差矩阵D为对角阵[，对角线元素为(O．040．040．04)。r=口oBJle：roL(2)3算例分析I【=H从推导过程可以看出，ITRF2005与ITRF2008可获取NEU坐标值。式中，B、的单位为弧度，ro下的转换参数，实质是求解式(6)

8、中的MME为余纬为0。时的纬圈半径，即r。=0。si