Some recent methods in non-rigid shape matching, with and without - - PowerPoint PPT Presentation

some recent methods in non rigid shape matching with and
SMART_READER_LITE
LIVE PREVIEW

Some recent methods in non-rigid shape matching, with and without - - PowerPoint PPT Presentation

Some recent methods in non-rigid shape matching, with and without learning GAMES 2019 webinar Maks Ovsjanikov Based on joint work with: E. Corman, Michael Bronstein, Emanuele Rodol, Justin Solomon, Adrien Butscher, Mirela Ben-Chen, Leonidas


slide-1
SLIDE 1

Some recent methods in non-rigid shape matching, with and without learning

GAMES 2019 webinar

Maks Ovsjanikov

Based on joint work with: E. Corman, Michael Bronstein, Emanuele Rodolà, Justin Solomon, Adrien Butscher, Mirela Ben-Chen, Leonidas Guibas, Frederic Chazal ….

slide-2
SLIDE 2

My Background

  • 2005 – 2010: PhD from Stanford University (advisor

Leonidas Guibas).

  • 2011: engineer at Google Inc.
  • Since 2012 Professor in the Computer Science Department at

Ecole Polytechnique in France.

slide-3
SLIDE 3

University

  • Ecole Polytechnique: located very near Paris.
  • Ranked 2nd best small university in the world in 2019.
  • Est. 1794. Students and professors such as Ampère, Cauchy,

Fourier, Hermite, Lagrange, Monge, Poincaré, Poisson…

  • Very international campus.
slide-4
SLIDE 4

Our group

  • Currently 5 PhD students and 2 PostDocs.
  • Part of the larger STREAM team dedicated to visual

computing with 3 other professors.

  • Many international collaborations: Stanford, MIT,

UCL, KAUST, Univ. Toronto, Univ. Rome, etc.

  • Funding for PhD students and postdocs !
slide-5
SLIDE 5

General Overview

Overall Objective: Create tools for computing and analyzing mappings between geometric objects.

slide-6
SLIDE 6

Talk Overview

Course Notes: [Related] Course Website:

http://www.lix.polytechnique.fr/~maks/fmaps_SIG17_course/

Linked from the website. Or use Attention: (significantly) more material than in the lectures

Sample Code:

See Sample Code link on the website.

  • r http://bit.do/fmaps2017

http://bit.do/fmaps2017_notes

Demo code for basic operations.

6

slide-7
SLIDE 7

Talk Overview

Motivation and Problem Taxonomy Rigid Matching: ICP Functional Map representation, properties Open problems, Q&A Basic pipeline for non-rigid matching Recent extensions, improvements

slide-8
SLIDE 8

What is a Shape?

Discrete: a graph embedded in 3D (triangle mesh). Continuous: a surface embedded in 3D.

  • Connected.
  • Manifold.
  • Without Boundary.

Common assumptions:

slide-9
SLIDE 9

What is a Shape?

5k – 200k triangles

Shapes from the FAUST, SCAPE, and TOSCA datasets

Discrete: a graph embedded in 3D (triangle mesh). Continuous: a surface embedded in 3D.

slide-10
SLIDE 10

Overall Goal

Given two shapes, find correspondences between them.

slide-11
SLIDE 11

Overall Goal

Given two shapes, find correspondences between them. Finding the best map between a pair of shapes.

slide-12
SLIDE 12

Problem Taxonomy

Local vs. Global

refinement (e.g. ICP) | alignment (search) . .

Rigid vs. Deformable

rotation, translation | general deformation.

Semi vs. Fully Automatic

given landmarks | a priori model

Learning-Based vs. Direct

known examples | unseen data

slide-13
SLIDE 13

Problem Taxonomy

Local vs. Global

refinement (e.g. ICP) | alignment (search) . .

Rigid vs. Deformable

rotation, translation | general deformation.

Semi vs. Fully Automatic

given landmarks | a priori model

Learning-Based vs. Direct

known examples | unseen data

slide-14
SLIDE 14

Why Shape Matching?

Given a correspondence, we can transfer:

texture and parametrization segmentation and labels deformation

Other applications: shape interpolation, reconstruction ...

Sumner et al. ‘04.

14

slide-15
SLIDE 15

Rigid Shape Matching

  • The unknowns are the rotation/translation

parameters of the source onto the target shape.

  • Given a pair of shapes, find the optimal Rigid

Alignment between them.

slide-16
SLIDE 16

Iterative Closest Point (ICP)

  • Classical approach: iterate between finding

correspondences and finding the transformation:

example in 2D

M

N

Given a pair of shapes, and , iterate:

  • 1. For each

find nearest neighbor .

  • 2. Find optimal transformation

minimizing:

arg min

R,t

  • i

∥Rxi + t − yi∥2

2

M

N

xi ∈ M yi ∈ N

Besl, McKay (1992). "A Method for Registration of 3-D Shapes".

slide-17
SLIDE 17

Iterative Closest Point

  • Classical approach: iterate between finding

correspondences and finding the transformation:

M

N

Given a pair of shapes, and , iterate:

  • 1. For each

find nearest neighbor .

  • 2. Find optimal transformation

minimizing:

arg min

R,t

  • i

∥Rxi + t − yi∥2

2

M

N

xi ∈ M yi ∈ N

slide-18
SLIDE 18

Iterative Closest Point

  • Classical approach: iterate between finding

correspondences and finding the transformation:

M

N

Given a pair of shapes, and , iterate:

  • 1. For each

find nearest neighbor .

  • 2. Find optimal transformation

minimizing:

arg min

R,t

  • i

∥Rxi + t − yi∥2

2

M

N

xi ∈ M yi ∈ N

slide-19
SLIDE 19

Iterative Closest Point

  • Classical approach: iterate between finding

correspondences and finding the transformation:

Given a pair of shapes, and , iterate:

  • 1. For each

find nearest neighbor .

  • 2. Find optimal transformation

minimizing:

arg min

R,t

  • i

∥Rxi + t − yi∥2

2

M

N

xi ∈ M yi ∈ N

M

N

slide-20
SLIDE 20

Iterative Closest Point

  • Classical approach: iterate between finding

correspondences and finding the transformation:

Given a pair of shapes, and , iterate:

  • 1. For each

find nearest neighbor .

  • 2. Find optimal transformation

minimizing:

arg min

R,t

  • i

∥Rxi + t − yi∥2

2

M

N

xi ∈ M yi ∈ N

M

N

slide-21
SLIDE 21
  • 1. Finding nearest neighbors: can be done with space-

partitioning data structures (e.g., KD-tree).

  • 2. Finding the optimal transformation

minimizing:

Iterative Closest Point

  • Classical approach: iterate between finding

correspondences and finding the transformation:

Can be done efficiently via SVD decomposition.

arg min

R,t

  • i

∥Rxi + t − yi∥2

2

M

N

Arun et al., Least- Squares Fitting of Two 3-D Point Sets

slide-22
SLIDE 22

Non-Rigid Shape Matching

Unlike rigid matching with rotation/translation, there is no compact representation to optimize for in non-rigid matching.

22

slide-23
SLIDE 23

Non-Rigid Shape Matching

What does it mean for a correspondence to be “good”? How to compute it efficiently in practice? Main Questions:

23

slide-24
SLIDE 24

Isometric Shape Matching

Good maps must preserve geodesic distances. Deformation Model:

Geodesic: length of shortest path lying entirely on the surface.

dM(x, y)

dN (T(x), T(y))

M

N

24

slide-25
SLIDE 25

Isometric Shape Matching

Approach:

Find the point mapping by minimizing the distance distortion: The unknowns are point correspondences.

Topt = arg min

T

  • x,y

∥dM(x, y) − dN (T(x), T(y))∥

dM(x, y)

dN (T(x), T(y))

M

N

slide-26
SLIDE 26

Isometric Shape Matching

Approach:

The space of possible solutions is highly non-linear, non-convex.

Problem:

Find the point mapping by minimizing the distance distortion:

Topt = arg min

T

  • x,y

∥dM(x, y) − dN (T(x), T(y))∥

dM(x, y)

dN (T(x), T(y))

M

N

slide-27
SLIDE 27

Functional Map Representation

We would like to define a representation of shape maps that is more amenable to direct optimization.

1. A compact representation for “natural” maps. 2. Inherently global and multi-scale. 3. Handles uncertainty and ambiguity gracefully. 4. Allows efficient manipulations (averaging, composition). 5. Leads to simple (linear) optimization problems.

27

slide-28
SLIDE 28

Functional Approach to Mappings

