Partial functional correspondence Michael Bronstein University of - - PowerPoint PPT Presentation

partial functional correspondence
SMART_READER_LITE
LIVE PREVIEW

Partial functional correspondence Michael Bronstein University of - - PowerPoint PPT Presentation

Partial functional correspondence Michael Bronstein University of Lugano Intel Corporation Lyon, 7 July 2016 1/59 Microsoft Kinect 2010 (Acquired by Intel in 2012) 4/59 Different form factor computers featuring Intel RealSense 3D camera


slide-1
SLIDE 1

1/59

Partial functional correspondence

Michael Bronstein

University of Lugano Intel Corporation Lyon, 7 July 2016

slide-2
SLIDE 2

Microsoft Kinect 2010

slide-3
SLIDE 3

4/59

(Acquired by Intel in 2012)

slide-4
SLIDE 4
slide-5
SLIDE 5

7/59

Different form factor computers featuring Intel RealSense 3D camera

slide-6
SLIDE 6

8/59

Deluge of geometric data

3D sensors Repositories 3D printers

slide-7
SLIDE 7

9/59

Applications

Deformable fusion Motion transfer Motion capture Texture mapping Dou et al. 2015; Sumner, Popovi´ c 2004; Faceshift; Cow image: Moore 2014

slide-8
SLIDE 8

10/59

Shape correspondence problem

Isometric

slide-9
SLIDE 9

10/59

Shape correspondence problem

Isometric Partial

slide-10
SLIDE 10

10/59

Shape correspondence problem

Isometric Partial Different representation

slide-11
SLIDE 11

10/59

Shape correspondence problem

Isometric Partial Different representation Non-isometric

slide-12
SLIDE 12

11/59

Computer Graphics Forum SGP 2016 Computer Graphics Forum SGP 2016 Best paper award

slide-13
SLIDE 13

12/59

Outline

Background: Spectral analysis on manifolds Functional correspondence Partial functional correspondence Non-rigid puzzles

slide-14
SLIDE 14

13/59

Riemannian geometry in one minute

Tangent plane TmM = local Euclidean representation of manifold (surface) M around m

m TmM M

slide-15
SLIDE 15

13/59

Riemannian geometry in one minute

Tangent plane TmM = local Euclidean representation of manifold (surface) M around m Riemannian metric ·, ·TmM : TmM × TmM → R depending smoothly on m

m TmM m′ Tm′M

slide-16
SLIDE 16

13/59

Riemannian geometry in one minute

Tangent plane TmM = local Euclidean representation of manifold (surface) M around m Riemannian metric ·, ·TmM : TmM × TmM → R depending smoothly on m Isometry = metric-preserving shape deformation

m TmM m′ Tm′M

slide-17
SLIDE 17

13/59

Riemannian geometry in one minute

Tangent plane TmM = local Euclidean representation of manifold (surface) M around m Riemannian metric ·, ·TmM : TmM × TmM → R depending smoothly on m Isometry = metric-preserving shape deformation Exponential map expm : TmM → M ‘unit step along geodesic’

m v TmM expm(v)

slide-18
SLIDE 18

14/59

Laplace-Beltrami operator

m f

Smooth field f : M → R

slide-19
SLIDE 19

14/59

Laplace-Beltrami operator

m f f ◦ expm

Smooth field f ◦ expm : TmM → R

slide-20
SLIDE 20

14/59

Laplace-Beltrami operator

Intrinsic gradient ∇Mf(m) = ∇(f ◦ expm)(0) Taylor expansion (f ◦ expm)(v) ≈ f(m) + ∇Mf(m), vTmM

m f f ◦ expm

slide-21
SLIDE 21

14/59

Laplace-Beltrami operator

Intrinsic gradient ∇Mf(m) = ∇(f ◦ expm)(0) Taylor expansion (f ◦ expm)(v) ≈ f(m) + ∇Mf(m), vTmM Laplace-Beltrami operator ∆Mf(m) = ∆(f ◦ expm)(0)

