Bind rows (features) of sf objects

Bind columns (variables) of sf objects

# S3 method for sf
rbind(..., deparse.level = 1)

# S3 method for sf
cbind(..., deparse.level = 1, sf_column_name = NULL)

st_bind_cols(...)

Arguments

...

objects to bind; note that for the rbind and cbind methods, all objects have to be of class sf; see dotsMethods

deparse.level

integer; see rbind

sf_column_name

character; specifies active geometry; passed on to st_sf

Value

cbind called with multiple sf objects warns about multiple geometry columns present when the geometry column to use is not specified by using argument sf_column_name; see also st_sf.

Details

both rbind and cbind have non-standard method dispatch (see cbind): the rbind or cbind method for sf objects is only called when all arguments to be binded are of class sf.

If you need to cbind e.g. a data.frame to an sf, use data.frame directly and use st_sf on its result, or use bind_cols; see examples.

st_bind_cols is deprecated; use cbind instead.

Examples

crs = st_crs(3857) a = st_sf(a=1, geom = st_sfc(st_point(0:1)), crs = crs) b = st_sf(a=1, geom = st_sfc(st_linestring(matrix(1:4,2))), crs = crs) c = st_sf(a=4, geom = st_sfc(st_multilinestring(list(matrix(1:4,2)))), crs = crs) rbind(a,b,c)
#> Simple feature collection with 3 features and 1 field #> geometry type: GEOMETRY #> dimension: XY #> bbox: xmin: 0 ymin: 1 xmax: 2 ymax: 4 #> epsg (SRID): 3857 #> proj4string: +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs #> a geom #> 1 1 POINT (0 1) #> 2 1 LINESTRING (1 3, 2 4) #> 3 4 MULTILINESTRING ((1 3, 2 4))
rbind(a,b)
#> Simple feature collection with 2 features and 1 field #> geometry type: GEOMETRY #> dimension: XY #> bbox: xmin: 0 ymin: 1 xmax: 2 ymax: 4 #> epsg (SRID): 3857 #> proj4string: +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs #> a geom #> 1 1 POINT (0 1) #> 2 1 LINESTRING (1 3, 2 4)
rbind(a,b)
#> Simple feature collection with 2 features and 1 field #> geometry type: GEOMETRY #> dimension: XY #> bbox: xmin: 0 ymin: 1 xmax: 2 ymax: 4 #> epsg (SRID): 3857 #> proj4string: +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs #> a geom #> 1 1 POINT (0 1) #> 2 1 LINESTRING (1 3, 2 4)
rbind(b,c)
#> Simple feature collection with 2 features and 1 field #> geometry type: GEOMETRY #> dimension: XY #> bbox: xmin: 1 ymin: 3 xmax: 2 ymax: 4 #> epsg (SRID): 3857 #> proj4string: +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs #> a geom #> 1 1 LINESTRING (1 3, 2 4) #> 2 4 MULTILINESTRING ((1 3, 2 4))
cbind(a,b,c) # warns
#> Simple feature collection with 1 feature and 3 fields #> Active geometry column: geom #> geometry type: POINT #> dimension: XY #> bbox: xmin: 0 ymin: 1 xmax: 0 ymax: 1 #> epsg (SRID): 3857 #> proj4string: +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs #> a a.1 a.2 geom geom.1 geom.2 #> 1 1 1 4 POINT (0 1) LINESTRING (1 3, 2 4) MULTILINESTRING ((1 3, 2 4))
if (require(dplyr)) dplyr::bind_cols(a,b)
#> Loading required package: dplyr
#> #> Attaching package: ‘dplyr’
#> The following object is masked from ‘package:testthat’: #> #> matches
#> The following objects are masked from ‘package:stats’: #> #> filter, lag
#> The following objects are masked from ‘package:base’: #> #> intersect, setdiff, setequal, union
#> Simple feature collection with 1 feature and 2 fields #> Active geometry column: geom #> geometry type: POINT #> dimension: XY #> bbox: xmin: 0 ymin: 1 xmax: 0 ymax: 1 #> epsg (SRID): 3857 #> proj4string: +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs #> a geom a1 geom1 #> 1 1 POINT (0 1) 1 LINESTRING (1 3, 2 4)
c = st_sf(a=4, geomc = st_sfc(st_multilinestring(list(matrix(1:4,2)))), crs = crs) cbind(a,b,c, sf_column_name = "geomc")
#> Simple feature collection with 1 feature and 3 fields #> Active geometry column: geomc #> geometry type: MULTILINESTRING #> dimension: XY #> bbox: xmin: 1 ymin: 3 xmax: 2 ymax: 4 #> epsg (SRID): 3857 #> proj4string: +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs #> a a.1 a.2 geom geom.1 geomc #> 1 1 1 4 POINT (0 1) LINESTRING (1 3, 2 4) MULTILINESTRING ((1 3, 2 4))
df = data.frame(x=3) st_sf(data.frame(c, df))
#> Simple feature collection with 1 feature and 2 fields #> geometry type: MULTILINESTRING #> dimension: XY #> bbox: xmin: 1 ymin: 3 xmax: 2 ymax: 4 #> epsg (SRID): 3857 #> proj4string: +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs #> a x geomc #> 1 4 3 MULTILINESTRING ((1 3, 2 4))
dplyr::bind_cols(c, df)
#> Simple feature collection with 1 feature and 2 fields #> geometry type: MULTILINESTRING #> dimension: XY #> bbox: xmin: 1 ymin: 3 xmax: 2 ymax: 4 #> epsg (SRID): 3857 #> proj4string: +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs #> a geomc x #> 1 4 MULTILINESTRING ((1 3, 2 4)) 3