dem_geocentric_grid — DEM `geodesy` op

• 数据种类:depth → pointmap

• 调用: import fullseye as fs; fs.ledger.dem_geocentric_grid(dem, lat0_deg, lon0_deg, cell_size, spherical=False)(要直接调用实现,import demops; demops.dem_geocentric_grid(dem, lat0_deg, lon0_deg, cell_size, spherical=False);从台账取用则 opsdem.get("dem_geocentric_grid"))

用法

将 DEM 的各格子单元转换为从地球中心看到的坐标。

把 (H, W) 的高程网格视为西北角位于 (lat0_deg, lon0_deg) 的局部平面, 将各单元映射到 ECEF(默认)或地心球坐标。

Args:

dem: (H, W) 的高程 [m]。

lat0_deg / lon0_deg: 网格西北角的经纬度 [度]。

cell_size: 单元尺寸 [m]。视纬度方向、经度方向相同

(Web 墨卡托的瓦片就是如此。

参见 dem_cell_size_webmercator)。

spherical: 为 True 时返回 (半径 r[m], 地心纬度[度], 经度[度])。

为 False(默认)时返回 ECEF (x, y, z) [m]。

Returns:

(H, W, 3)。

★ 地心纬度不是大地纬度。spherical=True 返回的纬度是从地球中心看的角度,与地图的纬度(椭球法线与赤道面所成的角) 最多相差 0.19 度。要回到地图时,请把 ECEF 一侧

传给 dem_ecef_to_geodetic。

详细使用指南

• dem_terrain_analysis 族使用指南

参考(示例数据、文献)

• 示例数据目录(下载 URL / 许可证) —— 2-D 用 skimage.data(BSD/公有领域)加合成图,3-D 给出真实数据源(Stanford/PDS 等)的下载 URL。

• 算子来历与参考文献 —— 该算子族所依据的研究/方法出处。

• 算法的正典(作者、年份)与用途见上面的族使用指南。

可运行的示例(实际调用该算子并已验证的样例)

• dem_geodesy_tour — py -3.11 examples/dem_geodesy_tour.py

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

—

同类别(`geodesy`)

dem_geodetic_to_ecef · dem_ecef_to_geodetic · dem_earth_curvature_drop · dem_cell_size_webmercator · dem_geodetic_slope · dem_geoid_height · dem_height_frame_convert · dem_height_frame_residual

*Provenance: demops.py — DEM 算子登记表。本条目由 tools/opdocs.py md 自动生成(请勿手工编辑)。*

dem_geocentric_grid — DEM geodesy op

用法