Preguntas etiquetadas con r

Me gustaría realizar una superposición espacial para identificar el polígono dentro del cual cae un conjunto de puntos. Sin embargo, primero quiero amortiguar y disolver los polígonos, de modo que los puntos que caen en cualquier lugar dentro de los polígonos fusionados (pero no dentro de ningún agujero) se etiquetarán …

9 r  sp 

He creado dos SpatialLinesobjetos en I: . Estos objetos fueron creados de esta manera: library(sp) xy <- cbind(x,y) xy.sp = sp::SpatialPoints(xy) spl1 <- sp::SpatialLines(list(Lines(Line(xy.sp), ID="a"))) Ahora quiero concluir de alguna manera que esta es la misma línea rotada y volteada, y que la diferencia entre ellas es igual a 0 …

9 r  geometry  similarity 

Estoy tratando de encontrar algunos centroides de un archivo de formas para los municipios daneses y luego encontrar el tiempo de conducción entre ellos. Uso Rla readShapePolyfunción de maptoolscombinado con la gCentroidfunción de rgeos, y todo funciona. Sin embargo, obtengo puntos espaciales como SpatialPoints: x y 1 571860.7 6225016 Coordinate …

9 shapefile  r  coordinate-system 

¿Cómo puedo encontrar un valor de celda máximo de un ráster apilado? Rmax <- maxValue(RAD1998.all[[1]]) funciona bien pero Rmax <- maxValue(RAD1998.all[[2]]) da NA. Por supuesto que no en una trama apilada. Aquí está mi código: RAD1998 <- raster(paste(getwd(), "/1998bil/1998ASC5min_ppt_spas1214_0001_19980202_0810_UTC.asc.bil", sep = "")) list.ras <- mixedsort(list.files(paste(getwd(), "/1998bil/", sep = ""), full.names …

9 raster  r 

He calculado el área de superficie de las distribuciones de especies (fusionando polígonos de archivos de forma), pero dado que esta área puede estar compuesta de polígonos bastante distantes, me gustaría calcular alguna medida de dispersión. Lo que he hecho hasta ahora es recuperar los centroides de cada polígono, calcular …

9 polygon  r  distance  centroids 

todos, disculpen las molestias, pero estoy siendo bastante nuevo con r enfrentado una dificultad crucial: quiero crear un mapa animado de Russin con cambios en el desempleo con diferentes años, como. En la imagen puedes ver datos de un año require(sp) require(maptools) require(RColorBrewer) require(rgdal) rus<-url("http://www.filefactory.com/file/4h1hb5c1cw7r/n/RUS_adm1_RData") print(load(rus)) unempl <- read.delim2(file="C:\\unempl1.txt", header …

9 cartography  r  animation 

Estoy tratando de eliminar la intersección de Poly2 en Poly1. Debajo están los coords de los 2 polígonos. > coords1 [,1] [,2] [1,] 992.0161 7462.531 [2,] 950.1962 7413.532 [3,] 902.3632 7421.207 [4,] 706.1985 7472.378 [5,] 654.0139 7711.760 [6,] 657.5960 7726.950 [7,] 786.2667 7844.380 [8,] 789.5935 7842.925 [9,] 981.7046 7498.659 [10,] …

9 r  polygon  intersection  symmetricdifference 

Estoy interesado en crear una superficie de trama que represente la intensidad solar promedio por unidad de área para una pequeña parcela de jardín. Me imagino que las entradas en el modelo serían las siguientes: Área de parcela: por simplicidad, supongamos que la parcela del jardín es circular. Altura media …

9 arcgis-desktop  python  arcpy  raster  r 

La imagen adjunta muestra una brecha forestal con pino rojo representado como círculos y pino blanco representado como cruces. Estoy interesado en determinar si existe una asociación positiva o negativa entre las dos especies de pinos (es decir, si crecen o no en las mismas áreas). Estoy al tanto de …

9 arcgis-desktop  python  arcpy  r  spatial-statistics 

En la documentación de spVi que el coordsargumento de la SpatialPolygonsfunción solo acepta 2 columnas. ¿No es posible crear polígonos tridimensionales? Estoy haciendo kriging espacio-temporal usando el paquete gstatagregando tiempo como una tercera dimensión espacial de a SpatialPointsDataFrame. Me gustaría extender eso para bloquear kriging con polígonos como bloques. ¿Hay …

9 r  interpolation  spatial-statistics  kriging  gstat 

Necesito una cuadrícula espacial como una cuadrícula maestra para diversos mapas temáticos. ¿Cómo produzco una cuadrícula espacial a partir de un ráster que descarta todos los píxeles de NA?

9 r 

Hay 100 SpatialPointsDataFrames en mi espacio de trabajo. Los cargué de esta manera: filenames <- list.files(path="", pattern="XYhectareTravelTimes_ez+.*shp") for(i in filenames){ filepath <- file.path("/",i) assign(i, readShapePoints(filepath)) Se llaman así: XYhectareTravelTimes_ez10.*shp XYhectareTravelTimes_ez11.*shp XYhectareTravelTimes_ez12.*shp etc. ¿Cómo puedo convertirlos en rásteres que recorren el espacio de trabajo? Soy un nuevo usuario de R y …

9 r  rasterization 

Refiriéndome a mi pregunta anterior , ¿cómo puedo "geocodificar" / convertir, usando R, una imagen ráster generada a un formato que pueda mostrarse como una capa en un mapa? En resumen, mi objetivo es recrear los resultados del software comercial SpatialKey ( captura de pantalla ). Estoy usando una instancia …

9 raster  r  georeferencing 

Probablemente una pregunta muy simple. Tengo una lista de aproximadamente mil ubicaciones geográficas potenciales (lat lon), y de esas necesito seleccionar las 200 ubicaciones que están "más dispersas". Creo que serían los 200 puntos con la distancia promedio total más alta. Piensa en las tiendas de una ciudad. ¿Hay un …

9 r  spatial-statistics  business 

Este es un poste cruzado de stackoverflow . Sé muy poco sobre el software SIG y estoy haciendo todo mi mapeo en R. Disculpas por adelantado si esta es una pregunta demasiado básica. Digamos que tengo dos shapefiles de diferentes fuentes pero con diferentes atributos. Digamos, uno es para límites …

9 r  overlay