Als «sp» getaggte Fragen

R-Paket mit Klassen und Methoden zum Erstellen und Bearbeiten von Geodaten

3
Eine Liste von räumlichen Polygonobjekten in R zusammenführen
Ich habe eine Liste der räumlichen Puffer (30000 Puffer), die ich mit der Funktion erstellt habe 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. …
16 polygon  r  sp  rgeos 
1
Entspricht dem Punkt des sp-Pakets in Polygon / Overing mit sf
Ich migriere Code vom sp-Paket zum neueren sf-Paket. In meinem vorherigen Code hatte ich ein Polygon-SpatialDataFrame (censimentoMap) und ein SpatialPointDataFrame (indirizzi.sp) und bekam die Polygonzellen-ID ("Cell110") für jeden Punkt, der darin liegt. points.data <- over(indirizzi.sp, censimentoMap[,"Cell110"]) Eigentlich habe ich zwei sf-Objekte erstellt: shape_sf <- st_read(dsn = shape_dsn) shape_sf <- st_transform(x=shape_sf, …
16 r  overlay  sp  sf 
2
Wie funktioniert das räumliche Polygon% über% Polygon, wenn Werte in r aggregiert werden?
Ich arbeite an einem Umweltepidemiologieprojekt, bei dem ich Punktexpositionen habe (~ 2.000 industrielle Schweineoperationen - IHOs). Diese IHOs sprühen auf Felder in der Nähe, aber die Fäkalien, Wassertropfen und Gerüche können kilometerweit wandern. Diese Punktbelichtungen erhalten also 3-mi-Puffer, und ich möchte die Anzahl der IHO-Belichtungen (verschiedener Arten - Summe der …
3
Anzahl der Punkte im Polygon mit R zählen?
Ich habe zwei Klassen, die dasselbe CRS (Latitute und Longitude) teilen: bolognaQuartieriMap: a SpatialPolygonDataFrameenthält Daten eines Stadtbezirks. crashPoints: a SpatialPointsDataFrameenthält Daten von Unfällen. Sie sind gut geplottet mit: plot(bolognaQuartieriMap) title("Crash per quartiere") plot(crashPoints, col="red",add=TRUE) Was ich brauche, ist die Anzahl der Punkte ( crashPoints) in jedem Polygon zu erhalten, die …
11 r  sp 
2
Puffer und Vereinigung in R, ohne auf ein Multi-Polygon zu reduzieren
Ich möchte eine räumliche Überlagerung durchführen, um das Polygon zu identifizieren, in das eine Reihe von Punkten fällt. Ich möchte jedoch zuerst die Polygone puffern und auflösen, sodass Punkte, die irgendwo in zusammengeführte Polygone (aber nicht in Löcher) fallen, durch das Überlagerungsverfahren ähnlich gekennzeichnet werden. Leider reduziert der von mir …
9 r  sp 
1
Entfernen Sie ein bestimmtes Polygon aus SpatialPolygonDataFrame
Ich verwende die R maptools- Bibliothek, um das Shapefile in eine Liste von Polygonen zu analysieren. Die Funktion readShapeSpatial shp <- readShapeSpatial("<path to my shapefile>") gibt mir ein SpatailPolygonDataFrame- Objekt. In meinem Beispiel hat mein SpatialPolygonsDataFrame die folgenden Spalten: > names(shp) [1] "AREA" "PERIMETER" "COMAREA_" "COMAREA_ID" "AREA_NUMBE" [6] "COMMUNITY" "AREA_NUM_1" …
9 shapefile  r  sp 
4
Verbinden Sie die CSV-Datei mit dem Shapefile
Ich wollte eine CSV-Datei mit einem Shapefile verbinden. Wenn ich ein Feld identifiziert habe, um den Join zu erleichtern, z. map <- spChFIDs(map, as.character(map$ID)) es kehrte zurück Fehler in spChFIDs (SP, x): Längen unterscheiden sich Kann mir jemand raten?
9 shapefile  r  csv  spatial-join  sp 
3
Verbinden Sie räumliche Punktdaten mit mehreren Polygondaten mit R.
Ich möchte eine räumliche Verknüpfung mit SpatialPolygonsDataFrame und SpatialPointsDataFrame herstellen. Erhalten Sie nur Einpunktdaten innerhalb mehrerer Polygone unter dem Code. sp::over(meuse, srdf) Einige Punktdaten in mehrere Polygone. Ich möchte Punktedaten erhalten, alle Polygondaten behalten. Zum Beispiel: Punkt A innerhalb des Polygons a und b. jetzt kann bekommen: | point | …
8 r  sp 
2
Raster in R neu projizieren: gibt eine Warnung aus, dass projizierte Punkte nicht endlich sind?
1. Frage Ich habe eine Warnung mit der Funktion projectRaster () im Rasterpaket in R festgestellt. Ein vollständig reproduzierbares Beispiel wird unten eingefügt. Warning message: In rgdal::rawTransform(projto_int, projfrom, nrow(xy), xy[, 1], : 33940 projected point(s) not finite Meine Frage ist: Ist diese Warnung ein Problem, das ich beheben muss, wenn …
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.