Ich arbeite an einem DEM in einem Gebiet, in dem es wenig Erleichterung gibt. Ich möchte das Flussnetz in das DEM "einbrennen", damit ich die Ansammlung und Länge der Strömung genau berechnen kann. Ich verwende ArcGIS Desktop
Ist es möglich, ein Shapefile pro Feature in Python zu teilen? (Am besten wäre eine Lösung, bei der ich die resultierenden Vektorobjekte vorübergehend im Speicher anstatt auf der Festplatte speichern kann). Der Grund: Ich möchte die Funktion gdal rasterizeLayer mit mehreren unterschiedlichen...
Ich habe ein Linien-Shapefile, das ein Straßennetz darstellt. Ich möchte diese Daten rastern, wobei die resultierenden Werte im Raster die Gesamtlänge der Linien anzeigen, die in die Rasterzelle fallen. Die Daten befinden sich in der British National Grid-Projektion, sodass die Einheiten Meter...
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...
Ich möchte ein Raster mit einer Auflösung von 25 Metern × 25 Metern erstellen, wobei jede Zelle den Abstand zur nächsten Küste enthält, der aus der Mitte der Zelle berechnet wird. Dazu habe ich nur ein Shapefile der neuseeländischen Küsten . Ich habe versucht, Dominic Royes Tutorial zu folgen, um...
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:...
Ich habe eine Punktebene, die die Verteilung einer Art zeigt. Ist es in QGIS möglich, diese Punktebene in ein Rasterraster zu konvertieren, in dem der Wert jeder Zelle der Anzahl der Punkte in dieser Zelle entspricht? Bisher konnte ich nur Attributwerte an die Rasterzellen anhängen. Ich habe meiner...
Ich konvertiere Vektor in Raster in R. Der Prozess war jedoch zu lang. Gibt es die Möglichkeit, das Skript in Multithread- oder GPU-Verarbeitung zu versetzen, um es schneller zu machen? Mein Skript zum gerasterten Vektor. r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do...