gamit学习总结 - 图文

. 打开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图片

-copt [file types]

–Compress remaining files for each dayserve when d [Default x k autcln.out] 压缩指定类型的每天文件 ? -aopt [file types]

–Copy the specified files for each day to the directory specified by –archivepth 将指定每天文件复制到archivepth文件夹。

? -dopt [ file types]

–Delete specified files after processing each day. [Default c] 删除指定类型文件

– -mailto

– Name of user and machine to send results summary. (e.g., simon@wegener.mit.edu) -snrres

–Calculate theoretical phase residuals from SNR. Requires raw data or RINEX with valid SNR values. ? -noacc

–Do not exclude satellites base on accuracy codes in the .sp3 file. [Default exclude] 不排除.sp3中基于精度码方式排除卫星。 ? -xver char

– 1-character x-file version (6th character of x-file). [Default is single-digit year] ? -noftp

–Use this option if no ftp download connections are available. 不联网

? -remakex Y/C/N

? Y - Remake x-files without checking.

? C - Remake x-file if orbit and session.info are incompatible and all RINEX files are available [Default C] ? N - Don't remake x-files under any circumstances

若tables中station.info消失，sh_gamit将会使用rinex的头文件信息。

guess_rcvant.dat 来将rinex的头文件信息翻译成gamit 可阅读形式。 Antenna height errors 在GLOBK可降低，但类型识别错误无法消除。

---------------------------------------------------------------------------------------------------------------------------------------------------

9.打开日期目录查看结算结果 判断解算质量主要有两个方面： 一是数据利用率

（即是否有足够多的可用于进行合理估计的数据参与计算）； 二是数据的噪声水平是否与选用的模型匹配。

（这一项主要有两个指标，基线分量误差的大小和 nrms（标准化的均方根误差)，对于基线分量误差的判断主要依据预先设定的站点坐标和轨道参数的先验约束，如果大于先验可以通过查看Q文件或者 autcln.sum 文件，是否有大量被autcln丢弃的数据。对于nrms一般应在0.25左右。如果大于0.5就意味着周跳没有消除，或出现了模型错误（固定站点坐标错误或者未模型化的卫星）等。解算结果文件autcln.sum.post中的rms残差为显示控制网站点的相对质量。

打开 sh_gamit_doy.summary，otesta.doy, qtesta.doy 文件查看解算总结和基线解算情况，解算成功与否判断： Postfit nrms: *****E+00 是否在0.25左右， 最差不能大于0.5，基线分量的误差大小是否小于预设值。 至此基线解算全部结束。

基线解算结果：sh_gamit_doy.summary

基线解算结果： qexpta.doy,oexpta.doy,hexpta.yrdoy

10.若软件限制了最大测站数据，最大卫星数目，最大对流层延迟参数，最大解算历元数目。更改于 Makefile.config文件。然后重新编。

11.tables表介绍。

12.格式问题

o文件名需要小写，并且内容格式符合UINX 文件名改小写：sh_casefold -dir d -files * 格式转UNIX: sh_dos2unix *

13.gamit模块以及结构

14.

Data Types

? L1 and L2 phase

– used for parameter estimation ? P1 and P2 range measurements – used to help clean phase data

? used to estimate receiver clock offset with tolerance of 1 microsecond (= 300 m)

? LC phase

– proportional to geometric range

– sensitive insensitive to ionospheric delay – linear combination of L1 and L2 ? LG phase

– measures ionospheric delay only – insensitive to geomery ? Widelane WL

– linear combination of L1, L2, P1 and P2

– Melbourne-Wubbena

– Should be constant (+/- noise)

15. Aim of GAMIT

Basic aim of GAMIT is to generate positions/orbits for one day of data. Standard product of GAMIT is a minimally constrained solution

16. gg/tables/rcvant.dat

File contains all receivers and antennas known to GAMIT

? Translates user-specified 6-character code to IGS SINEX standard.

? The file in templates directory provides mapping from RINEX header information to 6-char codes. ? Example lines

–rcv ASHTECH^Z-XII ASHZ12 –rcv Z-XII3 ASHZ12 –swv 5F 1.00

? Explanation given at top of file ? ^ used a represent space

? Sometimes need to be careful that pattern correctly gets correct type. Should not be a problem for IGS rinex files

? default lines allowed, but be careful

17. 模型说明：

1. 2. 3. 4.

otl 潮汐改正

vmfl GMF 投影函数

atml大气荷载模型，对高程影响较大，可消除周跳波动，可靠性需要进一步证实 atl大气抄袭荷载模型和met气象模型

下载地址：ftp://everest.mit.edu

1) 海洋潮汐。例如otl_FES2004.grid放在软件talbels目录下。链接到otl.grid。 2) 大气负荷。例如atmldisp_cm.2006，每年更新一次。连接到atmldisp_grid.2006 3) vmfl投影函数栅格，例如vmflgrd2006，连接到map.grid.2006。每年更新一次。