m f f ◦ expm

slide-22
SLIDE 22

14/59

Laplace-Beltrami operator

Intrinsic gradient ∇Mf(m) = ∇(f ◦ expm)(0) Taylor expansion (f ◦ expm)(v) ≈ f(m) + ∇Mf(m), vTmM Laplace-Beltrami operator ∆Mf(m) = ∆(f ◦ expm)(0)

m f f ◦ expm

Intrinsic (expressed solely in terms of the Riemannian metric)

slide-23
SLIDE 23

14/59

Laplace-Beltrami operator

Intrinsic gradient ∇Mf(m) = ∇(f ◦ expm)(0) Taylor expansion (f ◦ expm)(v) ≈ f(m) + ∇Mf(m), vTmM Laplace-Beltrami operator ∆Mf(m) = ∆(f ◦ expm)(0)

m f f ◦ expm

Intrinsic (expressed solely in terms of the Riemannian metric) Isometry-invariant

slide-24
SLIDE 24

14/59

Laplace-Beltrami operator

Intrinsic gradient ∇Mf(m) = ∇(f ◦ expm)(0) Taylor expansion (f ◦ expm)(v) ≈ f(m) + ∇Mf(m), vTmM Laplace-Beltrami operator ∆Mf(m) = ∆(f ◦ expm)(0)

m f f ◦ expm

Intrinsic (expressed solely in terms of the Riemannian metric) Isometry-invariant Self-adjoint ∆Mf, gL2(M)=f, ∆MgL2(M)

slide-25
SLIDE 25

14/59

Laplace-Beltrami operator

Intrinsic gradient ∇Mf(m) = ∇(f ◦ expm)(0) Taylor expansion (f ◦ expm)(v) ≈ f(m) + ∇Mf(m), vTmM Laplace-Beltrami operator ∆Mf(m) = ∆(f ◦ expm)(0)

m f f ◦ expm

Intrinsic (expressed solely in terms of the Riemannian metric) Isometry-invariant Self-adjoint ∆Mf, gL2(M)=f, ∆MgL2(M) ⇒ orthogonal eigenfunctions

slide-26
SLIDE 26

14/59

Laplace-Beltrami operator

Intrinsic gradient ∇Mf(m) = ∇(f ◦ expm)(0) Taylor expansion (f ◦ expm)(v) ≈ f(m) + ∇Mf(m), vTmM Laplace-Beltrami operator ∆Mf(m) = ∆(f ◦ expm)(0)

m f f ◦ expm

Intrinsic (expressed solely in terms of the Riemannian metric) Isometry-invariant Self-adjoint ∆Mf, gL2(M)=f, ∆MgL2(M) ⇒ orthogonal eigenfunctions Positive semidefinite ⇒ non-negative eigenvalues

slide-27
SLIDE 27

15/59

Discrete Laplacian

j i wij i αij βij ai αij

Undirected graph (V, E) (∆f)i ≈

  • (i,j)∈E

wij(fi − fj) Triangular mesh (V, E, F) (∆f)i ≈ 1 ai

  • (i,j)∈E

wij(fi − fj) wij =     

cot αij+cot βij 2

(i, j) ∈ Ei

1 2 cot αij

(i, j) ∈ Eb −

k=i wik

i = j else

ai = local area element

Tutte 1963; MacNeal 1949; Duffin 1959; Pinkall, Polthier 1993

slide-28
SLIDE 28

16/59

Fourier analysis (Euclidean spaces)

A function f : [−π, π] → R can be written as Fourier series f(x) =

  • ω

1 2π π

−π

f(ξ)eiωξdξ e−iωx

α1 α2 α3 = + + + . . .

slide-29
SLIDE 29

16/59

Fourier analysis (Euclidean spaces)

