Proponuję dwa rozwiązania: pierwsze za pomocą QGIS, drugie za pomocą Pythona (GDAL).
Rozwiązanie wykorzystujące QGIS
W QGIS możesz utworzyć mozaikę VRT.
Postępuj zgodnie z tą procedurą (patrz obrazek poniżej):
- Załaduj raster do panelu warstw ;
- Kliknij go prawym przyciskiem myszy i wybierz
Save As...;
- Sprawdź
Create VRTopcję;
- Wybierz folder, w którym zostaną zapisane dane wyjściowe;
- Ustaw zakres (jeśli chcesz pracować na całym rastrze, niczego nie modyfikuj);
- Wybierz, jeśli używasz bieżącej rozdzielczości (sugeruję pozostawić ją jako domyślną);
- Ustaw maksymalną liczbę kolumn i wierszy (w twoim przypadku powinno to być 20000 kolumn i 2000 wierszy);
- Naciśnij
OKprzycisk.
Na przykład użycie parametrów w powyższym oknie dialogowym na tym przykładowym rastrze (parametry, które ustawiam są wybierane losowo):
wygeneruje 100 płytek na ścieżce określonej w kroku 4:
Ładowanie ich do QGIS wygląda następująco:
Jak poprawnie napisano w komentarzach @ bugmenot123, wynik wygląda dziwnie tylko dlatego, że styl każdego obrazu pasuje do rozkładu wartości na obraz (ale dane są całkowicie w porządku).
Rozwiązanie wykorzystujące Python (GDAL)
Innym sposobem na uzyskanie tego samego wyniku jest użycie GDAL ( gdal_translate ).
W odniesieniu do tego samego przykładu opisanego powyżej możesz użyć tego skryptu:
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)
Oczywiście musisz dostosować wartości do konkretnego przypadku.