<div dir="ltr">Sorry for bringing this back up - I'm having problems converting tiff to .hgt. - how can I run a script that converts all geotiff (not 25m resolution but 30m) in folder /home/example/ to .hgt?<div>Basically I want to fill the voids in SRTM1v30 and need to convert .hgt to tiff as gdal and then get them back to hgt.</div><div><br></div><div>so first I simply run:</div><div><br></div><div>
<span style="font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-indent:0px;text-transform:none;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial;color:rgb(36,39,41);font-family:Consolas,Menlo,Monaco,"Lucida Console","Liberation Mono","DejaVu Sans Mono","Bitstream Vera Sans Mono","Courier New",monospace,sans-serif;font-size:13px;text-align:left;white-space:pre-wrap;background-color:rgb(239,240,241);float:none;display:inline"><br class="gmail-Apple-interchange-newline">for hgtfile in *.hgt;do gdal_fillnodata.py -md 500 $hgtfile $hgtfile.tif; done</span><span style="color:rgb(0,0,0);font-family:arial,sans-serif;font-size:12.8px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><span> </span></span>
<br></div><div><br></div><div>which will voildfill all .hgt files but gdal created rubbish (at least if I feed it to phyghtmap) if I try to simply voildfill hgt and have output also as hgt. So I need to convert them back to hgt before I can hand them over to pyhghtmap.</div><div>It should work vs a folder of the whole world 60N/60S.</div><div><br></div><div><br></div><div>Best in bash if possible, but perl is also fine if I understand the steps (I downloaded the SRTM1 from earthexplorer as .hgt - or better I will have them in a couple of days, download is pretty slow even though it is only 85GB zipped...)</div></div><div class="gmail_extra"><br><div class="gmail_quote">On 7 January 2018 at 13:25, franco_bez <span dir="ltr"><<a href="mailto:franco.bez@web.de" target="_blank">franco.bez@web.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Frank Stinner-2 wrote<br>
<span class="">> Hi Franco,<br>
> <br>
> i have seen an silly error in my area. Europa is a little bit greater ;)<br>
> Just i have "$lattop = 70". The split create 3608 tif's with 4,5 GB in 85<br>
> min. gdal4hgt create 1171 zipped hgt's with 5,8 GB in 3 h :((<br>
> And then i rebuild the tif's from the hgt's in 20 min. That are 4 GB.<br>
<br>
</span>OK<br>
<span class=""><br>
> Perhaps you have an to old gdal-version without lzw-compression? Or<br>
> perhaps<br>
<br>
</span>It's the current version of gdal<br>
<span class=""><br>
> -co \"COMPRESS=LZW\" -co \"TILED=YES\" -co \"PREDICTOR=2\"<br>
> <br>
> should be better<br>
> <br>
> -co COMPRESS=LZW -co TILED=YES -co PREDICTOR=2<br>
<br>
</span>That did the trick :-) <br>
<br>
Thanks a lot.<br>
<br>
now this is the new perl script for Ubuntu<br>
---<br>
<span class=""><br>
use strict;<br>
<br>
my $lonleft = -24;<br>
my $lonright = 50;<br>
my $latbottom = 27;<br>
</span>my $lattop = 70;<br>
<span class=""><br>
my $srcdir = "..";<br>
<br>
my $delta = 1 / 3600 / 2; # half pixel<br>
<br>
print "#!/bin/sh\n";<br>
</span><span class="">print "GDAL_TRANSLATE_OPT='-multi -wo NUM_THREADS=ALL_CPUS -t_srs EPSG:4326<br>
</span>-r near -ot Int16 -dstnodata -32768 -of GTiff -co COMPRESS=LZW -co TILED=YES<br>
-co PREDICTOR=2 -ts 3601 3601'\n";<br>
<span class="">for (my $lat = $latbottom; $lat != $lattop; $lat++) {<br>
for (my $lon = $lonleft; $lon != $lonright; $lon++) {<br>
printf "gdalwarp \$GDAL_TRANSLATE_OPT -te %.7f %.7f %.7f %.7f<br>
$srcdir/all.vrt %s%02d%s%03d.tif\n",<br>
$lon - $delta,<br>
$lat - $delta,<br>
$lon + 1 + $delta,<br>
$lat + 1 + $delta,<br>
'N',<br>
$lat >= 0 ? $lat : -$lat,<br>
$lon >= 0 ? 'E' : 'W',<br>
$lon >= 0 ? $lon : -$lon;<br>
}<br>
}<br>
<br>
</span>---<br>
<br>
Ciao,<br>
Franco<br>
<span class="im HOEnZb"><br>
<br>
<br>
--<br>
Sent from: <a href="http://gis.19327.n8.nabble.com/Mkgmap-Development-f5324443.html" rel="noreferrer" target="_blank">http://gis.19327.n8.nabble.<wbr>com/Mkgmap-Development-<wbr>f5324443.html</a><br>
</span><div class="HOEnZb"><div class="h5">______________________________<wbr>_________________<br>
mkgmap-dev mailing list<br>
<a href="mailto:mkgmap-dev@lists.mkgmap.org.uk">mkgmap-dev@lists.mkgmap.org.uk</a><br>
<a href="http://www.mkgmap.org.uk/mailman/listinfo/mkgmap-dev" rel="noreferrer" target="_blank">http://www.mkgmap.org.uk/<wbr>mailman/listinfo/mkgmap-dev</a><br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div><div>Felix Hartman - Openmtbmap.org & VeloMap.org<br></div>Schusterbergweg 32/8<br></div><div>6020 Innsbruck<br></div></div>Austria - Österreich</div></div></div></div>
</div>