gamit学习总结 - 图文

2020-02-21 21:05

. 打开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 to day directories. i.e. 035r 日期文件夹添加网的后缀 ? -yrext

? add year prefix to day directories. Ie 1997_053 日期文件夹添加年前缀 ? -pres

? Plot gamit residuals as skyplots . YES, NO, ELEV. [Default N] 绘制天空残差图 ? -nogifs

? Do NOT create skyplot GIF graphics files 不产生天空GIF图片


gamit学习总结 - 图文.doc 将本文的Word文档下载到电脑 下载失败或者文档不完整,请联系客服人员解决!

下一篇:职业经理人资质证书培训文件

相关阅读
本类排行
× 注册会员免费下载(下载后可以自由复制和排版)

马上注册会员

注:下载文档有可能“只有目录或者内容不全”等情况,请下载之前注意辨别,如果您已付费且无法下载或内容有问题,请联系我们协助你处理。
微信: QQ: