大地测量参考框架.ppt

《大地测量参考框架.ppt》由会员分享，可在线阅读，更多相关《大地测量参考框架.ppt（169页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、1 概论概论 大大地地基基准准（GeodeticDatum）：用以代表地球形体的旋转椭球，建立大地基准就是求定旋转椭球的参数及其定向（椭球旋转轴平行于地球的旋转轴，椭球的起始子午面平行于地球的起始子午面）和定位（旋转椭球中心与地球中心的关系）。ab12/26/20221大大地地测测量量参参考考系系统统（Geodetic ReferenceSystem）：坐标参考系统、高程参考系统、重力参考系统1）坐标参考系统：以旋转椭球为参照体建立的坐标系统，分为大地坐标系和空间直角坐标系两种形式。2）高程参考系统：以大地水准面为参照面的高程系统称为正高，以似大地水准面为参照面的高程系统称为正常高，以旋转椭球

2、面为参照面的高程系统称为大地高。3）重力参考系统：重力观测值的参考系统坐标系原点、坐标轴、尺度及其有关计算公式12/26/2022212/26/20223大地测量参考框架大地测量参考框架(GeodeticReferenceFrame)：是大地测量参考系统的具体实现，是通过大地测量手段确定的固定在地面上的控制网（点）所构建的，分为坐标参考框架、高程参考框架、重力参考框架。（http:/）12/26/20224国家平面控制网是按控制等级和施测精度分为一、二、三、四等网,含三角点、导线点共154348个。12/26/20225国家高程控制网按控制等级和施测精度分为一、二、三、四等网，共有水准点成果1

