<div dir="ltr"><div class="gmail_default" style="font-family:monospace,monospace">Hi all,</div><div class="gmail_default" style="font-family:monospace,monospace"><br></div><div class="gmail_default" style="font-family:monospace,monospace">I have been using the GDAL command line tools to pre-process various kinds of remote sensing imagery datasets in GeoTIFF format. This includes resampling 16 bit images (eg. from Landsat 8) to 8-bit and to combine bands into RBG composites. </div><div class="gmail_default" style="font-family:monospace,monospace"><br></div><div class="gmail_default" style="font-family:monospace,monospace">In GDAL 11.4 (and previous 11.x that I tried), this doesn't work any more:</div><div class="gmail_default" style="font-family:monospace,monospace"><br></div><div class="gmail_default" style=""><font face="monospace, monospace">gdal_translate -ot Byte -scale [sceneID]_Bn.tif </font><span style="font-family:monospace,monospace">[sceneID]_8bit_Bn.tif</span><br></div><div><br></div><div><div class="gmail_default" style="font-family:monospace,monospace">​I usually apply this to the already clipped RGB composites, but digging deeper, even for single-band images what seems to be happening is that the -scale option doesn't correctly calculates the source range:</div><div class="gmail_default" style="font-family:monospace,monospace"><br></div><div class="gmail_default" style=""><ul style=""><li style=""><font face="monospace, monospace">If I use the command above, or with -scale_1, the whole output file has a value 0 for each pixel</font></li><li style=""><font face="monospace, monospace">If I use the "-b 1 -scale", the whole output file has a value 255 at each pixel</font></li><li style=""><font face="monospace, monospace">If I use "-scale src_min src_max" with a manually inserted src_min and src_max, I get a reasonable result (though not necessarily identical to the GDAL result with automatically calculated min and max.</font></li></ul><div><font face="monospace, monospace"><br></font></div><div><span style="font-family:monospace,monospace">What is driving me absolutely bonkers, though, is that *occasionally* the command "-b 1 -scale" DOES produce a correct 8-bit file. It's not reproducible, though: if I delete all 8bit files and re-run the script over a whole set of 16 bit files, a whole different file may be resampled correctly. </span></div><div><span style="font-family:monospace,monospace"><br></span></div><div><span style="font-family:monospace,monospace">With GDAL 2.0.2, my old script (written when GDAL was at 1.8.x on my machine) seems to be working at the moment, but my main system is currently at 11.4, and I have a reason not to upgrade right away. Is there a way to make it work with 11.4? I could run gdalinfo first and extract min and max from the stats, and then feed this into gdal_translate, but I'd like to avoid this path if possible. A test file (a file subsetted with "</span><font face="monospace, monospace">gdal_translate -projwin... " from a whole Landsat 8 band) is here: <a href="http://www2.gi.alaska.edu/~cwaigl/images/LC80660142015170LGN00_B5_clip.tif">http://www2.gi.alaska.edu/~cwaigl/images/LC80660142015170LGN00_B5_clip.tif</a> (3.8 MB).</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">Thanks for your help,</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">Chris</font></div></div><br></div><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div style="color:rgb(0,0,0)"><font face="monospace, monospace"><span style="font-size:12.8px">Christine (Chris) Waigl - </span><a href="mailto:cwaigl@alaska.edu" style="font-size:12.8px;color:rgb(17,85,204)" target="_blank">cwaigl@alaska.edu</a><span style="font-size:12.8px"> -  </span><a value="+19074745483" style="font-size:12.8px;color:rgb(17,85,204)">+1-907-474-5483</a><span style="font-size:12.8px"> - Skype: cwaigl_work</span><br></font></div><div style="color:rgb(0,0,0)"><font face="monospace, monospace">Geophysical Institute, UAF, 903 Koyukuk Drive, Fairbanks, AK 99775-7320, USA</font></div></div></div></div></div>
</div>