30

Hi,

I just want to inform all who was following this thread that I succeeded to make contour from Aster v2 and Cgiar too.
Cgiar is the far best source of elevation data at the moment ( fixed all anomalies ).

I want to publish here perl script I did for data conversion, maybe they can be useful for someboby, maybe somebody can improve them.

Perl script for converting Aster dem.tif to hgt:

#!/usr/bin/perl -w

while (< *.tif>) {

$tile = $_;
$tile =~ s/dem.tif//;
$tile =~ s/ASTGTM2//;
$tile =~ m/([SN])(\d+)([EW])(\d+)/;

`gdalwarp -of GTiff -srcnodata 32767  -rcs -order 3 -ts 1201 1201 -multi ASTGTM2_${tile}_dem.tif ${tile}.tif` ;

`gdal_translate -of SRTMHGT ${tile}.tif ${tile}.hgt` ;

`rm ${tile}.tif ${tile}.hgt.aux.xml` ;

}

Perl script for converting cgiar 1 degree tiles to hgt:

#!/usr/bin/perl -w

while (< *.zip>) {
$file_asc = $_;
$file_asc=~ s/.zip//;

unless (-e “${file_asc}.asc”) {unzip $_};
open(READ, “${file_asc}.asc” ) or die $! ;
@lines = ;
close (READ);

  $x = $lines[2];
  $y = $lines[3];

  $x =~ s/xllcorner     //;
  $y =~ s/yllcorner     //;

if ($x > 0 ) {
$lat = “N”;
} else {
$lat = “S”;
}

if ($y > 0 ) {
$lon = “E”;
} else {
$lon = “W”;
}

$rx = sprintf(“%.0f”, $x);
$ry = sprintf(“%.0f”, $y);

$arx = abs $rx ;
$ary = abs $ry ;

$tile = “$lat” . “$ary” . “$lon” . “0” . “$arx” ;

`gdalwarp -of GTiff -rcs -order 3 -ts 1201 1201 -multi ${file_asc}.asc ${tile}.tif` ;

`gdal_translate -of SRTMHGT ${tile}.tif ${tile}.hgt` ;

rm ${file_asc}.asc ;

`rm ${tile}.tif ` ;
`rm ${tile}.hgt.aux.xml ` ;

}