ฉันเสนอวิธีแก้ปัญหาสองวิธี: อันแรกใช้ QGIS, อันที่สองใช้ Python (GDAL)
โซลูชันที่ใช้ QGIS
ใน QGIS คุณสามารถสร้างโมเสก VRT
โปรดปฏิบัติตามขั้นตอนนี้ (ดูภาพด้านล่าง):
- โหลดแรสเตอร์ในที่แผงเลเยอร์ ;
- คลิกขวาและเลือก
Save As...;
- ตรวจสอบ
Create VRTตัวเลือก;
- เลือกโฟลเดอร์ที่จะบันทึกผลลัพธ์ของคุณ
- กำหนดขอบเขต (ถ้าคุณต้องการทำงานกับแรสเตอร์ทั้งหมดอย่าแก้ไขอะไรเลย);
- เลือกว่าใช้ความละเอียดปัจจุบันหรือไม่ (ฉันแนะนำให้ปล่อยไว้เป็นค่าเริ่มต้น);
- กำหนดจำนวนคอลัมน์และแถวสูงสุด (ในกรณีของคุณควรเป็นคอลัมน์ 20,000 และ 2,000 แถว)
- กด
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)
เห็นได้ชัดว่าคุณต้องปรับค่าให้เข้ากับกรณีเฉพาะของคุณ