fullseye

ransac_cylinder — 3D `robust_fit` op

データ種: points × normals → primitive
呼び出し: import fullseye as fs; fs.ledger.ransac_cylinder(points, normals, thresh, iters=800, seed=0) (実装を直接呼ぶなら import ransac_fit; ransac_fit.ransac_cylinder(points, normals, thresh, iters=800, seed=0)、台帳から引くなら ops3d.get("ransac_cylinder"))

使い方

外れ値に頑健な RANSAC 円筒適合(点法線が必要)。

円筒表面の法線は軸に直交するので、2 点の法線の外積で軸方向を推定 → 軸に直交な平面へ全点を投影 → その平面内で円をフィット → |投影距離 - r| < thresh の inlier を最大化。最終 inlier ではより頑健に軸を再推定(法線群の SVD の最小特異方向 = 軸)し、投影円を最小二乗リフィットする。法線が無ければ呼び出し側で estimate してから渡す。

Args: points: (N,3) 点群。 normals: (N,3) 各点の(単位)法線。 thresh: inlier とみなす |投影距離-r| のしきい値。 iters: RANSAC 反復数。 seed: 乱数シード(決定論)。

Returns: (params, inlier_mask, info)。params = {“axis”: (3,) 単位軸, “point”: (3,) 軸上の一点, “radius”: float}。info には inlier 数 n_inliers / 比 inlier_ratio / iters。

背景知識ガイド(この op の手前にある物理・規約)

blas_threads_and_memory — 行列分解が遅い理由の知識 — BLAS スレッド・キャッシュ・メモリ配置

参考(サンプルデータ・文献)

サンプルデータカタログ(DL URL / ライセンス) — 2-D は skimage.data(BSD/public)+ 合成、3-D は実データ源(Stanford/PDS 等)の DL URL。
演算子の来歴・参考文献 — この op 族の元になった研究/手法の出典。

実行できる例(この op を実際に呼ぶ検証済みサンプル)

ransac_prim — py -3.11 examples_3d/ransac_prim.py

型が繋がる次の op(`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 · fit_cone · fit_torus · fit_ellipsoid

Provenance: ransac_fit.py — 3D operator registry. この per-op ノートは tools/opdocs.py md が自動生成(手編集しない)。

This site is open source. Improve this page.