A function f : [−π, π] → R can be written as Fourier series f(x) =

  • ω

1 2π π

−π

f(ξ)eiωξdξ

  • ˆ

f(ω)=f,e−iωxL2([−π,π])

e−iωx

α1 α2 α3 = + + + . . .

slide-30
SLIDE 30

16/59

Fourier analysis (Euclidean spaces)

A function f : [−π, π] → R can be written as Fourier series f(x) =

  • ω

1 2π π

−π

f(ξ)eiωξdξ

  • ˆ

f(ω)=f,e−iωxL2([−π,π])

e−iωx

α1 α2 α3 = + + + . . .

Fourier basis = Laplacian eigenfunctions: ∆e−iωx = ω2e−iωx

slide-31
SLIDE 31

17/59

Fourier analysis (non-Euclidean spaces)

A function f : M → R can be written as Fourier series f(m) =

  • k≥1
  • M

f(ξ)φk(ξ)dξ

  • ˆ

fk=f,φkL2(M)

φk(m)

= α1 + α2 + α3 + . . . f φ1 φ2 φ3

Fourier basis = Laplacian eigenfunctions: ∆Mφk = λkφk

slide-32
SLIDE 32

18/59

Outline

Background: Spectral analysis on manifolds Functional correspondence Partial functional correspondence Non-rigid puzzles

slide-33
SLIDE 33

19/59

Point-wise correspondence

m M n N t

Point-wise maps t: M → N

slide-34
SLIDE 34

19/59

Functional correspondence

f F(M) g F(N) T

Functional maps T: F(M) → F(N)

Ovsjanikov et al. 2012

slide-35
SLIDE 35

20/59

Functional correspondence

f g ↓ T ↓ Ovsjanikov et al. 2012

slide-36
SLIDE 36

20/59

Functional correspondence

f g ≈ a1 + a2 + · · · + ak ≈ b1 + b2 + · · · + bk ↓ T ↓ Ovsjanikov et al. 2012

slide-37
SLIDE 37

20/59

Functional correspondence

f g ≈ a1 + a2 + · · · + ak ≈ b1 + b2 + · · · + bk ↓ T ↓ ↓ C⊤ ↓ Translates Fourier coefficients from Φ to Ψ Ovsjanikov et al. 2012

slide-38
SLIDE 38

20/59

Functional correspondence

f g ≈ a1 + a2 + · · · + ak ≈ b1 + b2 + · · · + bk ↓ T ↓ ↓ C⊤ ↓ Translates Fourier coefficients from Φ to Ψ ≈ Ψk Φ⊤

k

g⊤Ψk = f ⊤ΦkC where Φk = (φ1, . . . , φk), Ψk = (ψ1, . . . , ψk) are truncated Laplace-Beltrami eigenbases

Ovsjanikov et al. 2012

slide-39
SLIDE 39

21/59

Functional correspondence in Laplacian eigenbases

For isometric simple spectrum shapes C is diagonal since ψi = ±Tφi

slide-40
SLIDE 40

22/59

Computing functional correspondence

Ovsjanikov et al. 2012

slide-41
SLIDE 41

22/59

Computing functional correspondence

f 1 f 2

· · ·

f q g1 g2

· · ·

gq

Given ordered set of functions f 1, . . . , f q on M and corresponding functions g1, . . . , gq on N (gi ≈ Tf i)

Ovsjanikov et al. 2012

slide-42
SLIDE 42

22/59

Computing functional correspondence

f 1 f 2

· · ·

f q g1 g2

· · ·

gq

Given ordered set of functions f 1, . . . , f q on M and corresponding functions g1, . . . , gq on N (gi ≈ Tf i) C found by solving a system of qk equations with k2 variables G⊤Ψk = F⊤ΦkC where F = (f 1, . . . , f q) and G = (g1, . . . , gq) are n × q and m × q matrices

Ovsjanikov et al. 2012

slide-43
SLIDE 43