3、14041个，水准路线长度为416619.1公里。12/26/20226国家重力基本网是确定我国重力加速度数值的参考框架。2000国家重力基本网包括21个重力基准点和126个重力基本点。12/26/202272000国家GPS控制网由国家测绘局高精度GPSA、B级网,总参测绘局GPS一、二级网，中国地壳运动观测网组成,共2609个点。12/26/20228讨论题：1大地基准、坐标系统、参考框架之间的关系。2“从整体到局部”的测量原则是如何通过坐标参考框架体现的？3大地原点、水准原点在建立大地测量参考框架中的作用是什么？12/26/20229参考椭球实例：贝塞尔椭球(1841年)，克拉克椭球(1

4、866年)，海福特椭球(1910年)和克拉索夫斯基椭球(1940年)等名称年代长半径(m)1/f贝塞尔椭球(Bessel)18416377397299.15克拉克椭球(Clarke)18666378206294.98海福特椭球(Hayford)19106378388297.00克拉索夫斯基椭球(Krassovsky)19406378245298.312/26/202210参考椭球大小、定位与定向参考椭球大小、定位与定向选择或求定椭球的几何参数(长半径 a和扁率)确定椭球短轴的指向(椭球定向)确定椭球中心的位置(椭球定位，建立大地原点）12/26/202211椭球定向椭球定向12/26/2022

5、12一点定位：椭球中心位置由大地原点的大地坐标所确定椭球定位12/26/202213多点定位：椭球中心位置由一组大地点的大地坐标所确定，大地原点的起算数据按下式求得。12/26/202214大地原点和大地起算数据大地原点和大地起算数据大地原点也叫大地基准点或大地起算点，大地原点的经纬度/大地高/至某一固定点的大地方位角称为大地起算数据。12/26/202215现代大地测量基准/卫星大地测量基准（几何特征+物理特征）：总地球椭球（椭球中心与地球质心重合，椭球旋转轴与地球旋转轴重合，椭球的起始子午面与地球的起始子午面重合，在全球范围内椭球面与地球表面最佳拟合）地球椭球的四个基本常数：地球椭球赤道半

6、径a，地心引力常数GM，地球重力场2阶带谐系数J2（由此导出椭球扁率f,）和地球自转角速度w。12/26/202216定义卫星大地测量基准，将涉及到地球重力场模型、地极运动模型、地球引力常数、地球自转速度等。不同大地测量基准的差异对坐标的影响，可根据公共点的大地观测数据求得，并进而求解出转换模型，实现不同基准下的坐标转换，但由于观测误差的存在，导致转换模型误差，其精度取决于公共点的数量和分布、观测精度、数据处理方法等。12/26/202217总地球椭球实例：WGS84,GRS80WGS84(GPS)GRS80(ITRS)a63781376378137fJ21/298.2572235631082

7、62.99910-81/298.25722210110826310-8(rad/s)729211510-11729211510-11GM(m3s-2)3986005108398600510812/26/202218WUHAN -2267749.162 5009154.325 3221290.762 BEIJING -2148743.784 4426641.236 4044655.935 SHANGHAI -2831733.268 4675666.039 3275369.521 KUNMING -1281255.473 5640746.079 2682880.117 URUMQI 193030.

8、873 4606851.324 4393311.421 LHASA -106937.669 5549269.591 3139215.762 WUHNBJFSSHAOKUNMURUMLHAS12/26/2022193坐标参考系统坐标参考系统 以参考椭球为基准的坐标系，叫做参心坐标系；以总地球椭球为基准的坐标系，叫做地心坐标系。无论是参心坐标系还是地心坐标系均可分为空间直角坐标系和大地坐标系两种，它们都与地球体固连在一起，与地球同步运动，因而又称为地固坐标系，以地心为原点的地固坐标系则称为地心地固坐标系(ECEF)，主要用于描述地面点的相对位置；另一类是空间固定的坐标系，与地球自转无关，称为惯性坐

9、标系或天球坐标系，主要用于描述卫星和地球的运行位置和状态。12/26/202220坐标系统是由坐标原点位置、坐标轴的指向、尺度和相关的计算模型所定义的。对于地固坐标系，坐标原点选在参考椭球中心或地心；坐标轴的指向具有一定的选择性，国际上通用的坐标系一般采用协议地极方向CTP(ConventionalTerrestrialPole)作为Z轴指向，因而称为协议坐标系；尺度采用国际标准长度单位；实现方式为大地测量理论、技术与方法。地球旋转轴的指向 1)空间指向的变化（岁差、章动）2)地球旋转轴相对于地球内部结构的变化（极移）12/26/20222112/26/202222空间指向的变化：岁差(pre

10、cession),章动(nutation)12/26/202223地球旋转轴相对于地球内部结构的变化：极点的地球旋转轴相对于地球内部结构的变化：极点的变化变化(极移极移,polar motion，国际协议原点，国际协议原点CIO)地球自转轴相对地球体的位置并不是固定的，地极点在地球表面上的位置是随时间而变化的，这种现象称为地极移动，简称极移。某一观测瞬间地球北极所在的位置称为瞬时极，某段时间内地极的平均位置称为平极。12/26/202224国际天文联合会(IAU)和国际大地测量与地球物理联合会(IUGG)在1967年于意大利共同召开的第32次讨论会上，建议采用国际上5个纬度服务(ILS)站以1

11、9001905年的平均纬度所确定的平极作为基准点，通常称为国际协议原点CIO(ConventionalInternationalOrigin)，它相对于19001905年平均历元1903.0。另外国际极移服务(IPMS)和国际时间局(BIH)等机构分别用不同的方法得到地极原点，与CIO相应的地球赤道面称为平赤道面或协议赤道面。12/26/20222512/26/2022263.1 1954年北京坐标系年北京坐标系 1954年北京坐标系可以认为是前苏联1942年坐标系的延伸。它的原点不在北京，而在前苏联的普尔科沃。相应的椭球为克拉索夫斯基椭球。椭球参数有较大误差。参考椭球面与我国大地水准面存在着

