La traduzione da GeoTIFF a ENVI GDAL provoca la rotazione
Ho un GeoTIFF ( puoi scaricarlo qui ~ 16 MB) che voglio convertire in ENVI raw . Sto usando questo comando GDAL:
gdal_translate -ot UInt16 -scale -of ENVI -outsize 1025 1025 dgm1_32372_5680_2_nw.tif dgm1_32372_5680_2_nw.raw
Il comando funziona bene, ma dà un avviso:
Warning 1: Geotransform matrix has non rotational terms
Quando confronto il file *.tif e il file *.raw in QGIS (3.14), il file *.raw è stato in qualche modo ruotato.
Quindi ho cercato su Google e ho trovato questa domanda . Ok, diamo un'occhiata:
grep "map info" dgm1_32372_5680_2_nw.hdr
dà:
map info = {Transverse Mercator, 1, 1, 32371999.5, 5679999.5, 1.95121951219512, 1.95121951219512, rotation=90}
Ma secondo questo post sul blog il bug è stato corretto con GDAL 2.20. Sto usando GDAL 3.0.4. Quindi qualcuno sa quale potrebbe essere il problema?
Risposte
Devi gdalwarp per correggere la rotazione:
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.