geodesic_mesh — 3D geodesic op

Data kinds: meshsignal

Call: import fullseye as fs; fs.ledger.geodesic_mesh(vertices: numpy.ndarray, faces: numpy.ndarray, source: int) -> numpy.ndarray (to call the implementation directly, import geodesic3d; geodesic3d.geodesic_mesh(vertices: numpy.ndarray, faces: numpy.ndarray, source: int) -> numpy.ndarray; from the registry, ops3d.get("geodesic_mesh"))

Usage

Geodesic distance from a source to each vertex by Dijkstra on the triangle mesh's edge graph. → (V,) float.

> The detailed description below is the original text — the summary and the headings are translated.

各三角形の 3 辺 (0,1),(1,2),(2,0) を無向辺として集め、(min,max) で一意化してから辺長 = 両端

頂点の Euclid 距離を重みにした疎グラフを作り、`dijkstra(directed=False)` で source からの

最短路長を返す。`d[source] = 0、source と辺で繋がらない頂点は inf`。単位は頂点座標の単位。

• `vertices: (V,3) 頂点座標(float64 に変換)。faces`: (M,3) 頂点添字(int に変換)。

• `source`: 始点頂点の添字。

• `faces が空のときは全頂点 inf`(source が範囲内ならそこだけ 0)を返し、例外は出さない。

実装上の要点: `csr_matrix` は同じ (i,j) を重複して渡すと重みを黙って加算するため、重複面・

非多様体・巻き順が不揃いなメッシュでも距離が膨らまないよう辺を一意化している。退化辺(i==j)は

捨てる。距離は辺に沿った折れ線長なので、粗いメッシュでは真の測地距離より長め(辺の走り方に

依存する異方性)になる。`voxel_to_meshconvex_hull` が返す (verts, faces) を

そのまま渡せる。

References (sample data, literature)

• Sample-data catalog (download URLs / licences) — 2-D uses skimage.data (BSD/public domain) plus synthetic images; 3-D lists download URLs for real data sources (Stanford, PDS, …).

• Operator provenance and references — the sources of the research/methods this op family came from.

Runnable examples (verified samples that actually call this op)

pcl_geodesicpy -3.11 examples_3d/pcl_geodesic.py

Ops the type connects to (they accept signal as input)

fuse_to_voxel

Same category (geodesic)

geodesic_distances · farthest_point_sampling · knn_graph


*Provenance: geodesic3d.py — 3D operator registry. This per-op note is generated by tools/opdocs.py md (do not hand-edit).*

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