dem_geoid_height — DEM `geodesy` op

• 資料種類:depth → signal

• 呼叫: import fullseye as fs; fs.ledger.dem_geoid_height(geoid, lat_deg, lon_deg, lat0_deg, lon0_deg, d_lat_deg, d_lon_deg)(要直接呼叫實作,import demops; demops.dem_geoid_height(geoid, lat_deg, lon_deg, lat0_deg, lon0_deg, d_lat_deg, d_lon_deg);從台帳取用則 opsdem.get("dem_geoid_height"))

用法

依緯度、經度對大地水準面高網格做雙線性內插，回傳 N [m] (signal)。

GNSS 回傳的是橢球高 h,地圖與設計圖使用的是正高 H。兩者之差是大地水準面高 N(H = h - N),在日本附近大致為 +30~+40 m。本 op 是讀取已公開的大地水準面模型網格(如 GSIGEO / EGM 系列那樣的等間隔網格)的工具,本身不包含模型 —— 使用哪個版本由呼叫端負責,版本一變高程會變動幾 cm。

Args:

geoid: (H, W) 的大地水準面高 [m]。geoid[0, 0] 對應

(lat0_deg, lon0_deg),列每增加一列緯度增加 d_lat_deg,

行每增加一行經度增加 d_lon_deg(符號正負皆可)。

lat_deg / lon_deg: 要查詢的點的緯度、經度 [度]。可為陣列。

lat0_deg / lon0_deg: 網格原點 [度]。

d_lat_deg / d_lon_deg: 網格步長 [度]。拒絕 0。

Returns:

(n,) 的大地水準面高 [m]。

可用閉合公式檢驗的部分:對 N 為緯度、經度的 一次式的網格, 雙線性內插是精確的(機器精度)。有二次項時誤差以 |∂²N| * d² / 8 為上界。

ValueError:網格不是 2-D / 小於 2x2 / 非有限,步長為 0 或非有限, 緯度超出 ±90,以及網格外的點(不會默默以邊緣值填補海域或範圍外 —— 外推的大地水準面高不是測量值)。

詳細使用指南

• 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

• poc_geodetic_height_frames — py -3.11 examples/poc_geodetic_height_frames.py

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

dem_height_frame_convert · 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_height_frame_convert · dem_height_frame_residual

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

dem_geoid_height — DEM geodesy op

用法