Given two shapes and a pointwise map The map induces a functional correspondence: TF (f) = g, where g = f ◦ T

T : N → M

M N T T

28

Functional maps: a flexible representation of maps between shapes, O., Ben-Chen, Solomon, Butscher, Guibas, SIGGRAPH 2012

slide-29
SLIDE 29

Functional Approach to Mappings

f : M → R

TF

TF (f) = g : N → R

The map induces a functional correspondence: T

29

TF (f) = g, where g = f ◦ T

Functional maps: a flexible representation of maps between shapes, O., Ben-Chen, Solomon, Butscher, Guibas, SIGGRAPH 2012

Given two shapes and a pointwise map T : N → M

slide-30
SLIDE 30

Functional Approach to Mappings

f : M → R

TF

TF (f) = g : N → R

The map induces a functional correspondence: T

30

TF (f) = g, where g = f ◦ T

Functional maps: a flexible representation of maps between shapes, O., Ben-Chen, Solomon, Butscher, Guibas, SIGGRAPH 2012

Given two shapes and a pointwise map T : N → M

slide-31
SLIDE 31

Functional Approach to Mappings

The induced functional correspondence is linear:

f : M → R

TF

TF (f) = g : N → R

TF (α1f1 + α2f2) = α1TF (f1) + α2TF (f2)

31

Given two shapes and a pointwise map T : N → M

slide-32
SLIDE 32

Functional Map Representation

The induced functional correspondence is complete.

f : M → R

TF

TF (f) = g : N → R

32

Given two shapes and a pointwise map T : N → M

slide-33
SLIDE 33

Observation

Express both and in terms of basis functions:

f

TF (f)

Since is linear, there is a linear transformation from to .

TF

{ai}

{bj}

M

f : M → R

g : N → R

TF

N

f =

  • i

aiφM

i

Assume that both: f ∈ L2(M), g ∈ L2(N)

g = TF (f) =

  • j

bjφN

j

slide-34
SLIDE 34

Functional Map Representation

Eigenfunctions of the Laplace-Beltrami operator: Generalization of Fourier bases to surfaces. Ordered by eigenvalues and provide a natural notion of scale.

λ0 = 0 λ1 = 2.6 λ2 = 3.4 λ3 = 5.1 λ4 = 7.6

∆φi = λiφi

Choice of Basis:

34

∆(f) = −div∇(f)

slide-35
SLIDE 35

Functional Map Representation

Eigenfunctions of the Laplace-Beltrami operator: Generalization of Fourier bases to surfaces. Form an orthonormal basis for . Ordered by eigenvalues and provide a natural notion of scale.

∆φi = λiφi

Choice of Basis:

Can be computed efficiently, with a sparse matrix eigensolver.

35

L2(M)

slide-36
SLIDE 36

Observation

Express both and in terms of basis functions:

f TF (f)

Since is linear, there is a linear transformation from to .

TF

{ai}

{bj}

M

f : M → R

g : N → R

TF

N

f =

  • i

aiφM

i

g = TF (f) =

  • i

biφN

i

slide-37
SLIDE 37

Functional Map Representation

Since the functional mapping TF is linear:

TF can be represented as a matrix C, given a choice of basis for

function spaces.

TF (α1f1 + α2f2) = α1TF (f1) + α2TF (f2)

37 Functional maps: a flexible representation of maps between shapes, O., Ben-Chen, Solomon, Butscher, Guibas, SIGGRAPH 2012

TF

=

slide-38
SLIDE 38

Functional Map Definition

Functional map: matrix C that translates coefficients from to .

ΦM

ΦN

38

slide-39
SLIDE 39

Example Maps in a Reduced Basis

Triangle meshes with pre-computed pointwise maps “Good” maps are close to being diagonal

39

(a) source (b) ground-truth map (c) left-to-right map (d) head-to-tail map

Try fmap_computation_demo on the course website

slide-40
SLIDE 40

Reconstructing from LB basis

Map reconstruction error using a fixed size matrix.

40

0.5 1 1.5 2 2.5 3 3.5 4 4.5

reconstruction error

Number of basis (eigen)-functions 27.9k vertices

Try fmap_reconstru ction_demo on the course website

slide-41
SLIDE 41

Functional Map algebra

  • 1. Map composition becomes matrix multiplication.
  • 2. Map inversion is matrix inversion (in fact, transpose).
  • 3. Algebraic operations on functional maps are possible.

E.g. interpolating between two maps with

C = αC1 +(1−α)C2.

41

slide-42
SLIDE 42

Talk Overview

Motivation and Problem Taxonomy Rigid Matching: ICP Functional Map Representation, properties Open Problems, Q&A Basic pipeline for non-rigid matching Extensions, Improvements

slide-43
SLIDE 43

In practice we do not know C. Given two objects our goal is to find the correspondence. How can the functional representation help to compute the map in practice?

Shape Matching

?

slide-44
SLIDE 44

Matching via Function Preservation

where Given enough pairs, we can recover C through a linear least squares system.

f =

i aiφM i ,

g =

i biφN i .

{a, b}

Suppose we don’t know C. However, we expect a pair of functions and to correspond. Then, C must be s.t.

Ca ≈ b

f : M → R

g : N → R

slide-45
SLIDE 45

Function preservation constraint is general and includes:

  • Attribute (e.g., color) preservation.
  • Descriptor preservation (e.g. Gauss curvature).
  • Landmark correspondences (e.g. distance to the point).
  • Part correspondences (e.g. indicator function).

Map Constraints

Suppose we don’t know C. However, we expect a pair of functions and to correspond. Then, C must be s.t.

Ca ≈ b

f : M → R

g : N → R

slide-46
SLIDE 46

Commutativity Constraints

Regularizations: Commutativity with other operators:

C

Note that the energy: is quadratic in C.

SM SN

CSM = SN C

∥CSM − SN C∥2

F

slide-47
SLIDE 47

Regularization

Lemma 1:

The mapping is isometric, if and only if the functional map matrix commutes with the Laplacian:

Implies that exact isometries result in diagonal functional maps.

C∆M = ∆N C

Functional maps: a flexible representation of maps between shapes, O., Ben-Chen, Solomon, Butscher, Guibas, SIGGRAPH 2012

Linking functional and point-to-point maps

slide-48
SLIDE 48

Basic Pipeline

Given a pair of shapes :

  • 1. Compute the first k (~80-100) eigenfunctions of the Laplace-

Beltrami operator. Store them in matrices:

  • 2. Compute descriptor functions (e.g., Wave Kernel Signature)
  • n . Express them in , as columns of :
  • 3. Solve
  • 4. Convert the functional map

to a point to point map T.

Copt

diagonal matrices of eigenvalues

  • f LB operator

M, N

ΦM, ΦN

ΦM, ΦN

∆M, ∆N :

Copt = arg min

C

∥CA − B∥2 + ∥C∆M − ∆N C∥2

A, B

M, N

slide-49
SLIDE 49

Recent Implementation

Recent implementation incorporating efficient spatial and spectral constraints.

https://github.com/llorz/SGA18_orientation_BCICP_code

Continuous and Orientation-preserving Correspondences via Functional Maps Jing Ren, Adrien Poulenard, Peter Wonka, Maks Ovsjanikov, SIGGRAPH Asia 2018

slide-50
SLIDE 50

Results

A very simple method that puts together many constraints and uses 100 basis functions gives reasonable results:

50 Functional maps: a flexible representation of maps between shapes, O., Ben-Chen, Solomon, Butscher, Guibas, SIGGRAPH 2012

slide-51
SLIDE 51

Results

radius 0.025 radius 0.05

Functional maps: a flexible representation of maps between shapes, O., Ben-Chen, Solomon, Butscher, Guibas, SIGGRAPH 2012

A very simple method that puts together many constraints and uses 100 basis functions gives reasonable results:

51

slide-52
SLIDE 52

Segmentation Transfer without P2P

To transfer functions we do not need to convert functional to pointwise maps. E.g. we can also transfer segmentations: for each segment, transfer its indicator function, and for each point pick the segment that gave the highest value.

slide-53
SLIDE 53

Talk Overview

Motivation and Problem Taxonomy Rigid Matching: ICP Functional Map Representation, properties Open Problems, Q&A Basic pipeline for non-rigid matching Extensions, Improvements

slide-54
SLIDE 54

Some Recent Extensions

Efficient Refinement Unsupervised Learning