18. 解算所需文件说明

类型文件名说明： 数据文件

bjfs0900.06o bjfs 为站点代码， 090 为年积日，0为时段号， 06 为2006 年，o 为观测文件 auto0900． 06n 广播星历文件，前4 个字节改为站点代码则为站点广播星历文件 igs13685． sp3 igs 为事后精密星历，1368 为GPS周，5为星期五 igr13685． sp3 igr 代表快速精密星历 igu13685． sp3 igu 代表预报精密星历

配置文件

lgigg6． 090 l 代表L 文件， gigg 代表工程名，6代表2006 年， 090 为年积日 station． Oct session． Info sestbl．

测站信息文件，记录各测站接收机和天线类型 时段信息文件，可以在数据处理过程中自动创建

解算参数控制文件，设置与解算有关的参数，达到不同的解算目的

sittbl． 表文件

设置测站坐标约束、大气模型等参数

ut1． 国际时间系统表 pole． 极移表 leap． sec 跳秒表 nutabl． 章动表 luntal． 月亮表 soltab． 太阳表

gdetic． dat 地球形状参数表

rcvant． dat 接收机及其天线型号对照表 svnav． dat 星号对照表

svs_exclude．Dat 坏卫星信息表 hi． dat 仪器高参数表 dcb． dat 差分码偏差表

autcln． cmd 周跳自动探测和修改命令表 antmod． dat 天线相位中心改正表

ot1． list /grid 全球格网或各测站的海洋潮汐参数表

19 .crx2rnx 压缩的d文件转换成o文件

Usage: crx2rnx [file] [-] [-f] [-s] [-h]

stdin and stdout are used if input file name is not given. - : output to stdout

-f : force overwrite of output file

-s : skip strange epochs (default:stop with error)

This option may be used for salvaging usable data when middle of the Compact RINEX file is missing. The data after the missing part, are, however, useless until the compression operation of all data are initialized at some epoch. Combination with use of -e option of RNX2CRX (ver.4.0 or after) may be effective.

Caution : It is assumed that no change in # of data types happens in the lost part of the data. -h : display help message

exit code = 0 (success) = 1 (error) = 2 (warning)

20．某个人的参数设置： 表2 解算控制参数的选择 控制选项控制参数 处理模式RELAX． 迭代方案1-ITER 观测量LC_AUTCLN

天顶延迟参数25 Number Zen = 25 ; number of zenith-delay parameters (default 1) 梯度参数3

--------------------Number gradients

海潮模型otl.list

光压模型BERNE 惯性框架J2000． 0 对流层误差模型SAAS

截止高度角 /( °) 10 ---------- Elevation Cutoff 测站约束( BJFS) /m 0． 05 /0． 05 /0． 1 测站约束( WUHN) /m 0． 025 /0． 01 /0． 1 测站约束( KUNM) /m 9． 999 /9． 999 /9． 999 9 其他默认

某篇论文写到;

不同截止高度角对基线精度的影响

表4 不同截止高度角对基线(BJFS_KUNM)精度的影响 精度指标 nrms σX /m σY /m σZ /m

高度角/( °)

20

25

± 0． 213 98 ± 0． 051 2 ± 0． 056 3 ± 0． 055 1

± 0． 222 52 ± 0． 048 5 ± 0． 049 1 ± 0． 051 1

5

10

15

± 0． 241 04 ± 0． 041 9 ± 0． 040 1 ± 0． 044 5

± 0． 241 04 ± 0． 041 9 ± 0． 040 1 ± 0． 044 5

± 0． 231 09 ± 0． 045 0 ± 0． 043 4 ± 0． 047 8

σL /L 4.90 × 10 E-9 4． 90 × 10 E-9 5． 49 × 10 E-9 6． 53 × 10 E-9 7． 67 × 10 E-9

不过高度角太低时路径延迟对水汽变化非常敏感，因此在一般的解算情况下选择10°的截止高 度角比较适宜。

不同天顶延迟参数对基线精度的影响

表5 不同天顶延迟参数对基线(BJFS_KUNM)精度的影响 精度指标 Nrms σX /m σY /m σZ /m σL /L

天顶延迟参数

25

± 0． 241 04 ± 0． 041 9 ± 0． 040 1 ± 0． 044 5

7

13

19

± 0． 239 47 ± 0． 042 5 ± 0． 040 8 ± 0． 045 9

± 0． 242 77 ± 0． 041 9 ± 0． 040 2 ± 0． 044 7

± 0． 241 66 ± 0． 041 8 ± 0． 039 8 ± 0． 044 2

4． 95 × 10 E-9 5． 05 × 10 E-9 4． 90 × 10 E-9 4． 90 × 10 E-9

