1、1、一共观测 6个时段,两个已知 WGS84点,其余为新点,基线长度在 4-25km,平均长度 11km。2、LGO 基线解算设置全部默认,导入数据为由 LGO生成的 RINEX数据,和直接导入DBX效果一样,包含天线相位偏差。解出的基线全是固定解。环路报告为默认设置。平差方案中约束两已知 WGS84点,自定义先验方差。3、TGO 基线解算设置除阀值由 3.5改为 2.5,其余全部默认,导入数据为 LGO生成的RINEX数据(O 文件的天线高手工+0.0644m,即加了一个 L1的相位偏差),对卫星进行了个别观测段的删除,解出的基线全是固定解。环路报告为默认设置。平差方案中约束两已知 WGS8
2、4点,加权方案自动。4、对比结果:4.1 基线LGO的解算速度明显高于 TGO,且使用 LEICA的静态数据,LGO 的基线结果和 TGO的基线结果基本符合相应等级的复测基线标准,但是基本上 TGO的基线长度要长过 LGO的结果,精度上 TGO比 LGO差很多。4.2 闭合环路LGO直接默认输出,仅生成 6个最小独立环,相对闭合差在 0-1ppm之间。TGO 也默认输出,生成 96个三边环,包换同步环和异步环。如果按照规范检核,TGO 有 9个同步环闭合差超限(超限很小,因为本身同步环闭合限差就很小,极易超限),而 LGO则无法进行检验,但是我并没有对 TGO中环路闭合差超限的相关基线进行删除
3、,而是依照规范的要求对“全部基线”进行了平差。4.3 无约束平差在 WGS84系统下进行,LGO 平差结果和 TGO平差结果对比如下:两种软件解算精度相当,中误差差值在 0-3mm,但是 TGO还是略逊 TGO,由于无约束平差中起算点的不同,两套结果的绝对值无法进行比较,但是如果起算点相同,平差结果的差值应当在 1cm以内。这里要提到的一点时,无论使用哪种软件,如果工区内没有 WGS84起算点,如果又是使用无约束平差的结果来求取转换参数的话,那么整个网,包括以后的 RTK工序,都必须基于同一个 WGS84起算点,而如果一切默认,那么使用LGO和 TGO求出的转换参数就会不同,且点位差值会有一系
4、统偏差,以此为例,则纬度和纬度差值 0.03秒也就是 0.9m,大地高差值 2.5m,这是比较可怕的局面。必须加以注意。而如果无约束平差的结果没有用途,仅用来检核粗差,那么就可以不用顾忌软件自己选取那一个点作为其算,直接默认即可,在约束平差的时候,加上必要的位置约束以后,相当于一个系统的平移,不会影响最终的平差结果。4.4 约束平差约束平差方案一致,都约束已知 84点 PCM016和 PCM023,绝对约束,默认中误差为0。平差结果对比如下:为了便于对比,平差结果按照 UTM49S投影显示,可以看出:两种软件平差点位北东坐标差值在 0-4cm之间,高程差较小,1cm 左右。但是这样的结果是不具
5、备对比性和统计性的。平面坐标的较大差值,在前面的基线长度较差已经有所反应。因为不可能做到LGO和 TGO进入平差阶段的基线解完全一样,所以就无法保证最终平差的极度吻合。但是,这样的差值在物探测量控制网阶段,是完全可以接受的。同时佐证了前面的莱卡天线在 TGO中的改正方案的正确性。不单单是物探规范,几乎所有的行业规范,对基线同步环闭合差的检核都过于严格,而异步环却又过于松懈,在这一点上,或许 LGO是做的更加合理的:不严格区分同步环异步环,只是统计最小独立环的闭合差。我也坚信,莱卡的经验要优于制定规范的专家们,所以:在监理或者甲方对闭合环路要求不完全照搬规范的话,少量极小的同步环超限是完全可以接
6、受的,且很多时候是必须的,但是异步环是绝对不能够超限的。期待我们的 GPS控制测量规范能做出更加合理的调整。经过近期对 TTC、Pinnacle、Waypoint 三个软件的摸索,感觉真的没有再做进一步研究的必要性。因为界面太丑,更新太慢,操作繁琐。我喜欢简单严谨。还是返回头重新拾起尘封的 LGO。这次重新开始 LGO,秉承的原则是:不急不躁,稳扎稳打,不浅尝辄止,主要是深入帮助文件里面去,解决以前遇到的仍不明晰的细节问题。ok, lets go!1、练习内容:GPS 静态后处理2、软件版本:7.0 中文特别版,全部模块可用。3、样本数据:LGO 自带的 static数据。4、地方坐标系统:U
7、TM32N,白赛尔 1841椭球。5、新建项目:learn-gps static,坐标系统 WGS84。6、导入原始的 DBX观测数据。发现由于是样本数据,所以各点观测比较随意,最短时间有 9min,最长时间有60min,且各观测段犬牙交错,基线解算决定采用手工方法。因为是练习,所以这次很详细的阅读了相关的帮助文件。LGO 郑重的两个重要提醒。值得仔细揣摩:如果在外业作业中,在同一个作业(JOB)内,对同一个测量点的不同测量时段的静态数据采用同一点号作为点名进行保存,则在使用 LGO进行数据处理时,要特别注意正确选择所需要的测量时段的数据。否则,可能会导致出现严重的数据处理误差。如果在外业作业
8、中,在同一台接收机的同一个作业(JOB)内,对在不同测量点上所采集的静态数据采用同一点名进行保存,则在使用 LGO进行数据处理时,要特别注意首先使用“单点定位解(SPP)”功能区分出不同测量点的数据,再正确选择所需要的测量点及测量时段的数据。否则,可能会导致出现严重的数据处理误差。也就是说,在实际的基线解算中,我们应当根据预先的计划来合理的选取可能出现的基线,由于野外调度和实际情况的需要,往往会在我们不希望两个测站形成同步观测的时候出现了同步观测,这时候,就应当采用手工模式或者设置形成基线的最小时间间隔,来达到想要的结果。7、GPS 处理参数的设置GPS 处理由两个主要的部分组成:选择观测时段
9、和选择处理参数、处理本身。用户需要干预的是选择时段和选择处理参数。因为处理本身无需任何干涉A、LGO 的基线处理模式分手工和自动两种。手工处理模式可以根据用户结合实际情况及需要进行设计如何计算数据。自动处理模式是从选择的时段中自动处理根据一组约束条件组合而成的所有的合理的基线。它只能选择流动站,不能选择参考站。LGO 会自动选择合适的参考站。处理顺序依赖于“自动处理参数”中定义的参数。 AAAA,自动处理参数选项卡:a、只有当处理模式设置成自动时,才可以使用自动处理参数。自动处理参数实际上是协助用户在自动处理的模式下进行基线计算的选取。“公共时间数据的最短时间”:同步观测时间的最短时间。同步时
10、间短于 300s,不予解算。这在剔除个别较短重叠时段时十分有用。“最大基线长度”:解算基线的最大长度。此选项,我能想到的大概就是批量处理个别较长基线时用到,如长度超过 200km的个别基线使用精密星历。批处理比较有用。另外:此项可输入的最大值是 1000,我曾尝试解算超过 1000km的基线,结果宣告失败,改成手工模式,可以没有限制。“处理方式”:选择“全部基线”的话,LGO 按照基线最短时间和最大基线长度为前提处理所有可能的相关基线。选择“独立设置”的话,LGO 只处理一组独立的基线集合。在 N个点的网中中有 N-1条独立基线。但是线性相关的基线也可以视为独立基线(例如三角形的三个边),前提
11、是:三条边没有同时观测。“坐标选择策略”:选择“距离”,第一个参考点的最短的基线首先被计算,然后计算次短基线。选择“时间”,具有最长公共时间的基线首先被计算,然后计算次长基线。如果“按时段”勾选,则在满足处理模式和选择策略的两个条件下,选最长时段的点做为第一个参考站。这个选项似乎意义不大,但是如果有一点的一个观测时段巨长,可以考虑勾选,如长时间长距离连测的点。“使用浮点解基准”,如果一个项目仅仅只有浮点解存在,则允许使用该点作为参考站做进一步处理。这个是极端情况。遇到的可能性不大。“重算已处理基线”,如果勾选,则原来已经计算了并且存储过的基线将再次被计算。这个选项比较任性,特别对于基线比较多需
12、要多次设置才能完成最终所有基线的解算,一次只能出来几条合格的基线,不勾选此项,可以保证已经成功处理的基线不丢失解算。那什么时候才勾选呢?一般不勾选,但是应该也有勾选的时候,那就是闭合环路出了很大问题,不过那个时候似乎可以直接删除结果了,更直接。呵呵。“计算控制点间的基线”,此项似乎一般不使用。似乎。但是似乎第二第三时段复测基线的时候为了和第一时段统一基准的时候。BBBB、“策略”选项卡确实应该把策略先进行阐述,LGO 似乎设置的不合理。-。“频率”:默认设置是“自动”,LGO 会自动为最终的解算选择最好的频率或者频率组合。软件自动处理的理论背景如下:由于 L1和 L2频率不同,在电离层中产生的
13、信号延迟也不同,利用这两种频率的线形组合可以计算消除电离层影响。然而,L3 的解算也同样破坏了整周模糊度。当模糊度还未固定时,采用浮点解计算。对于长基线来说(例如,长度大于 80M),使用浮点解是不可靠的(除非模糊度值固定)。如果观测时间足够长,则根据系统说明,L3 浮点解是足够准确的。如果可以事先解算 L1和 L2的模糊度,在无电离层影响的线形组合中采用 L1和 L2的整周模糊度进行第二次处理。在使用固定模糊度时可以消除电离层扰动。在可以解算模糊度但无法消除电离层影响时(如,大于 15M的基线)通常更优先采用这种策略。对于短基线来说,使用无电离层影响的线形组合会增加噪声,反而不好。最好采用标
14、准的 L1+L2解算。在自动频率下,双频数据基线长度超过 15KM,使用 L3(消除电离层)解算,根据情况可以有 L3浮点解和 L3固定解。如果基线长度小于 15KM,将处理 L1+L2。LGO 的这一点和规范是一致的,或许是规范参考这个而制定的。就是基线长度小于 15KM无条件必须是固定解,不然必须补测。选择 L1或者 L2将强制使用特定的一种频率计算一个解。选择 L1+L2将强制使用 L1和 L2进行计算,而不再限制基线的长度。似乎在长于15KM的基线,可以用 L1+L2弄出模糊度,但是又没有办法消除电离层影响的时候,这样设置是不合理的。而应当使用 L3。选择 L3使系统不限制基线长度使用
15、 L3解。“模糊度固定到:”:这个数值规定了解算模糊度的最长的基线距离,默认80KM,可以设置更高。在这里,LEICA 技术人员似乎转圈子了,如果频率是自动,那么距离小于 15KM没有问题,必须是固定解。长度超过 15KM的基线将使用 L3解,距离再长的,会出现 L3浮点解,他们又说长度较长的基线的模糊度固定是没有意义的,对于长基线,如果保证了观测时间够长,L3 浮点解是没有问题的。L3 浮点解足以满足系统的标称精度。那为何他们还要设置成 80KM呢?或许是一个经验值吧。或许浮点更好,在环路的检验中可以看出来。“采样率”:可以有 1 2 3 4 5 6 10 12 15 20 30 60秒,直
16、接设置成全部使用就可以了。“对流层模型”:对流层是高度约 30KM的大气部分,它可能导致电磁波的传播延迟。为了计算该延迟必须知道对流层的折射因子影响。LGO 提供如下的模型:Hopfied 模型、Simplified Hopfield 模型、Essen the first half computed from observations and the second half predicted orbit. The files are named according to the midpoint time in the file: 00, 06, 12, and 18 UT. Please
17、DO NOT write scripts to repeatedly check this area for new ephemeris files! For best reliability, you should use IGS Data Centers, automatically failing over to another one in case one is not available. Full details, including redundant access availability at IGS Global Data Centers, are available i
18、n the IGS Product table. This table indicates the most recent IGS ephemeris found by this server for recent days: IGS Final Orbit available for GPS weeks: 下面是我的理解:IGS共提供经过解算的三种类型的 GPS星历,钟差和地球定向。最终星历在 12天以后可用。快速星历在大约 17小时后可用。超快星历一天公布四次(UT 时间的 3点、9 点、15 点和 21点),包含 48小时的卫星轨道值,前 24小时的轨道值来自观测值后 24小时来自估算值
19、。文件名称依照该文件的中点时间,亦即:UT 时的 00点、06 点、12 点和 18点。请不要编写代码以执行重复性的在此区域检视新的星历文件。值得注意的是:你可以直接使用 IGS数据中心,以在其中一种星历不能使用的情况下可以快速转到另一种星历。全部的细节信息,包含其余的获取 IGS数据中心数据的途径,请参照“IGS 产品展台”这个表格显示了在本服务器中可以使用的最新的 IGS星历。(如上面的表格):今天是 2009年 11月 3日,计算出来的 GPS周和 GPS日分别是:1556 和 2(也就是星期三) 。从上表可以看出:今天可用的星历只有超快星历,如果今天在某一个时段进行了野外数据的采集,就可以参照时段进行星历的选择下载。前一天和再前一天,均只有快速星历可用。要想用到最终星历,那就要推到半个月之前了(也就是两周之前:1553 周)。下面给出了可以使用最终星历的 GPS周,打开 1553周: