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}。
以 scipy.optimize.least_squares 最小化子午面內點到母線的正交距離 a·sinα − ρ·cosα(a=軸向分量, ρ=半徑, α=半角)。初值把 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。
• blas_threads_and_memory — 矩陣分解為何慢 —— BLAS 執行緒、快取、記憶體配置
• 範例資料目錄(下載 URL / 授權) —— 2-D 用 skimage.data(BSD/公有領域)加合成圖,3-D 給出真實資料源(Stanford/PDS 等)的下載 URL。
• 運算子來歷與參考文獻 —— 該運算子族所依據的研究/方法出處。
• 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 運算子登記表。本條目由 tools/opdocs.py md 自動產生(請勿手動編輯)。*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.