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.
• blas_threads_and_memory — 행렬 분해가 느린 이유 — BLAS 스레드·캐시·메모리 배치
• 샘플 데이터 카탈로그(DL URL / 라이선스) —— 2-D 는 skimage.data(BSD/public)+ 합성, 3-D 는 실데이터 소스(Stanford/PDS 등)의 DL 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 연산자 레지스트리. 이 op 노트는 tools/opdocs.py md 가 자동 생성합니다(직접 편집하지 마세요).*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.