12.07.2015 Views

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 ...

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!