TIF Worldfile auslesen

Hallo morgen1,
zu 4 der 6 Zeilen kann ich dir ein bisschen was sagen.

Zeile 1 beschreibt die "Auflösung" oder den "Maßstab". Beide Begriffe passen nicht wirklich (deshalb die Gänsefüßchen). Die Zahl gibt an wieviele "Gittereinheiten" (je nach Koordinatensystem zb Meter oder Grad) einem Pixel des Kartenbilds (TIF oder was auch immer) auf der x-Achse entsprechen (also Rechtswert oder Längengrad pro Pixel).
Zeilen 2 und 3 habe ich bisher noch nicht so richtig kapiert. Hier geht es wohl irgendwie um eine eventuelle Rotation des Gitters.
Zeile 4 beschreibt (analog zu Zeile 1) die Auflösung entlang der y-Achse (Hochwert/Breitengrad pro Pixel). Auf der nördlichen Halbkugel ist der Wert negativ.
Zeile 5: UTM-Rechtswert oder Längengrad (in Dezimalgrad) des linken Bildrandes (bin mir jetzt nicht sicher ob der Wert bei westlichen Längengraden ebenfalls negativ sein muss oder ob das dann in Zeile 1 passiert).
Zeile 6: UTM-Hochwert oder Breitengrad des oberen Bildrandes.

Aus diesen Angaben lassen sich dann die Koordinaten der 4 Eckpunkte errechnen. Die Größe des Kartenausschnitts in Pixel muss natürlich auch bekannt sein.

Soweit mein leienhaftes Halbwissen. Hoffe du kommst damit schon mal ein Stückchen weiter.

Zitat von morgen1

Auch für GlobalMapper wäre es nützlich, wenn ich das endlich wüßte und dort ggf. die Projektionsangaben mit Hilfe des TFW berichtigen kann.
danke für Erklärungen

Erklärung des Worldfile-Inhalts und der Dateierweiterung hier, hier und hier. Worldfiles enthalten keinerlei Angaben zum verwendeten Referenzsystem und sind ohne begleitende Angabe von Projektion, Datum und Koordinatensystem (z.B. als PRJ-Datei oder als EPSG-Code) nutzlos, du kannst damit keine Projektionsangaben berichtigen.

Wenn du Geotiff-Dateien mit vollständigem Geotiff-Header als Quelle hast, schreibe vor der Bearbeitung mit Global Mapper Worldfile und PRJ-Datei raus. In der Regel geht die Georeferenzierung bei der Bearbeitung mit Grafikprogrammen verloren, falls du dabei auch die Bildgröße änderst, ist die gesamte Referenzierung hinfällig, weil pixelbezogen.

Was willst du noch mit Geotiff4, wenn du Global Mapper hast?

Zitat von morgen1

Jetzt erwarte ich, dass an der Landesgrenze die Teile, welche ausserhalb Östereichs liegen mit Null gefüllt sein müßten. Da werde ich wohl selbst ein miniprogramm schreiben müssen und die DEM-Werte für den Bereich ausserhalb der Grenze aus der Nasa einlesen müssen.

Such dir eine Shape-Datei mit dem Verlauf der österreichischen Grenze im Netz und nimm die in Global Mapper zum Maskieren, dann kannst du das auch automatisiert erledigen.

Zitat von morgen1

Die nächste Frage ist dann, ob die von GM aus den 10 m DEM erzeugten 1 Sec hgt besser sind als die der NASA ?

Angesichts der Lücken und Fehler in den SRTM-Daten würde ich das erwarten.

Zitat von morgen1

Wie könnte man das prüfen ?

Realistisch mit Hausmitteln gar nicht. Du kannst nur das Ergebnis deiner Konvertierung mit den SRTM-Daten überlagern und nachsehen, wo es auffällige Abweichungen gibt.

Die originalen SRTM-Daten verwende ich schon lange nicht mehr, sondern meistens die bearbeiteten Versionen von Jonathan de Ferranti. Wenn du viel Zeit und viel Rechenpower zur Verfügung hast, kannst du dich auch mal mit dem EU-DEM beschäftigen.

Zitat von morgen1

Ich habe mir ein Windows-Tool geschrieben, welches Leerstellen in HGT-Files auffüllen kann mit Daten aus anderen HGT Files.

Wäre m.M.n. nicht nötig gewesen. Weiter oben hatte ich schon einen einfacheren Weg aufgezeigt, es geht sogar noch einfacher:

1. Lade zuerst die SRTM-Daten der Fläche, in der deine österreichischen Höhendaten liegen -> erster Layer

2. Lade dann die österreichischen Höhendaten einfach dazu -> zweiter Layer

Der Stützpunktabstand und das Datenformat sind dabei völlig egal, solange Global Mapper es lesen kann.

3. Exportiere das Ganze dann in HGT-Dateien oder Geotiffs oder was auch immer du brauchst, Global Mapper exportiert für jeden Punkt den sichtbaren Layer, also auf der Fläche der österreichischen Höhendaten diese, auf den umgebenden Flächen die SRTM-Daten, wenn diese Flächen auf dem östereichischen Layer transparent gesetzt sind.

GM kann mit SRTM3-Basisdaten auch Pseudo-1"-HGTs erzeugen, die Zwischenwerte werden dann iirc linear interpoliert.