NAME
psxy - Plot lines, polygons, and symbols on maps
SYNOPSIS
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] ]
DESCRIPTION
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. If no
symbols are specified then the symbol code (see -S below)
must be present as last column in the input. Multiple seg-
ment 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 standard
output.
files
List one or more file-names. If no files are given,
psxy will read standard input.
-J Selects the map projection. Scale is UNIT/degree,
1:xxxxx, or width in UNIT (upper case modifier). UNIT
is cm, inch, or m, depending on the MEASURE_UNIT set-
ting in .gmtdefaults, but this can be overridden on the
command line by appending the c, i, or m to the
scale/width value.
CYLINDRICAL PROJECTIONS:
-Jclon0/lat0/scale (Cassini)
-Jjlon0/scale (Miller)
-Jmscale (Mercator - Greenwich and Equator as origin)
-Jmlon0/lat0/scale (Mercator - Give meridian and stan-
dard parallel)
-Joalon0/lat0/azimuth/scale (Oblique Mercator - point
and azimuth)
-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, with Equator
as y = 0)
-Jtlon0/lat0/scale (TM - Transverse Mercator, set ori-
gin)
-Juzone/scale (UTM - Universal Transverse Mercator)
-Jylon0/lats/scale (Basic Cylindrical Projection)
AZIMUTHAL PROJECTIONS:
-Jalon0/lat0/scale (Lambert).
-Jelon0/lat0/scale (Equidistant).
-Jflon0/lat0/horizon/scale (Gnomonic).
-Jglon0/lat0/scale (Orthographic).
-Jslon0/lat0/scale (General Stereographic)
CONIC PROJECTIONS:
-Jblon0/lat0/lat1/lat2/scale (Albers)
-Jdlon0/lat0/lat1/lat2/scale (Equidistant)
-Jllon0/lat0/lat1/lat2/scale (Lambert)
MISCELLANEOUS PROJECTIONS:
-Jhlon0/scale (Hammer)
-Jilon0/scale (Sinusoidal)
-Jk[f|s]lon0/scale (Eckert IV (f) and VI (s))
-Jnlon0/scale (Robinson)
-Jrlon0/scale (Winkel Tripel)
-Jvlon0/scale (Van der Grinten)
-Jwlon0/scale (Mollweide)
NON-GEOGRAPHICAL PROJECTIONS:
-Jpscale[/origin] (polar (theta,r) coordinates, option-
ally offset theta [0])
-Jxx-scale[l|ppow][/y-scale[l|ppow]] (Linear, log, and
power scaling)
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.
OPTIONS
No space between the option flag and the associated argu-
ments.
-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. When used with -S, lets
symbol color be determined by the z-value in the third
column. Additional fields are shifted over by one
column (optional size would be 4th rather than 3rd
field, etc.). If -S is not set, psxy expects the user
to supply a multisegment polygon file (requires -M) and
will look for -Zval strings in each multisegment
header. The val will control the color via the cpt
file.
-E Draw error bars. Append x and/or y to indicate which
bars you want to draw (Default is both x and y). 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
[0.25c (or 0.1i)]. Pen attributes for error bars may
also be set. [Defaults: width = 1, color = 0/0/0, tex-
ture = 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 system].
-P Selects Portrait plotting mode [GMT Default is
Landscape, see gmtdefaults to change this].
-S Plot symbols. If present, size is symbol size in the
unit set in .gmtdefaults (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 cir-
cumscribed 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
MEASURE_UNIT].
-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 plot-
distance units]. 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, 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 projec-
tion (-Se leaves the directions unchanged.) Further-
more, the axes lengths must be given in km instead of
plot-distance units.
-Sf fault. Give distance gap between ticks and ticklength.
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 (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 ANOT_FONT].
-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 must be found in columns 3 and
4. size, if present, will be interpreted as
arrowwidth/headlength/headwidth [Default 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.
-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 projec-
tion (-Sv leaves the directions unchanged.)
-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.
Optionally, append a label, or c (which will plot the
command string.)
-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 pen.
-X -Y
Shift origin of plot by (x-shift,y-shift). 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 chosen settings].
EXAMPLES
To plot solid red circles (diameter = 0.25 cm) at the posi-
tions listed in the file DSDP.xy on a Mercator map at 5
cm/degree of the area 150E to 154E, 18N to 23N, with tick-
marks every 1 degree and gridlines every 15 minutes, try:
psxy DSDP.xy -R150/154/18/23 -Jm5c -Sc0.25c -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.25 inch on the left side of the line,
spaced every 0.8 inch, use
psxy trench.xy -R150/200/20/50 -Jm0.15i -Sf0.8i/0.1iL -G255
-W -B10 | lpr br
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 mag-
nitude 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
BUGS
The -N option does not adjust the BoundingBox information so
you may have to post-process the PostScript outout with
epstool or ps2epsi to obtain a correct BoundingBox.
SEE ALSO
gmt(l), psbasemap(l), psxyz(l)