Bathymetry data

From QERM Wiki

This is rough, just trying to put it up there to be cleaned up later.

A slightly older attempt at doing something similar using the NOAA "Design-a-Grid" tool line is posted here: http://students.washington.edu/itaylor/R/#Maps .

Copy and paste the code from here and you get the plot above. The maps and mapdata packages are used to get the aspect ratio right and then add the coastline. Bath_map2.r

#grid data in ASCII raster format with no headers from  http://www.ngdc.noaa.gov/mgg/gdas/gd_designagrid.html
d = read.table('http://wiki.cbr.washington.edu/qerm/images/a/ad/Eli_bath.txt',header=F)
lonvec = seq(130,180,1/15) #4-minute grid
latvec = seq(30,70,1/15)
d2 = d[nrow(d):1,] #why it needs to be flipped over, who knows?

# new and fancier, extending to full range of topo.colors separately for land and water
dmin=min(d2)
dmax=max(d2)
NwaterColors = round(abs(dmin))
NlandColors = round(abs(dmax))

#taking first 1/3 of topo.colors of length NwaterColors for water
waterColorVec = topo.colors(NwaterColors*3)[1:NwaterColors]
#taking second 2/3 of topo.colors of length NlandColors for land
landColorVec = topo.colors(NlandColors*3/2)[round(NlandColors/2+1):round(1.5*NlandColors)]
#combining
allColorVec = c(waterColorVec,landColorVec)

#plot in R
library(maps)
library(mapdata)
map('worldHires',xlim=c(140,173),ylim=c(43,63), xaxs='i',yaxs='i')
image(lonvec,latvec,t(as.matrix(d2)), col=allColorVec, add=T)
map('worldHires',xlim=c(140,173),ylim=c(43,63), xaxs='i',yaxs='i', add=T)
box()
axis(1, seq(140,180,5),paste(seq(140,180,5),'°E', sep=''))
axis(2, seq(30,70,5),paste(seq(30,70,5),'°N', sep=''))