fullseye

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 ファミリガイド

参考(サンプルデータ・文献)

サンプルデータカタログ(DL URL / ライセンス) — 2-D は skimage.data(BSD/public)+ 合成、3-D は実データ源(Stanford/PDS 等)の DL URL。
演算子の来歴・参考文献 — この op 族の元になった研究/手法の出典。
アルゴリズムの正典(著者・年)と用途は上記ファミリ使い方ガイドに記載。

実行できる例(この op を実際に呼ぶ検証済みサンプル)

poc_geodetic_benchmarks_real — py -3.11 examples/poc_geodetic_benchmarks_real.py

型が繋がる次の op(`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 operator registry. この per-op ノートは tools/opdocs.py md が自動生成(手編集しない)。

This site is open source. Improve this page.