dem_geodetic_slope — DEM `geodesy` op

• 数据种类:depth → image2d

• 调用: import fullseye as fs; fs.ledger.dem_geodetic_slope(dem, lat0_deg, d_lat_deg, d_lon_deg, method='horn', units='degrees')(要直接调用实现,import demops; demops.dem_geodetic_slope(dem, lat0_deg, d_lat_deg, d_lon_deg, method='horn', units='degrees');从台账取用则 opsdem.get("dem_geodetic_slope"))

用法

经纬度网格(等角度间隔)DEM 的坡度。格子单元尺寸随纬度变化。

公开的 DEM 大多像“1 秒网格(mesh)”那样按角度等间隔。把这种网格

直接以固定单元尺寸传给 dem_slope,就会忽略经度方向的实际距离短了 cos(纬度) 倍这一点。若两轴都用纬度方向的尺寸, 东西向坡度会偏小(北纬 60 度时减半);若两轴都用经度方向的尺寸, 则反过来南北向偏大。偏向哪一边取决于把什么当作常数,所以正确做法不是记住“偏大/偏小”,而是按纬度计算尺寸。

★ 最初这个 docstring 写的是“东西向会偏大”,但用测试测量后 变成了一半(偏小)。关于符号或方向的主张要先测量再写。

这里按纬度计算东西向的单元尺寸,把网格修正为等距离后再测量。

Args:

dem: (H, W)。第 0 行为北。

lat0_deg: 西北角的纬度 [度]。

d_lat_deg / d_lon_deg: 每个单元的纬度、经度步长 [度]。正值。

method / units: 与 dem_slope 相同。

详细使用指南

• 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

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

—

同类别(`geodesy`)

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

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

dem_geodetic_slope — DEM geodesy op

用法