J'ai une liste de tampons spatiaux (30000 tampons) que j'ai construits avec la fonction 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 :...
Je migre le code du package sp vers le package sf plus récent. Mon code précédent, j'avais un polygone SpatialDataFrame (censimentoMap) et un SpatialPointDataFrame (indirizzi.sp) et j'ai obtenu l'ID de cellule de polygone ("Cell110") pour chaque point se trouvant à l'intérieur avec les instructions...
Je travaille sur un projet d'épidémiologie environnementale où j'ai des expositions ponctuelles (~ 2 000 porcs industriels - OHI). Ces OHI pulvérisent sur les champs voisins, mais les gouttelettes d'eau et les odeurs des matières fécales peuvent parcourir des kilomètres. Donc, ces expositions...
J'ai deux classes partageant le même CRS (Latitute et Longitude): bolognaQuartieriMap: un SpatialPolygonDataFramecontenant des données d'un arrondissement de la ville. crashPoints: un SpatialPointsDataFramecontenant des données d'accidents. Ils sont bien tracés en utilisant:
Je voulais joindre un fichier CSV à un fichier de formes. Lorsque j'ai identifié un champ pour faciliter la jointure, par exemple, map <- spChFIDs(map, as.character(map$ID)) il est revenu Erreur dans les spChFID (SP, x): les longueurs diffèrent Quelqu'un peut-il
Je voudrais effectuer une superposition spatiale pour identifier le polygone dans lequel se situe un ensemble de points. Cependant, je veux d'abord tamponner et dissoudre les polygones, de sorte que les points situés n'importe où dans les polygones fusionnés (mais pas dans les trous) seront...
J'utilise la bibliothèque R maptools pour analyser le fichier de formes dans une liste de polygones. La fonction readShapeSpatial shp <- readShapeSpatial("<path to my shapefile>") me donne un objet SpatailPolygonDataFrame . Dans mon exemple, mon SpatialPolygonsDataFrame a les colonnes...
J'ai un SpatialPointsDataFrame avec lequel je travaille dans R. Je veux ajouter un nouveau vecteur au dataframe qui contient, pour chaque point, la distance à l'autre point le plus proche dans le SpatialPointsDataFrame. J'ai regardé knearneigh dans le package spdep ainsi que spDistsN1 et spDists...
1. Question J'ai rencontré un avertissement en utilisant la fonction projectRaster () dans le package raster en R. Un exemple reproductible complet est collé ci-dessous. Warning message: In rgdal::rawTransform(projto_int, projfrom, nrow(xy), xy[, 1], : 33940 projected point(s) not finite Ma...
Je veux joindre spatialement avec SpatialPolygonsDataFrame et SpatialPointsDataFrame. obtenir une seule donnée ponctuelle dans plusieurs polygones, sous le code. sp::over(meuse, srdf) Certains points de données dans plusieurs polygones. Je veux que les données de points conservent toutes les...