Generic Mapping Tools GMT(1)
Total Page:16
File Type:pdf, Size:1020Kb
GMT(1) Generic Mapping Tools GMT(1) NAME GMT − The Generic Mapping Tools data processing and display software package INTRODUCTION GMT is a collection of public-domain Unix tools that allows you to manipulate x,y and x,y,z data sets (fil- tering, trend fitting, gridding, projecting, etc.) and produce PostScript illustrations ranging from simple x-y plots, via contour maps, to artificially illuminated surfaces and 3-D perspective views in black/white or full color.Linear,log10, and power scaling is supported in addition to 25 common map projections. The pro- cessing and display routines within GMT are completely general and will handle any(x,y) or (x,y,z) data as input. SYNOPSIS GMT is also a wrapper script that can start anyofthe programs: GMT module module-options where module is the name of a GMT program and the options are those that pertain to that particular pro- gram. GMT OVERVIEW The following is a summary of all the programs supplied with GMT and a very short description of their purpose. Detailed information about each program can be found in the separate manual pages. blockmean L2 (x,y,z) data filter/decimator blockmedian L1 (x,y,z) data filter/decimator blockmode Mode (x,y,z) data filter/decimator filter1d Filter 1-D data sets (time series) fitcircle Finds the best-fitting great circle to a set of points gmt2rgb Convert Sun rasterfile or grid to r,g,bgrids gmtconvert Convert between ASCII and binary 1-D tables gmtdefaults List the current default settings gmtmath Mathematical operations on data tables gmtset Set individual default parameters gmtselect Extract data subsets based on spatial criteria grdfilter Filter 2-D data sets in the space domain grd2cpt Makeacolor palette table from grid files grd2xyz Conversion from 2-D grid file to table data grdblend Blend several partially over-lapping grid files onto one grid grdclip Limit the z-range in gridded data grdcontour Contouring of 2-D gridded data grdcut Cut a sub-region from a grid file grdedit Modify header information in a 2-D grid file grdfft Operate on grid files in the wav enumber (or frequency) domain grdgradient Compute directional gradient from grid files grdhisteq Histogram equalization for grid files grdimage Produce images from 2-D gridded data grdinfo Get information about grid files grdlandmask Create mask grid file from shoreline data base grdmask Reset nodes outside a clip path to a constant grdmath Mathematical operations on grid files grdpaste Paste together grid files along a common edge grdproject Project gridded data onto a newcoordinate system grdreformat Converting between different grid file formats grdsample Resample a 2-D gridded data set onto a newgrid grdtrend Fits polynomial trends to grid files grdtrack Sampling of 2-D data set along 1-D track GMT 4.5.18 1July 2018 1 GMT(1) Generic Mapping Tools GMT(1) grdvector Plot vector fields from grid files grdview 3-D perspective imaging of 2-D gridded data grdvolume Volume calculations from 2-D gridded data greenspline Interpolation using Green’sfunctions for splines in 1-3 dimensions makecpt Makecolor palette tables mapproject Forward or inverse map projections of table data minmax Find extreme values in data tables nearneighbor Nearest-neighbor gridding scheme project Project data onto lines/great circles ps2raster Crop and convert PostScript files to raster images, EPS, and PDF psbasemap Create a basemap plot psclip Use polygon files to define clipping paths pscoast Plot coastlines and filled continents on maps pscontour Contour xyz-data by triangulation pshistogram Plot a histogram psimage Plot images (EPS or Sun raster files) on maps pslegend Plot legend on maps psmask Create overlay to mask out regions on maps psrose Plot sector or rose diagrams psscale Plot gray scale or color scale on maps pstext Plot text strings on maps pswiggle Drawtime-series along track on maps psxy Plot symbols, polygons, and lines on maps psxyz Plot symbols, polygons, and lines in 3-D sample1d Resampling of 1-D table data sets spectrum1d Compute various spectral estimates from time-series splitxyz Split xyz-files into several segments surface Acontinuous curvature gridding algorithm trend1d Fits polynomial or Fourier trends to y = f(x) data trend2d Fits polynomial trends to z = f(x,y) data triangulate Perform optimal Delaunay triangulation and gridding xyz2grd Convert equidistant xyz data to a 2-D grid file SEE ALSO Look up the individual man pages for more details and full syntax. Run GMT without options to list all GMT programs and to showall installation directories. Information is also available on the GMT home page gmt.soest.hawaii.edu REFERENCES Wessel, P., and W.H.F.Smith, 2018, The Generic Mapping Tools (GMT) version 4.5.18 Technical Refer- ence & Cookbook, SOEST/NOAA. Wessel, P., and W.H.F.Smith, 1998, New, ImprovedVersion of Generic Mapping Tools Released, EOS Trans., AGU, 79 (47), p. 579. Wessel, P., and W.H.F.Smith, 1995, NewVersion of the Generic Mapping Tools Released, EOS Trans., AGU, 76 (33), p. 329. Wessel, P., and W.H.F.Smith, 1995, NewVersion of the Generic Mapping Tools Released, http://www.agu.org/eos_elec/95154e.html, Copyright 1995 by the American Geophysical Union. Wessel, P., and W.H.F.Smith, 1991, Free Software Helps Map and Display Data, EOS Trans., AGU, 72 (41), p. 441. GMT 4.5.18 1July 2018 2 BLOCKMEAN(1) Generic Mapping Tools BLOCKMEAN(1) NAME blockmean − filter to block average (x,y,z)data by L2 norm SYNOPSIS blockmean [ xyz[w]file(s) ] −Ixinc[unit][=|+][/yinc[unit][=|+]] −Rxmin/xmax/ymin/ymax[r][−C ][−E ][ −F ][−H[i][nrec]][−S[w|z]][−V ][−W[io]][−:[i|o]][−b[i|o][s|S|d|D[ncol]|c[var1/...]] ] [ −f[i|o]col- info ] DESCRIPTION blockmean reads arbitrarily located (x,y,z)triples [or optionally weighted quadruples (x,y,z,w)] from stan- dard input [or xyz[w]file(s)]and writes to standard output a mean position and value for every non-empty block in a grid region defined by the −R and −I arguments. Either blockmean, blockmedian,or block- mode should be used as a pre-processor before running surface to avoid aliasing short wav elengths. These routines are also generally useful for decimating or averaging (x,y,z)data. You can modify the precision of the output format by editing the D_FORMAT parameter in your .gmtdefaults4 file, or you may choose binary input and/or output using single or double precision storage. xyz[w]file(s) 3[or 4] column ASCII file(s) [or binary,see −b]holding (x,y,z[,w]) data values. [w]isanoptional weight for the data. If no file is specified, blockmean will read from standard input. −I x_inc [and optionally y_inc]isthe grid spacing. Optionally,append a suffix modifier. Geographi- cal (degrees) coordinates:Append m to indicate arc minutes or c to indicate arc seconds. If one of the units e, k, i,or n is appended instead, the increment is assumed to be giveninmeter,km, miles, or nautical miles, respectively,and will be converted to the equivalent degrees longitude at the middle latitude of the region (the conversion depends on ELLIPSOID). If /y_inc is givenbut set to 0 it will be reset equal to x_inc;otherwise it will be converted to degrees latitude. All coor- dinates:If = is appended then the corresponding max x (east)or y (north)may be slightly adjusted to fit exactly the givenincrement [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 + 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 Appendix B for details. Note: if −Rgrdfile is used then grid spacing has already been initialized; use −I to override the values. −R xmin, xmax, ymin,and ymax specify the Region of interest. Forgeographic regions, these limits correspond to west, east, south, and north and you may specify them in decimal degrees or in [+-]dd:mm[:ss.xxx][W|E|S|N] format. Append r if lower left and upper right map coordinates are giveninstead of w/e/s/n. The twoshorthands −Rg and −Rd stand for global domain (0/360 and -180/+180 in longitude respectively,with -90/+90 in latitude). Alternatively,specify the name of an existing grid file and the −R settings (and grid spacing, if applicable) are copied from the grid. Forcalendar time coordinates you may either give (a) relative time (relative tothe selected TIME_EPOCH and in the selected TIME_UNIT;append t to −JX|x), or (b) absolute time of the form [date]T[clock](append T to −JX|x). At least one of date and clock must be present; the T is always required. The date string must be of the form [-]yyyy[-mm[-dd]] (Gregorian calendar) or yyyy[-Www[-d]] (ISO week calendar), while the clock string must be of the form hh:mm:ss[.xxx]. The use of delimiters and their type and positions must be exactly as indicated (however, input, output and plot formats are customizable; see gmtdefaults). OPTIONS −C Use the center of the block as the output location [Default uses the mean location]. −E Provide Extended report which includes s (the standard deviation of the mean), l,the lowest value, and h,the high value for each block. Output order becomes x,y,z,s,l,h[,w]. [Default outputs x,y,z[,w]. See −W for w output.