Dear all,
I am new to R and even newer to mapping with R. I have installed
PBSMapping, maptools and maps, and I am struggling to produce a map of
fish abundance in the Grand Banks. I have tried several approaches to
plot the bathymetry in my study area:
1) download the data from http://topex.ucsd.edu/cgi-bin/get_data.cgi
(42N46N313W307W)
To start with, the map I get is blank between lon -52W -53W.
Besides, it seems that the shallower isobath is 500 m. Or maybe I am
wrong here?
2) import a shapefile with the bathymetry, but I get the following error:
test <- read.shapefile("C:/data")
Error in seq.default(1, (num.points * 2), by = 2) :
wrong sign in 'by' argument
3) use a Rdata version of the above mentioned shapefile. This is the map
I like the most, but when I try to plot the study area (which is >200 nm
from the nearest shore), I get the error
map('worldHires',xlim=c(-53,-47),ylim=c(42,46))
Error in map.poly(database, regions, exact, xlim, ylim,
boundary, interior, : nothing to draw: all regions out of bounds
The next problem will be plotting fish abundance. I was thinking of
contours or scaled bubbles, but I am not sure I can do any of those with
PBSMapping. Am I right? Any advice?
Thanks for your help,
Elena
[[alternative HTML version deleted]]
