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] —— :func: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

用法