23/59

Key issues

How to recover point-wise correspondence with some guarantees (e.g. bijectivity)? How to automatically find corresponding functions F, G? Near isometric shapes: easy (a lot of structure!) Non-isometric shapes: hard Does not work well in case of missing parts and topological noise

slide-44
SLIDE 44

24/59

Partial Laplacian eigenvectors

φ2 φ3 φ4 φ5 φ6 φ7 φ8 φ9 ψ2 ψ3 ψ4 ψ5 ψ6 ψ7 ψ8 ψ9 ζ2 ζ3 ζ4 ζ5 ζ6 ζ7 ζ8 ζ9

Laplacian eigenvectors of a shape with missing parts (Neumann boundary conditions)

Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-45
SLIDE 45

25/59

Partial Laplacian eigenvectors

Functional correspondence matrix C

Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-46
SLIDE 46

26/59

Perturbation analysis: intuition

∆M ∆M ∆ ¯

M

φ1 φ2 φ3 φ1 φ2 φ3 ¯ φ1 ¯ φ2 ¯ φ3 M ¯ M

Ignoring boundary interaction: disjoint parts (block-diagonal matrix) Eigenvectors = Mixture of eigenvectors of the parts

Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-47
SLIDE 47

27/59

Perturbation analysis: eigenvalues

10 20 30 40 50 0.00 2.00 4.00 6.00 8.00 ·10−2 eigenvalue number r k N M

Slope r

k ≈ |M| |N | (depends on the area of the cut)

Consistent with Weil’s law

Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-48
SLIDE 48

28/59

Perturbation analysis: details

∆M ∆ ¯

M

∆M+tDM ∆ ¯ M+tD ¯ M

tE tE⊤ M ¯ M Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-49
SLIDE 49

28/59

Perturbation analysis: details

∆M ∆ ¯

M

∆M+tDM ∆ ¯ M+tD ¯ M

tE tE⊤ M ¯ M

“How would the Laplacian eigenvalues and eigenvectors of the red part change if we attached a blue part to it?”

Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-50
SLIDE 50

28/59

Perturbation analysis: details

PM P E DM n × n n × ¯ n M ¯ M

“How would the Laplacian eigenvalues and eigenvectors of the red part change if we attached a blue part to it?”

Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-51
SLIDE 51

29/59

Perturbation analysis: details

Denote ∆M + tPM = Φ(t)Λ(t)Φ(t)⊤, ∆ ¯

M = ¯

Φ ¯ Λ ¯ Φ⊤, Φ = Φ(0), and Λ = Λ(0). Theorem 1 (eigenvalues) The derivative of the non-trivial eigenvalues is given by d dtλi = φ⊤

i PMφi

PM = DM

  • Rodol`

a, Cosmo, B, Torsello, Cremers 2016

slide-52
SLIDE 52

29/59

Perturbation analysis: details

Denote ∆M + tPM = Φ(t)Λ(t)Φ(t)⊤, ∆ ¯

M = ¯

Φ ¯ Λ ¯ Φ⊤, Φ = Φ(0), and Λ = Λ(0). Theorem 1 (eigenvalues) The derivative of the non-trivial eigenvalues is given by d dtλi = φ⊤

i PMφi

PM = DM

  • Theorem 2 (eigenvectors) Assuming λi = λj for i = j and λi = ¯

λj for all i, j, the derivative of the non-trivial eigenvectors is given by d dtφi =

n

  • j=1

j=i

φ⊤

i PMφj

λi − λj φj +

¯ n

  • j=1

φ⊤

i P ¯

φj λi − ¯ λj ¯ φj P = E

  • Rodol`

a, Cosmo, B, Torsello, Cremers 2016

slide-53
SLIDE 53

30/59

Perturbation analysis: boundary interaction strength

Value of f

10 20

Eigenvector perturbation depends on length and position of the boundary Perturbation strength d

