|||
对于survey mode的数据36-48h 水平方向精度在0.7-1.0mm 垂直在2-4mm,超过三天经比较好,超过5天垂直方向的达到mm级别量级,在研究区域有很少的连续站点,应该运行一个或者两个站点覆盖整个的测量时间,在每一年的同一时间去观测数据,最好的观测方式处理流动模式的数据时,通常用流动gps数据一起提供连续GPS数据,并且绑定ITRF站点框架合适的相对加权在绑定连续GPS数据时,如果选用不同的gps天线,则天线高度不一样,同时一般情况下多路径影响与大气延迟成为最大的误差考虑问题
原理以及软件构建模式:
比起用伪距测量能达到厘米-分米的水平,选择精密定位用相位测量,可以达到毫米级水平,因为相位是模糊的,开始我们只能依赖随着时间变化显示台站位置或者卫星位置的误差的特征,如果我们有充足的测量充足的卫星充足的时间,可以反演估计所有站点的相对位置与大气延迟,对于跟踪时段8小时,我们可以达到毫米级的估计即便没有解决2pi的模糊度在相位上,解决模糊度非常重要,对于基线的组成部分代表着测量误差的精度,而不是相位的改变率,对于所有高纬度的台站,东西方向的部分就比较弱。
准备:
一、工程规划
处理GPS测点数据 (注意:至少需要处理两个点 否则会出错)以最新的参考框架itrf2008为参考点,目前中国境内有14个igs站点(包括台湾两个)选择参考点,在建立工程目录时
尽量选取方位IGS站点(GUAO KIT3 KUNM POL2 SELE TASH XIAN URUM CHUM LHAZ BISH)作为框架 bish eveb xian 没有下到rinex数据,查看igs数据缺少的情况,保持每天的数据包含至少三个站左右
或者利用lfile-search的脚本去寻找最近公里之内的igs站点(尽量最近)
二、更新tables
将gamit软件包里的tabels进行更新,
主要更新地球方向参数ut1.,wob.(当天需要更新);
跳秒文件leap.sec(按年更新),月球星历太阳星历及其章动表soltab. luntab. nutabl.(按年更新),双差误差文件(c1和p1伪距分别得到接收机的时间,之间的误差)dcb.dat(按照月更新),天线相位中心模型文件antmod.dat(当有新的天线加入时,需要更新)
更新eop参数 sh_update_eop -series usno
更新dcb误差 sh_update_dcb -o Y
还有一些文件包括星历有关的等等,在此不再一一说明。
三、数据转换
如果是天宝数据,那么直接用它的软件用convertto rinex转换(只选择gps数据,gamit只能处理看gps数据,注意rinex改名,用点位的名字然后+gps周+0+年+o,转换完检查天线的序列号的输入以及天线高的设置,并且接收机的序列改成严格意义上的TRIMBLE MNETR8 全部大写,station.info文件格式要求非常严格)或者用runpkr00+teqc编写脚本进行转换成rinex格式。runpkr -d GT*.T00 teqc -tr d -week **** GT*.dat>????```0.10o ,当然此时需要编写脚本进行处理大量的数据。rinex格式很重要,teqc必须加入week的形式 doy可以给出gps week的数目。
在 gamit中用sh-runpkr的命令以及sh-make-rinex的命令可以实现转化为原始到rinex格式的数据,利用sh_get_times -sort -f *.12o>time.sorted.2012 查看数据的间隔以及天数
当然如果是徕卡的接收机等数据类型,一般代理公司都自带转换软件或者考虑利用teqc进行转换等等。
准备先验坐标有助于良好的处理gps数据,可以用teqc得到
1、利用teqc生成rinex文件与nav文件
teqc +nav abcd3650.14n +obsabcd3650.14o -tr d 12343650.dat
2、利用teqc在qc模式观测联合nav文件得到伪距估计的近似坐标
teqc +qc -nav abcd3650.14n abcd3650.14o
处理近似坐标也可以用命令sh_rx2apr 得到
四、数据准备
建立工程目录 四个字母kash(注意:最好用4字母命名工程文件夹名,保持工程名与工程文件夹名一致,因为sites.defaults等文件与sh_gamit等命令要求工程名为四字母)
下一级建立 brdc(导航文件) igs(星历文件) rinex raw(保存原始数据文件包括ion文件/eph文件/dat文件/mes文件)
在工程目录下 sh_setup -yr2012,链接gg/tables,将下载的tables复制到里面每一年的数据进行链接。
brdc:sh_get_nav -archive sopac -yr 2012-doy 241 -ndays 26
igs:sh_get_orbits -archive sopac -yr2012 -doy 211 -ndays 26;
sh_get_orbits -archive sopac -yr 2012-doy 241 -ndays 26 -makeg no
得到gfile文件,gigsfy.ddd
rinex:sh_get_rinex -archive sopac -yr2012 -doy 241 -ndays 26 -sd ../tables/sites.defaults -expt kash 或者直接跟站点名 -sites 下载rinex数据,检查数据的天数
如果遇到更新下载文件失败,可以选择更换archive为cddis重新尝试
五、准备设置文件
1.sites.defaults (测站列表)当然还有自己的点位,都设置进去
l# ftprnx= sites to ftp from rinex data archives.
# ftpraw =sites to ftp from raw data archives.
# localrx= sites names used to search for rinex files on your local system.
# (requiredin conjunction with rnxfnd path variable set in process.defaults).
# xstinfo= sites to exclude from automatic station.infoupdating.
# xsite = sites to exclude from processing, all days or specifieddays
# GLOBK:
# glrepu =sites used in the GLRED repeatability solution (default is to use all)
# glreps =sites used for reference frame definition (stabilization) in
# GLORGfor the GLRED repeatabilitiy solution (default is IGS list)
# glts= sites to plot as timeseries from GLRED repeatability solution (default is all)
(注意:glreps表示为框架站;glrepu计算重复性,生成时间序列,igs站可以不用,但是自己的站点可以设置。xstinfo表示不再从原始数据头文件读取)
查看sites.defaults有没有xstinfo这样sh_gamit就不会尝试调用sh_upd_stnfo去改变station.info文件信息
# may use the following
all_sites kash xstinfo localrx
备选、假定已经将前面的rinex格式准备完毕,不需要gamit进行自行ftp文件,这样让gamit进行更快的运算,如果有需要选择进行重复性的站点或者作为参考站点的站点可以进行另行设置。
2.station.info (测站信息文件)
假如在site.defaults里都罗列好了所有的站点后,那么可以直接利用
利用 sh_upd_stnfo -ref station.info -w station.info.mit -l sd
得到新的station.info.mit,这样得到的文件放到rinex文件夹下,然后重新命名为station.info
然后直接运行sh_upd_stnfo-files *.12o -rxo
或者直接将station.info 放入到rinex文件夹下,直接运行sh_upd_stnfo-files *.12o
备注:更新站点天线高度以及model一集序列号码对每一段的采样时段在tables目录下直接使用:
sh_upd_stnfo -files ../rinex/*.12o
3.lfile.
选择的参考站点可以利用精确坐标得到:
在 tables目录下
grepGUAO itrf08.apr > coor.apr
grepLHAZ itrf08.apr >> coor.apr
grepBISH itrf08.apr >> coor.apr
grepKIT3 itrf08.apr >> coor.apr
grepKUNM itrf08.apr >> coor.apr
grepPOL2 itrf08.apr >> coor.apr
grepSELE itrf08.apr >> coor.apr
grepTASH itrf08.apr >> coor.apr
grepCHUM itrf08.apr >> coor.apr
grepXIAN itrf08.apr >> coor.apr
grepURUM itrf08.apr >> coor.apr
grepEVEB itrf08.apr >> coor.apr
但是对于itrfo8.apr没有的站点以及自己的站点,当然现在已经更新到itrf14的apr啦
生成不了精确坐标的参考站点可以选择剔除
(注意:选取的参照点要尽量靠近要计算的点,有的IGS站可能在计算的这一天缺失,所以换一天的数据就行了)
自己的点也要用到采用双差定位
sh_rx2apr-site at012150.09o -nav ../brdc/brdc2150.09n -ref guao2150.09o -apr../tables/itrf08.apr
直接cat *.apr > coor.apr
将算出来的apr文件中的数据复制到kash/tables/coor.apr,然后 mv coor.apr lfile.
利用sh—make—lfile直接生成lfile
备选、利用手持的GPS定位得到坐标,利用tform转换为xyz坐标,
或者利用SOPAC SCOUT自动运行gamit得到 近似坐标,或者利用APPS
最好的方法,至少有一个站点利用这个方法得到 http://sopac.ucsd.edu/scout.shtml
或者利用teqc的qc full模式,得到坐标:
sh_get_nav -archivesopac -yr 2012 -doy 253 -ndays 4
teqc+qc -nav brdc2530.12n rinex/anzc2530.12o
4.process.defaults 基本不动
itrf08.apr 改一下
在10.50中已经默认为itrf08_comb.apr
在10.60可以改为itrf14.apr
如果进行处理小时的或者不是按照天数进行处理,尤其在地震学中,要考虑session,当然这一条也可以在sh_gamit 中直接进行设置,所以这个文件可以不用变化
5.sestbl.
选择松弛模式,是为了更好的与全球解进行融合,如果处理火山周围数据,以及短基线或者超短基线数据,可以选择baseline的模式
Choiceof Experiment = RELAX.
选择投影模式VMF1 (注意:其实也可以不改,使用全球投影模式GMF)
利用大气负载,大气网格文件,海潮负载
DMap=VMF1 ; GMF(default)/NMFH/VMF1
WMap=VMF1 ; GMF(default)/NMFW/VMF1
Outputmet =Y ; write the a priori met values to a z-file (Y/N)
Usemet.grid =N ; not yet supported
Usemap.grid =Y ; VMF1 grid file with mapping functions and ZHD
Useotl.grid =Y ; Ocean tidal loading grid file, GAMIT-format converted from OSO
Useatl.grid =Y ; Atmospheric tides, grid file, not yet available
Useatml.grid =Y ; Atmospheric (non-tidal) loading grid file from LU, converted to GAMIT format
为了轨道的确定,选择CM框架,主要是用来用到IGS轨道,大气负载用CM
下载vmf1grd.2009放到tables目录下
(直接链接)
在 kash/tables目录下map文件一定对应好年份每一年的map投影文件对应每一年。
ftp://everest.mit.edu/pub/GRIDS
rmmap.grid 建立新的map.grid并且链接到vmf1grd.2009
rmotl.grid
ln -svmf1grd.2012 map.grid
ln -sotl_FES2004.grid otl.grid
6.sittbl. 基本不动
所有列表中的IGS站点在gamit中 计算松弛解,如果在globk中用作参考框架则用紧约束
将n和e的约束改为0.005 区域站改成约束为10. 10. 10.
grep*_gps sittbl.|awk ‘{print $2=0.005,$3=0.005}’
将没有在文件中的约束,全部加入约束信息。自己的站点用99.9 99.9 99.9 igs站点用0.005 0.005 0.05
将没有在表内的站点(选择的igs站点)进行编辑,输入0.05 0.05 0.05
挑选自己认为约束的站点进行紧约束
E.g.,P488 P488_GPS NNN 0.003 0.003 0.01
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-10-20 01:35
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社