Next Article in Journal
Entropy: An Inspiring Tool for Characterizing Turbulence–Combustion Interaction in Swirling Flames via Direct Numerical Simulations of Non-Premixed and Premixed Flames
Next Article in Special Issue
Patterns in Temporal Networks with Higher-Order Egocentric Structures
Previous Article in Journal
Towards Multifractality through an Ernst-Type Potential in Complex Systems Dynamics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Probabilistic Result on Impulsive Noise Reduction in Topological Data Analysis through Group Equivariant Non-Expansive Operators

Department of Mathematics, University of Bologna, 40126 Bologna, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Entropy 2023, 25(8), 1150; https://doi.org/10.3390/e25081150
Submission received: 5 May 2023 / Revised: 22 July 2023 / Accepted: 25 July 2023 / Published: 31 July 2023

Abstract

:
In recent years, group equivariant non-expansive operators (GENEOs) have started to find applications in the fields of Topological Data Analysis and Machine Learning. In this paper we show how these operators can be of use also for the removal of impulsive noise and to increase the stability of TDA in the presence of noisy data. In particular, we prove that GENEOs can control the expected value of the perturbation of persistence diagrams caused by uniformly distributed impulsive noise, when data are represented by L-Lipschitz functions from R to R .
MSC:
55N31; 62R40; 60-08; 65D18; 68T09; 68U05

1. Introduction

In the last thirty years, Topological Data Analysis (TDA) developed as a useful mathematical theory for analyzing data [1], benefiting from the dimensionality reduction guaranteed by topology [2,3,4]. Topology can construct low-dimensional representations of high-dimensional manifolds of data, thereby reducing the dimensionality of the parameter space and ultimately reducing the learning cost of machine learning models. One of the main tools in TDA is the concept of the persistence diagram, which is a collection of points in the real plane that describes the homological changes of the sublevel sets of suitable continuous functions. These changes give important information about the data of interest, focusing on some of their most relevant properties. Persistence diagrams can be used in the presence of noise, since a well-known stability theorem states that these topological descriptors change in a controlled way when we know that the functions expressing the filtrations we are interested in change in a controlled way with respect to the sup-norm. More precisely, the sup-norm distance between two functions is an upper bound for the bottleneck distance between the corresponding persistence diagrams [5]. Furthermore, L p -stability of persistence diagrams with respect to the sup-norm have been proved in [6]. Unfortunately, in many applications, the sup-norm of the noise is not guaranteed to be small, and hence these results cannot be directly applied. In particular, these results cannot be directly used when data are represented by functions from R to R and affected by impulsive noise, i.e., noise that carries a sudden sharp change of short duration, when the variable is seen as a time. It is indeed well known that persistence diagrams can completely change when the data are subject to impulsive noise.
Analogously, in the discrete setting of TDA, the presence of outliers in cloud points can drastically affect the corresponding persistence diagrams. The problem of managing outliers has been studied by several authors with different techniques. In [7], an approach based on confidence sets has been introduced. A method inspired by the k-nearest neighbors regression and the local median filtering has been used in [8], while the concept of the bagplot has been applied in [9]. In [10], an approach based on reproducing kernels has been proposed.
In our paper, we start exploring a different probabilistic approach, in the topological setting. The main idea is that studying the properties of data should primarily be based on the analysis of the observers that are asked to examine the data, since we cannot ignore that different observers can differently judge the same data. This approach has been initially proposed in [11] and requires both the definition of the space of group equivariant non-expansive operators (GENEOs) and the development of geometrical techniques to move around in this space [12,13]. In other words, in this model we should not wonder how we have to manage the data but rather how we have to manage the observers (i.e., GENEOs) analyzing the data. As a first step in this direction, in this paper we show how GENEOs could be used to obtain stability of persistence diagrams of 1D signals in the presence of impulsive noise. For the use of such operators in the comparison of 1D signals under the action of several transformation groups, we refer the interested reader to the paper [14]. This type of signal is important in many applications, such as those concerning EEG data and time series.
GENEOs have been studied in [14] as a new tool in TDA, since they allow for an extension of the theory that is not invariant under the action of every homeomorphism of the considered domain. This is important in applications where the invariance group is not the group of all homeomorphisms, such as the ones concerning shape comparison. Interestingly, GENEOs are also deeply related to the foliation method used to define the matching distance in two-dimensional persistent homology [15,16] and can be seen as a theoretical bridge between TDA and Machine Learning [11]. Furthermore, these operators make available lower bounds for the natural pseudo-distance d G ( φ 1 , φ 2 ) : = inf g G φ 1 φ 2 g , associated with a group G of self-homeomorphisms of the domain of the signals φ 1 , φ 2 [14]. For general information about the interest in the theory of GENEOs and its applications, we refer the reader to [17].
In our paper, we prove that GENEOs can control the expected value of the perturbation of persistence diagrams caused by uniformly distributed impulsive noise when data are represented by L-Lipschitz functions from R to R . In order to do that, we choose a “mother” bump function ψ , i.e., a continuous function that is non-negative, upper bounded by 1, and compactly supported, and we assume that our noise is made up of finite bumps each obtained by translating, heightening, and/or widening ψ . The function φ ^ = φ + R represents the corrupted data, where the noise is given by the function R = i = 1 k a i ψ ( b i ( x c i ) ) for some a i , b i , c i R , with b i > 0 , and some positive integer k.
In this situation, trying to use a convolution to approximate our starting data is not effective, because even if it does contract the bumps, it does not cut them and hence does not improve the sup-norm distance from the original data. A classical approach for the removal of impulsive noise is the one of using a median filter [18]. Although this approach would be quite efficient in the discrete case, let us remark that in our setting, we are considering continuous functions. The median for the continuous case is defined as the interval of those values m such that
m f ( x ) d x = m + f ( x ) d x ,
where f is a density of probability. However, this operator is not stable: a small alteration in the starting function could lead to a significant change in the median.
The operators we consider in this paper are F δ ( φ ) ( x ) = max ( φ ( x δ ) , φ ( x + δ ) ) and F ε ( φ ) ( x ) = min ( φ ( x ε ) , φ ( x + ε ) ) . The main idea is that F δ removes the noise “directed downwards” and F ε the noise “directed upwards”, and hence, their composition should be able to eliminate all the bumps. These operators are GENEOs with respect to isometries of the real line. We prove that moving in the space of GENEOs by taking suitable values for ε and δ , we can get quite close to restoring the original function φ , depending on how the bumps are positioned. The closer the bumps are to being Dirac delta functions and the further they are from each other, the better our approximation can be. On the ground of this result, we finally obtain an estimate of the expected value E ( F δ F ε ( φ ^ ) φ ) .
We conclude this introduction by observing that, while the operators we use in this paper are elementary, the theory of GENEOs makes available procedures to learn more complex operators by combining simple GENEOs through techniques belonging to the rising field of Geometric Deep Learning [11,12,13]. For this reason, we hope that our results can not only illustrate a practical use of some particular GENEOs but also be a first step towards the application of Machine Learning for increasing the stability of TDA in the presence of noise.
The paper is structured as follows. In Section 2, the mathematical background is laid. The case we consider and our notations are explained in Section 3. In Section 4, we prove the results that are needed in order to demonstrate our main results. In Section 5, the main theorems giving us probabilistic upper bounds are formulated. In Section 6, some examples and experiments are presented in order to better illustrate the use of our results. A brief discussion concludes the paper.

2. Mathematical Setting