12、自西向东明显的系统性的倾斜，在东部地区大地水准面差距最大达67m。定向不明确12/26/2022273.2 1980年国家大地坐标系年国家大地坐标系(1980(1980西安坐标系西安坐标系)1980年国家大地坐标系的特点是：采用1975年国际大地测量与地球物理联合会(IUGG)第16届大会上推荐的4个椭球基本参数。地球椭球长半径 a=6 378 140 m，地心引力常数 GM=3.986 0051014m3/s2,地球重力场二阶带球谐系数J2=1.082 6310-8,地球自转角速度=7.292 11510-5 rad/s。12/26/202228 参心大地坐标系是在1954年北京坐标系基础上

13、建立起来的。椭球面同似大地水准面在我国境内最为密合，是多点定位。定向明确。椭球短轴平行于地球质心指向地极原点JYD1968.0的方向大地原点地处我国中部，位于西安市以北60km处的泾阳县永乐镇，简称西安原点。大地高程基准采用1956年黄海高程系12/26/202229平差后提供的大地点成果属于1980年西安坐标系，它和原1954年北京坐标系的成果是不同的。这个差异除了由于它们各属不同椭球与不同的椭球定位、定向外，还因为前者是经过整体平差，而后者只是作了局部平差。不同坐标系统的控制点坐标可以通过一定的数学模型，在一定的精度范围内进行互相转换，使用时必须注意所用成果相应的坐标系统。12/26/20

14、22303.3 3.3 新新1954年北京坐标系年北京坐标系(BJ54新新)新1954年北京坐标系，是在GDZ80基础上，改变GDZ80相对应的IUGG1975椭球几何参数为克拉索夫斯基椭球参数，并将坐标原点(椭球中心)平移，使坐标轴保持平行而建立起来的。12/26/20223112/26/20223212/26/202233BJ54新的特点是：采用克拉索夫斯基椭球参数。是综合 GDZ80和BJ54建立起来的参心坐标系。采用多点定位，但椭球面与大地水准面在我国境内不是最佳拟合。定向明确，坐标轴与 GDZ80 相平行，椭球短轴平行于地球质心指向1968.0地极原点的方向 大地原点与 GDZ80

15、相同，但大地起算数据不同。12/26/202234 大地高程基准采用1956年黄海高程系。与 旧BJ54相比，所采用的椭球参数相同，其定位相近，但定向不同。旧BJ54的坐标是局部平差结果，而新BJ54是GDZ80 整体平差结果的转换值，两者之间无全国统一的转换参数，只能进行局部转换。12/26/2022353.4 3.4 地心地固坐标系地心地固坐标系 地心地固空间直角坐标系的定义是：原点O与地球质心重合，Z轴指向地球北极，X轴指向格林尼治平均子午面与地球赤道的交点，Y轴垂直于XOZ平面构成右手坐标系。地球北极是地心地固坐标系的基准指向点，地球北极点的变动将引起坐标轴方向的变化。12/26/20

16、2236地心地固大地坐标系的定义是：地球椭球的中心与地球质心重合，椭球面与大地水准面在全球范围内最佳符合，椭球的短轴与地球自转轴重合(过地球质心并指向北极)12/26/202237以协议地极CTP(ConventionalTerrestrialPole)为指向点的地球坐标系称为协议地球坐标系CTS(ConventionalTerrestrialSystem)，而以瞬时极为指向点的地球坐标系称为瞬时地球坐标系。在大地测量中采用的地心地固坐标系大多采用协议地极原点CIO为指向点，因而也是协议地球坐标系，一般情况下协议地球坐标系和地心地固坐标系代表相同的含义。3.5 3.5 协议地球坐标系协议地球坐

17、标系12/26/20223820世纪60年代以来，美国和原苏联等国家利用卫星观测等资料，开展了建立地心坐标系的工作。美国国防部曾先后建立过世界大地坐标系(WorldGeodeticSystem，简称为WGS)WGS60，WGS66和WGS72，并于1984年开始，经过多年修正和完善，建立起更为精确的地心坐标系统，称为WGS84。12/26/2022393.6 WGS84世界大地坐标系世界大地坐标系WGS84是一个协议地球参考系CTS。该坐标系的原点是地球的质心，Z 轴指向BIH1984.0定义的协议地球极CTP方向，X轴指向BIH1984.0零度子午面和CTP赤道的交点，Y轴和Z、X 轴构成右

18、手坐标系 12/26/202240WGS84坐标系统采用的4个基本参数是：a=6 378 137mGM=3 986 005108m3s-2C2,0=-484.166 8510-6=7 292 11510-11rad/s12/26/202241为了改善WGS84系统的精度，1994年6月，由美国国防制图局(DMA)将其和美国空军(AirForce)在全球的10个GPS跟踪站的数据加上部分IGS站的ITRF91数据，进行联合处理，并以IGS站在ITRF91框架下的站坐标为固定值，重新计算了这些全球跟踪站在1994.0历元的站坐标，更新为WGS84(G730)WGS84(G730)1996年，WGS

19、84坐标框架再次进行更新，得到了WGS84(G873)2004年进一步更新为WGS84(G1150)12/26/202242WGS84的体现与维持的体现与维持12/26/202243讨论题1地球旋转轴的指向变化特点及其与大地测量坐标系的关系。2岁差、章动、极移的不同点是什么？3WGS84的Z轴指向哪里？12/26/2022444 国际地球参考框架国际地球参考框架ITRFITRF（InternationalTerrestrialReferenceFrame）是由IERS（InternationalEarthRotationService）提供的国际地球参考框架，其构成是基于甚长基线干涉VLBI、

20、激光测月LLR、激光测卫SLR、GPS和卫星轨道跟踪和定位DORIS等空间大地测量技术的观测数据。这些观测数据首先由不同技术各自的分析中心进行处理，最后由IERS中心局（IERSCB）根据各分析中心的处理结果进行综合分析，得出ITRF的最终结果，并由IERS年度报告和技术备忘录向世界发布，提供各方面的应用。12/26/2022454.1 ITRF系列系列IERSCB每年将全球站的观测数据进行综合处理和分析，得到一个ITRF框架，并以IERS年报和IERS技术备忘录的形式发布。自1988年起，IERS已经发布ITRF88，89，90，91，92，93，94，96，97，ITRF2000，ITRF

21、2005等全球坐标参考框架。目前，IGS各种轨道产品的坐标参考基准采用的是ITRF2005参考框架。12/26/2022464.1.1 ITRF2000与其他框架的转换与其他框架的转换TRANSFORMATION PARAMETERS AND THEIR RATES FROM ITRF2000 TO PREVIOUS FRAMESSOLUTION T1 T2 T3 D R1 R2 R3 EPOCH Ref.UNITS-cm cm cm ppb .001 .001 .001 IERS Tech.Note#RATES T1 T2 T3 D R1 R2 R3UNITS-cm/y cm/y cm/y

22、ppb/y .001/y.001/y.001/y-ITRF97 0.67 0.61-1.85 1.55 0.00 0.00 0.00 1997.0 27 rates 0.00-0.06-0.14 0.01 0.00 0.00 0.02 ITRF96 0.67 0.61-1.85 1.55 0.00 0.00 0.00 1997.0 24 rates 0.00-0.06-0.14 0.01 0.00 0.00 0.02 ITRF94 0.67 0.61-1.85 1.55 0.00 0.00 0.00 1997.0 20 rates 0.00-0.06-0.14 0.01 0.00 0.00 0

23、.02 12/26/202247 ITRF93 1.27 0.65-2.09 1.95 -0.39 0.80 -1.14 1988.0 18 rates -0.29-0.02-0.06 0.01 -0.11 -0.19 0.07 ITRF92 1.47 1.35-1.39 0.75 0.00 0.00 -0.18 1988.0 15 rates 0.00-0.06-0.14 0.01 0.00 0.00 0.02 ITRF91 2.67 2.75-1.99 2.15 0.00 0.00 -0.18 1988.0 12 rates 0.00-0.06-0.14 0.01 0.00 0.00 0.

24、02 ITRF90 2.47 2.35-3.59 2.45 0.00 0.00 -0.18 1988.0 9 rates 0.00-0.06-0.14 0.01 0.00 0.00 0.02 ITRF89 2.97 4.75-7.39 5.85 0.00 0.00 -0.18 1988.0 6 rates 0.00-0.06-0.14 0.01 0.00 0.00 0.02 ITRF88 2.47 1.15-9.79 8.95 0.10 0.00 -0.18 1988.0 rates 0.00-0.06-0.14 0.01 0.00 0.00 0.02 12/26/202248Transfor

25、mation Parameters from ITRF2005 to ITRF2000at epoch 2000.0T1 T2 T3 D R1 R2 R3mm mm mm 10-9 mas mas mas0.1 -0.8 -5.8 0.40 0.000 0.000 0.000 12/26/202249已知wuhn的ITRF2000下1997.0参考历元的坐标(m)及变化率(m/y)为：-2267749.1625009154.3253221290.762-.0325-.0077-.0119求wuhn的ITRF2000下2009年3月19日为参考历元的坐标。12/26/2022504.1.2 IT

26、RF与与WGS84的转换的转换-Parameters from ITRF90 to WGS84-Doppler realized system:T1(m)T2(m)T3(m)D(ppm)R1(“)R2(“)R3(“)0.060 -0.517 -0.223 -0.011 0.0183 -0.0003 0.0070-New realizations of WGS84 based on GPS data,such as WGS84(G730 or G873):These new WGS84 realizations are coincident with ITRF at about 10-centi

27、meter level.For these realizations there are no official transformation parameters.This means that one can consider that ITRF coordinates are also expressed in WGS84 at 10 cm level.12/26/2022514.2 4.2 IERShttp:/TheIERSwasestablishedastheInternationalEarthRotationServicein1987bytheInternationalAstron

28、omicalUnionandtheInternationalUnionofGeodesyandGeophysicsanditbeganoperationon1January1988.In2003itwasrenamedtoInternationalEarthRotationandReferenceSystemsService.12/26/20225212/26/20225312/26/202254IERS的任务主要有以下几个方面：维持国际天球参考系统(ICRS)和框架(ICRF)；维持国际地球参考系统(ITRS)和框架(ITRF)；为当前应用和长期研究提供及时准确的地球自转参数(EOP)。12

29、/26/202255IERS ProductsThe IERS maintains the following main products:Earth orientation dataConventionsInternational Celestial Reference SystemInternational Celestial Reference FrameInternational Terrestrial Reference SystemInternational Terrestrial Reference Framehttp:/www.iers.org/iers/publication

30、s/tn/tn31/Geophysical fluids data12/26/202256ITRF2000STATIONPOSITIONS（m)ATEPOCH1997.0ANDVELOCITIES(m/y)BJFS-2148743.7844426641.2364044655.935-.0444.0141-.0013WUHN-2267749.1625009154.3253221290.762-.0325-.0077-.011912/26/2022574.3IGS4.3.1 HISTORYInternationalGlobalPositioningSystem(GPS)ServiceforGeod

