Proponho duas soluções: a primeira usando QGIS, a segunda usando Python (GDAL).
Solução usando QGIS
No QGIS, você pode criar um mosaico VRT.
Por favor, siga este procedimento (veja a imagem abaixo):
- Carregue a varredura no painel Camadas ;
- Clique com o botão direito do mouse e escolha
Save As...;
- Marque a
Create VRTopção;
- Escolha a pasta onde suas saídas serão salvas;
- Defina a extensão (se você quiser trabalhar em toda a varredura, não modifique nada);
- Escolha se estiver usando a resolução atual (sugiro deixá-la como padrão);
- Defina o número máximo de colunas e linhas (no seu caso, deve ser 20000 colunas e 2000 linhas);
- Pressione o
OKbotão
Por exemplo, o uso dos parâmetros na caixa de diálogo acima nesta amostra de varredura (os parâmetros que eu defini são escolhidos aleatoriamente):
irá gerar 100 blocos no caminho especificado na etapa 4:
Carregando-os no QGIS, eles ficam assim:
Como @ bugmenot123 disse corretamente nos comentários, o resultado parece estranho apenas porque o estilo de cada imagem se encaixa na distribuição de valores por imagem (mas os dados estão perfeitamente corretos).
Solução usando Python (GDAL)
Outra maneira de obter o mesmo resultado é o uso de GDAL ( gdal_translate ).
Com referência ao mesmo exemplo descrito acima, você pode usar este script:
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)
Obviamente, você precisa adaptar os valores ao seu caso específico.