Как преобразовать data.frame в список полигонов? - PullRequest
0 голосов
/ 02 июля 2018

У меня есть фрейм данных, подобный следующему:

long   lat   group  
1.0    1.0   1  
1.0    2.0   1  
2.0    2.0   1  
5.0    5.0   2  
5.0    6.0   2
6.0    6.0   2

Как я могу преобразовать его в список объектов многоугольника?

Спасибо!

1 Ответ

0 голосов
/ 02 июля 2018

Ваши последовательности точек не закрыты, поэтому в sf вам нужно будет немного прочесать их, чтобы превратить их в полигоны. Одним из вариантов является создание многоточечного объекта и приведение к многоугольнику:

library(tidyverse)
library(sf)

df <- data.frame(long = c(1, 1, 2, 5, 5, 6), 
                 lat = c(1, 2, 2, 5, 6, 6), 
                 group = c(1L, 1L, 1L, 2L, 2L, 2L))

df_sf <- df %>% 
    group_by(group) %>% 
    summarise(geometry = st_sfc(st_cast(st_multipoint(cbind(long, lat)), 'POLYGON'))) %>% 
    st_sf()

df_sf
#> Simple feature collection with 2 features and 1 field
#> geometry type:  POLYGON
#> dimension:      XY
#> bbox:           xmin: 1 ymin: 1 xmax: 2 ymax: 2
#> epsg (SRID):    NA
#> proj4string:    NA
#> # A tibble: 2 x 2
#>   group               geometry
#>   <int>              <POLYGON>
#> 1     1 ((1 1, 1 2, 2 2, 1 1))
#> 2     2 ((5 5, 5 6, 6 6, 5 5))

ggplot(df_sf) + geom_sf()

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...