dem_height_frame_convert — DEM `geodesy` op

• 資料種類:signal × signal → signal

• 呼叫: import fullseye as fs; fs.ledger.dem_height_frame_convert(height_m, geoid_height_m, frm='ellipsoidal', to='orthometric')(要直接呼叫實作,import demops; demops.dem_height_frame_convert(height_m, geoid_height_m, frm='ellipsoidal', to='orthometric');從台帳取用則 opsdem.get("dem_height_frame_convert"))

用法

轉換高度的基準：橢球高 h ↔ 正高 H(H = h - N)。signal。

同一個「高度」一詞下流通著 2 種不同的東西,而且兩者都不拋例外就登上地圖。本 op 強制寫明從哪一種轉換到哪一種 —— 雖然有預設值,但 frm 與 to 會留在程式碼中而不是 docstring 裡。

Args:

height_m: 高度 [m]。可為陣列。

geoid_height_m: 同一點的大地水準面高 N [m] —— 以 dem_geoid_height 查到的值。

純量或陣列皆可。

frm / to: "ellipsoidal"(橢球高)或 "orthometric"(正高)。

Returns:

(n,) 的轉換後高度 [m]。frm == to 時值不變。

可用閉合公式檢驗的部分:H = h - N 就是定義本身,所以真值由減法得出。往返(h→H→h)的誤差實測在 5e-13 m 以下(5,000 點,高度 −100~4000 m、N ±100 m)。並非逐位元一致 —— (h - N) + N 可能因捨入而變動 1 ulp, 實測中 5,000 點裡有 273 點如此。不以等號來寫「往返後恢復原值」,原因就在於此。

ValueError:未知的基準名、長度不一致、非有限。

詳細使用指南

• dem_terrain_analysis 族使用指南

參考(範例資料、文獻)

• 範例資料目錄(下載 URL / 授權) —— 2-D 用 skimage.data(BSD/公有領域)加合成圖,3-D 給出真實資料源(Stanford/PDS 等)的下載 URL。

• 運算子來歷與參考文獻 —— 該運算子族所依據的研究/方法出處。

• 演算法的正典(作者、年份)與用途見上面的族使用指南。

可執行的範例(實際呼叫該運算子並已驗證的樣例)

• poc_geodetic_benchmarks_real — py -3.11 examples/poc_geodetic_benchmarks_real.py

型別可銜接的下一個運算子(可接受 `signal` 作為輸入)

dem_height_frame_residual

同類別(`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_residual

*Provenance: demops.py — DEM 運算子登記表。本條目由 tools/opdocs.py md 自動產生(請勿手動編輯)。*

dem_height_frame_convert — DEM geodesy op

用法