Nonlinear diffusion filtering of data on the Earth’s surface
Tóm tắt
The paper deals with data filtering on closed surfaces using linear and nonlinear diffusion equations. We define a surface finite-volume method to approximate numerically parabolic partial differential equations on closed surfaces, namely on a sphere, ellipsoid or the Earth’s surface. The closed surface as a computational domain is approximated by a polyhedral surface created by planar triangles and we construct a dual co-volume grid. On the co-volumes we define a weak formulation of the problem by applying Green’s theorem to the Laplace–Beltrami operator. Then the finite-volume method is applied to discretize the weak formulation. Weak forms of elliptic operators are expressed through surface gradients. In our numerical scheme we use a piece-wise linear approximation of a solution in space and the backward Euler time discretization. Furthermore, we extend a linear diffusion on surface to the regularized surface Perona–Malik model. It represents a nonlinear diffusion equation, which at the same time reduces noise and preserves main edges and other details important for a correct interpretation of the real data. We present four numerical experiments. The first one has an illustrative character showing how an additive noise is filtered out from an artificial function defined on a sphere. Other three examples deal with the real geodetic data on the Earth’s surface, namely (i) we reduce a stripping noise from the GOCE satellite only geopotential model up to degree 240, (ii) we filter noise from the real GOCE measurements (the component
$$T_{zz})$$
, and (iii) we reduce a stripping noise from the satellite only mean dynamic topography at oceans. In all experiments we focus on a comparison of the results obtained by both the linear and nonlinear models presenting advantages of the nonlinear diffusion.
Tài liệu tham khảo
Alvarez L, Guichard F, Lions PL, Morel JM (1993) Axioms and fundamental equations of image processing. Arch Ration Mech 123(3):199–257
Barrett R, Berry M, Chan TF, Demmel J, Donato J, Dongarra J, Eijkhout V, Pozo R, Romine C, Van der Vorst H (1994) Templates for the solution of linear systems: building blocks for iterative methods, 2nd edn. SIAM, Philadelphia. http://www.netlib.org/linalg/html_templates/Templates.html
Bruinsma SL, Marty JC, Balmino G, Biancale R, Foerste C, Abrikosov O, Neumayer H (2010) GOCE gravity field recovery by means of the direct numerical method. Presented at the ESA living planet symposium, Bergen, Norway, 27th June–2nd July (2010). see also http://earth.esa.int/GOCE
Capuzzo Dolcetta I, Ferretti R (2001) Optimal stopping time formulation of adaptive image filtering. Appl Math Optim 43(3):245–258
Caselles V, Kimmel R, Sapiro G (1995) Geodesic active contours. In: Proceedings of the international conference on computer vision ’95, Boston, pp 694–699
Caselles V, Kimmel R, Sapiro G (1997) Geodesic active contours. Int J Comput Vis 22:61–79
Catté F, Lions PL, Morel JM, Coll T (1992) Image selective smoothing and edge detection by nonlinear diffusion. SIAM J Numer Anal 29:182–193
Chen Y, Vemuri BC, Wang L (2000) Image denoising and segmentation via nonlinear diffusion. Comput Math Appl 39:131–149
Čunderlík R, Mikula K, Mojzeš M (2008) Numerical solution of the linearized fixed gravimetric boundary-value problem. J Geod 82: 15–29
Dziuk G, Elliott CM (2007) Surface finite elements for parabolic equations. J Comput Math 25(4):385–407
ESA (2011) The GOCE user toolbox, version 2.0. http://earth.esa.int/gut
Kichenassamy S, Kumar A, Olver P, Tannenbaum A, Yezzi A (1995) Gradient flows and geometric active contours models. In: Proceedings of the international conference on computer vision’95, Boston
Kichenassamy S, Kumar A, Olver P, Tannenbaum A, Yezzi A (1996) Conformal curvature flows: from phase transitions to active vision. Arch Ration Mech Anal 134:275–301
Koenderink J (1984) The structure of images. Biol Cybern 50:363–370
Mikula K (2002) Image processing with partial differential equations. In: Bourlioux A, Gander M (eds) Modern methods in scientific computing and applications. NATO science series II, vol 75. Kluwer, Dordrecht, pp 283–322
Mikula K, Ramarosy N (2001) Semi-implicit finite volume scheme for solving nonlinear diffusion equations in image processing. Numerische Mathematik 89(3):561–590
Pavlis NK, Holmes SA, Kenyon SC, Factor JK (2008) An earth gravitational model to degree 2160: EGM2008. Presented at the 2008 general assembly of EGU, Vienna, Austria, April 13–18
Perona P, Malik J (1987) Scale space and edge detection using anisotropic diffusion. In: Proceedings of the IEEE society workshop on computer vision
Witkin AP (1983) Scale-space filtering. In: Proceedings of the 8th international joint conference on artificial intelligence, Karlsruhe, Germany, vol 2, pp 1019–1022