dtΦF ≤ c

  • ∂M f(m)dm, where

f(m) =

n

  • i,j=1

j=i

φi(m)φj(m) λi − λj 2

Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-54
SLIDE 54

31/59

Laplacian perturbation: typical picture

Plate Punctured plate Figure: Filoche, Mayboroda 2009

slide-55
SLIDE 55

32/59

Partial functional maps

Model shape M Query shape N Part M ⊆ M ≈ isometric to N Data F, G Partial functional map TG ≈ F(M)

M N M T Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-56
SLIDE 56

32/59

Partial functional maps

Model shape M Query shape N Part M ⊆ M ≈ isometric to N Data F, G Partial functional map G⊤ΨC ≈ F(M)⊤Φ

M N M C Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-57
SLIDE 57

32/59

Partial functional maps

Model shape M Query shape N Part M ⊆ M ≈ isometric to N Data F, G Partial functional map G⊤ΨC ≈ F⊤diag(v)Φ v ∈ F(M) indicator function of M

M N v C Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-58
SLIDE 58

32/59

Partial functional maps

Model shape M Query shape N Part M ⊆ M ≈ isometric to N Data F, G Partial functional map G⊤ΨC ≈ F⊤diag(η(v))Φ v ∈ F(M) indicator function of M η(t) = 1

2(tanh(2t − 1) + 1)

M N v C

Optimization problem w.r.t. correspondence C and part v min

C,v G⊤ΨC − F⊤diag(η(v))Φ2,1 + ρcorr(C) + ρpart(v)

Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-59
SLIDE 59

33/59

Partial functional maps

min

C,v G⊤ΨC − F⊤diag(η(v))Φ2,1 + ρcorr(C) + ρpart(v)

Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-60
SLIDE 60

33/59

Partial functional maps

min

C,v G⊤ΨC − F⊤diag(η(v))Φ2,1 + ρcorr(C) + ρpart(v)

Part regularization ρpart(v) = µ1

  • |N| −
  • M

η(v)dm 2 + µ2

  • M

ξ(v)∇Mvdm where ξ(t) ≈ δ

  • η(t) − 1

2

  • Rodol`

a, Cosmo, B, Torsello, Cremers 2016; BB 2008

slide-61
SLIDE 61

33/59

Partial functional maps

min

C,v G⊤ΨC − F⊤diag(η(v))Φ2,1 + ρcorr(C) + ρpart(v)

Part regularization ρpart(v) = µ1

  • |N| −
  • M

η(v)dm 2

  • area preservation

+ µ2

  • M

ξ(v)∇Mvdm

  • Mumford−Shah

where ξ(t) ≈ δ

  • η(t) − 1

2

  • Rodol`

a, Cosmo, B, Torsello, Cremers 2016; BB 2008

slide-62
SLIDE 62

33/59

Partial functional maps

min

C,v G⊤ΨC − F⊤diag(η(v))Φ2,1 + ρcorr(C) + ρpart(v)

Part regularization ρpart(v) = µ1

  • |N| −
  • M

η(v)dm 2

  • area preservation

+ µ2

  • M

ξ(v)∇Mvdm

  • Mumford−Shah

where ξ(t) ≈ δ

  • η(t) − 1

2

  • Correspondence regularization

ρcorr(C) = µ3C ◦ W2

F + µ4

  • i=j

(C⊤C)2

ij + µ5

  • i

((C⊤C)ii − di)2

Rodol` a, Cosmo, B, Torsello, Cremers 2016; BB 2008

slide-63
SLIDE 63

33/59

Partial functional maps

min

C,v G⊤ΨC − F⊤diag(η(v))Φ2,1 + ρcorr(C) + ρpart(v)

Part regularization ρpart(v) = µ1

  • |N| −
  • M

η(v)dm 2

  • area preservation

+ µ2

  • M

ξ(v)∇Mvdm

  • Mumford−Shah

where ξ(t) ≈ δ

  • η(t) − 1

2

  • Correspondence regularization

ρcorr(C) = µ3C ◦ W2

F

  • slant

+ µ4

  • i=j

(C⊤C)2

ij

  • ≈ orthogonality

+ µ5

  • i

((C⊤C)ii − di)2

  • rank≈r

Rodol` a, Cosmo, B, Torsello, Cremers 2016; BB 2008