In this section, we will recall some basic concepts we will use in this paper.
To avoid misunderstandings about pairs and intervals, we will use the symbols ( u , v ) , [ u , v ] , and ] u , v [ to denote a pair, a closed interval, and an open interval, respectively.

2.1. Representing Data as Real Functions

Let us consider a set Φ of bounded functions from a set X to R , which will represent the data we wish to take into account (e.g., functions describing sounds). We shall call Φ the set of admissible measurements on X. We endow Φ with the topology induced by the sup-norm · and the corresponding distance D Φ . A pseudo-metric D X can be defined on X by setting D X ( x 1 , x 2 ) = sup φ Φ | φ ( x 1 ) φ ( x 2 ) | for every x 1 , x 2 X . We recall that a pseudo-metric on a set X is a distance d without the property d ( x , y ) = 0 x = y (in other words, we do not have the so-called “identity of indiscernibles”). We will consider the topological space ( X , τ D X ) where τ D X is the topology induced by D X . A base for this topology is given by the set of open balls { B ( x , r ) : x X , r R } , where B ( x , r ) : = { x X : D X ( x , x ) < r } . The choice of this topology makes every function in Φ a continuous functions. As shown in [11], this fact enables us to use persistence diagrams in the study of Φ .

2.2. GENEOs as Operators Acting on Data

We are interested in considering transformations of data. Let Homeo Φ ( X ) be the group of Φ -preserving homeomorphisms from X to X with respect to the topology τ D X , meaning that every g in Homeo Φ ( X ) is a homeomorphism of X such as both φ g and φ g 1 belong to Φ for every φ in Φ . Let G be a subgroup of Homeo Φ ( X ) . G represents the set of transformations on data for which we will require equivariance to be respected.
Under the previously stated assumptions, we call the ordered pair ( Φ , G ) a perception pair. We can now introduce the concept of GENEO.
Definition 1. 
Let ( Φ , G ) and ( Ψ , H ) be perception pairs and assume that a homomorphism T : G H is given. A function F : Φ Ψ is called a Group Equivariant Non-Expansive Operator (GENEO) from ( Φ , G ) to ( Ψ , H ) with respect to T if the following properties hold:
1. 
(Group Equivariance) F ( φ g ) = F ( φ ) T ( g ) for every φ Φ , g G ;
2. 
(Non-Expansivity) D Ψ ( F ( φ 1 ) , F ( φ 2 ) ) D Φ ( φ 1 , φ 2 ) for every φ 1 , φ 2 Φ .
Let us now consider the set F all of all GENEOs from ( Φ , G ) to ( Ψ , H ) with respect to T : G H . The space F all is endowed with the extended pseudo-metric D F all , defined by setting D F all ( F 1 , F 2 ) = sup φ Φ D Ψ ( F 1 ( φ ) , F 2 ( φ ) ) for every F 1 , F 2 F all . The word extended refers to the possibility that D F all takes an infinite value.
The following result can be proven [11]:
Theorem 1. 
If ( Φ , D Φ ) , ( Ψ , D Ψ ) are compact and convex then the metric space ( F all , D F all ) is compact and convex.
If a non-empty set F F all is fixed, we can define the following pseudo-distance D F , Φ on Φ :
Definition 2. 
For any φ 1 , φ 2 in Φ we set
D F , Φ ( φ 1 , φ 2 ) = sup F F F ( φ 1 ) F ( φ 2 ) .
This pseudo-distance allows us to compare data by taking into account how agents operate on data. Notice how, if G becomes larger, the natural pseudo-distance d G becomes harder to compute but this new pseudo-distance D F , Φ becomes easier to evaluate. In order to find a lower bound for D F , Φ it is useful to introduce the notion of persistence diagram.

2.3. Persistence Diagrams

We will now recall some basic definitions and results in persistent homology. The interested reader can find more details in [4].
Let us consider an ordered pair ( X , φ ) where X is a topological space and φ : X R is a continuous function. For any t R we can set X t : = φ 1 ( , t ) . If u < v , the inclusion i u , v : X u X v induces a homomorphism i u , v k : H k ( X u ) H k ( X v ) between the kth homology groups of X u and X v . We can define the kth persistent homology group, with respect to φ and computed at the point ( u , v ) , as P H k ( u , v ) : = i u , v k ( H k ( X u ) ) . Moreover, we can define the kth persistent Betti numbers function r k ( u , v ) as the rank of P H k ( u , v ) .
The kth persistent Betti numbers function can be represented by the kth persistence diagram. This diagram is defined as the multi-set of all the ordered pairs ( u j , v j ) , where u j and v j are the times of birth and death of the jth k-dimensional hole in X, respectively. We call time of birth of a hole the first time at which the homology class appears, and time of death the first time at which the homology class merges with an older one. When a hole never dies, we set its time of death equal to . We also add to this set all points of the form ( w , w ) for w R .
In Figure 1 the filtration of the set X : = 0 , 3 4 π given by the function φ ( x ) = 2 sin x is illustrated. In this example, the topology on X is the one defined by the space Φ of admissible functions given by all functions a sin ( x + b ) with a , b R . The reader can easily check that this topology is the same as the Euclidean topology. The persistence diagram in degree k = 0 of the function φ is displayed in Figure 2.

2.4. Comparing Persistence Diagrams

Persistence diagrams can be efficiently compared by means of a suitable metric d match . In order to define it, we first define the pseudo-distance
δ ( ( x , y ) , ( x , y ) ) = min { max { | x x | , | y y | } , max { | x y | / 2 , | x y | / 2 } }
for all ( x , y ) , ( x , y ) { ( x , y ) R with x y } { ( x , ) with x R } by agreeing that y = , y = for y , = 0 , / 2 = , | ± | = , min { , c } = c , max { , c } = .
If two persistence diagrams, D , D , are given, we can set
d match ( D , D ) = inf σ Σ sup P D δ ( P , σ ( P ) )
where Σ represents the set of all bijections between the multisets D , D . This function is usually called the bottleneck distance between D and D .
For every degree k we can now define a new pseudo-metric:
D F , Φ match ( φ 1 , φ 2 ) = sup F F d match ( D F ( φ 1 ) , D F ( φ 2 ) )
where D F ( φ 1 ) , D F ( φ 2 ) are the persistence diagrams at degree k of the functions F ( φ 1 ) , F ( φ 2 ) , respectively.
In this paper we will limit ourselves to considering data represented as functions from R to R , and we recall that for this kind of data, persistence diagrams are non-trivial only in degree k = 0 (i.e., when persistent homology is used to count connected components). For this reason, in the following we will always assume k = 0 .

3. Our Model

In this paper, we will be mainly interested in the set Lip L of all L-Lipschitz functions from R to R , for some fixed constant L R , and in the set C 0 ( R ) of all functions from R to R that are continuous with respect to the Euclidean topology. We will set X = R and consider the perception pairs ( Lip L , G ) , ( C 0 ( R ) , G ) , where G is the group of the Euclidean isometries of R .
We will assume our noise originates from a finite number of copies of a “mother” non-negative continuous bump function ψ : R R , such that supp ( ψ ) σ , σ for some σ > 0 and ψ 1 . We recall that the support of a function is the closure of the set of points where the function is non-zero. After fixing two positive real numbers, η and β , the noise we will be adding is a function R belonging to the space R η , β that contains the null function and all functions of the form i = 1 k a i ψ ( b i ( x c i ) ) , where k is a positive integer, and a i , b i , c i are real numbers such that | c i c j | η for i j and b i β for every index i. For any R R η , β , we define S ( R ) : = i = 1 k c i σ β , c i + σ β and remark that if x S ( R ) , then R ( x ) = 0 .
Figure 3 shows how drastically persistence diagrams can change when the data are subject to impulsive noise. Our purpose will be to recover φ Lip L as well as possible from the function φ ^ = φ + R . An example of such a situation is depicted in Figure 4.
The following result will be of use.
Proposition 1. 
Let F 1 , F 2 be GENEOs from ( C 0 ( R ) , G ) to ( C 0 ( R ) , G ) with respect to the trivial homomorphism T = id : G G . Then F 1 F 2 is a GENEO from ( C 0 ( R ) , G ) to ( C 0 ( R ) , G ) with respect to T.
Proof. 
For every φ C 0 ( R ) , g G we have that
F 1 F 2 ( φ g ) = F 1 ( F 2 ( φ g ) ) = F 1 ( F 2 ( φ ) g ) = F 1 ( F 2 ( φ ) ) g = ( F 1 F 2 ) ( φ ) g .
Therefore, F 1 F 2 is G-equivariant. Moreover, for any φ 1 , φ 2 C 0 ( R )
D C 0 ( R ) ( F 1 F 2 ( φ 1 ) ) , F 1 F 2 ( φ 2 ) ) = D C 0 ( R ) ( F 1 ( F 2 ( φ 1 ) ) , F 1 ( F 2 ( φ 2 ) ) ) D C 0 ( R ) ( F 2 ( φ 1 ) ) , F 2 ( φ 2 ) ) D C 0 ( R ) ( φ 1 , φ 2 ) .
It follows that F 1 F 2 is non-expansive. □

