Verwendung von gdal2tiles für ein benutzerdefiniertes TIFF-Image, das von Anbietern zum Generieren von Kacheln erhalten wurde

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 der TMS-Bildanbieter die Kacheln bereitstellt. Bei Verwendung von gdalinfo sind hier einige Details des Bildes:

Driver: GTiff/GeoTIFF
Files: image.tif
Size is 52250, 56119
Coordinate System is:
PROJCS["WGS 84 / UTM zone 35S",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",27],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",10000000],
    UNIT["meters",1],
    AUTHORITY["EPSG","32735"]]
Origin = (606276.000000000000000,7197873.000000000000000)
Pixel Size = (0.500000000000000,-0.500000000000000)
Metadata:
  AREA_OR_POINT=Area
  TIFFTAG_MAXSAMPLEVALUE=13165
  TIFFTAG_MINSAMPLEVALUE=1
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
  TIFFTAG_SOFTWARE=ERDAS IMAGINE
  TIFFTAG_XRESOLUTION=1
  TIFFTAG_YRESOLUTION=1
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  606276.000, 7197873.000) ( 28d 3'21.59"E, 25d19'55.12"S)
Lower Left  (  606276.000, 7169813.500) ( 28d 3'29.55"E, 25d35' 7.17"S)
Upper Right (  632401.000, 7197873.000) ( 28d18'55.92"E, 25d19'47.60"S)
Lower Right (  632401.000, 7169813.500) ( 28d19' 5.85"E, 25d34'59.57"S)
Center      (  619338.500, 7183843.250) ( 28d11'13.23"E, 25d27'27.58"S)
Band 1 Block=512x512 Type=UInt16, ColorInterp=Gray
Band 2 Block=512x512 Type=UInt16, ColorInterp=Undefined
Band 3 Block=512x512 Type=UInt16, ColorInterp=Undefined
Band 4 Block=512x512 Type=UInt16, ColorInterp=Undefined

Mein erster Versuch war, das Bild mit gdal_translate zu georeferenzieren und dann mit gdalwarp die Projektion in EPSG: 3857 zu ändern, wie von Cäsium gefordert (siehe API-Referenz).

gdal_translate -of VRT -a_srs EPSG:4326 -gcp 606275 7197875 28.055987 -25.331974 -gcp 606275 7169814 28.058200 -25.585326 -gcp 632400.5 7197875 28.31553 -25.329876 -gcp 632400.5 7169814 28.318286 -25.583209 image.tif newImage1.vrt
gdalwarp -of VRT -t_srs EPSG:3857 newImage1.vrt newImage2.vrt

Ich erhalte jedoch viele der folgenden Fehler:

FEHLER 1: Breiten- oder Längengrad haben die Grenzwerte überschritten

Eine andere Methode, die ich versuchte, war, gdal2tiles direkt zu verwenden und die Kacheln zu generieren:

gdal2tiles.py image.tif

Dadurch wurde ein Ordner erstellt, wobei ein Unterordner (mit 18 gekennzeichnet) die einzige Zoomstufe ist, bei der Kacheln erstellt wurden. Die Bilder, die ich hier bekomme, sind jedoch völlig "falsch" und "verschwommen".

Ein Beispiel für eine der Kacheln:

Irgendwelche Vorschläge zum Generieren von Kacheln für dieses Bild Großes Bild eines bestimmten Bereichs mit gdal2tiles, damit ich es in Cäsium laden und anzeigen kann?

Aktualisieren

Nachdem ich den Vorschlag von @ iant ausprobiert hatte, verwendete ich die folgenden Befehle:

gdalwarp -co TILED=YES -co COMPRESS=DEFLATE -co BIGTIFF=YES -t_srs EPSG:3857 image.tif newImage.tif

Dieser Befehl funktionierte einwandfrei bis zum Ende, wo ich den folgenden Fehler bekam:

FEHLER 1: TIFFFillTile: Lesefehler in Zeile 43520, Spalte 47104; bekam 35788250 Bytes, erwartet 37421449

Ich bin mir nicht sicher, was dieser Fehler bedeutete. Ich habe ihn für den Moment verlassen und trotzdem ein endgültiges Bild "newImage.tif" erhalten, das vom gdalwarp-Schritt erstellt wurde. Damit habe ich gdal2tiles.py aufgerufen

gdal2tiles.py newImage.tif

Dies erzeugte einen Ordner mit den Unterordnern 10-18 (und nicht nur einer Zoomstufe 18, wie ich sie zuvor erhalten habe). Es liest sich auch perfekt in Cäsium, ohne Konsolenfehler, aber das Bild sieht immer noch "falsch" aus:

Ich denke, mein Problem könnte sein, wie @ user30184 vorgeschlagen hat "... Quelldaten passen nicht gut für gdal2tiles." Bis unser Anbieter uns jedoch etwas zur Verwendung mit gdal zur Verfügung stellen kann, ist dies alles, was ich habe.

Ich dachte darüber nach, vielleicht eine der Bands zu entfernen, um zu vermeiden, dass gdal die letzte Band als Alpha-Kanal unterbricht. Irgendwelche Vorschläge?

Ich denke, alles was Sie tun müssen, ist es neu zu projizieren mit:

gdalwarp -co TILED=YES -co COMPRESS=DEFLATE -t_srs EPSG:3857 newImage.tif image.tif

und dann kacheln:

gdal2tiles.py newImage.tif

Wenn Ihre Datei sehr groß ist, dauert es eine Weile.

Ich nehme an, Ihr Image ist eines der 4-Band-Produkte von Airbus DS:

http://www.intelligence-airbusds.com/de/4951-which-spectral-mode-do-i-choose

Gdal2tiles dient zum Aufteilen gängiger visueller Bilder in PNG-Kacheln. Solche Bilder verwenden 8 Bits pro Band und haben ein Band (Graustufen), 3 Bänder (Rot-Grün-Blau) von 4 Bändern (Reg-Grün-Blau + Alpha).

Ich würde sagen, dass Ihre Frage weitgehend irrelevant ist, da Ihre Quelldaten nicht gut für gdal2tiles geeignet sind. Sie können die unmittelbaren Probleme umgehen, die Sie gerade haben, aber das Endergebnis wird immer noch nicht gut sein, wenn Sie Ihre Daten nicht erneut verarbeiten.

Der Grund für eine nicht gut aussehende Kachel, die Sie Ihrer Frage beigefügt haben, kann sein, dass das vierte Datenband als Alphakanal interpretiert wird.