Aumentar a velocidade de cortar, mascarar e extrair varredura por muitos polígonos no R?
Estou extraindo a área e a porcentagem de cobertura de diferentes tipos de uso da terra de uma varredura baseada em vários milhares de limites de polígonos. Descobri que a função extrair funciona muito mais rápido se eu percorrer cada polígono individual e cortar, mascarar a varredura até o tamanho …
Unindo polígonos em R
Eu estou querendo saber como unir polígonos espaciais usando o código R? Estou trabalhando com dados do censo em que determinadas áreas mudam ao longo do tempo e desejo unir os polígonos e os dados correspondentes e simplesmente informar sobre as áreas unidas. Estou mantendo uma lista de polígonos que …
Lendo a classe de recurso no geodatabase de arquivo usando R?
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 …
Trabalhando com dados PostGIS no R?
Trabalho com R quase o tempo todo e agora estou usando-o para fazer mineração de dados espaciais. Eu tenho um banco de dados PostGIS com (obviamente) dados GIS. Se eu quiser fazer análises espaciais estatísticas e mapas de plotagem, é a melhor maneira de: exportar as tabelas como shapefiles ou; …
R converte o objeto sf novamente para SpatialPolygonsDataFrame
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
Como altero a cor de preenchimento do polígono e a cor da borda dos objetos SpatialPolygons?
Eu tenho SpatialPolygons ou SpatialPolygonsDataFrames que eu gostaria de plotar. Como altero a cor usando os três sistemas de plotagem em R (gráficos de base, gráficos de treliça e ggplot2)? Dados de exemplo: library(sp) Srs1 = Polygons(list(Polygon(cbind(c(2,4,4,1,2),c(2,3,5,4,2)))), "s1") Srs2 = Polygons(list(Polygon(cbind(c(5,4,2,5),c(2,3,2,2)))), "s2") SpDF <- SpatialPolygonsDataFrame( SpatialPolygons(list(Srs1,Srs2)), data.frame( z=1:2, row.names=c("s1","s2") ) …
25 polygon r
Obtendo TopologyException: o geom de entrada 1 é inválido devido à interseção automática em R?
O erro de auto-interseção 'TopologyException: Input geom 1 is invalid' que surge de geometrias de polígonos inválidas foi amplamente discutido. No entanto, não encontrei uma solução conveniente na Web que dependa apenas da funcionalidade R. Por exemplo, eu consegui criar um objeto 'SpatialPolygons' a partir da saída da map("state", ...)boa …
24 r polygon rgeos self-intersection
Como acelerar a plotagem de polígonos em R?
Eu quero plotar as fronteiras dos países da América do Norte sobre uma imagem rasterizada representando algumas variáveis ​​e, em seguida, sobrepor contornos no topo da plotagem usando R. muito lento! Ainda não o fiz no ggplot2, mas duvido que ele se sairá melhor em termos de velocidade. Eu tenho …
24 polygon r
Usando R para calcular a área de vários polígonos em um mapa que se cruza com outro polígono sobreposto
Eu tenho um shapefile baixado do Ordnance Survey que fornece limites da ala eleitoral (divisão) para um condado do Reino Unido. Utilizei com êxito o R para carregar o shapefile e plotei vários mapas usando ggplot2como descrito nesta pergunta . Está tudo funcionando muito bem. Agora, eu gostaria de criar …
22 shapefile r
Maneira correta de vincular SpatialPolygonsDataFrames com IDs de polígono idênticos?
Qual é o idioma R adequado para vincular SPDFs juntos quando os IDs se sobrepõem? Observe que aqui (como costuma ser o caso) os IDs são basicamente sem sentido, por isso é muito irritante que eu não possa fazer o rbind ignorá-los ... library(sp) library(UScensus2000) library(UScensus2000tract) data(state) # for state …
22 polygon r
Instanciando polígono espacial sem usar um shapefile em R
Portanto, a maneira usual de ler um shapefile no R é através do pacote maptools, assim: sfdata <- readShapeSpatial("/path/to/my/shapefile.shp", proj4string=CRS("+proj=longlat")) No entanto, tenho um caso de uso em que não tenho um shapefile.shp, mas tenho uma série de coordenadas poligonais 16.484375 59.736328125,17.4951171875 55.1220703125,24.74609375 55.0341796875,22.5927734375 61.142578125,16.484375 59.736328125 e sua projeção correspondente …