sdf_intersect — 3D sdf_csg op

データ種: sdf × sdfsdf

呼び出し: import fullseye as fs; fs.ledger.sdf_intersect(a, b) (実装を直接呼ぶなら import sdf_ops; sdf_ops.sdf_intersect(a, b)、台帳から引くなら ops3d.get("sdf_intersect"))

使い方

2 SDF の積集合 A∩B = 要素ごとの max(a, b)(両方の内側でのみ内側)。

ゼロ等値面は両境界の共通部分。外側で厳密 SDF、内側は保守的下界。

計算: `np.maximum(np.asarray(a, float64), np.asarray(b, float64))`。要素ごとの

max なので、ある点が A∩B の内側(負)になるのは `a < 0 かつ b < 0` のときだけ。

shape はブロードキャスト整合していればよく(不整合なら numpy の `ValueError`)、

それ以外の検査はしない。NaN は伝播、`±inf は厳密に伝播(max(a, -inf) = a`)。

入力: 同一グリッド上で評価した 2 つの SDF(`sphere_sdf / box_sdf` /

`esdf` の出力など、内側負・外側正)。距離の単位は入力と同じ。

返り値: ブロードキャスト後の shape の float64。`<= 0` が A∩B の占有。

注意:

外側の値は真の距離より小さく出うる: 交差の外側で、点が A の外かつ B の外の

とき `max(a, b)` は「遠い方の表面まで」の距離で、A∩B の表面はそれより遠い

ことがある(下界)。ゼロ等値面は厳密。

• 共通部分が無ければ全要素が正(占有 0)になり、エラーにはならない。

• 箱で球を切る・2 つの箱で角柱を作る、といった CSG の「切り出し」に使う。

A から B を抜くのは `sdf_subtract`。

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

• サンプルデータ カタログ(DL URL / ライセンス) — 2-D は skimage.data(BSD/public)+ 合成、3-D は実データ源(Stanford/PDS 等)の DL URL。

• 演算子の来歴・参考文献 — この op 族の元になった研究/手法の出典。

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

gear_metrologypy -3.11 examples_3d/gear_metrology.py

render_beautypy -3.11 examples_3d/render_beauty.py

sdf_csgpy -3.11 examples_3d/sdf_csg.py

型が繋がる次の op(sdf を入力に取れる)

sdf_to_occupancy · fuse_to_voxel · integrate · extract_surface_points · query_distance · sdf_union · sdf_subtract · sdf_smooth_union

同カテゴリ(sdf_csg)

grid_coords · sphere_sdf · box_sdf · plane_sdf · cylinder_sdf · torus_sdf · capsule_sdf · sdf_union


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

© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.