slide-64
SLIDE 64

34/59

Structure of partial functional correspondence

C W C⊤C

20 40 60 80 100 2 4

singular values Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-65
SLIDE 65

35/59

Alternating minimization

C-step: fix v∗, solve for correspondence C min

C G⊤ΨC − F⊤diag(η(v∗))Φ2,1 + ρcorr(C)

v-step: fix C∗, solve for part v min

v G⊤ΨC∗ − F⊤diag(η(v))Φ2,1 + ρpart(v)

Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-66
SLIDE 66

35/59

Alternating minimization

C-step: fix v∗, solve for correspondence C min

C G⊤ΨC − F⊤diag(η(v∗))Φ2,1 + ρcorr(C)

v-step: fix C∗, solve for part v min

v G⊤ΨC∗ − F⊤diag(η(v))Φ2,1 + ρpart(v)

Iteration 1 2 3 4 Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-67
SLIDE 67

36/59

Example of convergence

20 40 60 80 100 104 105 106 107 108 109 1010 Iteration Energy C-step v-step 5 10 15 20 25 Time (sec.) Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-68
SLIDE 68

37/59

Examples of partial functional maps

Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-69
SLIDE 69

37/59

Examples of partial functional maps

Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-70
SLIDE 70

37/59

Examples of partial functional maps

Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-71
SLIDE 71

37/59

Examples of partial functional maps

Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-72
SLIDE 72

38/59

Partial functional maps vs Functional maps

0.05 0.1 0.15 0.2 0.25 20 40 60 80 100 50 100 150 50 100 150 Geodesic error % Correspondences

PFM

  • Func. maps

Correspondence performance for different rank values k

Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-73
SLIDE 73

39/59

Partial correspondence performance

0.05 0.1 0.15 0.2 0.25 20 40 60 80 100

Geodesic Error % Correspondences

Cuts

0.05 0.1 0.15 0.2 0.25

Geodesic Error

Holes PFM RF IM EN GT

SHREC’16 Partial Matching benchmark Rodol` a et al. 2016; Methods: Rodol` a, Cosmo, B, Torsello, Cremers 2016 (PFM); Sahillio˘ glu, Yemez 2012 (IM); Rodol` a, Bronstein, Albarelli, Bergamasco, Torsello 2012 (GT); Rodol` a et al. 2013 (EN); Rodol` a et al. 2014 (RF)

slide-74
SLIDE 74

40/59

Partial correspondence performance

20 40 60 80 0.2 0.4 0.6 0.8 1

Partiality (%) Mean geodesic error

Cuts

20 40 60 80

Partiality (%)

Holes PFM RF IM EN GT

SHREC’16 Partial Matching benchmark Rodol` a et al. 2016; Methods: Rodol` a, Cosmo, B, Torsello, Cremers 2016 (PFM); Sahillio˘ glu, Yemez 2012 (IM); Rodol` a, Bronstein, Albarelli, Bergamasco, Torsello 2012 (GT); Rodol` a et al. 2013 (EN); Rodol` a et al. 2014 (RF)

slide-75
SLIDE 75

41/59

Deep learning + Partial functional maps

Correspondence Correspondence error

0.0 0.1

Boscaini, Masci, Rodol` a, B 2016

slide-76
SLIDE 76

42/59

Deep learning + Partial functional maps

Correspondence Correspondence error

0.0 0.1

Boscaini, Masci, Rodol` a, B 2016

slide-77
SLIDE 77

43/59

Outline

