Tag Archives: tutorial GAMIT

Tutorial GAMIT part 2

Posted on by
14

This article is continuing on GPS data processing using GAMIT. For reader who want to read previous article, please read this.

Assumed that you already done with your data processing using this step and you have at least more than 1 GPS data observation to combine. To make it clear, I will draw the assumptions as a picture below

Fig.1. GAMIT data processing

The main specific control for this step are :

  1. The list of sites to be used in defining the reference frame.
  2. The a priori coordinates for these sites.

for simplicity, user can use IGS station as sites for defining reference frame. It’s due to the assumption that IGS station is most stable station used in the processing. Using this step, reader will gain output of GAMIT/GLOBK processing which contain GPS observation sites and GPS reference sites, and get 3 epoch repeatibilities and velocities.

A. Create new folder.

create 2 new folder under your project folder, namely “h-files” and “h-gamit”. This folder is used to store h-files from global solution (under hfiles folder) and h-files from sh_gamit processing (under hgamit folder).

hint: to make it easy, you can link your h-files from sh_gamit processing which is originally stored in each DOY folder using this command (executed under hgamit folder)

ln -s ../???/h[expt] .

where : expt is the name of the project experiment.

B. Create list of H-file.

Create a list of H-files to be input to GLOBK for repeatabilities. Before that, convert H-file to glb file using this command.

htoglb glbf tables/[expt].svs hgamit/h* for your h-files from gamit processing, and

htoglb glbf tables/[expt].svs hfiles/h* for h-files from global solution.

After that, create list of H-files (under glbf folder) using this command :

