medial op• 資料種類:voxel → voxel
• 呼叫: import fullseye as fs; fs.ledger.distance_ridge(vol, min_radius=0.0)(要直接呼叫實作,import medial; medial.distance_ridge(vol, min_radius=0.0);從台帳取用則 ops3d.get("distance_ridge"))
• 經台帳呼叫的回傳值: fullseye.ledger.distance_ridge(...) 只回傳宣告的 out 型別 voxel 的值(底層函式還會回傳輔助資訊)。需要被捨棄的部分時,請用 fullseye.ledger.distance_ridge.raw(...),或直接呼叫 medial.distance_ridge。
• 底層回傳: (ridge, dist)
把 EDT 的脊(距離場的局部極大)提取為中軸。回傳 (ridge_mask, edt)。
計算各前景 voxel 的 EDT,把 26 鄰域的局部極大(自身 EDT 不小於周圍 26 個 voxel 的最大值) 視為 medial。此判據會依物體的局部維數自然地給出相應維數:
塊(球) -> EDT 單峰 -> 點狀 medial(中心 1 點)。
管(圓柱)-> 沿軸向平坦、沿徑向單峰 -> 線狀 medial(軸線)。
板(平板)-> 面內平坦、沿厚度方向單峰 -> 面狀 medial(中心面)。
邊界 voxel 的內側鄰居必然更大,因此不會成為極大,外殼自然被去除。平坦的脊(軸/面)因允許等值鄰接(>=)而作為連續的線/面保留下來。
Args:
vol: 二值 voxel(bool / 0-1 的 3D)。
min_radius: 捨棄 EDT 值不超過此值的弱脊的閾值(預設 0 = 不捨棄)。用於抑制雜訊。
Returns:
ridge_mask (bool 3D): medial voxel。
edt (float64 3D): 各 voxel 到背景的歐氏距離(= 局部半徑)。
• 範例資料目錄(下載 URL / 授權) —— 2-D 用 skimage.data(BSD/公有領域)加合成圖,3-D 給出真實資料源(Stanford/PDS 等)的下載 URL。
• 運算子來歷與參考文獻 —— 該運算子族所依據的研究/方法出處。
• pcl_geodesic — py -3.11 examples_3d/pcl_geodesic.py
voxel 作為輸入)voxel_to_mips · voxel_to_mesh · signed_distance_field · to_points · sobel3d · hessian3d · curvature_maps · edt_jfa
medial)skeletonize_vol · medial_axis_points · topology_signature · medial_match · skeleton_junctions3d · skeleton_endpoints3d · skeleton_prune3d · skeleton_branches3d
*Provenance: medial.py — 3D 運算子登記表。本條目由 tools/opdocs.py md 自動產生(請勿手動編輯)。*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.