dem_datum_shift_3param — DEM `geodesy` op

• Datenarten: keine → points (ein Operator, der allein durch seine Argumente bestimmt ist — er nimmt kein Bild und keine Daten entgegen)

• Aufruf: import fullseye as fs; fs.ledger.dem_datum_shift_3param(lat_deg, lon_deg, h_m, dx_m, dy_m, dz_m, a_from_m, f_from, a_to_m, f_to) (die Implementierung direkt: import demops; demops.dem_datum_shift_3param(lat_deg, lon_deg, h_m, dx_m, dy_m, dz_m, a_from_m, f_from, a_to_m, f_to); aus dem Register: opsdem.get("dem_datum_shift_3param"))

Verwendung

Wechsel des geodätischen Datums — geozentrische Translation mit 3 Parametern (points).

Selbst dasselbe "35 Grad 41 Minuten Nord" bezeichnet im alten japanischen geodätischen System und im Welt-Geodätischen System Punkte, die am Boden mehrere hundert Meter auseinanderliegen. Diese op führt diesen Wechsel als Translation in geozentrischen kartesischen Koordinaten aus (der Fall Molodensky-Badekas ohne Rotation und Maßstab)::

ECEF(a_from, f_from) + (dx, dy, dz) -> geodätisch(a_to, f_to)

★Keine Standardwerte für die Parameter. Von welchem Bezugssystem in welches, mit welcher Translation transformiert wird, unterscheidet sich nach Land, Region und Version; stillschweigende Annahmen liefern "Koordinaten ohne Ausnahme, aber um mehrere hundert Meter verschoben" —— das ist der Daseinszweck dieser op, daher muss der Aufrufer sie immer angeben.

Args:

lat_deg / lon_deg: Breite/Länge im Ausgangssystem [Grad]. Arrays erlaubt.

h_m: ellipsoidische Höhe im Ausgangssystem [m].

dx_m / dy_m / dz_m: geozentrische Translation [m] (Quelle → Ziel).

a_from_m / f_from: große Halbachse [m] und Abplattung des Ausgangsellipsoids.

a_to_m / f_to: große Halbachse [m] und Abplattung des Zielellipsoids.

Returns:

(n, 3) mit (Breite[Grad], Länge[Grad], ellipsoidische Höhe[m]).

Per geschlossener Form prüfbar: bei Translation 0 und gleichem Ellipsoid identische Transformation (die Hin-und-zurück-Untergrenze ist dieselbe wie bei dem_ecef_to_geodetic). Bei reiner Translation stimmt die Verschiebung am Boden mit den in der ECEF-Basis des Punktes zerlegten Komponenten überein.

ValueError: nicht-endliche Werte, Breite außerhalb ±90, nicht-positive große Halbachse, Abplattung außerhalb [0, 1).

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 `points` als Eingabe)

dem_ecef_to_geodetic · dem_geodetic_from_enu

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_geoid_height · dem_height_frame_convert

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

dem_datum_shift_3param — DEM geodesy op