Systèmes d'information géographique r

3

Créer un fichier de formes au niveau du comté à partir d'un ensemble de points à un espacement de 30 x 30 km à l'aide d'ArcGIS for Desktop?

J'ai un ensemble de données sous forme csv, comme: latitude, longitude, value -45, 45, 10 .... Je voudrais calculer les moyennes au niveau du comté à partir de ces données, en utilisant ArcGIS ou R + grass. J'ai la couche ArcGIS USA Counties, mais je n'ai pas réussi à terminer …

9 arcgis-desktop data geoprocessing r

1

Raster raster: masquer un raster par polygone supprime également les cellules partiellement couvertes?

Je veux masquer un raster en fonction des valeurs d'un polygone: c'est-à-dire définir à NA toutes les valeurs du raster qui ne sont pas couvertes par les polygones. Cependant, il semble que raster::mask()les cellules soient masquées (c'est-à-dire qu'elles soient définies sur NA) même si elles sont partiellement couvertes? Ceci est …

9 raster r masking

1

Tamponnage euclidien et géodésique en R

Pour comprendre la mise en mémoire tampon géodésique , l'équipe de développement du géotraitement d'Esri fait la distinction entre la mise en mémoire tampon euclidienne et géodésique. Ils concluent avec "la mise en mémoire tampon euclidienne effectuée sur les classes d'entités projetées peut produire des tampons trompeurs et techniquement incorrects. …

9 coordinate-system r buffer rgeos geodesic

1

Calculer l'aire de chaque polygone dans R

Je travaille avec un fichier de formes contenant de nombreux polygones. Comment puis-je ajouter un champ supplémentaire nommé "area_sqkm" et calculer la zone pour chaque polygone dans le fichier de formes?

9 shapefile r area

1

Supprimer certains polygones de SpatialPolygonDataFrame

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 suivantes: > names(shp) [1] "AREA" "PERIMETER" "COMAREA_" "COMAREA_ID" "AREA_NUMBE" [6] "COMMUNITY" …

9 shapefile r sp

4

Joignez le fichier CSV au fichier de formes

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 conseiller?

9 shapefile r csv spatial-join sp

1

Visualiser la direction de migration la plus probable étant donné les points avec les âges en utilisant R?

J'ai une question sur la meilleure façon de visualiser la direction de migration prévue compte tenu d'un ensemble de points avec des âges estimés. Imaginez que j'ai un ensemble de coordonnées géographiques représentant les communautés. Pour chaque point, j'ai une estimation de l'âge de cette communauté. Je veux dessiner une …

8 coordinate-system r vector

2

Le package R: sf pointe vers plusieurs lignes avec st_cast

Je veux créer plusieurs lignes à partir de points donnés en tant que sf objects. Si j'ai un certain nombre de points library(sf) pts <- st_multipoint(matrix(c(10, 10, 15, 20, 30, 30), nrow = 3, byrow = TRUE), dim = "XY") et j'utilise st_castpour en créer des lignes lines <- st_cast(pts, …

8 r sf

2

Sélection d'entités «au-dessus» ou «en dessous» d'une ligne à l'aide de R

Étant donné une ligne et un ensemble de points, je ne peux pas comprendre comment utiliser sfpour identifier de quel côté de la ligne chaque point tombe. Un petit exemple reproductible suit, adapté d'une question différente # Load Libraries ---------------------------------------------------------- library('sf') # Test data --------------------------------------------------------------- points.df <- data.frame( 'x' = …

8 r sf

1

Utiliser le rééchantillonnage par rapport à l'agrégation et étendre dans R pour avoir des rasters de résolution et d'étendue correspondantes

J'ai deux rasters de résolution et d'étendue différentes: > 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

3

Joignez des données de points spatiaux avec plusieurs données de polygones à l'aide de R

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 données de polygone. par exemple: point A dans les polygones a et …

8 r sp

1

Extraction de la moyenne de plusieurs couches raster à l'aide de R?

Je ne sais pas comment gérer les mégadonnées que j'ai reçues des radars. J'ai 144 images .tif par jour et je dois extraire une valeur quotidienne moyenne, pour chaque pixel, et l'exporter dans une image .tif. J'ai pensé empiler toutes les 144 images, mais je ne sais pas comment continuer …

8 r extract

4

Créer une carte en encart dans R

J'ai besoin de créer une carte de mon site d'échantillonnage, et le journal a demandé de faire une carte en encart montrant où dans le monde cette carte est, c'est mon code de carte. library(maps) library(GISTools) map('state', fill = FALSE, xlim = c(-125, -114), ylim = c(32.2, 42.5), xlab = …

8 r maps

1

Comment puis-je utiliser les routes à partir des données OSM dans R?

Je dois utiliser les routes pour comparer les positions GPS pour vérifier qu'il n'y a pas d'erreurs dans les données. J'ai téléchargé le map.osm à partir d'un site contenant des informations sur la ville de Rome et je l'ai transformé en lignes spatiales en R. src <- osmsource_file(file = "rome.osm") …

8 openstreetmap r gps distance road

2

Reprogrammer le raster en R: avertit que le ou les points projetés ne sont pas finis?

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 question est la suivante: cet avertissement est-il un problème que …

8 raster r coordinate-system sp rgdal

Questions marquées «r»