<div dir="ltr"><br><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Apr 14, 2014 at 11:16 AM, Nicole Stoffels <span dir="ltr"><<a href="mailto:stoffels@f2e.de" target="_blank">stoffels@f2e.de</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Gdal-List,<br>
<br>
I have a large ASTER file (downloaded from <a href="http://www.eea.europa.eu/data-and-maps/data/eu-dem#tab-european-data" target="_blank">http://www.eea.europa.eu/data-<u></u>and-maps/data/eu-dem#tab-<u></u>european-data</a>) in the GeoTIFF-format. It is compressed with LZW. What I want to do is:<br>

<br>
1. unpack it<br>
2. seperate a smaller region from it<br>
3. convert this regional GeoTIFF to a PNG<br>
<br>
What I am doing:<br>
<br>
$ gdal_translate -co COMPRESS=NONE EUD_CP-DEMS_6500035000-AA.tif tmp1.tif<br>
<br>
$ gdal_translate -srcwin 6500 3300 1000 1000 tmp1.tif tmp2.tif<br>
<br>
$ gdal_translate -of PNG tmp2.tif tmp2.png<br>
<br>
When I display the info with gdalinfo the first step seems to produce a decompressed file. The line COMPRESSION=LZW is not there anymore. </blockquote><div><br></div><div>You need to set the compression options on any gdal command, it does not use the compression of the input file. Also, you specifically asked for no compression, so it's doing what you told it to.</div>
<div>use "-co COMPRESS=LZW"</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">However, if I want to display the min and max with gdalinfo, I get the following error.<br>

<br>
ERROR 1: Failed to compute min/max, no valid pixels found in sampling.<br></blockquote><div><br></div><div>try without changing the extents and see if a single command (as suggested by Jean-Luc) works</div><div><br></div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
I thought that the file might be too large, but if I change the data_type from float 32 to e.g. Int16 I get 0 as min and max value. That's why I assume that already the first gdal_translate does not work as I actually want it to work.<br>
</blockquote><div><br></div><div>your input values are probably between 0 and 1 in this case</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
The other two transformations work without error messages. However, in the end I get a black png.<br>
<br>
I have tried all types of transformation with various options. But as I am new in the subject and do not know much about gdal, I have no further ideas how to generate a file from my large GeoTIFF to continue working with.<br>

<br>
I hope someone has an idea, what's going wrong with my gdal_translations. Thanks very much in advance.<br>
<br>
Best regards<br>
Nicole<br>
______________________________<u></u>_________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/<u></u>mailman/listinfo/gdal-dev</a><br>
</blockquote></div><br></div></div>