Tôi đề xuất hai giải pháp: giải pháp thứ nhất sử dụng QGIS, giải pháp thứ hai sử dụng Python (GDAL).
Giải pháp sử dụng QGIS
Trong QGIS bạn có thể tạo khảm VRT.
Vui lòng làm theo quy trình này (xem hình ảnh bên dưới):
- Tải raster trong Bảng điều khiển Lớp ;
- Nhấp chuột phải vào nó và chọn
Save As...;
- Kiểm tra
Create VRTtùy chọn;
- Chọn thư mục lưu kết quả đầu ra của bạn;
- Đặt phạm vi (nếu bạn muốn làm việc trên toàn bộ raster, đừng sửa đổi bất cứ điều gì);
- Chọn nếu sử dụng độ phân giải hiện tại (tôi đề nghị để mặc định như mặc định);
- Đặt số lượng cột và hàng tối đa (trong trường hợp của bạn, nó phải là 20000 cột và 2000 hàng);
- Nhấn
OKnút.
Ví dụ: việc sử dụng các tham số trong hộp thoại trên trên raster mẫu này (các tham số tôi đặt được chọn ngẫu nhiên):
sẽ tạo ra 100 ô trong đường dẫn được chỉ định ở bước 4:
Đang tải chúng trong QGIS, chúng trông như thế này:
Như @ bugmenot123 đã nói chính xác trong các bình luận, kết quả có vẻ kỳ lạ chỉ vì phong cách của mỗi hình ảnh phù hợp với phân phối giá trị trên mỗi hình ảnh (nhưng dữ liệu hoàn toàn tốt).
Giải pháp sử dụng Python (GDAL)
Một cách khác để có được kết quả tương tự là sử dụng GDAL ( gdal_translate ).
Với tham chiếu đến ví dụ tương tự được mô tả ở trên, bạn có thể sử dụng tập lệnh này:
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)
Bạn rõ ràng cần phải điều chỉnh các giá trị cho trường hợp cụ thể của bạn.