Preguntas etiquetadas con sp

Tengo una lista de buffers espaciales (30000 buffers) que construí con la función lapply: buff.pts <- lapply(1:nrow(pts.prj), FUN=function(l){ buff <- gBuffer(pts.prj[l,], width=1000) ## 1km return(buff) })) > head(buff.pts) [[1]] class : SpatialPolygons features : 1 extent : 307941.8, 311941.8, 4994518, 4998518 (xmin, xmax, ymin, ymax) coord. ref. : +proj=tmerc +lat_0=0 …

16 polygon  r  sp  rgeos 

Estoy migrando el código del paquete sp al paquete sf más nuevo. Mi código anterior tenía un polígono SpatialDataFrame (censimentoMap) y un SpatialPointDataFrame (indirizzi.sp) y obtuve el id. De celda de polígono ("Cell110") para cada punto con las instrucciones a continuación: points.data <- over(indirizzi.sp, censimentoMap[,"Cell110"]) En realidad, creé dos objetos …

16 r  overlay  sp  sf 

Estoy trabajando en un proyecto de epidemiología ambiental donde tengo exposiciones puntuales (~ 2,000 operaciones de cerdos industriales - OHI). Estas OHI rocían en los campos cercanos, pero las heces y el olor de las heces pueden viajar millas. Entonces, estas exposiciones puntuales obtienen amortiguadores de 3mi, y quiero saber …

12 r  point-in-polygon  sp  de-9im 

Tengo dos clases que comparten el mismo CRS (Latitute y Longitude): bolognaQuartieriMap: SpatialPolygonDataFramedatos que contienen los municipios de una ciudad. crashPoints: a que SpatialPointsDataFramecontiene datos de accidentes. Están bien trazados usando: plot(bolognaQuartieriMap) title("Crash per quartiere") plot(crashPoints, col="red",add=TRUE) Lo que necesito es obtener el número de puntos ( crashPoints) en cada …

11 r  sp 

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 

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 

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 

1 pregunta He encontrado una advertencia usando la función projectRaster () en el paquete ráster en R. A continuación se pega un ejemplo reproducible completo. Warning message: In rgdal::rawTransform(projto_int, projfrom, nrow(xy), xy[, 1], : 33940 projected point(s) not finite Mi pregunta es: ¿Es esta advertencia un problema que debo solucionar …

8 raster  r  coordinate-system  sp  rgdal 

Tengo un SpatialPointsDataFrame con el que estoy trabajando en R. Quiero agregar un nuevo vector al marco de datos que contiene, para cada punto, la distancia al otro punto más cercano en el SpatialPointsDataFrame. He mirado hasta las rodillas en el paquete spdep, así como spDistsN1 y spDists en el …

8 r  point  nearest-neighbor  sp