Hallo,
das wäre mein Ziel. Diese liegen als GeoTiff vor zusätzlich gibt es noch eine .tlw und eine .hdr Datei zu jedem tif.
Wie ich die Höhenlinien aus den hgt-Dateien bekomme ist mir klar. Gibt es eine Möglichkeit der Konvertierung der CIGAR-Daten nach hgt oder eine direkte Methode zur Konvertierung in ein osm-Format?
Die geotiffs mit MICRODEM
http://www.usna.edu/Users/oceano/pgu…e/microdem.htm
ins BIL format konvertieren und dann mit bilxstrm ins hgt
http://www.cplus.org/rmw/bilxsrtm.zip
konvertieren.
Möglichst große Bereiche überlappend (so ca. 20 Kacheln gehen auf einmal > 32 Bitgrenze!!) da sonst an den Kachelgrenzen hellere Linien erscheinen!
Hallo Henning,
das Thema Höhenlinien möchte ich für meine Karte dieses Jahr auch noch umsetzen.
Insofern wäre ich an deinem Ansatz und Vorgehensmodell sehr interessiert.
Gruß Klaus
PS: Meine Vorstellungen gehen derzeit in die Richtung, die Höhenlinien direkt in die Karte zu integrieren. Hierfür gibt es ja speziell die Linientypen 0x20, 0x21 und 0x22. Derzeit sind mir allerdings die Vor- und Nachteile der verschiedenen Ansätze (integrierte Höhenlinien, Höhenlinien als zusätzlicher, transparenter Layer, …) noch nicht klar.
ich wollte es mal versuchen aber bei mir scheitert es schon daran, dass ich zu blöd bin, den Download dort auf der Seite zu finden.
Ich weis nicht wie es früher war.
Zur Zeit gehst du auf die Hauptseite von CGIAR (nicht C_IG_AR bitte) und klickst dort auf den Link zu SRTM v4.1 (unter Quick Links). Dort findest du erst einmal viel Text und weit nach der Hälfte der Seite den Abschnitt Download. Und welche Überraschung 
dort findest du mehre Links zu verschiedenen Download-Seiten.
Die CGIAR Download-Seite benutzt zwar eine sehr kleine Weltkarte aber bei dem 5°x5° Raster funktioniert das noch.
Kurz vor dem Download-Abschnitt gibt es einen Abschnitt über alternative Quellen, der noch interessant sein könnte.
PS @aighes: Kannst du den Namen dieses Threads bitte korrigieren.
Edbert (EvanE)
Ich hab es versucht umzusetzen, allerdings ist es mit MICRODEM nicht wirklich praktikabel größere Gebiete umzusetzen. Von daher hab ich es erstmal mit diesem Weg nicht weiter probiert.
Ich sitze momentan einem Weg über imagemagick, geotiffcp und DEM2Topo. Wenn ich da weiter bin, melde ich mich.
Srtm dürfte aber deutlich weniger Widerstand bieten als die CGIAR-Daten. 
FTP: ftp://xftp.jrc.it/pub/srtmV4/tiff/
mit GUI: http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp
Danke für die Einladung, aighes.
Also ich mach das eigentlich schon seit Jahren… wundert mich, dass da sonst noch niemand drauf gekommen ist.
Und wenn man ganz penibel ist, schreibt man sich noch ein kleines Progrämmchen, was die .mp Dateien durchforstet und die überstehenden Randartefakte von DEM2Topo wegschneidet. Ist etwa ein halber Pixel GeoTIFF an jeder Kante, wo die Höhenlinien der benachbarten tiles sich überlappen. Aber das führt hier wohl etwas zu weit…
#!/usr/bin/perl
# Cut CGIAR SRTM data into smaller tiles of e.g. 1x1 degree
# usage:
# cgiar2tile.pl [out.tif] [top] [left] [bottom] [right]
#
# Just a quick and dirty hack that may very well not work for you. If it does,
# there should be a one pixel boundary left all around the output image, so that
# contour line extractors can work properly as well. The geo-registration
# matches the current SRTM 4.1 release and is correct, if used for whole degree
# tiles.
# e.g.:
# cgiar2tile.pl 5101.tif 52 -2 51 -1
#
# Currently, the CGIAR data is expected to reside in the current working
# directory upon execution of the script. If a tile is missing, you will be
# told about that.
#
# Please note the two variables $verbose and $execute, that allow to turn
# on/off the messages and performing a dry run without actually executing the
# commands.
#
#
# written by Olaf Kaehler
#
# original versions date back to 2006
# latest significant changes 2009-09-22
# first released version 2011-10-11
$im_convert="convert -define quantum:format=signed";
$gt_geotifcp="geotifcp";
$verbose=100;
$execute=1;
$outfilename = "out.tif";
if ( $#ARGV == 4 ) {
$outfilename = shift;
}
$top=shift;
$left=shift;
$bottom=shift;
$right=shift;
sub getTileFilename($;$) {
my $tile_x=shift;
my $tile_y=shift;
my $tile_0x = $tile_x;
my $tile_0y = $tile_y;
$tile_0x=~s/^(.)$/0$1/;
$tile_0y=~s/^(.)$/0$1/;
# version 3/4 files
my $filename = "srtm_".$tile_0x."_".$tile_0y.".tif";
if ( -r $filename ) { return $filename; }
# version 2 files
$filename = "z_".$tile_x."_".$tile_y.".tif";
if ( -r $filename ) { return $filename; }
# version 1 files:
$filename = "s_".$tile_0x."_".$tile_0y.".tif";
if ( -r $filename ) { return $filename; }
printf "tile $tile_x $tile_y does not exist...\n";
return $filename;
}
sub pos2pixel($;$) {
my $lat=shift;
my $lon=shift;
$lat = 60-$lat;
$lon = 180+$lon;
$lat = $lat*1200;
$lon = $lon*1200;
return $lat,$lon;
}
sub pixel2tile($;$) {
my $pix_x = shift;
my $pix_y = shift;
$tile_x = $pix_x/6000;
$tile_x =~ s/\..*//;
$pix_x = $pix_x - $tile_x*6000;
$tile_y = $pix_y/6000;
$tile_y =~ s/\..*//;
$pix_y = $pix_y - $tile_y*6000;
# tiles start at 1,1
$tile_x+=1;
$tile_y+=1;
return ($tile_x,$tile_y,$pix_x,$pix_y);
}
sub createTIF($;$;$;$) {
$top=shift;
$left=shift;
$bottom=shift;
$right=shift;
# get pixels on global SRTM data image
($gpix1_y,$gpix1_x) = pos2pixel($top,$left);
($gpix2_y,$gpix2_x) = pos2pixel($bottom,$right);
# we usually want the pixel *center* coordinates inside the given boundaries
$gpix1_x = $gpix1_x-1;
$gpix1_y = $gpix1_y-1;
$gpix2_x = $gpix2_x+1;
$gpix2_y = $gpix2_y+1;
if ($verbose>=3) {
print "extracting pixels from $gpix1_x,$gpix1_y to $gpix2_x,$gpix2_y\n";
}
# get tiles and local coordinates within the tiles
($tile1_x,$tile1_y,$lpix1_x,$lpix1_y) = pixel2tile($gpix1_x,$gpix1_y);
($tile2_x,$tile2_y,$lpix2_x,$lpix2_y) = pixel2tile($gpix2_x,$gpix1_y);
($tile3_x,$tile3_y,$lpix3_x,$lpix3_y) = pixel2tile($gpix1_x,$gpix2_y);
($tile4_x,$tile4_y,$lpix4_x,$lpix4_y) = pixel2tile($gpix2_x,$gpix2_y);
if ($verbose>=3) {
print " left top is on tile $tile1_x,$tile1_y at $lpix1_x,$lpix1_y\n";
print "(right top is on tile $tile2_x,$tile2_y at $lpix2_x,$lpix2_y)\n";
print "(left bottom is on tile $tile3_x,$tile3_y at $lpix3_x,$lpix3_y)\n";
print " right bottom is on tile $tile4_x,$tile4_y at $lpix4_x,$lpix4_y\n";
}
# crop the necessary patches from all tiles
my %patches;
$x_page=0;
foreach $tile_x ($tile1_x..$tile4_x) {
$x_start=0;
$x_end=5999;
if ($tile_x==$tile1_x) { $x_start=$lpix1_x; }
if ($tile_x==$tile4_x) { $x_end=$lpix4_x; }
$x_size = $x_end - $x_start + 1;
$y_page=0;
foreach $tile_y ($tile1_y..$tile4_y) {
$y_start=0;
$y_end=5999;
if ($tile_y==$tile1_y) { $y_start=$lpix1_y; }
if ($tile_y==$tile4_y) { $y_end=$lpix4_y; }
$y_size = $y_end - $y_start + 1;
$s_img=getTileFilename($tile_x,$tile_y);
if ($verbose>=2) {
print "from $s_img: ($x_start,$y_start) - ($x_end,$y_end)\n";
}
$patchname="tmp_c".scalar(keys %patches).".tif";
$cmdline = $im_convert." ".
$s_img." ".
"-crop ".$x_size."x".$y_size."+".$x_start."+".$y_start." ".
$patchname;
my $retval=0;
if ($verbose>=1) { print " $cmdline\n"; }
if ($execute>=1) { system $cmdline; $retval=$?; }
print "success: $retval\n";
if ($retval == 0) { $patches{$patchname} = $x_size."x".$y_size."+".$x_page."+".$y_page; }
$y_page += $y_size;
}
$x_page += $x_size;
}
# create mosaic of all the small patches
$cmdline = "$im_convert -page ".$x_page."x".$y_page."+0+0";
foreach $i (keys %patches) {
$cmdline = $cmdline." -page ".$patches{$i}." $i";
}
$cmdline = $cmdline." -mosaic tmp_c.tif";
if ($verbose>=1) { print " $cmdline\n"; }
if ($execute>=1) { system $cmdline; }
}
sub createGEO($;$) {
$top=shift;
$left=shift;
if ($execute>=1) { open($FILE, ">tmp_c.geo"); }
else { $FILE=stdout; }
print $FILE
"Geotiff_Information:\n".
" Version: 1\n".
" Key_Revision: 1.0\n".
" Tagged_Information:\n".
" ModelTiepointTag (2,3):\n".
" 0 0 0 \n".
#" ".($left-0.0004166)." ".($top+0.0004166)." 0 \n".
#" ".$left." ".$top." 0 \n".
" ".($left-0.000833333333)." ".($top+0.000833333333)." 0 \n".
" ModelPixelScaleTag (1,3):\n".
" 0.000833333333 0.000833333333 0 \n".
" End_Of_Tags.\n".
" Keyed_Information:\n".
" GTModelTypeGeoKey (Short,1): ModelTypeGeographic\n".
" GTRasterTypeGeoKey (Short,1): RasterPixelIsArea\n".
' GTCitationGeoKey (Ascii,218): "IMAGINE GeoTIFF Support\nCopyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved\n@(#)$RCSfile: egtf.c $ $Revision: 1.9.2.9 $ $Date: 2002/10/28 15:14:13EST $\nProjection Name = GCS_WGS_1984\nUnits = dd\nGeoTIFF Units = dd"'."\n".
" GeographicTypeGeoKey (Short,1): GCS_WGS_84\n".
" GeogAngularUnitsGeoKey (Short,1): Angular_Degree\n".
" End_Of_Keys.\n".
" End_Of_Geotiff.\n";
$cmdline = $gt_geotifcp." -g tmp_c.geo -c none tmp_c.tif $outfilename";
if ($verbose>=1) { print " $cmdline\n"; }
if ($execute>=1) { system $cmdline; }
}
if ($verbose>=2) {
print "extracting area from $top,$left to $bottom,$right\n";
}
createTIF($top,$left,$bottom,$right);
createGEO($top,$left);
Hallo Olaf,
könntest du deine finalen .mp-Dateien für einen Test zum Download bereit stellen?
Gruß Klaus
Hab gerade für einen Bekannten einen Satz Höhenlinien der Türkei erstellt, ca. 200MB als .img Dateien: http://island-olaf.de/tmp/turkey_contour.zip
Die .mp Dateien sind ja nur Zwischenprodukt und für den gesamten Türkei-Datensatz etwa 2 GB groß, ohne dabei zusätzliche Informationen zu enthalten. Aber wenn du unbedingt ein Beispiel haben willst: http://island-olaf.de/tmp/91314129.mp (ca 2.6MB für eine 1x1 Kachel)
Viel Spaß damit,
Olaf
(Natürlich könnte man auch andere Höhenlinien basteln, 20m Abstände, was auch immer. Muss man nur bei DEM2Topo so einstellen, wie man’s braucht)
Danke erstmal - ich werde dann mit Island testen.
Gruß Klaus (JAPH)
Viel Spaß beim Testen, aber Island wirst du in den SRTM Daten nicht finden. Vielleicht hilft dir http://www.viewfinderpanoramas.org/dem3.html oder vielleicht nimmst du auch gleich Tom’s Island-Karte http://ourfootprints.de/gps/mapsource-island.html . Der hat genau diese Arbeit nämlich schon gemacht allerdings mit cGPSmapper statt mit mkgmap.
Hallo,
allen vorangegangenen Beiträgen mein Dank!
Soeben habe ich völlig problemlos und auf Anhieb:
Grüsse,
Federico