İki çözüm öneriyorum: birincisi QGIS, ikincisi Python (GDAL) kullanıyor.
QGIS kullanarak çözüm
QGIS'te bir VRT mozaiği oluşturabilirsiniz.
Lütfen bu prosedürü izleyin (aşağıdaki resme bakın):
- Raster'i Katmanlar Paneline yükleyin ;
- Üzerine sağ tıklayın ve seçin
Save As...;
- Kontrol
Create VRTseçeneği;
- Çıktılarınızın kaydedileceği klasörü seçin;
- Kapsamı ayarlayın (tüm raster üzerinde çalışmak istiyorsanız, hiçbir şeyi değiştirmeyin);
- Mevcut çözünürlüğü kullanıp kullanmadığınızı seçin (varsayılan olarak bırakmanızı öneririm);
- Maksimum sütun ve satır sayısını ayarlayın (sizin durumunuzda, 20000 sütun ve 2000 satır olmalıdır);
- Düğmesine basın
OK.
Örneğin, bu örnek taramada yukarıdaki iletişim kutusundaki parametrelerin kullanılması (ayarladığım parametreler rastgele seçilir):
4. adımda belirtilen yolda 100 kutucuk oluşturur:
Onları QGIS'e yüklüyorlar, şöyle görünürler:
@ Bugmenot123 yorumlarda doğru bir şekilde söylediği gibi, sonuç her bir resmin stilinin görüntü başına değer dağılımına uyması nedeniyle garip görünüyor (ancak veriler mükemmel derecede iyi).
Python (GDAL) kullanarak çözüm
Aynı sonucu elde etmenin başka bir yolu da GDAL ( gdal_translate ) kullanılmasıdır.
Yukarıda açıklanan aynı örneğe istinaden, bu komut dosyasını kullanabilirsiniz:
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)
Değerleri özel durumunuza göre uyarlamanız gerekir.