总结sf:：st_intersection（）中几何图形重叠的属性

我想总结一组几何图形的属性，总结它们重叠的值

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$att求和，得到几何体#2（对应于重叠的那个）的att=2

---

任何帮助都将不胜感激。

您可以使用

原点

列表列检索“源自”相交的多边形，然后在其

附件

列上求和。类似的东西可以工作（至少在这个非常简单的用例中…）：

---

@我觉得它似乎工作得很好。我只是想知道如何避免一个循环，但我会接受它！谢谢当然，它可以通过*apply/purrr来完成，尽管有时我想知道为什么要这么做（除非有性能问题），因为For循环语法非常容易阅读和理解（至少对我来说…；-）

> 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

for (int in seq_along(d$id)) {
     d$att[int] = sum(sf[d$origins[[int]], ]$att, na.rm = TRUE)
}
> 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   2 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