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.