Shift or re-center geographical coordinates for a Pacific view
All longitudes < 0 are added to 360, to avoid for instance parts of Alaska
being represented on the far left and right of a plot because they have
values straddling 180 degrees. In general, using a projected
coordinate reference system is to be preferred, but this method permits a
geographical coordinate reference system to be used. This is the sf
equivalent of recenter
in the sp package and
ST_ShiftLongitude in PostGIS.
st_shift_longitude(x) ## S3 method for class 'sfc' st_shift_longitude(x, ...) ## S3 method for class 'sf' st_shift_longitude(x, ...)
x |
object of class sf or sfc |
... |
ignored |
## sfc pt1 = st_point(c(-170, 50)) pt2 = st_point(c(170, 50)) (sfc = st_sfc(pt1, pt2)) sfc = st_set_crs(sfc, 4326) st_shift_longitude(sfc) ## sf d = st_as_sf(data.frame(id = 1:2, geometry = sfc)) st_shift_longitude(d)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.