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 :...
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 <-...
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...
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
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
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...
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...
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...
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...
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...