4. Cutting Off the Noise by GENEOs

We start by introducing two families of GENEOs from ( C 0 ( R ) , G ) to ( C 0 ( R ) , G ) with respect to the identical homomorphism.
Definition 3. 
Let φ Lip L and ε > 0 . For all x R we define:
1. 
F ε ( φ ) ( x ) = max ( φ ( x ε ) , φ ( x + ε ) ) ;
2. 
F ε ( φ ) ( x ) = min ( φ ( x ε ) , φ ( x + ε ) ) .
Proposition 2. 
The maps F ε and F ε are GENEOs from ( C 0 ( R ) , G ) to ( C 0 ( R ) , G ) with respect to the identical homomorphism.
Proof. 
We start by proving that F ε is G-equivariant.
Let g be the translation x x + k with k R , then
F ε ( φ g ) = max { φ ( ( x + k ) ε ) , φ ( ( x + k ) + ε ) } = max { φ ( ( x ε ) + k ) , φ ( ( x + ε ) + k ) } = F ε ( φ ) g .
Let g be the symmetry x x , then
F ε ( φ g ) = max { φ ( ( x ) ε ) , φ ( ( x ) + ε ) } = max { φ ( ( x + ε ) ) , φ ( ( x ε ) ) } = max { φ ( x + ε ) , φ ( x ε ) } g = F ε ( φ ) g .
Since every isometry in G can be written as the composition of a symmetry and a translation, our statement follows.
Furthermore, for any x R
| F ε ( φ 1 ) ( x ) F ε ( φ 2 ) ( x ) | = = | max { φ 1 ( x ε ) , φ 1 ( x + ε ) } max { φ 2 ( x ε ) , φ 2 ( x + ε ) } | max { | φ 1 ( x ε ) φ 2 ( x ε ) | , | φ 1 ( x + ε ) φ 2 ( x + ε ) | } max { φ 1 φ 2 , φ 1 φ 2 } = φ 1 φ 2 .
It follows that F ε ( φ 1 ) F ε ( φ 2 ) φ 1 φ 2 , and hence F ε is non-expansive.
We observe that min { a , b } = max { a , b } , and hence F ε ( φ ) = F ε ( φ ) for every φ C 0 ( R ) . Moreover, the map that takes φ to φ is a GENEO from ( C 0 ( R ) , G ) to ( C 0 ( R ) , G ) with respect to the identical homomorphism. Therefore, Proposition 1 ensures that G-equivariance and non-expansivity also hold for F ε . □
Since Proposition 1 shows that the composition of GENEOs is still a GENEO, the operator F δ F ε is a GENEO from ( C 0 ( R ) , G ) to ( C 0 ( R ) , G ) with respect to the identical homomorphism.
We now want to prove that if a function φ ^ is obtained by adding impulsive noise to a function φ , then the value of
F δ F ε ( φ ^ ) φ
is bounded, and possibly small, provided that δ and ε are chosen appropriately. The main idea is that the operator F ε cuts the noise “directed upwards” and F δ cuts the noise “directed downwards”.
In order to proceed, we need two lemmas.
Lemma 1. 
Let R C 0 ( R ) , φ Lip L for some L R , and set φ ^ : = φ + R . Then for any ε > 0 and δ > 0
( i )
L ε + F ε ( R ) F ε ( φ ^ ) φ L ε + F ε ( R ) ;
( ii )
L δ + F δ ( R ) F δ ( φ ^ ) φ L δ + F δ ( R ) ;
( iii )
L ( δ + ε ) + F δ F ε ( R ) F δ F ε ( φ ^ ) φ L ( δ + ε ) + F δ F ε ( R ) .
Proof. 
Since φ is Lipschitz of constant L, we have that for any value x R   | φ ( x ε ) φ ( x ) | L ε and | φ ( x + ε ) φ ( x ) | L ε . Therefore,
F ε ( φ ^ ) ( x ) = F ε ( φ + R ) ( x ) = min { φ ( x ε ) + R ( x ε ) , φ ( x + ε ) + R ( x + ε ) } min { φ ( x ) + L ε + R ( x ε ) , φ ( x ) + L ε + R ( x + ε ) } = φ ( x ) + L ε + min { R ( x ε ) , R ( x + ε ) } = φ ( x ) + L ε + F ε ( R ) ( x ) .
Analogously, F ε ( φ ^ ) ( x ) φ ( x ) L ε + F ε ( R ) ( x ) . The same steps applied to F δ yield the second statement of the lemma. As for the last claim, we can see that:
F δ F ε ( φ ^ ) ( x ) = max { F ε ( φ ^ ) ( x δ ) , F ε ( φ ^ ) ( x + δ ) } max { φ ( x δ ) + L ε + F ε ( R ) ( x δ ) , φ ( x + δ ) + L ε + F ε ( R ) ( x + δ ) } max { φ ( x ) + L δ + L ε + F ε ( R ) ( x δ ) , φ ( x ) + L δ + L ε + F ε ( R ) ( x + δ ) } = φ ( x ) + L δ + L ε + max { F ε ( R ) ( x δ ) , F ε ( R ) ( x + δ ) } = φ ( x ) + L δ + L ε + F δ F ε ( R ) ( x ) .
Analogously, we can prove the lower bound. □
Henceforth, we will assume that any summation on an empty set of indexes is the null function.
Lemma 2. 
Let R R η , β and λ σ β . If λ ρ η 2 λ then
( a )
F ρ ( R ) ( x ) = a i < 0 a i ψ ( b i ( x ρ c i ) ) + a i ψ ( b i ( x + ρ c i ) ) 0
( b )
F ρ ( R ) ( x ) = a i > 0 a i ψ ( b i ( x ρ c i ) ) + a i ψ ( b i ( x + ρ c i ) ) 0
for all x R . Moreover F ρ ( R ) , F ρ ( R ) ( x ) R 2 λ , β .
Proof. 
We will suppose without loss of generality that c i < c i + 1 for all i = 1 , , k 1 and a i 0 for all indexes i.
We want to show that at least one of x ρ and x + ρ must always be outside of i = 1 k c i λ , c i + λ . If by contradiction both x ρ and x + ρ were in the same c i λ , c i + λ , then we would have ρ < λ , against our hypotheses. Moreover, x ρ and x + ρ cannot belong to different intervals c i λ , c i + λ and c j λ , c j + λ for i < j : if by contradiction x ρ c i λ , c i + λ and x + ρ c j λ , c j + λ for some i < j , then 2 ρ = x + ρ ( x ρ ) > ( c j λ ) ( c i + λ ) c j c i 2 λ η 2 λ , and hence we would have ρ > η 2 λ , against our hypotheses.
Since c i σ β , c i + σ β c i λ , c i + λ , it follows that at least one among the values R ( x ρ ) and R ( x + ρ ) must always be zero. Let us now set I i : = ( c i ρ ) λ , ( c i ρ ) + λ and I i + : = ( c i + ρ ) λ , ( c i + ρ ) + λ . These two intervals must be disjoint, since ( c i + ρ λ ) ( c i ρ + λ ) = 2 ρ 2 λ 0 .
Let us now consider { I 1 , I 1 + , , I k , I k + } . We will now prove that any two distinct elements from this set must be a disjoint. Since we have just proven that I i I i + = for i = 1 , , k , the following holds in the case i < j :
1.
I i + I j + = , since c j + ρ λ ( c i + ρ + λ ) c j c i 2 λ η 2 λ 0 ;
2.
I i I j = , since c j ρ λ ( c i ρ + λ ) c j c i 2 λ η 2 λ 0 ;
3.
I i + I j = , since c j ρ λ ( c i + ρ + λ ) c j c i 2 ρ 2 λ η 2 η 2 λ 2 λ = 0 ;
4.
I j + I i , since c j + ρ λ ( c i ρ + λ ) c j c i + 2 ρ 2 λ η + 2 λ 2 λ 0 .
Let us now fix k { 1 , , k } .
Since supp ( a k ψ ( b k ( x c k ) ) ) c k λ , c k + λ , then supp ( a k ψ ( b k ( x ρ c k ) ) ) I k and supp ( a k ψ ( b k ( x + ρ c k ) ) ) I k + .
Hence, for any x R we have that
F ρ ( R ) ( x ) = min { R ( x ρ ) , R ( x + ρ ) } = 0 i f   x i = 1 k ( I i I i + ) a j ψ ( b j ( x + ρ c j ) ) i f   x I j and a j < 0 a j ψ ( b j ( x ρ c j ) ) i f   x I j + and a j < 0 .
This means that
F ρ ( R ) ( x ) = a i < 0 a i ψ ( b i ( x ρ c i ) ) + a i ψ ( b i ( x + ρ c i ) ) 0 .
Now, given c i + ρ and c j ρ centers of bumps of F ρ ( R ) , we have that
| ( c i + ρ ) ( c j ρ ) | min { η 2 ρ , 2 ρ } 2 λ .
It follows that F ρ ( R ) R 2 λ , β .
By noting that F ρ ( R ) = F ρ ( R ) , we obtain the second part of the thesis. □
Let us remark that, in particular, if λ = 2 σ β , then F ρ ( R ) , F ρ ( R ) R 4 σ β , β . We observe that for there to exist a ρ that satisfies the hypotheses of Lemma 2, the inequality η 4 λ must hold.
We are now actually ready to prove a result that will be useful in the remainder of our paper.
Theorem 2. 
Given θ , β , L R , let R ¯ R θ , β , φ Lip L and set φ ^ : = φ + R ¯ . If 2 σ β ε θ 2 2 σ β then for any σ β δ 1 2 min { θ 2 ε , 2 ε } σ β the following inequality holds:
F δ F ε ( φ ^ ) φ L ( ε + δ ) .
Proof. 
Let x R . Firstly, let us notice that the condition 2 σ β ε θ 2 2 σ β implies that σ β 1 2 min { θ 2 ε , 2 ε } σ β . From Lemma 1, we know that
L ( δ + ε ) + F δ F ε ( R ¯ ) ( x ) F δ F ε ( φ ^ ) ( x ) φ ( x ) L ( δ + ε ) + F δ F ε ( R ¯ ) ( x ) .
Let us now prove that F δ F ε ( R ¯ ) ( x ) = 0 .
By applying Lemma 2 with η : = θ , λ : = 2 σ β , ρ : = ε and R : = R ¯ , we obtain
F ε ( R ¯ ) ( x ) = a i < 0 a i ψ ( b i ( x ε c i ) ) + a i ψ ( b i ( x + ε c i ) ) 0 .
Let us remark that F ε ( R ¯ ) R 4 σ β , β . Moreover,
1.
F δ F ε ( R ¯ ) ( x ) = max { F ε ( R ¯ ) ( x δ ) , F ε ( R ¯ ) ( x + δ ) } 0 since both terms are negative;
2.
F δ F ε ( R ¯ ) ( x ) = F δ ( F ε ( R ¯ ) ) ( x ) 0 .
These inequalities follow from Lemma 2, by setting η : = 4 σ β , λ : = σ β , ρ : = δ and R : = F ε ( R ¯ ) .
Therefore, we have proved that | F δ F ε ( φ ^ ) ( x ) φ ( x ) | L ( δ + ε ) for any x R . It follows that F δ F ε ( φ ^ ) φ L ( ε + δ ) . □
Let us remark that Theorem 2 works under the (only) implicit assumption that θ 8 σ β . In our setting, this should not be restrictive since it means that the noise added is made up of scattered, thin bumps, without any reference to the height of the bumps. This is what we expect when considering additive impulsive noise.
Corollary 1. 
Given θ , β , L R , let R R θ , β , φ Lip L , and set φ ^ : = φ + R . If θ 8 σ β , then
F σ β F 2 σ β φ ^ φ 3 L σ β .
Proof. 
The claim follows from Theorem 2 by taking ε : = 2 σ β and δ : = σ β .
Corollary 1 and the well-known stability of persistence diagrams with respect to the max-norm [5] immediately imply the following result, which is of interest in TDA (the symbol d match denotes the usual bottleneck distance between persistence diagrams).
Corollary 2. 
Given θ , β , L R , let R R θ , β , φ Lip L , and set φ ^ : = φ + R . If θ 8 σ β , and D and D are the persistence diagrams in degree 0 of the filtering functions φ and F σ β F 2 σ β ( φ ^ ) , respectively, then
d match ( D , D ) 3 L σ β .

