"plot"(x, xlab = "", ylab = "", showHemi = TRUE, asp, clongitude, clatitude, span, lonlabel = NULL, latlabel = NULL, sides = NULL, projection = NULL, expand = 1, mgp = getOption("oceMgp"), mar = c(mgp[1] + 1, mgp[1] + 1, 1, 1), bg, fill, border = NULL, col = "lightgray", axes = TRUE, cex.axis = par("cex.axis"), add = FALSE, inset = FALSE, geographical = 0, longitudelim, latitudelim, debug = getOption("oceDebug"), ...)
coastline
object, as read by read.coastline
or created by as.coastline
, or a list containing items named
longitude
and latitude
.plot,coastline-method
to set the aspect ratio to give natural latitude-longitude scaling somewhere
near the centre latitude on the plot. Often, it makes sense to set
asp
yourself, e.g. to get correct shapes at 45N, use
asp=1/cos(45*pi/180)
. Note that the land mass is not symmetric about
the equator, so to get good world views you should set asp=1
or set
ylim
to be symmetric about zero. Any given value of asp
is
ignored, if clongitude
and clatitude
are given (or
if the latter two are inferred from projection
.clongitude
and clatitude
are provided,
or alternatively if they can be inferred from substrings +lon_0
and +lat_0
in projection
, then
any provided value of asp
is ignored, and instead the plot aspect
ratio is computed based on the center latitude. If clongitude
and
clatitude
are known, then span
must also be provided, and
in this case, it is not permitted to also specify longitudelim
and
latitudelim
.span
must be supplied, if
clongitude
and clatitude
are supplied (or inferred from projection
).plot,coastline-method
. Using these arguments permits reasonably
simple customization. If they are are not provided, reasonable defaults
will be used.mapPlot
argument of the same name).
If set to FALSE
then no projection is used,
and the data are plotted in a cartesion frame, with aspect ratio set to
reduce distortion near the middle of the plot. This option is useful if the
coastline produces spurious horizontal lines owing to islands crossing the
plot edges (a problem that plagues map projections). If projection
is not set, a Mercator projection is used for latitudes below about 70
degrees, as if projection="+proj=merc"
had been supplied, or a
Stereopoloar one is used as if projection="+proj=stere"
. Otherwise,
projection
must be a character string identifying a projection
accepted by mapPlot
.expand
is ignored if either xlim
or ylim
is given.par(mgp)
, and also
for par(mar)
, computed from this. The default is tighter than the R
default, in order to use more space for the data and less for the axes.par("mar")
.polygon
, which does the drawing).polygon
,
which does the drawing.) This is ignored unless the coastline object is
fillable.TRUE
to plot axes.TRUE
to draw the coastline on an existing
plot. Note that this retains the aspect ratio of that existing plot, so it
is important to set that correctly, e.g. with asp=1/cos(lat * pi /
180)
, where clat
is the central latitude of the plot.geographical=0
, the axes are conventional, with decimal degrees as
the unit, and negative signs indicating the southern and western
hemispheres. If geographical=1
, the signs are dropped, with axis
values being in decreasing order within the southern and western
hemispheres. If geographical=2
, the signs are dropped and the axes
are labelled with degrees, minutes and seconds, as appropriate.latitudelim
provide a second way to
suggest plot ranges. Note that these may not be supplied if
clongitude
, clatitude
and span
are given.longitudelim
.TRUE
to get debugging information during
processing.yaxp=c(-90,90,4)
for a plot extending from pole to pole.plot,coastline-method
relied on a
now-defunct argument fill
to control colours; col
is to be
used now, instead. Also, in February, 2016, the arguments named
parameters
and orientation
were both removed, as they had
become nonfunctional about a year previously, in the transition to using
the rgdal
package to carry out map projections.longitudelim
, latitudelim
and projection
are all given,
then these arguments are passed to mapPlot
to produce the plot.
(The call uses bg
for col
, and uses col
, fill
and border
directly.) If the results need further customization,
users should use mapPlot
directly.If projection
is provided without longitudelim
or latitudelim
,
then mapPlot
is still called, but longitudelim
and
latitudelim
are computed from clongitude
, clatitude
and span
.
If projection
is not provided, much simpler plots are produced. These are
Cartesian, with aspect ratio set to minimize shape distortion at the central latitude.
Although these are crude, they have the benefit of always working, which cannot
be said of true map projections, which can be problematic in various ways owing
to difficulties in inverting projection calculations.
To get an inset map inside another map, draw the first map, do
par(new=TRUE)
, and then call plot,coastline-method
with a value of
mar
that moves the inset plot to a desired location on the existing
plot, and with bg="white"
.
coastline-class
explains the
structure of coastline objects, and also outlines the other functions
dealing with them.Other functions that plot oce
data: plot,adp-method
,
plot,adv-method
,
plot,amsr-method
,
plot,argo-method
,
plot,bremen-method
,
plot,cm-method
,
plot,ctd-method
,
plot,gps-method
,
plot,ladp-method
,
plot,lisst-method
,
plot,lobo-method
,
plot,met-method
,
plot,odf-method
,
plot,rsk-method
,
plot,satellite-method
,
plot,sealevel-method
,
plot,section-method
,
plot,tidem-method
,
plot,topo-method
,
plot,windrose-method
,
plotProfile
, plotScan
,
plotTS
, tidem-class
Other things related to coastline
data: [[,coastline-method
,
[[<-,coastline-method
,
as.coastline
,
coastline-class
,
coastlineBest
, coastlineCut
,
coastlineWorld
,
download.coastline
,
read.coastline.openstreetmap
,
read.coastline.shapefile
,
subset,coastline-method
,
summary,coastline-method
## Not run:
# library(oce)
# par(mar=c(2, 2, 1, 1))
# data(coastlineWorld)
# plot(coastlineWorld)
# plot(coastlineWorld, clongitude=-63.6, clatitude=44.6, span=1000)
#
# ## Canada in Lambert projection
# plot(coastlineWorld, clongitude=-95, clatitude=65, span=5500,
# grid=10, projection='+proj=laea +lon_0=-100 +lat_0=55')
# ## End(Not run)
Run the code above in your browser using DataLab