psxy - Plot lines, polygons, and symbols on maps
psxy files -Jparameters -Rwest/east/south/north[r] [ -A ] [ -Btickinfo
] [ -Ccptfile ] [ -E[x][y][cap][/pen] ] [ -Gfill ] [ -H[nrec] ] [ -K ]
[ -L ] [ -N ] [ -M[flag] ] [ -O ] [ -P ] [ -S[symbol][size] ] [
-U[/dx/dy/][label] ] [ -V ] [ -W[pen] ] [ -Xx-shift ] [ -Yy-shift ] [
-: ] [ -ccopies ] [ -bi[s][n] ]
psxy reads (x,y) pairs from files [or standard input] and generates
PostScript code that will plot lines, polygons, or symbols at those
locations on a map. If a symbol is selected and no symbol size given,
then psxy will interpret the third column of the input data as symbol
size in inch. If no symbols are specified then the symbol code (see
-S below) must be present as last column in the input. Multiple
segment files may be plotted using the -M option. If -S is not
selected, a line connecting the data points will be drawn instead. To
explicitly close polygons, use -L. Select a shade with -G. If -G is
set, -W will control whether the polygon outline is drawn or not. If
a symbol is selected, -G and -W determines the fill color and
outline/no outline, respectively. The PostScript code is written to
List one or more file-names. If no files are given, psxy will
read standard input.
-J Selects the map projection. Scale is inch/degree, 1:xxxxx, or
width in inch (upper case modifier).
-Jmscale (Mercator - Greenwich and Equator as origin)
-Jmlon0/lat0/scale (Mercator - Give meridian and standard
-Joalon0/lat0/azimuth/scale (Oblique Mercator - point and
-Joblon0/lat0/lon1/lat1/scale (Oblique Mercator - two points)
-Joclon0/lat0/lonp/latp/scale (Oblique Mercator - point and pole)
-Jqlon0/scale (Equidistant Cylindrical Projection (Plate Carree))
-Jtlon0/scale (TM - Transverse Mercator)
-Juzone/scale (UTM - Universal Transverse Mercator)
-Jylon0/lats/scale (Basic Cylindrical Projection)
-Jslon0/lat0/scale (General Stereographic)
-Jklon0/scale (Eckert VI)
-Jrlon0/scale (Winkel Tripel)
-Jpscale (Linear projection for polar (theta,r) coordinates)
-Jxx-scale[l|ppow][/y-scale[l|ppow]] (Linear, log, and power
More details can be found in the psbasemap manpages.
-R west, east, south, and north specify the Region of interest. To
specify boundaries in degrees and minutes [and seconds], use the
dd:mm[:ss] format. Append r if lower left and upper right map
coordinates are given instead of wesn.
No space between the option flag and the associated arguments.
-A Suppress drawing line segments as great circle Arcs. [Default
draws great circle arcs.]
-B Sets map boundary tickmark intervals. See psbasemap for details.
-C Give a color palette file and let symbol color be determined by
the z-value in the third column. Must be used with -S. If used,
then additional fields are shifted over by one column (optional
size would be 4th rather than 3rd field, etc.).
-E Draw error bars. Append x and/or y to indicate which bars you
want to draw. The x and/or y errors must be stored in the
columns after the (x,y) pair [or (x,y,size) triplet]. The cap
parameter indicates the length of the end-cap on the error bars
in inch [0.1i]. Pen attributes for error bars may also be set.
[Defaults: width = 1, color = 0/0/0, texture = solid].
-G Select filling of polygons and symbols. Append the shade
(0-255), color (r/g/b), or p|pdpi/pattern (polygons only)
[Default is no fill]. Note when -M is chosen, psxy will search
for -G and -W strings in all the subheaders and let any found
values over-ride the command line settings.
-H Input file(s) has Header record(s). Number of header records can
be changed by editing your .gmtdefaults file. If used, GMT
default is 1 header record.
-K More PostScript code will be appended later [Default terminates
the plot system].
-L Force closed polygons: connect the endpoints of the line-
segment(s) and draw polygons.
-M Multiple segment file. Segments are separated by a record whose
first character is flag. [Default is '>'].
-N Do NOT skip symbols that fall outside map border [Default plots
points inside border only].
-O Selects Overlay plot mode [Default initializes a new plot
-P Selects Portrait plotting mode [GMT Default is Landscape, see
gmtdefaults to change this].
-S Plot symbols. If present, size is symbol size in inch (unless c,
i, m, or p is appended). The uppercase symbols A, C, D, H, I, S,
T are normalized to have the same area as the circle, while the
corresponding lowercase symbols all are circumscribed by the
circle. Choose between these symbol codes:
-S Read symbol code (see below) from last column in the input data.
Cannot be used in conjunction with -b. Optionally, append c, i,
m, p to indicate that the size information in the input data is
in units of cm, inch, meter, or point, respectively. [Default is
-Sa star. size is diameter of circumscribing circle.
-Sb bar extending from base to y. size is bar width. Append u if
size is in x-units [Default is inch]. By default, base = 0.
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 counterclockwise from
horizontal), major_axis (in inch), and minor_axis (in inch) 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 inch.
-Sf fault. Give distance gap between ticks and ticklength in inch.
If gap is negative, it is interpreted to mean number of ticks
instead. Append l or r to draw tick on the left or right side of
line [Default is centered]. Upper case L or R draws a triangle
instead of line segment.
-Sh hexagon. size is diameter of circumscribing circle.
-Si inverted triangle. size is diameter of circumscribing circle.
-Sl letter or text string. Give size in inch, 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 *.
-Sp point. No size needs to be specified (1 pixel is used).
-Ss square. size is diameter of circumscribing circle.
-St triangle. size is diameter of circumscribing circle.
-Sv vector. Direction (in degrees counterclockwise from horizontal)
and length (in inch) must be found in columns 3 and 4. size, if
present, will be interpreted as arrowwidth/headlength/headwidth
(in inch) [Default is 0.03i/0.12i/0.1i inch]. 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.
-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
-Sx cross. size is diameter of circumscribing circle.
-U Draw Unix System time stamp on plot. User may specify where the
lower left corner of the stamp should fall on the page relative
to lower left corner of plot in inch [Default is (-0.75i/-
0.75i)]. Optionally, append a label, or c (which will plot the
-V Selects verbose mode, which will send progress reports to stderr
[Default runs "silently"].
-W Set pen attributes. [Defaults: width = 1, color = 0/0/0, texture
= solid]. Implicitly draws the outline of symbols with selected
Shift origin of plot by (x-shift,y-shift) inch [Default is
(a1i,a1i) for new plots, (0,0) for overlays]. Prepend a for
absolute coordinates; the default (r) will reset plot origin.
-: Toggles between (longitude,latitude) and (latitude,longitude)
input/output. [Default is (longitude,latitude)].
-c Specifies the number of plot copies. [Default is 1]
-bi Selects binary input. Append s for single precision [Default is
double]. Append n for the number of columns in the binary
file(s). [Default is the required number of columns given the
To plot solid red circles (diameter = 0.1 inch) at the positions
listed in the file DSDP.xy on a Mercator map at 2 inch/degree of the
area 150E to 154E, 18N to 23N, with tickmarks every 1 degree and
gridlines every 15 minutes, try:
psxy DSDP.xy -R150/154/18/23 -Jm2i -Sc0.1i -G255/0/0 -B1g15m | lpr
To plot the xyz values in the file quakes.xyzm as circles with size
given by the magnitude in the 4th column and color based on the depth
in the third using the color palette cpt on a linear map, try
psxy quakes.xyzm -R0/1000/0/1000 -JX6i -Sc -Ccpt -B200 > map.ps
To plot the file trench.xy on a Mercator map, with white triangles
with sides 0.1 inch on the left side of the line, spaced every 0.8
psxy trench.xy -R150/200/20/50 -Jm0.15i -Sf0.8i/0.1iL -G255 -W -B10 |
To plot the data in the file misc.d as symbols determined by the code
in the last column, and with size given by the magnitude in the 4th
column, and color based on the third column via the color palette cpt
on a linear map, try
psxy misc.d -R0/100/-50/100 -JX6i -S -Ccpt -B20 > t.ps
gmt, psbasemap, psxyz
Man(1) output converted with