5. Our Main Results

We are now ready to prove our main results. We start by recalling a simple statement concerning the probability p that any two distinct points in a randomly chosen set of cardinality k in an interval of length have a distance greater than η (cf., e.g., [19]). For the reader’s convenience, we report the proof here.
Lemma 3. 
Let X 1 , , X k , with k 2 , be independent random variables, uniformly distributed on the interval [ 0 , ] , for some > 0 . Let
M : = min 1 i , j k i j | X i X j |
be the minimal distance between two distinct random variables. Then, we have
P ( M > η ) = 1 if η 0 , 1 ( k 1 ) η k if 0 < η < k 1 , 0 if η k 1 .
Proof. 
It suffices to consider the case 0 < η < k 1 . By symmetry, we have
P ( M > η ) = k ! P ( ( M > η ) ( X 1 < X 2 < < X k ) )
(since X 1 , , X k are uniformly distributed)
= k ! Leb ( S ) k
where S = { x [ 0 , ] k x 1 < x 2 η < x 3 2 η < < x k ( k 1 ) η } , and Leb denotes the Lebesgue measure. Setting y i = x i ( i 1 ) η for i = 1 , , k , we have that Leb ( S ) = Leb ( S ) where
S = { y [ 0 , ( k 1 ) η ] k y 1 < y 2 < < y k } .
On the other hand, again by symmetry, we have
Leb ( S ) = Leb ( [ 0 , ( k 1 ) η ] k ) k ! = ( ( k 1 ) η ) k k !
and plugging this last identity into (1) we obtain the thesis. □
We can now prove the following result, concerning the expected value of the error F σ β F 2 σ β ( φ ^ ) φ . We remark explicitly that the error is random because it is a function of c 1 , , c k , which are independent random variables uniformly distributed on the interval [ 0 , ] .
Theorem 3. 
Let us choose a function φ Lip L , a non-negative continuous function ψ : R R with ψ 1 and supp ( ψ ) σ , σ for some σ > 0 , two positive numbers β and ℓ, and an integer k 2 . For i = 1 , , k , let us fix a i R and b i β , and set α ¯ : = max | a i | . Moreover, let c 1 , , c k be independent random variables, uniformly distributed on the interval [ 0 , ] . Let us consider the random variable φ ^ : = φ + R , where R ( x ) : = i = 1 k a i ψ ( b i ( x c i ) ) for any x R . If σ β < 8 ( k 1 ) , then
E F σ β F 2 σ β ( φ ^ ) φ 3 L σ β + k α ¯ 1 1 8 ( k 1 ) σ β k .
Proof. 
By setting δ = σ β and ε = 2 σ β in statement ( iii ) of Lemma 1, we have that
F σ β F 2 σ β ( φ ^ ) φ 3 L σ β + F σ β F 2 σ β ( R ) .
Since the operator F σ β F 2 σ β is non-expansive and F σ β F 2 σ β ( 0 ) = 0 , it follows that
F σ β F 2 σ β ( R ) R k α ¯ ψ k α ¯ .
Therefore, F σ β F 2 σ β ( φ ^ ) φ 3 L σ β + k α ¯ . If we apply Lemma 3 with η = 8 σ β , we obtain that R R 8 σ β , β with probability p : = 1 8 ( k 1 ) σ β k . If R R 8 σ β , β , we can apply Theorem 2 by setting δ = σ β , ε = 2 σ β and θ = 8 σ β , and hence we obtain that F σ β F 2 σ β ( φ ^ ) φ 3 L σ β with probability at least p. Since F σ β F 2 σ β ( φ ^ ) φ 3 L σ β + k α ¯ in any case, it follows that
E F σ β F 2 σ β ( φ ^ ) φ 3 L σ β p + 3 L σ β + k α ¯ 1 p = 3 L σ β + k α ¯ 1 p .
Theorem 3 and the well-known stability of persistence diagrams with respect to the max-norm [5] immediately imply the following result, which is of interest in TDA (the symbol d match denotes the usual bottleneck distance between persistence diagrams).
Corollary 3. 
Let us make the same assumptions of Theorem 3. Let D and D be the persistence diagrams in degree 0 of the filtering functions φ and F σ β F 2 σ β ( φ ^ ) , respectively. If σ β < 8 ( k 1 ) , then
E d match ( D , D ) 3 L σ β + k α ¯ 1 1 8 ( k 1 ) σ β k .
This result shows that the use of suitable GENEOs can make TDA (relatively) stable also in the presence of impulsive noise, under the assumptions we have considered in this paper.
Theorem 3 and Corollary 3 can be easily extended to the case that the number k of bumps is a random variable:
Theorem 4. 
Let k be a random variable that takes values in the subset S = { 2 , 3 , , K } of the positive integers, where we assume that the probability of the integer j S is p j . Let us choose a function φ Lip L , a non-negative continuous function ψ : R R with ψ 1 and supp ( ψ ) σ , σ for some σ > 0 , two positive numbers β and ℓ. Let c 1 k , , c k k be random variables that, conditioned to k, are independent and uniformly distributed on the interval [ 0 , ] . Moreover, for i = 1 , , k , let us fix a i k R and b i k β , and set α ¯ k : = max i | a i k | . Let us also consider the random variable φ ^ : = φ + R k , where R k ( x ) : = i = 1 k a i k ψ ( b i k ( x c i k ) ) for any x R . If σ β < 8 ( K 1 ) , then
E F σ β F 2 σ β ( φ ^ ) φ 3 L σ β + j = 2 K p j j α ¯ j 1 1 8 ( j 1 ) σ β j .
Proof. 
It follows immediately from Theorem 3. □
Corollary 4. 
Let us make the same assumptions of Theorem 4. Let D and D be the persistence diagrams in degree 0 of the filtering functions φ and F σ β F 2 σ β ( φ ^ ) , respectively. If σ β < 8 ( K 1 ) , then
E d match ( D , D ) 3 L σ β + j = 2 K p j j α ¯ j 1 1 8 ( j 1 ) σ β j .
Proof. 
It follows immediately from Corollary 3. □

