Preguntas etiquetadas con r

Tengo un conjunto de datos en formato csv, como: latitude, longitude, value -45, 45, 10 .... Me gustaría calcular las medias a nivel de condado a partir de estos datos, utilizando ArcGIS o R + grass. Tengo la capa de los condados de ArcGIS USA, pero no he podido completar …

9 arcgis-desktop  data  geoprocessing  r 

Quiero enmascarar un ráster basado en los valores de un polígono: es decir, establecer en NA todos los valores del ráster que no estén cubiertos por los polígonos. Sin embargo, parece que raster::mask()enmascara las células (es decir, las establece en NA) incluso si están parcialmente cubiertas. Esto es similar a …

9 raster  r  masking 

Al comprender el almacenamiento en búfer geodésico , el equipo de desarrollo de geoprocesamiento de Esri distingue entre el almacenamiento en búfer geodésico y euclidiano. Concluyen con "El almacenamiento en búfer euclidiano realizado en clases de entidad proyectadas puede producir amortiguadores engañosos y técnicamente incorrectos. Sin embargo, el almacenamiento en …

9 coordinate-system  r  buffer  rgeos  geodesic 

Estoy trabajando con un archivo de formas que tiene muchos polígonos. ¿Cómo agrego un campo más llamado "area_sqkm" y calculo el área para cada polígono en el shapefile?

9 shapefile  r  area 

Estoy usando la biblioteca R maptools para analizar el archivo de forma en una lista de polígonos. La función readShapeSpatial shp <- readShapeSpatial("<path to my shapefile>") me da un objeto SpatailPolygonDataFrame . En mi ejemplo, mi SpatialPolygonsDataFrame tiene las siguientes columnas: > names(shp) [1] "AREA" "PERIMETER" "COMAREA_" "COMAREA_ID" "AREA_NUMBE" [6] …

9 shapefile  r  sp 

Quería unir un archivo CSV a un shapefile. Cuando identifiqué un campo para facilitar la unión, por ejemplo, map <- spChFIDs(map, as.character(map$ID)) volvió Error en spChFID (SP, x): las longitudes difieren ¿Alguien puede aconsejar?

9 shapefile  r  csv  spatial-join  sp 

Gracias a la respuesta dada en esta pregunta , he podido subconjuntar y dibujar un mapa de divisiones electorales en parte del Reino Unido, en este caso Pembrokeshire. El marco de datos resultante es grande y contiene datos de Ordnance Survey, por lo que sería difícil publicarlo aquí, pero el …

9 shapefile  r 

Tengo una pregunta sobre la mejor manera de visualizar la dirección prevista de la migración dado un conjunto de puntos con edades estimadas. Imagine que tengo un conjunto de coordenadas geográficas que representan comunidades. Para cada punto, tengo una estimación de la edad de esa comunidad. Quiero dibujar una flecha …

8 coordinate-system  r  vector 

Quiero crear varias líneas a partir de puntos dados como sf objects. Si tengo varios puntos como library(sf) pts <- st_multipoint(matrix(c(10, 10, 15, 20, 30, 30), nrow = 3, byrow = TRUE), dim = "XY") y estoy usando st_castpara crear líneas de ellos lines <- st_cast(pts, "MULTILINESTRING") Siempre obtendré uno …

8 r  sf 

Dada una línea y un conjunto de puntos, no puedo entender cómo usar sfpara identificar de qué lado de la línea cae cada punto. Sigue un pequeño ejemplo reproducible, adaptado de una pregunta diferente # Load Libraries ---------------------------------------------------------- library('sf') # Test data --------------------------------------------------------------- points.df <- data.frame( 'x' = c(-53.50000, -54.15489, …

8 r  sf 

Tengo dos rásteres de diferente resolución y extensión: > 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 ymax …

8 raster  r  extents  resolution 

Quiero unirme espacialmente con SpatialPolygonsDataFrame y SpatialPointsDataFrame. obtener solo un punto de datos dentro de múltiples polígonos, bajo el código. sp::over(meuse, srdf) Algunos puntos de datos en múltiples polígonos. Quiero obtener datos de puntos y mantener todos los datos de polígonos. por ejemplo: punto A dentro del polígono ay b. …

8 r  sp 

No sé cómo lidiar con los grandes datos que tenía de los radares. Tengo 144 imágenes .tif por día y necesito extraer un valor medio diario, para cada píxel, y exportarlo en una imagen .tif. Pensé en apilar las 144 imágenes, pero no sé cómo continuar y cómo exportar el …

8 r  extract 

Necesito crear un mapa de mi sitio de muestreo, y la revista solicitó hacer un mapa insertado que muestre en qué parte del mundo está este mapa, este es mi código de mapa. library(maps) library(GISTools) map('state', fill = FALSE, xlim = c(-125, -114), ylim = c(32.2, 42.5), xlab = "lon", …

8 r  maps 

Necesito usar las carreteras para comparar las posiciones de GPS para verificar que no haya errores en los datos. Descargué el map.osm de un sitio que tiene información de la ciudad de Roma y lo transformé en Spatial Lines en R. src <- osmsource_file(file = "rome.osm") bb <- center_bbox(41.9, 12.3, …

8 openstreetmap  r  gps  distance  road