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 …
Ich habe zwei SpatialLinesObjekte in R erstellt : . Diese Objekte wurden folgendermaßen erstellt: library(sp) xy <- cbind(x,y) xy.sp = sp::SpatialPoints(xy) spl1 <- sp::SpatialLines(list(Lines(Line(xy.sp), ID="a"))) Jetzt möchte ich irgendwie zu dem Schluss kommen, dass dies dieselbe Linie ist, die gedreht und gespiegelt wird, und dass der Unterschied zwischen ihnen gleich …
Ich versuche, einige Schwerpunkte eines Shapefiles für dänische Gemeinden zu finden und anschließend die Fahrzeit zwischen ihnen zu ermitteln. Ich benutze Rdie readShapePolyFunktion von maptoolskombiniert mit der gCentroidFunktion von rgeosund alles funktioniert. Ich bekomme jedoch räumliche Punkte wie SpatialPoints: x y 1 571860.7 6225016 Coordinate Reference System (CRS) arguments: NA …
Wie finde ich einen maximalen Zellenwert aus einem gestapelten Raster? Rmax <- maxValue(RAD1998.all[[1]]) funktioniert aber gut Rmax <- maxValue(RAD1998.all[[2]]) gibt NA. Natürlich nicht in einem gestapelten Raster. Hier ist mein Code: RAD1998 <- raster(paste(getwd(), "/1998bil/1998ASC5min_ppt_spas1214_0001_19980202_0810_UTC.asc.bil", sep = "")) list.ras <- mixedsort(list.files(paste(getwd(), "/1998bil/", sep = ""), full.names = T, pattern = …
Ich habe die Oberfläche von Artenverteilungen berechnet (Zusammenführen von Polygonen aus Shapefiles), aber da diese Fläche aus ziemlich weit entfernten Polygonen bestehen kann, möchte ich ein Maß für die Dispersion berechnen. Was ich bisher getan habe, ist, die Schwerpunkte jedes Polygons abzurufen, den Abstand zwischen diesen zu berechnen und diese …
Alle, entschuldigen Sie die Störung, aber ich bin ziemlich neu mit einer entscheidenden Schwierigkeit: Ich möchte eine animierte Karte von Russin mit Veränderungen der Arbeitslosigkeit mit unterschiedlichen Jahren erstellen, wie z. Auf dem Bild sehen Sie Daten für ein Jahr require(sp) require(maptools) require(RColorBrewer) require(rgdal) rus<-url("http://www.filefactory.com/file/4h1hb5c1cw7r/n/RUS_adm1_RData") print(load(rus)) unempl <- read.delim2(file="C:\\unempl1.txt", header …
Ich bin daran interessiert, eine Rasterfläche zu erstellen, die die durchschnittliche Sonnenintensität pro Flächeneinheit für ein kleines Gartengrundstück darstellt. Ich stelle mir vor, die Eingaben in das Modell wären folgende: Grundstücksfläche - Nehmen wir der Einfachheit halber an, dass das Gartengrundstück kreisförmig ist. Durchschnittliche Baumhöhe des umgebenden Waldrandes Sonnenintensität bei …
Das beigefügte Bild zeigt eine Waldlücke mit roter Kiefer als Kreise und weißer Kiefer als Kreuze. Ich bin daran interessiert festzustellen, ob zwischen den beiden Kiefernarten eine positive oder negative Assoziation besteht (dh ob sie in denselben Gebieten wachsen oder nicht). Ich kenne Kcross und Kmulti im R-Spatstat-Paket. Da ich …
In der Dokumentation von habe spich gesehen, dass das coordsArgument der SpatialPolygonsFunktion nur 2 Spalten akzeptiert. Ist es nicht möglich, dreidimensionale Polygone zu erstellen? Ich mache räumlich-zeitliches Kriging mit Paket, gstatindem ich Zeit als dritte räumliche Dimension von a hinzufüge SpatialPointsDataFrame. Ich möchte das erweitern, um Kriging mit Polygonen als …
Ich benötige ein räumliches Raster als Master-Raster für verschiedene thematische Karten. Wie erstelle ich ein räumliches Gitter aus einem Raster, bei dem alle NA-Pixel verworfen werden?
In meinem Arbeitsbereich befinden sich 100 SpatialPointsDataFrames. Ich habe sie auf diese Weise geladen: filenames <- list.files(path="", pattern="XYhectareTravelTimes_ez+.*shp") for(i in filenames){ filepath <- file.path("/",i) assign(i, readShapePoints(filepath)) Sie heißen so: XYhectareTravelTimes_ez10.*shp XYhectareTravelTimes_ez11.*shp XYhectareTravelTimes_ez12.*shp etc. Wie kann ich sie in Raster konvertieren, die den Arbeitsbereich durchlaufen? Ich bin ein sehr neuer R-Benutzer …
Wie kann ich unter Bezugnahme auf meine vorherige Frage ein generiertes Rasterbild mit einem generierten Rasterbild in ein Format "geocodieren" / konvertieren, das als Layer auf einer Karte angezeigt werden kann? Kurz gesagt, ich möchte die Ergebnisse der kommerziellen SpatialKey- Software ( Screenshot ) neu erstellen . Ich verwende eine …
Wahrscheinlich eine sehr einfache Frage. Ich habe eine Liste von ungefähr tausend potenziellen geografischen Standorten (lat lon), und von denen muss ich die 200 Standorte auswählen, die "am weitesten verbreitet" sind. Ich denke, das wären die 200 Punkte mit der höchsten durchschnittlichen Gesamtentfernung. Denken Sie an Geschäfte in einer Stadt. …
Dies ist ein Crosspost von Stackoverflow . Ich weiß sehr wenig über GIS-Software und mache alle meine Mappings in R. Entschuldigung im Voraus, wenn dies eine zu grundlegende Frage ist. Angenommen, ich habe zwei Shapefiles aus verschiedenen Quellen, aber mit unterschiedlichen Attributen. Angenommen, einer ist für Verwaltungsgrenzen von Texas ( …
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.