ZoomOut: Spectral Upsampling for Efficient Shape Correspondence Melzi, Ren, Rodolà, Sharma, Wonka, Ovsjanikov, SIGGRAPH Asia 2019 Unsupervised Deep Learning for Structured Shape Matching Roufosse, Sharma, Ovsjanikov, ICCV, 2019 (oral).

slide-55
SLIDE 55

Main Question

55

What happens if the descriptors are bad?

slide-56
SLIDE 56

56

Learning approach to computing descriptors.

  • O. Litany, T. Remez, E. Rodolà, A. Bronstein, M. Bronstein: Deep functional maps:

Structured prediction for dense shape correspondence. In Proc. ICCV (2017).

D1

<latexit sha1_base64="/NqoDpuSU5l6AFS0+IbhLzW/c/k=">A CE3icbVDLSsNAFJ34rPFVdekmWAQRKYkIuizahRuxon1AE8pkepMOnUzCzEQo Z/gUv0Yd+LWD/Bb3Dhps9DWAwOHc+69c+/xE0alsu0vY2FxaXl tbRmrm9sbm2Xd3ZbMk4FgSaJWSw6PpbA KIemo pBJxGAI59B2x9e5X7 EYSkMX9QowS8CIecBpRgpaX7es/plSt21Z7AmidOQSqoQKNX/nb7MUkj4IowLGX sRPlZVgoShiMT eVkGAyxCF0NeU4Aul k1XH1qFW+lYQC/24sibq74 MR1KOIl9XRlgN5KyXi/953VQF 15GeZIq4GT6UZAyS8VWfrfVpwKIYiN MBFU72qRARaYKJ2O6dZB3yLgRs+9TUBgFYvjzMUijCgf69tC9yRnpk7Lmc1mnrROq45d e7OKrXLIrcS2kcH6Ag56BzV0DVqoCYiKERP6AW9Gs/Gm/FufExLF4yiZw/9gfH5AyFcnkg=</latexit> <latexit sha1_base64="/NqoDpuSU5l6AFS0+IbhLzW/c/k=">A CE3icbVDLSsNAFJ34rPFVdekmWAQRKYkIuizahRuxon1AE8pkepMOnUzCzEQo Z/gUv0Yd+LWD/Bb3Dhps9DWAwOHc+69c+/xE0alsu0vY2FxaXl tbRmrm9sbm2Xd3ZbMk4FgSaJWSw6PpbA KIemo pBJxGAI59B2x9e5X7 EYSkMX9QowS8CIecBpRgpaX7es/plSt21Z7AmidOQSqoQKNX/nb7MUkj4IowLGX sRPlZVgoShiMT eVkGAyxCF0NeU4Aul k1XH1qFW+lYQC/24sibq74 MR1KOIl9XRlgN5KyXi/953VQF 15GeZIq4GT6UZAyS8VWfrfVpwKIYiN MBFU72qRARaYKJ2O6dZB3yLgRs+9TUBgFYvjzMUijCgf69tC9yRnpk7Lmc1mnrROq45d e7OKrXLIrcS2kcH6Ag56BzV0DVqoCYiKERP6AW9Gs/Gm/FufExLF4yiZw/9gfH5AyFcnkg=</latexit> <latexit sha1_base64="/NqoDpuSU5l6AFS0+IbhLzW/c/k=">A CE3icbVDLSsNAFJ34rPFVdekmWAQRKYkIuizahRuxon1AE8pkepMOnUzCzEQo Z/gUv0Yd+LWD/Bb3Dhps9DWAwOHc+69c+/xE0alsu0vY2FxaXl tbRmrm9sbm2Xd3ZbMk4FgSaJWSw6PpbA KIemo pBJxGAI59B2x9e5X7 EYSkMX9QowS8CIecBpRgpaX7es/plSt21Z7AmidOQSqoQKNX/nb7MUkj4IowLGX sRPlZVgoShiMT eVkGAyxCF0NeU4Aul k1XH1qFW+lYQC/24sibq74 MR1KOIl9XRlgN5KyXi/953VQF 15GeZIq4GT6UZAyS8VWfrfVpwKIYiN MBFU72qRARaYKJ2O6dZB3yLgRs+9TUBgFYvjzMUijCgf69tC9yRnpk7Lmc1mnrROq45d e7OKrXLIrcS2kcH6Ag56BzV0DVqoCYiKERP6AW9Gs/Gm/FufExLF4yiZw/9gfH5AyFcnkg=</latexit> <latexit sha1_base64="/NqoDpuSU5l6AFS0+IbhLzW/c/k=">A CE3icbVDLSsNAFJ34rPFVdekmWAQRKYkIuizahRuxon1AE8pkepMOnUzCzEQo Z/gUv0Yd+LWD/Bb3Dhps9DWAwOHc+69c+/xE0alsu0vY2FxaXl tbRmrm9sbm2Xd3ZbMk4FgSaJWSw6PpbA KIemo pBJxGAI59B2x9e5X7 EYSkMX9QowS8CIecBpRgpaX7es/plSt21Z7AmidOQSqoQKNX/nb7MUkj4IowLGX sRPlZVgoShiMT eVkGAyxCF0NeU4Aul k1XH1qFW+lYQC/24sibq74 MR1KOIl9XRlgN5KyXi/953VQF 15GeZIq4GT6UZAyS8VWfrfVpwKIYiN MBFU72qRARaYKJ2O6dZB3yLgRs+9TUBgFYvjzMUijCgf69tC9yRnpk7Lmc1mnrROq45d e7OKrXLIrcS2kcH6Ag56BzV0DVqoCYiKERP6AW9Gs/Gm/FufExLF4yiZw/9gfH5AyFcnkg=</latexit>

D2

<latexit sha1_base64="4fQo4Jpdl8nw3AmI+ XyHse7Qqg=">A CE3icbVDLSsNAFJ3UV42vqks3wSKISEmKoMuiXbgRK9oHNKFMprfp0Mk zEyE voJLtWPcSdu/QC/xY2TNgt PTBwO feO/ceP2ZUKtv+MgpLy ura8V1c2Nza3untLvXklEiCDRJxCLR8bE Rjk0FVUMOrEAHPoM2v7oKvPbjyAkjfiDGsfghTjgdEAJVlq6r/eqvVLZrthTWIvEyUkZ5Wj0St9uPyJ CFwRhqXsOnasvBQLRQmDiekmEmJMRjiArqYchyC9dLrqxDrS t8aREI/rqyp+rsjxaGU49DXlSFWQznvZeJ/XjdRgwsvpTxOFHAy+2iQMEtFVna31acCiGJj T ARVO9qkSEWmCidjunWQd8i4EbPvY1BYBWJk9TFIg pn+jbAvc0Y6ZOy5nPZpG0qhXHrjh3Z+XaZ 5bER2gQ3SMH SOaugaNVATERSgJ/SCXo1n4814Nz5mpQUj79lHf2B8/gAjCZ5J</latexit> <latexit sha1_base64="4fQo4Jpdl8nw3AmI+ XyHse7Qqg=">A CE3icbVDLSsNAFJ3UV42vqks3wSKISEmKoMuiXbgRK9oHNKFMprfp0Mk zEyE voJLtWPcSdu/QC/xY2TNgt PTBwO feO/ceP2ZUKtv+MgpLy ura8V1c2Nza3untLvXklEiCDRJxCLR8bE Rjk0FVUMOrEAHPoM2v7oKvPbjyAkjfiDGsfghTjgdEAJVlq6r/eqvVLZrthTWIvEyUkZ5Wj0St9uPyJ CFwRhqXsOnasvBQLRQmDiekmEmJMRjiArqYchyC9dLrqxDrS t8aREI/rqyp+rsjxaGU49DXlSFWQznvZeJ/XjdRgwsvpTxOFHAy+2iQMEtFVna31acCiGJj T ARVO9qkSEWmCidjunWQd8i4EbPvY1BYBWJk9TFIg pn+jbAvc0Y6ZOy5nPZpG0qhXHrjh3Z+XaZ 5bER2gQ3SMH SOaugaNVATERSgJ/SCXo1n4814Nz5mpQUj79lHf2B8/gAjCZ5J</latexit> <latexit sha1_base64="4fQo4Jpdl8nw3AmI+ XyHse7Qqg=">A CE3icbVDLSsNAFJ3UV42vqks3wSKISEmKoMuiXbgRK9oHNKFMprfp0Mk zEyE voJLtWPcSdu/QC/xY2TNgt PTBwO feO/ceP2ZUKtv+MgpLy ura8V1c2Nza3untLvXklEiCDRJxCLR8bE Rjk0FVUMOrEAHPoM2v7oKvPbjyAkjfiDGsfghTjgdEAJVlq6r/eqvVLZrthTWIvEyUkZ5Wj0St9uPyJ CFwRhqXsOnasvBQLRQmDiekmEmJMRjiArqYchyC9dLrqxDrS t8aREI/rqyp+rsjxaGU49DXlSFWQznvZeJ/XjdRgwsvpTxOFHAy+2iQMEtFVna31acCiGJj T ARVO9qkSEWmCidjunWQd8i4EbPvY1BYBWJk9TFIg pn+jbAvc0Y6ZOy5nPZpG0qhXHrjh3Z+XaZ 5bER2gQ3SMH SOaugaNVATERSgJ/SCXo1n4814Nz5mpQUj79lHf2B8/gAjCZ5J</latexit> <latexit sha1_base64="4fQo4Jpdl8nw3AmI+ XyHse7Qqg=">A CE3icbVDLSsNAFJ3UV42vqks3wSKISEmKoMuiXbgRK9oHNKFMprfp0Mk zEyE voJLtWPcSdu/QC/xY2TNgt PTBwO feO/ceP2ZUKtv+MgpLy ura8V1c2Nza3untLvXklEiCDRJxCLR8bE Rjk0FVUMOrEAHPoM2v7oKvPbjyAkjfiDGsfghTjgdEAJVlq6r/eqvVLZrthTWIvEyUkZ5Wj0St9uPyJ CFwRhqXsOnasvBQLRQmDiekmEmJMRjiArqYchyC9dLrqxDrS t8aREI/rqyp+rsjxaGU49DXlSFWQznvZeJ/XjdRgwsvpTxOFHAy+2iQMEtFVna31acCiGJj T ARVO9qkSEWmCidjunWQd8i4EbPvY1BYBWJk9TFIg pn+jbAvc0Y6ZOy5nPZpG0qhXHrjh3Z+XaZ 5bER2gQ3SMH SOaugaNVATERSgJ/SCXo1n4814Nz5mpQUj79lHf2B8/gAjCZ5J</latexit>

