GAMIT/GLOBK长基线静态精密单点定位解算方法 |
您所在的位置：网站首页 › 华测一直单点解 › GAMIT/GLOBK长基线静态精密单点定位解算方法 |
GAMIT/GLOBK长基线静态精密单点定位解算方法
|
GAMIT/GLOBK长基线静态精密单点定位解算方法
GAMIT/GLOBK组合平差解算方式，通过利用全球分布的IGS台站观测文件、GPS卫星精密星历、精密钟差、对流层、电离层等各类改正模型可以获取测站高精度的三维坐标信息。其解算过程中采用IGS站作为已知控制站点，因此待解算测站点可以不考虑基线距离的限制、区域网站点间共模误差的影响，可以直接实现待解算测站精密单点解算，大大解决了常规GPS测量依赖公共站点、短基线解算的缺点。 分步处理与批处理方式进行解算，具体解算流程如下图（本人论文中的图）
工程目录用4个英文小字母建立：例如，test。本人习惯。 然后建立工程目录test，之后cd 进入到工程中文件夹中，分别建立年积日“100”，globk平差“glbf”、“gsoln”（这两个目录也可以不用新建，后续平差会自己生成），“brdc”，“igs”和“rinex”目录并依此下载所需要的导航文件n，精密星历文件sp3和观测文件o： 1.1下载brdc0940.18n（例）文件到brdc目录： for ((i=1;iccdoy}0.12n.Z fi done 1.2 下载igs19953.clk和clkigs19953.sp3到igs目录 for ((line=10;line doyfile gpsweek=`grep "GPS Week" doyfile | awk '{print $3}'` doyofweek=`grep "GPS Week" doyfile | awk '{print substr($7,1,1)}'` igufile="igs${gpsweek}${doyofweek}.sp3" if [ -f $igufile ];then echo "${igufile} is exist!" else wget -nd ftp://cddis.gsfc.nasa.gov/pub/gps/products/$gpsweek/${igufile}.Z uncompress * rm ${igufile}.Z rm doyfile fi done 1.3 下载自己观测文件和IGS站bjfs0940.18o观测文件到rinex目录 site="cusv irkj badg ulab sele chum po12 tash kit3 urum guao lhaz lcki lck2 kunm wuhn xian shao bjfs bjnm chan khaj osn1 daej tnml tcms" for ((line=144;linesitename}${line}0.17o.Z uncompress ${sitename}${ccdoy}0.17.Z sh_crx2rnx -f ${sitename}${ccdoy}0.17d fi done done 解压下载的rinex文件 ```c gunzip *.Z sh_crx2rnx -f *.17d上述步骤1就准备好了所要解算的目录及原始数据准备。 1.4tables更新表文件将所要处理观测文件的表文件根据观测时间进行更新，一般更新的表文件内容包括下述表内容，将这些更新完的表文件，本人建议放置在另外一个文件夹下，因为接下来要将gamit安装目录下的tables文件夹ln -s到text文件目录里，级别与“brdc”，“igs”和“rinex”同级，然后将更新后新的表文件替换至原来旧的tables文件夹中的文件。这样就可以保证gamit安装目录下的tables文件夹不被修改和破坏。 将gamit安装目录下的tables文件夹l拷贝到test目录下，并将新下载更新的表文件复制到tables文件夹下。 2.制作近似坐标文件lfile.lfile.是测站的先验坐标文件 IGS站的三维坐标可以通过IGS官网进行下载snx文件，直接获取，因为是官网公布的，其结果精度是非常可靠的，只不过是根据自己解算的时间，将速度变化量对应时间变化量进行修改，也可以自己提取，见下述步骤。本文提取的部分IGS站坐标如下： BJFS_GPS -2148744.08283 4426641.27265 4044655.92693 -0.03124 -0.00577 -0.00668 2005.000 0.0006 0.0007 0.0007 CHAN_GPS -2674427.27578 3757143.20242 4391521.65453 -0.02518 -0.00889 -0.00920 2006.785 0.0007 0.0007 0.0008 CHUM_GPS 1228950.79178 4508079.91433 4327868.50457 -0.02622 0.00631 0.00291 2005.000 0.0006 0.0008 0.0008 IRKJ_GPS -968328.57783 3794426.52419 5018167.26506 -0.02574 -0.00081 -0.00380 2005.000 0.0006 0.0007 0.0008 LHAZ_GPS -106941.43662 5549269.87743 3139215.00618 -0.04623 -0.00721 0.01407 2005.000 0.0006 0.0008 0.0007 SHAO_GPS -2831733.53698 4675665.93026 3275369.40443 -0.03116 -0.01045 -0.01020 2006.136 0.0007 0.0008 0.0007 TASH_GPS 1695945.20102 4487138.55096 4190140.69855 -0.02652 0.00670 0.00293 2004.620 0.0006 0.0008 0.0007 ULAB_GPS -1257408.54035 4099404.34411 4707992.62361 -0.02884 -0.00161 -0.00499 2005.000 0.0006 0.0007 0.0007 URUM_GPS 193030.63049 4606851.30008 4393311.46943 -0.03082 -0.00156 0.00723 2005.000 0.0006 0.0009 0.0009 所有lfile.提取完后要更新至tables文件夹下。 伪距单点定位提取伪距单点定位提取只是精度第一点，其他没什么影响。 lfile.分布提取 sh_rx2apr -site bjfs1460.17o -nav brdc1460.17n lfile. lfile.批直接提取a）提取观测值.o文件的先验XYZ坐标 grepPOSITION *.10o lfile.rnxb）将.rnx文件转化为.apr文件 rx2aprlfile.rnx 2010 056c）由.apr文件生成lfile.文件 gapr_to_llfile.rnx.apr lfile. 2010 056 载波相位差分提取 sh_rx2apr -site bmss1111.09O -nav brdc1111.09n -ref shao1111.09o -apr itrf05.apr 3.tables文件夹下sestbl.配置sittbl.的配置。nanj/tables下的sittbl.文件对各个测站的先验坐标（或钟差、大气模型等）进行约束。对高精度的已知坐标采取强约束，而对待求点采用松弛约束。如IGS站的坐标分量约束在较小的1～75px，对未知点的约束可以到5～10m。 sites.defaults文件用来控制需要参与解算的测站。在文档末尾可以根据提示编辑，来给定那些不参与解算的测站或是测站的某些天。 process.defaults文件用来控制处理过程中的很多细节，比如sampling interval, number of epochs, start time for processing，default globk .apr file等等，根据需求和提示进行编辑。 在这里均使用它们的默认值。如果更改后进行配置，配置完后要进行链接： ln -s ../tables/process.defaults 5、station.info文件制作将tables文件夹中的station.info文件复制到当前工作目录100中，并删除其中的测站信息，仅保留文件头，并执行下面语句 sh_upd_stnfo -files ../rinex/*.17o打开心生成的station.info文件，查看测站信息是否正确 gedit station.info 注意查看天线高，天线高量取方式等信息。如果有误，修改过了。 6.将tables、igs、brdc和rinex数据链接到工作目录在100目录下将test目录下igs、brdc和rinex数据进行链接 ln -s ../rinex/*.17o ./ ln -s ../igs/*.sp3 ./ ln -s ../brdc/*.17n ./在100目录下将test目录下tables文件数据进行链接 links.day 2017 146 jb80 ln -s ../tables/atl.grid ./atl.grid rm luntab. ln -s ../tables/luntab.2017.J2000 ./luntab. rm soltab. ln -s ../tables/soltab.2017.J2000 ./soltab. rm nutabl. ln -s ../tables/nutabl.2017 ./nutabl. ln -s ../tables/pmu.usno ./pmu. ln -s ../tables/pole.usno ./pole. ln -s ../tables/ut1.usno ./ut1. rm map.grid ln -s ../tables/vmf1grd.2017 ./map.grid ln -s ../tables/otl_FES2004.grid ./otl.grid ln -s ../tables/atl.grid ./atl.grid rm atml.grid ln -s ../tables/atmdisp_cm.2017 ./atml.grid ln -s ../tables/lfile. ./ ln -s ../tables/process.defaults 7.GAMIT解算基线进入test的项目文件夹，输入分布理命令进行解算： makexp sh_sp3fit -f igs19505.sp3 -o igsf -t sh_check_sess -sess 146 -type gfile -file gigsf7.146 makej brdc1460.17n jbrdc7.146 sh_check_sess -sess 146 -type jfile -file jbrdc7.146 makex jb80.makex.batch fixdrv djb807.146 csh bjb807.bat 8.解算结果结果文件将存放在名称为年积日的文件夹内，此例中，可供参考的结果文件为： sh_gamit_2010_006.summary 解算总结 qnanja.006 解算记录 onanja.006 解算记录的简略版，一般关注此文件 hnanja.10006 协防差矩阵、参数平差值 在上述文件中，可根据描述来判断解算结果是否符合相应的需求。其中，基线解算结果（O文件，即这里的onanja.006）中的postfit_nrms项优于0.3左右时最佳；如果大于1.0,则表示此解存在问题。 GAMIT批量处理解算命令：批处理之前，与分布处理一样，要把准备的数据，表文件都进行更新，整理。 1.文件根目录下（例如test）新建 igs/brdc/rinex目录文件，放置相应的文件。 2 文件根目录下（例如test）：sh_setup -yr 2017 -doy 146 链接了table文件，链接之前table文件参数必须全部配置好,更新好，具体更新参考分布处理方法（station.info/sittbl./l-file/sestl./sites.defaults/process.defaults。 3.tables目录（test的子目录）：sh_upd_stnfo -ref station.info -files ../rinex/*.17o 制作station.info文件，所有观测站文件全部放到rinex文件夹下 文件根目录下（例如test）：单天:sh_gamit -expt test -d 2017 146 -pres ELEV -orbit IGSF 多天：sh_gamit -expt test -d 2015 217 218 219 -orbit IGSF sh_gamit -expt test -s 2015 217 219 -orbit IGSF 解算完结果图
解算结果坐标在XYZcors中。 GLOBK批处理平差解算步骤 cors$ sh_glred -d 2017 160 -expt cors -opt H G E cors$ grep Unc. globk_0610_17161.org > XYZ0610 |
【本文地址】
今日新闻 |
推荐新闻 |