gpt4 book ai didi

r - 为什么 rbind() 工作而 bind_rows() 在组合这些 sf 对象时不起作用?

转载 作者:行者123 更新时间:2023-12-04 20:26:53 25 4
gpt4 key购买 nike

我有一些想要组合的 sf 对象。出于这个问题的目的,我将它们每个都剥离为一个多边形( rtmprtmp2 )。 rbind()do.call(rbind,...) 似乎工作正常,但 bind_rows() 没有。显然我可以解决这个问题,但我错过了什么?

(问题底部的数据。)

这些工作:

library(dplyr)  ## version 0.8.3
library(sf) ## version 0.8.0
r1 <- rbind(rtmp,rtmp2) ## works
r2 <- do.call(rbind,list(rtmp,rtmp2)) ## works

identical(r1,r2) ## TRUE
all.equal(r1,r2) ## fails??

Error: Can't join on 'geometry' x 'geometry' because of incompatible types (sfc_POLYGON/sfc / sfc_POLYGON/sfc)



这失败了:
bind_rows(list(rtmp,rtmp2))

Error in .subset2(x, i) : attempt to select less than one element in get1index
In addition: Warning messages: 1: In bind_rows_(x, .id) : Vectorizing 'sfc_POLYGON' elements may not preserve their attributes 2: In bind_rows_(x, .id) : Vectorizing 'sfc_POLYGON' elements may not preserve their attributes



数据
rtmp <- structure(list(PROVCODE = "ON", geometry = structure(list(structure(list(
structure(c(7201954.77714286, 7206422.6, 7206659.96857143,
7206183.55714286, 7206456.25714286, 7206918.87142857, 7207300.78857143,
7207337.03428571, 7207488.91714286, 7207400.89428571, 7208110.25714286,
7208360.52, 7208664.28571429, 7208710.88571429, 7208819.62285714,
7209207.96, 7209308.06285714, 7209665.33428571, 7209425.42857143,
7209568.68285714, 7209522.08285714, 7209870.72285714, 7210022.60571429,
7210017.42857143, 7210321.19428571, 7210507.59714286, 7210630.14,
7210407.49142857, 7210535.21142857, 7210378.15142857, 7210921.82285714,
7210906.87428571, 7211255.58285714, 7211310.45142857, 7211675.00857143,
7212183.14285714, 7208783.07714286, 7209760.22, 7208297.23714286,
7207629.9, 7201954.77714286, 935407.037142857, 937804.568571429,
937788.902857143, 937063.325714286, 936599.088571429, 936431.897142857,
936498.582857143, 936341.52, 936410.56, 936134.408571429,
936103.34, 935840.997142857, 935891.048571429, 935715.002857143,
935972.168571429, 936172.38, 935989.428571429, 935885.871428571,
935349.102857143, 935278.337142857, 934976.297142857, 934831.317142857,
934964.217142857, 934770.911428571, 934691.517142857, 934978.022857143,
934460.24, 934237.594285714, 933921.745714286, 933302.131428571,
933074.305714286, 932540.708571429, 932384.642857143, 931995.7,
932099.337142857, 931815.214285714, 929976.52, 928202.56,
927651.188571429, 927094.268571429, 935407.037142857), .Dim = c(41L,
2L))), class = c("XY", "POLYGON", "sfg"))), class = c("sfc_POLYGON",
"sfc"), precision = 0, bbox = structure(c(xmin = 7201954.77714286,
ymin = 927094.268571429, xmax = 7212183.14285714, ymax = 937804.568571429
), class = "bbox"), crs = structure(list(epsg = NA_integer_,
proj4string = NA_character_), class = "crs"), n_empty = 0L)), row.names = c(NA,
-1L), sf_column = "geometry", agr = structure(c(PROVCODE = NA_integer_), .Label = c("constant",
"aggregate", "identity"), class = "factor"), class = c("sf",
"tbl_df", "tbl", "data.frame"))

rtmp2 <- structure(list(PROVCODE = "ON", geometry = structure(list(structure(list(
structure(c(6914891.13617828, 7022510.44138576, 7130129.74659325,
7130129.74659325, 7022510.44138576, 6914891.13617828, 6914891.13617828,
896361.565655993, 958495.600487533, 896361.565655993, 772093.495992912,
709959.461161371, 772093.495992912, 896361.565655993), .Dim = c(7L,
2L), .Dimnames = list(NULL, c("x", "y")))), class = c("XY",
"POLYGON", "sfg"))), class = c("sfc_POLYGON", "sfc"), precision = 0, bbox = structure(c(xmin = 6914891.13617828,
ymin = 709959.461161371, xmax = 7130129.74659325, ymax = 958495.600487533
), class = "bbox"), crs = structure(list(epsg = NA_integer_,
proj4string = NA_character_), class = "crs"), n_empty = 0L)), row.names = 277L, class = c("sf",
"data.frame"), sf_column = "geometry", agr = structure(c(PROVCODE = NA_integer_), .Label = c("constant",
"aggregate", "identity"), class = "factor"))

最佳答案

我以前遇到过同样的问题,我发现 rbindlist确实有助于合并列表,但您必须将其转换回 sf对象使用 st_as_sf()
这对我有用:

p <- data.table::rbindlist(list(rtmp,rtmp2),
use.names = TRUE,
fill = TRUE,
idcol = NULL)
st_as_sf(p)

Simple feature collection with 2 features and 1 field
geometry type: POLYGON
dimension: XY
bbox: xmin: 7201955 ymin: 927094.3 xmax: 7212183 ymax: 937804.6
epsg (SRID): NA
proj4string: NA
PROVCODE geometry
1 ON POLYGON ((7201955 935407, 7...
2 ON POLYGON ((6914891 896361.6,...

关于r - 为什么 rbind() 工作而 bind_rows() 在组合这些 sf 对象时不起作用?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/58649456/

25 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com