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.

• 연산자의 내력·참고문헌 —— 이 연산자 족의 바탕이 된 연구/기법의 출처.

• 알고리즘의 정전(저자·연도)과 용도는 위의 패밀리 사용 가이드에 적혀 있습니다.

실행 가능한 예제(이 연산자를 실제로 호출하는 검증된 샘플)

• 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 연산자 레지스트리. 이 op 노트는 tools/opdocs.py md 가 자동 생성합니다(직접 편집하지 마세요).*

© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.