GAMIT软件Track模块使用帮助.doc

1、.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 % trackRTB where options are-m Host name of system with real

2、-time data port(Normally BNC would be running on this system)(Default local host)-p Port number on host supplying the data stream(Mandatory input).-f trackRT command file (similar to track command file)-r 4-char code for reference site.-d list of sites to be processed.If notspecified all sites in da

3、ta 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 Specifies a root name to which output will be directedWord 文档.(.out is appended for stdout output).Use ? in thepos_root, sum_file

4、 and csv_file names the ? will bereplaced by this string.Command file name must be given.% trackRTr Same list of option although -m, -p have no effect -d 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 tha

5、t 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 BNC 2.5 and greater that is now outp

6、utingascii rather than binary data (trackRTB handles BNC2.0) (110501)1.00- Initial version 100301 (2010 March 1)Word 文档.INTRODUCTION:-TrackRT/TrackRTB and trackRTr have many commands that allow it be controlledboth when started and during runs.Below the detailed commands are givenand here we give th

7、e 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 a

8、nd process noise. Typical site_statswould be where cit1 here is the reference site.site_statsall0.1 0.1 0.10.025 0.025 0.025cit10.0 0.0 0.00.0 0.0 0.0ATM_STATS: By default atmospheric delays are not estimated, and normallythese should be.Normally the reference site is set to zerowhen site separation

9、s are less the 500-1000km. Typical atm_statscommand would beatm_statsall0.20 0.00010 0.000! Unit m/sqrt(sec) - 0.0001 = 0.03 m/sqrt(day)cit10.00 0.00000 0.000Word 文档.If mixed antenna and receiver types are usedANTE_OFF: To specify the antenna heights and types and receiver typesANTMOD_FILE: Must be

10、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 controlle

11、d 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 identicalexcept trackRTB is expecting a binary data stream.COMMANDS-Input

12、/Output commands-Word 文档.SP3_DIR - Directory where sp3 files are stored- Center for orbits (default igs, igr and igu also tested)POS_ROOT Set the root part of the name for the output files.Root part of name.When ? included in the , the ? isreplaced with the -n string. Duration of data in each file.T

13、he 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).The resolution of the which is modulo the output intervaldepends on the output duration.For output durations greater than orequal to 1-day, t

14、he time is given as YYYYMMDD.For intervals shorterthan 1-day, it is YYYYMMDD:HHMN.Minium output interval is 1 minute.SUM_FILE 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).Fil

15、e names are time taggedWord 文档.according to the pos_root output interval.CSV_ROOT Set the root part of the name for comma separated values (CSV) outputfile.These files are used for AmCharts web plots.DCB_FILE Set the name of the data-code-bias (DCB) file.This file is part of theGAMIT tables director

16、y and should be updated regularly.It is used toremove biases in the Melbourne-Wubbena widelanes.The receiver type canbe specified with the RCV_TYPE or ANTE_OFF commands.ANTMOD_FILE Sets the name of a standard IGS antex file with phase center models forthe GPS ground antennas.Antenna types at specifi

17、c sites are given withthe ANTE_OFF command. This command my be used multiple times for sitespecific model with new models replacing previously read ones.UPDATE_FILE Allows new trackRT commands to be issued during a run.Once the file isread it needs to be deleted before the trackRT will re-read it.Fi

