Я пропоную два рішення: перше, що використовує QGIS, друге - за допомогою Python (GDAL).
Рішення за допомогою QGIS
У QGIS ви можете створити мозаїку VRT.
Дотримуйтесь цієї процедури (див. Зображення нижче):
- Завантажте растр на панель шарів ;
- Клацніть на ньому правою кнопкою миші та виберіть
Save As...;
- Перевірте
Create VRTваріант;
- Виберіть папку, в якій будуть збережені ваші результати;
- Встановіть масштаб (якщо ви хочете працювати над усією растрою, нічого не змінюйте);
- Виберіть, чи використовуєте поточну роздільну здатність (пропоную залишити її за замовчуванням);
- Встановіть максимальну кількість стовпців та рядків (у вашому випадку це має бути 20000 стовпців та 2000 рядків);
- Натисніть
OKкнопку.
Наприклад, використання параметрів у наведеному вище діалоговому вікні цього растрового зразка (параметри, які я встановлюю, вибираються випадковим чином):
буде генерувати 100 плиток на шляху, визначеному на кроці 4:
Завантажуючи їх у QGIS, вони виглядають приблизно так:
Як правильно сказано в коментарях @ bugmenot123, результат виглядає дивним лише тому, що стиль кожного зображення підходить саме до розподілу значень на зображення (але дані ідеально чудові).
Рішення за допомогою Python (GDAL)
Ще один спосіб отримати той же результат - використання GDAL ( gdal_translate ).
Посилаючись на той самий приклад, описаний вище, ви можете використовувати цей сценарій:
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)
Вам, очевидно, потрібно адаптувати значення до вашого конкретного випадку.