Neural net

  • Func. Map

T : S1 → S2

<latexit sha1_base64="jsg/YwaKstOMgPlQ2pkgRWaCgy8=">A CKHicbVBNTxsxEPVSPtMCoUhcerGIKlUIRbsICcQJAYdeKlJBACkbRbPOZGPhtVf2LCha8md6LPyY3qpc+Rtc EIOLfAkS8/vzdgzL8mVdBSGo2Dmw+zc/MLiUuXjp+WV1era5wtnCiuwKYwy 9ioBh0pqbJIkhVe5RcgShZfJ9fHYv7xB6 TR5zTIsZ1BqmVPCiAvdaob5wf8rBPx2Mq0T2CtufX3nU61FtbDCfhbEk1JjU3R6FSf4q4R Ya hALnWlGYU7sES1IoHFbiwmEO4hpSbHmqIUPXLifzD/lXr3R5z1h/NPGJ+m9HCZlzgyzxlRlQ3732xuJ7Xqug3n67lDov CLV4+ahXKE6Gj8PgXWlRkBp4AsJKPysXfbAgyEdWiU/Q72Lxh3/3NEcLZOxWGYN M6mHfrc03h6zik8rep3NW3KxU4/CevRzt3Z4NM1tkX1hm+wbi9geO2TfWYM1mWB37Be7Zw/B7+BP8DcYvZTOBNOedfYfgsdnmG lvQ= </latexit> <latexit sha1_base64="jsg/YwaKstOMgPlQ2pkgRWaCgy8=">A CKHicbVBNTxsxEPVSPtMCoUhcerGIKlUIRbsICcQJAYdeKlJBACkbRbPOZGPhtVf2LCha8md6LPyY3qpc+Rtc EIOLfAkS8/vzdgzL8mVdBSGo2Dmw+zc/MLiUuXjp+WV1era5wtnCiuwKYwy 9ioBh0pqbJIkhVe5RcgShZfJ9fHYv7xB6 TR5zTIsZ1BqmVPCiAvdaob5wf8rBPx2Mq0T2CtufX3nU61FtbDCfhbEk1JjU3R6FSf4q4R Ya hALnWlGYU7sES1IoHFbiwmEO4hpSbHmqIUPXLifzD/lXr3R5z1h/NPGJ+m9HCZlzgyzxlRlQ3732xuJ7Xqug3n67lDov CLV4+ahXKE6Gj8PgXWlRkBp4AsJKPysXfbAgyEdWiU/Q72Lxh3/3NEcLZOxWGYN M6mHfrc03h6zik8rep3NW3KxU4/CevRzt3Z4NM1tkX1hm+wbi9geO2TfWYM1mWB37Be7Zw/B7+BP8DcYvZTOBNOedfYfgsdnmG lvQ= </latexit> <latexit sha1_base64="jsg/YwaKstOMgPlQ2pkgRWaCgy8=">A CKHicbVBNTxsxEPVSPtMCoUhcerGIKlUIRbsICcQJAYdeKlJBACkbRbPOZGPhtVf2LCha8md6LPyY3qpc+Rtc EIOLfAkS8/vzdgzL8mVdBSGo2Dmw+zc/MLiUuXjp+WV1era5wtnCiuwKYwy 9ioBh0pqbJIkhVe5RcgShZfJ9fHYv7xB6 TR5zTIsZ1BqmVPCiAvdaob5wf8rBPx2Mq0T2CtufX3nU61FtbDCfhbEk1JjU3R6FSf4q4R Ya hALnWlGYU7sES1IoHFbiwmEO4hpSbHmqIUPXLifzD/lXr3R5z1h/NPGJ+m9HCZlzgyzxlRlQ3732xuJ7Xqug3n67lDov CLV4+ahXKE6Gj8PgXWlRkBp4AsJKPysXfbAgyEdWiU/Q72Lxh3/3NEcLZOxWGYN M6mHfrc03h6zik8rep3NW3KxU4/CevRzt3Z4NM1tkX1hm+wbi9geO2TfWYM1mWB37Be7Zw/B7+BP8DcYvZTOBNOedfYfgsdnmG lvQ= </latexit> <latexit sha1_base64="jsg/YwaKstOMgPlQ2pkgRWaCgy8=">A CKHicbVBNTxsxEPVSPtMCoUhcerGIKlUIRbsICcQJAYdeKlJBACkbRbPOZGPhtVf2LCha8md6LPyY3qpc+Rtc EIOLfAkS8/vzdgzL8mVdBSGo2Dmw+zc/MLiUuXjp+WV1era5wtnCiuwKYwy 9ioBh0pqbJIkhVe5RcgShZfJ9fHYv7xB6 TR5zTIsZ1BqmVPCiAvdaob5wf8rBPx2Mq0T2CtufX3nU61FtbDCfhbEk1JjU3R6FSf4q4R Ya hALnWlGYU7sES1IoHFbiwmEO4hpSbHmqIUPXLifzD/lXr3R5z1h/NPGJ+m9HCZlzgyzxlRlQ3732xuJ7Xqug3n67lDov CLV4+ahXKE6Gj8PgXWlRkBp4AsJKPysXfbAgyEdWiU/Q72Lxh3/3NEcLZOxWGYN M6mHfrc03h6zik8rep3NW3KxU4/CevRzt3Z4NM1tkX1hm+wbi9geO2TfWYM1mWB37Be7Zw/B7+BP8DcYvZTOBNOedfYfgsdnmG lvQ= </latexit>

FMNet

Neural net

slide-57
SLIDE 57

FMNet

57

Training loss: Learning approach to computing descriptors.

  • O. Litany, T. Remez, E. Rodolà, A. Bronstein, M. Bronstein: Deep functional maps:

Structured prediction for dense shape correspondence. In Proc. ICCV (2017).

slide-58
SLIDE 58

Our Goals

58

1. Avoid using ground truth correspondences

  • Replace supervised loss with unsupervised one

2. Avoid using geodesic distances

  • Perform all computations in the spectral domain

Main question: how to measure the quality of a map?

Note: related concurrent paper by Halimi et al. Unsupervised learning

  • f dense shape correspondence. In CVPR, 2019
slide-59
SLIDE 59

59

Replace supervised loss with unsupervised one

Unsupervised Deep Learning for Structured Shape Matching, J.-M. Rouffosse, A. Sharma, M. O., ICCV 2019

Our approach

FMNet FMNet

Ereg(C12, C21)

D1

D2

C12 =

C21 =

arg min

C

∥CAT (D1) − AT (D2)∥2, arg min

C

∥CAT (D2) − AT (D1)∥2

