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}。

用 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.