min_enclosing_sphere — 3D bounds op

データ種: pointsprimitive

呼び出し: import fullseye as fs; fs.ledger.min_enclosing_sphere(points, refine_iters: 'int' = 1000) -> 'Dict[str, object]' (実装を直接呼ぶなら import hull3d; hull3d.min_enclosing_sphere(points, refine_iters: 'int' = 1000) -> 'Dict[str, object]'、台帳から引くなら ops3d.get("min_enclosing_sphere"))

使い方

点群 (N,3) → 全点を含む(近似)最小包含球 {center(3), radius}。

`fit_sphere_3d(球面フィット)や ransac_sphere / hough_sphere_3d`(球検出)とは

異なり、全点を内包する最小の球(minimum enclosing ball, MEB)を解く。2 段構成で

「全点内包」を厳守しつつ半径を詰める:

1. Ritter (1990) 初期化 — 最遠の点対を粗く取り初期球にし、各点を走査して球外の点が

あれば「その点と既存球の両方を含む」最小の球へ 1 回膨らませる(膨張式

`new_r=(r+d)/2 / 中心を点方向へ (d-r)/(2d)` 進める)。新球が旧球を完全に含むため、

1 パスで全点内包を保証する。

2. Bădoiu–Clarkson (2003) core-set 反復による精緻化 — 反復 `i で最遠点 q` へ

中心を `1/(i+2)` だけ寄せる。真の最小包含球へ単調収束する(半径過大な Ritter の

ドリフトを詰める)。最後に半径を「中心からの最大距離」で確定するので、精緻化後も

必ず全点を内包(近似ゆえ半径が過小になり点が漏れることはない、安全側)。

精緻化した中心が Ritter より外接半径を縮められたときのみ採用する(常に Ritter 以下)。

真の最小球(厳密解は Welzl の乱択線形時間法)ではなく高速な (1+ε) 近似。

Parameters

----------

points : array_like (N,3)

入力点群(>= 1 点)。

refine_iters : int

Bădoiu–Clarkson 精緻化の反復数(既定 1000)。0 で Ritter のみ。

Returns

-------

dict

• `center`: (3,) float64 — 球中心(世界座標)。

• `radius`: float — 半径(全点を内包)。

Raises

------

ValueError

形状不正・非有限・点数 0、または `refine_iters` が負のとき(fail-closed)。

背景知識ガイド(この op の手前にある物理・規約)

blender_interop — Blender との併用 — 形を作って fullseye で測る(軸・単位・正解データの罠)

参考(サンプルデータ・文献)

• サンプルデータ カタログ(DL URL / ライセンス) — 2-D は skimage.data(BSD/public)+ 合成、3-D は実データ源(Stanford/PDS 等)の DL URL。

• 演算子の来歴・参考文献 — この op 族の元になった研究/手法の出典。

実行できる例(この op を実際に呼ぶ検証済みサンプル)

hull_boundspy -3.11 examples_3d/hull_bounds.py

型が繋がる次の op(primitive を入力に取れる)

fuse_to_voxel · angle_between_lines · angle_between_planes · angle_line_plane · distance_point_plane · distance_point_line · distance_line_line · distance_segment_segment

同カテゴリ(bounds)

convex_hull · aabb · obb


*Provenance: hull3d.py — 3D operator registry. この per-op ノートは tools/opdocs.py md が自動生成(手編集しない)。*

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