dem_datum_shift_3param — DEM `geodesy` op

• 数据种类:无 → points(仅由参数决定的算子 —— 不接受图像或数据输入)

• 调用: import fullseye as fs; fs.ledger.dem_datum_shift_3param(lat_deg, lon_deg, h_m, dx_m, dy_m, dz_m, a_from_m, f_from, a_to_m, f_to)(要直接调用实现,import demops; demops.dem_datum_shift_3param(lat_deg, lon_deg, h_m, dx_m, dy_m, dz_m, a_from_m, f_from, a_to_m, f_to);从台账取用则 opsdem.get("dem_datum_shift_3param"))

用法

大地基准(datum)的转换 —— 地心 3 参数平移(points)。

即使同为“北纬 35 度 41 分”，旧日本测地系与世界测地系所指的点在地面上相距数百米。本 op 将这一转换作为地心直角坐标中的平移(Molodensky-Badekas 中无旋转、无缩放的情形)来执行::

ECEF(a_from, f_from) + (dx, dy, dz) -> 大地坐标(a_to, f_to)

★参数不设默认值。从哪个成果到哪个成果、用哪个平移量转换，因国家、地区、版本而异，若默默假定，就会得到“不报异常但偏离数百米的坐标” —— 这正是本 op 存在的理由，所以必须让调用者写明。

Args:

lat_deg / lon_deg: 原成果中的纬度、经度 [度]。可为数组。

h_m: 原成果中的椭球高 [m]。

dx_m / dy_m / dz_m: 地心平移量 [m] (原 → 目标)。

a_from_m / f_from: 原椭球的长半径 [m] 与扁率。

a_to_m / f_to: 目标椭球的长半径 [m] 与扁率。

Returns:

(n, 3) 的 (纬度[度], 经度[度], 椭球高[m])。

可用闭式检验的内容：平移为 0 且为同一椭球时为恒等变换 (往返的下限与 dem_ecef_to_geodetic 相同)。只给出平移时的地面位移量，与该点 ECEF 基底下分解的分量一致。

ValueError：非有限值、纬度超出 ±90、长半径非正、扁率超出 [0, 1)。

详细使用指南

• 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

类型可衔接的下一个算子(可接受 `points` 作为输入)

dem_ecef_to_geodetic · dem_geodetic_from_enu

同类别(`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_datum_shift_3param — DEM geodesy op

用法