WANtaro雑記帳（Win_GMT：基本的なことなど[1]）

目次
GMT program summary (gmt_man.html より抜粋)
GMT Version 4.5.7 Technical Reference and Cookbook からの抜粋
>gmtdefaults -D による表示内容
プロット記号 (man の psxy より抜粋)

GMT program summary (gmt_man.html より抜粋)

No	Program	Function
1	blockmean	L2 (x,y,z) data filter/decimator
2	blockmedian	L1 (x,y,z) data filter/decimator
3	blockmode	Mode-estimating (x,y,z) data filter/decimator
4	filter1d	Filter 1-D data (time series)
5	fitcircle	Finds best-fitting great or small circles
6	gmt_shell_functions	Simplifying GMT bash scripting
7	gmt2rgb	Convert Sun raster or grdfile to red, green, blue component grids
8	gmtcolors	Information on how to specify colors in GMT
9	gmtconvert	Convert table data from one format to another
10	gmtdefaults	List the current default settings
11	gmtget	Get present value of a GMT defaults parameter
12	gmtlogo	Plot a GMT logo as an overlay
13	gmtmath	Reverse Polish calculator for table data
14	gmtselect	Select table subsets based on multiple spatial criteria
15	gmtset	Edit parameters in the .gmtdefaults file
16	grd2cpt	Make color palette table from grdfile
17	grd2xyz	Convert 2-D gridded data to table
18	grdblend	Blend several gridded data sets into one
19	grdclip	Limit the z-range in gridded data sets
20	grdcontour	Contouring of 2-D gridded data
21	grdcut	Cut a sub-region from a grd file
22	grdedit	Modify grd header information
23	grdfft	Operate on grdfiles in frequency domain
24	grdfilter	Filter 2-D data in space domain
25	grdgradient	Compute directional gradient from grdfiles
26	grdhisteq	Histogram equalization for grdfiles
27	grdimage	Produce images from 2-D gridded datar
28	grdinfo	Get information about grd files
29	grdlandmask	Creates mask grdfile from coastline database
30	grdmask	Set nodes outside a clip path to a constant
31	grdmath	Reverse Polish calculator for grdfiles
32	grdpaste	Paste together grdfiles along common edge
33	grdproject	Project gridded data onto new coordinate system
34	grdreformat	Convert from one grdformat to another
35	grdsample	Resample a 2-D gridded data onto new grid
36	grdtrack	Sampling of 2-D data along 1-D track
37	grdtrend	Fits polynomial trends to grdfiles (z = f(x,y))
38	grdvector	Plot vector fields from 2-D gridded data
39	grdview	3-D perspective imaging of 2-D gridded data
40	grdvolume	Calculating volume under a surface within a contour
41	greenspline	Gridding using Green's function splines
42	isogmt	Run GMT command or script in isolation mode
43	makecpt	Create GMT color palette tables
44	mapproject	Transformation of coordinate systems
45	minmax	Report extreme values in table datafiles
46	nearneighbor	Nearest-neighbor gridding scheme
47	project	Project data onto lines/great circles
48	ps2raster	Crop and convert PostScript to raster image, EPS and PDF
49	psbasemap	Create a basemap frame
50	psbbox	Replace BoundingBox line in PostScript files
51	psclip	Use polygon files as clipping paths
52	pscoast	Plot coastlines, filled continents, rivers, and political borders
53	pscontour	Direct contouring or imaging of xyz-data by triangulation
54	pshistogram	Plot a histogram
55	psimage	Plot Sun rasterfiles on a map
56	pslegend	Plot legend on a map
57	psmask	Create overlay to mask specified regions of a map
58	psrose	Plot sector or rose diagrams
59	psscale	Plot grayscale or colorscale
60	pstext	Plot textstrings
61	pswiggle	Draw anomalies along track
62	psxy	Plot symbols, polygons, and lines in 2-D
63	psxyz	Plot symbols, polygons, and lines in 3-D
64	sample1d	Resampling of 1-D data
65	spectrum1d	Compute spectral estimates from time-series
66	splitxyz	Split xyz files into several segments
67	surface	Continuous curvature gridding algorithm
68	trend1d	Fits polynomial or Fourier trends to y = f(x) series
69	trend2d	Fits polynomial trends to z = f(x,y) series
70	triangulate	Perform optimal Delauney triangulation on xyz data
71	triangulate	Perform optimal Delauney triangulation on xyz data
72	xyz2grd	Convert table to 2-D grd file