Background: Spectral analysis on manifolds Functional correspondence Partial functional correspondence Non-rigid puzzles

slide-78
SLIDE 78

Litani, BB 2012

slide-79
SLIDE 79

45/59

Partial correspondence

Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-80
SLIDE 80

45/59

Non-rigid puzzle

Litani, Rodol` a, BB, Cremers 2016

slide-81
SLIDE 81

46/59

Partial Laplacian eigenvectors

Functional correspondence matrix C

Rodol` a, Cosmo, B, Torsello, Cremers 2016

slide-82
SLIDE 82

47/59

Key observation

N N M M

CNN slant ∝ |N| |N| CMM slant ∝ |M| |M|

Litani, Rodol` a, BB, Cremers 2016

slide-83
SLIDE 83

47/59

Key observation

N N M M

CNM = CNN CN MCMM slant ∝ |N| |N| |M| |M|

Litani, Rodol` a, BB, Cremers 2016

slide-84
SLIDE 84

47/59

Key observation

N N M M

CNM = CNN CN MCMM slant ∝ |N| |N| |M| |M| = |N| |M|

Litani, Rodol` a, BB, Cremers 2016

slide-85
SLIDE 85

48/59

Non-rigid puzzles problem formulation

Input Model M Parts N1, . . . , Np Output Segmentation Mi ⊆ M Located parts Ni ⊆ Ni Correspondences Ci Clutter N c

i

Missing parts M0

Model Parts

M1 M2

C1 C2

N2 N1 N c

2

N c

1

M0 M N2 N1

Litani, Rodol` a, BB, Cremers 2016

slide-86
SLIDE 86

48/59

Non-rigid puzzles problem formulation

Data Fi, Gi Model basis Φ, Φ(Mi) Part bases Ψi, Ψi(Ni) Data term F⊤

i Φ(Mi) ≈ G⊤ i Ψi(Ni)Ci

Model Parts

M1 M2

C1 C2

N2 N1 N c

2

N c

1

M0 M N2 N1

Litani, Rodol` a, BB, Cremers 2016

slide-87
SLIDE 87

49/59

Non-rigid puzzles problem formulation

min

Ci Mi⊆M,Ni⊆Ni p

  • i=1

G⊤

i Ψ

Ψ Ψi(Ni)Ci − F⊤

i Φ

Φ Φ(Mi)2,1 + λM

p

  • i=0

ρpart(Mi) + λN

p

  • i=1

ρpart(Ni) + λcorr

p

  • i=1

ρcorr(Ci) s.t. Mi ∩ Mj = ∅ ∀i = j M0 ∪ M1 ∪ · · · ∪ Mp = M |Mi| = |Ni| ≥ α|Ni|,

Litani, Rodol` a, BB, Cremers 2016

slide-88
SLIDE 88

49/59

Non-rigid puzzles problem formulation

min

Ci ui,vi p

  • i=1

G⊤

i diag(η(ui))Ψ

Ψ ΨiCi − F⊤

i diag(η(vi))Φ

Φ Φ2,1 + λM

p

  • i=0

ρpart(η(vi)) + λN

p

  • i=1

ρpart(η(ui)) + λcorr

p

  • i=1

ρcorr(Ci) s.t.

p

  • i=1

η(ui) = 1 a⊤

Mui = a⊤ N vi ≥ αa⊤ Ni1

Litani, Rodol` a, BB, Cremers 2016

slide-89
SLIDE 89

50/59

Convergence example

Outer iteration 1

Litani, Rodol` a, BB, Cremers 2016

slide-90
SLIDE 90

50/59

Convergence example

Outer iteration 2

Litani, Rodol` a, BB, Cremers 2016

slide-91
SLIDE 91

50/59

Convergence example

Outer iteration 3

Litani, Rodol` a, BB, Cremers 2016

slide-92
SLIDE 92

51/59

Convergence example

80 90 100 110 120 130 140 150 160

Iteration number Time (sec)

