Meine Frage bezieht sich auf das Extrahieren von Werten aus Rastern an der Position von Punkten. Mit der Funktion extrahieren ist dies sehr einfach, und die Funktion gibt mir einen Datenrahmen mit den Werten aller Variablen in den Punkten. Ich möchte in diesem Datenrahmen die Koordinaten jedes...
Siehe diesen Link für weitere Details. Das Problem: Ich möchte ein kontinuierliches Raster (eines ohne Attributtabelle) Zelle für Zelle durchlaufen und den Wert der Zelle abrufen. Ich möchte diese Werte übernehmen und Bedingungen für sie ausführen und die unten beschriebenen Schritte der...
Ich habe versucht, ein Rasterbild aus einer Datenbank mit unregelmäßigen Abständen zu erstellen. Die Daten sehen aus wie- > head(s100_ras) x y z 1 267573.9 2633781 213.29545 2 262224.4 2633781 69.78261 3 263742.7 2633781 51.21951 4 259328.4 2633781 301.98413 5 264109.8 2633781 141.72414 6...
Ich möchte zwei Polygone erstellen. Eine davon ist die rechteckige Ausdehnung meines Rasters. Ich weiß, dass ich das Ausmaß des Rasters ermitteln kann, r <- raster("band5.tif") e <- extent(r) plot(e)wenn ich Folgendes verwende : Wie kann ich EIN Polygon von der Grenze des Rasters wie unten...
Ich habe Python-Code, mit dem Raster-Attributtabellen für GeoTiffs mit GDAL erstellt werden. In ArcGIS kann ich diese Attributtabellen gut sehen und sie verhalten sich wie ich es mag (als Klassifizierungen), aber in QGIS scheint das gesamte Konzept einer Attributtabelle zu fehlen. Gibt es eine...
Ich hatte nur begrenzten Erfolg beim Erstellen von Polygonen mit glatten Konturen aus Rasterebenen. Ich habe Konturlinien und Grenzpolylinien erstellt und dann QGIS Union verwendet. Dann benutze ich QGIS Polygonize. Verwenden Sie dann die Zonenstatistik für das ursprüngliche Raster, um den Wert für...
Ich wollte Raster in QGIS in Point Vector konvertieren. Ist es möglich? QGIS hat die Option, es in Polygon zu konvertieren, aber ich habe kein Werkzeug gefunden, um es in Punkt zu konvertieren. Kann mir jemand
Ich habe zwei LiDAR-Dateien (.las), eine ist original, sagen wir mit X Punkten. Und die andere ist eine Kopie der ersten .las-Datei, jedoch mit Y-Punkten, wobei Y kleiner als X ist. Jetzt möchte ich die Digital Elevation Models (DEMs) dieser beiden .las-Dateien vergleichen und veranschaulichen, wie...
Ich habe ein Arc / Info-Binärgitter - speziell ein ArcGIS-Flow-Akkumulations-Raster - und möchte alle Zellen mit einem bestimmten Wert (oder einem Wertebereich) identifizieren. Letztendlich hätte ich gerne ein Shapefile mit Punkten, die diese Zellen darstellen. Ich kann QGIS verwenden, um die...
Ich arbeite an einem Projekt, das sowohl "Rasterdaten" (Luftbilder und DEMs) als auch "Gitterdaten" umfasst, von denen ich verstehe, dass es sich um andere Arten von numerischen Messungen handelt, die an Gitterpunkten durchgeführt werden und die möglicherweise nicht bildbezogen sind (zB...
Ich habe ein Raster-Dataset im ASCII-Format. Mit Python muss ich die Werte minund maxim Dataset ermitteln. Mir wurde gesagt, dass die Kopfzeileninformationen der Schlüssel sind, der Dinge wie Anzahl der Zeilen / Spalten, Zellengröße usw. enthält. Können Sie nicht einfach die Header-Informationen...
Ich habe ein Raster von Lebensraumtypen für ein bestimmtes Gebiet in Schottland. Ich muss zukünftige Lebensraumszenarien mit Änderungen im Lebensraum erstellen, um die Lebensfähigkeit der Population einer Vogelart zu beurteilen. Zum Beispiel könnte es in Zukunft 10% mehr Forstwirtschaft in der...
Ich versuche herauszufinden, warum der Slot @ data @ values beim Einlesen eines NDVI-Rasters erst dann die tatsächlichen Werte enthält, wenn ich sie manuell einstelle. Beispielsweise: NDVI <- raster("./filename.tif", crs="+proj=longlat +datum=WGS84") NDVI@data@values ## returns: logical(0)...
Ich frage mich nur, ist es immer besser, mehr als 3 GCPs zu haben? Sollten die GCPs gleichmäßig verteilt sein? Wo finde ich einige allgemeine Ressourcen zu dieser Theorie? Anders als beim Ausprobieren eines Fehlers, der mir einfach beibringt, was zu diesem Zeitpunkt funktioniert und nicht, was in...
Ich habe einen Polygon-Layer aus Rastern erstellt und nur eine einzige Klassifizierungs-ID (und berechnete Acres) für jedes Polygon. Ich möchte alle berührenden Merkmale in größere, einzelne Polygone "auflösen". Die 3 rot eingekreisten Bereiche (unten) enthalten grüne Polygone / Pixel, die nur die...
Gibt es eine Möglichkeit, dieselbe Aufgabe wie das Dienstprogramm gdalbuildvrt mithilfe von GDAL-Python-Bindungen auszuführen? Bisher habe ich keine andere Möglichkeit gefunden, als ein vrt eines einzelnen Datasets zu erstellen und die XML manuell zu bearbeiten. Ich möchte ein vrt aus mehreren...
Ich habe den Quellcode mehrerer QGis-1.7.4-Rasterfilter gelesen, die Steigung, Aspekt und Krümmung berechnen. Der Filter enthält eine Formel zur Berechnung der Gesamtkrümmung, die mich verwundert. Die Quelldatei befindet sich in der aktuellen Version von QGis mit dem folgenden Pfad: qgis-1.7.4 /...
Ich muss es in eine UTM umwandeln, um den Puffer funktionsfähig zu machen. wets<-readOGR(dsn=".",layer="shapefile") r.raster <- raster() extent(r.raster) <- extent(wets) res(r.raster) <- 100 wets.r <- rasterize(wet,r.raster) plot(wets.r) wetsbuf<-buffer(wets.r,width=500) Während...
Ich versuche, Rasterdaten in Polygondaten zu konvertieren. Bisher verwendete ich shapefile.py, um Punktdaten für einzelne Lat / Lon im Raster zu erstellen. Aber das wird riesige Daten erzeugen und das Abfragen wird ein großes Problem sein. Es wäre also besser, wenn ich die Rasterdaten in...
Ich versuche, ein Shapefile mit GDALs RasterizeLayer in ein Raster zu brennen. Ich erstelle ein Interessenbereich-Raster aus einem anderen Shapefile bei einer bestimmten Pixelgröße. Dieser AOI dient dann als Basis für alle folgenden Rasterisierungen (gleiche Anzahl von Spalten und Zeilen, gleiche...