GMT Version 4.5.7 Technical Reference and Cookbook からの抜粋

Table 4.1: The 17 standardized GMT command line switches.

Option	Meaning
-B	Defines tickmarks, annotations, and labels for basemaps and axes
-H	Specifies that input/output tables have header record(s)
-J	Selects a map projection or coordinate transformation
-K	Allows more plot code to be appended to this plot later
-O	Allows this plot code to be appended to an existing plot
-P	Selects Portrait plot orientation [Default is landscape]
-R	Defines the extent of the map/plot region
-U	Plots a time-stamp, by default in the lower left corner of page
-V	Selects verbose operation; reporting on progress
-X	Sets the x-coordinate for the plot origin on the page
-Y	Sets the y-coordinate for the plot origin on the page
-b	Selects binary input and/or output
-c	Specifies the number of plot copies
-f	Specifies the data format on a per column basis
-g	Identify data gaps based on supplied criteria
-m	Specifies data in multiple segment format
-:	Assumes input geographic data are (lat,lon) and not (lon,lat)

Table 4.2: Interval type codes.

Flag	Description
a	Annotation tick spacing
f	Frame tick spacing
g	Grid tick spacing

Table 4.7: GMT text escape sequences.

Code	Effect
@~	Turns symbol font on or off
@+	Turns superscript on or off
@-	Turns subscript on or off
@#	Turns small caps on or off
@_	Turns underline on or off
@%fontno%	Switches to another font; @%% resets to previous font
@:size:	Switches to another font size; @:: resets to previous size
@;color;	Switches to another font color; @;; resets to previous color
@!	Creates one composite character of the next two characters
@@	Prints the @ sign itself

>gmtdefaults -D による表示内容

