We can eventually estimate the difference |k 1 (τ) − k 2 (τ)| by us<strong>in</strong>g e.g. the<strong>in</strong>equality|A 2 B 2 C 2 − A 1 B 1 C 1 | ≤ |A 1 B 1 | |C 2 − C 1 | + |A 1 C 2 | |B 2 − B 1 | + |B 2 C 2 | |A 2 − A 1 |on the difference <strong>of</strong> the <strong>in</strong>tegr<strong>and</strong>s( 12 − l 2(τ) − l 2 (θ))h 2 (θ) |c ′L} {{ 2 } {{ } }2(θ)|{{ }}B 2 C 2A 2( 1−2 − l 1(τ) − l 1 (θ))L 1} {{ }A 1to obta<strong>in</strong> that the above term is less or equal than‖h 1 ‖ 0 ‖c ′ 2 − c ′ 1‖ 0 + ‖c ′ 2‖ 0 ‖h 2 − h 0 ‖ 0 + ‖h 2 ‖ 0 ‖c ′ 2‖ 0∣ ∣∣∣ l 2 (θ) − l 2 (τ)L 2h 1 (θ) |c ′} {{ } }1(θ)|{{ }B 1 C 1− l ∣1(θ) − l 1 (τ) ∣∣∣L 1s<strong>in</strong>ce |A i | ≤ 1/2. In turn (s<strong>in</strong>ce the formulas def<strong>in</strong><strong>in</strong>g k 1 (τ), k 2 (τ) the parametersare bound by τ − 2π ≤ θ ≤ τ) then|A 2 − A 1 | =l 2 (τ) − l 2 (θ)∣− l ∣ 1(τ) − l 1 (θ) ∣∣∣ ∫ τ∫ τ∣ θ=|c′ 2(x)| dxθ−|c′ 1(x)| dx ∣∣∣∣≤L 2L 1 ∣ L 2L 1≤ π L 1 + L 2‖c ′ 1 − c ′L 1 L2‖ 0 + π ‖c′ 1‖ 0 + ‖c ′ 2‖ 0|L 1 − L 2 | ≤2 L 1 L 2≤ 4π 2 ‖c′ 1‖ 0 + ‖c ′ 2‖ 0L 1 L 2‖c ′ 1 − c ′ 2‖ 0(by equations (ii) <strong>and</strong> (i) <strong>in</strong> lemma 10.25). Summariz<strong>in</strong>g(|k 1 (τ) − k 2 (τ)| ≤ 2π ‖h 1 ‖ 0 + ‖h 2 ‖ 0 ‖c ′ 2‖ 0 4π 2 ‖c′ 1‖ 0 + ‖c ′ )2‖ 0‖c ′ 1 − c ′L 1 L2‖ 0 +2+ 2π‖c ′ 2‖ 0 ‖h 2 − h 0 ‖ 0 .If we derive, k i ′(θ) = h i(θ)|c ′ i (θ)|, so|k ′ 2(θ)−k ′ 2(θ)| ≤ |h 2(θ)| + |h 1 (θ)|2Symmetriz<strong>in</strong>g we prove (10.30).|c ′ 2(θ)−c ′ 1(θ)| + |c′ 2(θ)| + |c ′ 1(θ)||h ′22(θ)−h ′ 1(θ)| .Lemma 10.27 Conversely, let c 1 , c 2 be two C 1 immersed <strong>curves</strong>, <strong>and</strong> h 1 , h 2be two differentiable fields; then (by us<strong>in</strong>g once aga<strong>in</strong> eqn. (i) from the lemma10.25)‖D c1 h 1 −D c2 h 2 ‖ 0 ≤ ‖c 1‖ 1 + ‖c 2 ‖ 12ε 2 ‖h 1 −h 2 ‖ 1 + ‖h 1‖ 1 + ‖h 2 ‖ 12ε 2 ‖c 1 −c 2 ‖ 1 (10.35)where ε = m<strong>in</strong>(<strong>in</strong>f S 1 |c ′ 1|, <strong>in</strong>f S 1 |c ′ 2|).76
10.6.2 Existence <strong>of</strong> flow for the centroid energy (2.9)Theorem 10.29 Let us fix v ∈ lR n , let E(c) = 1 2 |avg c(c) − v| 2 ; let the <strong>in</strong>itialcurve c 0 ∈ C 1 , then the ˜H 1 gradient descent flow <strong>of</strong> E has an unique solutionC = C(t, θ), for all t ∈ lR, <strong>and</strong> C(t, ·) ∈ C 1 .Results <strong>of</strong> a numerical simulation <strong>of</strong> the above gradient descent flow areshown <strong>in</strong> figure 12 on the follow<strong>in</strong>g page.Before prov<strong>in</strong>g the theorem, let us comment on an <strong>in</strong>terest<strong>in</strong>g property <strong>of</strong>the above gradient flow.Remark 10.30 The length <strong>of</strong> the <strong>curves</strong> len(C(t, ·)) is constant dur<strong>in</strong>g theevolution <strong>in</strong> t.Pro<strong>of</strong>. Indeed it is easy to prove that ∂ t len(C(t, ·)) = 0: the Gâteaux differential<strong>of</strong> the length <strong>of</strong> a curve c is∫D len(c; h) = (D s c · D s h) ds ;substitut<strong>in</strong>g the value <strong>of</strong> D s C(t, ·)) from eqn. (10.25),c∂ t len(C(t, ·)) = D(len(C))(∂ t C)∫1 〈=λL 2 D s C, ( D s C〈(avg c (C) − v), (C − avg c (C))〉 + α )〉 ds =〈∫〉 ∫1=λL 2 (avg c (C) − v), (C − avg c (C)) ds + (D s C · α) ds = 0 .We now present the pro<strong>of</strong> <strong>of</strong> the theorem 10.29.Pro<strong>of</strong>. We will first <strong>of</strong> all prove existence <strong>and</strong> uniqueness <strong>of</strong> the gradient flow forsmall time, us<strong>in</strong>g the Cauchy–Lipschitz theorem 25 <strong>in</strong> the space M <strong>of</strong> immersed<strong>curves</strong>, seen as an open subset <strong>of</strong> the C 1 Banach space (see Def<strong>in</strong>ition 10.24).To this end we will prove that c ↦→ ∇˜H1E(c) is a locally Lipschitz functionalfrom C 1 <strong>in</strong>to itself. Afterward, we will directly prove that the solution exists forall times.So, let us fix c 1 , c 2 ∈ C 1 ∩ M that are two <strong>curves</strong> near c 0 . Let us fix ε > 0 byε := <strong>in</strong>fθ |c′ 0(θ)|/2 .By “near” we mean that, <strong>in</strong> all <strong>of</strong> the follow<strong>in</strong>g, we will require that ‖c i −c 0 ‖ 1 < ε,for i = 1, 2. We will need the follow<strong>in</strong>g (easy to prove) facts. (In the follow<strong>in</strong>g,the <strong>in</strong>dex i will represent both 1, 2).• <strong>in</strong>f θ |c ′ i (θ)| ≥ ε, so all <strong>curves</strong> <strong>in</strong> the neighborhood are immersed.25 A.k.a. as the Picard–L<strong>in</strong>delöf theorem77
- Page 1 and 2:
Metrics of curves in shape optimiza
- Page 3 and 4:
shape analysis where we study a fam
- Page 5 and 6:
• F (c) = F (c ◦ φ) for all cu
- Page 7 and 8:
κ > 0HNNHNHκ < 0Figure 1: Example
- Page 9 and 10:
In the case of planar curves c 1 ,
- Page 11 and 12:
(a) (b) (c) (d) (e)Figure 3: Segmen
- Page 13 and 14:
where φ may be chosen to beφ(x) =
- Page 15 and 16:
2.4.1 Example: geometric heat flowW
- Page 17 and 18:
2.4.5 Centroid energyWe will now pr
- Page 20 and 21:
shapes. Unfortunately, H 0 does not
- Page 22 and 23:
• If the second request is waived
- Page 24 and 25:
• The Fréchet space of smooth fu
- Page 26 and 27: Since φ k are homeomorphisms, then
- Page 28 and 29: 3.6.1 Riemann metric, lengthDefinit
- Page 30 and 31: Theorem 3.30 Suppose that M is a sm
- Page 32 and 33: Example 3.38 Let M = C ∞ ([−1,
- Page 34 and 35: • sometimes S 1 will be identifie
- Page 36 and 37: The proof is by direct computation.
- Page 38 and 39: The term preshape space is sometime
- Page 40 and 41: 5 Representation/embedding/quotient
- Page 42 and 43: 6.1.1 Length induced by a distanceI
- Page 44 and 45: 0101200000000000011111111111110000
- Page 46 and 47: 6.2.4 Applications in computer visi
- Page 48 and 49: of a small ball from A. The motion
- Page 50 and 51: CutΩ(The arrows represent the dist
- Page 52 and 53: 7.3 L 1 metric and Plateau problemI
- Page 54 and 55: Definition 8.3 (Flat curves) Let Z
- Page 56 and 57: Proof. Fix α 0 ∈ S \ Z. Let T =
- Page 58 and 59: 8.2.2 RepresentationThe Stiefel man
- Page 60 and 61: 9.3 Conformal metricsYezzi and Menn
- Page 62 and 63: 10.1.1 Related worksA family of met
- Page 64 and 65: Definition 10.7 (Convolution) A arc
- Page 66 and 67: and̂∇˜HjE(0) = ̂∇ H 0E(0),̂
- Page 68 and 69: and we simply integrate twice! More
- Page 70 and 71: Corollary 10.18 In particular, the
- Page 72 and 73: 10.6 Existence of gradient flowsWe
- Page 74 and 75: • The length functional (from C 1
- Page 78 and 79: 7.552.5-10 -7.5 -5 -2.5 2.5 5 7.5 1
- Page 80 and 81: By using (i) and (ii) from lemma 10
- Page 82 and 83: 10.6.3 Existence of flow for geodes
- Page 84 and 85: 10.7.1 Robustness w.r.to local mini
- Page 86 and 87: We already presented all the calcul
- Page 88 and 89: 10.9 New regularization methodsTypi
- Page 90 and 91: can be solved for k (this is not so
- Page 92 and 93: 2. Now• at t = 1/2 it achieves th
- Page 94 and 95: Definition 11.10 • The orbit is O
- Page 96 and 97: y eqn. (11.3), where the terms RHS
- Page 98 and 99: So Prop. 11.19 guarantees that the
- Page 100 and 101: Theorem 11.23 Suppose that the Riem
- Page 102 and 103: 2.2.3 Examples of geometric energy
- Page 104 and 105: 9 Riemannian metrics of immersed cu
- Page 106 and 107: contour continuation, 11convolution
- Page 108 and 109: normal vector, 6objective function,
- Page 110 and 111: References[1] Luigi Ambrosio, Giuse
- Page 112 and 113: [30] Serge Lang. Fundamentals of di
- Page 114 and 115: [57] Ganesh Sundaramoorthi, Anthony