dem_geoid_height — DEM `geodesy` op

• Datenarten: depth → signal

• Aufruf: import fullseye as fs; fs.ledger.dem_geoid_height(geoid, lat_deg, lon_deg, lat0_deg, lon0_deg, d_lat_deg, d_lon_deg) (die Implementierung direkt: import demops; demops.dem_geoid_height(geoid, lat_deg, lon_deg, lat0_deg, lon0_deg, d_lat_deg, d_lon_deg); aus dem Register: opsdem.get("dem_geoid_height"))

Verwendung

Wertet ein Geoidhöhengitter an Breite/Länge durch bilineare Interpolation aus und gibt N [m] zurück (signal).

GNSS liefert die ellipsoidische Höhe h, Karten und Konstruktionspläne verwenden die orthometrische Höhe H. Die Differenz ist die Geoidhöhe N (H = h - N), in der Umgebung Japans etwa +30 bis +40 m. Dieser op ist ein Werkzeug zum Lesen veröffentlichter Geoidmodell-Gitter (äquidistante Gitter wie die GSIGEO- / EGM-Familie) und enthält das Modell selbst nicht —— welche Version verwendet wurde, liegt in der Verantwortung des Aufrufers, und mit einer anderen Version verschieben sich die Höhen um einige cm.

Args:

geoid: (H, W)-Geoidhöhe [m]. geoid[0, 0] ist

(lat0_deg, lon0_deg); mit steigender Zeile wächst die Breite um d_lat_deg,

mit steigender Spalte die Länge um d_lon_deg (beide Vorzeichen zulässig).

lat_deg / lon_deg: Breite/Länge der abzufragenden Punkte [Grad]. Arrays erlaubt.

lat0_deg / lon0_deg: Ursprung des Gitters [Grad].

d_lat_deg / d_lon_deg: Schrittweite des Gitters [Grad]. 0 wird abgelehnt.

Returns:

(n,)-Geoidhöhen [m].

Mit einer geschlossenen Formel prüfbar: für ein Gitter, in dem N eine lineare Funktion von Breite und Länge ist, ist die bilineare Interpolation exakt (Maschinengenauigkeit). Bei quadratischen Termen ist der Fehler durch |∂²N| * d² / 8 nach oben beschränkt.

ValueError: Gitter nicht 2-D / kleiner als 2x2 / nicht endlich, Schrittweite 0 oder nicht endlich, Breite außerhalb ±90, und Punkte außerhalb des Gitters (Meeresgebiete oder Bereiche außerhalb werden nicht stillschweigend mit Randwerten gefüllt —— eine extrapolierte Geoidhöhe ist kein Messwert).

Ausführlicher Anwendungsleitfaden

• Leitfaden zur Familie dem_terrain_analysis

Referenzen (Beispieldaten, Literatur)

• Katalog der Beispieldaten (Download-URLs / Lizenzen) — 2-D nutzt skimage.data (BSD/Public Domain) plus synthetische Bilder, 3-D nennt Download-URLs echter Datenquellen (Stanford, PDS, …).

• Herkunft und Literatur der Operatoren — die Quellen der Forschung/Verfahren, auf denen diese Operatorfamilie beruht.

• Der kanonische Algorithmus (Autor, Jahr) und seine Anwendungen stehen im Familienleitfaden oben.

Ausführbare Beispiele (verifizierte Samples, die diesen Operator wirklich aufrufen)

• poc_geodetic_benchmarks_real — py -3.11 examples/poc_geodetic_benchmarks_real.py

• poc_geodetic_height_frames — py -3.11 examples/poc_geodetic_height_frames.py

Typkompatible Folge-Operatoren (nehmen `signal` als Eingabe)

dem_height_frame_convert · dem_height_frame_residual

Gleiche Kategorie (`geodesy`)

dem_geodetic_to_ecef · dem_ecef_to_geodetic · dem_geocentric_grid · dem_earth_curvature_drop · dem_cell_size_webmercator · dem_geodetic_slope · dem_height_frame_convert · dem_height_frame_residual

*Provenance: demops.py — DEM Operator-Registry. Diese Notiz wird von tools/opdocs.py md erzeugt (nicht von Hand bearbeiten).*

dem_geoid_height — DEM geodesy op