用包含橄榄酸间隔的列表柱的毫无疑问的SF数据框

Question

我从一个名为 polygon_history 的 sf 数据框开始，其中包含以下列：

polygon_id,valid_period,geometry

其中 valid_period 是使用 lubridate 包创建的两个日期之间的间隔。

我将 valid_period 细分为更小的间隔，并通过获取 valid_period 与日期范围列表 (sectionized_periods) 的交集将它们组合在一个列表中，我有兴趣单独分析，代码如下：

sectionized_periods
[1] 2021-09-01 UTC--2021-09-07 UTC 2021-09-07 UTC--2021-10-10 UTC 2021-10-10 UTC--2021-10-24 UTC 2021-10-24 UTC--2021-11-30 UTC

polygon_history %>%
  rowwise %>%
  mutate(period_section = list(c(lubridate::intersect(valid_period,sectionized_periods)))) 

现在我的表有以下列：

polygon_id, valid_period, period_sections, geometry

其中 period_sections 是一个列表列；每行都是 [lubridate] 间隔的列表。

period_sections 中的示例行：

polygon_history$period_section[1]
[[1]]
[1] 2021-09-01 UTC--2021-09-07 UTC 2021-09-07 UTC--2021-10-10 UTC 2021-10-10 UTC--2021-10-24 UTC 2021-10-24 UTC--2021-11-30 UTC

我尝试根据 period_sections 列表列取消嵌套此表，以便为列表的每个元素（每个间隔）保留一行，但是，当我运行下面的代码时，我只得到 period_section 中的数值而不是间隔。我相信这个数值代表间隔的持续时间，但我需要有间隔。

polygon_history <- polygon_history %>%
+     rowwise %>%
+     mutate(PERIOD_SECTION = list(c(lubridate::intersect(VALID_PERIOD,sectionized_periods)))) %>%
+     unnest(PERIOD_SECTION) %>%
+     sf::st_sf() 

polygon_history$PERIOD_SECTION
An object of class "vctrs:::common_class_fallback"
 [1]  518400 2851200 1209600 3196800  518400 2851200 1209600 3196800  518400 2851200 1209600 3196800  518400 2851200 1209600 3196800  518400 2851200
[19] 1209600

感谢任何有关如何解决此问题的见解。谢谢！

Answer 1

我能够使用unchop来完成我需要的事情。

 polygon_history <- polygon_history %>%
    +     rowwise %>%
    +     mutate(PERIOD_SECTION = list(c(lubridate::intersect(VALID_PERIOD,sectionized_periods)))) %>%
    +     unchop(PERIOD_SECTION) %>%
    +     sf::st_sf()