NAME¶
psrose - Plot (length, azimuth) as windrose diagram or polar histogram (sector
or rose diagram).
SYNOPSIS¶
psrose file [
-Asector_width[
r] ] [
-Btickinfo ] [
-C[
mode-file] ] [
-Eazimuth/elevation ] [
-Gfill ] [
-H[
nrec] ] [
-I ] [
-K ] [
-M[
parameters ] [
-O ] [
-P ] [
-Rr0/r1/az_0/az_1 ] [
-Sradial_scale[
n] ] [
-T ] [
-U[
/dx/dy/][
label] ] [
-V ] [
-W pen ] [
-Xx-shift ] [
-Yy-shift ]
[
-Zscale ] [
-ccopies ] [
-: ] [
-bi[
s][
n] ]
DESCRIPTION¶
psrose reads (length,azimuth) pairs from
file [or standard input]
and generates
PostScript code that will plot a windrose diagram.
Optionally (with
-A), polar histograms may be drawn (sector diagram or
rose diagram). Options include full circle and half circle plots. The
PostScript code is written to standard output.
- file
- Name of ASCII [or binary, see -b] data file. If no
file is given, psrose will read standard input.
OPTIONS¶
No space between the option flag and the associated arguments.
- -A
- Gives the sector width in degrees for sector and rose
diagram. [Default 0 means windrose diagram]. Append r to draw rose
diagram instead of sector diagram.
- -B
- Sets map boundary tickmark intervals. See psbasemap
for details. Remember that "x" here is radial distance and
"y" is azimuth. The ylabel may be used to plot a figure
caption.
- -C
- Plot vectors showing the principal directions given in the
modes file. If no file is given, compute and plot mean
direction.
- -E
- Sets the viewpoint's azimuth and elevation [180/90]'
- -G
- Selects shade or color for sector infill [Default is no
fill]. Specify the shade (0-255) or color (r/g/b, each in 0-255).
- -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.
- -I
- Inquire. Computes statistics needed to specify useful
-R. No plot is generated.
- -K
- More PostScript code will be appended later [Default
terminates the plot system].
- -M
- Specify new arrow attributes
tailwidth/headlength/headwidth/r/g/b to change the appearance of arrows
(Only if -C is set). [Default is
0.075c/0.3c/0.25c/0/0/0 (or 0.03
i/0.12i/0.1 i/0/0/0)].
- -bo
- Selects binary output. Append s for single precision
[Default is double].
- -P
- Selects Portrait plotting mode [GMT Default is
Landscape, see gmtdefaults to change this].
- -R
- Specifies the 'region' of interest in (r,azimuth) space. r0
is 0, r1 is max length in units. For azimuth, specify -90/90 for half
circle plot or 0/360 for full circle.
- -S
- Specifies radius of circle. Append n to normalize
input radii to go from 0 to 1.
- -T
- Specifies that the input data is orientation data (has a
180 degree ambiguity) instead of true 0-360 degree directions
[Default].
- -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.). The GMT parameters UNIX_TIME and
UNIX_TIME_POS can affect the appearance; see the gmtdefaults man
page for details.
- -V
- Selects verbose mode, which will send progress reports to
stderr [Default runs "silently"].
- -W
- Set pen attributes for sector outline or rose plot.
[Default is no outline].
- -X -Y
- Shift origin of plot by (x-shift,y-shift). Prepend
a for absolute coordinates; the default ( r) will reset plot
origin.
- -Z
- Multiply the data radii by scale. E.g., use
-Z0.001 to convert your data from m to km [Default is no
scaling].
- -c
- Specifies the number of plot copies. [Default is 1]
- -:
- Input file has (azimuth,radius) pairs rather than the
expected (radius,azimuth).
- -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 2 input columns].
EXAMPLES¶
To plot a half circle rose diagram of the data in the file fault_segments.az_r
(containing pairs of (azimuth, length in meters), using a 10 degree bin sector
width, on a circle of radius = 3 inch, grid going out to radius = 150 km in
steps of 25 km with a 30 degree sector interval, radial direction annotated
every 50 km, using a light blue shading (200/240/255) outlined by a solid red
pen (width = 0.75 points), draw the mean azimuth, and shown in Portrait
orientation, try:
psrose fault_segments.az_r
-R0/150/-90/90
-B50g25:"Fault
length":/g30:."Rose diagram":
-S3
i
-A10
r -G200/240/255
-W0.75
p/255/0/0
-Z0.001
-C -P -T -: | lpr
To plot a full circle wind rose diagram of the data in the file lines.r_az, on a
circle of radius = 5 cm, grid going out to radius = 500 units in steps of 100
with a 45 degree sector interval, using a solid pen (width = 0.5 point), and
shown in landscape [Default] orientation with UNIX timestamp and command line
plotted, try:
psrose lines.az_r
-R0/500/0/360
-S5
c
-Bg100/
g45:."Windrose diagram":
-W0.5
p
-Uc | lpr
BUGS¶
No default radial scale and grid settings for polar histograms. User must run
psrose -I to find max length in binned data set.
SEE ALSO¶
gmt(1gmt),
gmtdefaults(1gmt),
pshistogram(1gmt)