Ich schlage zwei Lösungen vor: die erste mit QGIS, die zweite mit Python (GDAL).
Lösung mit QGIS
In QGIS können Sie ein VRT-Mosaik erstellen.
Bitte befolgen Sie diese Prozedur (siehe Bild unten):
- Laden Sie das Raster im Ebenenbedienfeld .
- Klicken Sie mit der rechten Maustaste darauf und wählen Sie
Save As...;
- Aktivieren Sie die
Create VRTOption.
- Wählen Sie den Ordner, in dem Ihre Ausgaben gespeichert werden sollen.
- Legen Sie die Ausdehnung fest (wenn Sie mit dem gesamten Raster arbeiten möchten, ändern Sie nichts).
- Wählen Sie, ob Sie die aktuelle Auflösung verwenden möchten (ich empfehle, die Standardauflösung beizubehalten).
- Legen Sie die maximale Anzahl von Spalten und Zeilen fest (in Ihrem Fall sollten es 20000 Spalten und 2000 Zeilen sein).
- Drücken Sie die
OKTaste.
Zum Beispiel die Verwendung der Parameter im obigen Dialogfeld für dieses Beispiel-Raster (die von mir festgelegten Parameter werden nach dem Zufallsprinzip ausgewählt):
generiert 100 Kacheln in dem in Schritt 4 angegebenen Pfad:
Sie werden in QGIS geladen und sehen folgendermaßen aus:
Wie @ bugmenot123 in den Kommentaren richtig sagte, sieht das Ergebnis seltsam aus, nur weil sich der Stil jedes Bildes an die Verteilung der Werte pro Bild anpasst (aber die Daten sind vollkommen in Ordnung).
Lösung mit Python (GDAL)
Eine andere Möglichkeit, das gleiche Ergebnis zu erzielen, ist die Verwendung von GDAL ( gdal_translate ).
Mit Bezug auf dasselbe oben beschriebene Beispiel können Sie dieses Skript verwenden:
import os, gdal
in_path = 'C:/Users/Marco/Desktop/'
input_filename = 'dtm_5.tif'
out_path = 'C:/Users/Marco/Desktop/output_folder/'
output_filename = 'tile_'
tile_size_x = 50
tile_size_y = 70
ds = gdal.Open(in_path + input_filename)
band = ds.GetRasterBand(1)
xsize = band.XSize
ysize = band.YSize
for i in range(0, xsize, tile_size_x):
for j in range(0, ysize, tile_size_y):
com_string = "gdal_translate -of GTIFF -srcwin " + str(i)+ ", " + str(j) + ", " + str(tile_size_x) + ", " + str(tile_size_y) + " " + str(in_path) + str(input_filename) + " " + str(out_path) + str(output_filename) + str(i) + "_" + str(j) + ".tif"
os.system(com_string)
Sie müssen die Werte natürlich an Ihren speziellen Fall anpassen.