ls ../glbf/*.glx > [expt].gdl

where :

expt = name of experiment.

the *.gdl file may optionally contain additional parameters to indicate re-weighting and coupling. To do that, you can append “1.0 +” in the end of every line inside the *.gdl file.

C. Preparing control file.

before execute repeatabilities, don’t forget to copy globk_comb.cmd, glorg_comb.cmd and *.gdl file to GSOLN, and edit those 2 files (globk_com.cmd and glorg_comb.cmd). please refers to GAMIT/GLOBK manual book to help you to edit these files.

In addition, please copy stab_site.global, itrf08.apr, igs08_disc.eq, and pmu.usno from gg/tables to your project directory tables. this step may be vary, because these copies files are depend on your editing in globk_com.cmd and glorg_com.cmd)

D. Processing

To get repeatabilities, please type on your terminal under your project directory

glred 6 globk_com.prt globk_comb.log [expt].gdl globk.cmd

where = [expt] is the name of your *.gdl file.

This command will make globk_comb.org as output file solution for h-file in the *.gdl file. You should check this file for quality and uncertainties to the positions and velocities.

E. Plot result

Use this command for plotting your globk output

sh_plotcrd -f [file *.org] -expt [experiment name] -cols [number of column, type = 1 for single column or 2 for double colomn]

check your ps file to see the output plot of this command. If there are any outliers, they must be dealt with by exclusion or downweighting .

 

Kyoto, 22 Agustus 2014

Hidayat Panuntun

Tutorial GAMIT part 1

Posted on by
90

This example is discuss about GPS data processing using GAMIT/GLOBK

The example is set up to :

A. Conduct phase processing (sh_gamit).

B. Compute and plot daily repeatabilities (sh_glred).

So please note that this tutorial is only talked about phase processing using sh_gamit and compute daily repeatability using sh_glred, this tutorial is not (yet) talked about combine daily result using GLOBK processing.

As I said, using this tutorial will yield daily position of your data processing. Even your input rinex is multiple day, the command will generate daily .org file. Then, you can combine this daily result on one file repeatability using GLOBK (will be publish soon) (^_^)v.

enjoy!

–Begin Processing–

1. CREATE PROJECT DIRECTORY FOLDER.

You can done this work manually or by terminal typing on the path where you will work on GAMIT (e.g= Using, mkdir [your project]). Please note that the name of your project must be named with only 4 digit.

2. CREATE GPS DATA, IGS, and BRDC FOLDER.

Just same as before, just create folder and use “rinex”, “igs”, “brdc” as the name for the folder.

a. Rinex data.

Store your GPS data observation in this folder. If you wanna use RINEX from IGS station, you can download easily using GAMIT scritp named “sh_get_rinex”. just type in your terminal “sh_get_rinex -archive sopac cddis -yr [years] -doy [doy] -ndays [days] -sites [IGS sites]”

Noted:

  • years: 4 char year of nav data requested
  • doy: 3 char day of year of nav data requested
  • days: Number of consecutive days of data to retrieve
  • IGS sites: List of sites to be retrieved from the ftp archive (IGS websites)

3. INPUT DATA
firstly, set up GAMIT data environment using “sh_setup” command. Type the command on the terminal under your project folder.
e.g: “sh_setup -yr [year] -apr [itrf]”

Noted:

  • year: 4 char year of observation
  • itrf: ITRF that will used on processing

4. ESTABLISH INPUT FILE

Input files provided in advance are lfile, sistbl, sestbl, station.info, process default. Copy and edit those file from gg/tables to directory tables under your project folder.

  • lfile: input the coordinates apriory as initial coordinates (user can use coordinates from RINEX header).
  • sittbl: listed IGS stations used by GAMIT to compute loose constrained solution.
  • sestbl: list and setting of parameters used on GAMIT computation
  • station.info: contain all information related to station that will used on the processing.
  • process.default: file which contain the default process setting.
  • sites.default: file used for controlling the use of station in the processing.
*User is strongly recommended to read GAMIT/GLOBK manual book to make modification of those input file.

 5. AUTOMATIC BATCH PROCESSING (SH_GAMIT)

type

sh_gamit -s [year] d1 d2 -expt [expt] -pres [ress]

e.g: sh_gamit –s 2011 164 167 –expt regi -pres ELEV

noted:

  • -s d1 d2: where d1 is start day and d2 is stop day to be processed. e.g: 2011 164 167
  • -expt: 4 char name of experiment/s to run (check on your input file)
  • -pres: Plot gamit residuals as skyplots. Option: NO / YES (skyplot only) / ELEV (skyplot + phase vs elevation [Default N].

6. SH_GAMIT OUTPUT

Check the summary for number of stations, Postfit RMS , Postfit nrms, ambiguity resolution, and coordinate adjustments in your <doy> folder under your project directory. Look at the sky plots and phase vs elevation angle plots in the /gifs directory. If necessary, compare the q-files and sky plots with those in the check_files directory (Note that the sky plots in check_files are still in postscript, whereas those in /gifs have been converted to gif images.)

fig 1. phase vs elevation angle plot

fig 2. Sky plot

7. SH_GLRED

before start to use sh_glred, user must get template command files into /gsoln by typing (on the terminal) at the project level “sh_glred -cmd”

go to /gsoln directory under your project directory, edit globk_comb.cmd and glorg_comb.cmd to use translation-only stabilization, and to use for stablization the apr file and  stab_site list pre-generated for the example.

In globk_comb.cmd, comment out (add a character to column 1, user can put ‘x’ character)

apr_wob 10 10 0 0
apr_ut1 10 0

and uncomment (remove the ‘x’)

x apr_wob .25 .25 .1 .1
x apr ut1 .25 .1

In glorg_comb.cmd, comment out

pos_org xtran ytran ztran xrot yrot zrot
source ../tables/[stab_site.file]

and uncomment (remove the ‘x’)

x pos_org xtran ytran ztran
x source ../../tables/regional_stab_site

*please note, the user must copy [stab_site.file] from gg/tables to tables under the project directory. User can choose stab_site file that will use on the processing. Just make sure for change [stab_site.file] into stab_site name that will used in the processing.

After establishing path and other important file, then type on the terminal:

sh_glred -d [yr day] -s [yr1 d1 yr2 d2]  -expt [expt] -opt H G E LA

-s <yr1 d1 yr2 d2>:   start and stop dates, e.g. -s 2011 164 2011 167

-expt: 4 char name of experiment/s to run.

-opt: Optional input to set some output paramater.

H: Run htoglb on all ascii files present or linked within glfpth (usually procdir/glbf).

G: Run glred for combination or repeatabilities.

E: Run ensum and sh_baseline for plots

LA: call h-files from global global that locally store on your local directory.

*noted: for advance setting, please read GAMIT/GLOBK manual book.

8. SH_GLRED OUTPUT

The script as commanded will translate the GAMIT ascii h-files in each day directory to GLOBK binary h-files and put them into the /glbf directory ( H ); create a gdl file for each day listing the h-file; run GLOBK for each day ( G ) using globk_comb.cmd and glorg_comb.cmd ( G ); and generate time series plots ( E, from program ensum).

In /gsoln, use Ghostscript (‘gs’) to view the time series, which should show wrms repeatabilities at the level of 1 mm horizontal and 2 mm vertical, and nrms values < 1.0. If these are good, there is usually no need to inspect any other files.

Please don’t be hesitate to ask me if you face some difficulties on the processing.

-to be continued (insyaallah, for GLOBK processing).

Kyoto, 3 May 2014

Hidayat Panuntun