31、ynamics(IGS)formallybeganon1January1994.DuetotheexpansionofIGSobjectives,thenameoftheservicewaschangedtoInternationalGPSService(IGS)on1January1999.FollowingfurtherexpansionofIGS,integratingdatafromGLONASSsystemandplanningforthedeploymentofGalileosystem,thenamewaschangedtoInternationalGlobalNavigatio

32、nSatelliteSystem(GNSS)Service(IGS)on14March2005.12/26/2022584.3.2 GOALS AND OBJECTIVESThe IGS strives to:Provide the highest quality,reliable GNSS data and products,openly and readily available to all.Promote universal acceptance of IGS products,standards and conventions.Continuously innovate by att

33、racting leading-edge expertise and pursuing challenging projects and ideas.Seek and pursue new growth opportunities while responding to changing user needs.12/26/202259Sustain and nurture the IGS culture of collegiality,openness,inclusiveness,and cooperation.Maintain a voluntary organization with ef

34、fective leadership,governance,and management.12/26/2022604.3.3 Data&products Trackingdata(dailyRINEXfiles)(Global)Navigationmessagefiles(broadcastmessagesoftrackingstations)MeteorologicalFilesandSummaryfilesGPSsatelliteephemeridesEarthrotationparametersIGStrackingstationcoordinatesandvelocitiesGPSsa

