Diffshape-Younes.Pdf

Total Page:16

File Type:pdf, Size:1020Kb

Diffshape-Younes.Pdf Deformation analysis for shape and image processing Laurent Younes Contents General notation 5 Chapter 1. Elements of Hilbert space theory 7 1. Definition and notation 7 2. Examples 8 3. Orthogonal spaces and projection 9 4. Orthonormal sequences 10 5. The Riesz representation theorem 11 6. Embeddings and the duality paradox 12 7. Weak convergence in a Hilbert space 15 8. Spline interpolation 16 9. Building V and its kernel 22 Chapter 2. Ordinary differential equations and groups of diffeomorphisms 35 1. Introduction 35 2. A class of ordinary differential equations 37 3. Groups of diffeomorphisms 44 Chapter 3. Matching deformable objects 49 1. General principles 49 2. Matching terms and greedy descent procedures 49 3. Regularized matching 70 Chapter 4. Groups, manifolds and invariant metrics 95 1. Introduction 95 2. Differential manifolds 95 3. Submanifolds 98 4. Lie Groups 99 5. Group action 102 6. Riemannian manifolds 103 7. Gradient descent on Riemannian manifolds 110 Bibliography 113 3 General notation The differential of a function ϕ :Ω → Rd, Ω being an open subset of Rd, is denoted Dϕ : x 7→ Dϕ(x). It takes values in the space of linear operators from Rd to Rd so that Dϕ(x) can be identified to d × k matrix (containing the partial derivatives of the coordinates of ϕ). d Idk is the identity matrix in R . ∞ CK (Ω) is the set of infinitely differentiable functions on Ω, an open subset of Rd. 5 CHAPTER 1 Elements of Hilbert space theory 1. Definition and notation A set H is a (real) Hilbert space if i) H is a vector space over R (i.e., addition and scalar multiplication are defined on H). 0 0 0 ii) H has an inner product denoted (h, h ) 7→ hh , h iH , for h, h ∈ H. This inner product is a bilinear form, which is positive definite and the associ- p ated norm is kh]kH = hh , hiH . iii) H is a complete space for the topology associated to the norm. Here are some explanations of iii). Converging sequences for the norm topology are sequences hn for which there exists h ∈ H such that kh − hnkH → 0. Property iii) means that if a sequence (hn, n ≥ 0) in h is a Cauchy sequence, i.e., it collapses in the sense that, for every positive ε there exists n0 > 0 such that khn − hn0 kH ≤ ε for n ≥ n0, then it necessarily has a limit: there exists h ∈ H such that khn − hk → 0 when n tends to infinity. If condition ii) is weakened to the fact that k.kH is a norm (not necessarily induced by an inner product), one says that H is a Banach space. If H satisfies i) and ii), it is called a pre-Hilbert space. On pre-Hilbert spaces, the Schwartz inequality holds: Proposition 1 (Schwartz inequality). If H is pre-Hilbert, and h, h0 ∈ H, then 0 0 hh , h iH ≤ khkH kh kH The first consequence of this property is Proposition 2. The inner product on H is continuous for the norm topology. Proof. The inner product is a function H × H → R. Letting ϕ(h, h0) = 0 hh , h iH , we have, by Schwartz inequality, and introducing a sequence (hn) which converges to h 0 0 0 |ϕ(h, h ) − ϕ(hn, h )| ≤ kh − hnkH kh kH → 0 which proves the continuity with respect to the first coordinate, and also with respect to the second coordinate by symmetry. ut Working with a complete normed vector space is essential when dealing with infinite sums of elements: if (hn) is a sequence in H, and if khn+1 + ··· + hn+kk can be made arbitrarily small for large n and any k > 0, then completeness implies P∞ that the series n=0 hn has a limit in h. In particular, absolutely converging series in H converge: X X (1) khnkH < ∞ ⇒ hn converges. n≥0 n≥0 7 8 1. ELEMENTS OF HILBERT SPACE THEORY We shall add a fourth condition to our definition of a Hilbert space iv) H is separable for the norm topology: there exists a countable subset S in H such that, for any h ∈ H and any ε > 0, there exists h0 ∈ S such that kh − h0k ≤ ε. In the following, Hilbert spaces will always be a separable Hilbert spaces (this is not always the case in the literature). A Hilbert space isometry between two Hilbert spaces H and H0 is a an invertible 0 0 0 0 linear map F : H → H such that, for all h, h ∈ H, hF h , F h iH0 = hh , h iH . A sub-Hilbert space of H is a subspace H0 of H (i.e., a non-empty subset closed 0 by linear combination) which is closed for the norm topology: if (hn) is a sequence 0 0 0 in H and khn − hkH → 0 for some h ∈ H, then h ∈ H . Topological closedness implies that H0 is itself a Hilbert space, since Cauchy sequences in H0 also are Cauchy sequences in H, hence converge in H, hence in H0 since H0 is closed. The next proposition shows that every finite dimensional subspace of a Hilbert space is closed, and hence a Hilbert subspace. Proposition 3. If K is a finite dimensional subspace of H, then K is closed in H Proof. let e1, . , ep be a basis of K. Let (hn) be a sequence in K which Pp converges to some h ∈ H: then hn may be written hn = k=1 aknek and for all l = 1, . , p: p X hhn , eliH = hek , eliH akn k=1 p p Let an be the vector in R with coordinates (akn, k = 1, . , p) and un ∈ R with coordinates (hhn , ekiH , k = 1, . , p). Let also S be the matrix with coefficients skl = hek , eliH , so that the previous system may be written: un = San. The matrix S is invertible: if b belongs to the null space of S, a quick computation shows that p 2 X T biei = b Sb = 0 i=1 H which is only possible when b = 0, since (e1, . , en) is a basis. We therefore have −1 an = S un. Since un converges to u with coordinates (hh , ekiH , k = 1, . , p) (by −1 continuity of the inner product), we obtain the fact that an converges to a = S u. Pp Pp But this implies that k=1 aknek → k=1 akek and since the limit is unique, we Pp have h = k=1 akek ∈ K. ut 2. Examples n 0 2.1. Finite dimensional Euclidean spaces. H = R with hh , h i n = Pn 0 R i=1 hihi is the standard example of a finite dimensional Hilbert space. 2.2. `2 space of real sequences. Let H be the set of real sequences h = P∞ 2 (h1, h2,...) such that i=1 hi < ∞. Then H is a Hilbert space. 3. ORTHOGONAL SPACES AND PROJECTION 9 2.3. L2 space of functions. Let k and d be two integers. Let Ω be an open subset of Rd. We define L2(Ω, Rd) as the set of all square integrable functions h :Ω → Rd, with inner-product Z 0 0 hh , h iL2 = h(x)h (x)dx Ω Integrals are taken with respect to the Lebesgue measure on Ω, and two functions which coincide everywhere except on a set of null Lebesgue measure are consid- ered as equal. The fact that L2(Ω, Rd) is a Hilbert space is a standard result in integration theory. 3. Orthogonal spaces and projection Let O be a subset of H. The orthogonal space to O is defined by ⊥ O = {h : h ∈ H, ∀o ∈ O, hh , oiH = 0} Theorem 1. O⊥ is a sub-Hilbert space of H. Proof. Stability by linear combination is obvious, and closedness is a conse- quence of the continuity of the inner product. ut When K is a sub-Hilbert space of H (a closed subspace) and h ∈ H, one can define the variational problem: 0 0 (PK (h)) : find k ∈ K such that kk − hkH = min {kh − k kH : k ∈ K} The following theorem is fundamental Theorem 2. If K is a closed subspace of H and h ∈ H, (PK (h)) has a unique solution, characterized by the property k ∈ K and h − k ∈ K⊥ Definition 1. The solution of problem (PK (h)) in the previous theorem is called the orthogonal projection of h on K and will be denoted πK (h). Proposition 4. πK : H → K is a linear, continuous function and πK⊥ = id − πK . 0 0 Proof. Let d = inf {kh − k kH : k ∈ K}. The proof relies on the following construction: if k, k0 ∈ K, a direct computation shows that 0 2 k + k 0 2 2 0 2 h − + kk − k kH /4 = kh − kkH + kh − k kH /2 2 H 0 2 0 k+k and the fact that (k + k )/2 ∈ K implies h − 2 ≥ d so that H 0 2 2 0 2 2 kk − k kH ≤ kh − kkH + kh − k kH /2 − d . Now, from the definition of the infimum, one can find a sequence kn in K such that 2 2 −n kkn − hkH ≤ d + 2 for each n. The previous inequality implies that 2 −n −m kkn − kmkH ≤ 2 + 2 /2 which implies that kn is a Cauchy sequence, and therefore converges to a limit k 0 which belongs to K and kk − hkH = d. If the minimum is attained for another k in K, we have, by the same inequality 0 2 2 2 2 kk − k kH ≤ d + d /2 − d = 0 .
Recommended publications
  • Computer Science at Kent
    Computer Science at Kent . Constructive Potential Theory: Foundations and Applications. Marius-Constantin BUJORIANU Computing Laboratory, University of Kent at Canterbury, UK . Manuela-Luminita BUJORIANU Department of Engineering, University of Cambridge, UK . ––––––––––––––––––––––––––––– Technical Report No: 6-02 Date June 2002 . Copyright c 2001 University of Kent at Canterbury Published by the Computing Laboratory, University of Kent, Canterbury, Kent CT2 7NF, UK. Abstract Stochastic analysis is now an important common part of computing and mathematics. Its applications are impressive, ranging from stochastic concurrent and hybrid systems to finances and biomedicine. In this work we investigate the logical and algebraic foundations of stochastic analysis and possible applications to computing. We focus more concretely on functional analysis theoretic core of stochastic analysis called potential theory. Classical potential theory originates in Gauss and Poincare’s work on partial differential equations. Modern potential theory now study stochastic processes with their adjacent theory, higher order differential operators and their combination like stochastic differential equations. In this work we consider only the axiomatic branches of modern potential theory, like Dirichlet forms and harmonic spaces. Due to the inherently constructive character of axiomatic potential theory, classical logic has no enough ability to offer a proper logical foundation. In this paper we propose the weak commutative linear logics as a logical framework for reasoning about the processes described by potential theory. The logical approach is complemented by an algebraic one. We construct an algebraic theory with models in stochastic analysis, and based on this, and a process algebra in the sense of computer science. Applications of these in area of hybrid systems, concurrency theory and biomedicine are investigated.
    [Show full text]
  • Towards a Simple Mathematical Theory of Citation Distributions
    Katchanov SpringerPlus (2015) 4:677 DOI 10.1186/s40064-015-1467-8 RESEARCH Open Access Towards a simple mathematical theory of citation distributions Yurij L. Katchanov* *Correspondence: [email protected] Abstract National Research University The paper is written with the assumption that the purpose of a mathematical theory Higher School of Economics (HSE), 20 Myasnitskaya Ulitsa, of citation is to explain bibliometric regularities at the level of mathematical formalism. 101000 Moscow, Russia A mathematical formalism is proposed for the appearance of power law distributions in social citation systems. The principal contributions of this paper are an axiomatic characterization of citation distributions in terms of the Ekeland variational principle and a mathematical exploration of the power law nature of citation distributions. Apart from its inherent value in providing a better understanding of the mathematical under- pinnings of bibliometric models, such an approach can be used to derive a citation distribution from first principles. Keywords: Bibliometrics, Citation distributions, Power law distribution, Wakeby distribution Mathematics Subject Classification: 91D30, 91D99 Background Scholars have been investigating their own citation practice for too long. Bibliometrics already forced considerable changes in citation practice Michels and Schmoch (2014). Because the overwhelming majority of bibliometric studies focus on the citation sta- tistics of scientific papers (see, e.g., Adler et al. 2009; Albarrán and Ruiz-Castillo 2011; De Battisti and Salini 2013; Nicolaisen 2007; Yang and Han 2015), special attention is devoted to citation distributions (see, inter alia, Radicchi and Castellano 2012; Ruiz-Cas- tillo 2012; Sangwal 2014; Thelwall and Wilson 2014; Vieira and Gomes 2010). However, the fundamentals of the citation distribution (or CD for convenience) are far from being well established and the universal law of CD is still unknown (we do not go into details, and refer the reader to Bornmann and Daniel 2009; Eom and Fortunato 2011; Peterson et al.
    [Show full text]
  • Smoothing - Strichartz Estimates for Dispersive Equations Perturbed by a First Order Differential Operator
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Electronic Thesis and Dissertation Archive - Università di Pisa UNIVERSITÀ DI PISA SCUOLA DI DOTTORATO "GALILEO GALILEI" DIPARTIMENTO DI MATEMATICA "L. TONELLI" DOTTORATO DI RICERCA IN MATEMATICA - ANNO 2003 CODICE SETTORE SCIENTIFICO DISCIPLINARE MAT/05 Ph. D. Thesis Smoothing - Strichartz Estimates for Dispersive Equations Perturbed by a First Order Differential Operator Mirko Tarulli Di Giallonardo Advisor Candidate Prof.VladimirGeorgiev MirkoTarulliDiGiallonardo Contents Contents ii Introduction 1 1. Background of Mathematical Physics 1 1.1. General Setting 1 1.2. Problem Setting 2 2. A-Priori Estimates 3 2.1. Resolvent Estimates 3 2.2. Dispersive Estimates 4 2.3. Smoothing Estimates 5 3. Plan of the Thesis 6 4. Acknowledgements 6 Chapter 1. Functional Analysis Background 9 1. Operator and Spectral Theory 9 2. Linear operators in Banach spaces 9 3. Symmetric strictly monotone operators on Hilbert space 14 4. Spectral Families 19 5. Integration with Respect to a Spectral Family 21 6. Some facts about holomorphic functions 24 7. Spectral Theorem for Self-Adjoint Operators 27 8. Some Interpolation Results 35 8.1. Interpolation for sequences with values in Banach spaces 37 9. Hardy-Littlewood-Sobolev Inequality 37 10. TT ∗ Method 38 11. Paley-Littlewood Partition of Unity 39 Chapter 2. Resolvent Estimates 45 1. The Limiting Absorption Principle 45 1.1. The Free Case 45 1.2. The Perturbed Case 49 2. The Resonances 52 3. The Resolvent Estimates 52 2 4. Resolvent estimate for R± λ 56 ∇ 0 Chapter 3. Applications to a Class of Dispersive Equations 59 iii Contents 1.
    [Show full text]
  • Download Article (PDF)
    DEMONSTRATE MATHEMATICA Vol. XXXIV No 3 2001 Joanna Napiôrkowska REGULARITY OF SOLUTIONS OF NONLINEAR PARABOLIC EQUATIONS 1. Introduction One of the main topics considered by mathematicians dealing with evo- lution partial differential equations is the problem of regularity of solu- tions. This problem was investigated by many authors, in particular by A. Friedman [FR], O. A. Ladyzhenskaya, V. A. Solonnikov and N. N. Ural- ceva [L-S-U], S. Smale [SM], C. Bardos [BA]. In the context of semi-linear evolution equations of parabolic type, the regularity was studied, e.g., by R. Temam [TE]. Consider a nonlinear evolution problem of the form ù + Au = F(u), t > 0, (1.1) u(0)I = u0. In this paper we study the smoothness of solutions for the Cauchy prob- lem of the type (1.1) with A being an abstract counterpart of an elliptic op- erator. Our main result is a necessary and sufficient condition of the higher order regularity of solutions for such nonlinear partial differential equation. In Section 2 we introduce some preliminaries. Section 3 is devoted to the regularity result of Temam [TE] in the case of linear equations. Section 4 contains the main result of the present paper. Finally, the applications and examples are described in Section 5. 2. Preliminaries Let if be a real Hilbert space with scalar product (•, •) and the norm || • ||. Let (A,D(A)) be a linear operator acting from D(A) into H, where D(A) C H is the domain of A. Assume that (A,D(A)) is a positive (i.e., (Au, u) > 0 for all u € D(A) \ {0}) and self-adjoint (see, e.g., [PA]) operator with D(A) dense in H.
    [Show full text]
  • Spectral Analysis of Elastic Waveguides
    BEYKENT ÜNİVERSİTESİ FEN VE MÜHENDİSLİK BİLİMLERİ DERGİSİ CİLT SAYI:13/1 www.dergipark.gov.tr SPECTRAL ANALYSIS OF ELASTIC WAVEGUIDES Mahir HASANSOY * ABSTRACT This paper deals with the spectral analysis of elastic waveguides shaped like an infinite cylinder with a bounded or unbounded cross section. Note that waveguides with bounded cross sections were studied in the framework of the operator and operator pencil theory. In this paper we extend this theory on elastic waveguides with a unbounded cross section. The problem reduces to the spectral theory of a one-parameter family of unbounded operators. The spectral structure, the asymptotics of the eigenvalues, comparisons between the solutions of different problems and the existence of special guided modes for these operators are the main questions that we study in this paper. We use an operator approach to solve these problems, and on the basis of this approach, we suggest alternative methods to solve spectral problems arising in the theory of both closed elastic waveguides and elastic waveguides with unbounded cross sections. Keywords: Eigenvalue, spectrum, operator, wave solution, elastic waveguide, perturbation. *Makale Gönderim Tarihi: 15.05.2020 ; Makale Kabul Tarihi : 30.05.2020 Makale Türü: Araştırma DOI: 10.20854/bujse.738083 *Beykent University, Faculty of Arts and Sciences, Department of Mathematics, Ayazağa Mahallesi, Hadım Koruyolu Cd. No:19, Sariyer Istanbul, Turkey ([email protected]) (ORCID ID: 0000-0002-9080-4242) AMS subject Classifications: Primary: 47A56, 47A10, 49R50; Secondary: 34L15. 43 BEYKENT ÜNİVERSİTESİ FEN VE MÜHENDİSLİK BİLİMLERİ DERGİSİ CİLT SAYI:13/1 www.dergipark.gov.tr ELASTİK DALGA KLAVUZLARININ SPEKTRAL ANALİZİ Mahir HASANSOY * ÖZ Bu makalede, sınırlı veya sınırsız kesite sahip sonsuz bir silindir biçiminde şekillendirilmiş elastik dalga kılavuzlarının spektral analizi ele alınmıştır.
    [Show full text]
  • Table of Contents
    A Comprehensive Course in Geometry Table of Contents Chapter 1 - Introduction to Geometry Chapter 2 - History of Geometry Chapter 3 - Compass and Straightedge Constructions Chapter 4 - Projective Geometry and Geometric Topology Chapter 5 - Affine Geometry and Analytic Geometry Chapter 6 - Conformal Geometry Chapter 7 - Contact Geometry and Descriptive Geometry Chapter 8 - Differential Geometry and Distance Geometry Chapter 9 - Elliptic Geometry and Euclidean Geometry Chapter 10 WT- Finite Geometry and Hyperbolic Geometry _____________________ WORLD TECHNOLOGIES _____________________ A Course in Boolean Algebra Table of Contents Chapter 1 - Introduction to Boolean Algebra Chapter 2 - Boolean Algebras Formally Defined Chapter 3 - Negation and Minimal Negation Operator Chapter 4 - Sheffer Stroke and Zhegalkin Polynomial Chapter 5 - Interior Algebra and Two-Element Boolean Algebra Chapter 6 - Heyting Algebra and Boolean Prime Ideal Theorem Chapter 7 - Canonical Form (Boolean algebra) Chapter 8 -WT Boolean Algebra (Logic) and Boolean Algebra (Structure) _____________________ WORLD TECHNOLOGIES _____________________ A Course in Riemannian Geometry Table of Contents Chapter 1 - Introduction to Riemannian Geometry Chapter 2 - Riemannian Manifold and Levi-Civita Connection Chapter 3 - Geodesic and Symmetric Space Chapter 4 - Curvature of Riemannian Manifolds and Isometry Chapter 5 - Laplace–Beltrami Operator Chapter 6 - Gauss's Lemma Chapter 7 - Types of Curvature in Riemannian Geometry Chapter 8 - Gauss–Codazzi Equations Chapter 9 -WT Formulas
    [Show full text]
  • $ W $-Sobolev Spaces: Theory, Homogenization and Applications
    W -SOBOLEV SPACES: THEORY, HOMOGENIZATION AND APPLICATIONS ALEXANDRE B. SIMAS AND FABIO´ J. VALENTIM Abstract. Fix strictly increasing right continuous functions with left limits d d Wi : R → R, i = 1,...,d, and let W (x) = Pi=1 Wi(xi) for x ∈ R . We construct the W -Sobolev spaces, which consist of functions f having weak ∇ generalized gradients W f = (∂W1 f,...,∂Wd f). Several properties, that are analogous to classical results on Sobolev spaces, are obtained. W -generalized elliptic and parabolic equations are also established, along with results on ex- istence and uniqueness of weak solutions of such equations. Homogenization results of suitable random operators are investigated. Finally, as an applica- tion of all the theory developed, we prove a hydrodynamic limit for gradient processes with conductances (induced by W ) in random environments. 1. Introduction The space of functions that admit differentiation in a weak sense has been widely studied in the mathematical literature. The usage of such spaces provides a wide application to the theory of partial differential equations (PDE), and to many other areas of pure and applied mathematics. These spaces have become associated with the name of the late Russian mathematician S. L. Sobolev, although their origins predate his major contributions to their development in the late 1930s. In theory of PDEs, the idea of Sobolev space allows one to introduce the notion of weak solutions whose existence, uniqueness, regularities, and well-posedness are based on tools of functional analysis. In classical theory of PDEs, two important classes of equations are: elliptic and parabolic PDEs. They are second-order PDEs, with some constraints (coerciveness) in the higher-order terms.
    [Show full text]
  • ¢¡¤£¦ ¥¨§ ©§ ©§ ¡¤ © ! "$# % %£ %'&) (10 #¡¤20 ©"30£ ! "©¡¤ "¡¤ 465¤
    ¢¡¤£¦¥¨§ © § © §¡¤ © ! £ ¡¤ © £ © ¡¤ ¡¤ "$#¤%'&'(*)¨+,+.-*)%0/21 3 %4/5)6#7#78¤9*+287:;)<&'(*=6+28¤%?>@ ¢H. I ¡7 © ¨A B 0© ¡¤C D ¡6C ¡7E ¡7¡'F G © £ KJ B L<MON¢P QSR TP U¦V!WYX7Z P U ¢¡¤£¦¥¢§©¨ §¢ ¥ £¦¨¤¤¤ ¡¤!¥ "# $£¦¥ £¦¨¨% &('*),+.-0/21436587:9;74'=<>3@?¤A3CBED Acknowledgments I wish to express my deep gratitude to my thesis advisor Professor V.Matsaev for his valuable guidance, essential support and permanent attention. It is also a pleasure to thank to Professor N.Kopachevsky who introduced me to the wonderful world of Mathematics and support me during many years. I am also indebted to Professor I.Gohberg for many suggestions and encouragement. I thank the staff of Tel Aviv University for the permanent help in technical problems. Last but not least I want to thank my wife Natalie, my parents and my children for patience and warm moral support. i Abstract The Floquet approach is a main tool of the theory of linear ordinary differential equations (o.d.e.) with periodic coefficients. Such equations arise in many physi- cal and technical applications. At the same time, a lot of theoretical and applied problems that appear in quantum mechanics, hydrodynamics, solid state physics, theory of wave conductors, parametric resonance theory, etc. lead to periodic partial differential equations. The case of ordinary periodical differential equations is in a sharp contrast with a partial one, because a space of solutions of partial differential equation is, generally, an infinite-dimensional space. Methods of the classical Flo- quet theory (i.e. method of monodromy operator and method of substitution) in general are not applicable and new tools and technique are necessary.
    [Show full text]