T(D1)

T(D2)

slide-60
SLIDE 60

Our approach

Bijectivity Area-preservation Functional map close to pointwise one.

All penalties are in the reduced basis. 50x faster than FMNet

Near-isometry

E1(C12, C21) = ∥C12C21 − Id∥2 E1(C12, C21) = ∥C21C12 − Id∥2

<latexit sha1_base64="WgkVSqF5wLBrfFa0rp5a+7DQM5k=">A CpHicfVFta9RAEN6kvtT4dtaPg we6in1SI5C/SIUW0FBsaLXFm6vYbOZ5JZuNmF3Ixwxf8J/5z8p+MXNXUBtxYFlH56Z 2bn2aS wtgw/OH5G1euXru+eSO4e v2nbuDe1tHpqw1xykvZalP EmZQCoVTK6zEk0ojKxKJx8nZfpc/ oraiFJ9scsK5wXLlcgEZ9ZR8eA7lZhZ2gBNMBeqYVqzZdvI9k0cjfbjJpq02+DuSdQ+g1dAv8Ga7Dl4Ae9Sx5 OACiF4H+qr wX/1ZRVGk/FKgW+cKO48EwHIergMsg6sGQ9HEYD37StOR1gcpy YyZRWFl56 rFVxiG9DaYMX4 Gctx5qBiBZp5s/KuhceOS ErtTvKwor9U9Gw phlkbjKgtmFuZjryH/lZrXNXs4boaraouLrQVktwZbQfQSkQiO3cukA41q4twJfM 24d 8V0AN0u2j84Pp+rFAzW+rnDWU6L4Rq3W453e5Q4NyKLnpzGRxNxlE4j 7tDPde975tkgfkERmRiOySPfKWHJIp4eTce+g9 9Ub+E/+9/9mfrkt9r9fcJ3+Ff/oLK7PJCw= </latexit> <latexit sha1_base64="WgkVSqF5wLBrfFa0rp5a+7DQM5k=">A CpHicfVFta9RAEN6kvtT4dtaPg we6in1SI5C/SIUW0FBsaLXFm6vYbOZ5JZuNmF3Ixwxf8J/5z8p+MXNXUBtxYFlH56Z 2bn2aS wtgw/OH5G1euXru+eSO4e v2nbuDe1tHpqw1xykvZalP EmZQCoVTK6zEk0ojKxKJx8nZfpc/ oraiFJ9scsK5wXLlcgEZ9ZR8eA7lZhZ2gBNMBeqYVqzZdvI9k0cjfbjJpq02+DuSdQ+g1dAv8Ga7Dl4Ae9Sx5 OACiF4H+qr wX/1ZRVGk/FKgW+cKO48EwHIergMsg6sGQ9HEYD37StOR1gcpy YyZRWFl56 rFVxiG9DaYMX4 Gctx5qBiBZp5s/KuhceOS ErtTvKwor9U9Gw phlkbjKgtmFuZjryH/lZrXNXs4boaraouLrQVktwZbQfQSkQiO3cukA41q4twJfM 24d 8V0AN0u2j84Pp+rFAzW+rnDWU6L4Rq3W453e5Q4NyKLnpzGRxNxlE4j 7tDPde975tkgfkERmRiOySPfKWHJIp4eTce+g9 9Ub+E/+9/9mfrkt9r9fcJ3+Ff/oLK7PJCw= </latexit> <latexit sha1_base64="WgkVSqF5wLBrfFa0rp5a+7DQM5k=">A CpHicfVFta9RAEN6kvtT4dtaPg we6in1SI5C/SIUW0FBsaLXFm6vYbOZ5JZuNmF3Ixwxf8J/5z8p+MXNXUBtxYFlH56Z 2bn2aS wtgw/OH5G1euXru+eSO4e v2nbuDe1tHpqw1xykvZalP EmZQCoVTK6zEk0ojKxKJx8nZfpc/ oraiFJ9scsK5wXLlcgEZ9ZR8eA7lZhZ2gBNMBeqYVqzZdvI9k0cjfbjJpq02+DuSdQ+g1dAv8Ga7Dl4Ae9Sx5 OACiF4H+qr wX/1ZRVGk/FKgW+cKO48EwHIergMsg6sGQ9HEYD37StOR1gcpy YyZRWFl56 rFVxiG9DaYMX4 Gctx5qBiBZp5s/KuhceOS ErtTvKwor9U9Gw phlkbjKgtmFuZjryH/lZrXNXs4boaraouLrQVktwZbQfQSkQiO3cukA41q4twJfM 24d 8V0AN0u2j84Pp+rFAzW+rnDWU6L4Rq3W453e5Q4NyKLnpzGRxNxlE4j 7tDPde975tkgfkERmRiOySPfKWHJIp4eTce+g9 9Ub+E/+9/9mfrkt9r9fcJ3+Ff/oLK7PJCw= </latexit> <latexit sha1_base64="WgkVSqF5wLBrfFa0rp5a+7DQM5k=">A CpHicfVFta9RAEN6kvtT4dtaPg we6in1SI5C/SIUW0FBsaLXFm6vYbOZ5JZuNmF3Ixwxf8J/5z8p+MXNXUBtxYFlH56Z 2bn2aS wtgw/OH5G1euXru+eSO4e v2nbuDe1tHpqw1xykvZalP EmZQCoVTK6zEk0ojKxKJx8nZfpc/ oraiFJ9scsK5wXLlcgEZ9ZR8eA7lZhZ2gBNMBeqYVqzZdvI9k0cjfbjJpq02+DuSdQ+g1dAv8Ga7Dl4Ae9Sx5 OACiF4H+qr wX/1ZRVGk/FKgW+cKO48EwHIergMsg6sGQ9HEYD37StOR1gcpy YyZRWFl56 rFVxiG9DaYMX4 Gctx5qBiBZp5s/KuhceOS ErtTvKwor9U9Gw phlkbjKgtmFuZjryH/lZrXNXs4boaraouLrQVktwZbQfQSkQiO3cukA41q4twJfM 24d 8V0AN0u2j84Pp+rFAzW+rnDWU6L4Rq3W453e5Q4NyKLnpzGRxNxlE4j 7tDPde975tkgfkERmRiOySPfKWHJIp4eTce+g9 9Ub+E/+9/9mfrkt9r9fcJ3+Ff/oLK7PJCw= </latexit>

E4(C) =

  • i

∥CXfi − YgiC∥2

<latexit sha1_base64="hKJXk5HeOm7C6uQBrmfQ19aPFMk=">A CQXicbZB T9swGIadskHXDejYcRdrFVI3QUmqSuyCVNFN4jIBEqWdmhI57pfUqu1EtoNUhf6M/RqObD9iP2E3tGsvOKWHUfZJlh6/rz/b3xumnGnjur+d0tqLl+sb5VeV1282t7arb3cudZIp Cl2a8ET1Q6KBMwldw yHfq AiJBDL5x0Cr93DUqzRF6YaQpDQWLJIkaJsVJQPfgatOqdj/gI+zoTQc5m2L/BHdwP8i wm3 8Pcj gjrWuGoG1ZrbcBeFn4O3hBpa1l QnfujhGYCpKGcaD3w3NQMc6IMox mFT/TkBI6ITEMLEoiQA/zxWAzvGuVEY4SZ c0eKH+25ET
  • fVUhPakIGasV71C/J83yEz0eZgzmWYGJH18KMo4NgkuUsIjpoAaPrVAqGL2r5iOiSLU2Cwr/hewsyj4Zu89TUERk6hPuU9ULJic2dlif6+gik3LW83mOVw2G57b8M5btfbxMrcyeo8+oDry0CFqoxN0hrqIoh/oFv1Ev5w7549z7/x9PFpylj3v0JNy5g/uta5G</la
texit> <latexit sha1_base64="hKJXk5HeOm7C6uQBrmfQ19aPFMk=">A CQXicbZB T9swGIadskHXDejYcRdrFVI3QUmqSuyCVNFN4jIBEqWdmhI57pfUqu1EtoNUhf6M/RqObD9iP2E3tGsvOKWHUfZJlh6/rz/b3xumnGnjur+d0tqLl+sb5VeV1282t7arb3cudZIp Cl2a8ET1Q6KBMwldw yHfq AiJBDL5x0Cr93DUqzRF6YaQpDQWLJIkaJsVJQPfgatOqdj/gI+zoTQc5m2L/BHdwP8i wm3 8Pcj gjrWuGoG1ZrbcBeFn4O3hBpa1l QnfujhGYCpKGcaD3w3NQMc6IMox mFT/TkBI6ITEMLEoiQA/zxWAzvGuVEY4SZ c0eKH+25ET
  • fVUhPakIGasV71C/J83yEz0eZgzmWYGJH18KMo4NgkuUsIjpoAaPrVAqGL2r5iOiSLU2Cwr/hewsyj4Zu89TUERk6hPuU9ULJic2dlif6+gik3LW83mOVw2G57b8M5btfbxMrcyeo8+oDry0CFqoxN0hrqIoh/oFv1Ev5w7549z7/x9PFpylj3v0JNy5g/uta5G</la
