# gpsgridder¶

Interpolate GPS strains using Green’s functions for elastic deformation

## Synopsis¶

**gmt gpsgridder** [ *table* ]
**-G***outfile*
[ **-I***increment* ]
[ **-R***region* ]
[ **-C**[**n**]*value*[**+f***file*] ]
[ **-E**[*misfitfile*] ]
[ **-F**[**d**|**f**]*fudge* ]
[ **-L** ]
[ **-N***nodefile* ]
[ **-S***nu* ]
[ **-T***maskgrid* ]
[ **-V**[*level*] ]
[ **-W**[**w**] ]
[ **-b**binary ]
[ **-d**nodata ]
[ **-e**regexp ]
[ **-f**flags ]
[ **-h**headers ]
[ **-o**flags ]
[ **-qi**flags ]
[ **-x**[[-]*n*] ]
[ **-:**[**i**|**o**] ]
[ **--PAR**=*value* ]

**Note:** No space is allowed between the option flag and the associated arguments.

## Description¶

**gpsgridder** grids 2-D vector data such as GPS velocities by using a coupled
model based on 2-D elasticity. The degree of coupling can be tuned by adjusting
the effective Poisson’s ratio. The solution field can be tuned to extremes such
as incompressible (1), typical elastic (0.5) or even an unphysical value of -1
that basically removes the elastic coupling of vector interpolation. Smoothing
is offered via the optional elimination of small eigenvalues.

## Required Arguments¶

*table*- table with GPS strain rates at discrete locations. We expect the input
format to be
*x y u v*[*du dv*] (see**-W**to specify data uncertainties or weights). If*lon lat*is given you must supply**-fg**and we will use a flat Earth approximation in the calculation of distances.

**-G***outfile*- Name of resulting output file. (1) If options
**-R**,**-I**, and possibly**-r**are set we produce two equidistant output grids. In this case,*outfile*must be a name template containing the C format specifier %s, which will be replaced with u and v, respectively. (2) If option**-T**is selected then**-R**,**-I**cannot be given as the*maskgrid*determines the region and increments. Again, the*outfile*must be a name template for the two output grids. (3) If**-N**is selected then the output is a single ASCII (or binary; see**-bo**) table written to*outfile*; if**-G**is not given then this table is written to standard output. The**-G**option is ignored if**-C**or**-C**0 is given.

## Optional Arguments¶

**-C**[**n**]*value*[**+f***file*]- Find an approximate surface fit: Solve the linear system for the
spline coefficients by SVD and eliminate the contribution from all
eigenvalues whose ratio to the largest eigenvalue is less than
*value*[Default uses Gauss-Jordan elimination to solve the linear system and fit the data exactly]. If*value*is in 0–1 range the we assume it is the fraction of eigenvalues to keep. Optionally, append**+f***file*to save the eigenvalue ratios to the specified file for further analysis. If a negative*value*is given then**+f***file*is required and execution will stop after saving the eigenvalues, i.e., no surface output is produced. Specify**-Cn***value*to retain only the*value*largest eigenvalues.**Note**: 1/4 of the total number of data constraints is a good starting point for further experiments.

**-E**[*misfitfile*]- Evaluate the spline exactly at the input data locations and report
statistics of the misfit (mean, standard deviation, and rms) for
*u*and*v*separately and combined. Optionally, append a filename and we will write the data table, augmented by two extra columns after each of the*u*and*v*columns holding the spline estimates and misfits. If**-W**is given we also add two more columns with the chi^2 values.

**-F**[**d**|**f**]*fudge*- The Green’s functions are proportional to terms like 1/r^2 and log(r)
and thus blow up for r == 0. To prevent that we offer two fudging schemes:
**-Fd***del_radius*lets you add a constant offset to all radii and must be specified in the user units. Alternatively, use**-Ff***factor*which will compute*del_radius*from the product of the shortest inter-point distance and*factor*[0.01].

