# load NW Atlantic data and convert to class bathy
data(nw.atlantic)
atl <- as.bathy(nw.atlantic)
library("ggplot2")
head(fortify(atl))
# one can now use bathy objects with ggplot directly
ggplot(atl) + geom_contour(aes(x=x, y=y, z=z))
# which allows complete plot configuration
atl.df <- fortify(atl)
ggplot(atl.df, aes(x=x, y=y)) +
geom_raster(aes(fill=z), data=atl.df[atl.df$z <= 0,]) +
geom_contour(aes(z=z),
breaks=c(-100, -200, -500, -1000, -2000, -4000),
colour="white", size=0.1
)
Run the code above in your browser using DataLab