NAME
pshistogram - Bin data and plot histograms
SYNOPSIS
pshistogram file -Jxxscale[/yscale] -Wbin_width [ -2 ] [ -A ] [ -Btick-
info ] [ -C ] [ -Eazimuth/elevation ] [ -Gfill ] [ -H[nrec] ] [ -I[o] ]
[ -K ] [ -Lpen ] [ -O ] [ -P ] [ -Q ] [ -Rwest/east/south/north[r] ] [
-S ] [ -U[/dx/dy/][label] ] [ -V ] [ -Xx-shift ] [ -Yy-shift ] [ -Ztype
] [ -ccopies ] [ -bi[s][n] ] [ -f[i|o]colinfo ]
DESCRIPTION
pshistogram reads the first column from file [or standard input] and
calculates histogram parameters based on the bin-width provided. Using
these parameters, scaling, and optional range parameters it will gener-
ate PostScript code that plots a histogram. A cumulative histogram may
also be specified.
file ASCII [or binary, see -b] datafile. If no file is given, pshis-
togram will read standard input.
-Jx xscale[/yscale] (Linear scale(s) in distance unit/data unit).
-W Sets the bin width used for histogram calculations.
OPTIONS
No space between the option flag and the associated arguments.
-2 Read second rather than first column.
-A Plot the histogram horizontally from x = 0 [Default is verti-
cally from y = 0].
-B Sets map boundary annotation and tickmark intervals; see the
psbasemap man page for all the details.
-C Center bin on each value. [Default is left edge].
-E Sets the viewpoint’s azimuth and elevation (for perspective
view) [180/90].
-G Select filling of bars. [Default is no fill]. Specify the shade
(0-255) or color (r/g/b, each in 0-255, h/s/v, with h in 0-360,
s, and v in 0-1, and set COLOR_MODEL = hsv, or c/m/y/k, each in
0-100 %). Alternatively, specify -Gpdpi/pattern, where pattern
gives the number of the built-in pattern (1-90) or the name of a
Sun 1-, 8-, or 24-bit raster file. The dpi sets the resolution
of the image. For 1-bit rasters: use -GP for inverse video, or
append :Fr/g/b[B[r/g/b]] to specify fore- and background colors
(use r/g/b = - for transparency). See GMT Cookbook & Technical
Reference Appendix E for information on individual patterns.
-H Input file(s) has Header record(s). Number of header records can
be changed by editing your .gmtdefaults4 file. If used, GMT
default is 1 header record. Use -Hi if only input data should
have header records [Default will write out header records if
the input data have them].
-I Inquire about min/max x and y after binning. No plotting is
done. Append o to output an ASCII table of the resulting x,y
data.
-K More PostScript code will be appended later [Default terminates
the plot system].
-L Draw bar outline using the specified pen thickness. [Default is
no outline].
-O Selects Overlay plot mode [Default initializes a new plot sys-
tem].
-P Selects Portrait plotting mode [GMT Default is Landscape, see
gmtdefaults to change this].
-Q Draw a cumulative histogram.
-R xmin, xmax, ymin, and ymax specify the Region of interest. For
geographic 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 given instead of wesn. The
two shorthands -Rg -Rd stand for global domain (0/360 or
-180/+180 in longitude respectively, with -90/+90 in latitude).
For calendar time coordinates you may either give relative time
(relative to the selected TIME_EPOCH and in the selected
TIME_UNIT; append t to -JX|x), or 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 as indicated (however, input/output
and plotting formats are flexible). If not given, pshistogram
will automatically find reasonable values for the region.
-S Draws a stairs-step diagram instead of histogram.
-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"].
-X -Y Shift origin of plot by (x-shift,y-shift). Prepend a for abso-
lute coordinates; the default (r) will reset plot origin. Give
c to center plot using current page size.
-Z Choose between 6 types of histograms: 0 = counts [Default], 1 =
frequency_percent, 2 = log (1.0 + count), 3 = log (1.0 + fre-
quency_percent), 4 = log10 (1.0 + count), 5 = log10 (1.0 + fre-
quency_percent).
-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 2 input columns].
-f Special formatting of input and output columns (time or geo-
graphical data) Specify i(nput) or o(utput) [Default is both
input and output]. Give one or more columns (or column ranges)
separated by commas. Append T (Absolute calendar time), t (time
relative to chosen TIME_EPOCH), x (longitude), y (latitude), g
(geographic coordinate), or f (floating point) to each column or
column range item.
EXAMPLES
To draw a histogram of the data v3206.t containing seafloor depths,
using a 250 meter bin width, center bars, and draw bar outline, use:
pshistogram v3206.t -JXh -W250 -C -L0.5p -V > plot.ps
If you know the distribution of your data, you may explicitly specify
range and scales. E.g., to plot a histogram of the y-values (2nd col-
umn) in the file errors.xy using a 1 meter bin width, plot from -10 to
+10 meters @ 0.75 cm/m, annotate every 2 m and 100 counts, and use
black bars, run:
cut -f2 errors.xy | pshistogram -W1 -R-10/10/0/0 -Jx0.75c/0.01c
-B2:Error:/100:Counts: -G0 -V > plot.ps
Since no y-range was specified, pshistogram will calculate ymax in even
increments of 100.
SEE ALSO
gmt(l), psbasemap(l), psrose(l), psxy(l)
GMT4.0 1 Oct 2004 PSHISTOGRAM(l)
Man(1) output converted with
man2html