Metrics of curves in shape optimization and analysis - Andrea Carlo ...
Metrics of curves in shape optimization and analysis - Andrea Carlo ...
Metrics of curves in shape optimization and analysis - Andrea Carlo ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
10.6 Existence <strong>of</strong> gradient flowsWe recall this def<strong>in</strong>ition (that was already presented <strong>in</strong>formally <strong>in</strong> the <strong>in</strong>troduction).Def<strong>in</strong>ition 10.22 (Gradient descent flow) Given a differentiable energy E :M → lR, <strong>and</strong> a metric 〈, 〉 c , let ∇E(c) be the gradient. Let us fix moreoverc 0 : S 1 → lR n , c 0 ∈ M. The gradient descent flow <strong>of</strong> E is the solutionC = C(t, θ) <strong>of</strong> the <strong>in</strong>itial value P.D.E.{∂t C = −∇E(C)C(0, θ) = c 0 (θ)We present an example computation for the geodesic active contour modelon a radially symmetric “image”.Example 10.23 LetC(t, θ) = r(t) (cos θ, s<strong>in</strong> θ) ,φ(x) = d(|x|)with r, d : lR → lR + ; the energy takes the form∫E(C) = φ(C(t, s)) ds = 2π r(t) d(r(t))Then C is the ˜H 1 gradient descent iffCr ′ (t) = − 1 (d(r(t)) + r(t) d ′ (r(t)) ) .λ2πwhereas C is the H 0 gradient descent iffr ′ (t) = −2π ( d(r(t)) + r(t)d ′ (r(t)) ) ,where we use the def<strong>in</strong>ition (10.4) <strong>of</strong> H 0 .Note also thatPro<strong>of</strong>. Indeed,∂ t E(C) = r ′ (t) ( d(r(t)) + r(t)d ′ (r(t)) ) .ds = r(t) dθD s =1r(t) D θlen(C) = 2πr(t)φ(C) = d(r(t))∇φ(x) = x|x| d′ (|x|) for x ≠ 0∇φ(C) = d ′ (r(t))(cos θ, s<strong>in</strong> θ)avg c (∇φ(C)) = 072