**-I***xinc*[**+e**|**n**][/*yinc*[**+e**|**n**]]*x_inc*[and optionally*y_inc*] is the grid spacing. Optionally, append a suffix modifier.**Geographical (degrees) coordinates**: Append**m**to indicate arc minutes or**s**to indicate arc seconds. If one of the units**e**,**f**,**k**,**M**,**n**or**u**is appended instead, the increment is assumed to be given in meter, foot, km, Mile, nautical mile or US survey foot, respectively, and will be converted to the equivalent degrees longitude at the middle latitude of the region (the conversion depends on PROJ_ELLIPSOID). If*y_inc*is given but set to 0 it will be reset equal to*x_inc*; otherwise it will be converted to degrees latitude.**All coordinates**: If**+e**is appended then the corresponding max*x*(*east*) or*y*(*north*) may be slightly adjusted to fit exactly the given increment [by default the increment may be adjusted slightly to fit the given domain]. Finally, instead of giving an increment you may specify the*number of nodes*desired by appending**+n**to the supplied integer argument; the increment is then recalculated from the number of nodes and the domain. The resulting increment value depends on whether you have selected a gridline-registered or pixel-registered grid; see GMT File Formats for details.**Note**: If**-R***grdfile*is used then the grid spacing (and registration) have already been initialized; use**-I**(and**-r**) to override the values.

**-L**- Leave trend alone. Do
*not*remove a planer (2-D) trend from the data before fitting the spline. [Default removes least squares plane, fits normalized residuals, and restores plane].

**-N***nodefile*- ASCII file with coordinates of desired output locations
**x**in the first column(s). The resulting*w*values are appended to each record and written to the file given in**-G**[or stdout if not specified]; see**-bo**for binary output instead. This option eliminates the need to specify options**-R**,**-I**, and**-r**.

**-R***xmin*/*xmax*/*ymin*/*ymax*[**+r**][**+u***unit*] (more …)- Specify the region of interest.

**-S***nu*- Specify Poisson’s ratio to use for this 2-D elastic sheet [0.5].
**Note**: 1.0 is incompressible in a 2-D formulation while -1 removes all coupling between the two directions.

**-T***maskgrid*- Only evaluate the solutions at the nodes
in the
*maskgrid*that are not set to NaN. This option eliminates the need to specify options**-R**,**-I**(and**-r**).

**-W**[**w**]- One-sigma data uncertainties for
*u*and*v*are provided in the last two columns. We then compute least squares weights that are inversely proportional to the square of the uncertainties. Append**w**if weights are given instead of uncertainties, in which case we just use the weights as provided (no squaring). This results in a weighted least squares fit. Note that**-W**only has an effect if**-C**is used. [Default uses no weights or uncertainties].

**-V**[*level*] (more …)- Select verbosity level [
**w**].

**-d**[**i**|**o**]*nodata*(more …)- Replace input columns that equal
*nodata*with NaN and do the reverse on output.

**-e**[**~**]*”pattern”*|**-e**[**~**]/*regexp*/[**i**] (more …)- Only accept data records that match the given pattern.
**-fg**- Geographic grids (dimensions of longitude, latitude) will be converted to meters via a “Flat Earth” approximation using the current ellipsoid parameters.

**-h**[**i**|**o**][*n*][**+c**][**+d**][**+m***segheader*][**+r***remark*][**+t***title*] (more …)- Skip or produce header record(s). Not used with binary data.

**-i***cols*[**+l**][**+s***scale*][**+o***offset*][,*…*][,**t**[*word*]] (more …)- Select input columns and transformations (0 is first column,
**t**is trailing text, append*word*to read one word only).

**-qi**[~]*rows*[**+c***col*][**+a**|**f**|**s**] (more …)- Select input rows or data range(s) [all].

**-r**[**g**|**p**] (more …)- Set node registration [gridline].

**-:**[**i**|**o**] (more …)- Swap 1st and 2nd column on input and/or output.

**-^**or just**-**- Print a short message about the syntax of the command, then exit (NOTE: on Windows just use
**-**). **-+**or just**+**- Print an extensive usage (help) message, including the explanation of any module-specific option (but not the GMT common options), then exit.
**-?**or no arguments- Print a complete usage (help) message, including the explanation of all options, then exit.
**--PAR**=*value*- Temporarily override a GMT default setting; repeatable. See gmt.conf for parameters.

## Units¶

For map distance unit, append *unit* **d** for arc degree, **m** for arc
minute, and **s** for arc second, or **e** for meter [Default], **f**
for foot, **k** for km, **M** for statute mile, **n** for nautical mile,
and **u** for US survey foot. By default we compute such distances using
a spherical approximation with great circles (**-jg**) using the authalic radius
(see PROJ_MEAN_RADIUS). You can use **-jf** to perform
“Flat Earth” calculations (quicker but less accurate) or **-je** to perform
exact geodesic calculations (slower but more accurate; see
PROJ_GEODESIC for method used).

## Examples¶

To compute the *u* and *v* strain rate grids from the GPS data set *gps.txt*,
containing *x y u v du dv*, on a 2x2 arc minute grid for California, try

gmt gpsgridder gps.txt -R-125/-114/31/41 -I2m -fg -W -r -Ggps_strain_%s.nc -V

## References¶

Haines, A. J. et al., 2015, *Enhanced Surface Imaging of Crustal Deformation*, SpringerBriefs in Earth Sciences,
doi:10.1007/978-3-319-21578-5_2.

Sandwell, D. T. and P. Wessel, 2016, Interpolation of 2-D Vector Data Using Constraints from Elasticity,
*Geophys. Res. Lett., 43*, 10,703-10,709,
http://dx.doi.org/10.1002/2016GL070340