18、le isonly read if it exists.NOTE: File should be removed before trackRT isWord 文档.run or else it will be read when the command file is read (ie., it willoverwrite the commands in the command file.OUT_TYPE Specifies types of output coordinates.All types can specified in asingle string with no spaces.

19、The types areNEU- North, East, Up differences from the reference site or from thecoordinates given in the REF_NEU command.GEOD- Geodetic latitude, longitude and height (in the GEOD format, thetotal atmosheric delay is given, while in the other formats theadjustment to the apriori delay is given).XYZ

20、- Cartersian XYZ coordinatesDHU- Delta horizontal and Up coordinates from the apriori coordinatesof each site (default output type)OUT_SIG_LIMIT 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

21、meters.Analysis commandsWord 文档.-DATA_NOISE PRNAllows 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 .These values affect the sigmas printed forthe positio

22、n 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 Specifies type of data to process.Types supported are L1 L2 LC P1 P1 PCwhich can be combined i

23、n a contiguous string. Files names include thisstring 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 Set the used of the GPT temperature pressure model and the GPT dry and wetWord 文档.mapping functions. Optional relat

24、ive 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 theinformation by site on the subsequent linesSITE_POSSite Site is the four character name of the site (more characters can

25、 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 containspositon and velocity and the epoch in deciminal years to which the positionrefers.ANTE_OFFSite Specifies the type of antenna and its position

26、 of antena reference point (ARP)at each site.The antenna name including radome should be specified with theWord 文档.official 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 station.inf

27、o long antenna and this extra characterbefore the radome name must be removed.The antenna name and radomecan 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 option

28、ally specified here as well (SeeRCV_TYPEcommand).RCV_TYPESiteSpecifies 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 i

29、n gamit/tables/rcvant.dat.These codes can also be given in theante_off command.An up-to-date DCB_FILE command must be used tospecify theDCB biases. The files are available from the MIT ftp site and update once permonth.Word 文档.SITE_STATSSiteGives statistics to assign to the kinematic station positio

30、ns.The are the three sigmas in XYZ for the initial postion and 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 appliedto all kinematic sites (NOTE: the fixed site do not change position).

31、e.g.site_statsall1.01.0 1.00.010 0.010 0.010base0.00.0 0.00.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 otherwiseallstation positions will be estimated and noramlly drift by large amounts.The fix

32、ed station does NOT need to be the reference site.ATM_STATSSite Word 文档.Gives the statistics for the atmospheric delays by site. A random walk processsmoise 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(adde

33、dverion 1.2) a dH/dt variance term so that during rapid height changes moreprocess noise can be added to zenith delay estimate.The process noisevariance is(2 per sqrt(sec) where dh/dt is m/s.The defaultvaleis 0.00023 m2/sqrt(s). (Term only needed for aircraft processing).e.g.,atm_statsall0.100.00010

34、.00000t39a0.100.00010.00023base0.000.00000.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 de

35、lay is not estimated (other stations absorbe the atmosphericdelays at the base station).With kinematic positioning over few hundred kmsiteseparations, the atmospheric delay at one site is normally fixed.Word 文档.Ambiguity resolution/data control commands.-trackRT uses a combination of the Melbourne-W

36、ubbena widelane (MW-WL), theextrawidelane (EX-WL) and the floating point estimates of the ionospheric freeambiguity (LC)to resolve integer ambiguities.If we denote the number of integer cycleambiguitiesat L1 and L2 by N1 and N2,The MW-WL is an estimate of N1-N2 based onphase and rangedata; EX-WL = N

37、1 - f1/f2 N2 and is an integer for L1 cycles, but 1.283 N2 for L2cycles;LC = 2.546 N1 - 1.984 N2.The EX-WL is unaffected by geometric rangeschanges, but doesdepend on the ionospheric delay.For short baselines, the EX-WL should benear zero forcorrect choices of N1 and N2.The LC residual should also b

38、e near zero whenN1 and N2are correct.The problem in ambiguity resolution is that difference choices ofN1 andWord 文档.N2 can make different linear combinations small.For example, errors in N1 andN2 of3 and 4 cycles will change LC by 0.298 cycles (56.6 mm), MW-WL by 1 cycle andEX-WL by2.132 cycles (405

39、 mm).On long baselines, at low elevation angles, ionosphericdelaysof 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 del

40、aysare commonThe follow commands are use to control the ambiguity resolution.AMB_SET Sets parameters for ambiguity resolution.The input parameters are:Minumum of values need to allow bias fixingMaximum number to be used in computing sigma ofmean MW-WLMinimum sigma for LC and MW-WL for ambiguity fixi

41、nfWord 文档.Weighting factor for MW-WL in chi*2Weighting factor for extra-wide langeMinimum sigma to assigned to float estimates ofambiguities.Max chi*2 value allowed for ambiguity to be resolved.EXWL_SET - magnitude of jump in EX-WL to have cycle slip added (default0.10) - minimum sigma for mean ex-w

42、l (cycles, default 0.02 cycles)- Scaling factor for length. Scale 0.1 results in 0.1 cyclesover 100 km (default) - Elevation angle factor that increases sigma as (1 +factor/sin(elev)MWWL_SET - magnitude of jump in MW-WL to have cycle slip added(default 5.0) - minimum sigma for mean MW-WL (cycles, de

43、fault 0.10 cycles) - Maximum number of values to use to computemean sigmaWord 文档. - Minimum number needed to resolve ambiguityDD_SET Sets parameters for double difference processing. - Magnitude of jump in double differences on bias fixed datathatwill introduce a cycle slip- Minumum number of double

44、 differences for epoch to beprocessed.(If too few than errors in the data can be detected than this cancause large position errors, default of 4 double differencesallowsredunancy).RMS_EDIT_TOL is an n-sigma condition where sigma is based on datanoise model. Minimum phase sigma to use to that no phas

45、e residual less than* are deleted number of sequential delete data, before ambigity and cycleslipsare reset (assumed missed cycle slip).Word 文档.Testing and evaluting commandsSTATUS Writes status information to the current summary file at intervalsThe types of reports are given by typeP - Parameter estimatesA - Ambiguity resolution report (shows resolved and unresolved)W