Deformable Sequence Reconstruction 1 EG 2012 Tutorial: Dynamic - - PowerPoint PPT Presentation

deformable sequence reconstruction
SMART_READER_LITE
LIVE PREVIEW

Deformable Sequence Reconstruction 1 EG 2012 Tutorial: Dynamic - - PowerPoint PPT Presentation

Deformable Sequence Reconstruction 1 EG 2012 Tutorial: Dynamic Geometry Processing Eurographics 2012, Cagliari, Italy Deformable Shape Matching Basic Principle Eurographics 2012, Cagliari, Italy Example Correspondences? Eurographics 2012,


slide-1
SLIDE 1

Deformable Sequence Reconstruction

EG 2012 Tutorial: Dynamic Geometry Processing

1

slide-2
SLIDE 2

Eurographics 2012, Cagliari, Italy

Deformable Shape Matching

Basic Principle

slide-3
SLIDE 3

Eurographics 2012, Cagliari, Italy

Example

Correspondences?

slide-4
SLIDE 4

Eurographics 2012, Cagliari, Italy

What are We Looking for?

Problem Statement: Given:

  • Two surfaces S1, S2 ⊆ ℝ3

We are looking for:

  • A reasonable deformation function f : S1 → ℝ3

that brings S1 close to S2

?

S1 S2

f

slide-5
SLIDE 5

Eurographics 2012, Cagliari, Italy

Example

?

correspondences? no shape match too much deformation

  • ptimum
slide-6
SLIDE 6

Eurographics 2012, Cagliari, Italy

This is a Trade-Off

Deformable Shape Matching is a Trade-Off:

  • We can match any two shapes

using a weird deformation field

  • We need to trade-off:

– Shape matching (close to data) – Regularity of the deformation field (reasonable match)

slide-7
SLIDE 7

Eurographics 2012, Cagliari, Italy

Components: Matching Distance: Deformation / rigidity:

Variational Model

slide-8
SLIDE 8

Eurographics 2012, Cagliari, Italy

Variational Model

Variational Problem:

  • Formulate as an energy minimization problem:
slide-9
SLIDE 9

Eurographics 2012, Cagliari, Italy

Part 1: Shape Matching

Data Term:

  • Objective Function:
  • Distance measures:

– Least-squares (L2) – Reweighted (robustness) – Hausdorf distance – Lp-distances, etc.

  • L2 measure is frequently used (models Gaussian noise)

– Reweighting/truncation for robustness

S2 f (S1)

slide-10
SLIDE 10

Eurographics 2012, Cagliari, Italy

Surface Approximation

Basic: Closest point matching

  • “Point-to-point” energy
  • Usually iterated: “Iterated Closest Points (ICP)”

– Establish nearest-neighbor correspondences – Minimize energy (with regularizer)

f (S1)

S2

slide-11
SLIDE 11

Eurographics 2012, Cagliari, Italy

Surface Approximation

Improvement: Linear approximation

  • “Point-to-plane” energy
  • Fit plane to k-nearest neighbors

S2

f (S1)

slide-12
SLIDE 12

Eurographics 2012, Cagliari, Italy

Robust Least-Squares

Robustness: Reweighting

  • Ignore Outliers

– Large distance – Connection to normal at large angle – Many matches to one point

slide-13
SLIDE 13

Eurographics 2012, Cagliari, Italy

Variational Model

Variational Problem:

  • Formulate as an energy minimization problem:
slide-14
SLIDE 14

Eurographics 2012, Cagliari, Italy

What is a “nice” deformation field?

  • Elastic deformation

– Volumetric elasticity – Thin shell model (more complex)

  • Intrinsic

– Isometric matching

  • Smooth deformations

– “Thin-plate-splines” (TPS) – Allowing strong deformations, but keep shape

Deformation Model

slide-15
SLIDE 15

Eurographics 2012, Cagliari, Italy

What is a “nice” deformation field?

  • Elastic deformation

– Volumetric elasticity – Thin shell model (more complex)

  • Intrinsic

– Isometric matching

  • Smooth deformations

– “Thin-plate-splines” (TPS) – Allowing strong deformations, but keep shape

Deformation Model

slide-16
SLIDE 16

Eurographics 2012, Cagliari, Italy

How to Detect Deformations?

S1 V1 f S2 ∇f f (V1)

slide-17
SLIDE 17

Eurographics 2012, Cagliari, Italy

How to Detect Deformations?

f

slide-18
SLIDE 18

Eurographics 2012, Cagliari, Italy

Elastic Volume Model

Extrinsic Volumetric “As-Rigid-As Possible”

  • Measure orthogonality
  • Integrate over deviation from orthogonality

S1 V1 f S2 ∇f

f (V1)

slide-19
SLIDE 19

Eurographics 2012, Cagliari, Italy

Deformable ICP