35、telliteandIGStrackingstationclockinformationZenithpathdelayestimates12/26/2022614.3.4 IGS组织主要由以下几部分组成组织主要由以下几部分组成：1)跟踪站网（Tracking Stations）IGS跟踪站网由全球24小时全天候观测的GPS跟踪站组成，IGS跟踪站划分为全球跟踪站（GlobalSites）、区域跟踪站（RegionalSites）和地方跟踪站（LocalSites）。建立连续观测跟踪站网的目的是计算卫星轨道，确定地球参考框架及地球自转参数等。12/26/202262Thereare382stat

36、ionsasof3Apr2006.Eachstationisassigneda4-charactersiteID.TherawstationdataisconvertedtoRINEXformatandaccessedthroughInternet.12/26/2022632)操作中心(OperationalCenters)该中心负责一个地方跟踪站网的运作，进行GPS原始观测数据的格式转换，以及为网内各跟踪站提供管理和技术支持。3）数据中心(Data Centers)全球数据中心（GlobalDataCenter）主要任务是为分析中心及外部用户提供数据服务。包括：CDDIS:CrustalDy

37、namicsDataInformationSystemSIO:ScrippsInstitutionofOceanographyIGN:InstitutGeographiqueNational12/26/202264区域数据中心（ReginalDataCenter）区域数据中心的目的是从数个地方数据中心收集数据，满足本区域的需求，减少数据传输流量，并且将数据传输至全球数据中心，包括：BKG:BundesamtfuerKartographieundGeodaesieAUSLIG:AustralianSurveyingandLandInformationNRCAN:NaturalResourceso

