Hi All, I'm trying to use the 'coastsCoarse' dataset from library 'rworldmap' as a base map for plotting a global dataset. It works fine when I use unprojected data. However, if I project the data and base map, coastsCoarse ends up with one horizontal line segment that spans the whole map. You can see what I mean like this: # --------------------------------------- library(rgdal) library(rworldmap) data(coastsCoarse) plot(coastsCoarse) # Looks good mercator.proj4 = "+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs" coastsMerc = spTransform(coastsCoarse, CRS(mercator.proj4)) spplot(coastsMerc) # horizontal line across plot connects east and west Beringia... # --------------------------------------- It seems like the boundary of northern Russia didn't get cut off correctly. I imagine the fix is straightforward, but I've been banging away at it to no avail. Does anyone have a simple solution that either (a) locates and removes the offending pan-global segment or (b) doesn't produce it to begin with? Thanks much, Ryan
Projected SpatialLinesDataFrame has a segment spanning the globe
1 message · Ry Guy