Re: [R-sig-Geo] Rgshhs for getting maps around longitude 0

2010-10-18 Thread Karl Ove Hufthammer
Roger Bivand wrote: Is it necessary to fetch the map twice? If so, which xlim values should be used, and how should the two maps be merged into one SpatialPolygons? Note that if I change ‘xlim’ to c(0,30), Great Britain is missing ... :-( As things stand, I suspect that there is no

Re: [R-sig-Geo] Rgshhs for getting maps around longitude 0

2010-10-17 Thread Roger Bivand
On Thu, 14 Oct 2010, Karl Ove Hufthammer wrote: Dear list members The GSHHS map data uses longitudes from 0 to 360, which causes problems fetching maps around longitude 0. Example: library(maptools) gpclibPermit() gshhs.c.b - system.file(share/gshhs_c.b, package=maptools) xl=c(-30,30)

[R-sig-Geo] Rgshhs for getting maps around longitude 0

2010-10-14 Thread Karl Ove Hufthammer
Dear list members The GSHHS map data uses longitudes from 0 to 360, which causes problems fetching maps around longitude 0. Example: library(maptools) gpclibPermit() gshhs.c.b - system.file(share/gshhs_c.b, package=maptools) xl=c(-30,30) yl=c(50,80) rmap=Rgshhs(gshhs.c.b, xlim=xl, ylim=yl)