medial_axis_points — 3D 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)` で最大内接球の中心を取る。

참고(샘플 데이터·문헌)

• 샘플 데이터 카탈로그(DL URL / 라이선스) —— 2-D 는 skimage.data(BSD/public)+ 합성, 3-D 는 실데이터 소스(Stanford/PDS 등)의 DL 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 연산자 레지스트리. 이 op 노트는 tools/opdocs.py md 가 자동 생성합니다(직접 편집하지 마세요).*

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