Perguntas com a marcação «r»

Eu usei a integração do R com o QGIS e achei uma ferramenta muito útil. Tendo examinado a extensão da caixa de ferramentas R para 9 e a extensão para 10 , parece que ambos oferecem apenas agrupamento de pontos para demonstrar a prova de conceito. Outros métodos de análise …

22 arcgis-desktop  r  spatstat 

Eu criei 2 polígonos no QGIS. Utilizando-os em R, os polígonos se tornam automaticamente SpatialPolygonsDataFrame (SPDF). Gostaria de mesclá-los em um único SPDF (como é super fácil no ArcGis usando a ferramenta Mesclar ). Estou certo de que deve haver uma maneira simples de concluir isso em R, mas não …

22 r  spatial-analyst  merge  maptools 

Esta pergunta foi migrada do Stack Overflow porque pode ser respondida no Geographic Information Systems Stack Exchange. Migrou há 4 anos . Estou tentando executar uma junção espacial entre dados de ponto e dados de polígono. Eu tenho dados que indicam as coordenadas espaciais de um evento no meu arquivo …

21 r  spatial-join  point-in-polygon 

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 

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 

Dada uma amostra grande (~ 1 milhão) de pontos desigualmente distribuídos - é possível gerar uma grade irregular (em tamanho, mas também pode ter uma forma irregular se for possível?) Que conterá uma quantidade mínima especificada de n pontos? É menos importante para mim se as 'células' geradas dessa grade …

20 r  clustering  quadtree 

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 traçando algumas geolocalizações (lon, lat)de tweets coletados usando a API do Stream com 7mi. raio. A figura abaixo mostra dois mapas dos mesmos dados, mas usando alfa / transparência diferentes para cada observação. Olhando para o mapa à esquerda, notei que as geolocalizações seguem principalmente uma espécie de "grade" …

19 r  geolocation  twitter 

coast<-readShapeSpatial("coastline.shp") landc<-readShapeSpatial("landcover.shp") ro<-readShapeSpatial("roads.shp") bc<-gBuffer(ro,width=100) landc$ratings=1 landc$ratings[landc$LANDUSE_ID==4]=0 Acima, eu pego qualquer categoria que tenha 4 e, na nova coluna, coloque-a como 0. Nesse ponto, quero a coluna também nomeada ratingspara bc, onde será necessário 0 se estiver dentro do buffer e 1 se estiver fora. O problema é que o bcé …

19 r  cartography 

Eu tenho um conjunto de imagens rasterizadas que representam um mês específico ao longo dos anos e quero fazer uma linha do tempo dos valores médios de uma área usando um shapefile. Como extraio os valores das rasters e os importo em R de maneira que eu possa usá-los?

19 r  spatial-statistics 

O I de Moran , uma medida de autocorrelação espacial, não é uma estatística particularmente robusta (pode ser sensível a distribuições distorcidas dos atributos de dados espaciais). Quais são algumas técnicas mais robustas para medir a autocorrelação espacial? Estou particularmente interessado em soluções prontamente disponíveis / implementáveis ​​em uma linguagem …

19 algorithm  r  spatial-statistics 

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 

Estou tentando realizar uma união em um campo comum após mesclar dois shapefiles adjacentes. Os arquivos de forma acabam com pelo menos um pedaço fino de espaço entre eles. Quando tento uma união, recebo o seguinte erro de orfão órfão: Erro no createPolygonsComment (p): rgeos_PolyCreateComment: furo órfão, não é possível …

18 polygon  r  merge 

Meu objetivo é modificar um arquivo de forma existente mesclando certos polígonos. Depois de importar o shapefile e usar o comando UnionSpatialPolygons, recebo o contorno do polígono que desejo. No entanto, agora este é um objeto SpatialPolygons e não um SpatialPolygonsDataFrame, portanto, não consigo exportá-lo para um shapefile usando o …

18 shapefile  cartography  r  ogr