A PDE approach to data-driven sub-Riemannian geodesics in SE(2)

E.J. Bekkers, R. Duits, A. Mashtakov, G.R. Sanguinetti

Onderzoeksoutput: Bijdrage aan tijdschriftTijdschriftartikelAcademicpeer review

22 Citaties (Scopus)

Uittreksel

We present a new flexible wavefront propagation algorithm for the boundary value problem for sub-Riemannian (SR) geodesics in the roto-translation group $SE(2) = \mathbb{R}^2 \rtimes S^1$ with a metric tensor depending on a smooth external cost $\mathcal{C}:SE(2) \to [\delta,1]$, $\delta>0$, computed from image data. The method consists of a first step where an SR-distance map is computed as a viscosity solution of a Hamilton--Jacobi--Bellman system derived via Pontryagin's maximum principle (PMP). Subsequent backward integration, again relying on PMP, gives the SR-geodesics. For $\mathcal{C}=1$ we show that our method produces the global minimizers. Comparison with exact solutions shows a remarkable accuracy of the SR-spheres and the SR-geodesics. We present numerical computations of Maxwell points and cusp points, which we again verify for the uniform cost case $\mathcal{C}=1$. Regarding image analysis applications, trackings of elongated structures in retinal and synthetic images show that our line tracking generically deals with crossings. We show the benefits of including the SR-geometry.
TaalEngels
Pagina's2740-2770
Aantal pagina's31
TijdschriftSIAM Journal on Imaging Sciences
Volume8
Nummer van het tijdschrift4
DOI's
StatusGepubliceerd - 1 dec 2015

Vingerafdruk

Maximum principle
Data-driven
Geodesic
Pontryagin Maximum Principle
Wavefronts
Sub-Riemannian Geometry
Image analysis
Boundary value problems
Tensors
Costs
Global Minimizer
Hamilton-Jacobi
Viscosity Solutions
Viscosity
Cusp
Image Analysis
Wave Front
Numerical Computation
Geometry
Tensor

Trefwoorden

    Citeer dit

    Bekkers, E.J. ; Duits, R. ; Mashtakov, A. ; Sanguinetti, G.R./ A PDE approach to data-driven sub-Riemannian geodesics in SE(2). In: SIAM Journal on Imaging Sciences. 2015 ; Vol. 8, Nr. 4. blz. 2740-2770
    @article{44410feb24344d44b14afd9a6334e8ba,
    title = "A PDE approach to data-driven sub-Riemannian geodesics in SE(2)",
    abstract = "We present a new flexible wavefront propagation algorithm for the boundary value problem for sub-Riemannian (SR) geodesics in the roto-translation group $SE(2) = \mathbb{R}^2 \rtimes S^1$ with a metric tensor depending on a smooth external cost $\mathcal{C}:SE(2) \to [\delta,1]$, $\delta>0$, computed from image data. The method consists of a first step where an SR-distance map is computed as a viscosity solution of a Hamilton--Jacobi--Bellman system derived via Pontryagin's maximum principle (PMP). Subsequent backward integration, again relying on PMP, gives the SR-geodesics. For $\mathcal{C}=1$ we show that our method produces the global minimizers. Comparison with exact solutions shows a remarkable accuracy of the SR-spheres and the SR-geodesics. We present numerical computations of Maxwell points and cusp points, which we again verify for the uniform cost case $\mathcal{C}=1$. Regarding image analysis applications, trackings of elongated structures in retinal and synthetic images show that our line tracking generically deals with crossings. We show the benefits of including the SR-geometry.",
    keywords = "roto-translation group, Hamilton-Jacobi equations, vessel tracking, sub-Riemannian geometry, morphological scale spaces",
    author = "E.J. Bekkers and R. Duits and A. Mashtakov and G.R. Sanguinetti",
    year = "2015",
    month = "12",
    day = "1",
    doi = "10.1137/15M1018460",
    language = "English",
    volume = "8",
    pages = "2740--2770",
    journal = "SIAM Journal on Imaging Sciences",
    issn = "1936-4954",
    publisher = "Society for Industrial and Applied Mathematics (SIAM)",
    number = "4",

    }

    A PDE approach to data-driven sub-Riemannian geodesics in SE(2). / Bekkers, E.J.; Duits, R.; Mashtakov, A.; Sanguinetti, G.R.

    In: SIAM Journal on Imaging Sciences, Vol. 8, Nr. 4, 01.12.2015, blz. 2740-2770.

    Onderzoeksoutput: Bijdrage aan tijdschriftTijdschriftartikelAcademicpeer review

    TY - JOUR

    T1 - A PDE approach to data-driven sub-Riemannian geodesics in SE(2)

    AU - Bekkers,E.J.

    AU - Duits,R.

    AU - Mashtakov,A.

    AU - Sanguinetti,G.R.

    PY - 2015/12/1

    Y1 - 2015/12/1

    N2 - We present a new flexible wavefront propagation algorithm for the boundary value problem for sub-Riemannian (SR) geodesics in the roto-translation group $SE(2) = \mathbb{R}^2 \rtimes S^1$ with a metric tensor depending on a smooth external cost $\mathcal{C}:SE(2) \to [\delta,1]$, $\delta>0$, computed from image data. The method consists of a first step where an SR-distance map is computed as a viscosity solution of a Hamilton--Jacobi--Bellman system derived via Pontryagin's maximum principle (PMP). Subsequent backward integration, again relying on PMP, gives the SR-geodesics. For $\mathcal{C}=1$ we show that our method produces the global minimizers. Comparison with exact solutions shows a remarkable accuracy of the SR-spheres and the SR-geodesics. We present numerical computations of Maxwell points and cusp points, which we again verify for the uniform cost case $\mathcal{C}=1$. Regarding image analysis applications, trackings of elongated structures in retinal and synthetic images show that our line tracking generically deals with crossings. We show the benefits of including the SR-geometry.

    AB - We present a new flexible wavefront propagation algorithm for the boundary value problem for sub-Riemannian (SR) geodesics in the roto-translation group $SE(2) = \mathbb{R}^2 \rtimes S^1$ with a metric tensor depending on a smooth external cost $\mathcal{C}:SE(2) \to [\delta,1]$, $\delta>0$, computed from image data. The method consists of a first step where an SR-distance map is computed as a viscosity solution of a Hamilton--Jacobi--Bellman system derived via Pontryagin's maximum principle (PMP). Subsequent backward integration, again relying on PMP, gives the SR-geodesics. For $\mathcal{C}=1$ we show that our method produces the global minimizers. Comparison with exact solutions shows a remarkable accuracy of the SR-spheres and the SR-geodesics. We present numerical computations of Maxwell points and cusp points, which we again verify for the uniform cost case $\mathcal{C}=1$. Regarding image analysis applications, trackings of elongated structures in retinal and synthetic images show that our line tracking generically deals with crossings. We show the benefits of including the SR-geometry.

    KW - roto-translation group

    KW - Hamilton-Jacobi equations

    KW - vessel tracking

    KW - sub-Riemannian geometry

    KW - morphological scale spaces

    U2 - 10.1137/15M1018460

    DO - 10.1137/15M1018460

    M3 - Article

    VL - 8

    SP - 2740

    EP - 2770

    JO - SIAM Journal on Imaging Sciences

    T2 - SIAM Journal on Imaging Sciences

    JF - SIAM Journal on Imaging Sciences

    SN - 1936-4954

    IS - 4

    ER -