6. Examples and Experiments

We will now validate our approach based on GENEOs by giving two examples and illustrating some experimental results.

6.1. Examples

In order to verify how our approach works, we will set τ n : = 1 1 n 2 σ β + 1 n θ 2 2 σ β and consider the upper bound F τ n 2 F τ n φ ^ φ 3 2 L τ n , obtained by applying Corollary 1. We observe that τ n 2 σ β for every index n, and lim n + τ n = 2 σ β . We will examine two examples that use the GENEOs F τ n 2 F τ n and show how our method based on such operators and the method based on convolutions differ, as for their capability in removing additive impulsive noise. Moreover, we will compare the actual error F τ n 2 F τ n φ ^ φ to its upper bound 3 2 L τ n , by running several simulations. The convolutions that will be applied in our examples use the functions T h : R R defined by setting
T h ( x ) : = h 2 i f 1 h x 1 h 0 o t h e r w i s e
for h > 0 . We will see that, although the convolution with such functions is also a GENEO, it will not be able to efficiently remove the noise.
Our noise function R will be constructed starting from the mother function ψ defined by setting ψ ( x ) : = e 1 1 1 x 2 for x 1 , 1 and ψ ( x ) : = 0 for x 1 , 1 . Using the notation introduced in the previous sections, we will set σ = 1.1 , thus satisfying the condition supp ( ψ ) σ , σ . The impulsive noise will be added in an interval , . The following parameters are considered, with these respective uniform distributions:
  • k Unif { 1 , , 10 }
  • a i Unif 100 , 100 for i = 1 , , k
  • b i Unif 0 , 100 for i = 1 , , k
  • c i Unif + σ β , σ β for i = 1 , , k .
We set β : = min i = 1 , , N b i , α ¯ : = max i = 1 , , N | a i | , and η : = min i j | c i c j | . After producing random values for the parameters N , a i , b i , c i , our algorithm checks whether η > 8 σ β ; otherwise it generates another set of parameters.

6.1.1. First Example

Let us consider the function
φ ( x ) : = sin x i f 4 π x 4 π 0 o t h e r w i s e
for x R . We observe that φ Lip L , for L = 1 . We will add noise in the interval [ 4 π , 4 π ] (meaning = 4 π ) and visualize the results in such an interval. Figure 5 illustrates how the function φ ^ looks compared to φ .
We will start by considering how well the convolution φ ^ T n can approximate the original function φ , when n goes from 3 to 100. From Figure 6, it is immediately apparent that the max-norm distance between φ ^ T n and φ remains quite large.
If we apply a convolution with T 1 n , for 3 n 100 , we obtain the results displayed in Figure 7, showing that all information represented by the function φ is progressively destroyed.
In contrast, if we apply the operator F τ n 2 F τ n φ ^ , for n = 3 , 5 , 20 , 100 , we obtain the results displayed in Figure 8. As we can see, this operator is much more efficient in removing the bumps and restoring the function φ . We recall that as n increases, τ n tends to 2 σ β , and our operator becomes more effective at removing the noise.
As a matter of fact, when we apply a convolution with the function T h and check the corresponding errors via the sup-norm, we obtain the results displayed in Figure 9. Since lim h + φ ^ T h = φ ^ and lim h + φ ^ T 1 h = 0 , we obtain that the errors tend to φ ^ φ = max i = 1 , , N | a i | = α ¯ and φ , respectively.
In contrast, if we apply the operator F τ n 2 F τ n , we obtain the results displayed in Figure 10. This shows that the upper bound for the error stated in Corollary 1 is quite tight. As we can expect, the best denoising is achieved when we replace τ n with 2 σ β (see Figure 11).
Finally, we executed 1000 simulations. In each simulation, we added random impulsive noise to the function φ to produce a function φ ^ . We then applied F σ β F 2 σ β to φ ^ to see how close the upper bound in Corollary 1 is to the actual error F σ β F 2 σ β ( φ ^ ) φ . The same parameters as in the beginning of this example were used. As we can see in Figure 12, the overestimation committed by our upper bound is often quite close to zero relative to the Lipschitz constant L of the function φ .
This suggests that this upper bound is quite accurate.

