Home - VicenteYanez/GFA GitHub Wiki

Introduction (in construction)

My name is Vicente Yañez and I am a geologist interested in computer science. GNSS Field Analysis (GFA) is a compilation of python functions that I wrote for my internship in Centro Sismológico Nacional (CSN) and during my undergraduate thesis (spanish).

This software has the functions to calculate vectors from a GNSS time series data (from the trajectory model of Bevis & Brown, 2014 or from a simple lineal fitting), to obtain a velocity field that represent the cortical deformation. Then, from the velocity field this program include the tools to calculate the velocity tensor, and from this, the vorticity and stretching tensors (from the equation of Davis & Titus, 2011).

The Python Modules requires

GFA depends on Numpy, scipy, matplotlib and optionaly, cartopy.

Structure

GFA is divided in two principal sections. The GNSS Analysis and the Field Analysis. The first is focused on the calculation of the velocity vector and the second one, use the output of GNSS Analysis to analyze the velocity field in different ways.

GNSS Analysis

Classes

  • TimeSeriesControl
  • ModelControl
  • ModeloTrayectoria

Functions

  • loadGPS
  • lineal_model
  • hvsd
  • fun_vector

Field Analysis

Functions

  • field
  • geometry

Individual Scripts

Finally, GFA includes various scripts that use this classes and functions. Also, GFA includes examples of use, and its results.You can see the complete list writing gfa --help on bash.

Scripts

  • select: Cut and save the time series giving a specific time and location area.
  • buildmodelall: Calculate the trajectory model for all the stations in the db using the same parameters for all of them.
  • build
  • plot
  • vector

Examples

  • example_vorticity: Calculate the vorticity in the southern Andes from velocity vectors between the years 2007-2010.

To do list...

  • To improve the documentation and translate the doc string in spanish.
  • To create a user interface for the field analysis
  • To calculate geostadistical information from the GNSS station and their time series
  • To relate the station with geodatabase information like, rock type, morphology, distance to the trench, etc. The objective is display the data when you are calculating the model.