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 sfc
st_shift_longitude(x, ...)

# S3 method for sf
st_shift_longitude(x, ...)

Arguments

x

object of class sf or sfc

...

ignored

Examples

## sfc pt1 = st_point(c(-170, 50)) pt2 = st_point(c(170, 50)) (sfc = st_sfc(pt1, pt2))
#> Geometry set for 2 features #> geometry type: POINT #> dimension: XY #> bbox: xmin: -170 ymin: 50 xmax: 170 ymax: 50 #> epsg (SRID): NA #> proj4string: NA
#> POINT (-170 50)
#> POINT (170 50)
sfc = st_set_crs(sfc, 4326) st_shift_longitude(sfc)
#> Geometry set for 2 features #> geometry type: POINT #> dimension: XY #> bbox: xmin: 170 ymin: 50 xmax: 190 ymax: 50 #> epsg (SRID): 4326 #> proj4string: +proj=longlat +datum=WGS84 +no_defs
#> POINT (190 50)
#> POINT (170 50)
## sf d = st_as_sf(data.frame(id = 1:2, geometry = sfc)) st_shift_longitude(d)
#> Simple feature collection with 2 features and 1 field #> geometry type: POINT #> dimension: XY #> bbox: xmin: 170 ymin: 50 xmax: 190 ymax: 50 #> epsg (SRID): 4326 #> proj4string: +proj=longlat +datum=WGS84 +no_defs #> id geometry #> 1 1 POINT (190 50) #> 2 2 POINT (170 50)