<div>We use geotif with units=seconds quite a lot.</div><div><br></div>I found that running :<div>gdalwarp -t_srs <srs-def> <input path> <output path></div><div><br></div><div>1. If target srs_def is a prj file with EPSG code for GCS WGS84 and units= degree and my input file is a geotiff with <font><span style="background-color:rgba(255,255,255,0)">GCS WGS84 and units= Seconds - the output file is a "normal" geotiff with GCS WGS84 and units= degree.</span></font></div>
<div><font><span style="background-color:rgba(255,255,255,0)">2. If target srs_def in a prj file with EPSG code for GCS WGS84 and units= SECONDS and my input file is a geotiff with GCS WGS84 and units= degree - the output is a geotiff with the units undetected correctly.</span></font></div>
<div><font><span style="background-color:rgba(255,255,255,0)">In order to correct the output file I used ArcCatalog and set the units to seconds.</span></font></div><div><font><span style="background-color:rgba(255,255,255,0)">After the correction the file (and the header was ok).</span></font></div>
<div><font><span style="background-color:rgba(255,255,255,0)">Note: the prj must include authority = EPSG and the code (in this case 4326) otherwise the geotiff is produced as Custom CS and it m</span></font><span style="background-color:rgba(255,255,255,0);font-size:small">ay or may not be identified correctly and fully by apps.</span></div>
<div><font><span style="background-color:rgba(255,255,255,0)">Also note that with default units you can run gdalwarp with EPSG:4326 instead of the prj file.</span></font></div><div><font><span style="background-color:rgba(255,255,255,0)"><br>
</span></font></div><div><font><span style="background-color:rgba(255,255,255,0)">It seems like there is a functional program bug since GDAL reads non-default units OK but cannot produce a geotiff with seconds. Something similar will probable happen with UTM PCSs if the target is feet instead of meter (the default).</span></font></div>
<div><font><span style="background-color:rgba(255,255,255,0)"><br></span></font></div><div><font><span style>I would like to know if somebody else found this problem and how he solved it.</span></font></div><div><font><span style>Of course, I wait to get confirmation it is a bug or maybe I did something wrong.</span></font></div>
<div><font><span style>If it is a "new bug" I woułd like to know what can I do to advance the solution for it.</span></font></div><div><font><span style><br></span></font></div><div><font><span style>This bug reminds me of a very old ArcMap bug that manifests itself when you open a polygon or a raster with non-default units like second. If this is the first item you open, the units are not read correctly until you reload the GCS into the Data Frame. And than you get Degrees...the logic (or the lack of it) is similar.</span></font></div>
<div><font><span style><br></span></font></div><div><font><span style>Best regards,</span></font></div><div><font><span style>Moty</span></font></div>