Preguntas etiquetadas con sf

¿Existe una forma nativa de SF (es decir, "correcta") de extraer todo excepto la columna de geometría de un objeto de características simples? Esto funciona df <- dplyr::select(as.data.frame(sf), -geometry) pero el select( , -geometry)paso se siente innecesario. Además, no elimina los atributos de geometría.

44 r  sf 

¿Cómo convierto un objeto sf de nuevo a un SpatialPolygonsDataFrame? Estoy ejecutando R 3.4.0 y sf versión 0.4.1. La documentación de SF menciona una as_Spatialfunción, pero esta no está disponible en mi instalación (por alguna razón).

27 r  sf 

¿Existe una función para recortar el objeto de mapa sf, similar a la maptools::pruneMap(lines, xlim= c(4, 10), ylim= c(10, 15))utilizada para SpatialPolygon o SpatialLine? Estoy considerando st_intersection()pero puede haber una manera adecuada.

20 r  clip  sf 

¿Cómo puedo convertir los puntos lon-lat en características simples (sfg) y luego ponerlos en una colección de características simples (sfc)? Aquí hay un MWE que no funciona, pero es el más cercano al que he llegado. library(data.table) library(sf) # The DT data.table is the data I have (but 10,000s of …

20 r  sf 

Estoy migrando el código del paquete sp al paquete sf más nuevo. Mi código anterior tenía un polígono SpatialDataFrame (censimentoMap) y un SpatialPointDataFrame (indirizzi.sp) y obtuve el id. De celda de polígono ("Cell110") para cada punto con las instrucciones a continuación: points.data <- over(indirizzi.sp, censimentoMap[,"Cell110"]) En realidad, creé dos objetos …

16 r  overlay  sp  sf 

¿Cómo simplifico un sfpolígono sin introducir huecos y astillas? Con un shapefile, por ejemplo, usaría rmapshaper::ms_simplify(): library("pryr") library("rgdal") library("rmapshaper") download.file("https://borders.ukdataservice.ac.uk/ukborders/easy_download/prebuilt/shape/England_gor_2011.zip", destfile = "regions.zip") unzip("regions.zip") regions <- readOGR(".", "england_gor_2011") object_size(regions) # ~13MB regions <- ms_simplify(regions) object_size(regions) # < 1MB He intentado lo sf::st_cast()que, desde las páginas de manual, dice: Convertir geometría …

14 r  simplify  sf 

Usando el paquete R sf, ¿cómo se combinan los sfcobjetos? Por ejemplo, dado el siguiente código, ¿cómo se puede crear un solo sfcobjeto sfc12que incluya las geometrías de ambos sfc1y sfc2? ( length(sfc12)debe ser 2.) library(sf) pt1 = st_point(c(0,1)) pt2 = st_point(c(1,1)) sfc1 = st_sfc(pt1) # An sfc object sfc2 …

12 r  sf 

Estoy intentando escribir varias capas a la misma GeoPackage en R, pero me da un error, Creation of output file failed. He intentado buscar documentación de lectura y escritura en archivos .gpkg con RGDAL, por ejemplo, para averiguar si writeOGR()realmente admite varias capas, con poco éxito. ¿Es esto posible, si …

11 gdal  r  rgdal  geopackage  sf 

Quiero crear varias líneas a partir de puntos dados como sf objects. Si tengo varios puntos como library(sf) pts <- st_multipoint(matrix(c(10, 10, 15, 20, 30, 30), nrow = 3, byrow = TRUE), dim = "XY") y estoy usando st_castpara crear líneas de ellos lines <- st_cast(pts, "MULTILINESTRING") Siempre obtendré uno …

8 r  sf 

Dada una línea y un conjunto de puntos, no puedo entender cómo usar sfpara identificar de qué lado de la línea cae cada punto. Sigue un pequeño ejemplo reproducible, adaptado de una pregunta diferente # Load Libraries ---------------------------------------------------------- library('sf') # Test data --------------------------------------------------------------- points.df <- data.frame( 'x' = c(-53.50000, -54.15489, …

8 r  sf