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

用法