fit_ellipsoid — 3D `robust_fit` op

• 데이터 종류: points → primitive

• 호출: import fullseye as fs; fs.ledger.fit_ellipsoid(points) -> 'dict'(구현을 직접 호출하려면 import fit_primitives_ext; fit_primitives_ext.fit_ellipsoid(points) -> 'dict', 원장에서 가져오려면 ops3d.get("fit_ellipsoid"))

사용법

점군에 임의 자세의 삼축 타원체를 대수 적합해 `{center, axes, radii, residual}` 을 돌려줍니다.

> 아래 상세 설명은 원문입니다 —— 요약과 제목은 번역되어 있습니다.

一般二次曲面 `xᵀA x + b·x + c = 0 を、楕円体を保証する拘束 4J − I² = 1` の下で

一般化固有問題 `Sr v1 = λ C v1(scipy.linalg.eig`)として解く(Li & Griffiths 2004)。

各実固有ベクトルを楕円体へ復元し、正定値(= 実在する楕円体)へ復元でき残差 RMS が最小

のものを採用する(復元時の正定値検査そのものが厳密な楕円体判定)。初期値不要・決定論・

大域解。数値安定化のため点群を

重心と RMS 半径で無次元化してから解き、パラメータを world 座標へ戻す。

Args:

points: (N,3) 点群(最低 10 点)。外れ値には無防備(必要なら事前に inlier 選別)。

Returns:

dict: ``{"center": (3,), "axes": (3,3) 列=主軸(半径降順), "radii": (3,) 半径(降順),

"residual": float Taubin 近似の点-面距離 RMS}``。

Raises:

ValueError: 形状不正/点数不足/正定値な楕円体解が得られない(平面状の退化・

非楕円面・被覆不足)など 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_cone · fit_torus

*Provenance: fit_primitives_ext.py — 3D 연산자 레지스트리. 이 op 노트는 tools/opdocs.py md 가 자동 생성합니다(직접 편집하지 마세요).*

fit_ellipsoid — 3D robust_fit op

사용법