plot.btrack: Plots an estimated track by month.

In galuardi/analyzepsat: Microwave Telemtry PSAT analysis

Description

Plots an estimated track by month. May include confidence intervals and bathymetric contours

Usage

plot.btrack(btrack, map, cex = 1.5, ci = F, bathy = NULL, add = F, bathlevels = c(-100, -200),  alpha  =  0.15, bymonth = T, pch = 21, bg = 4)

Arguments

btrack	Estimated track. Track is typically estimated using make.btrack
map	Map object to be used. May be gshhs object created through maptools, or a gshhs coastline extracted through plotmap.r in the kftrack family of libraries
cex	Size of points
ci	Should confidence interval be plotted?
bathy	Should a bathymtric contour be added? The default is to plot the -100m and -200m contours. Color will be drawn from the bath.colors palette
add	Should this plot be added to the existing plot?
bathlevels	Bathymetric contours to be drawn
alpha	If ci = T, then this refers to the degree of transparency in the error ellipses. This is useful for eyeballing heavily used areas.
bymonth	Indicates if the colors of points plotted represent months. This will use plot.by.month
pch	The point type
bg	Background color for the point. This is only used if bymonth = F

Value

A plot of an estimated track.

Author(s)

Benjamin Galuardi

See Also

plot.by.month

Examples

data(ATL)
data(fake)
data(myramps)

bathy = get.bath.data(-80, -50, 30, 50, folder = getwd(), seaonly = T, res = 1)
fit = kftrack(fake[,1:5])

par(mar=c(8,4,4,4))
fmat = prepb(fit, fake)
fmat$Lon_E = fmat$Lon_E-360
plot.btrack(fmat, map = map2, bathy = bathy)

#add scales
plot.btrack(fmat, map = map2)
blue.colors  = colorRampPalette(c("#00007F", "blue", "#007FFF", "cyan"))
#blueg.colors = colorRampPalette(c("grey10", "gray40", "gray50", "gray70", "lightsteelblue4","lightsteelblue3","lightsteelblue2","lightsteelblue1","lemonchiffon1"))
image(bathy$lon+360, bathy$lat, t(bathy$data), add=T, col = blue.colors(100), zlim = c(-10000,0))
plot.btrack(fmat, map = map2, add=T)
plot(map2$SP, col = 'khaki', pbg = 'lightblue', add=T)
.add.month.scale()
.add.bathy.scale('blue') # or 'grey'

galuardi/analyzepsat documentation built on May 17, 2019, 3:25 p.m.