38、fCanadaGODC(formerlyCIGNET):GeoscienceLaboratory,NOAAJPL:JetPropulsionLaboratory12/26/202265地方数据中心（LocalDataCenter）地方数据中心与跟踪站直接相联系，他们的任务是进行数据格式转换、压缩数据、备份数据及传输数据至区域数据中心。4）分析中心(AnalysisCenters)分析中心每天的基本任务是从全球数据中心获取GPS跟踪站观测数据，独立地进行数据处理和分析，估计得到GPS卫星轨道、地球自转参数和站坐标及其速度，并将各自的分析成果返回给数据中心。分析中心包括NRCan、GFZ、JPL、

39、CODE、ESA、SIO和NGS等七个。12/26/2022665）协调分析中心（Analysis Centers Coordinator）协调分析中心负责监督分析中心的工作，并且对各个分析中心的产品进行质量分析、评价，最后进行加权平均，得到最终的IGS产品。并在数据采集后十天内，发送给IGS数据中心和IGS中心局信息系统CBIS。6）中心局(Central Bureau)中心局主要负责IGS日常工作，包括组织会议，制定标准及出版相关出版物等。由中心局维护的中心局信息系统CBIS，包含IGS的各种产品和信息，设有web网站和ftp服务器。位于法国巴黎的全球数据中心IGN设有其ftp服务的镜像站

40、。用户通过web网站或ftp服务器都可以访问IGS跟踪站数据及IGS轨道产品。12/26/2022674.3.5 获取获取IGS观测数据和产品观测数据和产品1)文件命名规则文件命名规则ThetrackingdataismadeavailableindailyRINEXfiles.AlltrackingsitesgenerateatleastadailyRINEXobservationandnavigationmessagefile.Thefilesareorganizedbytype,year,dayofyear,GPSweek,dayofweek,siteID.12/26/202268Loc

41、ating IGS data,products,and format definitionsKey to directory and file name variablesddayofweek(0-6)ssss4-characterIGSsiteIDor4-characterLEOIDddddayofyear(1-366)ww2-digitweekofyear(00-51)hh2-digithourofday(00-23)wwww4-digitGPSweekhsingleletterforhourofday(a-x=0-23)yy2-digityearmmminuteswithinhouryy

42、yy4-digityearleonameoflow-earthorbitingsatelliteSoftwaretouncompressfilesendingin.Zisavailablefromhttp:/www.gzip.org(ThiswillassistyouinlocatingandusingIGSdataandproducts.ItisnotanexhaustivecompilationofallfilesanddirectoriesavailableattheDataCenters)12/26/202269Observation Files Naming conventionss

43、sssdddf.yyD.Z,ssssdddf.yyO.Z(ftp:/)Navigation Message Files Naming conventions ssssdddf.yyN.Z(ftp:/)Global Navigation Message Filebrdcdddf.yyN.Z(ftp:/)Convert software between.yyD and.yyO:CRX2RNX.exeandRNX2CRX.exe(ftp:/igscb.jpl.nasa.gov/pub/software/RNXCMP_2.4.2/dos/exe_win2000/)12/26/2022702)How t

44、o Access IGS Tracking DataIdentification of the Tracking Siteshttp:/Data DownloadConnecttothecenterbyftporhttpanddownloadthefiles.Forexample:ftp:/pubftp:/ftp:/ftp:/pub/product/Example:wuhn0230.05o,wuhn0230.05n12/26/2022713）DownloadPreciseEphemerisigs13170.sp3(3Apr.2005)4)GPSCalendarforMarch,200812/2

