Я хочу выполнить то, что я бы назвал объединением 2 объектов sf.
У меня есть следующий код:
east.west.sf <- st_sfc(st_polygon(list(cbind(c(1,2,2,1,1),c(0,0,2,2,0)))),
st_polygon(list(cbind(c(0,1,1,0,0),c(0,0,2,2,0)))))
east.west.df <- data.frame(var1=c("east", "west"), var2=c(1,2))
east.west <- st_sf(east.west.df, geom=east.west.sf)
north.south.sf <- st_sfc(st_polygon(list(cbind(c(0,2,2,0,0),c(1,1,2,2,1)))),
st_polygon(list(cbind(c(0,2,2,0,0),c(0,0,1,1,0)))))
north.south.df <- data.frame(var3=c("north", "south"), var4=c(FALSE, TRUE))
north.south <- st_sf(north.south.df, geom=north.south.sf)
map.overlay <- st_union(east.west, north.south)
st_area(map.overlay)
Я ожидал, что все четыре области будут по 1 единице. Почему они 3 единицы?
Аналогичная операция union () с sp-объектами дает мне то, что я ожидаю.
Как я могу получить четыре объекта по 1 единице с моими двумя начальными объектами sf?