Terjemahan GeoTIFF ke ENVI GDAL menyebabkan rotasi
Saya memiliki GeoTIFF ( Anda dapat mengunduhnya di sini ~ 16MB) yang ingin saya ubah menjadi ENVI mentah . Saya menggunakan perintah GDAL ini:
gdal_translate -ot UInt16 -scale -of ENVI -outsize 1025 1025 dgm1_32372_5680_2_nw.tif dgm1_32372_5680_2_nw.raw
Perintah berjalan dengan baik, tetapi memberikan peringatan:
Warning 1: Geotransform matrix has non rotational terms
Ketika saya membandingkan file * .tif dan * .raw di QGIS (3.14), file * .raw entah bagaimana bisa diputar.
Jadi saya melakukan googling dan menemukan pertanyaan ini . Oke, mari kita lihat:
grep "map info" dgm1_32372_5680_2_nw.hdr
memberikan:
map info = {Transverse Mercator, 1, 1, 32371999.5, 5679999.5, 1.95121951219512, 1.95121951219512, rotation=90}
Tetapi menurut posting blog ini , bug telah diperbaiki dengan GDAL 2.20. Saya menggunakan GDAL 3.0.4. Jadi, apakah ada yang tahu apa masalahnya?
Jawaban
Anda perlu melakukan gdalwarp untuk memperbaiki rotasi:
gdalwarp -of VRT dgm1_32372_5680_2_nw.tif dgm1_32372_5680_2_nw.vrt
Creating output file that is 2000P x 2000L.
Processing input file dgm1_32372_5680_2_nw.tif.
gdal_translate -ot UInt16 -scale -of ENVI -outsize 1025 1025 dgm1_32372_5680_2_nw.vrt dgm1_32372_5680_2_nw.raw
Input file size is 2000, 2000
0...10...20...30...40...50...60...70...80...90...100 - done.