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 denke, diese srs ist …
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 möchte. Das zweite Bild ist das resultierende Raster, aber ich möchte …
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 numpy def reproject_point(point, srs, target_srs): ''' Reproject a pair of coordinates from …
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: "Bei sorgfältiger Prüfung besteht der Unterschied in der fraglichen …
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 beende es. Ich fand dann heraus, dass …
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 -crop_to_cutline input.tif output.tif Der tatsächliche Bereich, in den ich einen Clip einfügen möchte, wird anfangs jedoch immer …
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 GDAL-Werkzeugen in der Verarbeitungs-Toolbox identisch, was sinnvoller ist, aber …
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 Bildes mit der graduellen …
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 Files: S16W180.hgt Size is 1201, 1201 [...] Lower …
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. Gdalwarp gibt die netcdf-Datei aus, aber die Ausgabedaten sind alle Nullen …
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 Prozess, den ich verwende: import os os.system('gdal_translate -of GTiff …
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 und anzeigen, wobei …
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 to newproj <- "+init=epsg:4714" # using raster package to reproject pr1 <- projectRaster(r, …
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, qgis,utils QgsApplication.setPrefixPath("C:/OSGeo4W64/apps/qgis", True) …
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.