一連のジオメトリの属性を要約し、重複する値を合計したいと思います。
library(devtools)
install_github("r-spatial/sf")
library(sf)
m = rbind(c(0,0), c(1,0), c(1,1), c(0,1), c(0,0))
p = st_polygon(list(m))
n = rbind(c(0.5,0.5), c(1.5,0.5), c(1.5,1.5), c(0.5,1.5), c(0.5,0.5))
q = st_polygon(list(n))
s = st_sfc(list(p, q))
sf = st_sf(s, att=c(1,1))
d = st_intersection(sf)
d$id <- 1:nrow(d)
plot(d['att'])
plot(st_centroid(d['att']), add = TRUE, col = 'red')
> d
#Simple feature collection with 3 features and 4 fields
#geometry type: POLYGON
#dimension: XY
#bbox: xmin: 0 ymin: 0 xmax: 1.5 ymax: 1.5
#epsg (SRID): NA
#proj4string: NA
# att geometry n.overlaps origins id
#1 1 POLYGON ((1 0.5, 1 0, 0 0, ... 1 1 1
#2 1 POLYGON ((0.5 1, 1 1, 1 0.5... 2 1, 2 2
#3 1 POLYGON ((0.5 1, 0.5 1.5, 1... 1 2 3
上記の最小限の例では、d$att を合計して、ジオメトリ #2 (オーバーラップに対応するもの) の att=2 を取得したいと考えています。
どんな助けでも大歓迎です。