6.1.2. Second Example

Let us consider the function
φ ( x ) : = 1 1000 ( x 5 ) ( x 3 ) ( x + 1 ) ( x + 4 ) ( x + 5 ) i f 5 x 5 0 o t h e r w i s e
for x R . The coefficient 1 1000 was chosen so that the Lipschitz constant L would be comparable to the one in the previous example. In this example, L = 27 25 . Figure 13 illustrates how the function φ ^ looks compared to φ .
We will start by considering how well the convolution φ ^ T n can approximate the original function φ , when n goes from 3 to 100. From Figure 14, it is immediately apparent that the max-norm distance between φ ^ T n and φ remains quite large.
If we apply a convolution with T 1 n , for 3 n 100 , we obtain the results displayed in Figure 15, showing that all information represented by the function φ is progressively destroyed.
In contrast, if we apply the operator F τ n 2 F τ n φ ^ , for n = 3 , 5 , 20 , 100 , we obtain the results displayed in Figure 16. As we can see, this operator is much more efficient in removing the bumps and restoring the function φ .
When we apply a convolution with the function T h and check the corresponding errors via the sup-norm, we obtain the results displayed in Figure 17. As we have already seen, since lim h + φ ^ T h = φ ^ and lim h + φ ^ T 1 h = 0 , we obtain that the errors tend to φ ^ φ = max i = 1 , , k | a i | = α ¯ and φ , respectively.
In contrast, if we apply the operator F τ n 2 F τ n , we obtain the results displayed in Figure 18. This shows that the upper bound for the error stated in Corollary 1 is quite tight. As we can expect, the best denoising is achieved when we replace τ n with 2 σ β (see Figure 19).
Finally, we again executed 1000 simulations. We used the same methodology as in the previous case, but this time we considered the polynomial presented at the beginning of this example. As we can see in Figure 20, the overestimation committed by our upper bound is often quite close to zero relative to the Lipschitz constant L of the polynomial φ .

6.2. Experiments

