. 打开LINUX系统,登陆。(选择安装语言包,关键更新)重启。语言包的设置在system settings-language support- install languages 选择chinese(simplified)点击apply changes 等待系统自动安装,重启,重复步骤进入support 选择系统和菜单语言 汉语 ok(文件夹尽量保持英文状态,以便在终端中输入);
2.GAMIT安装包中没有包含海洋潮汐模型otl.grid 和map——vmf1函数模型需要自己下载连接。 otl.grid在ftp://garner.ucsd.edu/archive/garner/gamit/tables/ 中下载 (约696M) vmf1在ftp://everest.mit.edu//pub/GRIDS/ 中下载(约220M) 如没有下载otl.grid模型可以通过设置禁用该模型。 模型的禁用:
打开/opt/gamit/tables/sestbl. 将use otl.grid的之改为N,表示禁用海洋潮汐模型,同时将 Tides spplied 改为23(原始为31)。
独立下载模型的连接:
进入/opt/gamit/tables 运行终端
运行:ls -al 进行连接 查看otl与map的连接情况 rm otl.grid 移除已有的连接
ln -s /opt/otlmap/otl_FES2004.grid otl.grid 进行连接 (/opt/otlmap/otl_FES2004.grid为模型的路径) rm map.grid
ln -s /opt/otlmap/vmf1grd.2012_063 map.grid (map.grid 目前有vmf提供模型) ls -al 连接并查看是否已连接成功。 连接结束。(连接完成后安装更新时会自动连接)
其它otl模型: ", otl_FES99.grid 182.7 Mb ", otl_CSR4.grid 45.6 Mb ", otl_GOT00.grid 45.6 Mb ", otl_NAO99b.grid 87.1 Mb
Atmosphere模型也许年更新: atml files 185 Mb/year
3.sites.defaults 准备介绍
sites.defaults文件是站处理控制文件。格式如下: all_sites expt xstinfo
brus_gps expt ftprnx ftpraw localrx xstinfo xsite glrepu glreps glts
其中 ftprnx ftpraw localrx xstinfo xsite为gamit选项, glrepu glreps glts为GLOBK选项,格式具体含义如下: expt 为4个字符的项目名称(工程名称)
ftprnx 表示该站点的rinex数据从FTP站点下载 ftpraw 表示该站点的raw数据从FTP站点下载 localrx 表示站点的rinex数据位本地数据
xstinfo 表示该站点不自动更新station.info表中的信息 xsite 该站点的所有天或指定日期的数据不处理 glrepu 表示站点用于计算重复解(默认全部)
glreps 用于定义站点的在GLRED处理中的参考框架(默认为IGS站点列表) glts 在GLRED处理中以时间序列绘制站点图(默认全部)
一般通过自己编制的sites.defaults文件可以设置所有站点不自动更新station.info表中的信息。示例如下: all_sites test xstinfo (test为工程名) LHAS_GPS test localrx glrepu glreps
SHAO_GPS test localrx glrepu
4.station.info文件准备
station.info文件可以采用sh_upd_stnfo 命令生成也可以人工编辑得到。
对于部分不规范的rinex文件建议采用人编辑的方法,并参照gamit的rcvant(天线与接收)文件进行编辑,确保后期处理的顺利。
在更新前应将sites.defaults文件中的 all_sites expt xstinfo注释掉,否则容易报错,无法更新。
采用sh_upd_stnfo 命令生成station.info文件方法:
将tables中的station.info文件拷贝至rinex中,在rinex文件夹下运行sh_upd_stnfo -l sd 提取IGS站点的信息,得到station.info.new
然后见station.info.new改名为station.info,在运行 sh_upd_stnfo -files *12o -expt test 即可。
采用手工编辑的办法:
将tables中的station.info文件拷贝至rinex中,在rinex文件夹下运行sh_upd_stnfo -l sd 提取IGS站点的信息。会自动生成station.info.new 文件。
打开station.info.new删除无关条目,增加其他站点的信息。此时应打开rinex文件进行核对,确保无误后将其保存为station.info。应注意将天线高参考设为DHARP(天线参考中心——一般为天线底部)
格式:
– Fixed format file with station information
? *TRCK SITE Station Name Ant Ht Ant N Ant E \\ ? 7ODM 7ODM Seven Oaks Dam 0.0083 0.0000 0.0000 \\ ??
Rcvr AntCod HtCod Vers \\ ? ASHZ12 ATDMRB DHPAB 9.10 \\ ? Year Doy SN Start Stop ? 2000 110 0 00 00 00 24 00 00 ? (\\ = line break for this display)
– This file must be sorted by TRCK, then epoch
TRCK and SITE – 4-character code
– match the 4-char code in RINEX file name – TRCK = SITE for static analysis 在静态分析中
? Change for kinematic analysis 对于动态分析要改变
– File must be sorted by this field, then epoch ? station_name
– is a long (12-char) description
AntHt
– Raw antenna height as measured 原始测量高度 – Critical measurement! 至关重要 ? HtCod
– Defines geometry of AntHt measurement 定义几何类型
– DHPAB is RINEX standard is vertical height to antenna reference point (ARP) (DHPAB指相对于站点的垂直高度) ? N and E offsets 基准偏移 – e.g., uncalibrated tribrach – rarely used
– occupied point w.r.t. benchmark
RCVR
– 6-character codes for the receiver type – not critical except for some old receivers – or \? Vers
– version number of firmware in the receiver 接收机固件版本号 – Should be correct to handle SNR coding properly
AntCod and HtCod are critical 关键
– Sets phase center model 涉及相位中心模型
– conversion of antenna height measurement 天线高度的解译 – SLBGP -- Slant bottom ground plane 倾斜的底部大地平面
? year and doy specify year and day of year when the values apply
? times set part of day
– (if it changes during UTC day)
? This representation will be augmented in future version of GAMIT
5.sittbl.编辑
打开sittble文件删除无关站点,增加工程站点设置各站点的约束属性。
对于IGS框架点或你需要约束的点设置为0.050 0.050 0.050,其他点100. 100. 100. 一般来说约束点为IGS站点且不少于3个。 示例:
SITE FIX --COORD.CONSTR.-- BJFS BJFS_GPS NNN 100. 100. 100. GUAO GUAO_GPS NNN 100. 100. 100. LHAZ LHAZ_GPS NNN 100. 100. 100. SHAO SHAO_GPS NNN 100. 100. 100. URUM URUM_GPS NNN 0.050 0.050 0.05 WUHN WUHN_GPS NNN 0.050 0.050 0.05 KUNM KUNM_GPS NNN 0.050 0.050 0.05
下载sittbl.refined 文件,里面有最近的各站的可靠性的一个评测。
要想点位置精度高则设置IGS站的约束小一些,即认为IGS站已经很准,强约束,若目的是求基线后面平差则可以将IGS站的约束设置大些。
6.process.defaults 文件可修改部分的解释
process.default文件中需要修改的有采样间隔,最大历元数 ,开始时间 等,示例如下: set sint = '30' 采样间隔依据需要和rinex的情况可修改为15 30 60
set nepc = '2880' 最大历元数 依据观测时长计算调整。
set stime = '0 0' 开始时间(小时 分)依据rinex文件修改,一般要求时间与rinex时间差不大于0.1h 否则报错。 set use_rxc = \设置为N 不依据rinex更新lfile文件。
set aprf = itrf08.apr 采用08框架数据(默认itrf2000) 如不采用GLOBK平差可不设置。
set minxf = '300' Xfile文件大小的最小值,小于该值的将被抛弃。依据历元数设置。历元数少时应减小该值。 set mailto = '***@163.com' 设置为有效邮箱以便接受gamit结算结果。
7.sestbl.文件依据实际需要进行修改。示例如下:
hoice of Experiment = BASELINE (计算高精度基线)
BASELINE/RELAX./ORBIT默认为BASELINE,但一般修改为relax.(BASLINE为固定轨道)(如只需要基线文件可以选择BASELINE)
Choice of Observable = LC_AUTCLN 观测值类型选择
LC_AUTCLN (default), LC_HELP (codeless L2), L1_ONLY (L1 soln from dual freq),L2_ONLY (L2 soln from dual freq),L1,L2_INDEPENDENT (L1 + L2 from dual freq),
L1&L2 (same as L1,L2_INDEPENDENT but with ion constraint);L1_RECEIVER (must add 'L1only' in autcln.cmd) 对于15km以下的短基线应选择L1,L2_INDEPENDENT 或L1_ONLY
Type of analysis =0-iter 迭代多少次[计算结果作为初始值再次计算一般设置1-iter]
Zenith delay =N 天顶延迟 对于短基线选N
Interval zen = 2 ;2 hrs = 13 knots/day (default is 1 ZD per day) 天顶延迟修改为1h,默认2h
Update T/L files = L_ONLY ;T_AND_L (default), T_ONLY, L_ONLY, NONE 人工编辑时可设为NONE
# DMap = GMF ; GMF(default)/NMFH/VMF1 如果已下载VMF1 这可设为VMF1 --重新连接。
# WMap = GMF ; GMF(default)/NMFW/VMF1 如果已下载VMF1 这可设为VMF1 Use map.grid = Y ;VMF1 grid file with mapping functions and ZHD
Inertial frame = J2000 ; B1950/J2000 (default = J2000)
8.sh_gamit 批处理解释
> & ! sh_gamit_test.log 生成过程详解文件sh_gamit_test.log 当采用快速星历处理时需要设置 -orbit IGSU/IGSR
###要处理最近的数据需要将gamit/gamit/fixdrv文件夹下,打开armake.f文件,将2010改为2100.然后再gamit文件夹下运行./install_software.即重新编译安装。(未验证) -mailto lyzlrz@163.com 联网后发送
---------------------------------------------------------------------------------------------------------------------------------------------------
Basic example
? sh_gamit -d 1998 240 -orbit SIOF -expt emed –eops usno
? process GPS data from day 240 of 1998 (Aug 28) ? use any data in the local rinex directory
? ftp data for sites specified by the user ? optionally scan multiple directories to find data
? use Final orbit estimate from Scripps Institution of Oceanography (SIO) 轨道估计数据来源 ? use EOP data from US Naval Observatory (USNO) [IERS Bulletin A] EOP数据来源 ? Name the run “emed\ 工程名 –identify files –select data
List options
– sh_gamit -help
? Select days of data to process
? -d yr days
– where yr and days are the data to be processed e.g., -d 1997 153 156 179 ? -s yr d1 d2
–where d1 is start day and d2 is stop day to be processed e.g. 1997 153 178 ? -r day
–where day is number of days before current date.
sh_gamit options
? -expt
– 4 char name of experiment/solution being run [Default expt] ? -orbit
–Type of orbit to use. [Default SIOF]. Optional second argument gives the day-of-year for the orbit if not the same as processing day ? -orbt
–4 char ID for orbit files. [Default 4 char extracted from the orbit selected] ? -eops
– Name of EOP series to use. (bull_b, bull_a, etc.) [Default usno]
-sessinfo
–Session - sampling interval, number of epochs, start time (default 30 2880 0 0) ? -rx_doy_plus n
–Search RINEX file names this many days into the future for current-day data 寻找往前推迟n天的rinex数据 ? -rx_doy_minus n
– Search RINEX file names this many days into the past for current-day data 寻找往后n天的rinex数据 ? -xsite list
–List of sites to be excluded from the processing. 解算排除的站
? -netext
–add network suffix
? add year prefix to day directories. Ie 1997_053 日期文件夹添加年前缀 ? -pres
? Plot gamit residuals as skyplots
? Do NOT create skyplot GIF graphics files 不产生天空GIF图片