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 軸 +
半径 ρ の軸成分 t に対する線形回帰(`ρ = m·t + b` の傾き m=tanα, 切片ゼロ点=頂点)。
軸の向きは「頂点から離れるほど ρ が増える(+方向に開く)」に正規化する。
Args:
points: (N,3) 点群(最低 6 点)。
Returns:
dict: ``{"apex": (3,), "axis": (3,) 単位軸(開く向き), "half_angle": float [rad],
"residual": float 点-面距離の RMS}``。
Raises:
ValueError: 形状不正/点数不足/半角 ~0(円柱へ縮退)など fail-closed。
• blas_threads_and_memory — 行列分解が遅い理由の知識 — BLAS スレッド・キャッシュ・メモリ配置
• サンプルデータ カタログ(DL URL / ライセンス) — 2-D は skimage.data(BSD/public)+ 合成、3-D は実データ源(Stanford/PDS 等)の DL URL。
• 演算子の来歴・参考文献 — この op 族の元になった研究/手法の出典。
• fit_primitives_ext — py -3.11 examples_3d/fit_primitives_ext.py
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 が自動生成(手編集しない)。*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.