Hello,

I have a shapefile with polygons located in Oceania, some of which straddle longitude =-180/0. When plotted, horizontal scratches result. I assumed that use of 'nowrapSpatialPolygons' would do the trick, but the plot looks the same before and after. Clearly, I'm not utilizing it properly. Could someone look at the example code below Ultimately, I need to successfully rasterize the data onto a -180/180 grid, so maybe another solution exists...

Thanks!
Kevin

The folder (WDPA Oceania) with the example shapefile can be downloaded here: http://dl.dropbox.com/u/14314000/WDPA%20Oceania.zip

orig=readShapeSpatial("sql_statement",CRS("+proj=longlat"))

plot(orig)

out=nowrapSpatialPolygons(orig,offset=180)

plot(out)






	[[alternative HTML version deleted]]

