gamit说明文档-个人总结
gamit基线解算
任务名称：Gamit基线解算简介Gamit（Globk+Gamit）是一款用于全球GPS观测数据处理的软件，用来进行GPS的基线解算。
基线解算是指通过测量不同GPS接收站点之间的距离、方向角和倾角，确定不同站点之间的相对位置和变形情况。
Gamit基线解算是一种常用的地球物理测量技术，用于地壳变形、地震活动和大气延迟等领域的研究。
基线解算原理基线解算是通过GPS观测数据中的载波相位和伪距进行计算，并使用全球参考框架来确定地球表面不同站点之间的相对位置变化。
基线是指两个或多个GPS接收站点之间的距离和方向。
基线解算的原理是将GPS观测数据转换为坐标，然后通过数学模型和计算方法计算出站点之间的基线向量。
评估基线向量的精度和可靠性，可以帮助测量者判断地壳变形、地震活动以及大气延迟等现象。
基线解算的关键步骤包括：1.数据预处理：对原始GPS观测数据进行编辑、过滤和修正，消除掉仪器误差、信号传播误差和大气延迟等因素。
2.伪距差分：通过差分GPS观测数据，计算出不同站点之间的伪距差分，以获得更准确的GPS观测数据。
3.载波相位差分：通过差分GPS观测数据的载波相位，计算出不同站点之间的载波相位差分，以获得更高精度的GPS观测数据。
4.基线解算：根据差分后的GPS观测数据，使用数学模型和计算方法计算出不同站点之间的基线向量。
5.网运动学分析：通过基线解算结果，评估基线的精度和可靠性，判断地壳变形、地震活动和大气延迟等现象。
6.结果分析：对基线解算结果进行分析和解读，为相关领域的研究和应用提供依据。
Gamit基线解算流程Gamit基线解算流程包括以下主要步骤：1.数据准备：收集并准备全球GPS观测数据，在线或离线获取必要的测站坐标和地球物理模型数据。
2.数据预处理：对原始GPS观测数据进行编辑、过滤和修正，消除掉系统误差和环境影响因素。
3.伪距差分和载波相位差分：通过差分GPS观测数据的伪距和载波相位，计算出不同站点之间的伪距差分和载波相位差分。
gamit软件安装使用
三章 GAMIT在GPS反演大气水汽中的应用3.1 GPS数据处理软件GPS定位技术的普遍使用和深入研究使GPS已广泛应用到各个领域，其软件科学的发展和定位导航的需要促进了GPS 定位软件研发，同时GPS 精度的要求使精密解算GPS 观测数据的软件不断改进创新。
现在国际上普遍使用且精度较高的大地测量数据处理软件有如下软件：德国GFZ 地学中心的EPOS.P.V3 软件，瑞士BERNE 大学研制的Bernese 软件，美国麻省理工学院和加州大学圣地亚哥分校Scripps 海洋研究所联合研制开发的GAMIT/GLOBK 软件和美国宇航局喷气推进实验室研制的GIPSY/OASIS 软件。
针对工程应用而研制的商业软件主要包括TGO、TBC、Pinnacle、LGO、SKI-Pro 、Kiss以及国内各GPS 厂商自带的随机基线解算软件[22]。
本文研究采用高精度数据解算软件，下面对几种GPS高精度软件做简单的介绍。
（1）Bernese 软件Bernese 软件的功能非常强大，除了能定轨、定位、估计地球自转参数之外，还大量吸收融合各种有效改善定轨、定位精度的方法。
Bernese 软件能处理GPS 的两种数据，即非差和双差，用非差方法可以进行严密单点定位，也可用双差方法进行基线解算和整网平差。
此外，处理GPS 数据的同时还能处理GLONASS 及SLR 数据，重要的是它还能对GPS 数据和GLONASS 数据同时处置。
相比于其他高精度软件，Bernese的最大特点就是它能利用GPS 数据估计接收机天线的相位中心偏差及变化，并能够处理SLR 和GLONASS 卫星观测数据。
该软件由约一千个数据处理程序和百来个菜单程序通过文件有机地结合而组成。
经实验证明BPE 具有自动批处理功能且能够满足高精度定位的作用。
此软件运算速度快并且解算精度高，在大批观测量的数据解算中更能体现出一定的速度优势。
该软件的主体源程序由FORTRAN、Perl 语言写成并可在Windows、Unix 等多种计算平台上使用。
GAMIT安装步骤_修订版
GAMIT安装步骤_修订版1.系统要求在安装GAMIT之前，确保你的计算机符合以下最低配置要求：- 操作系统：支持Windows、Linux或Mac OS X。
-处理器：推荐使用64位处理器。
-内存：推荐至少8GB的内存。
-存储空间：至少需要10GB的可用磁盘空间。
- 软件依赖：安装了Perl、Python、Fortran和Gfortran等软件依赖。
3.安装依赖软件- Linux：使用包管理器安装Perl、Python和Gfortran等软件依赖。
例如，在Ubuntu系统上，可以使用以下命令进行安装：```sudo apt-get updatesudo apt-get install perl python gfortran```- Mac OS X：使用Homebrew或MacPorts等包管理器安装Perl、Python和Gfortran等软件依赖。
4.解压GAMIT软件包5.配置环境变量为了方便在任意位置运行GAMIT命令，需要将GAMIT安装目录添加到系统的环境变量中。
根据你的操作系统进行如下配置：```export PATH=$PATH:/path/to/gamit```其中，/path/to/gamit是GAMIT安装目录的路径。
6.配置GAMIT参数文件-GAMIT_BIN：GAMIT可执行文件所在的目录。
-GAMIT_GEN_DIR：GAMIT输出文件所在的目录。
-GAMIT_EXE_DIR：各个GAMIT模块的可执行文件所在的目录。
- GAMIT_SHELL：使用的命令行解释器，一般为bash。
-GAMIT_ARCHIVE：GAMIT源代码库所在的目录。
7.添加观测数据8.运行GAMIT完成上述步骤后，可以使用终端进入GAMIT安装目录，并执行以下命令来运行GAMIT：```./rungamit```总结：。
GAMIT添加自定义接收机天线
GAMIT添加自定义接收机天线第一步：在rcvant.dat 添加接收机和天线信息1.1 添加接收机信息特别说明，格式必须严格执行！！！在文件中RECEIVERS 部分添加接收机的简码(RECCOD)、IGS代码(IGS 20-char code)、描述(Description)下面为英文格式说明：# The (new) flags in column 37 for receivers indicate how differential code biases are to # be treated:# 'P' : receiver is cross-correlating and requires correction of P2' and C1# Rogue SNR, Trimble 4000, etc.# 'C' : receiver is non-cross-correlating but reports C1 instead of P1# Trimble 4700, 5700, Leica RS500, CRS1000, SR9600, etc. unless AS is off# 'N' : receiver is non-cross-correlating and reports true P1, P2# Blank indicates that we have not determined the correct flag and GAMIT will stop until # this file is modified to reflect the correct description of the receiver.##RECCOD IGS 20-char code | Description (comment field)#format(1x,a6,8x,a20) |#示例：添加天宝NETR9（假设不存在，其实存在的）TRNTR9 TRIMBLE NETR9 C | L1/L2+L2C/L5 GLONASS L1/L2 with 2 Maxwell-6 ASIC, eth + SBAS, 440 channel1.2 添加天线信息特别说明，格式必须严格执行！！！在文件中ANTENNAS 部分添加天线的简码(ANTCOD)、天线IGS代码(IGS 20-char code)、描述(Description)下面为英文格式说明：# Note: The 20-character RINEX code for antennas includes as of July 4, 1999, a designation# of 15 characters for the antenna, the 16th charcater is blank, and characters 17-20 are # reserved for antenna cover (radome) codes. The codes are tabulated at the end of thisfile.#ANTCOD ____FULL_NAME_______ |#format(1x,a6,8x,a15,1x,A4) |#示例：添加天宝TRM57970.00（假设不存在，其实存在的）TRZG2R * TRM57970.00 | Zephyr GNSS II - RoHS compliant solder L1/L2/L5/G1/G2/G3/E1/E2/E5ab/E6/Compass以上步骤是在rcvant.dat中添加天线和接收机信息，主要是将接收机全名和天线全名分别对应一个6位的简码！！！，文件格式须严格一致。
GAMIT软件Track模块使用帮助
GAMIT软件Track模块使用帮助TrackRT : Realtime GPS kinematic processing programTrackRTB : Realtime GPS kinematic processing program for use with Pre2.5 BNCTrackRTr : Rinex file emulator of trackRT for post-processing evaluationRunstring:% trackRT <options>% trackRTB <options>where options are-m <machine> Host name of system with real-time data port(Normally BNC would be running on this system)(Default local host)-p <port> Port number on host supplying the data stream(Mandatory input).-f <command file> trackRT command file (similar to track command file) -r <ref code> 4-char code for reference site.-d <list of 4-character codes> list of sites to be processed. If not specified all sites in data stream will be processed (there is a limit on maximum number of sites that can be processed depending on max_site in track. The referencecode should NOT be given here.)-n <root name> Specifies a root name to which output will be directed (.out is appended for stdout output). Use ? in the pos_root, sum_file and csv_file names the ? will be replaced by this string.Command file name must be given.% trackRTr <options> Same list of option although -m, -p have no effect -d <list of rinex files> list of rinex files (site names are generated fromthe first 4-characters of rinex file names)VERSIONS:---------1.13 -- Fixed time-equate iin SaveObsA so that a differnce of morethan 0.1 seconds in needed for a new epoch (trackRT only)Added reference relative humidity for GPT. (120103)1.12 -- Handle bad PRN from BNC where G23 == 203.1.11 -- Version to handle missing satellites in SP3 files (110505)1.10 -- Version to handler BNC2.5 and greater that is now outputing ascii rather than binary data (trackRTB handles BNC2.0) (110501) 1.00 -- Initial version 100301 (2010 March 1)INTRODUCTION:-------------TrackRT/TrackRTB and trackRTr have many commands that allow it be controlledboth when started and during runs. Below the detailed commands are given and here we give the minimum commands needed.SP3_DIR : (unless sp3 files are in the current directory)SITE_POS: Site positions must be given for every site to be processed. SITE_STATS: To set the apriori sigma and process noise for the sites. NOTE: The reference site position should be set to zerofor the apriori sigma and process noise. Typical site_stats would be where cit1 here is the reference site.site_statsall 0.1 0.1 0.1 0.025 0.025 0.025cit1 0.0 0.0 0.0 0.0 0.0 0.0ATM_STATS: By default atmospheric delays are not estimated, and normally these should be. Normally the reference site is set to zero when site separations are less the 500-1000km. Typicalatm_statscommand would beatm_statsall 0.20 0.00010 0.000 ! Unit m/sqrt(sec) -> 0.0001 = 0.03m/sqrt(day)cit1 0.00 0.00000 0.000If mixed antenna and receiver types are usedANTE_OFF: To specify the antenna heights and types and receiver types ANTMOD_FILE: Must be given to get the antenna phase center modelsRCV_TYPE: Needed for mixed receiver types (entry can be specified in ANTE_OFF also and this command would not be needed).DCB_FILE: Up to date, data code bias (DCB) file (part of GAMIT ftp area).IncludingUPDATE_FILE is useful so that trackRT can be controlled on the fly.TrackRTB should be used for older versions of BNC. If no files are created by trackRT or trackRTB (list of files should be listed shortly after theprogram is started) then try using the other version. The program are identicalexcept trackRTB is expecting a binary data stream.COMMANDS--------Input/Output commands---------------------@ SP3_DIR <Directory> <center><Directory> -- Directory where sp3 files are stored<center> -- Center for orbits (default igs, igr and igu also tested)@ POS_ROOT <root> <duration>Set the root part of the name for the output files.<root> Root part of name. When ? included in the <root>, the ? is replaced with the -n string.<duration> Duration of data in each file. The designations of d, h, or mmay be used to specify the units of days, hours, minutes. Default is d.The position file names take the form: (see commands below)<root>.<outtype>.<site>.<start time>.<datatype>The resolution of the <start time> which is modulo the output interval depends on the output duration. For output durations greater than or equal to 1-day, the time is given as YYYYMMDD. For intervals shorter than 1-day, it is YYYYMMDD:HHMN. Minium output interval is 1 minute.@ SUM_FILE <root>Sets the root part of the summary file name. Using ? in the namewith be replaced by the -n string. (Default if command is not givenis trackRT or the -n string when -n used). File names are time tagged according to the pos_root output interval.@ CSV_ROOT <root>Set the root part of the name for comma separated values (CSV) output file. These files are used for AmCharts web plots.@ DCB_FILE <file name>Set the name of the data-code-bias (DCB) file. This file is part of the GAMIT tables directory and should be updated regularly. It is used to remove biases in the Melbourne-Wubbena widelanes. The receiver type can be specified with the RCV_TYPE or ANTE_OFF commands.@ ANTMOD_FILE <file name>Sets the name of a standard IGS antex file with phase center models for the GPS ground antennas. Antenna types at specific sites are given with the ANTE_OFF command. This command my be used multiple times for site specific model with new models replacing previously read ones.@ UPDATE_FILE <file name>Allows new trackRT commands to be issued during a run. Once the file is read it needs to be deleted before the trackRT will re-read it. File is only read if it exists. NOTE: File should be removed before trackRT is run or else it will be read when the command file is read (ie., it will overwrite the commands in the command file.@ OUT_TYPE <NEU+GEOD+XYZ+DHU>Specifies types of output coordinates. All types can specified in a single string with no spaces. The types areNEU -- North, East, Up differences from the reference site or from the coordinates given in the REF_NEU command.GEOD -- Geodetic latitude, longitude and height (in the GEOD format, the total atmosheric delay is given, while in the other formats the adjustment to the apriori delay is given).XYZ -- Cartersian XYZ coordinatesDHU -- Delta horizontal and Up coordinates from the apriori coordinates of each site (default output type)@ OUT_SIG_LIMIT <sigma (m)>Sets the maximum sigma of a position estimate for it to be output.If pseudorange data types are used, the default value of 1 m needsto increased to 10-100 meters.Analysis commands-----------------@ DATA_NOISE <L1 (m)> <L2 (m)> <P1 (m)> <P2 (m)> <Elev Weight> [PRN] Allows specificiation of the noise in the L1 phase, L2 phase,P1 range and P2 range, and the weight given to elevation angle depepence (at ver 1.20); variance is scale by (1+(W/sin(el))^2)where W is the <Elev Weight>. These values affect the sigmas printed for the position determinationsOptional: PRN may be added and noise assigned to that PRN (ifnon-PRN form is used, this will replace all PRN specific valuesso use the non-PRN first followed by specific PRN values@ DATA_TYPE <type>Specifies type of data to process. Types supported are L1 L2 LC P1 P1 PCwhich can be combined in a contiguous string. Files names include this string at the end.Vers 1.0: Only LC, LCPC and PC have been tested. Specifically L1 only datahave not been fully implemented yet.@ USE_GPTGMF <Relative humidity (0-1)>Set the used of the GPT temperature pressure model and the GPT dry and wetmapping functions. Optional relative humity added vers 1.27; default 0.00.Default is the older MTT atmospheric model.Site and Antenna informationSite information is entered by first giving the command and then the information by site on the subsequent lines@ SITE_POS@ Site <X (m)> <Y (m)> <Z (m)> <Vx (m)> <Vy (m)> <Vz (m)> <Epoch (yrs)>Site is the four character name of the site (more characters can be includedbut only the first 4 are checked). Site names that do not appeat in the listof sites to be processed are ignored). The remainder of the line contains positon and velocity and the epoch in deciminal years to which the position refers.@ ANTE_OFF@ Site <ARP dN (m)> <ARP dE (m)> <ARP dU (m)> <Antenna Name> <Receiver Code>Specifies the type of antenna and its position of antena reference point (ARP)at each site. The antenna name including radome should be specified with theofficial IGS name for a standard ANTEX file or with a unique name that appearsin the ANTEX file for site specific calibrations. (Note: There is oneadditional character in the long antenna and this extra characterbefore the radome name must be removed. The antenna name and radome can becopied directly from the rinex file if present. One more ANTEX files must bespecified with the ANTMOD_FILE command for the antenna names to be useful. The receiver type DCB code can be optionally specified here as well (See RCV_TYPEcommand).@ RCV_TYPE@ Site <Receiver code N/P/C>Specifies the type of data-code-bias (DCB) correction needed for the receiver.Code specifies the type of L1 and L2 ranges being measures. The choices areP -- Pcode, C -- C/A and N C/A with cross corelation for L2 range. The codescan be found in gamit/tables/rcvant.dat. These codes can also be given in theante_off command. An up-to-date DCB_FILE command must be used to specify theDCB biases. The files are available from the MIT ftp site and update once permonth.@ SITE_STATS@ Site <Apriori Sigmas in XYZ (m)> <RW noises in XYZ (m/sqrt(sec))> Gives statistics to assign to the kinematic station positions. The <Apriori Sigma in XYZ> are the three sigmas in XYZ for the initial postion and<RW noise in XYZ> are the three sigmas for the random walk in position. Unitsof the ransom walk are m/ssqrt(s).ALL can be used for the station name and the same statistics will be applied to all kinematic sites (NOTE: the fixed site do not change position).e.g.site_statsall 1.0 1.0 1.0 0.010 0.010 0.010base 0.0 0.0 0.0 0.000 0.000 0.000(1 meters apriori sigmas and changes 10 mm/sqrt(s)=600 mm/sqrt(hr)).NOTE: One station should always set to zero sigmas and random walk otherwise allstation positions will be estimated and noramlly drift by large amounts. The fixed station does NOT need to be the reference site.@ ATM_STATS@ Site <Apriori Zenith delay sigma> <RW noise in Zenith delay> <RW dH/dt noise>Gives the statistics for the atmospheric delays by site. A random walk process'smoise increases as sqrt(time) where time is time in seconds. The values are the initial sigma in meters, RW changes in meters per sqrt(second) and (addedverion 1.2) a dH/dt variance term so that during rapid height changes more process noise can be added to zenith delay estimate. The process noise variance is(<RW dH/dt noise)*abs(dh/dt)>^2 per sqrt(sec) where dh/dt is m/s. The default valeis 0.00023 m^2/sqrt(s). (Term only needed for aircraft processing). e.g.,atm_statsall 0.10 0.0001 0.00000t39a 0.10 0.0001 0.00023base 0.00 0.0000 0.00000Sets the apriori sigma as 10cm and allows the delay to change 6 mm/sqrt(hr) and2.3 mm/sqrt(s) when height is changing at 10 m/s (fast ascent or desent) for thesite t39a (aircraft). In the case, here "base" is a reference station at whichthe atmospheric delay is not estimated (other stations absorbe the atmosphericdelays at the base station). With kinematic positioning over few hundred km siteseparations, the atmospheric delay at one site is normally fixed.Ambiguity resolution/data control commands.-------------------------------------------trackRT uses a combination of the Melbourne-Wubbena widelane (MW-WL), the extrawidelane (EX-WL) and the floating point estimates of the ionospheric free ambiguity (LC)to resolve integer ambiguities. If we denote the number of integer cycle ambiguitiesat L1 and L2 by N1 and N2, The MW-WL is an estimate of N1-N2 based on phase and rangedata; EX-WL = N1 - f1/f2 N2 and is an integer for L1 cycles, but 1.283 N2 for L2 cycles;LC = 2.546 N1 - 1.984 N2. The EX-WL is unaffected by geometric ranges changes, but doesdepend on the ionospheric delay. For short baselines, the EX-WL should be near zero forcorrect choices of N1 and N2. The LC residual should also be near zero when N1 and N2are correct. The problem in ambiguity resolution is that difference choices of N1 andN2 can make different linear combinations small. For example, errors in N1 and N2 of3 and4 cycles will change LC by 0.298 cycles (56.6 mm), MW-WL by 1 cycle and EX-WL by2.132 cycles (405 mm). On long baselines, at low elevation angles, ionospheric delaysof 400 mm are common. The most common error is a N1=N2=1 cycle error. For thiscombination, the MW-WL is unaffected and LC changes by 0.562 cyc (107 mm) and EX-WL by0.283 cyc (54 mm). Even on relative short baselines, 54 mm ionospheric delays are commonThe follow commands are use to control the ambiguity resolution.@ AMB_SET <RelRank> <FloatSigma (2)> <MWWL Fact> <EXWL Fact> <Min AmbSig> <MaxChi>Sets parameters for ambiguity resolution. The input parameters are: <RelRank><WL min> Minumum of values need to allow bias fixing<WL avN> Maximum number to be used in computing sigma of mean MW-WL<FloatSigma (2)> Minimum sigma for LC and MW-WL for ambiguity fixinf <MWWL Fact> Weighting factor for MW-WL in chi**2<EXWL Fact> Weighting factor for extra-wide lange<Min AmbSig> Minimum sigma to assigned to float estimates of ambiguities.<MaxChi> Max chi**2 value allowed for ambiguity to be resolved. @ EXWL_SET <Jump> <Min Sigma> <Scale> <Elev Fact><Jump> -- magnitude of jump in EX-WL to have cycle slip added (default 0.10)<Min Sigma> -- minimum sigma for mean ex-wl (cycles, default 0.02 cycles) <Scale> -- Scaling factor for length. Scale 0.1 results in 0.1 cycles over 100 km (default)<Elev Factor> -- Elevation angle factor that increases sigma as (1 + factor/sin(elev))@ MWWL_SET <Jump> <Min Sigma> <Max Averaging number> <Min number> <Jump> -- magnitude of jump in MW-WL to have cycle slip added (default 5.0)<Min Sigma> -- minimum sigma for mean MW-WL (cycles, default 0.10 cycles) <Max Averaging number> -- Maximum number of values to use to compute mean sigma<Min number> -- Minimum number needed to resolve ambiguity@ DD_SET <Jump (cycle)> <Min Number>Sets parameters for double difference processing.<Jump (cycle)> -- Magnitude of jump in double differences on bias fixed data thatwill introduce a cycle slip<Min Number> -- Minumum number of double differences for epoch to be processed.(If too few than errors in the data can be detected than this cancause large position errors, default of 4 double differences allowsredunancy).@ RMS_EDIT_TOL <n-sigma Tolerance> <min sigma> <Reset number><n-sigma tolerance> is an n-sigma condition where sigma is based on data noise model.<min sigma> Minimum phase sigma to use to that no phase residual less than <<min sigma>*<n-sigma tolerance> are deleted<Reset number> number of sequential delete data, before ambigity and cycle slipsare reset (assumed missed cycle slip).Testing and evaluting commands@ STATUS <type> <# epoch>Writes status information to the current summary file at <# epoch> intervalsThe types of reports are given by "type"P -- Parameter estimatesA -- Ambiguity resolution report (shows resolved and unresolved)W -- Widelanes (Melbourne-Wubbena and Extra-widelanes)R -- Postfit residuals are current epochC -- Report current A and W entries onlytype PAWR will output all reports.@ NUM_EPOCHS <number><number> is the number of epochs of data to be processed before stopping. When shortduraton files are output (POS_ROOT command) the number of epochs may not be reachedbecuase the epoch counter is reset with each new file.@ START_TIME <yy mm dd hh min sec.>Used to set the start time. Useful with trackRTr to have the processing start at thesame time as the realtime stream.@ DEBUG <Up to 10 epoch numbers in pairs>Sets which epochs will report detailed debug and status information (see STATUS commandas well). The pairs are used as:1,2 -- Mostly model information. Useful if data are generating large residuals3,4 -- Parameter estimates, widelanes, residuals5,6 -- One-way OMC and single differences7,8 -- Not Used9,10 -- When 9 is non-zero, antenna model and SP3 information.@ RESET <ALL/list of sites>Reset command. Resets the filter state vector and resolved ambiquities for alist of sites or ALL sites (generally used in the update_file to fix problems)@ EXCLUDE_SVS <list of PRN numbers to be excluded>Excludes satellites from being processed. Useful when a satellite is not in an SP3 file or not available during the times there are data. Example: exclude_svs 26 22 13Output types------------The output position files have two header lines that explain the columnsin thefile. The output in the summary and output files contain the followingbasicblocks are written to the summary file. The time in these files is givenbyan epoch number which is a count of the number of epochs of data (at the referencesite) since the start of the current data files. The output positionfiles givethis counter and the corresponding GPS times and deciminal days.CSLIP line are cycle slip detectionsCSLIP Ep 540 Site THMG G 29 DMW-WL/Tol 4.31 4.00 DEX-WL/Tol 1.78 0.20 cyc, Elev 10.00 degIndicates the jump in the MW-WL and EW-WL (and the tolerance for the detection)that caused the slip to be detected.CSLIP Ep 961 Site THMG G 12 DD RESID 2.18 Tol 0.50 cyc, Elev63.85 degIndicates a slip added due to repeated double difference residuals.When ambiguities are resolved, three lines are output giving thestatisticsAMBFIX WMAP PRN29 EP 599 RelRank 101.76 FC ------ dL12 4 3Dchi 6.54 691.02 AMB 36AMBWLS WMAP PRN29 EP 599 RG 133 599 FX 3 iL12 0 0Means -0.16 -0.05 RMS 0.74 0.06 # 467 eN12 -0.54 -0.38AzEl 274.32 12.35AMBCON WMAP PRN29 EP 599 NCont 3 MW Res -0.159 0.185 Chi2 0.74 EXRes -0.051 0.068 Chi2 0.55 LC Res -0.046 0.020 Chi2 5.25AMBFIX gives the epoch number, relative rank, the Fix Code (always ------when fixed), thechanges to the number of L1 and L2 cycles from the intially asigned values,thechanges in chi**2 for the best and second best choices of ambiguities andthenumber of the ambiguity resolved. When the status feature is used, ambiguitiesnot resolved yet are listed as AMBFREE entries.AMBWLS gives the mean values and RMS scatters of the MW-WL and EX-WL (cycles).The iL12 estimate here is the changes in L1 and L2 cycles implied by justthe widelanes with no contribution from the LC estimates. On short baselines,these values are expected to be zero.AMBCON gives the contributiond to chi**2 from each of the input type. Ncontis the number of contributions (3 for LC data type). Values given are the residal (Res) with its sigmas and Chi**2 (Chi2) contribution,When STATUS output is requested the following output types are possible. STATUS REPORT Epoch 1800 Date 2010 05 06 08 36 47.000 Type PAWR PARAMETER ESTIMATES Epoch 1800 for NP 16 Type P STATUS REPORT...PARAMETER Ep 1800 NP 1 MSCG dATM -0.0466 +- 0.0001 m Var 0.681E-08PARAMETER Ep 1800 NP 2 PSDM dATM -0.0001 +- 0.0001 m Var 0.436E-08...PARAMETER Ep 1800 NP 7 LJRN G 27 1 0.1281 +- 0.0253 cyc Var 0.639E-03PARAMETER Ep 1800 NP 8 LJRN G 05 1 6.1172 +- 0.0011 cyc Var 0.112E-05Parameter values for position adjustments (XYZ), atmospheric delay adjustementsand float ambiguities are given.Ambiguity report for both both fixed and free ambiguities. ResEpoch isthe epoch whenthe ambiquity was fixed (if fixed in previous data file block, may be greater thancurrent count).AMBIGUITY Report Ep 1800 Number of ambiquities 59AMBFIXD # 11 MSCG PRN10 RelRank 616.27 FCode ------ dL12 -1-1 Dchi 1.01 776.94 Elev 72.02 deg; ResEpoch 156AMBFIXD # 12 MSCG PRN27 RelRank 135.91 FCode ------ dL12 -1-1 Dchi 5.39 766.54 Elev 29.86 deg; ResEpoch 156..AMBFREE # 41 LJRN PRN27 RelRank 11.37 FCode RW---- dL12 00 Dchi 25.69 294.89 Elev 30.60 deg; ResEpoch 0AMBFREE # 42 LJRN PRN05 RelRank 11.76 FCode R----C dL12 76 Dchi 106.15 1251.17 Elev 36.16 deg; ResEpoch 0The FCode gives the reason for not being fixed the entries are: RWSSSCR - Relative rankW - Widelane sigma too large or not enough data yet (see MWWWL_set command)S - LC sigma too largeS - L1 sigma too large when L1+L2 data type used (not tested)S - L2 sigma too large when L1+L2 data type used (not tested)C - Chi**2 increment too large.Widelane report shows status of average values of widelanes and current LC estimate. X means fixed, R ambiguity is still free.WIDELANE Report Ep 1800 Number of ambiquities 59WIDELANE 11 MSCG PRN10 EP 1800 Range 2 1800 # 1799 NC 3 X MW Res 0.06 0.10 Chi2 0.30 EX Res -0.18 0.19 Chi2 0.82 LC Res -0.01 0.02 Chi2 0.10...WIDELANE 19 MSCG PRN29 EP 1800 Range 421 1800 # 1380 NC 3 R MW Res -0.09 0.15 Chi2 0.35 EX Res 0.07 0.26 Chi2 0.07 LC Res 0.17 0.02 Chi2 74.30Postfit residuals for data types being used.POSTFIT RESIDUAL Report Ep 1800 Number of DD 76POST 1 Ep 1800 MSCG PRN10 - CIT1 PRN02 Res -0.0945 +- 0.0381 LC AMB 11 Elev 37.31 degPOST 2 Ep 1800 MSCG PRN27 - CIT1 PRN02 Res -0.1040 +- 0.0441 LC AMB 12 Elev 22.17 deg...POST 41 Ep 1800 MSCG PRN04 - CIT1 PRN02 Res 2.9712 +- 1.8988 PC AMB 13 Elev 30.60 degPOST 42 Ep 1800 MSCG PRN30 - CIT1 PRN02 Res -0.2485 +- 1.8116 PC AMB 14 Elev 36.16 deg..++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++EXAMPLE COMMAND FILE:--------------------* FILE trackRT_pbo.cmd* Example TrackRT command file to process some PBO sites.* Run BNC with as the TCP/IP port (login information * can be obtained from rtgps@.) and select sites P496, P497, * P498, P505 to be cast on specific port (port 3765 for example).* Run trackRT as* trackRT -p 3765 -r p497 -d p496 p498 p505 -f trackRT_pbo.cmd -n P497 &* (if BNC is run on another machine used -m <machine name> as well.* A relative directory ../sp3_files is assumed to contain up to date* sp3 file.sp3_dir ../sp3_files* These positions are post April 4, Baja earthquake (extracted from * globk output with sh_exglk -f <globk org file> -apr <apr file>* entries from apr_file* Nore: Additional sites can be given heresite_posP498_GGU -2313650.15153 -4835550.43945 3444474.55150 -0.00658 0.02008 0.01478 2010.259 0.0033 0.0059 0.0042 -1.0000 -1.0000 -1.0000P497_GGU -2315937.74208 -4838712.09671 3438545.10828 -0.00945 0.01927 0.01884 2010.259 0.0035 0.0063 0.0044 -1.0000 -1.0000 -1.0000P496_GGU -2319722.10265 -4842511.39340 3430709.79068 -0.01697 0.02370 0.02587 2010.259 0.0050 0.0089 0.0062 -1.0000 -1.0000 -1.0000P503_GGU -2325056.31018 -4826758.67413 3449209.09091 -0.01167 0.01951 0.01660 2010.259 0.0033 0.0059 0.0042 -1.0000 -1.0000 -1.0000P505_AGU -2309739.55681 -4802072.50955 3493258.58085 -0.00377 0.00608 0.00160 2010.404 0.0033 0.0057 0.0040* PBO sites antenna and receiver information (all are NetRS and thus C type)* Again more sites can given). (Note: If antenna information extracted from* gamit , one space must be removed before radome string). * Antenna/Radom combination must be in antmod_file for model to be used). ante_offp475 0.00 0.00 0.0083 TRM29659.00 SCIT Cp066 0.00 0.00 0.0083 TRM29659.00 SCIT Cp472 0.00 0.00 0.0083 TRM29659.00 SCIT Cp478 0.00 0.00 0.0083 TRM29659.00 SCIT Cp494 0.00 0.00 0.0083 TRM41249.00 SCIT Cp496 0.00 0.00 0.0083 TRM29659.00 SCIT Cp497 0.00 0.00 0.0083 TRM29659.00 SCIT Cp498 0.00 0.00 0.0083 TRM41249.00 SCIT Cp500 0.00 0.00 0.0083 TRM29659.00 SCIT Cp503 0.00 0.00 0.0083 TRM29659.00 SCIT Cp505 0.00 0.00 0.0083 TRM29659.00 SCIT Cp510 0.00 0.00 0.0083 TRM29659.00 SCIT C* Change path as needed.antmod_file /home/tah/gg/tables/antmod.datdcb_file /home/tah/gg/tables/dcb.datdata_type LCPCsite_statsall 0.1 0.1 0.1 0.025 0.025 0.025p497 0.0 0.0 0.0 0.0 0.0 0.0atm_statsall 0.20 0.00010 0.000 ! Unit m/sqrt(sec) -> 0.0001 = 0.03m/sqrt(day)p497 0.00 0.00000 0.000pos_root ? 1h# Used at MIT for web output. Results can be viewed at# /kmeduna/# csv_root /net/chandler/var/www/kmeduna/trackrt_view/data/P497# Useful to have this option. Note file should be removed before trackRT # is run and created when needed.update_file upd_app.cmd+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ +++++++++INSTALLATION------------trackRT requires the BKG NTRIP Client (BNC) and QT libraries and include files.Ntrip is the Networked Transport of RTCM via Internet Protocol and the clientprogram BNC provides access to the realtime data streams and casts the datevia an internet port GNSS raw data in a structure that is decoded by trackRT.Currently trackRT only uses the GPS data streams. TheNTRIP system is discussedat http://igs.bkg.bund.de/ntrip/. The BNC program can be downloaded fromhttp://igs.bkg.bund.de/ntrip/download. The BNC documentation constains theinstructions for obtaining the QT libraries and includes needed for its installation from source files. trackRT only needs the executable version ofBNC. The QT software can be obtained directly from Noia at the site /downloadsFor trackRT, the libraries and include files are needed and these are containedin Frameworks libraries 4.6.2 for your platform (linux or Mac). Earlier versionshould also work. For Linux and Max OSX, QT can be downloaded directly from./downloads/linux-x11-cpp/downloads/qt-for-open-source-cpp-development-on-m ac-os-xQT takes several hours to install.On linux systems, other tools that can be needed are:sudo app-get install gfortransudo app-get install g++sudo app-get build-dep qt4-qmake(The last app-get may be need if the qt install can not find X11 libraries andincludes.)For MacOSXThe basic qt install steps aresudo port install qt4-macIf you don't have MacPorts, it's a very easy way to get lots of *nix tools for the Mac. You can download MacPorts from .Once QT is installed, the makefile in the trackRT directory should be modifiedto include the directories for the libraries and includes. Example entries aregiven in the makefile. The actual entries are for the MIT QT installation (whichis not quite standard).NOTE: There is a file makefile (not Makefile) for making trackRT series of programs.To run trackRT, bnc is run first and the desired stations added the mountpoints.。
GAMIT使用
GAMIT使用
GAMIT的使用需要依赖于多个输入数据源。
首先，需要GPS观测数据
来确定地球表面的位置。
这些观测数据通常通过全球分布的GPS测量站获得。
其次，需要卫星轨道测量数据来确定卫星位置和运动的精度。
最后，
还需要一个大地坐标系以及与卫星测量数据相关的各种参数。
这些数据最
终用于计算地壳形变的模型估计。
GAMIT的工作流程主要分为数据处理和模型估计两个部分。
在数据处
理方面，数据首先通过一个预处理阶段进行数据检查和修正，以去除可能
的误差和异常。
接下来，数据经过一个精确的观测数据处理流程，包括时
钟和轨道参数的建模和估计。
最后，通过使用所有这些数据和参数，计算
出地壳的运动模型。
在模型估计方面，GAMIT使用单差估计方法，将所有的观测数据与参
考站数据进行比较。
通过不同站点之间的差异，可以进一步估计出地壳运
动的速率和方向。
此外，GAMIT还可以提供地壳变形模型的参数，如弹性
滞后和地震震级。
这些估计结果对于地球物理学家来说非常重要，可以帮
助他们更好地理解地壳运动的动力学过程。
GAMIT具有一些突出的特点和优势。
首先，它是一个开源的软件工具，可以免费使用。
其次，它具有较高的计算精度和稳定性，可以处理大规模
的观测数据。
此外，GAMIT还具有较好的可扩展性，可以与其他地球物理
研究工具进行集成和扩展。
最后，GAMIT提供了用户友好的图形界面和详
细的文档，以帮助用户更好地理解和使用该工具。
GAMIT软件Track模块使用帮助
信息系统安全与风险管理的重要性与策略概述：随着信息技术的快速发展，信息系统安全与风险管理变得愈发重要。
在数字化时代，大量的数据和信息储存在电子设备和网络中，这使得信息系统更容易受到攻击和威胁。
本文将探讨信息系统安全和风险管理的重要性以及一些策略，以帮助组织确保其信息系统的安全性和可靠性。
一、信息系统安全的重要性信息系统安全指的是保护系统中的信息和数据免受未经授权的访问、窃听、篡改或破坏。
信息系统安全至关重要，主要原因如下：1. 保护数据和隐私：信息系统中存储着大量的敏感数据和个人信息，包括财务数据、客户信息和商业秘密等。
如果这些数据泄露或被攻击，将对组织和个人造成严重损失，包括财务损失、声誉受损以及法律风险等。
2. 确保商业连续性：当信息系统受到攻击或遭受故障时，组织的业务运转可能会受到影响，导致生产中断、服务停止或流程受阻。
确保信息系统的安全性可以减少这些风险，保证组织的商业连续性。
3. 遵守法规和合规要求：根据国家法规和行业标准，组织有责任保护其信息系统的安全性。
未能遵守法规和合规要求将面临巨大的法律和经济风险。
4. 保护知识产权：信息系统安全对于保护组织的知识产权至关重要。
公司的重要机密、专利和商业计划等必须受到安全保护，以防止未经授权的披露或泄露。
二、信息系统风险管理的重要性信息系统风险管理旨在识别、评估和应对潜在的风险和威胁。
以下是信息系统风险管理的重要性：1. 提升整体安全性：通过风险管理流程，组织可以更好地理解其信息系统面临的威胁和脆弱性。
这将有助于采取适当的措施来保护信息系统，提高整体的安全性。
2. 减轻损失：风险管理活动有助于预见和减少潜在的风险和威胁。
通过及时应对和预防，组织可以减轻潜在损失的影响。
3. 合理分配资源：风险管理还可以帮助组织识别和优先处理最高风险的区域和问题。
这有助于合理分配有限的资源，确保最大的安全回报。
4. 遵守法规要求：风险管理是组织履行合规要求的一种方式。
GAMIT使用笔记
大气层研究和空间空间电离层研究使用到是GAMIT模块，精密定位还GAMIT、GLOBK两个模块都需要。
安装完成后的几个重要文件：gg/gamit（基线平差）和gg/kf（Kalman Filter）两个目录下到模块是用fortran编写的。
gg/com是cshell编写到脚本，重要用于gamit和kf目录下的模块的组织。
gg/tables是表文件。
年更新LUNTAB、SOLTAB 、NUTABL、LEAP.SEC周更新UT1、POLEsh_gamit批处理要求工程目录下至少有rinex brdc gfiles三个目录。
分别放O文件，N文件，卫星轨道文件g文件，这样做的目的是把文件分类，最后这些文件都会被link到单天的目录之下。
注意：需要将所有观测文件和表文件都link到单天目录下的，sh_gamit能自动完成link功能。
模型说明：1.otl 潮汐改正2.vmfl GMF 投影函数3.atml大气荷载模型，对高程影响较大，可消除周跳波动，可靠性需要进一步证实4.atl大气抄袭荷载模型和met气象模型星历文件：e/n, sp3, g,te/n为广播星历，主要用来你和卫星和接收机的种差g文件是根据sp3文件拟合的某天的圆形轨道参数t文件是根据观测文件和g文件求出的卫星位置，是gamit专用格式gcc编译器作用：将常见的编程语言转化为c语言。
安装gcc需要把原来到gcc覆盖。
在/usr目录下，具体怎么做，不是很清楚。
软件中的栅格文件：下载地址：ftp://1)海洋潮汐。
例如otl_FES2004.grid放在软件talbels目录下。
链接到otl.grid。
2)大气负荷。
例如atmldisp_cm.2006，每年更新一次。
连接到atmldisp_grid.20063)vmfl投影函数栅格，例如vmflgrd2006，连接到map.grid.2006。
每年更新一次。
以example为例作一个实例：1）在/media/Tool/TOOL/专业工具/GAMIT下新建文件夹10-05-18-EXAMPLE，在该目录下建立tables目录。
解算软件GAMIT处理流程文档
解算软件GAMIT处理流程文档.txt小时候觉得父亲不简单，后来觉得自己不简单，再后来觉得自己孩子不简单。
越是想知道自己是不是忘记的时候，反而记得越清楚。
gamit 处理流程流程概览：1、规划工程2、更新table文件夹3、创建工程目录，连接table4、编辑sites.defaults5、手动建立其余必须目录6、下载、准备数据7、制作准备文件8、sh_gamit批处理解算9、具体介绍：1、工程规划：处理中国大陆IGS站8个（wuhn、shao、chan、bjfs、kunm、guao、urum、lhaz、），台湾IGS站2个(tnml、twtf)数据跨度：2008年10月1日~2008年10月5日，2009年10月1日~2009年10月5日框架站（6个）：wuhn、chan、bjfs、kunm、lhaz、twtf2、更新~/gg/tables文件夹(需要一定权限)可到/processing/gamit/ 更新所需表文件，或者到ftp:///archive/garner/gamit/tables/ 下载（用户名：ananymous 密码：邮箱）通常来自IERS的地球自转参数表ut1. 、pole. 需要每日或周更新章动表nutabl. 太阳星历表soltab. 月亮星历表luntab. 跳秒表leap.sec 需要每年更新接收机/天线列表rcvant.dat 在有新仪器使用时需要更新卫星列表svnav.dat 在有新卫星使用时更新码偏文件(P1-C1、P1-P2)dcb.dat 需要月更新如果想方便的话，每次解算时全部更新3、建立工程目录，例如china2008,同时运行sh_setup 命令链接talbes文件夹可以查看ut1. 文件，是否在你处理时间内，处理的是2008年数据，显然在范围内。
4、编辑工程china2008文件夹下的talbes中的sites.defaults文件，把需要处理的测站加进去，并规划参考固定站(globk需要)。
GAMIT软件Track模块使用帮助
GAMIT软件Track模块使用帮助TrackRT : Realtime GPS kinematic processing programTrackRTB : Realtime GPS kinematic processing program for use with Pre2.5 BNC TrackRTr : Rinex file emulator of trackRT for post-processing evaluationRunstring:% trackRT <options>% trackRTB <options>where options are-m <machine> Host name of system with real-time data port(Normally BNC would be running on this system)(Default local host)-p <port> Port number on host supplying the data stream(Mandatory input).-f <command file> trackRT command file (similar to track command file) -r <ref code> 4-char code for reference site.-d <list of 4-character codes> list of sites to be processed. If notspecified all sites in data stream will be processed(there is a limit on maximum number of sites that can beprocessed depending on max_site in track. The referencecode should NOT be given here.)-n <root name> Specifies a root name to which output will be directed(.out is appended for stdout output). Use ? in thepos_root, sum_file and csv_file names the ? will bereplaced by this string.Command file name must be given.% trackRTr <options> Same list of option although -m, -p have no effect -d <list of rinex files> list of rinex files (site names are generated fromthe first 4-characters of rinex file names)VERSIONS:---------1.13 -- Fixed time-equate iin SaveObsA so that a differnce of morethan 0.1 seconds in needed for a new epoch (trackRT only)Added reference relative humidity for GPT. (120103)1.12 -- Handle bad PRN from BNC where G23 == 203.1.11 -- Version to handle missing satellites in SP3 files (110505)1.10 -- Version to handler BNC2.5 and greater that is now outputingascii rather than binary data (trackRTB handles BNC2.0) (110501) 1.00 -- Initial version 100301 (2010 March 1)INTRODUCTION:-------------TrackRT/TrackRTB and trackRTr have many commands that allow it be controlled both when started and during runs. Below the detailed commands are given and here we give the minimum commands needed.SP3_DIR : (unless sp3 files are in the current directory)SITE_POS: Site positions must be given for every site to be processed.SITE_STATS: To set the apriori sigma and process noise for the sites.NOTE: The reference site position should be set to zerofor the apriori sigma and process noise. Typical site_statswould be where cit1 here is the reference site.site_statsall 0.1 0.1 0.1 0.025 0.025 0.025cit1 0.0 0.0 0.0 0.0 0.0 0.0ATM_STATS: By default atmospheric delays are not estimated, and normally these should be. Normally the reference site is set to zerowhen site separations are less the 500-1000km. Typical atm_statscommand would beatm_statsall 0.20 0.00010 0.000 ! Unit m/sqrt(sec) -> 0.0001 = 0.03 m/sqrt(day) cit1 0.00 0.00000 0.000If mixed antenna and receiver types are usedANTE_OFF: To specify the antenna heights and types and receiver types ANTMOD_FILE: Must be given to get the antenna phase center modelsRCV_TYPE: Needed for mixed receiver types (entry can be specified inANTE_OFF also and this command would not be needed).DCB_FILE: Up to date, data code bias (DCB) file (part of GAMIT ftp area).IncludingUPDATE_FILE is useful so that trackRT can be controlled on the fly.TrackRTB should be used for older versions of BNC. If no files are createdby trackRT or trackRTB (list of files should be listed shortly after theprogram is started) then try using the other version. The program are identical except trackRTB is expecting a binary data stream.COMMANDS--------Input/Output commands---------------------SP3_DIR <Directory> <center><Directory> -- Directory where sp3 files are stored<center> -- Center for orbits (default igs, igr and igu also tested)POS_ROOT <root> <duration>Set the root part of the name for the output files.<root> Root part of name. When ? included in the <root>, the ? is replaced with the -n string.<duration> Duration of data in each file. The designations of d, h, or m may be used to specify the units of days, hours, minutes.Default is d.The position file names take the form: (see commands below)<root>.<outtype>.<site>.<start time>.<datatype>The resolution of the <start time> which is modulo the output interval depends on the output duration. For output durations greater than or equal to 1-day, the time is given as YYYYMMDD. For intervals shorter than 1-day, it is YYYYMMDD:HHMN. Minium output interval is 1 minute.SUM_FILE <root>Sets the root part of the summary file name. Using ? in the namewith be replaced by the -n string. (Default if command is not givenis trackRT or the -n string when -n used). File names are time taggedaccording to the pos_root output interval.CSV_ROOT <root>Set the root part of the name for comma separated values (CSV) outputfile. These files are used for AmCharts web plots.DCB_FILE <file name>Set the name of the data-code-bias (DCB) file. This file is part of theGAMIT tables directory and should be updated regularly. It is used to remove biases in the Melbourne-Wubbena widelanes. The receiver type can be specified with the RCV_TYPE or ANTE_OFF commands.ANTMOD_FILE <file name>Sets the name of a standard IGS antex file with phase center models forthe GPS ground antennas. Antenna types at specific sites are given withthe ANTE_OFF command. This command my be used multiple times for site specific model with new models replacing previously read ones.UPDATE_FILE <file name>Allows new trackRT commands to be issued during a run. Once the file is read it needs to be deleted before the trackRT will re-read it. File isonly read if it exists. NOTE: File should be removed before trackRT isrun or else it will be read when the command file is read (ie., it will overwrite the commands in the command file.OUT_TYPE <NEU+GEOD+XYZ+DHU>Specifies types of output coordinates. All types can specified in a single string with no spaces. The types areNEU -- North, East, Up differences from the reference site or from the coordinates given in the REF_NEU command.GEOD -- Geodetic latitude, longitude and height (in the GEOD format, the total atmosheric delay is given, while in the other formats theadjustment to the apriori delay is given).XYZ -- Cartersian XYZ coordinatesDHU -- Delta horizontal and Up coordinates from the apriori coordinates of each site (default output type)OUT_SIG_LIMIT <sigma (m)>Sets the maximum sigma of a position estimate for it to be output.If pseudorange data types are used, the default value of 1 m needsto increased to 10-100 meters.Analysis commands-----------------DATA_NOISE <L1 (m)> <L2 (m)> <P1 (m)> <P2 (m)> <Elev Weight> [PRN] Allows specificiation of the noise in the L1 phase, L2 phase,P1 range and P2 range, and the weight given to elevation angledepepence (at ver 1.20); variance is scale by (1+(W/sin(el))^2)where W is the <Elev Weight>. These values affect the sigmas printed for the position determinationsOptional: PRN may be added and noise assigned to that PRN (ifnon-PRN form is used, this will replace all PRN specific valuesso use the non-PRN first followed by specific PRN valuesDATA_TYPE <type>Specifies type of data to process. Types supported are L1 L2 LC P1 P1 PC which can be combined in a contiguous string. Files names include thisstring at the end.Vers 1.0: Only LC, LCPC and PC have been tested. Specifically L1 only data have not been fully implemented yet.USE_GPTGMF <Relative humidity (0-1)>Set the used of the GPT temperature pressure model and the GPT dry and wetmapping functions. Optional relative humity added vers 1.27; default 0.00. Default is the older MTT atmospheric model.Site and Antenna informationSite information is entered by first giving the command and then the information by site on the subsequent linesSITE_POSSite <X (m)> <Y (m)> <Z (m)> <Vx (m)> <Vy (m)> <Vz (m)><Epoch (yrs)>Site is the four character name of the site (more characters can be includedbut only the first 4 are checked). Site names that do not appeat in the listof sites to be processed are ignored). The remainder of the line contains positon and velocity and the epoch in deciminal years to which the position refers.ANTE_OFFSite <ARP dN (m)> <ARP dE (m)> <ARP dU (m)> <Antenna Name><Receiver Code>Specifies the type of antenna and its position of antena reference point (ARP) at each site. The antenna name including radome should be specified with theofficial IGS name for a standard ANTEX file or with a unique name that appearsin the ANTEX file for site specific calibrations. (Note: There is oneadditional character in the long antenna and this extra character before the radome name must be removed. The antenna name and radome can becopied directly from the rinex file if present. One more ANTEX files must be specified with the ANTMOD_FILE command for the antenna names to be useful. The receiver type DCB code can be optionally specified here as well (SeeRCV_TYPEcommand).RCV_TYPESite <Receiver code N/P/C>Specifies the type of data-code-bias (DCB) correction needed for the receiver. Code specifies the type of L1 and L2 ranges being measures. The choices are P -- Pcode, C -- C/A and N C/A with cross corelation for L2 range. The codes can be found in gamit/tables/rcvant.dat. These codes can also be given in the ante_off command. An up-to-date DCB_FILE command must be used to specify theDCB biases. The files are available from the MIT ftp site and update once per month.SITE_STATSSite <Apriori Sigmas in XYZ (m)> <RW noises in XYZ (m/sqrt(sec))> Gives statistics to assign to the kinematic station positions. The<Apriori Sigma in XYZ> are the three sigmas in XYZ for the initial postion and<RW noise in XYZ> are the three sigmas for the random walk in position. Units of the ransom walk are m/ssqrt(s).ALL can be used for the station name and the same statistics will be appliedto all kinematic sites (NOTE: the fixed site do not change position).e.g.site_statsall 1.0 1.0 1.0 0.010 0.010 0.010base 0.0 0.0 0.0 0.000 0.000 0.000(1 meters apriori sigmas and changes 10 mm/sqrt(s)=600 mm/sqrt(hr)).NOTE: One station should always set to zero sigmas and random walk otherwise allstation positions will be estimated and noramlly drift by large amounts.The fixed station does NOT need to be the reference site.ATM_STATSSite <Apriori Zenith delay sigma> <RW noise in Zenith delay> <RW dH/dt noise>Gives the statistics for the atmospheric delays by site. A random walk process's moise increases as sqrt(time) where time is time in seconds. The valuesare the initial sigma in meters, RW changes in meters per sqrt(second) and (addedverion 1.2) a dH/dt variance term so that during rapid height changes more process noise can be added to zenith delay estimate. The process noise variance is(<RW dH/dt noise)*abs(dh/dt)>^2 per sqrt(sec) where dh/dt is m/s. The default valeis 0.00023 m^2/sqrt(s). (Term only needed for aircraft processing).e.g.,atm_statsall 0.10 0.0001 0.00000t39a 0.10 0.0001 0.00023base 0.00 0.0000 0.00000Sets the apriori sigma as 10cm and allows the delay to change 6 mm/sqrt(hr) and 2.3 mm/sqrt(s) when height is changing at 10 m/s (fast ascent or desent) for the site t39a (aircraft). In the case, here "base" is a reference station at whichthe atmospheric delay is not estimated (other stations absorbe the atmospheric delays at the base station). With kinematic positioning over few hundred kmsiteseparations, the atmospheric delay at one site is normally fixed.Ambiguity resolution/data control commands.-------------------------------------------trackRT uses a combination of the Melbourne-Wubbena widelane (MW-WL), the extrawidelane (EX-WL) and the floating point estimates of the ionospheric free ambiguity (LC)to resolve integer ambiguities. If we denote the number of integer cycle ambiguitiesat L1 and L2 by N1 and N2, The MW-WL is an estimate of N1-N2 based on phase and rangedata; EX-WL = N1 - f1/f2 N2 and is an integer for L1 cycles, but 1.283 N2 for L2 cycles;LC = 2.546 N1 - 1.984 N2. The EX-WL is unaffected by geometric ranges changes, but doesdepend on the ionospheric delay. For short baselines, the EX-WL should be near zero forcorrect choices of N1 and N2. The LC residual should also be near zero when N1 and N2are correct. The problem in ambiguity resolution is that difference choices ofN1 andN2 can make different linear combinations small. For example, errors in N1 and N2 of3 and4 cycles will change LC by 0.298 cycles (56.6 mm), MW-WL by 1 cycle and EX-WL by2.132 cycles (405 mm). On long baselines, at low elevation angles, ionospheric delaysof 400 mm are common. The most common error is a N1=N2=1 cycle error.For thiscombination, the MW-WL is unaffected and LC changes by 0.562 cyc (107 mm) and EX-WL by0.283 cyc (54 mm). Even on relative short baselines, 54 mm ionospheric delays are commonThe follow commands are use to control the ambiguity resolution.AMB_SET <RelRank> <FloatSigma (2)> <MWWL Fact> <EXWL Fact> <Min AmbSig> <MaxChi>Sets parameters for ambiguity resolution. The input parameters are:<RelRank><WL min> Minumum of values need to allow bias fixing<WL avN> Maximum number to be used in computing sigma ofmean MW-WL<FloatSigma (2)> Minimum sigma for LC and MW-WL for ambiguity fixinf<MWWL Fact> Weighting factor for MW-WL in chi**2<EXWL Fact> Weighting factor for extra-wide lange<Min AmbSig> Minimum sigma to assigned to float estimates of ambiguities.<MaxChi> Max chi**2 value allowed for ambiguity to be resolved.EXWL_SET <Jump> <Min Sigma> <Scale> <Elev Fact><Jump> -- magnitude of jump in EX-WL to have cycle slip added (default 0.10)<Min Sigma> -- minimum sigma for mean ex-wl (cycles, default 0.02 cycles)<Scale> -- Scaling factor for length. Scale 0.1 results in 0.1 cyclesover 100 km (default)<Elev Factor> -- Elevation angle factor that increases sigma as (1 +factor/sin(elev))MWWL_SET <Jump> <Min Sigma> <Max Averaging number> <Min number><Jump> -- magnitude of jump in MW-WL to have cycle slip added (default 5.0)<Min Sigma> -- minimum sigma for mean MW-WL (cycles, default 0.10 cycles) <Max Averaging number> -- Maximum number of values to use to compute mean sigma<Min number> -- Minimum number needed to resolve ambiguityDD_SET <Jump (cycle)> <Min Number>Sets parameters for double difference processing.<Jump (cycle)> -- Magnitude of jump in double differences on bias fixed data thatwill introduce a cycle slip<Min Number> -- Minumum number of double differences for epoch to be processed.(If too few than errors in the data can be detected than this cancause large position errors, default of 4 double differences allowsredunancy).RMS_EDIT_TOL <n-sigma Tolerance> <min sigma> <Reset number><n-sigma tolerance> is an n-sigma condition where sigma is based on data noise model.<min sigma> Minimum phase sigma to use to that no phase residual less than <<min sigma>*<n-sigma tolerance> are deleted<Reset number> number of sequential delete data, before ambigity and cycle slipsare reset (assumed missed cycle slip).Testing and evaluting commandsSTATUS <type> <# epoch>Writes status information to the current summary file at <# epoch> intervals The types of reports are given by "type"P -- Parameter estimatesA -- Ambiguity resolution report (shows resolved and unresolved)W -- Widelanes (Melbourne-Wubbena and Extra-widelanes)R -- Postfit residuals are current epochC -- Report current A and W entries onlytype PAWR will output all reports.NUM_EPOCHS <number><number> is the number of epochs of data to be processed before stopping. When shortduraton files are output (POS_ROOT command) the number of epochs may not be reachedbecuase the epoch counter is reset with each new file.START_TIME <yy mm dd hh min sec.>Used to set the start time. Useful with trackRTr to have the processing start at thesame time as the realtime stream.DEBUG <Up to 10 epoch numbers in pairs>Sets which epochs will report detailed debug and status information (see STATUS commandas well). The pairs are used as:1,2 -- Mostly model information. Useful if data are generating large residuals3,4 -- Parameter estimates, widelanes, residuals5,6 -- One-way OMC and single differences7,8 -- Not Used9,10 -- When 9 is non-zero, antenna model and SP3 information.RESET <ALL/list of sites>Reset command. Resets the filter state vector and resolved ambiquities for a list of sites or ALL sites (generally used in the update_file to fix problems)EXCLUDE_SVS <list of PRN numbers to be excluded>Excludes satellites from being processed. Useful when a satellite is notin an SP3 file or not available during the times there are data.Example: exclude_svs 26 22 13Output types------------The output position files have two header lines that explain the columns in the file. The output in the summary and output files contain the following basic blocks are written to the summary file. The time in these files is given byan epoch number which is a count of the number of epochs of data (at the referencesite) since the start of the current data files. The output position files givethis counter and the corresponding GPS times and deciminal days.CSLIP line are cycle slip detectionsCSLIP Ep 540 Site THMG G 29 DMW-WL/Tol 4.31 4.00DEX-WL/Tol 1.78 0.20 cyc, Elev 10.00 degIndicates the jump in the MW-WL and EW-WL (and the tolerance for the detection)that caused the slip to be detected.CSLIP Ep 961 Site THMG G 12 DD RESID 2.18 Tol 0.50 cyc, Elev 63.85 degIndicates a slip added due to repeated double difference residuals.When ambiguities are resolved, three lines are output giving the statisticsAMBFIX WMAP PRN29 EP 599 RelRank 101.76 FC ------ dL12 4 3 Dchi 6.54 691.02 AMB 36AMBWLS WMAP PRN29 EP 599 RG 133 599 FX 3 iL12 0 0 Means -0.16 -0.05 RMS 0.74 0.06 # 467 eN12 -0.54 -0.38 AzEl 274.32 12.35AMBCON WMAP PRN29 EP 599 NCont 3 MW Res -0.159 0.185 Chi20.74 EX Res -0.051 0.068 Chi2 0.55 LC Res -0.046 0.020 Chi2 5.25 AMBFIX gives the epoch number, relative rank, the Fix Code (always ------ when fixed), thechanges to the number of L1 and L2 cycles from the intially asigned values, the changes in chi**2 for the best and second best choices of ambiguities and the number of the ambiguity resolved. When the status feature is used, ambiguities not resolved yet are listed as AMBFREE entries.AMBWLS gives the mean values and RMS scatters of the MW-WL and EX-WL (cycles).The iL12 estimate here is the changes in L1 and L2 cycles implied by justthe widelanes with no contribution from the LC estimates. On short baselines, these values are expected to be zero.AMBCON gives the contributiond to chi**2 from each of the input type. Ncont is the number of contributions (3 for LC data type). Values given arethe residal (Res) with its sigmas and Chi**2 (Chi2) contribution,When STATUS output is requested the following output types are possible. STATUS REPORT Epoch 1800 Date 2010 05 06 08 36 47.000 Type PAWR PARAMETER ESTIMATES Epoch 1800 for NP 16 Type P STATUS REPORT ...PARAMETER Ep 1800 NP 1 MSCG dATM -0.0466 +- 0.0001 m Var 0.681E-08PARAMETER Ep 1800 NP 2 PSDM dATM -0.0001 +- 0.0001 m Var 0.436E-08...PARAMETER Ep 1800 NP 7 LJRN G 27 1 0.1281 +- 0.0253 cyc Var 0.639E-03PARAMETER Ep 1800 NP 8 LJRN G 05 1 6.1172 +- 0.0011 cyc Var 0.112E-05Parameter values for position adjustments (XYZ), atmospheric delay adjustementsand float ambiguities are given.Ambiguity report for both both fixed and free ambiguities. ResEpoch is the epoch whenthe ambiquity was fixed (if fixed in previous data file block, may be greater than current count).AMBIGUITY Report Ep 1800 Number of ambiquities 59AMBFIXD # 11 MSCG PRN10 RelRank 616.27 FCode ------ dL12-1 -1 Dchi 1.01 776.94 Elev 72.02 deg; ResEpoch 156 AMBFIXD # 12 MSCG PRN27 RelRank 135.91 FCode ------ dL12-1 -1 Dchi 5.39 766.54 Elev 29.86 deg; ResEpoch 156..AMBFREE # 41 LJRN PRN27 RelRank 11.37 FCode RW---- dL120 0 Dchi 25.69 294.89 Elev 30.60 deg; ResEpoch 0 AMBFREE # 42 LJRN PRN05 RelRank 11.76 FCode R----C dL12 7 6 Dchi 106.15 1251.17 Elev 36.16 deg; ResEpoch 0The FCode gives the reason for not being fixed the entries are: RWSSSCR - Relative rankW - Widelane sigma too large or not enough data yet (see MWWWL_set command)S - LC sigma too largeS - L1 sigma too large when L1+L2 data type used (not tested)S - L2 sigma too large when L1+L2 data type used (not tested)C - Chi**2 increment too large.Widelane report shows status of average values of widelanes and current LC estimate. X means fixed, R ambiguity is still free.WIDELANE Report Ep 1800 Number of ambiquities 59WIDELANE 11 MSCG PRN10 EP 1800 Range 2 1800 # 1799 NC 3 X MW Res 0.06 0.10 Chi2 0.30 EX Res -0.18 0.19 Chi2 0.82 LC Res -0.01 0.02 Chi2 0.10...WIDELANE 19 MSCG PRN29 EP 1800 Range 421 1800 # 1380 NC 3 R MW Res -0.09 0.15 Chi2 0.35 EX Res 0.07 0.26 Chi2 0.07 LC Res 0.17 0.02 Chi2 74.30Postfit residuals for data types being used.POSTFIT RESIDUAL Report Ep 1800 Number of DD 76POST 1 Ep 1800 MSCG PRN10 - CIT1 PRN02 Res -0.0945 +- 0.0381 LC AMB 11 Elev 37.31 degPOST 2 Ep 1800 MSCG PRN27 - CIT1 PRN02 Res -0.1040 +- 0.0441 LC AMB 12 Elev 22.17 deg...POST 41 Ep 1800 MSCG PRN04 - CIT1 PRN02 Res 2.9712 +-1.8988 PC AMB 13 Elev 30.60 degPOST 42 Ep 1800 MSCG PRN30 - CIT1 PRN02 Res -0.2485 +-1.8116 PC AMB 14 Elev 36.16 deg..++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++EXAMPLE COMMAND FILE:--------------------* FILE trackRT_pbo.cmd* Example TrackRT command file to process some PBO sites.* Run BNC with as the TCP/IP port (login information* can be obtained from .) and select sites P496, P497,* P498, P505 to be cast on specific port (port 3765 for example).* Run trackRT as* trackRT -p 3765 -r p497 -d p496 p498 p505 -f trackRT_pbo.cmd -n P497 & * (if BNC is run on another machine used -m <machine name> as well.* A relative directory ../sp3_files is assumed to contain up to date* sp3 file.sp3_dir ../sp3_files* These positions are post April 4, Baja earthquake (extracted from* globk output with sh_exglk -f <globk org file> -apr <apr file>* entries from apr_file* Nore: Additional sites can be given heresite_posP498_GGU -2313650.15153 -4835550.43945 3444474.55150 -0.00658 0.02008 0.01478 2010.259 0.0033 0.0059 0.0042 -1.0000 -1.0000 -1.0000P497_GGU -2315937.74208 -4838712.09671 3438545.10828 -0.00945 0.01927 0.01884 2010.259 0.0035 0.0063 0.0044 -1.0000 -1.0000 -1.0000P496_GGU -2319722.10265 -4842511.39340 3430709.79068 -0.01697 0.02370 0.02587 2010.259 0.0050 0.0089 0.0062 -1.0000 -1.0000 -1.0000P503_GGU -2325056.31018 -4826758.67413 3449209.09091 -0.01167 0.01951 0.01660 2010.259 0.0033 0.0059 0.0042 -1.0000 -1.0000 -1.0000P505_AGU -2309739.55681 -4802072.50955 3493258.58085 -0.00377 0.00608 0.00160 2010.404 0.0033 0.0057 0.0040* PBO sites antenna and receiver information (all are NetRS and thus C type) * Again more sites can given). (Note: If antenna information extracted from* gamit , one space must be removed before radome string).* Antenna/Radom combination must be in antmod_file for model to be used).ante_offp475 0.00 0.00 0.0083 TRM29659.00 SCIT Cp066 0.00 0.00 0.0083 TRM29659.00 SCIT C p472 0.00 0.00 0.0083 TRM29659.00 SCIT C p478 0.00 0.00 0.0083 TRM29659.00 SCIT C p494 0.00 0.00 0.0083 TRM41249.00 SCIT C p496 0.00 0.00 0.0083 TRM29659.00 SCIT C p497 0.00 0.00 0.0083 TRM29659.00 SCIT C p498 0.00 0.00 0.0083 TRM41249.00 SCIT C p500 0.00 0.00 0.0083 TRM29659.00 SCIT C p503 0.00 0.00 0.0083 TRM29659.00 SCIT C p505 0.00 0.00 0.0083 TRM29659.00 SCIT C p510 0.00 0.00 0.0083 TRM29659.00 SCIT C* Change path as needed.antmod_file /home/tah/gg/tables/antmod.datdcb_file /home/tah/gg/tables/dcb.datdata_type LCPCsite_statsall 0.1 0.1 0.1 0.025 0.025 0.025p497 0.0 0.0 0.0 0.0 0.0 0.0atm_statsall 0.20 0.00010 0.000 ! Unit m/sqrt(sec) -> 0.0001 = 0.03 m/sqrt(day) p497 0.00 0.00000 0.000pos_root ? 1h# Used at MIT for web output. Results can be viewed at# /kmeduna/# csv_root /net/chandler/var//kmeduna/trackrt_view/data/P497# Useful to have this option. Note file should be removed before trackRT# is run and created when needed.update_file upd_app.cmd++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++INSTALLATION------------trackRT requires the BKG NTRIP Client (BNC) and QT libraries and include files. Ntrip is the Networked Transport of RTCM via Internet Protocol and the client program BNC provides access to the realtime data streams and casts the date via an internet port GNSS raw data in a structure that is decoded by trackRT. Currently trackRT only uses the GPS data streams. TheNTRIP system is discussedat igs.bkg.bund.de/ntrip/. The BNC program can be downloaded fromigs.bkg.bund.de/ntrip/download. The BNC documentation constains the instructions for obtaining the QT libraries and includes needed for itsinstallation from source files. trackRT only needs the executable version of BNC. The QT software can be obtained directly from Noia at the siteqt.nokia./downloadsFor trackRT, the libraries and include files are needed and these are containedin Frameworks libraries 4.6.2 for your platform (linux or Mac). Earlier version should also work. For Linux and Max OSX, QT can be downloaded directly from qt.nokia..qt.nokia./downloads/linux-x11-cppqt.nokia./downloads/qt-for-open-source-cpp-development-on-mac-os-xQT takes several hours to install.On linux systems, other tools that can be needed are:sudo app-get install gfortransudo app-get install g++sudo app-get build-dep qt4-qmake(The last app-get may be need if the qt install can not find X11 libraries and includes.)For MacOSXThe basic qt install steps aresudo port install qt4-macIf you don't have MacPorts, it's a very easy way to get lots of *nix toolsfor the Mac. You can download MacPorts from .Once QT is installed, the makefile in the trackRT directory should be modifiedto include the directories for the libraries and includes. Example entries are given in the makefile. The actual entries are for the MIT QT installation (which is not quite standard).NOTE: There is a file makefile (not Makefile) for making trackRT series of programs.To run trackRT, bnc is run first and the desired stations added the mountpoints. Once bnc streaming is started, trackRT is started with the machine name (if different from the bnc machine) and port (given in the bnc FEED tab) given. From。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
解压
# tar xvzf .......
1、# yum install tcsh
# yum install gcc
# yum install gcc-gfortran
# yum install gmt
2、更改csh路径
# chsh
/bin/csh
#echo $PATH(测试路径是否正确)
3、安装前将instal_software中的libx11_a修改为libx11_so
# vim install_software (查找到libx11_a)
将stdio.h中的getline改为getline1(stdio.h显示的路径指令：
# whereis stdio.h)
4、# vim Makefile.config
将其中的-------for linux from 0.0.1 to 3.0.0------那一行下面的内容改掉
改成：OS_ID Linux 0001 3107
5 # ./install_software
安装成功。
6、配置csh.cshrc
输入指令 # vim /etc/csh.cshrc
# /etc/csh.cshrc: system-wide .cshrc file for csh(1) and tcsh(1)
# set prompt
set prompt = "[`echo $USER@$cwd`]% "
set history = 50
#set system path
set path = (. /bin /sbin /etc /usr/bin /usr/sbin /usr/bin/X11 /usr/X11R6/bin)
set path = (/usr/local/bin /usr/local/sbin /usr/local/lib /usr/lib $path)
setenv PATH "/usr/sbin:/sbin:${PATH}"
setenv MANPATH /usr/local/man:/usr/X11R6/man
setenv LIBPATH /usr/lib:/usr/local/lib
#set Gamit/Globk path
set gamitpath = /home/czf/share/gamit
set path = ($gamitpath/gamit/bin $gamitpath/kf/bin $gamitpath/com $path)
setenv HELP_DIR $gamitpath/help/
setenv _POSIX2_VERSION 199209
#set GMT path
set gmtpath = /usr/lib/gmt
set path = ($gmtpath/bin $gmtpath/include $gmtpath/lib $gmtpath/man $gmtpath/share
$gmtpath/www $path)
set path = (/etc/gmt /usr/include/gmt /usr/lib/gmt /usr/share/gmt /usr/share/doc/gmt $path)
setenv MANPATH "$gmtpath/man:$MANPATH"
#set full domain for GAMIT shell-scripts used to transfer data (sh_get_nav, sh_get_rinex,
#sh_get_orbits, and
# sh_get_hfiles)
set host = `\hostname | \awk -F. '{print $1}'`
set dom = `\hostname | \awk -F. '{print $2}'`
if( $dom == '' ) then
alias hostname 'echo $host.mit.edu'
endif
# set alias
# alias cd 'cd \!*; set prompt= '[`echo $USER@$cwd`]%''
alias pwd 'echo $cwd'
alias ls 'ls --classify'
alias ll 'ls -l'
alias la 'ls -a'
alias rm 'rm -i'
alias mv 'mv -i'
alias cls 'clear'