如何根据 R 中的相邻多边形填充缺失值？答案

【问题标题】：How to fill missing values based on neighboring polygons in R?如何根据 R 中的相邻多边形填充缺失值？
【发布时间】：2020-07-15 22:29:54
【问题描述】：

假设如下数据集并添加一些缺失值（仅用于说明）：

library(dplyr)
library(sf)

demo(nc, ask = FALSE, verbose = FALSE)
nc$AREA[c(30, 45)] <- NA

我可以得到每个县的所有相邻多边形：

nc %>% mutate(
   INTERSECT = purrr::map(.x = geometry, .f = st_intersects, y = st_geometry(nc))
)

这给了我每行相邻县的索引列表。现在我想用相邻多边形的平均面积填充缺失的面积值。我将如何使用这些索引来获取相应行的平均值？

【问题讨论】：

nc长什么样子？

标签： r purrr sf

【解决方案1】：

index <- st_touches(nc, nc)

output <- nc %>% 
  mutate(AREA = ifelse(is.na(AREA),
                             apply(index, 1, function(i){mean(.$AREA[i])}),
                             AREA))

output$AREA[c(30, 45)]

[1] 0.1510 0.1335

检查答案：

两个多边形邻居的索引。

index[c(30, 45)]

[[1]]
[1] 13 14 29 37 48

[[2]]
[1] 44 87

手动查找区域。

mean(output$AREA[index[[30]]])

[1] 0.151

mean(output$AREA[index[[45]]])

[1] 0.1335

【讨论】：