一般情况下，天顶延迟参数取为13，即每2 h 产生一个天顶延迟估值是比较合适的。天顶延迟参数变化对长基线的精度影响不是非常明显。

选取不同观测量对基线精度的影响

表6 选取不同观测量对基线(BJFS_KUNM)精度的影响 精度指标 nrms

观测量

LC_HELP

LC_ONLY ± 0． 239 27 ± 0． 047 9 ± 0． 044 6 ± 0． 055 3

± 0． 239 27 ± 0． 047 9 ± 0． 044 6 ± 0． 055 3

LC_AUTCLN

σX /m σY /m σZ /m σL /L

不同测站坐标约束对基线(BJFS_KUNM)精度的影响 精度指标

测站约束/m

0.1/0.1/0.1

0.5/0.5/0.5

1/1/1

0.01/0.01/0.02

0.05/0.05/0.06

L1＆L2

± 0． 241 04 ± 0． 041 9 ± 0． 040 1 ± 0． 044 5

± 0． 239 28 ± 0． 047 3 ± 0． 044 1 ± 0． 053 8

4． 90 × 10 E-9 6． 53 × 10 E-9 6． 73 × 10 E-9 6． 73 × 10 E-9

LC_AUTCLN( 观测值是无电离层的线性组合)

nrms ± 0． 239 38 ± 0． 241 18 ± 0． 037 2 ± 0． 034 8 ± 0． 036 0

± 0． 241 07 ± 0． 047 4 ± 0． 040 1 ± 0． 046 1

± 0． 240 64 ± 0． 119 4 ± 0． 081 1 ± 0． 119 9

± 0． 240 59 ± 0． 223 9 ± 0． 146 2 ± 0． 225 9

σX /m ± 0． 013 3 σY /m ± 0． 019 9 σZ /m ± 0． 015 3 σL /L

21. 使用观测值类型的选择： 测站约束越严格，精度越高。

3． 66 × 10 E-9 4． 41 × 10 E-9 4． 75 × 10 E-9 8． 02 × 10 E-9 9． 11 × 10 E-9

适合长基线的LC_AUTCLN 解( 基于码求解宽巷模糊度的LC 解)和LC_HELP 解( 基于电离层约束求解宽巷模糊度 的LC 解) ; 适合短基线的L1_ONLY 解( 双频观测中的L1 解) 、L2_ONLY 解( 双频观测中的L2 解) 、L1L2_INDENT 解( 双频观测中独立的L1、L2 组合解) 、L1＆L2 解( 双频观测中电离层约束的L1、L2组合解) 、L1 _ RECEIVER 解( 单频观测的L1解) 等。

对于L1L2 _INDENT和L1 _ONLY 和LC_AUTCLN三种解，由于短基线两端观测环境相关性极强，短基线宜采用差分模式消除观测误差的L1L2 _INDENT和L1 _ONLY 解，而不宜采用模型改正的LC_AUTCLN解。

22.基线长度对基准点的选择 基线处理时应注意如下3点:

1)对于短基线(小于10 ｋｍ), 其基准站坐标精度应优于

1ｍ, 只需将基准站与国家Ａ、Ｂ级网点联测0.5ｈ即可得到。 2)对于中长基线(小于100 ｋｍ), 其基准站坐标精度应优于 0.1 ｍ, 应取国家Ａ级网点作为基准点, 并与其同步联测24 ｈ 以上。

3)对于长基线(大于100 ｋｍ), 其基准站坐标精度应优于0.05 ｍ, 应取全球站作为基准站, 并与其同步联测48 ｈ以上。

22.周跳修复

23. gapr_to_l

sh_gamit calls gapr_to_l

? Creates an L-file from (time-dependent) coordinates and velocities from globk apr file

sh_gamit merges entries from existing lfile for sites not in apriori file

? If script still can’t find coordinates: --lfile坐标来源其它方法 – Uses pseudorange data in rinex file to generate point position or differential range position relative

to site in sites.defaults –rinex文件之中的伪距值产生，或者相对于sites.defaults中站的差分距离位置 – Use Rinex header coordinates (can block this by

setting use_rxc N in process.defaults) – Rinex头文件产生 – During the sh_gamit run, the coordinates are

updated if they are in error by > 30 cm ---若误差超过30cm，坐标将更新

各测站的近似坐标。GAMIT要求是球面坐标，可以先在测站的O文件里提 取近似坐标XYZ，生成lfile.apr文件，然后在/tables目录下输入命令： gapr_to_l lfile.apr lfile. “” yr doy

后续. 疑问： One.

guess_rcvant.dat

Two.If systematic patterns appear in sky map residual

plots, antenna type may be wrong

本文来源：https://www.bwwdw.com/article/i9j5.html