Ich schreibe ein einfaches Dienstprogramm, um Stapel von Multiband-Geotiff-Rasterdateien auf denselben (kleineren) Bereich zuzuschneiden. Mit gdalwarp kann ich eine Datei mithilfe eines Shapefiles mit einem einzelnen Polygonausschnitt ganz einfach zuschneiden: gdalwarp -cutline clipper.shp...
Ich bin sehr neu in GIS. Ich benutze gdal, um eine Landnutzungs- / Landbedeckungskarte einzulesen, und ich muss die Länge bestimmter Landbedeckungstypen auswählen, um sie in einem anderen Datensatz zu indizieren, der nur in Länge ausgedrückt wird. Unfortuantely, ich verstehe nicht , die Form der x...
Ich wollte ungefähr 20 Dem-Tiffs mit jeweils ungefähr 100 MB zusammenführen und habe gdalwarp verwendet (ich kenne die Option für virtuelle Dateien): gdalwarp $(list_of_tiffs) merged.tiff Dieser Befehl dauerte jedoch sehr lange und verarbeitete nach 15 Minuten das fünfte Bild aus dem Satz. Ich...
Ich möchte ein Webportal entwickeln, das die von QGIS dynamisch erstellten Bilder anzeigt, aber von QGIS bereitgestellte Tiff-Bilder sind zu groß und können mit Bildbetrachtern nicht auf einer Website angezeigt werden. Aus diesem Grund habe ich mich entschlossen, diese tiff-Bilder in jpg / jpeg...
Ich versuche, meinen Python-Code mit pypy zu optimieren . Während ich die meisten Bibliotheken importieren und erfolgreich ausführen kann, kann ich keine GDAL-Bibliotheken importieren. Dieselben Bibliotheken werden beim Import in Python 2.7 problemlos ausgeführt. Ich habe versucht, GDAL-Speicherort...
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 füge einige Shapefiles zusammen und ich hatte einige Probleme damit in QGIS, also verwende ich ogr2ogr direkt. Ich mache das (in einem Batch): ogr2ogr -overwrite %destination% %n1% ogr2ogr -update -append %destination% %n2% -nln all_new ogr2ogr -update -append %destination% %n3% -nln all_new...
Meistens bin ich ein Vektorkerl, aber für ein aktuelles Projekt muss ich mit ein paar alten gescannten Papierkarten arbeiten (ww2-Bombenschadenskarten für London, wenn Sie interessiert sind!) Wir haben die Karten gescannt und georeferenziert und möchten nun eine gekachelte Verbundebene für die...
Gibt es eine Möglichkeit, alle 0-Wert-Pixel in einem Raster mithilfe von GDAL auf Nodata zu setzen? Ich möchte Fillnodata verwenden, um einige Löcher zu füllen, aber die Löcher sind auf 0 gesetzt, daher muss ich sie zuerst
Ich habe ein ganzzahliges Raster, für das ich eine Raster-Attributtabelle mit Python und GDAL erstellen möchte. Ich kann eine GDAL-Rasterattributtabelle in Python wie folgt erstellen: >>> rat = gdal.RasterAttributeTable() Das funktioniert gut, wie wir sehen können: >>> rat...
Ich exportiere gerade eine Postgis-Tabelle mit diesen Tipps nach shp, kann aber mit derselben Bibliothek (OGR) keine shp nach Postgis importieren. Irgendeine Idee? Vielen Dank
In der Dokumentation für gdal_calc wird der Befehlszeilen-Rasterrechner mit Numpy-Syntax angegeben . Später gibt es einige Beispiele, wo in einem von ihnen: gdal_calc.py -A input.tif --outfile = result.tif --calc = "A * (A> 0)" --NoDataValue = 0 - bedeutet , dass Werte von Null und darunter auf...
Ich erstelle regelmäßig meine eigenen GeoTIFF-Raster mit GDAL in Python, zB: from osgeo import gdal from numpy import random data = random.uniform(0, 10, (300, 200)) driver = gdal.GetDriverByName('GTiff') ds = driver.Create('MyRaster.tif', 200, 300) band = ds.GetRasterBand(1) band.WriteArray(data)...
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 versuche, einige Rasterdaten mit ogr / gdal zu verarbeiten, und ich kann anscheinend nicht alle Kerne auf meinem Computer voll ausnutzen. Wenn ich den Prozess nur auf einem einzelnen Kern ausführe, wird dieser Kern zu 100% genutzt. Wenn ich versuche, mich in mehrere Kerne aufzuteilen (im...
Ich habe GDAL installiert und in der Befehlszeile überprüft, es funktioniert und wenn ich die Datei gdal2tiles.py in der Befehlszeile ausführe, funktioniert es, aber wenn ich dieselbe Datei gdal2tiles.py mit Python IDLE ausführe, erhalte ich "Kein Modul mit dem Namen" osgeo " '" Error. Kannst du...
Ich habe eine TIF-Bilddatei und eine .tfw-Weltdatei (erstellt aus OSM-Daten mit Mapertive). Wie kann ich das (unter Verwendung der Linux-Befehlszeile) in eine GeoTIFF-Datei konvertieren? (dh ich möchte eine Datei haben, nicht 2). Ich bin sicher, dass es mit GDAL oder so möglich
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 Rasterdaten in einer Postgres-Tabelle, die ich als Numpy-Array in Python aufnehmen möchte. Ich benutze psycopg2, um eine Verbindung zur Datenbank herzustellen. Ich kann die Daten herunterladen, sie werden jedoch als Zeichenfolge zurückgegeben (wahrscheinlich als serialisierte Binärdatei)....
Ich habe Probleme bei der Verwendung der Polygonize-Funktion in Python. Das Kochbuchbeispiel hierfür finden Sie hier . Der relevante Teil meines Codes ist: sourceRaster = gdal.Open('myraster.tif') band = sourceRaster.GetRasterBand(1) bandArray = band.ReadAsArray() outShapefile = "polygonized"...