45、6/2022725 ITRF框架坐标的传递框架坐标的传递方法：与IGS连续运行跟踪站的数据进行联合处理，得到待求点与IGS站之间的准确的基线向量，从而求得测区内一点的ITRF坐标，再以该点作为已知点进行控制网平差，得到全部其他点的ITRF坐标。特殊性：由于在我国IGS站较少，测区与IGS站的距离往往很长，甚至达到1000km以上，解算这种超长基线，一般需要特殊软件，我国用的较多的是Gamit.应用实例(某工程控制网)12/26/20227312/26/202274WUHN-2267749.4617 5009154.2928 3221290.6716BJFS-2148744.1162 44266

46、41.2776 4044655.8947LHAS -106938.0821 5549269.5274 3139215.8675SHAO-2831733.5363 4675665.9494 3275369.397900G2-2081460.5626 5097545.5610 3208420.9766 0.31 0.55 0.3712/26/20227512/26/20227612/26/202277为什么要用为什么要用GAMIT(1)在国家GPS控制网和城市GPS骨架网中,网的边长在几十至几百公里之间,要求精度达到10-7或更高；10-710-9的边长相对精度和毫米级的点位绝对精度测定地壳形变、

47、板块运动、地球自转和极移；精度要求是一般的商用软件系统所不能达到的,因此需要一个方便使用的、高精度的GPS 定位软件系统。12/26/202278为什么要用为什么要用GAMIT(2)美国麻省理工学院(MIT)等研制的GAMIT/GLOBK软件；瑞士伯尔尼大学天文研究所研制的BERNESE软件；美国宇航局(NASA)喷气推进实验室(JPL)研制的GIPSY/OASIS数据处理软件。这些软件都具有确定GPS卫星轨道的功能，都同时可处理多种型号接收机的观测数据。此外还有JPL的综合处理软件QOCA等12/26/202279GAMITGAMIT的介绍的介绍用来分析GPS载波数据，对一系列参数进行估计。

48、-测站的坐标(假定各个历元是不动的)-卫星轨道参数(初始条件、推估参数、相位中心的偏移)-地球定向参数(EOP)-大气延迟参数(Time-dependent天顶方向和倾斜方向的延迟)-载波的整周模糊度12/26/202280GLOBK的介绍的介绍对GPS基线、侧视雷达(SLR)、甚长基线测量(VLBI)数据的进行参数估计；联合这些数据，确立协方差阵来估计其他各种参数-测站速度（每个位置以mm/yr的时间积分）-卫星轨道的调整参数-EOP参数估计-地动方面的参数的估计(co-seismicdisplacements)12/26/202281GAMIT的发展里程的发展里程1970晚期，MIT正在开

49、发GPS接收机时期开始初始代码起源于1960-1970计算行星星历和VLBI软件1987年是Unix系统1992年IGS服务使得开始软件的自动计算处理计划mid-1990s完全自动处理，包括IGS持续跟踪站和GPS测量大会站的数据12/26/202282GAMIT主要模块主要模块Arc(GPS卫星轨道积分)Model-计算载波的理论值-计算载波的相关的部分参数的估值Autcln-清理GPS观测数据（粗差剔除）和修复周跳Solve-利用最小二乘来进行参数估计-计算合适解12/26/202283GAMIT次要模块次要模块计算GPS卫星的偏航表海洋潮汐加载系数-直接的测站观测值-格网的插值估计钟差模

50、型-卫星种的频率差-地面钟的模型载波残差显示-按照时间系列-按照星空图12/26/202284GAMIT/GLOBK 数据处理流程图数据处理流程图12/26/202285数据准备数据准备建立工作目录,如处理qingdao的数据就建立目录/qingdao；在此目录中再以年月日为目录名建立子目录/qingdao/322同时在这一级建立/qingdao/tables 目录。从相关机构网站上找到对应年和时段的各种信息文件放入tables表中 12/26/202286需要准备的文件需要准备的文件(1)GAMIT控制文件初始卫星轨道-广播星历-IGS站的精密星历（SP3）-GAMITg-files地球方向