#
#	GMT-SYSTEM 4.5.7 Defaults file
#
#-------- Plot Media Parameters -------------
PAGE_COLOR		= white
PAGE_ORIENTATION	= landscape
PAPER_MEDIA		= a4
#-------- Basemap Annotation Parameters ------
ANNOT_MIN_ANGLE		= 20
ANNOT_MIN_SPACING	= 0
ANNOT_FONT_PRIMARY	= Helvetica
ANNOT_FONT_SIZE_PRIMARY	= 14p
ANNOT_OFFSET_PRIMARY	= 0.2c
ANNOT_FONT_SECONDARY	= Helvetica
ANNOT_FONT_SIZE_SECONDARY	= 16p
ANNOT_OFFSET_SECONDARY	= 0.2c
DEGREE_SYMBOL		= ring
HEADER_FONT		= Helvetica
HEADER_FONT_SIZE	= 36p
HEADER_OFFSET		= 0.5c
LABEL_FONT		= Helvetica
LABEL_FONT_SIZE		= 24p
LABEL_OFFSET		= 0.3c
OBLIQUE_ANNOTATION	= 1
PLOT_CLOCK_FORMAT	= hh:mm:ss
PLOT_DATE_FORMAT	= yyyy-mm-dd
PLOT_DEGREE_FORMAT	= ddd:mm:ss
Y_AXIS_TYPE		= hor_text
#-------- Basemap Layout Parameters ---------
BASEMAP_AXES		= WESN
BASEMAP_FRAME_RGB	= black
BASEMAP_TYPE		= fancy
FRAME_PEN		= 1.25p
FRAME_WIDTH		= 0.2c
GRID_CROSS_SIZE_PRIMARY	= 0c
GRID_PEN_PRIMARY	= 0.25p
GRID_CROSS_SIZE_SECONDARY	= 0c
GRID_PEN_SECONDARY	= 0.5p
MAP_SCALE_HEIGHT	= 0.2c
POLAR_CAP		= 85/90
TICK_LENGTH		= 0.2c
TICK_PEN		= 0.5p
X_AXIS_LENGTH		= 25c
Y_AXIS_LENGTH		= 15c
X_ORIGIN		= 2.5c
Y_ORIGIN		= 2.5c
UNIX_TIME		= FALSE
UNIX_TIME_POS		= BL/-2c/-2c
UNIX_TIME_FORMAT	= %Y %b %d %H:%M:%S
#-------- Color System Parameters -----------
COLOR_BACKGROUND	= black
COLOR_FOREGROUND	= white
COLOR_NAN		= 128
COLOR_IMAGE		= adobe
COLOR_MODEL		= rgb
HSV_MIN_SATURATION	= 1
HSV_MAX_SATURATION	= 0.1
HSV_MIN_VALUE		= 0.3
HSV_MAX_VALUE		= 1
#-------- PostScript Parameters -------------
CHAR_ENCODING		= ISOLatin1+
DOTS_PR_INCH		= 300
GLOBAL_X_SCALE		= 1
GLOBAL_Y_SCALE		= 1
N_COPIES		= 1
PS_COLOR		= rgb
PS_IMAGE_COMPRESS	= lzw
PS_IMAGE_FORMAT		= ascii
PS_LINE_CAP		= butt
PS_LINE_JOIN		= miter
PS_MITER_LIMIT		= 35
PS_VERBOSE		= FALSE
TRANSPARENCY		= 0
#-------- I/O Format Parameters -------------
D_FORMAT		= %.12lg
FIELD_DELIMITER		= tab
GRIDFILE_FORMAT		= nf
GRIDFILE_SHORTHAND	= FALSE
INPUT_CLOCK_FORMAT	= hh:mm:ss
INPUT_DATE_FORMAT	= yyyy-mm-dd
IO_HEADER		= FALSE
N_HEADER_RECS		= 1
NAN_RECORDS		= pass
OUTPUT_CLOCK_FORMAT	= hh:mm:ss
OUTPUT_DATE_FORMAT	= yyyy-mm-dd
OUTPUT_DEGREE_FORMAT	= D
XY_TOGGLE		= FALSE
#-------- Projection Parameters -------------
ELLIPSOID		= WGS-84
MAP_SCALE_FACTOR	= default
MEASURE_UNIT		= cm
#-------- Calendar/Time Parameters ----------
TIME_FORMAT_PRIMARY	= full
TIME_FORMAT_SECONDARY	= full
TIME_EPOCH		= 2000-01-01T12:00:00
TIME_IS_INTERVAL	= OFF
TIME_INTERVAL_FRACTION	= 0.5
TIME_LANGUAGE		= us
TIME_UNIT		= d
TIME_WEEK_START		= Sunday
Y2K_OFFSET_YEAR		= 1950
#-------- Miscellaneous Parameters ----------
HISTORY			= TRUE
INTERPOLANT		= akima
LINE_STEP		= 0.025c
VECTOR_SHAPE		= 0
VERBOSE			= FALSE

プロット記号( man の psxy より抜粋)

注意：Sq は省略してあります．

The uppercase symbols A, C, D, G, H, I, N, S, T are normalized to have the same area as a circle with diameter size, while the size of the corresponding lowercase symbols refers to the diameter of a circumscribed circle. Choose between these symbol codes:

−S-	x-dash (-). size is the length of a short horizontal line segment.
−S+	plus (+). size is diameter of circumscribing circle.
−Sa	star. size is diameter of circumscribing circle.
−Sb	Vertical bar extending from base to y. size is bar width. Append u if size is in x-units [Default is plot-distance units]. By default, base = ymin. Append bbase to change this value.
−SB	Horizontal bar extending from base to x. size is bar width. Append u if size is in y-units [Default is plot-distance units]. By default, base = xmin. Append bbase to change this value.
−Sc	circle. size is diameter of circle.
−Sd	diamond. size is diameter of circumscribing circle.
−Se	ellipse. Direction (in degrees counter-clockwise from horizontal), major_axis, and minor_axis must be found in columns 3, 4, and 5.
−SE	Same as −Se, except azimuth (in degrees east of north) should be given instead of direction. The azimuth will be mapped into an angle based on the chosen map projection (−Se leaves the directions unchanged.) Furthermore, the axes lengths must be given in km instead of plot-distance units. An exception occurs for a linear projection in which we assume the ellipse axes are given in the same units as −R.
−Sf	front. −Sfgap/size[dir][type][:offset]. Supply distance gap between symbols and symbol size. If gap is negative, it is interpreted to mean the number of symbols along the front instead. Append dir to plot symbols on the left or right side of the front [Default is centered]. Append type to specify which symbol to plot: box, circle, fault, slip, or triangle. [Default is fault]. Slip means left-lateral or right-lateral strike-slip arrows (centered is not an option). Append :offset to offset the first symbol from the beginning of the front by that amount [Default is 0].
−Sg	octagon. size is diameter of circumscribing circle.
−Sh	hexagon. size is diameter of circumscribing circle.
−Si	inverted triangle. size is diameter of circumscribing circle.
−Sj	Rotated rectangle. Direction (in degrees counter-clockwise from horizontal), x-dimension, and y-dimension must be found in columns 3, 4, and 5.
−SJ	Same as −Sj, except azimuth (in degrees east of north) should be given instead of direction. The azimuth will be mapped into an angle based on the chosen map projection (−Sj leaves the directions unchanged.) Furthermore, the dimensions must be given in km instead of plot-distance units. An exception occurs for a linear projection in which we assume the dimensions are given in the same units as −R.
−Sk	kustom symbol. Append /size, and we will look for a definition file called .def in (1) the current directory or (2) in ~/.gmt or (3) in $GMT_SHAREDIR/custom. The symbol as defined in that file is of size 1.0 by default; the appended size will scale symbol accordingly. Users may add their own custom *.def files; see CUSTOM SYMBOLS below.
−Sl	letter or text string (less than 64 characters). Give size, and append /string after the size. Note that the size is only approximate; no individual scaling is done for different characters. Remember to escape special characters like *. Optionally, you may append %font to select a particular font [Default is ANNOT_FONT_PRIMARY].
−Sm	math angle arc, optionally with one or two arrow heads. The size is the radius of the arc. Start and stop directions (in degrees counter-clockwise from horizontal) for arc must be found in columns 3 and 4. Use −Smf to add arrow head at first angle, −Sml for arrow head at last angle, and −Smb for both [Default is no arrow heads].
−Sn	pentagon. size is diameter of circumscribing circle.
−Sp	point. No size needs to be specified (1 pixel is used).
−Sr	rectangle. No size needs to be specified, but the x- and y-dimensions must be found in columns 3 and 4.
−Ss	square. size is diameter of circumscribing circle.
−St	triangle. size is diameter of circumscribing circle.
−Sv	vector. Direction (in degrees counter-clockwise from horizontal) and length must be found in columns 3 and 4. size, if present, will be interpreted as arrowwidth/headlength/headwidth [Default unit is 0.075c/0.3c/0.25c (or 0.03i/0.12i/0.1i)]. By default arrow attributes remains invariant to the length of the arrow. To have the size of the vector scale down with decreasing size, append nnorm, where vectors shorter than norm will have their attributes scaled by length/norm. To center vector on balance point, use −Svb; to align point with the vector head, use −Svh; to align point with the vector tail, use −Svt [Default]. To give the head point’s coordinates instead of direction and length, use −Svs. Upper case B, H, T, S will draw a double-headed vector [Default is single head].
−SV	Same as −Sv, except azimuth (in degrees east of north) should be given instead of direction. The azimuth will be mapped into an angle based on the chosen map projection (−Sv leaves the directions unchanged.)
−Sw	pie wedge. Start and stop directions (in degrees counter-clockwise from horizontal) for pie slice must be found in columns 3 and 4.
−SW	Same as −Sw, except azimuths (in degrees east of north) should be given instead of the two directions. The azimuths will be mapped into angles based on the chosen map projection (−Sw leaves the directions unchanged.)
−Sx	cross (x). size is diameter of circumscribing circle.
−Sy	y-dash (\|). size is the length of a short vertical line segment.