Saya mengusulkan dua solusi: yang pertama menggunakan QGIS, yang kedua menggunakan Python (GDAL).
Solusi menggunakan QGIS
Di QGIS Anda dapat membuat mosaik VRT.
Silakan ikuti prosedur ini (lihat gambar di bawah):
- Muat raster di Panel Layers ;
- Klik kanan padanya dan pilih
Save As...;
- Periksa
Create VRTopsi;
- Pilih folder tempat output Anda akan disimpan;
- Tetapkan batasnya (jika Anda ingin bekerja di seluruh raster, jangan modifikasi apa pun);
- Pilih jika menggunakan resolusi saat ini (saya sarankan untuk membiarkannya sebagai default);
- Tetapkan jumlah maksimum kolom dan baris (dalam kasus Anda, itu harus 20.000 kolom dan 2000 baris);
- Tekan
OKtombolnya.
Misalnya, penggunaan parameter dalam dialog di atas pada raster sampel ini (parameter yang saya set dipilih secara acak):
akan menghasilkan 100 ubin di jalur yang ditentukan pada langkah 4:
Memuatnya di QGIS, mereka terlihat seperti ini:
Seperti @ bugmenot123 dengan benar katakan di komentar, hasilnya terlihat aneh hanya karena gaya setiap gambar cocok dengan distribusi nilai per gambar (tetapi data baik-baik saja).
Solusi menggunakan Python (GDAL)
Cara lain untuk mendapatkan hasil yang sama adalah penggunaan GDAL ( gdal_translate ).
Dengan merujuk pada contoh yang sama yang dijelaskan di atas, Anda dapat menggunakan skrip ini:
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)
Anda jelas perlu menyesuaikan nilai-nilai dengan kasus spesifik Anda.