NAME¶
pshistogram - Bin data and plot histograms
SYNOPSIS¶
pshistogram file
-Jxxscale[/yscale]
-Wbin_width [ -2 ] [
-A ] [ -Btickinfo ] [
-C ] [ -Eazimuth/elevation ] [
-Gfill ] [ -H[nrec] ]
[ -I[o] ] [ -K ] [
-L pen ] [ -O ] [ -P ]
[ -Q ] [
-Rwest/east/south/north[ r] ] [
-S ] [
-U[/dx/dy/][label] ] [
-V ] [ -Xx-shift ] [
-Y y-shift ] [ -Ztype
] [ -ccopies ] [
-bi[s ][n] ]
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 generate
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,
pshistogram 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.
- -B
- Sets map boundary tickmark intervals. See psbasemap for
details.
- -2
- Read second rather than first column.
- -A
- Plot the histogram horizontally from x = 0 [Default is vertically from y =
0].
- -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
- 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 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].
- -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].
- -Q
- Draw a cumulative histogram.
- -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. 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
absolute coordinates; the default ( r) will reset plot origin.
- -Z
- Choose between 6 types of histograms: 0 = counts [Default], 1 =
frequency_percent, 2 = log (1.0 + count), 3 = log (1.0 +
frequency_percent), 4 = log10 (1.0 + count), 5 = log10 (1.0 +
frequency_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].
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, try:
pshistogram v3206.t
-JXh
-W250
-C -L0.5
p -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 column) 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, try:
cut -f2 errors.xy | pshistogram
-W1
-R-10/10/0/0
-Jx0.75
c/0.01
c -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(1gmt),
psbasemap(1gmt),
psrose(1gmt),
psxy(1gmt)