Perguntas com a marcação «r»

Eu tenho um modelo de elevação plotado em R r <- raster("example.dem") plot(r, col = topo.colors(20)) Alguns dos valores estão abaixo do nível do mar (0) e eu gostaria que esses fossem exibidos em vermelho. Como atribuir intervalos específicos a cores específicas na plotagem ()?

30 raster  r  legend  color 

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 …

29 raster  r  performance  masking  clip 

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 …

29 polygon  r  spatial-statistics  attribute-joins  census 

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 

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; …

27 postgis  postgresql  r 

Ao pesquisar na web, as soluções para encontrar centróides de polígonos aparecem com bastante frequência. O que me interessa é encontrar um centróide de um conjunto de pontos. Uma média ponderada das sortes. Eu apreciaria se alguém pudesse fornecer alguns ponteiros, pseudo-código (ou melhor ainda, um pacote R que já …

27 r  centroids 

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 

Estou tendo problemas para localizar qualquer discussão ou comparação comparativa de diferentes formatos de arquivo raster (por exemplo, para uso na análise de dados em R). Alguém tem alguma idéia de por que determinados formatos podem ser mais rápidos ou mais lentos? Ou as diferenças devem ser mínimas? Especificamente, estou …

25 raster  geotiff-tiff  r  file-formats  netcdf 

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 

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 

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 

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 

Eu tenho uma imagem rasterizada como esta Quero excluir o ponto de ruído isolado e suavizar a borda (digamos os pontos únicos brancos e quero que as bordas da entidade pareçam suaves). Como posso fazer isso no ArcGIS ou no R? O tamanho da célula é 30 * 30.

22 arcgis-desktop  raster  r  generalization  smoothing 

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 

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 …

22 shapefile  polygon  r