<div dir="ltr">Please disregard my previous message with the gdalinfo output. Here is the correct one:<br><br><font face="monospace">$ gdalinfo "dem.tif"<br>Driver: GTiff/GeoTIFF<br>Files: dem.tif<br>Size is 1541, 1335<br>Coordinate System is:<br>GEOGCRS["WGS 84",<br>    ENSEMBLE["World Geodetic System 1984 ensemble",<br>        MEMBER["World Geodetic System 1984 (Transit)"],<br>        MEMBER["World Geodetic System 1984 (G730)"],<br>        MEMBER["World Geodetic System 1984 (G873)"],<br>        MEMBER["World Geodetic System 1984 (G1150)"],<br>        MEMBER["World Geodetic System 1984 (G1674)"],<br>        MEMBER["World Geodetic System 1984 (G1762)"],<br>        MEMBER["World Geodetic System 1984 (G2139)"],<br>        MEMBER["World Geodetic System 1984 (G2296)"],<br>        ELLIPSOID["WGS 84",6378137,298.257223563,<br>            LENGTHUNIT["metre",1]],<br>        ENSEMBLEACCURACY[2.0]],<br>    PRIMEM["Greenwich",0,<br>        ANGLEUNIT["degree",0.0174532925199433]],<br>    CS[ellipsoidal,2],<br>        AXIS["geodetic latitude (Lat)",north,<br>            ORDER[1],<br>            ANGLEUNIT["degree",0.0174532925199433]],<br>        AXIS["geodetic longitude (Lon)",east,<br>            ORDER[2],<br>            ANGLEUNIT["degree",0.0174532925199433]],<br>    USAGE[<br>        SCOPE["Horizontal component of 3D system."],<br>        AREA["World."],<br>        BBOX[-90,-180,90,180]],<br>    ID["EPSG",4326]]<br>Data axis to CRS axis mapping: 2,1<br>Origin = (76.738472222222228,43.403749999999995)<br>Pixel Size = (0.000277777777778,-0.000277777777778)<br>Metadata:<br>  AREA_OR_POINT=Area<br>Image Structure Metadata:<br>  COMPRESSION=DEFLATE<br>  INTERLEAVE=BAND<br>Corner Coordinates:<br>Upper Left  (  76.7384722,  43.4037500) ( 76d44'18.50"E, 43d24'13.50"N)<br>Lower Left  (  76.7384722,  43.0329167) ( 76d44'18.50"E, 43d 1'58.50"N)<br>Upper Right (  77.1665278,  43.4037500) ( 77d 9'59.50"E, 43d24'13.50"N)<br>Lower Right (  77.1665278,  43.0329167) ( 77d 9'59.50"E, 43d 1'58.50"N)<br>Center      (  76.9525000,  43.2183333) ( 76d57' 9.00"E, 43d13' 6.00"N)<br>Band 1 Block=1541x2 Type=Int16, ColorInterp=Gray<br>  NoData Value=32767</font><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Sep 1, 2024 at 12:58 AM Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><u></u>

  
    
  
  <div>
    <p>What is the gdalinfo output on dem.tif ?<br>
    </p>
    <br>
    <blockquote type="cite">
      <div dir="ltr">
        <div><span style="box-sizing:border-box;color:rgb(31,35,40);font-size:14px">Is this
            a bug? It appears that the data type change might be
            influenced by the nodata value.<font face="monospace"><br>
            </font></span></div>
      </div>
    </blockquote>
    Yes, there is some logic to take the "union" of the largest data
    types involved in the warping (but it should probably ignore nodata
    values that are outside of the validity range of the input data
    type. I didn't check what it does currently)<br>
    <blockquote type="cite">
      <div dir="ltr">
        <div><span style="box-sizing:border-box;color:rgb(31,35,40);font-size:14px"><br>
          </span></div>
      </div>
      <br>
      <fieldset></fieldset>
      <pre>_______________________________________________
gdal-dev mailing list
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a>
</pre>
    </blockquote>
    <pre cols="72">-- 
<a href="http://www.spatialys.com" target="_blank">http://www.spatialys.com</a>
My software is free, but my time generally not.</pre>
  </div>

</blockquote></div>