Dear List,

I have been struggling all the day long with this issue but I didn't find
any workaround.

I would like to have exactly the map that come out from the code below but
without the ugly line crossing one side to the other.

world.map <- map(fill=T, plot=F)
IDs <- sapply(strsplit(world.map$names, ":"), function(x) x[1])
poly.map <- map2SpatialPolygons(world.map, IDs=IDs,
proj4string=CRS("+proj=longlat"))
merc.map <- spTransform(poly.map, CRS("+proj=merc"))

# Set the map to -80 to 80
bb <- cbind(c(-180, 180), c(-80, 80))
bb.poly <- SpatialPoints(bb, proj4string = CRS("+proj=longlat"))
bb.merc <- spTransform(bb.poly, CRS("+proj=merc"))
bb.merc <- as.data.frame(bb.merc)

# Here the ugly plot
plot(merc.map, col="gray", border="gray", xlim=bb.merc[,1],
ylim=bb.merc[,2])

Any help is very appreciated,
/Gis

	[[alternative HTML version deleted]]

