dem_height_frame_residual — DEM `geodesy` op

• 資料種類:signal × signal × signal → table

• 呼叫: import fullseye as fs; fs.ledger.dem_height_frame_residual(h_ellipsoidal_m, h_orthometric_m, geoid_height_m, tol_m=0.1)(要直接呼叫實作,import demops; demops.dem_height_frame_residual(h_ellipsoidal_m, h_orthometric_m, geoid_height_m, tol_m=0.1);從台帳取用則 opsdem.get("dem_height_frame_residual"))

用法

> 該運算子的說明尚無譯文,以下照原文給出。

`h - H - N の残差 —— 高さの取り違えを**静かにさせない**検出器(table`)。

楕円体高・標高・ジオイド高は 1 つの恒等式で結ばれている(`h = H + N`)。

3 つそろった点でその残差を測れば、**どれかが別の基準・別のモデル・別の版で

作られている**ことが数値で出る。測量成果どうしなら残差は cm 級に収まり、

楕円体高をそのまま標高として使っていれば残差は `-N`(日本付近で −30〜−40 m)

に張り付く。

Args:

h_ellipsoidal_m: 楕円体高 [m] (GNSS の返り)。配列可。

h_orthometric_m: 同じ点の標高 [m] (地図・水準測量)。

geoid_height_m: 同じ点のジオイド高 [m]。

tol_m: 「合っている」とみなす閾値 [m]。既定 0.10 m。

Returns:

``{"residual_m": (n,), "rms_m", "max_abs_m", "median_m", "n",

"n_over_tol", "tol_m", "fraction_over_tol"}``。

閉じた式で検査できること: 恒等式どおりに作った 3 つ組では残差は 0 になる

(実測 rms 1.2e-14 m。丸めの分だけ厳密な 0 ではない)。標高の列に楕円体高を

そのまま入れると残差は −N に張り付く —— 取り違えの量がそのまま出る。

ValueError: 長さの食い違い、非有限、負の *tol_m*。

詳細使用指南

• 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

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

—

同類別(`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 運算子登記表。本條目由 tools/opdocs.py md 自動產生(請勿手動編輯)。*

dem_height_frame_residual — DEM geodesy op

用法