texit> <latexit sha1_base64="hKJXk5HeOm7C6uQBrmfQ19aPFMk=">A CQXicbZB T9swGIadskHXDejYcRdrFVI3QUmqSuyCVNFN4jIBEqWdmhI57pfUqu1EtoNUhf6M/RqObD9iP2E3tGsvOKWHUfZJlh6/rz/b3xumnGnjur+d0tqLl+sb5VeV1282t7arb3cudZIp Cl2a8ET1Q6KBMwldw yHfq AiJBDL5x0Cr93DUqzRF6YaQpDQWLJIkaJsVJQPfgatOqdj/gI+zoTQc5m2L/BHdwP8i wm3 8Pcj gjrWuGoG1ZrbcBeFn4O3hBpa1l QnfujhGYCpKGcaD3w3NQMc6IMox mFT/TkBI6ITEMLEoiQA/zxWAzvGuVEY4SZ c0eKH+25ET
  • fVUhPakIGasV71C/J83yEz0eZgzmWYGJH18KMo4NgkuUsIjpoAaPrVAqGL2r5iOiSLU2Cwr/hewsyj4Zu89TUERk6hPuU9ULJic2dlif6+gik3LW83mOVw2G57b8M5btfbxMrcyeo8+oDry0CFqoxN0hrqIoh/oFv1Ev5w7549z7/x9PFpylj3v0JNy5g/uta5G</la
texit> <latexit sha1_base64="hKJXk5HeOm7C6uQBrmfQ19aPFMk=">A CQXicbZB T9swGIadskHXDejYcRdrFVI3QUmqSuyCVNFN4jIBEqWdmhI57pfUqu1EtoNUhf6M/RqObD9iP2E3tGsvOKWHUfZJlh6/rz/b3xumnGnjur+d0tqLl+sb5VeV1282t7arb3cudZIp Cl2a8ET1Q6KBMwldw yHfq AiJBDL5x0Cr93DUqzRF6YaQpDQWLJIkaJsVJQPfgatOqdj/gI+zoTQc5m2L/BHdwP8i wm3 8Pcj gjrWuGoG1ZrbcBeFn4O3hBpa1l QnfujhGYCpKGcaD3w3NQMc6IMox mFT/TkBI6ITEMLEoiQA/zxWAzvGuVEY4SZ c0eKH+25ET
  • fVUhPakIGasV71C/J83yEz0eZgzmWYGJH18KMo4NgkuUsIjpoAaPrVAqGL2r5iOiSLU2Cwr/hewsyj4Zu89TUERk6hPuU9ULJic2dlif6+gik3LW83mOVw2G57b8M5btfbxMrcyeo8+oDry0CFqoxN0hrqIoh/oFv1Ev5w7549z7/x9PFpylj3v0JNy5g/uta5G</la
texit>

E3(C) = ∥Λ2C − CΛ1∥2

<latexit sha1_base64="CzTN9/B8D Ce+3FqU2comsGDPO0=">A CQnicbVBdaxNBFJ2NtcZU2 gf RkaCqlo2E0FfRFC04IPlkYwH5CNy93JT pkdnaZmRXCun/DX+Nj9T/0L/hW+mgfOvmg2KQXLhzOPf OnBMmgmvjupdO4dHG480nxaelrWfPt3fKL152dJwq hm0Wi1j1QtAouMS24UZgL1EIUSiwG06as3n3OyrNY/nVTBMcRDCWfMQZGEsFZfckOKw2D+hH6v+g/me7OIQgq+e0Sd/avmO83Aq+1SkNyhW35s6LrgNvCSpkWa2g/M8fxiyNUBomQOu+5yZmkIEynAnMS36qMQE2gTH2LZQ oR5kc2c53bfMkI5iZVsaOmf/38g 0noa hVYZgTnXq7MZ+dCsn5rRh0HGZ IalGzx0CgV1MR0FhMdcoXMiKkFwBS3f6XsHBQwY8Ms+cdovSg8tXfPElRgYvU680GNIy5z623sv5mhk 3LW81mHXTqNc+teV/eVRpHy9yK5BXZI1XikfekQT6RFmkTRn6SX+Q3+eNcOH+dK+d6IS04y51dcq+cm1vE0a2N</latexi t> <latexit sha1_base64="CzTN9/B8D Ce+3FqU2comsGDPO0=">A CQnicbVBdaxNBFJ2NtcZU2 gf RkaCqlo2E0FfRFC04IPlkYwH5CNy93JT pkdnaZmRXCun/DX+Nj9T/0L/hW+mgfOvmg2KQXLhzOPf OnBMmgmvjupdO4dHG480nxaelrWfPt3fKL152dJwq hm0Wi1j1QtAouMS24UZgL1EIUSiwG06as3n3OyrNY/nVTBMcRDCWfMQZGEsFZfckOKw2D+hH6v+g/me7OIQgq+e0Sd/avmO83Aq+1SkNyhW35s6LrgNvCSpkWa2g/M8fxiyNUBomQOu+5yZmkIEynAnMS36qMQE2gTH2LZQ oR5kc2c53bfMkI5iZVsaOmf/38g 0noa hVYZgTnXq7MZ+dCsn5rRh0HGZ IalGzx0CgV1MR0FhMdcoXMiKkFwBS3f6XsHBQwY8Ms+cdovSg8tXfPElRgYvU680GNIy5z623sv5mhk 3LW81mHXTqNc+teV/eVRpHy9yK5BXZI1XikfekQT6RFmkTRn6SX+Q3+eNcOH+dK+d6IS04y51dcq+cm1vE0a2N</latexi t> <latexit sha1_base64="CzTN9/B8D Ce+3FqU2comsGDPO0=">A CQnicbVBdaxNBFJ2NtcZU2 gf RkaCqlo2E0FfRFC04IPlkYwH5CNy93JT pkdnaZmRXCun/DX+Nj9T/0L/hW+mgfOvmg2KQXLhzOPf OnBMmgmvjupdO4dHG480nxaelrWfPt3fKL152dJwq hm0Wi1j1QtAouMS24UZgL1EIUSiwG06as3n3OyrNY/nVTBMcRDCWfMQZGEsFZfckOKw2D+hH6v+g/me7OIQgq+e0Sd/avmO83Aq+1SkNyhW35s6LrgNvCSpkWa2g/M8fxiyNUBomQOu+5yZmkIEynAnMS36qMQE2gTH2LZQ oR5kc2c53bfMkI5iZVsaOmf/38g 0noa hVYZgTnXq7MZ+dCsn5rRh0HGZ IalGzx0CgV1MR0FhMdcoXMiKkFwBS3f6XsHBQwY8Ms+cdovSg8tXfPElRgYvU680GNIy5z623sv5mhk 3LW81mHXTqNc+teV/eVRpHy9yK5BXZI1XikfekQT6RFmkTRn6SX+Q3+eNcOH+dK+d6IS04y51dcq+cm1vE0a2N</latexi t> <latexit sha1_base64="CzTN9/B8D Ce+3FqU2comsGDPO0=">A CQnicbVBdaxNBFJ2NtcZU2 gf RkaCqlo2E0FfRFC04IPlkYwH5CNy93JT pkdnaZmRXCun/DX+Nj9T/0L/hW+mgfOvmg2KQXLhzOPf OnBMmgmvjupdO4dHG480nxaelrWfPt3fKL152dJwq hm0Wi1j1QtAouMS24UZgL1EIUSiwG06as3n3OyrNY/nVTBMcRDCWfMQZGEsFZfckOKw2D+hH6v+g/me7OIQgq+e0Sd/avmO83Aq+1SkNyhW35s6LrgNvCSpkWa2g/M8fxiyNUBomQOu+5yZmkIEynAnMS36qMQE2gTH2LZQ oR5kc2c53bfMkI5iZVsaOmf/38g 0noa hVYZgTnXq7MZ+dCsn5rRh0HGZ IalGzx0CgV1MR0FhMdcoXMiKkFwBS3f6XsHBQwY8Ms+cdovSg8tXfPElRgYvU680GNIy5z623sv5mhk 3LW81mHXTqNc+teV/eVRpHy9yK5BXZI1XikfekQT6RFmkTRn6SX+Q3+eNcOH+dK+d6IS04y51dcq+cm1vE0a2N</latexi t>

