[mkgmap-dev] DEM-file generation from geo-tiff sources ?

Sun Jan 7 07:59:14 GMT 2018

I don't know why my "raw-text" blocks got eaten by the nabble.com web
interface

franco_bez wrote
> I got the following warnings

 gdalwarp $GDAL_TRANSLATE_OPT -te -24.0001389 26.9998611 -22.9998611
28.0001389 ../all.vrt N27W024.tif
Creating output file that is 3601P x 3601L.
Warning 6: Driver GTiff does not support "COMPRESS creation option
Warning 6: Driver GTiff does not support "TILED creation option
Warning 6: Driver GTiff does not support "PREDICTOR creation option

> Here's my Ubuntu-script

use strict;

my $lonleft = -24;
my $lonright = 50;
my $latbottom = 27;
my $lattop = 64;

my $srcdir = "..";

my $delta = 1 / 3600 / 2;     # half pixel

print "#!/bin/sh\n";
#print "GDAL_TRANSLATE_OPT='-multi -wo NUM_THREADS=ALL_CPUS -t_srs EPSG:4326
-r near -ot Int16 -dstnodata -32768 -of GTiff -co \"COMPRESS=LZW\" -co
\"TILED=YES\" -co \"PREDICTOR=2\" -ts 3601 3601'\n";
print "GDAL_TRANSLATE_OPT='-multi -wo NUM_THREADS=ALL_CPUS -t_srs EPSG:4326
-r near -ot Int16 -dstnodata -32768 -of GTiff -ts 3601 3601'\n";
for (my $lat = $latbottom; $lat != $lattop; $lat++) {
    for (my $lon = $lonleft; $lon != $lonright; $lon++) {
       printf "gdalwarp \$GDAL_TRANSLATE_OPT -te %.7f %.7f %.7f %.7f
$srcdir/all.vrt %s%02d%s%03d.tif\n",
              $lon - $delta,
              $lat - $delta,
              $lon + 1 + $delta,
              $lat + 1 + $delta,
              'N',
              $lat >= 0 ? $lat : -$lat,
              $lon >= 0 ? 'E' : 'W',
              $lon >= 0 ? $lon : -$lon;
    }
}

--
Sent from: http://gis.19327.n8.nabble.com/Mkgmap-Development-f5324443.html