dem_slope — DEM `surface` op

• 数据种类:depth → image2d

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

用法

> 该算子的说明尚无译文,以下照原文给出。

傾斜角。`units は "degrees" / "radians" / "percent"`。

平面なら閉形式 `atan(|grad|)` と一致する(実測 1e-13 以下)。

手順: 3x3 近傍から `dz/dx(東向き)と dz/dy`(北向き = 行が減る向き)を

取り、`g = hypot(dz/dx, dz/dy)。"degrees" は atan(g)` を度で、

`"radians" はラジアン、"percent" は 100 g`(45 度 = 100 %)。

縁は端の値を複製して埋める(`pad mode="edge"`)ので、外周 1 セルは内側より

緩めに出る。

• `dem: (H, W) の標高 [m]、3x3 以上、実数、inf 不可。欠測は nan` で

渡す ―― `-9999 のような番兵値がそのまま入っている(<= -9000` があり

`nan が無い)と拒否する。nan` は 3x3 の範囲に伝播する。

• `cell_size`: セル辺長 [m]、正の有限値。bool / 文字列は拒否。緯度経度格子は

先に `dem_geodetic_slope` 側を使う。

• `method: "horn"(Horn 1981 の 3x3 重み付き差分、既定)/ "central"`

(Zevenbergen–Thorne の中央差分)。

• 返り値: `(H, W) float64。degrees は [0, 90)`。

• 失敗はすべて `ValueError(形・番兵値・cell_size`・選択肢)。

`dem_aspect と対で使う。dem_hillshade` はこの 2 つから陰影を作る。

详细使用指南

• 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

• poc_crop_phenotyping — py -3.11 examples/poc_crop_phenotyping.py

• poc_dem_terrain — py -3.11 examples/poc_dem_terrain.py

• poc_geodetic_height_frames — py -3.11 examples/poc_geodetic_height_frames.py

• poc_multibeam_bathymetry — py -3.11 examples/poc_multibeam_bathymetry.py

• poc_stockpile_volume — py -3.11 examples/poc_stockpile_volume.py

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

—

同类别(`surface`)

dem_aspect · dem_curvature · dem_roughness · dem_tpi

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

dem_slope — DEM surface op

用法