Topological derivative

The topological derivative is, conceptually, a derivative of a shape functional with respect to infinitesimal changes in its topology, such as adding an infinitesimal hole or crack. When used in higher dimensions than one, the term topological gradient is also used to name the first-order term of the topological asymptotic expansion, dealing only with infinitesimal singular domain perturbations. It has applications in shape optimization, topology optimization, image processing and mechanical modeling.

Definition

Let ${\displaystyle \Omega }$ be an open bounded domain of ${\displaystyle \mathbb {R} ^{d}}$, with ${\displaystyle d\geq 2}$, which is subject to a nonsmooth perturbation confined in a small region ${\displaystyle \omega _{\varepsilon }({\tilde {x}})={\tilde {x}}+\varepsilon \omega }$ of size ${\displaystyle \varepsilon }$ with ${\displaystyle {\tilde {x}}}$ an arbitrary point of ${\displaystyle \Omega }$ and ${\displaystyle \omega }$ a fixed domain of ${\displaystyle \mathbb {R} ^{d}}$. Let ${\displaystyle \Psi }$ be a characteristic function associated to the unperturbed domain and ${\displaystyle \Psi _{\varepsilon }}$ be a characteristic function associated to the perforated domain ${\displaystyle \Omega _{\varepsilon }=\Omega \backslash {\overline {\omega _{\varepsilon }}}}$. A given shape functional ${\displaystyle \Phi (\Psi _{\varepsilon }({\tilde {x}}))}$ associated to the topologically perturbed domain, admits the following topological asymptotic expansion:

${\displaystyle \Phi (\Psi _{\varepsilon }({\tilde {x}}))=\Phi (\Psi )+f(\varepsilon )g({\tilde {x}})+o(f(\varepsilon ))}$

where ${\displaystyle \Phi (\Psi )}$ is the shape functional associated to the reference domain, ${\displaystyle f(\varepsilon )}$ is a positive first order correction function of ${\displaystyle \Phi (\Psi )}$ and ${\displaystyle o(f(\varepsilon ))}$ is the remainder. The function ${\displaystyle g({\tilde {x}})}$ is called the topological derivative of ${\displaystyle \Phi }$ at ${\displaystyle {\tilde {x}}}$.

Applications

Structural mechanics

The topological derivative can be applied to shape optimization problems in structural mechanics.[1] The topological derivative can be considered as the singular limit of the shape derivative. It is a generalization of this classical tool in shape optimization.[2] Shape optimization concerns itself with finding an optimal shape. That is, find ${\displaystyle \Omega }$ to minimize some scalar-valued objective function, ${\displaystyle J(\Omega )}$. The topological derivative technique can be coupled with level-set method.[3]

In 2005, the topological asymptotic expansion for the Laplace equation with respect to the insertion of a short crack inside a plane domain had been found. It allows to detect and locate cracks for a simple model problem: the steady-state heat equation with the heat flux imposed and the temperature measured on the boundary.[4] The topological derivative had been fully developed for a wide range of second-order differential operators and in 2011, it had been applied to Kirchhoff plate bending problem with a fourth-order operator.[5]

Image processing

In the field of image processing, in 2006, the topological derivative has been used to perform edge detection and image restoration. The impact of an insulating crack in the domain is studied. The topological sensitivity gives information on the image edges. The presented algorithm is non-iterative and thanks to the use of spectral methods has a short computing time.[6] Only ${\displaystyle O(Nlog(N))}$ operations are needed to detect edges, where ${\displaystyle N}$ is the number of pixels.[7] During the following years, other problems have been considered: classification, segmentation, inpainting and super-resolution.[7][8][9][10][11] This approach can be applied to gray-level or color images.[12] Until 2010, isotropic diffusion was used for image reconstructions. The topological gradient is also able to provide edge orientation and this information can be used to perform anisotropic diffusion.[13]

In 2012, a general framework is presented to reconstruct an image ${\displaystyle u\in L^{2}(\Omega )}$ given some noisy observations ${\displaystyle Lu+n}$ in a Hilbert space ${\displaystyle E}$ where ${\displaystyle \Omega }$ is the domain where the image ${\displaystyle u}$ is defined.[11] The observation space ${\displaystyle E}$ depends on the specific application as well as the linear observation operator ${\displaystyle L:L^{2}(\Omega )\rightarrow E}$. The norm on the space ${\displaystyle E}$ is ${\displaystyle \|.\|_{E}}$. The idea to recover the original image is to minimize the following functional for ${\displaystyle u\in H^{1}(\Omega )}$:

${\displaystyle \|C^{1/2}\nabla u\|_{L^{2}(\Omega )}^{2}+\|Lu-v\|_{E}^{2}}$

where ${\displaystyle C}$ is a positive definite tensor. The first term of the equation ensures that the recovered image ${\displaystyle u}$ is regular, and the second term measures the discrepancy with the data. In this general framework, different types of image reconstruction can be performed such as[11]

