Als «gdalwarp» getaggte Fragen

Ich habe ein Geotiff-Rasterbild mit einem Koordinatensystem mit Längen von 0 bis 360. Die horizontale Bildmitte beträgt 180 Längen. Siehe Bild unten: Ich möchte es in EPSG umwandeln: 4326 SRS mit -180 180 Längenbereich. Und ich möchte, dass die Bildmitte auf dem Greenwich-Meridian (0) liegt. Ich...

Ich versuche, meine Rasterdatei gemäß einer Shapefile-Geometrie zu beschneiden. Ich benutze den folgenden Code gdalwarp -cutline INPUT.shp INPUT.tif OUTPUT.tif Es entsteht jedoch eine schwarze Farbe außerhalb der Shapefile-Geometrie. Ich gebe hier das Beispiel. Im ersten Bild, das ich ausschneiden...

Ich versuche, die GDAL-Pythonbindungen neu zu projizieren / zu resampeln, erhalte jedoch geringfügig andere Ergebnisse als mit dem Befehlszeilendienstprogramm gdalwarp. Siehe Update unten für ein kürzeres Beispiel Dieses Skript veranschaulicht den Python-Ansatz: from osgeo import osr, gdal import...

Nachdem ich eine Reihe von Rastern mit gdalwarpprojiziert und am Raster (über -tap) ausgerichtet hatte, stellte ich fest, dass die Ausgabe-Raster erheblich größer waren als die ursprünglichen Raster. Eine ziemlich gründliche Websuche ergab dieses Trac-Problem : Frank Warmerdam erklärte den Grund:...

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 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...

In QGIS 2.XX hatten die Werkzeuge der Dropdown-Menüleiste unten eine bearbeitbare Aufrufbox, in der Befehle, die in der Benutzeroberfläche nicht verfügbar waren, manuell eingegeben werden konnten (durch Klicken auf eine Bleistift-Schaltfläche). In 3.0 sind die Werkzeuge in der Menüleiste mit den...

Ich habe einige Bilder mit GDAL mosaikiert und möchte das Endergebnis verbessern, indem ich eine verblassende / allmähliche Alpha-Ebene in Richtung der Bildkante verwende, um die scharfen Kanten in der Mitte des Mosaiks zu entfernen. Das Problem, das ich habe, ist, dass der Teil jedes einzelnen...

Ich habe versucht, den Code von QGis Save Raster als gerendertes Bild zu verwenden . iist eine Rasterebene: pipelayer = i pipeextent = pipelayer.extent() pipewidth, pipeheight = (pipelayer.width(), pipelayer.height()) piperenderer = pipelayer.renderer() pipeprovider = pipelayer.dataProvider() crs =...

Ich verwende gdalwarp, um SRTM-Kacheln in der Nähe der Datenlinie zu manipulieren (dh 180 °, auch bekannt als Antimeridian). SRTM-Kacheln haben eine sehr geringe Überlappung (1/2 Pixel) mit dem Meridian. Sie können dies mit gdalinfo sehen: gdalinfo S16W180.hgt Driver: SRTMHGT/SRTMHGT File Format...

Ich versuche, mit gdalwarp von Geostationary auf Lambert-konform umzuordnen. Meine Eingabedaten sind in netcdf und in geografischen Koordinaten (Grad) und ich möchte die neu zugeordneten Daten in netcdf ausgeben. Ich habe eine entsprechende vrt-Datei für die eingegebenen netcdf-Daten erstellt....

Ich möchte ein Raster mit pythonund georeferenzieren GDAL. Mein aktueller Ansatz ist es, eine hässliche Liste von Bodenkontrollpunkten aufzurufen gdal_translateund zu gdalwarpverwenden os.system. Ich würde wirklich gerne einen Weg finden, dies von Haus aus zu tun python. Dies ist der aktuelle...

Ich habe ein bisschen Mühe gehabt, Kacheln für ein hochauflösendes Bild zu erzeugen, das wir haben. Das aktuelle Bild ist ein sehr großes (+ 20 GB) Bild, das als GeoTiff-Datei gespeichert wird. Ich möchte die Kacheln mit dem Befehlszeilenprogramm gdal2tiles generieren und sie dann in Cäsium öffnen...

Ich benutze Python und QGIS 2.0. Ich versuche, Raster in einem Ordner nach einem Polygon-Feature zu schneiden. Es ist das erste Mal, dass ich "PyQGIS" benutze (sagen wir mal). Wie auch immer, ich bekomme mein einfaches Skript nicht zum Laufen, jeder Vorschlag wäre sehr dankbar! import qgis.core,...

Ich versuche ein Raster zu projizieren. In R gibt es dazu die projectRaster()Funktion (unter einem vollständig reproduzierbaren Beispiel): # example Raster require(raster) r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40) r <- setValues(r, 1:ncell(r)) projection(r) # project...

Ich versuche, einige Rasterdaten (alle .tif) mit gdalwarp in Python bis zu einer gröberen Auflösung (von 0,05 bis 0,25 Grad) zusammenzufassen, aber der Befehl funktioniert nicht. Anstatt eine Ausgabe mit einem großen Wertebereich zu erhalten, sind alle Werte der Ausgabe 0. Die Auflösung und die...