¿Cómo intersecar líneas y polígonos en R?

Question

Siguiendo una respuesta acerca de intersectar polígonos con líneas para dividir un polígono en unidades de polígono más pequeñas (en QGIS), quería intentar lo mismo en R. Sin embargo, ¡no puedo encontrar un método que funcione!

over() no tiene un método para polígonos que se intersectan con líneas; encontré gIntersection() de rgeos, pero falla:

require(sp)
require(rgeos)

poly <- readShapePoly("polygon.shp")
lines <- readShapeLines("lines.shp")

chopped <- gIntersection(poly, lines)

Dando como resultado:

Error en RGEOSBinTopoFunc(spgeom1, spgeom2, byid, id, "rgeos_intersection") : 
  UnsupportedOperationException: GeometryGraph::add(Geometry *): unknown geometry type: N4geos4geom18GeometryCollectionE

Actualización: Aquí hay un enlace a los archivos en cuestión.

Actualización 2: PaulG señala que funciona y después de actualizar rgeos y R me deshice del error anterior. Gracias PaulG ...

Sin embargo, gIntersection da como resultado un objeto SpatialLines independientemente de si ingreso (poly, lines) o (lines, poly) - mientras que la operación que hice en QGIS (o Arc, en los viejos tiempos) dividirá un polígono con las líneas y dará como resultado un objeto de polígono, no de línea.

Entonces, ¿cómo divido mi polígono con las líneas y obtengo polígonos como resultado?

Answer 1

He enfrentado el mismo problema varias veces. Me enfocaré en dos métodos que me funcionan

1) La función intersect del paquete raster. Devuelve un SpatialLinesDataframe que incluye el id de cada polígono.

library(raster)
newLines <- raster::intersect(SpatialLinesDataframe, SpatialPolygonDataframe)

2) Llamando a QGIS desde R con el paquete RQGIS, como se muestra aquí. Creo que este método es más rápido que otros, pero hay que tener cuidado con la versión de QGIS. Aquí está el github del creador. Además, con este paquete se pueden utilizar cualquiera de los geoalgoritmos de QGIS.

library(RQGIS)    
qgis_env <- set_env("PATH")
params <- get_args_man(alg = "qgis:intersection", 
                           qgis_env = qgis_env)
params$INPUT <- line_1
params$INPUT2 <- poly_1
out <- run_qgis(alg = "qgis:intersection",
                params = params,
                load_output = params$OUTPUT,
                qgis_env = qgis_env)

EDIT en 2018: Estaba pensando que puedes intersectar un polígono con otro para obtener polígonos más pequeños: Echa un vistazo a este enlace para el Manual del nuevo paquete sf de R (Pebesma, 2018).

• Edzer Pebesma (2018). sf: Simple Features for R. Versión del paquete 0.6-3. https://CRAN.R-project.org/package=sf
• Pebesma, E., 2018. Simple Features for R: Soporte Estandarizado para Datos Vectoriales Espaciales. The R Journal.

Answer 2

Puedes intentar usando RSAGA. No estoy muy familiarizado con ello, pero el comando sería algo así:

rsaga.geoprocessor("libshapes_polygons", "Intersección de Polígono-Línea", list(POLYGONS="polygonshape.shp",LINES="lineshape.shp",INTERSECT="result.shp"))