medial op• 資料種類:voxel → points
• 呼叫: import fullseye as fs; fs.ledger.medial_axis_points(vol, min_radius=0.0)(要直接呼叫實作,import medial; medial.medial_axis_points(vol, min_radius=0.0);從台帳取用則 ops3d.get("medial_axis_points"))
• 經台帳呼叫的回傳值: fullseye.ledger.medial_axis_points(...) 只回傳宣告的 out 型別 points 的值(底層函式還會回傳輔助資訊)。需要被捨棄的部分時,請用 fullseye.ledger.medial_axis_points.raw(...),或直接呼叫 medial.medial_axis_points。
• 底層回傳: (points, radii)
把 medial 體素的座標與局部半徑(= 該點的 EDT 值)變成點群。回傳 (points, radius)。
把粗的部分呈面狀、細的部分呈線狀分布的 medial 點,以 (M,3) 座標(z,y,x)及各自的局部半徑 (M,) 回傳。半徑最大的點指向形狀最「內側」= 中心。
Args:
vol: 二值 voxel(bool / 0-1 的 3D)。
min_radius: 排除半徑不超過此值的點(抑制雜訊)。
Returns:
points (float64, (M,3)): medial voxel 座標(z, y, x)。
radius (float64, (M,)): 各點的 EDT 值(= 局部內切半徑)。
步驟: 回傳由 distance_ridge(vol, min_radius) 得到的 ridge_mask 的 np.argwhere (陣列 index,列依 z-major 的字典序,具確定性)以及該位置的 edt 值。座標是 voxel index,不乘 spacing(需要物理座標時由呼叫端計算 points * (sz, sy, sx),半徑同樣乘以等向 spacing)。
引數與檢驗: vol 為 3-D(非零 = 前景)。不是 3-D、空、NaN/Inf 時為 ValueError,min_radius < 0 也為 ValueError。radius 只含 > min_radius 的點(邊界 voxel 的 EDT 不超過 1,因此用 min_radius=1 幾乎可以去掉外殼雜訊)。
邊界行為: 沒有前景,或全部點都不超過 min_radius 時,points 為 (0, 3)、radius 為 (0,) 的空陣列(不報錯)。
用途: medial_match 的半徑分布、通往點雲 op(smallest_sphere3 / fit_line3 等)的橋樑,用 np.argmax(radius) 取最大內切球的中心。
• 範例資料目錄(下載 URL / 授權) —— 2-D 用 skimage.data(BSD/公有領域)加合成圖,3-D 給出真實資料源(Stanford/PDS 等)的下載 URL。
• 運算子來歷與參考文獻 —— 該運算子族所依據的研究/方法出處。
• medial_topology — py -3.11 examples_3d/medial_topology.py
points 作為輸入)points_to_voxel · gaussians_to_voxel · estimate_point_normals · to_points · match_points_ncc · match_pca · moment_axes · icp_point2point_3d
medial)distance_ridge · skeletonize_vol · 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.