Skip to content
Prev 29166 / 29559 Next
Original Raw

st_intersection produce Geometry type: GEOMETRY instead of Geometry type: POLYGON

Manuel Spínola 
Dear list members,

I am trying to "crop" a polygon (grid) with a polygon, but the result is an
sf object Geometry type: GEOMETRY, instead of an sf object Geometry type:
POLYGON.

How can I obtain an sf POLYGON?

nc = st_read(system.file("shape/nc.shp", package="sf"))



nc <- st_transform(nc, 3857)



g_50000 <- st_make_grid(nc, cellsize = 50000) |> st_as_sf()



g_50000 <- g_50000[nc, ]



g_50000_d <- st_union(g_50000)



g_25000 = st_make_grid(g_50000_d, cellsize = 25000) |> st_as_sf()



g_25000 # Geometry type: POLYGON



g_25000_c <- st_intersection(g_25000, g_50000_d)



g_25000_c # Geometry type: GEOMETRY

Thread (7 messages)

Manuel Spínola st_intersection produce Geometry type: GEOMETRY instead of Geometry type: POLYGON Feb 10 Bede-Fazekas Ákos st_intersection produce Geometry type: GEOMETRY instead of Geometry type: POLYGON Feb 10 Josiah Parry st_intersection produce Geometry type: GEOMETRY instead of Geometry type: POLYGON Feb 10 Manuel Spínola st_intersection produce Geometry type: GEOMETRY instead of Geometry type: POLYGON Feb 10 Manuel Spínola st_intersection produce Geometry type: GEOMETRY instead of Geometry type: POLYGON Feb 10 Marcelino de la Cruz Rot st_intersection produce Geometry type: GEOMETRY instead of Geometry type: POLYGON Feb 10 Manuel Spínola st_intersection produce Geometry type: GEOMETRY instead of Geometry type: POLYGON Feb 10