E2(C) = ∥CT C − Id∥2

<latexit sha1_base64="SY1mjogmSReAQo1WPj+rMid7iWw=">A CL3icbVBJSwMxGM3Urdat6lGEYBGqaJkpgl6EYhX0IFawC3Tak m/1mAmMyQZoYw9+Ws8qj9GvIhXf4MX0+Xg9iDw8r4lec8LOVPatl+txMTk1PRMcjY1N7+wuJReXqmoIJIUyjTg ax5RAFn Asqa Q61UALxPQ5V76Y4qFdvQSoWiCvdC6Hhk65gHUaJNlIrvX7SymeLW/gQu3e42LzCRbyLz9rm1sxj3Epn7Jw9BP5LnDHJoDFKrfSn2w5o5IPQlBOl6o4d6kZMpGaUQz/lRgpCQm9IF+qGCuKDasRDG328aZQ27gTSHKHxUP0+ERNfqZ7vmU6f6Gv1uzYQ/6vVI905 aMRMhJEGQUcPdSKOdYAHmeA2k0A17xlCqGTmr5heE0moNsml3GMwXiScm70XIUi A7kdu0R2fSb6xlvX3RmwlEnL+Z3NX1LJ5xw751zuZQpH49ySaA1toCxy0D4qoFNUQmVE0T16QE/o2Xq0Xqw363 UmrDGM6voB6yPL0p4pcI=</latexit> <latexit sha1_base64="SY1mjogmSReAQo1WPj+rMid7iWw=">A CL3icbVBJSwMxGM3Urdat6lGEYBGqaJkpgl6EYhX0IFawC3Tak m/1mAmMyQZoYw9+Ws8qj9GvIhXf4MX0+Xg9iDw8r4lec8LOVPatl+txMTk1PRMcjY1N7+wuJReXqmoIJIUyjTg ax5RAFn Asqa Q61UALxPQ5V76Y4qFdvQSoWiCvdC6Hhk65gHUaJNlIrvX7SymeLW/gQu3e42LzCRbyLz9rm1sxj3Epn7Jw9BP5LnDHJoDFKrfSn2w5o5IPQlBOl6o4d6kZMpGaUQz/lRgpCQm9IF+qGCuKDasRDG328aZQ27gTSHKHxUP0+ERNfqZ7vmU6f6Gv1uzYQ/6vVI905 aMRMhJEGQUcPdSKOdYAHmeA2k0A17xlCqGTmr5heE0moNsml3GMwXiScm70XIUi A7kdu0R2fSb6xlvX3RmwlEnL+Z3NX1LJ5xw751zuZQpH49ySaA1toCxy0D4qoFNUQmVE0T16QE/o2Xq0Xqw363 UmrDGM6voB6yPL0p4pcI=</latexit> <latexit sha1_base64="SY1mjogmSReAQo1WPj+rMid7iWw=">A CL3icbVBJSwMxGM3Urdat6lGEYBGqaJkpgl6EYhX0IFawC3Tak m/1mAmMyQZoYw9+Ws8qj9GvIhXf4MX0+Xg9iDw8r4lec8LOVPatl+txMTk1PRMcjY1N7+wuJReXqmoIJIUyjTg ax5RAFn Asqa Q61UALxPQ5V76Y4qFdvQSoWiCvdC6Hhk65gHUaJNlIrvX7SymeLW/gQu3e42LzCRbyLz9rm1sxj3Epn7Jw9BP5LnDHJoDFKrfSn2w5o5IPQlBOl6o4d6kZMpGaUQz/lRgpCQm9IF+qGCuKDasRDG328aZQ27gTSHKHxUP0+ERNfqZ7vmU6f6Gv1uzYQ/6vVI905 aMRMhJEGQUcPdSKOdYAHmeA2k0A17xlCqGTmr5heE0moNsml3GMwXiScm70XIUi A7kdu0R2fSb6xlvX3RmwlEnL+Z3NX1LJ5xw751zuZQpH49ySaA1toCxy0D4qoFNUQmVE0T16QE/o2Xq0Xqw363 UmrDGM6voB6yPL0p4pcI=</latexit> <latexit sha1_base64="SY1mjogmSReAQo1WPj+rMid7iWw=">A CL3icbVBJSwMxGM3Urdat6lGEYBGqaJkpgl6EYhX0IFawC3Tak m/1mAmMyQZoYw9+Ws8qj9GvIhXf4MX0+Xg9iDw8r4lec8LOVPatl+txMTk1PRMcjY1N7+wuJReXqmoIJIUyjTg ax5RAFn Asqa Q61UALxPQ5V76Y4qFdvQSoWiCvdC6Hhk65gHUaJNlIrvX7SymeLW/gQu3e42LzCRbyLz9rm1sxj3Epn7Jw9BP5LnDHJoDFKrfSn2w5o5IPQlBOl6o4d6kZMpGaUQz/lRgpCQm9IF+qGCuKDasRDG328aZQ27gTSHKHxUP0+ERNfqZ7vmU6f6Gv1uzYQ/6vVI905 aMRMhJEGQUcPdSKOdYAHmeA2k0A17xlCqGTmr5heE0moNsml3GMwXiScm70XIUi A7kdu0R2fSb6xlvX3RmwlEnL+Z3NX1LJ5xw751zuZQpH49ySaA1toCxy0D4qoFNUQmVE0T16QE/o2Xq0Xqw363 UmrDGM6voB6yPL0p4pcI=</latexit>

Replace supervised loss with unsupervised one

slide-61
SLIDE 61

Datasets

61

Datasets released as part of: Continuous and Orientation-preserving Correspondences via Functional Maps, J. Ren, A. Poulenard, P. Wonka, M. O, SIGGRAPH Asia 2018

slide-62
SLIDE 62

62

Comparison to unsupervised methods

Results

Unsupervised Deep Learning for Structured Shape Matching, J.-M. Rouffosse,

  • A. Sharma, M. O., ICCV 2019

State-of-the art among unsupervised methods.

slide-63
SLIDE 63

63

State-of-the art among unsupervised methods.

Results

Unsupervised Deep Learning for Structured Shape Matching, J.-M. Rouffosse,

  • A. Sharma, M. O., ICCV 2019
slide-64
SLIDE 64

64

Comparison to supervised methods Comparable results even to supervised methods

Results

Unsupervised Deep Learning for Structured Shape Matching, J.-M. Rouffosse,

  • A. Sharma, M. O., ICCV 2019
slide-65
SLIDE 65

65

Comparable results even to supervised methods

Results

Unsupervised Deep Learning for Structured Shape Matching, J.-M. Rouffosse,

  • A. Sharma, M. O., ICCV 2019
slide-66
SLIDE 66

66

Results

Original vs. learned descriptors.

Unsupervised Deep Learning for Structured Shape Matching, J.-M. Rouffosse,

  • A. Sharma, M. O., ICCV 2019
slide-67
SLIDE 67

67

  • 1. How can we build up a functional map progressively?
  • 2. Given a small functional map, can we use it to transfer

high frequency functions?

  • 3. Simplify and speed-up functional map refinement?

Several related questions

slide-68
SLIDE 68

ZoomOut

68

A two-lines-of-code algorithm:

1) Given a functional map C1 of size k x k convert it to a p2p map T. 2) Convert T to C2 of size (k+1) x (k+1)

Repeat for progressively larger k

ZoomOut: Spectral Upsampling for Efficient Shape Correspondence, S. Melzi, J. Ren, A. Sharma, E. Rodolà, P. Wonka, M. O., SIGGRAPH Asia 2019

slide-69
SLIDE 69

69

Upsampling vs. computing directly:

ZoomOut: Spectral Upsampling for Efficient Shape Correspondence, S. Melzi, J. Ren, A. Sharma, E. Rodolà, P. Wonka, M. O., SIGGRAPH Asia 2019

ZoomOut

slide-70
SLIDE 70

70

Extreme case, from 2x2 to 100x100

Dataset provided by the Natural History Museum in Paris.

ZoomOut – Results

slide-71
SLIDE 71

ZoomOut – Results

71

From 5x5 to 50x50

ZoomOut: Spectral Upsampling for Efficient Shape Correspondence, S. Melzi, J. Ren, A. Sharma, E. Rodolà, P. Wonka, M. O., SIGGRAPH Asia 2019

