fullseye

fit_cone — 3D `robust_fit` op

データ種: points → primitive
呼び出し: import fullseye as fs; fs.ledger.fit_cone(points) -> 'dict' (実装を直接呼ぶなら import fit_primitives_ext; fit_primitives_ext.fit_cone(points) -> 'dict'、台帳から引くなら ops3d.get("fit_cone"))

使い方

点群に無限円錐を当てはめ {apex, axis, half_angle, residual} を返す。

子午面での点-母線直交距離 a·sinα − ρ·cosα(a=軸成分, ρ=半径, α=半角)を scipy.optimize.least_squares で最小化する。初期値は PCA の 3 主軸をそれぞれ軸候補にし、半径 ρ の軸成分 t に対する線形回帰 ρ = m·t + b から取る (半角 = atan(m)、頂点 = ρ が 0 になる軸上の点 t = -b/m)。残差最小の候補を採る。軸の向きは「頂点から離れるほど ρ が増える(+方向に開く)」に正規化する。

Args: points: (N,3) 点群(最低 6 点)。

Returns: dict: {"apex": (3,), "axis": (3,), "half_angle": float, "residual": float}。 axis は単位軸(開く向き)、half_angle は半角 [rad]、residual は点-面距離の RMS。

Raises: ValueError: 形状不正/点数不足/半角 ~0(円柱へ縮退)など fail-closed。

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

blas_threads_and_memory — 行列分解が遅い理由の知識 — BLAS スレッド・キャッシュ・メモリ配置

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

サンプルデータカタログ(DL URL / ライセンス) — 2-D は skimage.data(BSD/public)+ 合成、3-D は実データ源(Stanford/PDS 等)の DL URL。
演算子の来歴・参考文献 — この op 族の元になった研究/手法の出典。

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

fit_primitives_ext — py -3.11 examples_3d/fit_primitives_ext.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

同カテゴリ(`robust_fit`)

ransac_plane · ransac_sphere · ransac_line · ransac_cylinder · fit_torus · fit_ellipsoid

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

This site is open source. Improve this page.