Sub-Riemannian Snakes on the Projective Line Bundle with Applications to Segmentation of SEM Images
Pith reviewed 2026-05-10 15:16 UTC · model grok-4.3
The pith
A new snake model on the projective line bundle switches between fast spatial tracking and full geodesics to segment overlapping structures in SEM images without cusps or asymmetry.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We provide a practical solution to global tracking on the projective line bundle by introducing a snake model where we only compute the distance map where needed. Our method introduces a geometric criterion for switching between fast spatial snakes and computing minimizing geodesics of a new projective line bundle model. The new pseudo-distance underlying our geometric model is both symmetric and cusp-free, in contrast to previous geodesic sub-Riemannian models on R² × P¹. Our pseudo-distance satisfies the triangle inequality on a large set that we characterize, and includes a connected-component-informed cost function, which is highly advantageous in applications. Experiments on Scanning E
What carries the argument
The sub-Riemannian snake model on R² × P¹ that switches via a geometric criterion between fast spatial snakes and full minimizing geodesics, driven by a new symmetric cusp-free pseudo-distance that incorporates connected-component costs.
If this is right
- The pseudo-distance being symmetric and cusp-free removes two common failure modes of earlier sub-Riemannian models on the same bundle.
- The connected-component-informed cost term improves handling of separate but nearby objects in segmentation tasks.
- The characterized set on which the triangle inequality holds provides a concrete domain where the pseudo-distance behaves metrically.
- Automatic segmentation of overlapping electronic structures becomes feasible without manual intervention or prohibitive run times.
Where Pith is reading between the lines
- The switching mechanism could be adapted to other oriented tracking problems such as vessel or fiber segmentation in medical imaging.
- The explicit characterization of the triangle-inequality set opens the possibility of proving stronger metric properties or constructing true distances by restriction.
- Replacing the current cost function with data-driven terms learned from annotated SEM images might further reduce the need for the geometric switch.
- The same snake construction could be tested on non-SEM modalities that exhibit similar overlapping oriented structures to check generality beyond the reported experiments.
Load-bearing premise
The geometric criterion for switching between fast spatial snakes and full minimizing geodesics correctly identifies the regions where expensive full computation is required without missing critical features or causing segmentation errors on real SEM images.
What would settle it
A collection of SEM images containing overlapping structures in which the switching criterion selects fast snakes yet the resulting contours contain breaks, merges, or missed overlaps that disappear when full geodesic computation is forced everywhere.
read the original abstract
Geodesic tracking on the projective line bundle $\R^2 \times P^1 $ has many uses, including the segmentation of objects in images. However, global tracking requires expensive distance map computations. We provide a practical solution to this problem by introducing a snake model on $\R^2 \times P^1$, where we only compute the distance map where needed. Our method introduces a geometric criterion for switching between fast spatial snakes and computing minimizing geodesics of a new projective line bundle model. The new pseudo-distance underlying our geometric model is both symmetric and cusp-free, in contrast to previous geodesic sub-Riemannian models on $\R^2 \times P^1$. Our pseudo-distance satisfies the triangle inequality on a large set that we characterize, and includes a connected-component-informed cost function, which is highly advantageous in applications. Experiments on Scanning Electron Microscopy (SEM) images demonstrate our method's robust, automatic segmentation of overlapping electronic structures.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a practical snake model on the projective line bundle R² × P¹ for geodesic tracking in images. It introduces a geometric switching criterion that alternates between fast spatial snakes and full minimizing geodesics of a new pseudo-distance; this pseudo-distance is asserted to be symmetric and cusp-free (unlike prior sub-Riemannian models), to obey the triangle inequality on a characterized large set, and to incorporate a connected-component-informed cost. The approach is applied to automatic segmentation of overlapping electronic structures in SEM images, avoiding the expense of global distance-map computation.
Significance. If the switching criterion is reliable and the pseudo-distance properties hold, the work would deliver a computationally tractable sub-Riemannian snake method that retains the advantages of orientation-aware tracking while mitigating cusps and symmetry issues. The connected-component cost and the explicit characterization of the triangle-inequality domain are potentially useful extensions for segmentation tasks involving overlapping objects.
major comments (3)
- [§4.2] §4.2 (switching criterion): The geometric criterion for deciding when to invoke the full pseudo-distance geodesics rather than spatial snakes is presented only descriptively; no proof or exhaustive case analysis is given that it detects precisely the loci where the cheaper approximation deviates from the true minimizer. This is load-bearing for the central efficiency claim, as failure on global topology of overlapping structures could produce under-segmentation or spurious boundaries despite the claimed cusp-freeness.
- [§3.1] §3.1, Definition of the pseudo-distance: The symmetry and cusp-freeness are asserted by construction, yet the explicit formula and the argument that it eliminates the cusps of earlier R² × P¹ models are not supplied; without them the contrast with previous work cannot be verified and the triangle-inequality characterization on the “large set” remains uncheckable.
- [§5] §5 (experiments): The SEM segmentation results are described as “robust and automatic,” but no quantitative metrics (Dice, Hausdorff, or boundary error) or ablation on the switching threshold appear; the single qualitative figure does not substantiate that the criterion avoids missing critical features across varied overlapping configurations.
minor comments (2)
- [§2] Notation for the projective line bundle is introduced without a clear statement of the identification P¹ ≅ S¹/∼; a short paragraph clarifying the quotient would aid readers unfamiliar with sub-Riemannian geometry on bundles.
- [§3.3] The connected-component cost function is mentioned as “highly advantageous” but its precise additive term and weighting parameter are not written explicitly; adding the formula would make the model reproducible.
Simulated Author's Rebuttal
We thank the referee for the careful and constructive review of our manuscript. We address each major comment point by point below, indicating the revisions we will make to improve clarity and substantiation.
read point-by-point responses
-
Referee: [§4.2] §4.2 (switching criterion): The geometric criterion for deciding when to invoke the full pseudo-distance geodesics rather than spatial snakes is presented only descriptively; no proof or exhaustive case analysis is given that it detects precisely the loci where the cheaper approximation deviates from the true minimizer. This is load-bearing for the central efficiency claim, as failure on global topology of overlapping structures could produce under-segmentation or spurious boundaries despite the claimed cusp-freeness.
Authors: The switching criterion is derived directly from the pseudo-distance properties and the geometry of R² × P¹: it triggers full geodesic computation precisely when the spatial snake path would incur an orientation cost exceeding the connected-component term, which by construction occurs at loci where the approximation deviates due to potential cusps or topology changes in overlaps. While a complete exhaustive case analysis over all continuous configurations is intractable, we will expand §4.2 with a formal justification based on the triangle-inequality domain and include illustrative examples of overlapping structures to demonstrate correct detection. revision: partial
-
Referee: [§3.1] §3.1, Definition of the pseudo-distance: The symmetry and cusp-freeness are asserted by construction, yet the explicit formula and the argument that it eliminates the cusps of earlier R² × P¹ models are not supplied; without them the contrast with previous work cannot be verified and the triangle-inequality characterization on the “large set” remains uncheckable.
Authors: We will supply the explicit formula for the pseudo-distance in the revised §3.1: it is the infimum of integrated costs over admissible curves in R² × P¹ using a symmetric sub-Riemannian metric that treats θ and θ+π equivalently, augmented by the connected-component cost. Cusp-freeness follows because the projective identification removes the need for orientation reversal penalties present in prior asymmetric models; symmetry is immediate from the metric definition. The triangle-inequality domain will be characterized explicitly as the set where the connected-component term dominates local orientation costs, with a proof sketch added. revision: yes
-
Referee: [§5] §5 (experiments): The SEM segmentation results are described as “robust and automatic,” but no quantitative metrics (Dice, Hausdorff, or boundary error) or ablation on the switching threshold appear; the single qualitative figure does not substantiate that the criterion avoids missing critical features across varied overlapping configurations.
Authors: We agree that quantitative metrics and ablation would strengthen the experimental validation. In the revised manuscript we will add Dice coefficients, Hausdorff distances, and boundary error statistics computed against manual annotations on a collection of SEM images exhibiting varied overlaps. We will also include an ablation study varying the switching threshold to quantify its effect on avoiding under-segmentation. revision: yes
Circularity Check
No significant circularity; new pseudo-distance and switching criterion introduced independently
full rationale
The paper defines a new pseudo-distance on R² × P¹ with asserted properties (symmetry, cusp-freeness, triangle inequality on a characterized set, plus connected-component cost) and a geometric switching criterion between spatial snakes and full geodesics. These are presented as novel constructions whose properties are stated and verified separately from inputs. No equations or definitions reduce by construction to fitted parameters, self-citations, or prior ansatzes from the same authors. The derivation chain relies on explicit new model elements and experimental validation on SEM images rather than tautological renaming or load-bearing self-reference. This is the expected non-finding for a paper introducing an applied geometric model.
Axiom & Free-Parameter Ledger
axioms (1)
- standard math Sub-Riemannian geometry on the projective line bundle R² × P¹ provides a suitable structure for direction-aware geodesic tracking in images
invented entities (1)
-
Symmetric cusp-free pseudo-distance on R² × P¹
no independent evidence
Reference graph
Works this paper leans on
-
[2]
Journal of Mathematical Imaging and Vision24, 307–326 (2006) https://doi.org/10
Citti, G., Sarti, A.: A Cortical Based Model of Perceptual Completion in the Roto-Translation Space. Journal of Mathematical Imaging and Vision24, 307–326 (2006) https://doi.org/10. 1007/s10851-005-3630-2
work page 2006
-
[3]
Journal of Mathematical Imaging and Vision14, 225–236 (2001) https://doi.org/10.1023/A:1011281928379
Cohen, L.D.: Multiple Contour Finding and Per- ceptual Grouping Using Minimal Paths. Journal of Mathematical Imaging and Vision14, 225–236 (2001) https://doi.org/10.1023/A:1011281928379
-
[4]
Benmansour, F., Cohen, L.D.: Tubular Structure Segmentation Based on Minimal Path Method and Anisotropic Enhancement. International Journal of Computer Vision92, 192–210 (2011) https: //doi.org/10.1007/s11263-010-0331-0
-
[5]
International Journal of Computer Vision122, 458–483 (2017) https://doi.org/10
Chen, D., Mirebeau, J.-M., Cohen, L.D.: Global Minimum for a Finsler Elastica Minimal Path Approach. International Journal of Computer Vision122, 458–483 (2017) https://doi.org/10. 1007/s11263-016-0975-5
work page 2017
-
[6]
IJCV132, 349–391 (2024) https://doi.org/10.1137/15M1018460
Chen, D., Mirebeau, J.M., Shu, H., Cohen, L.D.: A Region-Based Randers Geodesic Approach for Image Segmentation. IJCV132, 349–391 (2024) https://doi.org/10.1137/15M1018460
-
[7]
Journal of Mathematical Imaging and Vision60, 816–848 (2018) https://doi.org/10
Duits, R., Meesters, S.P.L., Mirebeau, J.-M., Portegies, J.M.: Optimal Paths for Variants of the 2D and 3D Reeds–Shepp Car with Applications in Image Analysis. Journal of Mathematical Imaging and Vision60, 816–848 (2018) https://doi.org/10. 1007/s10851-018-0795-z
work page 2018
-
[8]
Differential Geometry and its Applications 65, 55–77 (2019) https://doi.org/10.1016/j.difgeo
Franceschiello, B., Mashtakov, A., Citti, G., Sarti, A.: Geometrical optical illusion via sub- Riemannian geodesics in the roto-translation group. Differential Geometry and its Applications 65, 55–77 (2019) https://doi.org/10.1016/j.difgeo. 2019.03.007
-
[9]
In: Nielsen, F., Bar- baresco, F
Bekkers, E.J., Duits, R., Mashtakov, A., Sachkov, Y.: Vessel Tracking via Sub-Riemannian Geodesics on the Projective Line Bundle. In: Nielsen, F., Bar- baresco, F. (eds.) Geometric Science of Informa- tion, pp. 773–781. Springer, Cham (2017). https: //doi.org/10.1007/978-3-319-68445-1 89
-
[10]
Boscain, U., Charlot, G., Rossi, F.: Existence of Planar Curves Minimizing Length and Curvature. Proceedings of the Steklov Institute of Mathemat- ics270, 43–56 (2010) https://doi.org/10.1134/ S0081543810030041
work page 2010
-
[11]
https://arxiv.org/abs/2409.18002
van den Berg, N.J., Mula, O., Vis, L., Duits, R.: Connected Components on Lie Groups and Applications to Multi-Orientation Image Analysis (2024). https://arxiv.org/abs/2409.18002
-
[12]
International Journal of Computer Vision1(4), 321–331 (1988) https:// doi.org/10.1007/BF00133570
Kass, M., Witkin, A., Terzopoulos, D.: Snakes: Active contour models. International Journal of Computer Vision1(4), 321–331 (1988) https:// doi.org/10.1007/BF00133570
-
[13]
In: Nielsen, M., Johansen, P., Olsen, O.F., Weickert, J
Chan, T., Vese, L.: An Active Contour Model without Edges. In: Nielsen, M., Johansen, P., Olsen, O.F., Weickert, J. (eds.) Scale-Space The- ories in Computer Vision, pp. 141–151. Springer, Berlin, Heidelberg (1999). https://doi.org/10. 1007/3-540-48236-9 13
work page 1999
-
[14]
Numerische Mathematik66, 1–31 (1993) https://doi.org/10.1007/BF01385685
Caselles, V., Catt´ e, F., Coll, T., Dibos, F.: A Geo- metric Model for Active Contours in Image Pro- cessing. Numerische Mathematik66, 1–31 (1993) https://doi.org/10.1007/BF01385685
-
[15]
Malladi, R., Sethian, J.A., Vemuri, B.C.: Shape modeling with front propagation: a level set approach. IEEE Transactions on Pattern Analy- sis and Machine Intelligence17(2), 158–175 (1995) https://doi.org/10.1109/34.368173
-
[16]
Cambridge University Press, Cambridge (2001)
Sapiro, G.: Geometric Partial Differential Equations and Image Analysis. Cambridge University Press, Cambridge (2001)
work page 2001
-
[17]
International Journal of Com- puter Vision22, 61–79 (1997) https://doi.org/10
Caselles, V., Kimmel, R., Sapiro, G.: Geodesic Active Contours. International Journal of Com- puter Vision22, 61–79 (1997) https://doi.org/10. 27 8 12 16 width in pixels 0.50 0.75 1.00 1.25 1.50 1.75 2.00dMASD in pixels New geodesic tracking model dc with new located at bottom top 8 12 16 width in pixels 0.50 0.75 1.00 1.25 1.50 1.75 2.00dMASD in pixels P...
work page 1997
-
[18]
Jost, J.: Riemannian Geometry and Geomet- ric Analysis, 7th edn. Universitext, p. 697. Springer, Cham (2017). https://doi.org/10.1007/ 978-3-319-61860-9
work page 2017
-
[19]
Grad- uate studies in mathematics, vol
Evans, L.C.: Partial Differential Equations. Grad- uate studies in mathematics, vol. 19. American Mathematical Society, Providence, R.I. (2010)
work page 2010
-
[20]
Journal of Computational Physics169(2), 463–502 (2001) https://doi.org/10.1006/jcph.2000.6636
Osher, S., Fedkiw, R.P.: Level Set Methods: An Overview and Some Recent Results. Journal of Computational Physics169(2), 463–502 (2001) https://doi.org/10.1006/jcph.2000.6636
-
[21]
Overcoming catastrophic forgetting in neural networks
Sethian, J.A., Vladimirsky, A.: Fast methods for the Eikonal and related Hamilton–Jacobi equations on unstructured meshes. In: Proceedings of the National Academy of Sciences, vol. 97, pp. 5699–5703 (2000). https://doi.org/10.1073/pnas. 090060097
-
[22]
Osher, S., Sethian, J.A.: Fronts propagating with curvature-dependent speed: Algorithms based on Hamilton-Jacobi formulations. Journal of Compu- tational Physics79(1), 12–49 (1988) https://doi. org/10.1016/0021-9991(88)90002-2
-
[23]
Vladimir Vovk, Alex Gammerman, and Glenn Shafer.Algorithmic learning in a random world
Tsitsiklis, J.N.: Efficient algorithms for globally optimal trajectories. IEEE Transactions on Auto- matic Control40(9), 1528–1538 (1995) https:// doi.org/10.1109/9.412624
-
[24]
In: Proceedings of the National Academy of Sciences, vol
Kimmel, R., Sethian, J.A.: Computing geodesic paths on manifolds. In: Proceedings of the National Academy of Sciences, vol. 95, pp. 8431– 8435 (1998). https://doi.org/10.1073/pnas.95.15. 8431
-
[25]
PhD thesis, University of Bath (2017) 28
Chen, D.: The visual analysis of complex natural phenomena. PhD thesis, University of Bath (2017) 28
work page 2017
-
[26]
SIAM Journal on Numerical Analysis52(4), 1573– 1599 (2014) https://doi.org/10.1137/120861667
Mirebeau, J.-M.: Anisotropic Fast-Marching on Cartesian Grids Using Lattice Basis Reduction. SIAM Journal on Numerical Analysis52(4), 1573– 1599 (2014) https://doi.org/10.1137/120861667
-
[27]
Image Processing On Line9, 47–93 (2019) https://doi.org/10.5201/ ipol.2019.227
Mirebeau, J.-M., Portegies, J.: Hamiltonian Fast Marching: A Numerical Solver for Anisotropic and Non-Holonomic Eikonal PDEs. Image Processing On Line9, 47–93 (2019) https://doi.org/10.5201/ ipol.2019.227
work page 2019
-
[28]
In: Bubba, T.A., Papafitsoros, K., Sch¨ onlieb, C.-B
van den Berg, N.J., Sherry, F.M., Berendschot, T.T.J.M., Duits, R.: Crossing-Preserving Geodesic Tracking on Spherical Images. In: Bubba, T.A., Papafitsoros, K., Sch¨ onlieb, C.-B. (eds.) Scale Space and Variational Methods in Computer Vision, pp. 192–204. Springer, Cham (2025). https: //doi.org/10.1007/978-3-031-92369-2 15
-
[29]
Sanguinetti, G., Bekkers, E., Duits, R., Janssen, M., Mashtakov, A., Mirebeau, J.-M.: Sub- Riemannian fast marching inSE(2). In: Pardo, A., Kittler, J. (eds.) Progress in Pattern Recog- nition, Image Analysis, Computer Vision, and Applications : CIARP 2015: Proceedings. LNCS, pp. 366–374. Springer, Cham (2015). https://doi. org/10.1007/978-3-319-25751-8 44
-
[30]
van den Berg, N.J., Smets, B.M.N., Pai, G., Mire- beau, J.-M., Duits, R.: Geodesic Tracking via New Data-Driven Connections of Cartan Type for Vascular Tree Tracking. Journal of Mathemati- cal Imaging and Vision66, 198–230 (2024) https: //doi.org/10.1007/s10851-023-01170-x
-
[31]
Houben, T., Huisman, T., Pisarenco, M., van der Sommen, F., de With, P.: Training proce- dure for scanning electron microscope 3D sur- face reconstruction using unsupervised domain adaptation with simulated data. Journal of Micro/Nanopatterning, Materials, and Metrology 22(3) (2023) https://doi.org/10.1117/1.JMM.22. 3.031208
-
[32]
Lindeberg, T.: Edge Detection and Ridge Detec- tion with Automatic Scale Selection. International Journal of Computer Vision30(2), 117–156 (1998) https://doi.org/10.1023/A:1008097225773
-
[33]
Pattern Recognition Let- ters7(3), 181–190 (1988) https://doi.org/10.1016/ 0167-8655(88)90063-3
Sj¨ oberg, F., Bergholm, F.: Extraction of diffuse edges by edge focusing. Pattern Recognition Let- ters7(3), 181–190 (1988) https://doi.org/10.1016/ 0167-8655(88)90063-3
work page 1988
-
[34]
Computational imaging and vision
ter Haar Romeny, B.M.: Front-end vision and multi-scale image analysis : multi-scale computer vision theory and applications, written in Mathe- matica. Computational imaging and vision. CIVI. Kluwer Academic Publishers, Netherlands (2003)
work page 2003
-
[35]
Lecture Notes in Morphogenesis
Petitot, J.: Elements of Neurogeometry: Func- tional Architectures of Vision, 1st edn. Lecture Notes in Morphogenesis. Springer, Cham (2017). https://doi.org/10.1007/978-3-319-65591-8
-
[36]
Duits, R., Boscain, U., Rossi, F., Sachkov, Y.: Association Fields via Cuspless Sub-Riemannian Geodesics inSE(2). Journal of Mathematical Imaging and Vision49, 384–417 (2014) https:// doi.org/10.1007/s10851-013-0475-y
-
[37]
Journal of Mathematical Imaging and Vision65(6), 819–843 (2023) https://doi.org/10
Bellaard, G., Bon, D.L., Pai, G., Smets, B.M., Duits, R.: Analysis of (sub-)Riemannian PDE- G-CNNs. Journal of Mathematical Imaging and Vision65(6), 819–843 (2023) https://doi.org/10. 1007/s10851-023-01147-w
work page 2023
-
[38]
ESAIM: COCV20(3), 748–770 (2014) https://doi.org/10.1051/cocv/2013082
Boscain, U., Duits, R., Rossi, F., Sachkov, Y.: Curve cuspless reconstruction via sub-Riemannian geometry. ESAIM: COCV20(3), 748–770 (2014) https://doi.org/10.1051/cocv/2013082
-
[39]
Quarterly of Applied mathematics, AMS68(2), 255–331 (2010)
Duits, R., Franken, E.M.: Left-invariant parabolic evolutions onSE(2) and contour enhancement via invertible orientation scores. Quarterly of Applied mathematics, AMS68(2), 255–331 (2010)
work page 2010
-
[40]
IJCV72, 79–102 (2007) https: //doi.org/10.1007/s11263-006-8894-5
Duits, R., Felsberg, M., Granlund, G., ter Haar Romeny, B.: Image Analysis and Recon- struction Using a Wavelet Transform Constructed from a Reducible Representation of the Euclidean Motion Group. IJCV72, 79–102 (2007) https: //doi.org/10.1007/s11263-006-8894-5
-
[41]
In: Nielsen, F., Barbaresco, F
Sherry, F.M., van de Geijn, C., Bekkers, E.J., Duits, R.: Orientation Scores Should Be a Piece of Cake. In: Nielsen, F., Barbaresco, F. (eds.) Geometric Science of Information, pp. 224–233. Springer, Cham (2025). https://doi.org/10.1007/ 978-3-032-03924-8 23
work page 2025
-
[42]
PhD thesis, Technis- che Universiteit Eindhoven (2005)
Duits, R.: Perceptual Organization in Image Anal- ysis: A Mathematical Approach Based on Scale, Orientation and Curvature. PhD thesis, Technis- che Universiteit Eindhoven (2005). https://doi. org/10.6100/IR595308
-
[43]
Soft- wareX12, 100605 (2020) https://doi.org/10.1016/ j.softx.2020.100605
van Kessel, L., Hagen, C.W.: Nebula: Monte Carlo simulator of electron–matter interaction. Soft- wareX12, 100605 (2020) https://doi.org/10.1016/ j.softx.2020.100605
-
[44]
Journal of Mathematical Imag- ing and Vision65, 209–239 (2023) https://doi
Smets, B.M.N., Portegies, J., Bekkers, E.J., Duits, R.: PDE-Based Group Equivariant Convolutional Neural Networks. Journal of Mathematical Imag- ing and Vision65, 209–239 (2023) https://doi. org/10.1007/s10851-022-01114-x
-
[45]
Moiseev, I., Sachkov, Y.L.: Maxwell strata in sub- Riemannian problemon the group of motions of a plane. ESAIM: Control, Optimisation and Cal- culus of Variations16(2), 380–399 (2010) https: //doi.org/10.1051/cocv/2009004
-
[46]
Lee, J., Haralick, R., Shapiro, L.: Morphologic edge detection. IEEE Journal on Robotics and Automation3(2), 142–156 (1987) https://doi.org/ 10.1109/JRA.1987.1087088
-
[47]
CVGIP: Graphical Models and Image Process- ing56(1994) https://doi.org/10.1006/cgip.1994.1042
Lee, T.C., Kashyap, R.L., Chu, C.N.: Building Skeleton Models via 3-D Medial Surface Axis Thinning Algorithms. CVGIP: Graphical Models and Image Processing56(6), 462–478 (1994) https: //doi.org/10.1006/cgip.1994.1042 29
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.