Ich habe die Integration von R mit QGIS verwendet und fand es ein sehr nützliches Werkzeug. Nach einem Blick auf die R-Toolbox- Erweiterung für 9 und die Erweiterung für 10 bieten beide nur Punkt-Clustering an, um den Proof-of-Concept zu demonstrieren. Können andere in R verfügbare Analysemethoden in ArcGIS Desktop über …
Ich habe 2 Polygone in QGIS erstellt. Wenn Sie sie in R verwenden, werden Polygone automatisch zu SpatialPolygonsDataFrame (SPDF). Ich möchte sie in einzelne SPDF zusammenführen (wie es in ArcGis mit Tool Merge sehr einfach ist ). Ich bin mir sicher, dass es eine einfache Möglichkeit geben sollte, dies in …
Diese Frage wurde von Stack Overflow migriert, da sie in Geographic Information Systems Stack Exchange beantwortet werden kann. Vor 4 Jahren migriert . Ich versuche, eine räumliche Verknüpfung zwischen Punktdaten und Polygondaten durchzuführen. Ich habe Daten, die die räumlichen Koordinaten eines Ereignisses in meiner CSV-Datei A angeben, und eine andere …
Ich versuche, eine Tabelle direkt aus einer ESRI-File-Geodatabase in R zu lesen. Eine Beispieldatendatei kann hier heruntergeladen werden . Die Datenbank enthält eine Point-Feature-Class (Zone9_2014_01_Broadcast) und zwei verknüpfte Tabellen (Zone9_2014_01_Vessel und Zone9_2014_01_Voyage). Sie können das Shapefile in R mit dem folgenden Befehl readOGRaus dem rgeosPaket lesen : library(rgeos) library(downloader) download("https://coast.noaa.gov/htdata/CMSP/AISDataHandler/2014/01/Zone9_2014_01.zip", …
Gibt es eine Funktion zum Zuschneiden von SF-Kartenobjekten, ähnlich wie maptools::pruneMap(lines, xlim= c(4, 10), ylim= c(10, 15))sie für SpatialPolygon oder SpatialLine verwendet werden? Ich überlege, st_intersection()aber es kann einen richtigen Weg geben.
Ist es bei einer großen Stichprobe (~ 1 Million) ungleichmäßig verteilter Punkte möglich, ein unregelmäßiges Raster (in der Größe, aber möglicherweise auch unregelmäßig in der Form?) Zu generieren, das eine festgelegte Mindestanzahl von n Punkten enthält? Es ist für mich weniger wichtig, wenn generierte "Zellen" eines solchen Gitters genau n …
Wie kann ich Lon-Lat-Punkte in einfache Features (sfg) konvertieren und sie dann in eine einfache Feature-Sammlung (sfc) einfügen? Hier ist ein MWE, der nicht funktioniert, dem ich aber am nächsten gekommen bin. library(data.table) library(sf) # The DT data.table is the data I have (but 10,000s of rows, each row is …
Ich zeichne einige Geolokationen (lon, lat)von Tweets, die mit der Stream-API gesammelt wurden, mit einer Länge von 11 km. Radius. Das folgende Bild zeigt zwei Karten mit denselben Daten, wobei jedoch für jede Beobachtung ein anderes Alpha / eine andere Transparenz verwendet wird. Auf der Karte auf der linken Seite …
coast<-readShapeSpatial("coastline.shp") landc<-readShapeSpatial("landcover.shp") ro<-readShapeSpatial("roads.shp") bc<-gBuffer(ro,width=100) landc$ratings=1 landc$ratings[landc$LANDUSE_ID==4]=0 Oben nehme ich eine Kategorie mit 4 und setze sie in der neuen Spalte als 0. An dieser Stelle möchte ich die Spalte ratingsauch für die benennen bc, wo es 0 nehmen wird, wenn es innerhalb des Puffers ist und 1, wenn es außerhalb …
Ich habe eine Reihe von Rasterbildern, die einen bestimmten Monat über die Jahre darstellen, und ich möchte mithilfe eines Shapefiles eine Zeitleiste der Durchschnittswerte eines Bereichs erstellen. Wie extrahiere ich die Werte aus den Rastern und importiere sie in R, sodass ich sie verwenden kann?
Morans I , ein Maß für die räumliche Autokorrelation, ist keine besonders robuste Statistik (es kann empfindlich auf verzerrte Verteilungen der räumlichen Datenattribute reagieren). Was sind robustere Techniken zur Messung der räumlichen Autokorrelation? Ich interessiere mich besonders für Lösungen, die in einer Skriptsprache wie R leicht verfügbar / implementierbar sind. …
In Zillow gibt es eine Reihe von Shapefiles für verschiedene Stadtteile großer US-amerikanischer Städte. Ich wollte mit R überprüfen, ob in bestimmten Stadtteilen bestimmte Gebäude vorhanden sind: library(rgeos) library(sp) library(rgdal) df <- data.frame(Latitude =c(47.591351, 47.62212,47.595152), Longitude = c(-122.332271,-122.353985,-122.331639), names = c("Safeco Field", "Key Arena", "Century Link")) coordinates(df) <- ~ Latitude …
Ich habe zwei Polygone. Eines enthält Felder (X, Y, Z) und das andere enthält Bodentypen (A, B, C, D). Ich möchte wissen, in welchem Gebiet jedes Feldes sich welcher Bodentyp befindet. Ich habe folgendes versucht: library(rgdal) library(rgeos) Field<-readOGR("./","Field") Soil<-readOGR("./","Soil") Results<-gIntersects(Soil,Field,byid=TRUE) rownames(Results)<-Field@data$FieldName colnames(Results)<-Soil@data$SoilType > Results A B C D Z TRUE …
Ich versuche, eine Vereinigung auf einem gemeinsamen Feld durchzuführen, nachdem zwei benachbarte Shapefiles zusammengeführt wurden. Zwischen den Shapefiles befindet sich mindestens ein dünner Streifen. Wenn ich versuche, eine Gewerkschaft zu gründen, wird der folgende Fehler angezeigt: Fehler in createPolygonsComment (p): rgeos_PolyCreateComment: Verwaiste Bohrung, enthält kein Polygon für Bohrung bei Index …
Mein Ziel ist es, ein vorhandenes Shapefile durch Zusammenführen bestimmter Polygone zu ändern. Nachdem ich das Shapefile importiert und den Befehl UnionSpatialPolygons verwendet habe, erhalte ich den gewünschten Polygonumriss. Dies ist jedoch jetzt ein SpatialPolygons-Objekt und kein SpatialPolygonsDataFrame, sodass ich es nicht mit writeOGR in ein Shapefile exportieren kann. Wie …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.