ฉันใช้ gdalwarp เพื่อจัดการกระเบื้อง SRTM ใกล้กับ dateline (เช่น 180 °หรือที่เรียกว่า antimeridian) ไทล์ SRTM มีการทับซ้อนกันเล็กน้อย (1/2 พิกเซล) เล็กน้อยกับเมริเดียน คุณสามารถดูได้โดยใช้ gdalinfo:
gdalinfo S16W180.hgt
Driver: SRTMHGT/SRTMHGT File Format
Files: S16W180.hgt
Size is 1201, 1201
[...]
Lower Left (-180.0004167, -16.0004167) (180d 0' 1.50"W, 16d 0' 1.50"S)
Upper Right (-178.9995833, -14.9995833) (178d59'58.50"W, 14d59'58.50"S)
[...]
ดังนั้นแหล่งขยายช่วงวันที่ด้วยจำนวนเล็กน้อย
สิ่งนี้ทำให้เกิดปัญหากับ gdalwarp ซึ่งจบลงด้วยการสร้างผลลัพธ์ที่ครอบคลุมทั่วโลก
gdalwarp -t_srs "epsg:900913" S16W180.hgt test.tif
gdalinfo test.tif
Driver: GTiff/GeoTIFF
Files: test.tif
Size is 1703, 5
[...]
Lower Left (-20037508.330,-1806798.473) (180d 0' 0.00"W, 16d 7'13.00"S)
Upper Right (20032839.451,-1689152.120) (179d57'29.01"E, 15d 5'45.84"S)
โปรดสังเกตว่าเส้นลองจิจูดยาว (เกือบ) ทั่วทั้งโลกและจำนวนบรรทัดมีขนาดเล็กโดยไม่คาดคิด (5)
นี่เป็นข้อผิดพลาดใน gdalwarp หรือไม่? ถ้าไม่ใช่ตัวเลือกที่ถูกต้องในการส่งผ่านไปยัง gdalwarp เพื่อรับเอาต์พุตที่เหมาะสมคืออะไร