gamit-工程处理流程－详细.txt
工程说明：工程名为：db05。数据来自中国大陆7个IGS站数据，其中kunm、urum、lhaz三个作为已知站，其他4站有：wuhn,shao,bjfs,guao,计算日期从2009年5月1日到2009年5月4日。年纪日（121～124），以下为具体操作流程：
1、在工程文件夹db05下，首先建立tables文件夹的链接，命令如下。完成后一定要检查utl. pole. soltbl. luntab. nutabl.等文件，主要是处理的数据日期是否在这些文件有效范围内，否则到网站http://sopac.ucsd.edu/processing/gamit/ 下载
sh_setup -yr 2009
2、在工程文件夹下建立各种工作必须目录，命令如下：
mkdir rinex igs brdc
3、在rinex文件夹下，运行下面命令，下载4天的IGS的观测文件，主要是中国大陆的IGS站，共用7站，其中kunm、urum、lhaz3站作为已知站
sh_get_rinex -archive sopac -yr 2009 -doy 121 -ndays 4 -sites kunm urum lhaz khaj wuhn shao bjfs guao
Information extracted from ftp_info
########################################
ftpsite garner.ucsd.edu
ftplogin anonymous zhao@linux
ftpdir /pub/rinex/YYYY/DDD
ftpcmd ftp -inv garner.ucsd.edu
wgetsite garner.ucsd.edu
wlogin anonymous zhao@linux
########################################
directory = /pub/rinex/2009/121
Try 1: Getting rinex data for day 121 from sopac archive, garner.ucsd.edu
Attempting to download sites: kunm urum lhaz khaj wuhn shao bjfs guao
ftp -inv garner.ucsd.edu
Successfully downloaded: kunm1210.09d.Z urum1210.09d.Z lhaz1210.09d.Z khaj1210.09d.Z wuhn1210.09d.Z shao1210.09d.Z bjfs1210.09d.Z guao1210.09d.Z
directory = /pub/rinex/2009/122
Try 1: Getting rinex data for day 122 from sopac archive, garner.ucsd.edu
Attempting to download sites: kunm urum lhaz khaj wuhn shao bjfs guao
ftp -inv garner.ucsd.edu
Successfully downloaded: kunm1220.09d.Z urum1220.09d.Z lhaz1220.09d.Z khaj1220.09d.Z wuhn1220.09d.Z shao1220.09d.Z bjfs1220.09d.Z guao1220.09d.Z
directory = /pub/rinex/2009/123
Try 1: Getting rinex data for day 123 from sopac archive, garner.ucsd.edu
Attempting to download sites: kunm urum lhaz khaj wuhn shao bjfs guao
ftp -inv garner.ucsd.edu
Successfully downloaded: kunm1220.09d.Z kunm1230.09d.Z urum1230.09d.Z lhaz1230.09d.Z khaj1230.09d.Z wuhn1230.09d.Z shao1230.09d.Z bjfs1230.09d.Z guao1230.09d.Z
directory = /pub/rinex/2009/124
Try 1: Getting rinex data for day 124 from sopac archive, garner.ucsd.edu
Attempting to download sites: kunm urum lhaz khaj wuhn shao bjfs guao
ftp -inv garner.ucsd.edu
Successfully downloaded: kunm1240.09d.Z urum1240.09d.Z lhaz1240.09d.Z khaj1240.09d.Z wuhn1240.09d.Z shao1240.09d.Z bjfs1240.09d.Z guao1240.09d.Z
4、在SP3文件夹下运行如下命令，获取4天精密星历文件。
sh_get_orbits -archive sopac -yr 2009 -doy 121 -ndays 4 -makeg no
ftp_info = sh_get_ftp_info -archive sopac -type sp3
Information extracted from ftp_info
########################################
ftpsite garner.ucsd.edu
ftplogin anonymous zhao@linux
ftpdir /pub/products/GPSW
ftpcmd ftp -inv garner.ucsd.edu
########################################
Getting requested orbit files
5、在rinex文件夹下运行如下命令，获取4天的广播星历，这个星历包含所有卫星，文件以brdc开头，而不是测站出观测到的广播星历。
sh_get_nav -archive sopac -yr 2009 -doy 121 -ndays 4
Information extracted from ftp_info
########################################
ftpsite garner.ucsd.edu
ftplogin anonymous zhao@linux
ftpdir /pub/rinex/YYYY/DDD
ftpcmd ftp -inv garner.ucsd.edu
########################################
Getting auto or brdc files for 2009 121 for 4 days from garner.ucsd.edu
6、准备lfile.文件。这里分为两部分，已知约束点坐标（IGS）从tables目录下的itrf00.apr文件中拷贝，求解的未知点利用观测o文件手动生成（脚本如下），
制作批处理脚本genLfile(注意修改年份，这里是09o)放在rinex文件夹下，
#!/bin/bash
for file in *09o
do sh_rx2apr -site $file
done
cat *.apr >>./lfile.
rm *.apr
把上面文件保存在rinex文件夹下，执行方法如下：
（1）、修改可执行行：chmod +x genLfile
执行：./genLfile
或者利用如下命令：
（2）、sh genLfile
制作好的lfile. 如下替换掉tables文件夹中的相应文件：
BJFS_GPS -2148743.784 4426641.236 4044655.935 -0.0444 0.0141 -0.0013 1997.0 21601M001 BEIJING .014 .019 .017 .0040 .0055 .0049
GUAO_GPS 228380.790000 4631962.420000 4367041.080000 0.0000 0.0000 0.0000 2009.34
KUNM_GPS -1281255.473 5640746.079 2682880.117 -0.0318 -0.0024 -0.0203 1997.0 21609M001 KUNMING .010 .022 .012 .0037 .0083 .0043
LHAS_GPS -106937.669 5549269.591 3139215.762 -0.0462 -0.0076 0.0121 1997.0 21613M001 LHASA .001 .003 .002 .0005 .0011 .0007
SHAO_GPS -2831733.583000 4675665.958000 3275369.410000 0.0000 0.0000 0.0000 2009.34
URUM_GPS 193030.873 4606851.324 4393311.421 -0.0306 -0.0055 -0.0003 1997.0 21612M001 URUMQI .003 .006 .006 .0009 .0022 .0020
WUHN_GPS -2267749.162 5009154.325 3221290.762 -0.0325 -0.0077 -0.0119 1997.0 21602M001 WUHAN .002 .003 .002 .0007 .0012 .0009
7、准备sites.defaults 如下：（kunm,urum,lhaz做为已知点）:xstinfo表示不更新station.info文件，glrepu表示重复解，glreps表示作为参考框架
BJFS_GPS db05 ftprnx glrepu
GUAO_GPS db05 ftprnx glrepu
KUNM_GPS db05 ftprnx xstinfo glrepu glreps
LHAS_GPS db05 ftprnx xstinfo glrepu glreps
SHAO_GPS db05 ftprnx glrepu
URUM_GPS db05 ftprnx xstinfo glrepu glreps
WUHN_GPS db05 ftprnx glrepu
8、制作station.info文件，利用在rinex文件夹下用命令手动生成，命令如下：
sh_upd_stnfo -files *.09o -expt db05
最终如下：
BJFS BJFS 2009 121 0 0 30 2009 121 23 59 30 0.0460 DHARP 0.0000 0.0000 ASHZ12 9.20 ATDMCB mstinf2: bjfs1210.09o
BJFS BJFS 2009 122 0 0 30 2009 122 23 59 30 0.0460 DHARP 0.0000 0.0000 ASHZ12 9.20 ATDMCB mstinf2: bjfs1220.09o
BJFS BJFS 2009 123 0 0 30 2009 123 23 59 30 0.0460 DHARP 0.0000 0.0000 ASHZ12 9.20 ATDMCB mstinf2: bjfs1230.09o
BJFS BJFS 2009 124 0 0 30 2009 124 23 59 30 0.0460 DHARP 0.0000 0.0000 ASHZ12 9.20 ATDMCB mstinf2: bjfs1240.09o
.
GUAO GUAO 2009 121 0 0 0 2009 121 23 59 30 0.0800 DHARP 0.0000 0.0000 ATUZ12 9.93 ATDMRB mstinf2: guao1210.09o
GUAO GUAO 2009 122 0 0 0 2009 122 23 59 30 0.0800 DHARP 0.0000 0.0000 ATUZ12 9.93 ATDMRB mstinf2: guao1220.09o
GUAO GUAO 2009 123 0 0 0 2009 123 23 59 30 0.0800 DHARP 0.0000 0.0000 ATUZ12 9.93 ATDMRB mstinf2: guao1230.09o
GUAO GUAO 2009 124 0 0 0 2009 124 23 59 30 0.0800 DHARP 0.0000 0.0000 ATUZ12 9.93 ATDMRB mstinf2: guao1240.09o
.
KUNM KUNM 2009 121 0 0 0 2009 121 23 59 30 0.0793 DHARP 0.0000 0.0000 TR8000 0.00 TRBROG mstinf2: kunm1210.09o?
KUNM KUNM 2009 122 0 0 0 2009 122 23 59 30 0.0793 DHARP 0.0000 0.0000 TR8000 0.00 TRBROG mstinf2: kunm1220.09o?
KUNM KUNM 2009 123 0 0 0 2009 123 23 59 30 0.0793 DHARP 0.0000 0.0000 TR8000 0.00 TRBROG mstinf2: kunm1230.09o?
KUNM KUNM 2009 124 0 0 0 2009 124 23 59 30 0.0793 DHARP 0.0000 0.0000 TR8000 0.00 TRBROG mstinf2: kunm1240.09o?
.
LHAZ LHAZ 2009 121 0 38 0 2009 121 23 59 30 0.1330 DHARP 0.0000 0.0000 TRMSST 0.00 ATGGDB mstinf2: lhaz1210.09o??
LHAZ LHAZ 2009 122 0 38 30 2009 122 23 59 30 0.1330 DHARP 0.0000 0.0000 TRMSST 0.00 ATGGDB mstinf2: lhaz1220.09o??
.LHAZ LHAZ 2009 123 0 38 30 2009 123 23 59 30 0.1330 DHARP 0.0000 0.0000 TRMSST 0.00 ATGGDB mstinf2: lhaz1230.09o??
LHAZ LHAZ 2009 124 0 38 30 2009 124 23 59 30 0.1330 DHARP 0.0000 0.0000 TRMSST 0.00 ATGGDB mstinf2: lhaz1240.09o??
SHAO SHAO 2009 121 0 0 0 2009 121 23 59 30 0.0814 DHARP 0.0000 0.0000 ATUZ12 0.00 TRBROG mstinf2: shao1210.09o?
SHAO SHAO 2009 122 0 0 0 2009 122 23 59 30 0.0814 DHARP 0.0000 0.0000 ATUZ12 0.00 TRBROG mstinf2: shao1220.09o?
SHAO SHAO 2009 123 0 0 0 2009 123 23 59 30 0.0814 DHARP 0.0000 0.0000 ATUZ12 0.00 TRBROG mstinf2: shao1230.09o?
SHAO SHAO 2009 124 0 0 0 2009 124 23 59 30 0.0814 DHARP 0.0000 0.0000 ATUZ12 0.00 TRBROG mstinf2: shao1240.09o?
.
URUM URUM 2009 121 0 0 0 2009 121 23 59 30 0.0460 DHARP 0.0000 0.0000 TRMSST 3.30 TPSC3D mstinf2: urum1210.09o?
URUM URUM 2009 122 0 0 0 2009 122 23 59 30 0.0460 DHARP 0.0000 0.0000 TRMSST 3.30 TPSC3D mstinf2: urum1220.09o?
URUM URUM 2009 123 0 0 0 2009 123 23 59 30 0.0460 DHARP 0.0000 0.0000 TRMSST 3.30 TPSC3D mstinf2: urum1230.09o?
URUM URUM 2009 124 0 0 0 2009 124 23 59 30 0.0460 DHARP 0.0000 0.0000 TRMSST 3.30 TPSC3D mstinf2: urum1240.09o?
.
WUHN WUHN 2009 121 0 0 30 2009 121 23 59 30 2.3610 DHARP -0.0094 -0.0022 ASHZ12 9.20 ATDMCE mstinf2: wuhn1210.09o
WUHN WUHN 2009 122 0 0 30 2009 122 23 59 30 2.3610 DHARP -0.0094 -0.0022 ASHZ12 9.20 ATDMCE mstinf2: wuhn1220.09o
WUHN WUHN 2009 123 0 0 30 2009 123 23 59 30 2.3610 DHARP -0.0094 -0.0022 ASHZ12 9.20 ATDMCE mstinf2: wuhn1230.09o
WUHN WUHN 2009 124 0 0 30 2009 124 23 59 30 2.3610 DHARP -0.0094 -0.0022 ASHZ12 9.20 ATDMCE mstinf2: wuhn1240.09o
9、利用gamit批处理模式解算基线，在工程目录下执行如下命令：
sh_gamit -d 2009 121 122 123 124 -expt db05 -noftp > & ! sh_gamit_db05.log
可用gedit打开 sh_gamit_db05.log查看进行到了那里，是否出了错误等。
运行完毕后查看各年纪日文件夹下的sh_gamit_***.summary、qdb05a.***等文件，确定是否解算成功、精度符合规范等。此工程 Postfit nrms均在0.2附近，说明解算成功。
其中一个总结文件sh_gamit_121.summary如下：
Input options -d 2009 121 122 123 124 -expt db05 -noftp
Processing 2009 121 GPS week 1529 5 Using node: linux-zhao Started at: 09_06_01_16:09:35
Processing directory: /home/zhao/project/db05/121
Disk Usage: 82428.7 Free 17969.6 Mbyte. Used 83%
Number of stations used 7 Total xfiles 7
X-files below minxf (300KB) threshold: 法访问 法访问 法访问 法访问 法访问 法访问 法访问
Sites excluded by xsite command
Postfit RMS rms, total and by satellite
RMS IT Site All 01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32
RMS 16 ALL 8.0 69 66 113 69 60 84 85 76 73 66 88 89 78 87 66 82 70 68 77 95 74 82 94 73 77 69 73 78 76 82 88 90
Best and Worst two sites:
RMS 16 GUAO 6.9 7 5 8 7 6 8 6 8 6 6 8 9 8 6 6 6 6 5 6 8 7 7 8 6 5 6 7 6 8 6 7 8
RMS 16 SHAO 6.9 0 6 7 6 0 6 7 7 7 7 6 8 8 8 5 8 7 7 6 8 5 6 6 8 7 8 5 7 7 8 7 8
RMS 16 WUHN 8.6 0 6 10 8 0 8 8 9 10 8 7 10 9 9 7 9 8 8 8 11 7 9 9 8 7 8 8 8 9 10 11 9
RMS 16 KUNM 10.8 8 9 22 7 6 7 15 8 7 6 14 8 6 12 8 11 8 6 13 15 8 9 16 10 10 7 12 11 6 12 10 12
Double difference statistics
Prefit nrms: 0.15731E+01 Postfit nrms: 0.19834E+00
Prefit nrms: 0.15688E+01 Postfit nrms: 0.20417E+00
Prefit nrms: 0.15731E+01 Postfit nrms: 0.18983E+00
Prefit nrms: 0.15688E+01 Postfit nrms: 0.19448E+00
Number of double differences: 27955
Phase ambiguities (Total WL-fixed NL-fixed): 180 171 151 AUTCLN (Total Expected) 180 180
Phase ambiguities WL fixed 95.0% NL fixed 83.9%
Processing 2009 121 GPS week 1529 5 Using node: linux-zhao Finished at: 09_06_01_16:14:55
基线解算就此完成。平差可以采用其他软件（如武大的PowerAdj）或采用gamit相应的globk。本工程采用globk平差。过程如下：
10、利用如下脚本生成要求的的文件夹链接，在工程目录下执行这个脚本，方法同前：
for file in 1?? # this is for current directory
do
ln -s $file 2009_$file
done
11、运行如下命令，实现批处理平差模式：
sh_glred -s 2009 121 2009 124 -yrext -local db05 -opt H G E
因为上面命令没有指定工程名，因此采用了默认的expr工程名。运行结束后查看gsoln文件夹下的sum.expr文件，详细结果在globk_expt_09***.org等4个文件中。
其中sum.expr 如下：
ENSUM: GLOBK NEU component summary: includes the following solutions
File tmp.ensum_input: Combination of ALL networks
Stations enu S# # mean length sig wrms nrms slope sig wrms nrms dur mean line #
(m) (m) (mm) (mm/yr) (mm/yr) (mm) (yrs) (yrs)
SHAO_GPS- N 1 4 3461996.3402 0.0355 18.3 0.26 1418.45 11634.00 21.5 0.30 0.01 2009.33 1- 7
SHAO_GPS- E 1 4 11552738.9544 1.2501 225.6 0.09 19993.08 ******** 262.4 0.10 0.01 2009.33 12- 18
SHAO_GPS- U 1 4 21.5172 0.4081 202.6 0.25 ******** ******** 238.5 0.29 0.01 2009.33 23- 29
BJFS_GPS- N 1 4 4409209.3915 0.1631 52.8 0.16 5482.05 53351.20 60.1 0.18 0.01 2009.33 34- 40
BJFS_GPS- E 1 4 9939262.4549 1.2346 195.4 0.08 18376.46 ******** 225.8 0.09 0.01 2009.33 45- 51
BJFS_GPS- U 1 4 86.9263 0.4167 169.2 0.20 ******** ******** 189.5 0.23 0.01 2009.33 56- 62
WUHN_GPS- N 1 4 3398768.0109 0.1700 23.1 0.07 2624.29 55538.72 25.9 0.08 0.01 2009.33 67- 73
WUHN_GPS- E 1 4 10965014.4971 1.2472 198.7 0.08 18470.75 ******** 229.8 0.09 0.01 2009.33 78- 84
WUHN_GPS- U 1 4 25.3512 0.4111 235.7 0.29 ******** ******** 276.5 0.34 0.01 2009.33 89- 95
KUNM_GPS- N 1 4 2786275.3981 0.4423 17.2 0.02 3132.77 ******** 16.1 0.02 0.01 2009.33 100- 106
KUNM_GPS- E 1 4 10368693.2207 1.2557 144.6 0.06 14968.96 ******** 164.8 0.07 0.01 2009.33 111- 117
KUNM_GPS- U 1 4 1985.7652 0.4130 310.3 0.38 ******** ******** 367.3 0.44 0.01 2009.33 122- 128
LHAZ_GPS- N 1 4 3301439.1381 0.7038 59.8 0.04 5547.11 ******** 69.2 0.05 0.01 2009.33 133- 139
LHAZ_GPS- E 1 4 8813208.7511 1.2231 86.0 0.04 8353.36 ******** 99.0 0.04 0.01 2009.33 144- 150
LHAZ_GPS- U 1 4 3624.1771 0.4132 305.0 0.37 ******** ******** 353.8 0.43 0.01 2009.33 155- 161
URUM_GPS- N 1 4 4876678.6982 0.7840 124.8 0.08 12312.75 ******** 143.2 0.09 0.01 2009.33 166- 172
URUM_GPS- E 1 4 7037402.0994 1.1504 98.7 0.04 4520.03 ******** 119.3 0.05 0.01 2009.33 177- 183
URUM_GPS- U 1 4 858.4000 0.4196 196.7 0.23 ******** ******** 206.9 0.25 0.01 2009.33 188- 194
GUAO_GPS- N 1 4 4839181.6651 0.7923 124.1 0.08 12314.34 ******** 142.3 0.09 0.01 2009.33 199- 205
GUAO_GPS- E 1 4 7042905.1263 1.1501 96.0 0.04 4482.59 ******** 115.9 0.05 0.01 2009.33 210- 216
GUAO_GPS- U 1 4 2028.2399 0.4193 198.1 0.24 ******** ******** 207.0 0.25 0.01 2009.33 221- 227
对比各站的IGS坐标，可见符合要求！
至此整个工程的解算过程就完成了，当然这里只是示范性的解算，有很多参数可以在process.defaults、sestbl.、sittbl.中具体设置，而不采用默认值。
gamit\globk软件虽然采用可以批处理模式，使得相对简化了处理的复杂性，但只用在不断的反复解算中才能掌握各种处理的技巧、问题解决的方法！
爱情
(1)
安装
(2)
北斗
(1)
毕业
(1)
变量
(1)
测绘局
(1)
测量
(2)
插件
(1)
查询
(1)
常用
(1)
成果转化
(1)
词汇
(1)
慈善
(1)
答辩
(2)
代码
(2)
电台
(1)
发泄
(1)
感悟
(4)
高程
(1)
搞笑
(5)
共产党
(1)
古诗词
(1)
管理
(1)
函数
(3)
绘图
(1)
加密
(1)
交际
(1)
教程
(4)
教育
(2)
解决
(4)
解密
(1)
精度
(1)
酒桌
(2)
开源
(2)
科技
(1)
科学
(1)
刻录
(1)
老外
(1)
励志
(5)
连续剧
(1)
恋爱
(1)
列表
(1)
领导
(1)
美食
(2)
名人
(3)
命令
(4)
区别
(1)
日记
(2)
软件
(12)
商业
(3)
时政
(1)
视频
(1)
数据
(1)
算法
(1)
投影
(2)
图论
(1)
网络
(1)
网站
(1)
卫星
(3)
未成年
(1)
慰问
(1)
文本
(1)
文件
(2)
下载
(4)
笑话
(2)
学习
(9)
遥感
(1)
疑问
(5)
营销
(1)
娱乐
(3)
源代码
(2)
政策
(1)
指导
(3)
智慧
(5)
主成分分析
(2)
抓图
(1)
专家
(1)
资料
(5)
字符串
(1)
最短路径
(1)
坐标
(1)
baidu
(3)
Bernese
(1)
blog
(2)
c#
(9)
China
(3)
Dijkstra
(1)
DNS
(1)
doris
(1)
dos
(1)
excel
(1)
firefox
(3)
GAMIT
(8)
gcc
(1)
GIS
(3)
GMT
(1)
GPS
(5)
ITRF
(1)
linux
(5)
mapx
(1)
matlab
(6)
movie
(4)
music
(3)
oracle
(2)
pic
(1)
PPT
(3)
PROJ.4
(2)
python
(2)
QQ
(2)
rinex
(1)
shell
(2)
sql
(1)
teqc
(3)
tools
(1)
tps
(1)
ubuntu
(5)
USA
(1)
website
(1)
博客归档
-
▼
2010
(118)
-
▼
六月
(33)
- psvelomeca.c ----c语言源代码
- QQ默认的表情 恢复
- 5 分钟平静你的内心：佛家箴言130条
- 包菜炒辣椒
- excel 下拉菜单怎么做，可以进行选择数据
- 番茄炒包菜
- 科学家眼中的世界--图2
- how to 10个网站（教你怎么做）
- 两篇“另类”高考作文
- 卖苹果
- 人生三大陷阱----大意、轻信、贪婪。
- 这些话好帅，经典得让人心痛。。。
- oracle 监听程序 启动
- gamit-工程处理流程－详细
- oracle常用sql语句
- bernese 错误.log
- 笑话----几则
- 答辩通知
- 笑话-考试--斗志
- 期待看的电影---忠奸人
- 把文件刻录为镜像文件iso---软件nero
- 跑步时在想什么？
- 全国门类、学科、专业代码查询
- 完美新娘---烂片一部
- 如何变成交际大师
- 不喝酒的后果如此严重
- 我国成功发射第四颗北斗导航卫星
- firefox 关闭安全更新
- 使用 Scrapbook作为个人资料管理系统
- FWTools 是一套处理GIS和遥感数据的开源工具
- Proj.4 简明文档 中文
- 计算椭球上两点间的距离 geod
- 在ubuntu上安装proj 投影转换程序
-
▼
六月
(33)
订阅：
博文评论 (Atom)
没有评论:
发表评论