In order to check how good the upper bound stated in Theorem 3 is, we have made the following experiment. In the first step, we fixed = 20 and σ = 1.1 , and assumed the following parameters to be given: L > 0 , N N , α > 0 , β > 0 , k N , a 1 , , a k ] 0 , α [ , b 1 , , b k > β , c 1 , , c k ] 0 , [ .
Firstly, we used the parameters L > 0 , N N to generate a random L-Lipschitz function φ in the following way. We randomly chose and sorted in ascending order N points x 1 , , x N in the open interval ] 0 , [ , with uniform distribution. Hence, we obtained the following decomposition: ] 0 , [ = ] 0 , x 1 ] ] x 1 , x 2 ] ] x N , [ . We defined our Lipschitz function φ to be 0 outside ] 0 , [ . After setting ( x 0 , y 0 ) = ( 0 , 0 ) and ( x N + 1 , y N + 1 ) = ( , 0 ) for i { 1 , , N + 1 } , the value y i of the function at x i was randomly chosen, with uniform distribution in an interval that allows for an L-Lipschitz extension to [ 0 , ] of the function, i.e.,
max { y i 1 L ( x i x i 1 ) , L ( x i ) } , min { y i 1 + L ( x i x i 1 ) , L ( x i ) } .
Finally, the graph of the Lipschitz function φ on [ 0 , ] was obtained by connecting each point ( x i 1 , y i 1 ) to ( x i , y i ) with a segment for i { 1 , , N + 1 } . We observe that φ constructed this way is an L-Lipschitz function.
Secondly, we used the parameters α > 0 , β > 0 , k N , a 1 , , a k ] 0 , α [ , b 1 , , b k > β , c 1 , , c k ] 0 , [ to generate a noise function as follows. We considered the mother function ψ defined by setting ψ ( x ) : = e 1 1 1 x 2 for x 1 , 1 and ψ ( x ) : = 0 for x 1 , 1 . For each L-Lipschitz function φ produced in the previously described way, we considered the function φ ^ = φ + i = 1 k a i ψ ( b i ( x c i ) ) .
In Figure 21 and Figure 22, some examples of the functions we have produced are displayed for N = 3 and N = 7 , respectively. In each figure, the functions are displayed without noise (left) and with added noise (right).
In the second step, we fixed = 20 and σ = 1.1 again, and considered a probabilistic model where α , β , L are given and the values N, k, a i , b i , c i are random variables. In this setting, we compared the noise φ ^ φ with the probabilistic upper bound stated in Theorem 3 and the value F σ β F 2 σ β ( φ ^ ) φ , which represents the reduced noise that we can obtain by applying our method. In order to average our results, for each triplet ( α , β , L ) with α { 50 , 55 , 60 , , 100 } , β { 3 , 4 , 5 , , 13 } , and L { 1 , 2 , , 10 } , we randomly generated 100 examples of an L-Lipschitz function φ and its noisy version φ ^ by randomly choosing the parameters N, k, a i , b i , c i according to the following distributions:
  • N Unif { 1 , , 10 }
  • k Unif { 1 , , 10 }
  • a i Unif 0 , α for i = 1 , , k
  • b i Unif β , 20 for i = 1 , , k
  • c i Unif 3 σ β , 3 σ β for i = 1 , , k .
Then, for the chosen values of each one of the three variables α , β , L , we computed the mean, with respect to the other two variables, of the average of the noise φ ^ φ and the average of the reduced noise F σ β F 2 σ β ( φ ^ ) φ obtained by our method, both evaluated for ( φ , φ ^ ) varying in the set of cardinality 100 that we produced. We also computed the mean of the probabilistic upper bound stated in Theorem 3 with respect to the same two variables. The results are displayed in Figure 23, Figure 24 and Figure 25. We remind the reader that α and β , respectively, express the maximum height and the thinness of the noise bumps, while L is a Lipschitz constant for each function φ that we are interested in. These results illustrate the effectiveness of using GENEOs in the reduction of impulsive noise.

7. Conclusions

In our paper, we have proved a stability property for persistence diagrams of functions from R to R in the presence of impulsive noise. This property shows that TDA can also be useful when noise drastically changes the topology of the sublevel sets of the filtering functions we are considering, and highlights some new potential interactions between TDA and the theory of GENEOs. The experimental section shows that our approach is indeed able to remove impulsive noise in most cases. It would be interesting to investigate the possibility of extending our method to real-valued functions defined on n-dimensional domains and allowing the function ψ to be a random variable, by selecting suitable GENEOs. We observe that the operators F ε and F ε introduced in Definition 3 can be easily adapted to the case of functions from R n to R , by setting F ε ( φ ) ( x ) = max { φ ( x + v ) | v R n , v = ε } and F ε ( φ ) ( x ) = min { φ ( x + v ) | v R n , v = ε } for all x R n . Applying a suitable operator F δ F ε to real-valued filtering functions defined on R n should allow us to control the expected value of the perturbation of persistence diagrams caused by uniformly distributed impulsive noise, even for n > 1 and in degrees greater than 0. Our long-term goal is to realize Machine Learning methods in which optimal operators for impulse noise removal can be learned and selected in the convex hull of a given finite set of GENEOs, using a theorem which guarantees that the convex combination of GENEOs is still a GENEO (cf. [11]). We plan to devote our research to this topic in the future.

Author Contributions

P.F. devised the project. All authors contributed equally to the manuscript. The authors of this paper have been listed in alphabetical order. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by the European Union in the context of the H2020 EU Marie Curie Initial Training Network project named WAKEUPCALL, and by INdAM-GNSAGA.

Data Availability Statement

Data available on request from the authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Edelsbrunner, H.; Morozov, D. Persistent homology: Theory and practice. In European Congress of Mathematics; European Mathematical Society: Zürich, Germany, 2013; pp. 31–50. [Google Scholar]
  2. Biasotti, S.; Floriani, L.D.; Falcidieno, B.; Frosini, P.; Giorgi, D.; Landi, C.; Papaleo, L.; Spagnuolo, M. Describing shapes by geometrical-topological properties of real functions. ACM Comput. Surv. 2008, 40, 12:1–12:87. [Google Scholar] [CrossRef]
  3. Carlsson, G. Topology and data. Bull. Amer. Math. Soc. 2009, 46, 255–308. [Google Scholar] [CrossRef] [Green Version]
  4. Edelsbrunner, H.; Harer, J. Persistent homology—A survey. In Surveys on Discrete and Computational Geometry; American Mathematical Society: Providence, RI, USA, 2008; Volume 453, pp. 257–282. [Google Scholar]
  5. Cohen-Steiner, D.; Edelsbrunner, H.; Harer, J. Stability of persistence diagrams. Discrete Comput. Geom. 2007, 37, 103–120. [Google Scholar] [CrossRef] [Green Version]
  6. Cohen-Steiner, D.; Edelsbrunner, H.; Harer, J.; Mileyko, Y. Lipschitz functions have Lp-stable persistence. Found. Comput. Math. 2010, 10, 127–139. [Google Scholar] [CrossRef]
  7. Fasy, B.T.; Lecci, F.; Rinaldo, A.; Wasserman, L.; Balakrishnan, S.; Singh, A. Confidence sets for persistence diagrams. Ann. Stat. 2014, 42, 2301–2339. [Google Scholar] [CrossRef]
  8. Buchet, M.; Chazal, F.; Dey, T.K.; Fan, F.; Oudot, S.Y.; Wang, Y. Topological analysis of scalar fields with outliers. In Proceedings of the 31st International Symposium on Computational Geometry (SoCG 2015), Eindhoven, The Netherlands, 22–25 June 2015; Leibniz International Proceedings in Informatics (LIPIcs). Arge, L., Pach, J., Eds.; Schloss Dagstuhl–Leibniz-Zentrum für Informatik: Dagstuhl, Germany, 2015; Volume 34, pp. 827–841. [Google Scholar]
  9. Adler, R.J.; Agami, S. Modelling persistence diagrams with planar point processes, and revealing topology with bagplots. J. Appl. Comput. Topol. 2019, 3, 139–183. [Google Scholar] [CrossRef] [Green Version]
  10. Vishwanath, S.; Fukumizu, K.; Kuriki, S.; Sriperumbudur, B.K. Robust persistence diagrams using reproducing kernels. In Advances in Neural Information Processing Systems; Larochelle, H., Ranzato, M., Hadsell, R., Balcan, M.F., Lin, H., Eds.; Curran Associates, Inc.: New York, NY, USA, 2020; Volume 33, pp. 21900–21911. [Google Scholar]
  11. Bergomi, M.G.; Frosini, P.; Giorgi, D.; Quercioli, N. Towards a topological–geometrical theory of group equivariant non-expansive operators for data analysis and machine learning. Nat. Mach. Intell. 2019, 1, 423–433. [Google Scholar] [CrossRef] [Green Version]
  12. Conti, F.; Frosini, P.; Quercioli, N. On the construction of Group Equivariant Non-Expansive Operators via permutants and symmetric functions. Front. Artif. Intell. 2022, 5, 786091. [Google Scholar] [CrossRef] [PubMed]
  13. Bocchi, G.; Botteghi, S.; Brasini, M.; Frosini, P.; Quercioli, N. On the finite representation of linear group equivariant operators via permutant measures. Ann. Math. Artif. Intell. 2023, 1–23. [Google Scholar] [CrossRef]
  14. Frosini, P.; Jabłoński, G. Combining persistent homology and invariance groups for shape comparison. Discrete Comput. Geom. 2016, 55, 373–409. [Google Scholar] [CrossRef] [Green Version]
  15. Cerri, A.; Fabio, B.D.; Ferri, M.; Frosini, P.; Landi, C. Betti numbers in multidimensional persistent homology are stable functions. Math. Methods Appl. Sci. 2013, 36, 1543–1557. [Google Scholar] [CrossRef] [Green Version]
  16. Cerri, A.; Ethier, M.; Frosini, P. On the geometrical properties of the coherent matching distance in 2D persistent homology. J. Appl. Comput. Topol. 2019, 3, 381–422. [Google Scholar] [CrossRef] [Green Version]
  17. Micheletti, A. A new paradigm for artificial intelligence based on group equivariant non-expansive operators. Eur. Math. Soc. Mag. 2023, 128, 4–12. [Google Scholar] [CrossRef]
  18. Vaseghi, S.V. Impulsive Noise: Modelling, Detection and Removal; John Wiley and Sons, Ltd.: Chichester, UK, 2008; Chapter 13; pp. 341–358. [Google Scholar]
  19. Earnest, M. Average Minimum Distance between n Points Generate i.i.d. with Uniform Dist. Mathematics Stack Exchange. Available online: https://math.stackexchange.com/q/2001026 (accessed on 13 July 2021).
Figure 1. How the sublevel sets change with respect to the filtration induced by φ .
Figure 1. How the sublevel sets change with respect to the filtration induced by φ .
Entropy 25 01150 g001
Figure 2. Persistence diagram of the function 2 sin x on 0 , 3 4 π . The point ( 0 , ) describes the existence of a connected component that is born at zero and never dies. The point ( 2 , 2 ) claims that there is a connected component born at 2 that dies (merges with the other one) at 2. The trivial points on the diagonal u = v are not displayed.
Figure 2. Persistence diagram of the function 2 sin x on 0 , 3 4 π . The point ( 0 , ) describes the existence of a connected component that is born at zero and never dies. The point ( 2 , 2 ) claims that there is a connected component born at 2 that dies (merges with the other one) at 2. The trivial points on the diagonal u = v are not displayed.
Entropy 25 01150 g002
Figure 3. How drastically impulsive noise can influence persistence diagrams.
Figure 3. How drastically impulsive noise can influence persistence diagrams.
Entropy 25 01150 g003
Figure 4. In the figure on the left, we have our original function φ . In the middle we have a noise function R, and in the right figure, we have the corrupted function φ ^ : = φ + R .
Figure 4. In the figure on the left, we have our original function φ . In the middle we have a noise function R, and in the right figure, we have the corrupted function φ ^ : = φ + R .
Entropy 25 01150 g004
Figure 5. Example 1: Comparison between φ and φ ^ .
Figure 5. Example 1: Comparison between φ and φ ^ .
Entropy 25 01150 g005
Figure 6. Example 1: Denoising via convolution with T h ( h = n ).
Figure 6. Example 1: Denoising via convolution with T h ( h = n ).
Entropy 25 01150 g006
Figure 7. Example 1: Denoising via convolution with T h ( h = 1 n ).
Figure 7. Example 1: Denoising via convolution with T h ( h = 1 n ).
Entropy 25 01150 g007
Figure 8. Example 1: Denoising via the proposed GENEO F τ n 2 F τ n .
Figure 8. Example 1: Denoising via the proposed GENEO F τ n 2 F τ n .
Entropy 25 01150 g008
Figure 9. Example 1: Error made by applying a convolution with T n (left) or T 1 n (right).
Figure 9. Example 1: Error made by applying a convolution with T n (left) or T 1 n (right).
Entropy 25 01150 g009
Figure 10. Example 1: Error made by using the GENEO F τ n 2 F τ n .
Figure 10. Example 1: Error made by using the GENEO F τ n 2 F τ n .
Entropy 25 01150 g010
Figure 11. Example 1: Denoising via F σ β F 2 σ β .
Figure 11. Example 1: Denoising via F σ β F 2 σ β .
Entropy 25 01150 g011
Figure 12. Example 1: Histogram counting the number of cases in each of ten bins concerning the overestimation value 3 L σ β F σ β F 2 σ β ( φ ^ ) φ obtained in our simulations. Most of the cases belong to the first bin, specifically 665 cases, containing the functions φ for which 0 3 L σ β F σ β F 2 σ β ( φ ^ ) φ 0.02 .
Figure 12. Example 1: Histogram counting the number of cases in each of ten bins concerning the overestimation value 3 L σ β F σ β F 2 σ β ( φ ^ ) φ obtained in our simulations. Most of the cases belong to the first bin, specifically 665 cases, containing the functions φ for which 0 3 L σ β F σ β F 2 σ β ( φ ^ ) φ 0.02 .
Entropy 25 01150 g012
Figure 13. Example 2: Comparison between φ and φ ^ .
Figure 13. Example 2: Comparison between φ and φ ^ .
Entropy 25 01150 g013
Figure 14. Example 2: Denoising via convolution with T h ( h = n ).
Figure 14. Example 2: Denoising via convolution with T h ( h = n ).
Entropy 25 01150 g014
Figure 15. Example 2: Denoising via convolution with T h ( h = 1 n ).
Figure 15. Example 2: Denoising via convolution with T h ( h = 1 n ).
Entropy 25 01150 g015
Figure 16. Example 2: Denoising via the proposed GENEO F τ n 2 F τ n .
Figure 16. Example 2: Denoising via the proposed GENEO F τ n 2 F τ n .
Entropy 25 01150 g016
Figure 17. Example 2: Error made by applying a convolution with T n (left) or T 1 n (right).
Figure 17. Example 2: Error made by applying a convolution with T n (left) or T 1 n (right).
Entropy 25 01150 g017
Figure 18. Example 2: Error made by using the GENEO F τ n 2 F τ n .
Figure 18. Example 2: Error made by using the GENEO F τ n 2 F τ n .
Entropy 25 01150 g018
Figure 19. Example 2: Denoising via F σ β F 2 σ β .
Figure 19. Example 2: Denoising via F σ β F 2 σ β .
Entropy 25 01150 g019
Figure 20. Example 2: Histogram counting the number of cases in each of ten bins concerning the overestimation value 3 L σ β F σ β F 2 σ β ( φ ^ ) φ obtained in our simulations: 431 cases belong to the first bin and 154 belong to the second one; hence, most of the cases belong to the first two bins containing the functions φ for which 0 3 L σ β F σ β F 2 σ β ( φ ^ ) φ 0.24 .
Figure 20. Example 2: Histogram counting the number of cases in each of ten bins concerning the overestimation value 3 L σ β F σ β F 2 σ β ( φ ^ ) φ obtained in our simulations: 431 cases belong to the first bin and 154 belong to the second one; hence, most of the cases belong to the first two bins containing the functions φ for which 0 3 L σ β F σ β F 2 σ β ( φ ^ ) φ 0.24 .
Entropy 25 01150 g020
Figure 21. Five examples of the functions we produced for N = 3 . In each figure, the functions are displayed without noise (left) and with added noise (right).
Figure 21. Five examples of the functions we produced for N = 3 . In each figure, the functions are displayed without noise (left) and with added noise (right).
Entropy 25 01150 g021
Figure 22. Five examples of the functions we produced for N = 7 . In each figure, the functions are displayed without noise (left) and with added noise (right).
Figure 22. Five examples of the functions we produced for N = 7 . In each figure, the functions are displayed without noise (left) and with added noise (right).
Entropy 25 01150 g022
Figure 23. Plots of the means of the averaged values of φ ^ φ (yellow), the means of 3 L σ β + k α 1 1 8 ( k 1 ) σ β k (brown), and the means of the averaged values of F σ β F 2 σ β ( φ ^ ) φ (blue) for β { 3 , 4 , 5 , , 13 } and L { 1 , 2 , , 10 } , when α varies in the set { 50 , 55 , 60 , , 100 } .
Figure 23. Plots of the means of the averaged values of φ ^ φ (yellow), the means of 3 L σ β + k α 1 1 8 ( k 1 ) σ β k (brown), and the means of the averaged values of F σ β F 2 σ β ( φ ^ ) φ (blue) for β { 3 , 4 , 5 , , 13 } and L { 1 , 2 , , 10 } , when α varies in the set { 50 , 55 , 60 , , 100 } .
Entropy 25 01150 g023
Figure 24. Plots of the means of the averaged values of φ ^ φ (yellow), the means of 3 L σ β + k α 1 1 8 ( k 1 ) σ β k (brown), and the means of the averaged values of F σ β F 2 σ β ( φ ^ ) φ (blue) for α { 50 , 55 , 60 , , 100 } and L { 1 , 2 , , 10 } , when β varies in the set { 3 , 4 , 5 , , 13 } .
Figure 24. Plots of the means of the averaged values of φ ^ φ (yellow), the means of 3 L σ β + k α 1 1 8 ( k 1 ) σ β k (brown), and the means of the averaged values of F σ β F 2 σ β ( φ ^ ) φ (blue) for α { 50 , 55 , 60 , , 100 } and L { 1 , 2 , , 10 } , when β varies in the set { 3 , 4 , 5 , , 13 } .
Entropy 25 01150 g024
Figure 25. Plots of the means of the averaged values of φ ^ φ (yellow), the means of 3 L σ β + k α 1 1 8 ( k 1 ) σ β k (brown), and the means of the averaged values of F σ β F 2 σ β ( φ ^ ) φ (blue) for α { 50 , 55 , 60 , , 100 } and β { 3 , 4 , 5 , , 13 } , when L varies in the set { 1 , 2 , , 10 } .
Figure 25. Plots of the means of the averaged values of φ ^ φ (yellow), the means of 3 L σ β + k α 1 1 8 ( k 1 ) σ β k (brown), and the means of the averaged values of F σ β F 2 σ β ( φ ^ ) φ (blue) for α { 50 , 55 , 60 , , 100 } and β { 3 , 4 , 5 , , 13 } , when L varies in the set { 1 , 2 , , 10 } .
Entropy 25 01150 g025
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Frosini, P.; Gridelli, I.; Pascucci, A. A Probabilistic Result on Impulsive Noise Reduction in Topological Data Analysis through Group Equivariant Non-Expansive Operators. Entropy 2023, 25, 1150. https://doi.org/10.3390/e25081150

AMA Style

Frosini P, Gridelli I, Pascucci A. A Probabilistic Result on Impulsive Noise Reduction in Topological Data Analysis through Group Equivariant Non-Expansive Operators. Entropy. 2023; 25(8):1150. https://doi.org/10.3390/e25081150

Chicago/Turabian Style

Frosini, Patrizio, Ivan Gridelli, and Andrea Pascucci. 2023. "A Probabilistic Result on Impulsive Noise Reduction in Topological Data Analysis through Group Equivariant Non-Expansive Operators" Entropy 25, no. 8: 1150. https://doi.org/10.3390/e25081150

APA Style

Frosini, P., Gridelli, I., & Pascucci, A. (2023). A Probabilistic Result on Impulsive Noise Reduction in Topological Data Analysis through Group Equivariant Non-Expansive Operators. Entropy, 25(8), 1150. https://doi.org/10.3390/e25081150

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop