par(mar=c(2,2,2,2))
image(Bathymetry$x,Bathymetry$y,Bathymetry$z,col=femmecol(100),
asp=TRUE,xlab="dg",ylab="dg")
contour(Bathymetry$x,Bathymetry$y,Bathymetry$z,asp=TRUE,add=TRUE)
# remove land
zz <- Bathymetry$z
zz[zz>0]<-0
image(Bathymetry$x,Bathymetry$y,zz,col=c(femmecol(100),"black"),asp=TRUE)
contour(Bathymetry$x,Bathymetry$y,zz,asp=TRUE,add=TRUE)
Run the code above in your browser using DataLab