Perguntas com a marcação «r»

Eu tenho um conjunto de dados no formato csv, como: latitude, longitude, value -45, 45, 10 .... Eu gostaria de calcular as médias de nível de município a partir desses dados, usando grama ArcGIS ou R +. Eu tenho a camada ArcGIS USA Counties, mas não consegui concluir com êxito …

9 arcgis-desktop  data  geoprocessing  r 

Quero mascarar uma varredura com base nos valores de um polígono: ou seja, defina como NA todos os valores da varredura que não são cobertos pelos polígonos. No entanto, parece que raster::mask()oculta as células (ou seja, define-as para NA), mesmo que estejam parcialmente cobertas? Isso é semelhante à pergunta: Cortar …

9 raster  r  masking 

Em Entendendo o buffer geodésico , a equipe de desenvolvimento de geoprocessamento da Esri distingue entre buffer euclidiano e geodésico. Eles concluem com "O buffer euclidiano realizado em classes de recursos projetados pode produzir buffers enganosos e tecnicamente incorretos. No entanto, o buffer geodésico sempre produzirá resultados geograficamente precisos porque …

9 coordinate-system  r  buffer  rgeos  geodesic 

Estou trabalhando com um shapefile que possui muitos polígonos. Como adiciono mais um campo chamado "area_sqkm" e calculo a área para cada polígono no shapefile?

9 shapefile  r  area 

Estou usando a biblioteca maptools do R para analisar o shapefile em uma lista de polígonos. A função readShapeSpatial shp <- readShapeSpatial("<path to my shapefile>") me fornece um objeto SpatailPolygonDataFrame . No meu exemplo, meu SpatialPolygonsDataFrame possui as seguintes colunas: > names(shp) [1] "AREA" "PERIMETER" "COMAREA_" "COMAREA_ID" "AREA_NUMBE" [6] "COMMUNITY" …

9 shapefile  r  sp 

Eu queria associar um arquivo CSV a um shapefile. Quando identifiquei um campo para facilitar a junção, por exemplo, map <- spChFIDs(map, as.character(map$ID)) voltou Erro em spChFIDs (SP, x): os comprimentos diferem Alguém pode aconselhar?

9 shapefile  r  csv  spatial-join  sp 

Graças à resposta dada nesta pergunta , fui capaz de agrupar e desenhar um mapa das divisões eleitorais em parte do Reino Unido, neste caso Pembrokeshire. O quadro de dados resultante é grande e contém dados do Ordnance Survey, pelo que seria difícil postar aqui, mas o quadro de dados …

9 shapefile  r 

Eu tenho uma pergunta sobre a melhor forma de visualizar a direção prevista da migração, considerando um conjunto de pontos com idades estimadas. Imagine que eu tenho um conjunto de coordenadas geográficas, representando comunidades. Para cada ponto, tenho uma estimativa da idade dessa comunidade. Quero desenhar uma seta no mapa …

8 coordinate-system  r  vector 

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 

Eu tenho dois rasters de resolução e extensão diferentes: > res(Elevation) [1] 0.002083333 0.002083333 > res(Ann_precip) [1] 0.008333333 0.008333333 > extent(Elevation) class : Extent xmin : -15.07722 xmax : -7.641806 ymin : 7.193611 ymax : 12.67694 > extent(Ann_precip) class : Extent xmin : -15.075 xmax : -7.641667 ymin : 7.191667 …

8 raster  r  extents  resolution 

Quero ingressar espacialmente com SpatialPolygonsDataFrame e SpatialPointsDataFrame. obtenha apenas dados de um ponto em vários polígonos, sob o código sp::over(meuse, srdf) Alguns dados de pontos em vários polígonos. Eu quero obter dados de pontos manter todos os dados de polígono. por exemplo: ponto A dentro do polígono a e b. …

8 r  sp 

Não sei como lidar com grandes dados que eu tinha dos radares. Tenho 144 imagens .tif por dia e preciso extrair um valor médio diário para cada pixel e exportá-lo em uma imagem .tif. Pensei em empilhar todas as 144 imagens, mas não sei como continuar e como exportar o …

8 r  extract 

Preciso criar um mapa do meu site de amostragem, e o diário solicitou que fizesse um mapa embutido mostrando em que local do mundo esse mapa é, esse é o código do meu mapa. library(maps) library(GISTools) map('state', fill = FALSE, xlim = c(-125, -114), ylim = c(32.2, 42.5), xlab = …

8 r  maps 

Preciso usar as estradas para comparar as posições do GPS e verificar se não há erros nos dados. Eu baixei o map.osm de um site que tinha informações da cidade de Roma e o transformou em Spatial Lines in R. src <- osmsource_file(file = "rome.osm") bb <- center_bbox(41.9, 12.3, 1000, …

8 openstreetmap  r  gps  distance  road