dem_tpi — DEM surface op

• 数据种类:depth → image2d

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

用法

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

地形位置指数 TPI —— 自セルと 8 近傍平均の差 [m]。正が尾根、負が谷。

式: `tpi = z - mean(8 近傍の z)`(Weiss の TPI を 3x3 の最小近傍で取ったもの)。

縁は端の値を複製して埋める。近傍が 1 セル固定なので、**スケールは

`cell_size` 1 つぶんに固定**で、より広い尾根/谷を見たいときは先に

`dem` を粗くする(この op に半径の引数は無い)。

• `dem: (H, W) [m]、3x3 以上、実数、inf 不可。欠測は nan`(番兵値

`<= -9000 が nan 無しで入っていれば拒否)。nan` は 3x3 に伝播する。

• `cell_size`: 結果には影響しない(差分に距離は入らない)が、単位が [m] である

ことを呼び出し側に意識させるため受け取り、正の有限値かを検査する。

• 返り値: `(H, W)` float64 [m]。0 付近が斜面の途中、正が凸(尾根・頂)、

負が凹(谷・窪地)。

• 失敗: `ValueError(形、番兵値、cell_size`)。

`dem_roughness`(同じ近傍の RMS)と組み合わせると地形分類の特徴になる。

详细使用指南

• dem_terrain_analysis 族使用指南

参考(示例数据・文献)

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

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

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

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

• dem_terrain_analysis_tour — py -3.11 examples/dem_terrain_analysis_tour.py

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

—

同类别(surface)

dem_slope · dem_aspect · dem_curvature · dem_roughness

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

© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.