How to build a deformable ICP algorithm

  • Pick a surface distance measure
  • Pick an deformation model / regularizer

) ( ) ( ) (

) ( ) (

f E f E f E

r regularize match

+ =

slide-20
SLIDE 20

Eurographics 2012, Cagliari, Italy

Deformable ICP

Deformable ICP Algorithm

  • Select model: E(match), E(regularizer)
  • Initialize f (S1) with S1 (i.e., f = id)
  • (Non-linear) optimization:

– Newton, Gauss Newton – LBGFS (quick & effective)

slide-21
SLIDE 21

Eurographics 2012, Cagliari, Italy

Animation Reconstruction

Reconstructing Sequences of Deformable Shapes

slide-22
SLIDE 22

Eurographics 2012, Cagliari, Italy

“Factorization” Approach

t = 0 t = 1 t = 2

data urshape

S f f f

deformation

slide-23
SLIDE 23

Eurographics 2012, Cagliari, Italy

Hierarchical Merging S f(S) data f

slide-24
SLIDE 24

Eurographics 2012, Cagliari, Italy

Hierarchical Merging S f(S) data f

slide-25
SLIDE 25

Eurographics 2012, Cagliari, Italy

Initial Urshapes S f(S) data f

slide-26
SLIDE 26

Eurographics 2012, Cagliari, Italy

Initial Urshapes S f(S) data f

slide-27
SLIDE 27

Eurographics 2012, Cagliari, Italy

Alignment S f(S) data f

slide-28
SLIDE 28

Eurographics 2012, Cagliari, Italy

Alignment S f(S) data f

slide-29
SLIDE 29

Eurographics 2012, Cagliari, Italy

Hierarchical Alignment S f(S) data f

slide-30
SLIDE 30

Eurographics 2012, Cagliari, Italy

Hierarchical Alignment S f(S) data f

slide-31
SLIDE 31

Eurographics 2012, Cagliari, Italy

Global Optimization Energy Function

E(f, S) = Edata + Edeform + Esmooth

Components

  • Edata(f, S)

– data fitting

  • Edeform(f)

– elastic deformation, smooth trajectory

  • Esmooth(S)

– smooth surface

Final Optimization

  • Minimize over all frames
slide-32
SLIDE 32

Eurographics 2012, Cagliari, Italy

Deformation Field S f

  • Elastic energy
  • Smooth trajectories
slide-33
SLIDE 33

Eurographics 2012, Cagliari, Italy

Additional Terms More Regularization

  • Acceleration:

Eacc = ∫T ∫V|∂t

2 f|2

 Smooth trajectories

  • Velocity (weak):

Evel = ∫T ∫V|∂t f|2

 Damping

slide-34
SLIDE 34

Eurographics 2012, Cagliari, Italy

Surface Reconstruction

Data fitting

  • Smooth surface
  • Fitting to noisy data

34

S

fi

  • 1(Di)

S fi

  • 1(Di)

S

f

slide-35
SLIDE 35

Eurographics 2012, Cagliari, Italy

Results

(Joint work with: Bart Adams, Maksim Ovsjanikov, Alexander Berner, Martin Bokeloh, Philipp Jenke, Leonidas Guibas, Hans-Peter Seidel, Andreas Schilling)

slide-36
SLIDE 36

Eurographics 2012, Cagliari, Italy Eurographics 2012, Cagliari, Italy

slide-37
SLIDE 37

Eurographics 2012, Cagliari, Italy Eurographics 2012, Cagliari, Italy

slide-38
SLIDE 38

Eurographics 2012, Cagliari, Italy Eurographics 2012, Cagliari, Italy

slide-39
SLIDE 39

Eurographics 2012, Cagliari, Italy Eurographics 2012, Cagliari, Italy

slide-40
SLIDE 40
slide-41
SLIDE 41

Eurographics 2012, Cagliari, Italy

Edeform(f)

Elastic Deformation Energy S Di f

Regularization

  • Elastic energy
  • Smooth trajectories
slide-42
SLIDE 42

Eurographics 2012, Cagliari, Italy

geometry

Discretization Example Approach:

  • Full resolution geometry
  • Subsample deformation

deformation

slide-43
SLIDE 43

Eurographics 2012, Cagliari, Italy

Discretization “Subspace” Approach:

  • Sample volume
  • Place basis functions
  • Decouple from resolution of geometry

geometry deformation

slide-44
SLIDE 44

Eurographics 2012, Cagliari, Italy

Surface Reconstruction Data fitting

  • Smooth surface
  • Fitting to noisy data

S Di

Esmooth(S)

fi

  • 1(Di)

S fi

  • 1(Di)

S

f

slide-45
SLIDE 45

Eurographics 2012, Cagliari, Italy

Example

?

correspondences? no shape match too much deformation

  • ptimum