• image denoising with ${\displaystyle E=L^{2}(\Omega )}$ and ${\displaystyle Lu=u}$,
• image denoising and deblurring with ${\displaystyle E=L^{2}(\Omega )}$ and ${\displaystyle Lu=\phi \ast u}$ with ${\displaystyle \phi }$ a motion blur or Gaussian blur,
• image inpainting with ${\displaystyle E=L^{2}(\Omega \backslash \omega )}$ and ${\displaystyle Lu=u|_{\Omega \backslash \omega }}$, the subset ${\displaystyle \omega \subset \Omega }$ is the region where the image has to be recovered.

In this framework, the asymptotic expansion of the cost function ${\displaystyle J_{\Omega }(u_{\Omega })={\frac {1}{2}}\int _{\Omega }u_{\Omega }^{2}}$ in the case of a crack provides the same topological derivative ${\displaystyle g(x,n)=-\pi c(\nabla u_{0}.n)(\nabla p_{0}.n)-\pi (\nabla u_{0}.n)^{2}}$ where ${\displaystyle n}$ is the normal to the crack and ${\displaystyle c}$ a constant diffusion coefficient. The functions ${\displaystyle u_{0}}$ and ${\displaystyle p_{0}}$ are solutions of the following direct and adjoint problems.[11]

${\displaystyle -\nabla (c\nabla u_{0})+L^{*}Lu_{0}=L^{*}v}$ in ${\displaystyle \Omega }$ and ${\displaystyle \partial _{n}u_{0}=0}$ on ${\displaystyle \partial \Omega }$ ${\displaystyle -\nabla (c\nabla p_{0})+L^{*}Lp_{0}=\Delta u_{0}}$ in ${\displaystyle \Omega }$ and ${\displaystyle \partial _{n}p_{0}=0}$ on ${\displaystyle \partial \Omega }$

Thanks to the topological gradient, it is possible to detect the edges and their orientation and to define an appropriate ${\displaystyle C}$ for the image reconstruction process.[11]

In image processing, the topological derivatives have also been studied in the case of a multiplicative noise of gamma law or in presence of Poissonian statistics.[14]

Inverse problems

In 2009, the topological gradient method has been applied to tomographic reconstruction.[15] The coupling between the topological derivative and the level set has also been investigated in this application.[16]

References

1. J. Sokolowski and A. Zochowski, 44On topological derivative in shape optimization44, 1997
2. Topological Derivatives in Shape Optimization, Jan Sokołowski, May 28, 2012. Retrieved November 9, 2012
3. G. Allaire, F. Jouve, Coupling the level set method and the topological gradient in structural optimization, in IUTAM symposium on topological design optimization of structures, machines and materials, M. Bendsoe et al. eds., pp3-12, Springer (2006).
4. S. Amstutz, I. Horchani, and M. Masmoudi. Crack detection by the topological gradient method. Control and Cybernetics, 34(1):81–101, 2005.
5. S. Amstutz, A.A. Novotny, Topological asymptotic analysis of the Kirchhoff plate bending problem. ESAIM: COCV 17(3), pp. 705-721, 2011
6. L. J. Belaid, M. Jaoua, M. Masmoudi, and L. Siala. Image restoration and edge detection by topological asymptotic expansion. CRAS Paris, 342(5):313–318, March 2006.
7. D. Auroux and M. Masmoudi. Image processing by topological asymptotic analysis. ESAIM: Proc. Mathematical methods for imaging and inverse problems, 26:24–44, April 2009.
8. D. Auroux, M. Masmoudi, and L. Jaafar Belaid. Image restoration and classification by topological asymptotic expansion, pp. 23–42, Variational Formulations in Mechanics: Theory and Applications, E. Taroco, E.A. de Souza Neto and A.A. Novotny (Eds), CIMNE, Barcelona, Spain, 2007.
9. D. Auroux and M. Masmoudi. A one-shot inpainting algorithm based on the topological asymptotic analysis. Computational and Applied Mathematics, 25(2-3):251–267, 2006.
10. D. Auroux and M. Masmoudi. Image processing by topological asymptotic expansion. J. Math. Imaging Vision, 33(2):122–134, February 2009.
11. S. Larnier, J. Fehrenbach and M. Masmoudi, The topological gradient method: From optimal design to image processing, Milan Journal of Mathematics, vol. 80, issue 2, pp. 411–441, December 2012.
12. D. Auroux, L. Jaafar Belaid, and B. Rjaibi. Application of the topological gradient method to color image restoration. SIAM J. Imaging Sci., 3(2):153–175, 2010.
13. S. Larnier and J. Fehrenbach. Edge detection and image restoration with anisotropic topological gradient. In 2010 IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP), pages 1362–1365, March 2010.
14. A. Drogoul, G. Aubert, The topological gradient method for semi-linear problems and application to edge detection and noise removal.
15. D. Auroux, L. Jaafar Belaid, and B. Rjaibi. Application of the topological gradient method to tomography. In ARIMA Proc. TamTam'09, 2010.
16. T. Rymarczyk, P. Tchórzewski, J. Sikora, Topological Approach to Image Reconstruction in Electrical Impedance Tomography, ADVCOMP 2014 : The Eighth International Conference on Advanced Engineering Computing and Applications in Science

Books

A. A. Novotny and J. Sokolowski, Topological derivatives in shape optimization, Springer, 2013.

External links

• Allaire and al. Structural optimization using topological and shape sensitivity via a level set method