slide-72
SLIDE 72

72

From 20x20 to 120x120

ZoomOut – Results

slide-73
SLIDE 73

73

Ours is 50-300x faster than state-of-the-art with higher accuracy Evaluated on:

  • Intrinsic symmetry detection
  • Complete matching
  • Partial matching
  • Function transfer

… Compared against 14 baselines

ZoomOut – Results

slide-74
SLIDE 74

74

Consider the optimization problem:

Theorem:

if and only if the point-to-point map is an isometry. ZoomOut can be derived as a iterative method for solving this

  • ptimization problem.

: functional map arising from some pointwise map. : leading principal submatrix of .

E(C) = 0

<latexit sha1_base64="K2LMUGTHjne5kmz0VeSMUsRlsWY=">A CI3icbVDLSsNAFJ34tr5SXboZLIK lEQE3QjiA9yICtYWmlIm05t26GQSZiZKCfkUl+rHuBM3LvwSN07aL R6YOBw7mPOPX7MmdKO82FNTE5Nz8zOzZcWFpeWV+zy6p2KEkmhRiMeyYZPFHAm
  • KaZ5tCIJZDQ51D3+6d5vX4PUrFI3OpBDK2QdAULGCXaSG27fL7lhUT3/CA9zb xEXbadsWpOkPgv8QtSAUVuG7bX14nok IQlNOlGq6TqxbKZGaUQ5ZyUsUxIT2SReahgoSgmqlQ+sZ3jRKBweRNE9oPFR/TqQkVGoQ+qYzt6nGa7n4X62Z6OCwlTIRJxoEHX0UJBzr
COc54A6TQDUfGEKoZMYrpj0iCdUmrZJ3BuYWCZdm71UMkuhI7qQekd2Qiczc1vV2c1Yyabnj2fwld3tV16m6N/uV45Mitzm0j bQFnLRATpGF+ga1RBFD+gRPaMX68l6td6s91HrhFXMrKFfsD6/AVEio30=</latexit> <latexit sha1_base64="K2LMUGTHjne5kmz0VeSMUsRlsWY=">A CI3icbVDLSsNAFJ34tr5SXboZLIK lEQE3QjiA9yICtYWmlIm05t26GQSZiZKCfkUl+rHuBM3LvwSN07aL R6YOBw7mPOPX7MmdKO82FNTE5Nz8zOzZcWFpeWV+zy6p2KEkmhRiMeyYZPFHAm
  • KaZ5tCIJZDQ51D3+6d5vX4PUrFI3OpBDK2QdAULGCXaSG27fL7lhUT3/CA9zb xEXbadsWpOkPgv8QtSAUVuG7bX14nok IQlNOlGq6TqxbKZGaUQ5ZyUsUxIT2SReahgoSgmqlQ+sZ3jRKBweRNE9oPFR/TqQkVGoQ+qYzt6nGa7n4X62Z6OCwlTIRJxoEHX0UJBzr
COc54A6TQDUfGEKoZMYrpj0iCdUmrZJ3BuYWCZdm71UMkuhI7qQekd2Qiczc1vV2c1Yyabnj2fwld3tV16m6N/uV45Mitzm0j bQFnLRATpGF+ga1RBFD+gRPaMX68l6td6s91HrhFXMrKFfsD6/AVEio30=</latexit> <latexit sha1_base64="K2LMUGTHjne5kmz0VeSMUsRlsWY=">A CI3icbVDLSsNAFJ34tr5SXboZLIK lEQE3QjiA9yICtYWmlIm05t26GQSZiZKCfkUl+rHuBM3LvwSN07aL R6YOBw7mPOPX7MmdKO82FNTE5Nz8zOzZcWFpeWV+zy6p2KEkmhRiMeyYZPFHAm
  • KaZ5tCIJZDQ51D3+6d5vX4PUrFI3OpBDK2QdAULGCXaSG27fL7lhUT3/CA9zb xEXbadsWpOkPgv8QtSAUVuG7bX14nok IQlNOlGq6TqxbKZGaUQ5ZyUsUxIT2SReahgoSgmqlQ+sZ3jRKBweRNE9oPFR/TqQkVGoQ+qYzt6nGa7n4X62Z6OCwlTIRJxoEHX0UJBzr
COc54A6TQDUfGEKoZMYrpj0iCdUmrZJ3BuYWCZdm71UMkuhI7qQekd2Qiczc1vV2c1Yyabnj2fwld3tV16m6N/uV45Mitzm0j bQFnLRATpGF+ga1RBFD+gRPaMX68l6td6s91HrhFXMrKFfsD6/AVEio30=</latexit> <latexit sha1_base64="K2LMUGTHjne5kmz0VeSMUsRlsWY=">A CI3icbVDLSsNAFJ34tr5SXboZLIK lEQE3QjiA9yICtYWmlIm05t26GQSZiZKCfkUl+rHuBM3LvwSN07aL R6YOBw7mPOPX7MmdKO82FNTE5Nz8zOzZcWFpeWV+zy6p2KEkmhRiMeyYZPFHAm
  • KaZ5tCIJZDQ51D3+6d5vX4PUrFI3OpBDK2QdAULGCXaSG27fL7lhUT3/CA9zb xEXbadsWpOkPgv8QtSAUVuG7bX14nok IQlNOlGq6TqxbKZGaUQ5ZyUsUxIT2SReahgoSgmqlQ+sZ3jRKBweRNE9oPFR/TqQkVGoQ+qYzt6nGa7n4X62Z6OCwlTIRJxoEHX0UJBzr
COc54A6TQDUfGEKoZMYrpj0iCdUmrZJ3BuYWCZdm71UMkuhI7qQekd2Qiczc1vV2c1Yyabnj2fwld3tV16m6N/uV45Mitzm0j bQFnLRATpGF+ga1RBFD+gRPaMX68l6td6s91HrhFXMrKFfsD6/AVEio30=</latexit>

ZoomOut – Rationale

ZoomOut: Spectral Upsampling for Efficient Shape Correspondence, S. Melzi, J. Ren, A. Sharma, E. Rodolà, P. Wonka, M. O., SIGGRAPH Asia 2019

slide-75
SLIDE 75

75

In some cases also works for non-isometric shapes

ZoomOut – Non-isometric

slide-76
SLIDE 76

Other Extensions

Maps in Collections Promoting Pointwise Maps

Informative Descriptor Preservation via Commutativity for Shape Matching, Nogneng, O., Eurographics 2017 Functional map networks for analyzing and exploring large shape collections Huang, Wang, Guibas, SIGGRAPH 2014

Manifold Optimization

MADMM: A generic algorithm for non-smooth optimization on manifolds. Kovnatsky, Glashoff, M. Bronstein, ECCV, 2016.

Huang et al., SIGGRAPH 2014

slide-77
SLIDE 77

Consistency via Latent Space Optimization

Application to Co-segmentation:

Image Co-Segmentation via Consistent Functional Maps Wang, Huang, Guibas, CVPR 2013

slide-78
SLIDE 78

Other Extensions

Measuring Differences between shapes Tangent Vector Field processing Maps Between Partial shapes

Partial Functional Correspondence, Rodolà, Cosmo, A. Bronstein, Torsello, Cremers, CGF 2017 Map-Based Exploration of Intrinsic Shape Differences and Variability Rustamov, Ovsjanikov, Azencot, Ben-Chen, Chazal, Guibas, SIGGRAPH 2014 An Operator Approach to Tangent Vector Field Processing Azencot, Ben-Chen, Chazal, Ovsjanikov, SGP, 2013.

Azencot et al., SGP 2014

slide-79
SLIDE 79

Some Open Problems

What is the optimal choice of basis? How to guarantee a continuous pointwise map? What are better deformation models? Shape interpolation without converting to p2p?

slide-80
SLIDE 80

Conclusions

Functional maps provide an efficient way to encode “generalized” mappings. Can be computed in practice with simple (least squares) optimization. Many different constraints can be incorporated: pointwise maps, consistency in collections, etc. Recent work incorporating learning of descriptors.

slide-81
SLIDE 81

Questions?

Acknowledgements:

  • A. Poulenard, M.-J. Rakotosaona, Y. Ponty, J.-M. Rouffosse, A. Sharma, S.

Melzi, E. Rodolà, J. Ren, P. Wonka …. Work supported by KAUST OSR Award No. CRG-2017-3426, a gift from Nvidia and the ERC Starting Grant StG-2017-758800 (EXPROTEA)

Thank You