Perguntas com a marcação «sf»

Existe uma maneira nativa sf (isto é, "correta") de extrair tudo, exceto a coluna geometria, de um objeto de recursos simples? Isso funciona df <- dplyr::select(as.data.frame(sf), -geometry) mas o select( , -geometry)passo parece desnecessário. Além disso, ele não remove os atributos de geometria.

43 r  sf 

Como converter um objeto sf novamente em um SpatialPolygonsDataFrame? Estou executando o R 3.4.0 e a versão 0.4.1 do sf. A documentação do sf menciona uma as_Spatialfunção, mas isso não está disponível na minha instalação (por algum motivo).

27 r  sf 

Existe uma função para cortar o objeto sf map, semelhante ao maptools::pruneMap(lines, xlim= c(4, 10), ylim= c(10, 15))usado para SpatialPolygon ou SpatialLine? Estou considerando, st_intersection()mas pode haver maneira adequada.

20 r  clip  sf 

Como converter pontos long-lat em recursos simples (sfg) e depois colocá-los em uma coleção de recursos simples (sfc)? Aqui está um MWE que não funciona, mas é o mais próximo que eu cheguei. library(data.table) library(sf) # The DT data.table is the data I have (but 10,000s of rows, each row …

20 r  sf 

Estou migrando código do pacote sp para o pacote sf mais recente. Meu código anterior tinha um polígono SpatialDataFrame (censimentoMap) e um SpatialPointDataFrame (indirizzi.sp) e obtive o ID da célula do polígono ("Cell110") para cada ponto dentro da instrução abaixo: points.data <- over(indirizzi.sp, censimentoMap[,"Cell110"]) Na verdade, eu criei dois objetos …

16 r  overlay  sp  sf 

Como simplifico um sfpolígono sem introduzir espaços e lascas? Com um shapefile, por exemplo, eu usaria 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 Eu tentei o sf::st_cast()que, nas páginas de manual, afirma: Transmitir geometria …

14 r  simplify  sf 

Usando o pacote R sf, como combinar sfcobjetos? Por exemplo, dado o código a seguir, como alguém criaria um único sfcobjeto sfc12que inclui as geometrias de ambos sfc1e sfc2? ( length(sfc12)deve ser 2.) library(sf) pt1 = st_point(c(0,1)) pt2 = st_point(c(1,1)) sfc1 = st_sfc(pt1) # An sfc object sfc2 = st_sfc(pt2) …

12 r  sf 

Estou tentando gravar várias camadas no mesmo GeoPackage no R, mas recebo um erro Creation of output file failed,. Eu tentei procurar documentação de leitura e gravação em arquivos .gpkg com RGDAL, por exemplo, para descobrir se writeOGR()realmente suporta várias camadas, com pouco sucesso. Isso é possível, se sim, como …

11 gdal  r  rgdal  geopackage  sf 

Eu quero criar várias linhas a partir de pontos determinados como sf objects. Se eu tiver vários pontos como library(sf) pts <- st_multipoint(matrix(c(10, 10, 15, 20, 30, 30), nrow = 3, byrow = TRUE), dim = "XY") e eu estou usando st_castpara criar linhas deles lines <- st_cast(pts, "MULTILINESTRING") Sempre …

8 r  sf 

Dada uma linha e um conjunto de pontos, não consigo descobrir como usar sfpara identificar qual lado da linha cada ponto cai. Um pequeno exemplo reproduzível segue, adaptado de uma pergunta diferente # Load Libraries ---------------------------------------------------------- library('sf') # Test data --------------------------------------------------------------- points.df <- data.frame( 'x' = c(-53.50000, -54.15489, -54.48560, -52.00000, …

8 r  sf