30 32 34 36 38 40 42 44 46 48

Litani, Rodol` a, BB, Cremers 2016

slide-93
SLIDE 93

52/59

“Perfect puzzle” example

Model/Part Synthetic (TOSCA) Transformation Isometric Clutter No Missing part No Data term Dense (SHOT)

Litani, Rodol` a, BB, Cremers 2016

slide-94
SLIDE 94

52/59

“Perfect puzzle” example

Model/Part Synthetic (TOSCA) Transformation Isometric Clutter No Missing part No Data term Dense (SHOT) Segmentation

Litani, Rodol` a, BB, Cremers 2016

slide-95
SLIDE 95

52/59

“Perfect puzzle” example

Model/Part Synthetic (TOSCA) Transformation Isometric Clutter No Missing part No Data term Dense (SHOT) Correspondence

Litani, Rodol` a, BB, Cremers 2016

slide-96
SLIDE 96

53/59

Overlapping parts example

Model/Part Synthetic (FAUST) Transformation Near-isometric Clutter Yes (overlap) Missing part No Data term Dense (SHOT) Segmentation

Litani, Rodol` a, BB, Cremers 2016

slide-97
SLIDE 97

53/59

Overlapping parts example

Model/Part Synthetic (FAUST) Transformation Near-isometric Clutter Yes (overlap) Missing part No Data term Dense (SHOT) Correspondence

Litani, Rodol` a, BB, Cremers 2016

slide-98
SLIDE 98

53/59

Overlapping parts example

Model/Part Synthetic (FAUST) Transformation Near-isometric Clutter Yes (overlap) Missing part No Data term Dense (SHOT)

0.0 0.1

Correspondence error

Litani, Rodol` a, BB, Cremers 2016

slide-99
SLIDE 99

54/59

Missing parts example

Model/Part Synthetic (TOSCA) Transformation Isometric Clutter Yes (extra part) Missing part Yes Data term Dense (SHOT)

Litani, Rodol` a, BB, Cremers 2016

slide-100
SLIDE 100

54/59

Missing parts example

Model/Part Synthetic (TOSCA) Transformation Isometric Clutter Yes (extra part) Missing part Yes Data term Dense (SHOT) Segmentation

Litani, Rodol` a, BB, Cremers 2016

slide-101
SLIDE 101

54/59

Missing parts example

Model/Part Synthetic (TOSCA) Transformation Isometric Clutter Yes (extra part) Missing part Yes Data term Dense (SHOT) Correspondence

Litani, Rodol` a, BB, Cremers 2016

slide-102
SLIDE 102

55/59

Scanned data example

Model/Part Synthetic (TOSCA) / Scan Transformation Non-Isometric Clutter No Missing part No Data term Sparse deltas

Litani, Rodol` a, BB, Cremers 2016

slide-103
SLIDE 103

55/59

Scanned data example

Model/Part Synthetic (TOSCA) / Scan Transformation Non-Isometric Clutter No Missing part No Data term Sparse deltas Segmentation

Litani, Rodol` a, BB, Cremers 2016

slide-104
SLIDE 104

56/59

Non-rigid puzzle vs Partial functional map

Partial functional map (pair-wise) Non-rigid puzzle

Rodol` a, Cosmo, B, Torsello, Cremers 2016; Litani, Rodol` a, BB, Cremers 2016

slide-105
SLIDE 105

57/59

Summary

New insights about spectral properties of Laplacians Extension of functional correspondence framework to the partial setting Practically working methods for challenging shape correspondence settings Code available (SGP Reproducibility Stamp) Some over-engineering - can be done simpler! (stay tuned...)

slide-106
SLIDE 106

58/59

  • D. Cremers

Supported by

  • E. Rodolà
  • O. Litany
  • L. Cosmo
  • A. Torsello
  • A. Bronstein
slide-107
SLIDE 107

59/59

Thank you!