MB-System Unix Manual Page
mbm_plot
Section: MB-System 5.0 (1)
Updated: 3 June 2013
Index
NAME
mbm_plot - Create an executable shellscript which will generate
a GMT plot of swath sonar swath data.
VERSION
Version 5.0
SYNOPSIS
mbm_plot -Fformat -Ifile
[-A[magnitude[/azimuth]
-C[cont_int/col_int/tic_int/lab_int/tic_len/lab_hgt]
-D[flipcolor/flipshade]
-Gcolor_mode[F] -H
-N[ttick/tannot/dannot/tlen[/nhgt] | F | FP]
-Oroot -Ppagesize
-S[color/shade] -T
-Uorientation -V
-W[color_style[/palette[ncolors]] | cptfile] ]
Additional Options:
[-Btickinfo
-Jprojection[/scale | width]
-Ltitle[:scale_label] -Mmisc
-Q -Rw/e/s/n -X -Y
-Zmin/max ]
Miscellaneous Options:
[-MGDgmtdef/value
-MGFscale_loc
-MGL[f][x]lon0/lat0/slat/length[m]
-MGQdpi
-MGTx/y/size/angle/font/just/text
-MGU[/dx/dy/][label]
-MMAfactor/mode/depth
-MMByr/mo/da/hr/mn/sc
-MMDmode/scale[/min/max]
-MMEyr/mo/da/hr/mn/sc
-MMLlonflip
-MMNnplot -MMPpings
-MMSspeedmin
-MMTtimegap -MMZalgorithm
-MNA[nhgt[/P] | P]
-MNP[pingnumber_tick/pingnumber_annot/pingnumber_tick_len]
-MTCfill -MTDresolution
-MTGfill -MTIriver[/pen]
-MTNborder[/pen] -MTSfill
-MTWpen
-MXGfill -MXIxy_file
-MXM
-MXSsymbol/size -MXWpen]
DESCRIPTION
mbm_plot is a macro to generate a shellscript of MB-System and GMT commands
which, when executed, will generate a Postscript plot of the
specified swath sonar data. The plot may include bathymetry color
fill (-G1), bathymetry color shaded relief (-G2), bathymetry
shaded with amplitudes (-G3), greyshade fill amplitude (-G4),
greyshade fill sidescan (-G5), contoured bathymetry (-C),
or annotated navigation. The plot may also include navigation tracks, text
labels, xy data in lines or symbols, and coastlines.
Five different color
schemes are included. The plot will be scaled to fit on
the specified page size or, if the scale is user defined,
the page size will be chosen in accordance with the plot
size. The primary purpose of this macro is to allow the
simple, semi-automated production of nice looking maps with
a few command line arguments. For users seeking more control
over the plot appearance, a number of additional optional
arguments are provided. Truly ambitious users may edit the
plot shellscript to take advantage of MB-System and GMT
capabilities not supported by this macro.
AUTHORSHIP
David W. Caress (caress@mbari.org)
Monterey Bay Aquarium Research Institute
Dale N. Chayes (dale@ldeo.columbia.edu)
Lamont-Doherty Earth Observatory
Suzanne H. O'Hara (sohara@ldeo.columbia.edu)
Lamont-Doherty Earth Observatory
SIMPLE DESCRIPTION OF BASIC OPTIONS
- -A
-
magnitude/azimuth or magnitude/median
Sets the parameters which control how mbswath generates
simulated illumination of bathymetry, which can be either
shaded relief bathymetry or bathymetry draped with amplitude data.
If mode is set to 2 (shaded relief bathymetry) using the
-G option, then the value magnitude
is an effective vertical exageration which modulates the intensity of
the shading; typical values are in the 1-5 range. The value azimuth
is the azimuth from which the bathymetry is illuminated.
If mode is set to 3 (bathymetry shaded using amplitudes) using the
-G option, then the value magnitude
modulates the intensity of the shading; the value median sets the
amplitude value which serves as the zero or neutral level.
- -C
-
Given by itself, the -C option produces a contour plot
of the bathymetry data with a contour interval chosen according
to the data in the file or files. If the -G option is used,
the default contours will be drawn in black; otherwise, the
default contours will be drawn in four colors (black, red,
green, and blue) with color changes, annotations, and tickmarks
every fourth contour interval. Additional optional parameters are
described in the COMPLETE DESCRIPTION OF OPTIONS section below.
- -D
-
[flipcolor/flipshade]
This option flips the color and shading conventions used by
mbm_plot.
Normally, the color or grayscale tables used for color bathymetry
maps run from
cool colors (or dark grays) for large depth values
to hot colors (or light grays) for small depth values.
In contrast, sidescan and beam amplitude data is normally plotted
using light grays (or hot colors) for small amplitudes and
dark grays (or cool colors) for large amplitudes.
If -D is given alone or with flipcolor = 1,
it applies to the color table used
for color or gray fill plots, shaded or unshaded. If the plot
is to be shaded, either by synthetic illumination (-G2)
or by overlaying amplitude data (-G3 option), then
setting flipshade = 1 will cause the shading convention
to be reversed (e.g. high amplitudes overlaid as light shading).
Using -D0/1 will flip the shading convention
but leave the default color convention.
- -F
-
format
Sets the data format for the input data.
If format < 0, then the input file specified
with the -I option will actually contain a list of input swath sonar
data files. This program uses the MBIO library
and will read or write any swath sonar
format supported by MBIO. A list of the swath sonar data formats
currently supported by MBIO and their identifier values
is given in the MBIO manual page. Default: format = -1.
- -G
-
mode[F]
Turns on color fill swath plot and sets the style of the plot.
mode = 1: Color fill of bathymetry data.
mode = 2: Color shaded relief bathymetry.
mode = 3: Bathymetry shaded using amplitude data.
mode = 4: Grayscale fill of amplitude data.
mode = 5: Grayscale fill of sidescan data.
If "F" is appended to mode, then mbm_plot will attempt
to plot amplitude or sidescan data that have been filtered with mbfilter.
If the desired filtered data files do not exist, then the plotting script
generated by mbm_plot will fail when
mbswath exits with an error message. Filtered amplitude
data are stored in ancilliary files ending with ".ffa", and filtered
sidescan files end in ".ffs". Filtering of bathymetry data is not supported,
and so appending "F" to mode values of 1 or 2 will have no effect.
- -H
-
This "help" flag cause the program to print out a description
of its operation and then exit immediately.
- -I
-
filename
Sets the input filename. If format > 0 (set with the
-f option) then the swath sonar data contained in infile
is read and processed. If format < 0 (the default),
then infile
is assumed to be an ascii file containing a list of the input swath sonar
data files to be processed and their formats. The program will read
the data in each one of these files.
In the infile file, each
data file should be followed by a data format identifier, e.g.:
datafile1 11
datafile2 24
This program uses the MBIO library and will read or write any swath sonar
format supported by MBIO. A list of the swath sonar data formats
currently supported by MBIO and their identifier values
is given in the MBIO manual page. The default input filename is
"datalist.mb-1".
- -N
-
Given by itself, this option causes a navigation track plot to be generated.
Additional optional parameters, including annotation control, are
described in the COMPLETE DESCRIPTION OF OPTIONS section below.
- -O
-
root
Sets the root used to construct the filename of the output shellscript
(root.cmd) and names of files created when the shellscript is
run. By default, the
name of the input data file or list file is used as the root.
- -P
-
pagesize
This option sets the size of the page the plot will be centered
on. If the user does not set the plot scale, the plot will be
sized as large as will fit on the designated page. If the user
sets the plot scale such that the plot will not fit on the
designated page, a larger page will be used.
The supported page sizes include ANSI A, B, C, D, E,
F, and E1, as well as most metric page sizes. See the
COMPLETE DESCRIPTION OF OPTIONS section
below for a complete list of
the supported page sizes. The default page size is A.
- -S
-
[color/shade]
This option enables effective histogram equalization of the
color and/or shading of thedata. The equalization is
not achieved by changing the data values, but rather by
constructing the color or shading tables so that
the boundaries in the tables encompass equal fractions of
the datapoints. This serves to focus color or shading contrasts
in value ranges corresponding to the bulk of the data values,
and is particularly useful for
enhancing the contrast of sidescan and beam amplitude plots.
If -S is given alone or with color = 1,
it enables equalization of the color table used
for color or gray fill plots, shaded or unshaded. If the plot
is to be shaded, either by synthetic illumination (-G2)
or by overlaying amplitude data (-G3 option), then
setting shade = 1 will cause the shading to be equalized.
Using -S0/1 will equalize the shading without
equalizing the color table.
- -T
-
If -T is given, it causes a coastline to be drawn
on the map. The default is to draw the coastline and shade
all dry land a uniform gray. To exercise greater control of
the coastline plotting, use the
-MTC, -MTD, -MTG, -MTI,
-MTN, -MTS, and -MTW options described in the
COMPLETE DESCRIPTION OF OPTIONS section below.
- -U
-
orientation
Normally the orientation of the plot (portrait or landscape)
is selected automatically so as to maximize the plot scale.
The -U option allows the user to set the plot orientation. If
orientation = 1, a portrait plot will be produced; if
orientation = 2, a landscape plot will be produced.
- -V
-
Causes mbm_plot to operate in "verbose" mode
so that it outputs
more information than usual.
- -W
-
[color_style[/palette[ncolors]] | cptfile]
This option controls the color scheme used for color
fill plots.
If color_style = 1 [default], then
the color scheme used will be a continuous grading
of colors. If color_style = 2, the color scheme
will be a set of discrete color intervals. The color
palette used is set using palette. Five palettes
are available:
palette = 1: Haxby colors [default]
palette = 2: high Intensity colors
palette = 3: low Intensity colors
palette = 4: grayscale
palette = 5: uniform grayscale
A complete description of the color palettes is given
in the COMPLETE DESCRIPTION OF OPTIONS section below.
The ncolors parameter sets the number of color
values used in plotting, whether the colors are
represented in a continuous color scale or a
stepped, discrete color scale [default is 11].
If the option argument is the path to an existing GMT
color palette (CPT) file, then that CPT file and its
color scheme will be used for the plot
COMPLETE DESCRIPTION OF OPTIONS
- -A
-
magnitude/azimuth or magnitude/median
Sets the parameters which control how mbm_plot generates
simulated illumination of bathymetry, which can be either
shaded relief bathymetry or bathymetry draped with amplitude data.
If mode is set to 2 (shaded relief bathymetry) using the
-G option, then the value magnitude
is an effective vertical exageration which modulates the intensity of
the shading; typical values are in the 1-5 range. The value azimuth
is the azimuth from which the bathymetry is illuminated.
If mode is set to 3 (bathymetry shaded using amplitudes) using the
-G option, then the value magnitude
modulates the intensity of the shading; the value median sets the
amplitude value which serves as the zero or neutral level.
- -B
-
tickinfo
Sets map boundary tickmark intervals. See the psbasemap
manual page for details. By default the program chooses
basemap annotations based on the map boundaries.
- -C
-
[cont_int/col_int/tic_int/lab_int/tic_len/lab_hgt]
Given by itself, the -C option produces a contour plot
of the bathymetry data with a contour interval chosen according
to the data in the file or files. If the -G option is used,
the default contours will be drawn in black; otherwise, the
default contours will be drawn in four colors (black, red,
green, and blue) with color changes, annotations, and tickmarks
every fourth contour interval. If any of the optional
parameters are appended, these values will control the contour interval
and other contour characteristics.
Contours will be
generated at invervals of cont_int meters. Color changes
will occur at intervals of col_int meters. Contours will have
downhill facing tickmarks tic_len inches long every tic_int
meters. Contours will have annotations
lab_hgt inches high every lab_int meters.
- -D
-
[flipcolor/flipshade]
This option flips the color and shading conventions used by
mbm_plot.
Normally, the color or grayscale tables used for color bathymetry
maps run from
cool colors (or dark grays) for large depth values
to hot colors (or light grays) for small depth values.
In contrast, sidescan and beam amplitude data is normally plotted
using (light grays (or hot colors) for small amplitudes and
dark grays (or cool colors) for large amplitudes.
If -D is given alone or with flipcolor = 1,
it applies to the color table used
for color or gray fill plots, shaded or unshaded. If the plot
is to be shaded, either by synthetic illumination (-G2)
or by overlaying amplitude data (-G3 option), then
setting flipshade = 1 will cause the shading convention
to be reversed (e.g. high amplitudes overlaid as light shading).
Using -D0/1 will flip the shading convention
but leave the default color convention.
- -F
-
format
Sets the data format for the input data.
If format < 0, then the input file specified
with the -I option will actually contain a list of input swath sonar
data files. This program uses the MBIO library
and will read or write any swath sonar
format supported by MBIO. A list of the swath sonar data formats
currently supported by MBIO and their identifier values
is given in the MBIO manual page. Default: format = -1.
- -G
-
mode
Turns on color fill swath plot and sets the style of the plot.
mode = 1: Color fill of bathymetry data.
mode = 2: Color shaded relief bathymetry.
mode = 3: Bathymetry shaded using amplitude data.
mode = 4: Grayscale fill of amplitude data.
mode = 5: Grayscale fill of sidescan data.
- -H
-
This "help" flag cause the program to print out a description
of its operation and then exit immediately.
- -I
-
filename
Sets the input filename. If format > 0 (set with the
-f option) then the swath sonar data contained in infile
is read and processed. If format < 0 (the default),
then infile
is assumed to be an ascii file containing a list of the input swath sonar
data files to be processed and their formats. The program will read
the data in each one of these files.
In the infile file, each
data file should be followed by a data format identifier, e.g.:
datafile1 11
datafile2 24
This program uses the MBIO library and will read or write any swath sonar
format supported by MBIO. A list of the swath sonar data formats
currently supported by MBIO and their identifier values
is given in the MBIO manual page.
- -J
-
projection[/scale | /width]
Selects the map projection. By default the map projection is
Mercator and the plot scale is chosen to fit on the selected
page size (see -P option). The user may specify a
different projection to be used, in which case the plot scale
is still automatically chosen to fit the page. The user may
also specify both the projection and the plot scale. If
the projection specifying character is upper case, a plot
width rather than a plot scale is used.
The scale values are specified in inch/degree or in 1:xxxxx
ratios. Plot widths are specified in inches. If the user
specifies a plot scale such that the plot will not fit
on the default A size page, a appropriately larger page
size will be chosen.
CYLINDRICAL PROJECTIONS:
-Jclon0/lat0/scale (Cassini)
-Jmscale (Mercator)
-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)
-Juzone/scale (UTM - Universal Transverse Mercator)
-Jylon0/lats/scale (Basic Cylindrical Projection)
AZIMUTHAL PROJECTIONS:
-Jalon0/lat0/scale (Lambert).
-Jelon0/lat0/scale (Equidistant).
-Jglon0/lat0/scale (Orthographic).
-Jslon0/lat0/scale (General Stereographic)
CONIC PROJECTIONS:
-Jblon0/lat0/lat1/lat2/scale (Albers)
-Jllon0/lat0/lat1/lat2/scale (Lambert)
MISCELLANEOUS PROJECTIONS:
-Jhlon0/scale (Hammer)
-Jilon0/scale (Sinusoidal)
-Jklon0/scale (Eckert VI)
-Jnlon0/scale (Robinson)
-Jrlon0/scale (Winkel Tripel)
-Jwlon0/scale (Mollweide)
NON-GEOGRAPHICAL PROJECTIONS:
-Jpscale (Linear projection for polar (theta,r) coordinates)
-Jxx-scale[l|ppow][/y-scale[l|ppow]] (Linear, log, and power scaling)
More details can be found in the psbasemap manpages.
- -L
-
title:scalelabel
Sets the title and the label for the colorscale (if used) of
the plot. Note that a colon (:) rather than a slash (/) is
used to separate the labels. Colons cannot be used in the
labels themselves. If this option is not used, then a default title
and colorscale label are provided. If the title is supplied
alone, a default colorscale label will be provided. To force
no title use -L" "; to force no title or colorscale
label use -L" : ".
- -M
-
A series of "miscellaneous" options are provided which are
given as -M followed by a two character identifier, followed
by any other parameters associated with that option.
The -M options may be strung together separated by
colons, e.g. "-MGQ100:GU", which is equivalent to
"-MGQ -MGU".
- -MGD
-
gmtdef/value
Allows the user to set the GMT default values used as
the plot is constructed. This command may be given repeatedly
to set as many GMT defaults as required. For example, to
set the basemap annotation font to Courier, use
"-MGDANOT_FONT/Courier".
- -MGF
-
scale_loc
Sets the location of the color scale. The possible values
of scale_loc are:
scale_loc = b: bottom of plot
scale_loc = t: top of plot
scale_loc = l: left of plot
scale_loc = r: right of plot
[Default scale_loc = b]
- -MGL
-
[f][x]lon0/lat0/slat/length[m]
Draws a simple map scale centered on lon0/lat0. Use -Lx to specify position in inch instead.
Scale is calculated at latitude slat, length is in km [miles if m is appended].
Use -Lf to get a "fancy" scale [Default is plain].
- -MGQ
-
dpi
Sets the resolution in dots per inch of the raster image used
for color fill maps. Larger values of dpi produce larger
Postscript plot files. [Default is 100].
- -MGT
-
x/y/size/angle/font/just/text
Causes a text label to plotted on the map.
size is text size in points, angle is
measured in degrees counter-clockwise from horizontal,
fontno sets the font type, justify sets the alignment.
If fontno starts with a leading hyphen, then
the remainder of fontno is taken to be a
textstring with the desired fontname. See the
gmtdefaults man page for names and numbers of available fonts
(or run pstext -L). The alignment
number refers to the part of the textstring that will be mapped
onto the (x,y) point: 1 = Lower Left corner,
2 = Lower Center, 3 = Lower Right, 5 = Mid Left, 6 = Mid Center,
7 = Mid Right, 9 = Upper Left, 10 = Upper Center,
11 = Upper Right. This option may be given as many times as
needed.
- -MGU
-
[/dx/dy/][label]
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.75,-0.75)]. Optionally,
append a label, or c (which will plot the command string.)
- -MMAfactor/mode/depth
-
This option determines how the along-track dimension of the
beam or pixel footprints is calculated. If mode = 1,
then the fore-aft beam angle width of the sonar is used so that
the width increases towards the outer parts of the swath.
The fore-aft beam angle width (MB-System internally stores
a value for each format/sonar) is multiplied by the factor
value; a factor < 1.0 can be useful if the data highly
oversamples the seafloor and a factor > 1.0 can fill in
plots of data which undersample the seafloor. If the data
stream does not include depth values (e.g. one is plotting
pure sidescan data), then the depth value sets the
depth value in meters used in the footprint calculations.
If mode = 2, then the along-track dimension of the beam
or pixel footprints is just the along-track distance between
pings multiplied by the factor value.
If mode = 3, then each data point is
plotted as a point, and the factor and depth parameters
are ignored.
Default: factor = 1.0, mode = 1, depth = 3000.0.
- -MMByr/mo/da/hr/mn/sc
-
Sets the starting time for data allowed in the input data; pings
with times before the starting time will be ignored.
Default: yr/mo/da/hr/mn/sc = 1962/2/21/10/30/0.
- -MMDmode/scale[/min/max]
-
Sets scaling of beam amplitude or sidescan pixel values which
can be applied before plotting. If mode = 1 or 2, then
a linear scaling of the form:
scaled_value = scale * (value - min) / (max - min)
is applied. If mode = 3 or 4, then a log10 scaling of
the form:
scaled_value = scale * (20 * log10(value) - min) / (max - min)
is applied. If mode = 2 or 4, then the value (or 20*log10(value))
will be clipped to min if it is smaller than min or max
if it is greater than max; this clipping happens prior to the
multiplication by scale. Default: mode = 1, scale = 1.0,
min = 0.0, max = 1.0.
- -MMEyr/mo/da/hr/mn/sc
-
Sets the ending time for data allowed in the input data; pings
with times after the ending time will be ignored.
Default: yr/mo/da/hr/mn/sc = 2062/2/21/10/30/0.
- -MMLlonflip
-
If the -R option is not used to explicitly set the plot bounds, then
the lonflip value sets the range of the longitude values used for
calculating the desired bounds. If lonflip=-1
then the longitude values will be in the range from -360 to 0
degrees. If lonflip=0 then the longitude values will be in the
range from -180 to 180 degrees. If lonflip=1 then the longitude
values will be in the range from 0 to 360 degrees. Default:
mbm_plot uses the lonflip value set by mbdefaults.
- -MMNnplot
-
Sets the number of pings to be read in before each contouring
episode. See the description of the -MMZalgorithm option
for advice on reasonable values
Default: nplot = 50 unless -MMZ1 is specified,
in which case the default is nplot = 5.
- -MMPpings
-
Sets the ping averaging of the input swath sonar data. If pings = 1, then
no ping averaging is performed. If pings > 0, then
that number of input pings will be averaged to produce one output
ping. If pings = 0, then the ping averaging will automatically
be done so that the along-track ping spacing is equal to the across-track
beam spacing.
Default: pings = 1 (no ping averaging).
- -MMSspeedmin
-
Sets the minimum speed in km/hr (5.5 kts ~ 10 km/hr) allowed in
the input data; pings associated with a smaller ship speed will not be
processed. Default: speed = 0.
- -MMT
-
timegap
Sets the maximum time gap in minutes between adjacent pings allowed before
the data is considered to have a gap. Default: timegap = 1.
- -MMZ
-
algorithm
Sets the contouring algorithm to be used. If algorithm=0,
a simple ping to ping contouring approach is used; this algorithm
is fast but produces poor looking contours when used with data
where beams from one ping may lie "behind" beams from previous
pings (this happens for sonars that ping at nonnull pitch
angles or for the "inside" beams when ships make sharp turns).
If algorithm=1 then a triangular network is constructed from
the available soundings and this network is in turn contoured;
this algorithm is slow but produces good looking contours in
most cases. It is important to note that the time required for
"triangle" algorithm increases with the square of the number
of beams to be contoured; thus it is sensible to keep the number
of pings contoured at a time small (e.g. use -N5).
The time required for the "ping to ping" algorithm varies linearly
with the number of pings contoured; thus larger numbers of pings
may be reasonably contoured at a time (e.g. use -N50).
Default: algorithm = 0 unless format = 41.
- -MNA
-
[nhgt[/P] | P]
Turns on filename annotation of navigation tracks. If -MNA is given
without specifying any controlling parameters, then the lettering height
nhgt is 0.15 and the filenames are plotted parallel to the
navigation track from the start of the track. The lettering height can be
specified using either -MNAnhgt or -MNAnhgt/P.
If -MNAP or -MNAnhgt/P is specified,
the filename will be plotted perpendicular to the navigation track.
Filename annotation can also be specified using the -N option.
Defaults: Filename annotation off.
- -MNP
-
[pingnumber_tick/pingnumber_annot/pingnumber_tick_len]
Turns on ping number (or shot number) annotation of navigation tracks.
Tick marks are made along the shiptrack at pingnumber_tick intervals; these
are tlen inches long. Longer tick marks are made along the
shiptrack at pingnumber_annot intervals; these are 1.5 times tlen
inches long.
Defaults: Pingnumber annotation off. If the -MNP option is given without
specifying the controlling parameters, then pingnumber_tick = 50,
pingnumber_annot = 100, and pingnumber_tick_len = 0.1.
- -MTC
-
fill
Coastline plotting option.
Set the shade (0-255), color (r/g/b), or pattern
(p|Pdpi/pattern; see -MTG) for lakes [Default is the
fill chosen for "wet" areas (-S)].
- -MTD
-
resolution
Coastline plotting option.
Selects the resolution of the coastline data set to use ((f)ull,
(h)igh, (i)ntermediate, (1)ow, and (c)rude). The
resolution drops off by 80% between data sets. [Default
is l].
- -MTG
-
fill
Coastline plotting option.
Select painting or clipping of "dry" areas. Append a
shade, color, pattern, or c for clipping. Specify the
shade (0-255) or color (r/g/b), or -MTGpdpi/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. dpi sets the resolution of the image. See GMT
Cookbook & Technical Reference Appendix E for
information on individual patterns.
- -MTI
-
river[/pen]
Coastline plotting option.
Draw rivers. Specify the type of rivers and
[optionally] append pen attributes [Default pen:
width = 1, color = 0/0/0, texture = solid]. Choose
from the list of river types below. Repeat option -I
as often as necessary.
1 = Permanent major rivers
2 = Additional major rivers
3 = Additional rivers
4 = Minor rivers
5 = Intermittent rivers - major
6 = Intermittent rivers - additional
7 = Intermittent rivers - minor
8 = Major canals
9 = Minor canals
10 = Irrigation canals
a = All rivers and canals (1-10)
r = All permanent rivers (1-4)
i = All intermittent rivers (5-7)
c = All canals (8-10)
- -MTN
-
border[/pen]
Coastline plotting option.
Draw political boundaries. Specify the type of
boundary and [optionally] append pen attributes
[Default pen: width = 1, color = 0/0/0, texture =
solid]. Choose from the list of boundaries below.
Repeat option -MTN as often as necessary.
1 = National boundaries
2 = State boundaries within the Americas
3 = Marine boundaries
a = All boundaries (1-3)
- -MTS
-
fill
Coastline plotting option.
Select painting or clipping of "wet" areas. Append the
shade (0-255), color (r/g/b), pattern (see -MTG), or c
for clipping.
- -MTW
-
pen
Coastline plotting option.
Append pen attributes [Defaults: width = 1, color = 0/0/0,
texture = solid].
- -MXG
-
fill
Select filling of symbols for xy plotting.
Set the shade (0-255) or color
(r/g/b) [Default is no fill]. To reset no fill,
use fill = "N".
For polygons, you may optionally specify
-Gpicon_size/pattern, where
pattern gives the number of the
image pattern (1-32) OR the name of a
icon-format file. icon_size sets
the unit size in inch.
To invert black and white pixels, use
-GP instead of -Gp. See
GMTs Cookbook & Technical Reference
Appendix E for information on individual patterns.
- -MXI
-
xy_file
Specifies a file containing (x,y) pairs to be plotted
as lines or symbols. The line and symbol characteristics
are set using the last -MXG, -MXS, and -MXW
options used. All of the -MX commands can be
given multiple times, so by stringing series of these
commands together the user can plot different files
using different line or symbol characteristics.
[Default is a solid black line].
- -MXM
-
Toggles expectation for xy data files having multiple
segments, in which each segment is to be plotted
separately. Segments are separated by a
record whose first character is '>'. By default,
unsegmented files are expected. Users may give this
command multiple times, allowing some input files to
be handled as segmented and others not.
- -MXS
-
symbol/size
Selects symbol to be used for plotting the next xy data
file. Setting symbol = "N" causes line plotting.
Choose between:
- -MXSa
-
star. size is radius of circumscribing circle.
- -MXSb
-
bar extending from base to y. size is bar width. By default,
base = 0. Append /base to change this value. Append u if size
is in x-units [Default is inch].
- -MXSc
-
circle. size is diameter of circle.
- -MXSd
-
diamond. size is side of diamond.
- -MXSe
-
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.
- -MXSf
-
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.
- -MXSh
-
hexagon. Give side in inch.
- -MXSi
-
inverted triangle. Give side in inch.
- -MXSl
-
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 *.
- -MXSp
-
point. No size needs to be specified (1 pixel is used).
- -MXSs
-
square. Give side in inch.
- -MXSt
-
triangle. Give side in inch.
- -MXSv
-
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.03/0.12/0.1 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.
- -MXSV
-
Same as -MXSv, 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 (-MXSv leaves the directions
unchanged.)
- -MXSx
-
cross. Give length in inch.
- -MXW
-
pen
Set pen attributes for xy plotting. See chapter 4.12 in the
GMT Technical reference for a discussion of GMT pen values.
[Defaults: width = 1, color = 0/0/0,
texture = solid].
- -N
-
[ttick/tannot/dannot/tlen[/nhgt/nperp]] | F | FP]
This option causes a navigation track plot to be generated, and
can also set the start of each swath file to be annotated with the
filename. If the -N option is given alone, then the
navigation track will be plotted without any annotation, The optional
parameters allow users to control the details of the navigation
track annotation. Time marks are made with "X" marks along the
shiptrack; annotated time marks show the time in HH:MM format
next to the time mark and annotated date marks show the time
and julian day in HH:MM/DDD format. The "X" marks are
tlen inches high for normal time marks and
1.5 times tlen inches high for annotated time or date
marks. The interval of time ticks, annotated time ticks, and
annotated date ticks are given in hours by ttick,
tannot, and dannot, respectively.
If the nhgt parameter is not given when the other parameters
are specified, then no filename annotation will be done. If given,
nhgt sets the height in inches of the filename annotation
and turns that annotation on. If given as 1, nperp causes
the filename annotation to be perpendicular to the shiptrack rather
than parallel (the default). If the -NF is given, then
a navigation track will be generated using the default parameters
and also with filename annotation along the shiptrack. If the
-NFP is given, then a navigation track will be generated
with the default parameters and also with filename annotation
perpendicular to the shiptrack.
Defaults if annotation is enabled: ttick = 0.25;
tannot = 1.0; dannot = 4.0;
tlen = 0.1; nhgt = 0.1; nperp = 0.
- -O
-
root
Sets the root used to construct the filename of the output shellscript
(root.cmd) and names of files created when the shellscript is
run. By default, the
name of the input data file or list file is used as the root.
- -P
-
pagesize
This option sets the size of the page the plot will be centered
on. If the user does not set the plot scale, the plot will be
sized as large as will fit on the designated page. If the user
sets the plot scale such that the plot will not fit on the
designated page, a larger page will be used.
The supported page sizes are:
American ANSI sizes:
A 8.5 x 11.0 in. ( 215.9 x 279.4 mm)
B 11.0 x 17.0 in. ( 279.4 x 431.8 mm)
C 17.0 x 22.0 in. ( 431.8 x 558.8 mm)
D 22.0 x 34.0 in. ( 558.8 x 863.6 mm)
E 34.0 x 44.0 in. ( 863.6 x 1117.6 mm)
F 28.0 x 40.0 in. ( 711.2 x 1016.0 mm)
E1 44.0 x 68.0 in. (1117.6 x 1727.2 mm)
Metric ISO A sizes:
A0 841.0 x 1189.0 mm (33.11 x 46.81 in.)
A1 594.0 x 841.0 mm (23.39 x 33.11 in.)
A2 420.0 x 594.0 mm (16.54 x 23.39 in.)
A3 297.0 x 420.0 mm (11.69 x 16.54 in.)
A4 210.0 x 297.0 mm ( 8.27 x 11.69 in.)
A5 148.0 x 210.0 mm ( 5.83 x 8.27 in.)
A6 105.0 x 148.0 mm ( 4.13 x 5.83 in.)
A7 74.0 x 105.0 mm ( 2.91 x 4.13 in.)
A8 52.0 x 74.0 mm ( 2.05 x 2.91 in.)
A9 37.0 x 52.0 mm ( 1.46 x 2.05 in.)
A10 26.0 x 37.0 mm ( 1.02 x 1.46 in.)
Metric ISO B sizes:
B0 1000.0x 1414.0 mm (39.37 x 55.67 in.)
B1 707.0 x 1000.0 mm (27.83 x 39.37 in.)
B2 500.0 x 707.0 mm (19.68 x 27.83 in.)
B3 353.0 x 500.0 mm (13.90 x 19.68 in.)
B4 250.0 x 353.0 mm ( 9.84 x 13.90 in.)
B5 176.0 x 250.0 mm ( 6.93 x 9.84 in.)
B6 125.0 x 176.0 mm ( 4.92 x 6.93 in.)
B7 88.0 x 125.0 mm ( 3.46 x 4.92 in.)
B8 62.0 x 88.0 mm ( 2.44 x 3.46 in.)
B9 44.0 x 62.0 mm ( 1.73 x 2.44 in.)
B10 31.0 x 44.0 mm ( 1.22 x 1.73 in.)
Metric ISO C sizes:
C0 914.4 x 1300.5 mm (36.00 x 51.20 in.)
C1 650.2 x 914.4 mm (25.60 x 36.00 in.)
C2 457.2 x 650.2 mm (18.00 x 25.60 in.)
C3 325.1 x 457.2 mm (12.80 x 18.00 in.)
C4 228.6 x 325.1 mm ( 9.00 x 12.80 in.)
C5 162.6 x 228.6 mm ( 6.40 x 9.00 in.)
C6 114.3 x 162.6 mm ( 4.50 x 6.40 in.)
C7 81.3 x 114.3 mm ( 3.20 x 4.50 in.)
MB-System large format sizes:
m1 1371.6 x 1828.8 mm (54.00 x 72.00 in.)
m2 1371.6 x 2133.6 mm (54.00 x 84.00 in.)
m3 1371.6 x 2438.4 mm (54.00 x 96.00 in.)
m4 1524.0 x 1828.8 mm (60.00 x 72.00 in.)
m5 1524.0 x 2133.6 mm (60.00 x 84.00 in.)
m6 1524.0 x 2438.4 mm (60.00 x 96.00 in.)
The default page size is A.
- -Q
-
Normally, the output plot generation shellscript
includes lines which execute
a program to display the Postscript image on the screen.
This option causes those lines to be commented out so
that executing the shellscript produces a Postscript plot
but does not attempt to display it on the screen.
The program
to be used to display the Postscript is set
using mbdefaults;
the default value can be overridden by setting the environment
variable $MB_PS_VIEWER.
- -R
-
west/east/south/north
Sets the longitude and latitude bounds within which swath sonar
data will be read. Normally the bounds are automatically chosen
to include all of the input data.
- -S
-
[color/shade]
This option enables effective histogram equalization of the
color and/or shading of thedata. The equalization is
not achieved by changing the data values, but rather by
constructing the color or shading tables so that
the boundaries in the tables encompass equal fractions of
the datapoints. This serves to focus color or shading contrasts
in value ranges corresponding to the bulk of the data values,
and is particularly useful for
enhancing the contrast of sidescan and beam amplitude plots.
If -S is given alone or with color = 1,
it enables equalization of the color table used
for color or gray fill plots, shaded or unshaded. If the plot
is to be shaded, either by synthetic illumination (-G2)
or by overlaying amplitude data (-G3 option), then
setting shade = 1 will cause the shading to be equalized.
Using -S0/1 will equalize the shading without
equalizing the color table.
- -T
-
If -T is given, it causes a coastline to be drawn
on the map. The default is to draw the coastline and shade
all dry land a uniform gray. To exercise greater control of
the coastline plotting, use the
-MTC, -MTD, -MTG, -MTI,
-MTN, -MTS, and -MTW options described in the
- -U
-
orientation
Normally the orientation of the plot (portrait or landscape)
is selected automatically so as to maximize the plot scale.
The -U option allows the user to set the plot orientation. If
orientation = 1, a portrait plot will be produced; if
orientation = 2, a landscape plot will be produced.
- -V
-
Causes mbm_plot to operate in "verbose" mode so that it
outputs more information than usual.
- -W
-
[color_style[/palette[ncolors]] | cptfile]
This option controls the color scheme used for color
fill plots.
If color_style = 1 [default], then
the color scheme used will be a continuous grading
of colors. If color_style = 2, the color scheme
will be a set of discrete color intervals. The color
palette used is set using palette. Seven palettes
are available:
palette = 1: Haxby colors [default]
palette = 2: high Intensity colors
palette = 3: low Intensity colors
palette = 4: grayscale
palette = 5: uniform grayscale
palette = 6: uniform black
palette = 7: uniform white
The RGB definitions of the color palettes are:
color palette 1 - Haxby Color Table
red: 255 255 255 255 240 205 138 106 50 40 37
green: 255 186 161 189 236 255 236 235 190 127 57
blue: 255 133 68 87 121 162 174 255 255 251 175
color palette 2 - High Intensity Colors
red: 255 255 255 255 128 0 0 0 0 128 255
green: 0 64 128 255 255 255 255 128 0 0 0
blue: 0 0 0 0 0 0 255 255 255 255 255
color palette 3 - Low Intensity Colors
red: 200 194 179 141 90 0 0 0 0 90 141
green: 0 49 90 141 179 200 141 90 0 0 0
blue: 0 0 0 0 0 0 141 179 200 179 141
color palette 4 - Grayscale
red: 255 230 204 179 153 128 102 77 51 26 0
green: 255 230 204 179 153 128 102 77 51 26 0
blue: 255 230 204 179 153 128 102 77 51 26 0
color palette 5 - Uniform Grayscale
red: 128 128 128 128 128 128 128 128 128 128 128
green: 128 128 128 128 128 128 128 128 128 128 128
blue: 128 128 128 128 128 128 128 128 128 128 128
color palette 6 - Uniform Black
red: 0 0 0 0 0 0 0 0 0 0 0
green: 0 0 0 0 0 0 0 0 0 0 0
blue: 0 0 0 0 0 0 0 0 0 0 0
color palette 7 - Uniform White
red: 255 255 255 255 255 255 255 255 255 255 255
green: 255 255 255 255 255 255 255 255 255 255 255
blue: 255 255 255 255 255 255 255 255 255 255 255
The Haxby colors have been adapted from a palette
developed by Dr. William Haxby of the Lamont-Doherty
Earth Observatory; this palette is pleasing to the
eye and well suited for shading. The high intensity
colors describe linear paths through RGB space from
red to blue to green to purple; because the colors are high
intensity they are not well suited to shading.
The low intensity colors are similar to the high
intensity, but muted and thus well suited to shading.
The grayscale palette runs linearly from white to
black and is commonly used for plots of sidescan and amplitude
data. The uniform grayscale is useful for non-color
shaded relief plots.
The ncolors parameter sets the number of color
values used in plotting, whether the colors are
represented in a continuous color scale or a
stepped, discrete color scale [default is 11].
If the option argument is the path to an existing GMT
color palette (CPT) file, then that CPT file and its
color scheme will be used for the plot
- -X
-
Normally, mbm_plot creates an executable shellscript and
then exits. This option will cause the shellscript to be executed
in the background before mbm_plot exits.
- -Y
-
Normally, mbm_plot generates nicely rounded numbers
for the boundaries of the color palette. Often, the resulting
color bounds extend well outside the range of the gridded data.
This option causes
the color boundaries to be uniformly distributed between the
minimum and maximum values of the grid.
- -Z
-
min/max
This option overrides the minimum and maximum values of
bathymetry data, affecting the color palette and the
contour interval if those parameters are not specified
by the user.
EXAMPLES
Suppose we have obtained a swath sonar data file called
sb2112_example.mb41 collected using a SeaBeam 2112 sonar.
This file contains bathymetry, beam amplitude, and
sidescan data. In order to obtain initial views of the
data in the file, we use mbm_plot to generate
shellscripts which in turn generate plots when executed.
The following five commands generate plotting shellscripts
for color fill bathymetry overlaid with contours, color shaded
relief bathymetry, color fill bathymetry overlaid with
amplitudes, grayscale amplitudes, and grayscale sidescan,
respectively:
mbm_plot -F41 -I sb2112_example.mb41 -G1 -C \
-N -V -Obathcont
mbm_plot -F41 -I sb2112_example.mb41 -G2 \
-N -V -Obathshade
mbm_plot -F41 -I sb2112_example.mb41 -G3 -S0/1 \
-N -V -Obathamp
mbm_plot -F41 -I sb2112_example.mb41 -G4 -S \
-N -V -Oamp
mbm_plot -F41 -I sb2112_example.mb41 -G5 -S \
-N -V -Oss
When the following shellscripts are executed, each will
generate a postscript plot file and then display the plot
on the screen:
bathcont.cmd
bathshade.cmd
bathamp.cmd
amp.cmd
ss.cmd
Note that we use the -S option to apply histogram
equalization to the amplitude and sidescan data, but not
the bathymetry data. Also note that by specifying -N
we obtain a track plot of the ship's navigation overlaid
on the color or grayscale file plots.
Now suppose we have a set of SeaBeam 2112 data files
comprising a short survey and that we want a plot
of all the data together. We create an ASCII text file which has
a list of the filenames, each followed by the appropriate MBIO
format id number, e.g.:
sb2112_example_1.mb41 41
sb2112_example_2.mb41 41
sb2112_example_3.mb41 41
If the name of the data file list is "datalist", then
using "-F-1 -Idatalist" will
cause the macro to operate on all of the files
together. We desire a plot including color fill bathymetry
overlaid with 25 meter contours and the ship's navigation.
For this plot, we want a very bright colortable and we would like
the colors to be discretely stepped rather than continuous;
thus we use -W2/2. We also choose to use the
-X option, which causes mbm_plot to execute
the shellscript it creates in the background before
exiting. The mbm_plot command is:
mbm_plot -F-1 -Idatalist -G1 -C25 -N \
-X -V -Obathtest
As an example, the contents of the plotting shellscript
"bathtest.cmd" are:
#
# Shellscript to create Postscript plot of swath sonar data
# Created by macro mbm_plot
#
# This shellscript created by following command line:
# mbm_plot -F-1 -Idatalist -G1 -C25 -N -V -Obathtest
#
# Save existing GMT defaults
echo Saving GMT defaults...
gmtdefaults -L > gmtdefaults$$
#
# Set new GMT defaults
echo Setting new GMT defaults...
gmtset ANOT_FONT Helvetica
gmtset LABEL_FONT Helvetica
gmtset HEADER_FONT Helvetica
gmtset ANOT_FONT_SIZE 8
gmtset LABEL_FONT_SIZE 8
gmtset HEADER_FONT_SIZE 10
gmtset FRAME_WIDTH 0.074999999999999997
gmtset TICK_LENGTH 0.074999999999999997
gmtset PAGE_ORIENTATION LANDSCAPE
gmtset COLOR_BACKGROUND 0/0/0
gmtset COLOR_FOREGROUND 255/255/255
gmtset COLOR_NAN 255/255/255
#
# Make color palette table file
echo Making color palette table file...
echo 2975 255 255 255 3150 255 186 133 > bathtest.cpt
echo 3150 255 186 133 3325 255 161 68 >> bathtest.cpt
echo 3325 255 161 68 3500 255 189 87 >> bathtest.cpt
echo 3500 255 189 87 3675 240 236 121 >> bathtest.cpt
echo 3675 240 236 121 3850 205 255 162 >> bathtest.cpt
echo 3850 205 255 162 4025 138 236 174 >> bathtest.cpt
echo 4025 138 236 174 4200 106 235 255 >> bathtest.cpt
echo 4200 106 235 255 4375 50 190 255 >> bathtest.cpt
echo 4375 50 190 255 4550 40 127 251 >> bathtest.cpt
echo 4550 40 127 251 4725 37 57 175 >> bathtest.cpt
#
# Run mbswath
echo Running mbswath...
mbswath -f-1 -Idatalist \
-Jm22.198543775528325 \
-R114.210795/114.430905/-31.91322/-31.62458 \
-Cbathtest.cpt \
-p1 -A1 -Z1 \
-p1 \
-P -X1.8069392647842304 -Y2 -K -V > bathtest.ps
#
# Run mbcontour
echo Running mbcontour...
mbcontour -f-1 -Idatalist \
-Jm22.198543775528325 \
-R114.210795/114.430905/-31.91322/-31.62458 \
-A50/100000/100000/100000/0.01/0.1 \
-D0.25/1/4/0.15 \
-p1 \
-P -K -O -V >> bathtest.ps
#
# Make color scale
echo Running psscale...
psscale -Cbathtest.cpt \
-D2.4431/-0.5000/4.8861/0.1500h \
-B":.Depth (meters):" \
-P -K -O -V >> bathtest.ps
#
# Make basemap
echo Running psbasemap...
psbasemap -Jm22.198543775528325 \
-R114.210795/114.430905/-31.91322/-31.62458 \
-B5m/5m:."Data List File datalist": \
-P -O -V >> bathtest.ps
#
# Delete surplus files
echo Deleting surplus files...
rm -f bathtest.cpt
#
# Reset GMT default fonts
echo Resetting GMT fonts...
mv gmtdefaults$$ .gmtdefaults
#
# Run xpsview
echo Running xpsview in background...
xpsview -ps a -maxp 4m bathtest.ps &
#
# All done!
echo All done!
SEE ALSO
mbsystem(1), mbcontour(1),
mbswath(1), mbdefaults(1),
mbm_grdplot(1), mbm_grd3dplot(1), mbfilter(1)
BUGS
By making this macro more useful, we have also made it
more complex.
Index
- NAME
-
- VERSION
-
- SYNOPSIS
-
- DESCRIPTION
-
- AUTHORSHIP
-
- SIMPLE DESCRIPTION OF BASIC OPTIONS
-
- COMPLETE DESCRIPTION OF OPTIONS
-
- EXAMPLES
-
- SEE ALSO
-
- BUGS
-
Last Updated: 3 June 2013
Return to list of MB-System manual pages...
Back
to MB-System Home Page...