skeleton_branches3d — 3D medial op

数据种类:voxelvoxel

调用: import fullseye as fs; fs.ledger.skeleton_branches3d(vol, min_length=0)(要直接调用实现,import medial; medial.skeleton_branches3d(vol, min_length=0);从台账取用则 ops3d.get("skeleton_branches3d"))

用法

把 3D 骨架在分叉点处切成枝(线段)。2D r2_split_skeleton_lines 的 3D 版。

> 以下的详细说明为原文 —— 摘要与标题已翻译。

分岐点 voxel を除いた残りが枝。min_length > 0 なら、26 連結成分の voxel 数が

それ未満の断片を除去する。

手順: `_ensure_skeleton で bool 化(interior voxel があれば skeletonize_vol`

で細線化)→ 26 近傍次数 `>= 3 の voxel(分岐点)を取り除く → min_length > 0`

なら `scipy.ndimage.label`(3x3x3 全 1 構造 = 26 連結)で成分に分け、voxel 数が

`min_length` 未満の成分を落とす。

返り値: 入力と同形の bool 配列(枝 voxel = True)。枝ごとのラベルは返さない

— 枝を個別に扱うには返り値を `vol_label(branches, 26)` に通す(成分数 = 枝数)。

分岐点 voxel そのものは結果に含まれないので、枝の両端は分岐点の 1 voxel 手前で

終わる。

引数と検証: `min_length は voxel 数(int()` で切り捨て。0 なら除去しない)。

入力が 3-D でない・空・NaN/Inf は `ValueError`。細線化が必要で

scikit-image が無ければ `ImportError`。前景が無ければ全 False。

注意: 26 近傍次数は分岐の対角隣接で 3 以上になりやすく、分岐点が数 voxel の

塊として除かれるため、枝が実際より短く出ることがある。閉ループだけの骨格は

分岐点が無く、全体が 1 本の枝として残る。

参考(示例数据・文献)

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

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

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

medial_topologypy -3.11 examples_3d/medial_topology.py

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

voxel_to_mips · voxel_to_mesh · signed_distance_field · to_points · sobel3d · hessian3d · curvature_maps · edt_jfa

同类别(medial)

distance_ridge · skeletonize_vol · medial_axis_points · topology_signature · medial_match · skeleton_junctions3d · skeleton_endpoints3d · skeleton_prune3d


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

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