Perguntas com a marcação «rgdal»

Eu tenho um recurso contido em um geodatabase maior que 2 GB como um shapefile exportado. Eu preciso executar uma função de extração em R para atribuir os polígonos com dados de um arquivo raster. Exportar o recurso como uma tabela não é uma solução. Como posso ler as classes …

27 r  features  file-geodatabase  feature-class  rgdal 

Estou tentando ler uma tabela diretamente de um geodatabase de arquivo ESRI para R. Um arquivo de dados de exemplo pode ser baixado aqui . O banco de dados contém uma classe de recurso de ponto (Zone9_2014_01_Broadcast) e duas tabelas vinculadas (Zone9_2014_01_Vessel e Zone9_2014_01_Voyage). Você pode ler o shapefile no …

21 r  file-geodatabase  table  rgdal 

Zillow tem um conjunto de arquivos de forma para diferentes bairros das principais cidades dos EUA. Eu queria verificar se certos edifícios estavam presentes em determinados bairros usando R: library(rgeos) library(sp) library(rgdal) df <- data.frame(Latitude =c(47.591351, 47.62212,47.595152), Longitude = c(-122.332271,-122.353985,-122.331639), names = c("Safeco Field", "Key Arena", "Century Link")) coordinates(df) <- …

19 r  rgdal  zillow 

Eu tenho dois polígonos. Um contém campos (X, Y, Z) e o outro contém tipos de solo (A, B, C, D). Quero saber qual área de cada campo contém qual tipo de solo. Eu tentei o seguinte: library(rgdal) library(rgeos) Field<-readOGR("./","Field") Soil<-readOGR("./","Soil") Results<-gIntersects(Soil,Field,byid=TRUE) rownames(Results)<-Field@data$FieldName colnames(Results)<-Soil@data$SoilType > Results A B C D …

19 shapefile  polygon  r  intersection  rgdal 

Eu quero criar dois polígonos. Um deles é da extensão retangular da minha varredura. Eu sei que posso obter a extensão da varredura usando r <- raster("band5.tif") e <- extent(r) plot(e)me dá o seguinte: Como posso criar UM polígono do limite da varredura, como mostrado abaixo? É assim que minha …

13 raster  r  rgdal 

Estou usando gSimplify(pacote rgeos) para simplificar as geometrias de um shapefile. A função funciona bem, mas agora não consigo gravar a saída em um novo shapefile. Eu tentei algumas maneiras: writeOGR(simplyshape, file, driver="ESRI Shapefile", layer='test') eu recebo obj deve ser um SpatialPointsDataFrame, SpatialLinesDataFrame ou SpatialPolygonsDataFrame e com: writePolyShape(simplyshape, file) Eu …

12 shapefile  r  simplify  rgdal 

Eu tenho o NE lat lnge SW lat lng. Meu objetivo é simples, estou usando o RStudio e quero criar uma caixa delimitadora a partir dos dois lats e longos acima e, em seguida, gerar Polygon e quero adicionar pontos aleatórios no polígono quadrado. Eu tenho NELat / Lng e …

11 r  rgdal 

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 tenho um arquivo de formato ESRI muito grande (> 8000 polígonos, o arquivo .shp é> 32 MB), que eu quero ler no R no Windows 7. Estou familiarizado com rgdal e readOGR (). No entanto, com esse tamanho de arquivo de forma, a leitura de todo o arquivo de …

9 shapefile  r  rgdal 

Estou convertendo vetor para varredura em R. No entanto, o processo foi muito longo. Existe a possibilidade de colocar o script no processamento multithread ou GPU para fazer isso mais rapidamente? Meu script para vetor rasterizado. r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do raster res(r.raster) = 10 …

9 raster  r  rasterization  rgdal  parallel 

1. Pergunta Encontrei um aviso usando a função projectRaster () no pacote raster em R. Um exemplo completo e reproduzível é colado abaixo. Warning message: In rgdal::rawTransform(projto_int, projfrom, nrow(xy), xy[, 1], : 33940 projected point(s) not finite Minha pergunta é: Esse aviso é um problema que preciso corrigir se estiver …

8 raster  r  coordinate-system  sp  rgdal