I am using newobj <- nowrapRecenter(obj)
, from the package maptools
to recenter a SpatialPolygonsDataFrame
such that longitude runs 0 - 360
(pacific view) rather than -180
to 180
(Atlantic view). Unfortunately, the resulting newobj
is a SpatialPolygon
only --- the data is lost. I expect I can merge the data back on, but I'm not sure how... generally I would merge by a data column, and this new polygon has no data.
I am using this GEOLEV1 shapefile:
https://international.ipums.org/international/gis_harmonized_1st.shtml