Propongo dos soluciones: la primera con QGIS, la segunda con Python (GDAL).
Solución usando QGIS
En QGIS puede crear un mosaico VRT.
Siga este procedimiento (vea la imagen a continuación):
- Cargue el ráster en el Panel de capas ;
- Haga clic derecho sobre él y elija
Save As...;
- Marque la
Create VRTopción;
- Elija la carpeta donde se guardarán sus resultados;
- Establezca la extensión (si desea trabajar en todo el ráster, no modifique nada);
- Elija si usa la resolución actual (sugiero dejarla como predeterminada);
- Establezca el número máximo de columnas y filas (en su caso, debe ser 20000 columnas y 2000 filas);
- Presione el
OKbotón
Por ejemplo, el uso de los parámetros en el cuadro de diálogo anterior en este ráster de muestra (los parámetros que configuro se eligen al azar):
generará 100 mosaicos en la ruta especificada en el paso 4:
Al cargarlos en QGIS, se ven así:
Como @ bugmenot123 dijo correctamente en los comentarios, el resultado parece extraño solo porque el estilo de cada imagen se ajusta a la distribución de valores por imagen (pero los datos están perfectamente bien).
Solución usando Python (GDAL)
Otra forma de obtener el mismo resultado es el uso de GDAL ( gdal_translate ).
Con referencia al mismo ejemplo descrito anteriormente, puede 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, necesita adaptar los valores a su caso específico.