Hallo Joern,
danke für die Blumen, ich bin einfach Deiner Ermutigung gefolgt und habe die Programme ausprobiert, die weoli erwähnt hatte. Ich habe ein paar Berechnungen angestellt zum Pixelabstand. Wäre schön, wenn Du kurz überprüfen könntest, ob ich mich da nicht verrechnet habe:
Ich habe einmal mit GPSMapEdit zwei Polylines erzeugt mit Kassel in der Mitte, eine Polyline in Nord-Süd Richtung und eine in Ost-West Richtung, jeweils mit einer Distanz von 5 Grad. Mit rechte Maustaste/Properties habe ich die Länge abgelesen und erhalte für die SRTM3 Daten folgende Werte:
347.300 Meter Distanz Ost-West / 6.000 Pixel -->57,8 Meter Pixelabstand in Ost-West Richtung
255.600 Meter Distanz Nord-Süd / 6.000 Pixel -->42,6 Meter Pixelabstand in Nord-Süd Richtung
347.300 Meter Distanz Ost-West / 5 Grad -->69.460 Meter Gradabstand in Ost-West Richtung
255.600 Meter Distanz Nord-Süd / 5 Grad -->51.120 Meter Gradabstand in Nord-Süd Richtung
Wenn ich die SRTM3 v4 Daten von CGIAR in MICRODEM einlese (Data Manipulation/Edit/DEM Header), zeigt der Data Header einen Pixelabstand von jeweils 0,000 833 354 000 ° Länge und Breite an. Daraus ergibt sich, daß die eingelesene Datei um 8,61 Meter (Längengrad) bzw. 6,34 Meter (Breitengrad) über die Größe von 5° mal 5° hinausragt.
Wenn ich die SRTM3 v4.1 Daten einlese, zeigt der Data Header einen Pixelabstand von jeweils 0,000 833 333 354 ° an. Daraus ergibt sich, daß die eingelesene Datei um 0,86 Zentimeter (Längengrad) bzw. 0,63 Zentimeter (Breitengrad) über die Größe von 5° mal 5° hinausragt.
Wenn ich nun den Pixelabstand über die Funktion Data Manipulation/Edit/DEM Header auf jeweils 0,000 833 333 333 ° ändere, dann erhalte ich 0,14 Millimeter (Längengrad) bzw. 0,10 Millimeter (Breitengrad) weniger als 5° mal 5°.
5 Grad / 6.000 Pixel = 0,000 833 333 333 333 ... Grad/Pixel
0,000 833 354 000 Grad/Pixel * 6.000 Pixel = 5,000 124 000 Grad --> 0,000 124 000 Grad zuviel
0,000 833 333 354 Grad/Pixel * 6.000 Pixel = 5,000 000 124 Grad --> 0,000 000 124 Grad zuviel
0,000 833 333 333 Grad/Pixel * 6.000 Pixel = 4,999 999 998 Grad --> 0,000 000 002 Grad zuwenig
0,000 124 000 Grad * 69.460 Meter / 1 Grad = 8,613 040 m
0,000 124 000 Grad * 51.120 Meter / 1 Grad = 6,338 880 m
0,000 000 124 Grad * 69.460 Meter / 1 Grad = 0,861 304 cm
0,000 000 124 Grad * 51.120 Meter / 1 Grad = 0,633 888 cm
-0,000 000 002 Grad * 69.460 Meter / 1 Grad = -0,138 920 mm
-0,000 000 002 Grad * 51.120 Meter / 1 Grad = -0,102 240 mm
Wenn ich die v4.1 Daten nun so nehme, wie sie sind (und ohne den Pixelabstand im Data Header zu ändern), dann weist das nordöstlichste Pixel in der Matrix einen Fehler von weniger als einem Zentimeter auf. (Bei den v4 Daten sind es immerhin mehr als 5 Meter und hier kommt noch die Problematik des südwestlichsten Pixels dazu).
Ergebnis: mit den v4.1 Daten bekomme ich jetzt genau 25 (1°mal1°) *.hgt Dateien aus 1 (5°mal5°) *.ASC Datei, wobei die südwestliche Ecke jeweils genau auf N*.000E*.000 liegt.
Ich wiederhole mich nur, wenn ich schreibe, daß mein vorrangiges Ziel war, auf schnelle und einfache Weise die v4 Daten von CGIAR für GPS-Track-Analyse.net verfügbar zu machen, um eine einheitliche Datenbasis für den Vergleich verschiedener Tracks zu haben.
Deine Ausführungen zu den Höhenmodellen sind im Moment zuviel für mich. Zu gegebener Zeit, werde ich aber dankbar darauf zurückkommen. Ich erinnere mich noch, mal aufgeschnappt zu haben, daß im Harz die Landkarten nicht richtig zusammenpaßten, als aus zwei deutschen Staaten wieder einer wurde. Und das wiederum erinnert mich an die verschiedenen Auffassungen und Unstimmigkeiten zur Datierung der Ägyptischen Geschichte. Und wenn man möchte, kann man dann auch noch über die kategorische Ablehnung/Ignoranz einiger Wissenschaftler gegenüber den "biochemische[n] Einwände[n] gegen die Evolutionstheorie" von Michael Behe den Kopf schütteln, aber das ist nun wirklich abseits des Themas.
Ich habe oft mit Kopfschütteln die zahlreichen Forderungen hier im Forum nach absolut genauen Höhendaten gelesen. Und Du hast da schon einige Beiträge zu geschrieben, wenn ich es recht überblicke.
Tipp: die v4.1 Daten im ArcInfoASCII Format gibt es vom FTP Server des JRC (Italien). Verknüpfung zur Internetseite von CGIAR-CSI. Die 872 Dateien haben eine Gesamtgröße von 14,1GB. Zur Beachtung: die Datein sind nur für den nichtkommerziellen Gebrauch. Wer die Daten kommerziell nutzen möchte, muss sich an Dr. Andy Jarvis (A.Jarvis(at)cgiar.org) wenden.
Gruß, Martin