*Article* **The Riemann-Lebesgue Integral of Interval-Valued Multifunctions**

#### **Danilo Costarelli 1, Anca Croitoru 2, Alina Gavrilu¸t <sup>2</sup> and Alina Iosif <sup>3</sup> and Anna Rita Sambucini 1,\***


Received: 10 November 2020; Accepted: 16 December 2020; Published: 20 December 2020

**Abstract:** We study Riemann-Lebesgue integrability for interval-valued multifunctions relative to an interval-valued set multifunction. Some classic properties of the *RL* integral, such as monotonicity, order continuity, bounded variation, convergence are obtained. An application of interval-valued multifunctions to image processing is given for the purpose of illustration; an example is given in case of fractal image coding for image compression, and for edge detection algorithm. In these contexts, the image modelization as an interval valued multifunction is crucial since allows to take into account the presence of quantization errors (such as the so-called round-off error) in the discretization process of a real world analogue visual signal into a digital discrete one.

**Keywords:** Riemann-Lebesgue integral; interval valued (set) multifunction; non-additive set function; image processing

**MSC:** 28B20; 28C15; 49J53

#### **1. Introduction**

The theory of multifunctions is an important field of research. Since interval arithmetic, introduced by Moore in [1], it appears a natural option for handling the uncertainty in data and in sensor measurements, particular attention was addressed to the study of interval-valued multifunctions and multimeasures because of their applications in statistics, biology, theory of games, economics, social sciences and software, to keep track of rounding errors in calculations and of uncertainties in the knowledge of the exact values of physical and technical parameters (see for example [2–5]). In fact, since the uncertainty of information could affect an expert's opinion, the ability to consider the uncertainty information during the process could be very important, see for example [2–4,6–11] and the references therein.

However, in some recent papers, interval-valued multifunctions have been applied also to some new directions, involving signal and image processing. Digital images are in fact the result of a discretization of the reality; namely sampled version of a continuous signal. Hence, there are different sources of uncertainty and ambiguity to be considered when performing image processing tasks, see for example [12,13]. For instance, the applications of fractal image coding for image compression [14,15] is one of the topic in which interval-valued multifunctions have been applied. Clearly, image compression techniques [16] are very useful in order to speed up the processes of digital image transmission and to improve the efficiency of image storage for high dimensional databases [17]. Further, applications of interval-valued multifunctions to the implementation of edge detection algorithms can also be found (see e.g., [13,18]).

In the literature several methods of integration for functions and multifunctions have been studied extending the Riemann and Lebesgue integrals. In this framework a generalization of Riemann sums was given in [19–37] while another generalization is due to Kadets and Tseytlin [38], who introduced the absolute Riemann-Lebesgue |*RL*| and unconditional Riemann-Lebesgue *RL* integrability, for Banach valued functions with respect to countably additive measures. They proved that in finite measure space, the Bochner integrability implies |*RL*| integrability which is stronger than *RL* integrability that implies Pettis integrability. Regarding this last extension contributions are given also in [21,23,34,39].

In the last decade the study of non-additive set functions and multifunctions has recently received a wide recognition, (see also [3,9,10,40–46]). In this paper, motivated by the large number of fields in which the interval-valued multifunction can be applied, we introduce a new type of integral of an interval-valued multifunction *G* with respect to an interval-valued submeasure *M* with respect to the weak interval order relation introduced in [4] by Guo and Zhang. Although the construction procedure of the integral is similar to the one given in [34,38,39], the integral proposed is a generalization of it since we are concerned with the study of a Riemann-Lebesgue set-valued integrand with respect to an arbitrary interval-valued set function, not necessarily countably additive. So the novelty of this construction concerns not only the codomain of the integrands but also the non-additivity of the measure with respect to which they are integrated. The main results on this subject are Theorem 1, in which the additivity of the integral is proved even if the pair (*G*, *M*) does not satisfy this property; the monotonicity and the order continuity are established in Theorems 2 and 4 and a convergent result given in Theorem 5.

The paper is organized as follows: in Section 2 the basic concepts and terminology are introduced together with some remarks. In Section 3 we introduce the RL-integral of an interval-valued multifunction with respect to an interval valued subadditive multifunction and we provide a comprehensive treatment of the integration theory together with a comparison with other integrals defined in the same setting (Remark 8). An example of an application in image processing is given in Section 3.1. The applications concerning image processing discussed in the present paper is given for the purpose of illustration and is new. The main reason for which we discuss the above application is to provide examples and justifications of the uses of interval-valued multifunctions to concrete applications in Image Processing. The advantage of using the notion of interval-valued multifunction in signal analysis is that this formalism allows to include in a unique framework possible uncertainty or the noise on the evaluation of an image at any given pixel.

#### **2. Preliminaries**

Let *S* be a nonempty at least countable set, P(*S*) the family of all subsets of *S* and A a *σ*-algebra of subsets of *S*. The symbol R<sup>+</sup> <sup>0</sup> denotes, as usual, the set of non negative real numbers.

#### **Definition 1** ([34], Definition 2.1)**.**


We denote by P the class of all the countable partitions of *S* and if *A* ∈ A is fixed, by P<sup>A</sup> we denote the class of all the countable partitions of the set *A*.

**Definition 2** ([34], Definition 2.2)**.** *Let m* : A → [0, +∞) *be a non-negative function, with m*(∅) = 0. *A set A* ∈ A *is said to be an atom of m if m*(*A*) > 0 *and for every B* ∈ A*, with B* ⊂ *A, it is m*(*B*) = 0 *or m*(*A*\*B*) = 0*.*

*m is said to be:*


*Moreover m satisfies property* (*σ*) *if the ideal of m-zero sets is stable under countable unions (see for example [34], Definition 2.3).*

We denote by the symbol *ck*(R) the family of all non-empty convex compact subsets of R, by convention, {0} = [0, 0]. We consider on *ck*(R) the Minkowski addition (*A* + *B* : = {*a* + *b* : *a* ∈ *A*, *b* ∈ *B*}) and the standard multiplication by scalars. *A* := sup{|*x*|: *x* ∈ *A*}. *dH* is the Hausdorff distance in *ck*(R), while *e*(*A*, *B*) = sup{*d*(*x*, *B*), *x* ∈ *A*} and *dH*(*A*, *B*) = max{*e*(*A*, *B*), *e*(*B*, *A*)}.

(*ck*(R), *dH*) is a complete metric space ([48,49]), but is not a linear space since the subtraction is not well defined.

If *A* = [*a*, *b*] then *A* = max{|*a*|, |*b*|}. Moreover

$$\begin{aligned} d\_H([a,b],[c,d]) &= \max\{|a-c|, |b-d|\}, \qquad \forall \ a, b, c, d \in \mathbb{R} \\ d\_H([0,a],[0,b]) &= |b-a| \qquad \forall \ a, b \in \mathbb{R}\_0^+. \end{aligned}$$

In the family *ck*(R) the following operations are also considered, for every *a*, *b*, *c*, *d* ∈ R:

(**i**) [*a*, *b*] *·* [*c*, *d*]=[*ac*, *bd*];

(**ii**) [*a*, *b*] ⊆ [*c*, *d*] if and only if *c* ≤ *a* ≤ *b* ≤ *d*;

(**iii**) [*a*, *b*] ! [*c*, *d*] if and only if *a* ≤ *c* and *b* ≤ *d*; (weak interval order)

(**iv**) [*a*, *b*] ∧ [*c*, *d*]=[min{*a*, *c*}, min{*b*, *d*}];

(**v**) [*a*, *b*] ∨ [*c*, *d*]=[max{*a*, *c*}, max{*b*, *d*}].

In general there is no relation between ! (**iii**) and ⊆ (**ii**); they only coincide on the subfamily {[0, *<sup>a</sup>*], *<sup>a</sup>* <sup>≥</sup> <sup>0</sup>}. Let *ck*(R<sup>+</sup> <sup>0</sup> ) := {[*a*, *b*], *a*, *b* ∈ R and 0 ≤ *a* ≤ *b*}.

In this paper we consider (*ck*(R<sup>+</sup> <sup>0</sup> ), *dH*, !), namely the space *ck*(R<sup>+</sup> <sup>0</sup> ) is endowed with the Hausdorff distance and the weak interval order. As a particular case of [20] (Definition 2.1) we have:

**Definition 3.** *Let* (*an*)*n*,(*bn*)*<sup>n</sup> be two sequences of real numbers so that* 0 ≤ *an* ≤ *bn*, ∀*n* ∈ N*.*

*The series* ∑<sup>∞</sup> *<sup>n</sup>*=0[*an*, *bn*] :<sup>=</sup> {∑<sup>∞</sup> *<sup>n</sup>*=<sup>0</sup> *yn* : *an* ≤ *yn* ≤ *bn*, ∀*n* ∈ N } *is called convergent if the sequence of partial sums Sn* := [∑*<sup>n</sup> <sup>k</sup>*=<sup>0</sup> *ak*, <sup>∑</sup>*<sup>n</sup> <sup>k</sup>*=<sup>0</sup> *bk*] *is dH-convergent to it.*

**Remark 1.** *It is easy to see that* ∑<sup>∞</sup> *<sup>n</sup>*=0[*an*, *bn*]=[*u*, *<sup>v</sup>*]*, with* <sup>0</sup> <sup>≤</sup> *<sup>u</sup>* <sup>≤</sup> *<sup>v</sup>* <sup>&</sup>lt; <sup>∞</sup>*, if and only if* <sup>∑</sup><sup>∞</sup> *<sup>n</sup>*=<sup>0</sup> *an* = *u and* ∑<sup>∞</sup> *<sup>n</sup>*=<sup>0</sup> *bn* = *v*.

We recall the following definition for the integrable Banach-valued functions *f* : *S* → *X* with respect to non-negative measures given in [38,39]:

**Definition 4.** *A function f is called unconditional Riemann-Lebesgue (RL ) m-integrable (on S) if there exists b* ∈ *X such that for every ε* > 0*, there exists a countable partition P<sup>ε</sup> of S, so that for every countable partition <sup>P</sup>* = {*An*}*n*∈<sup>N</sup> *of <sup>S</sup> with <sup>P</sup>* ≥ *<sup>P</sup>ε, <sup>f</sup> is bounded on every An, with <sup>m</sup>*(*An*) > 0 *and for every tn* ∈ *An, <sup>n</sup>* ∈ N*, the series* ∑+<sup>∞</sup> *<sup>n</sup>*=<sup>0</sup> *f*(*tn*)*m*(*An*) *is unconditional convergent and*

$$\left\| \sum\_{n=0}^{+\infty} f(t\_n) m(A\_n) - b \right\| < \varepsilon.$$

The vector *b* (necessarily unique) is called *the Riemann-Lebesgue m-integral of f on S* and it is denoted by (*RL*) *S f dm*. The *RL* definition of the integrability on a subset *A* ∈ A is given in the classical manner.

**Remark 2.** *We remember that, in the countably additive case, unconditional RL-integrability is stronger than Birkhoff integrability (in the sense of Fremlin), see Ref. [23] and the references therein; while the notion of unconditional Riemann-Lebesgue integrability coincides with Birkhoff's one given in [21] (Definition 1, Proposition 2.6 and note at p. 8).*

For the properties of this integral with respect to a submeasure we refer to the results given in [34]. Moreover we have that

**Proposition 1.** *Let gn* : *<sup>S</sup>* <sup>→</sup> <sup>R</sup><sup>+</sup> <sup>0</sup> *be an increasing sequence of bounded RL integrable function with respect to a submeasure <sup>μ</sup>* : A → <sup>R</sup><sup>+</sup> <sup>0</sup> *of bounded variation. If there exists a g* : *<sup>S</sup>* <sup>→</sup> <sup>R</sup><sup>+</sup> <sup>0</sup> *such that*


*then g is RL integrable with respect to μ and*

$$\lim\_{n \to \infty} \left( \mathbb{R} \mathbb{L} \right) \int\_{S} \mathbb{g}\_{n} \, d\mu = (\mathbb{R} \mathbb{L}) \int\_{S} \mathbb{g} \, d\mu.$$

**Proof.** Since *gn* ↑, by the monotonicity we have that (*RL*) *S gndμ* ↑ so sup*<sup>n</sup>* (*RL*) *S gn dμ* = lim*n*→<sup>∞</sup> (*RL*) *S gn <sup>d</sup><sup>μ</sup>* <sup>=</sup> *<sup>u</sup>* <sup>∈</sup> <sup>R</sup><sup>+</sup> <sup>0</sup> . Thanks to uniform convergence *g* is bounded; let *L* > 0 an upper bound for *g*.

Let *ε* > 0 be fixed and consider *k*(*ε*) ∈ N be such that

$$\begin{aligned} \left| \mathcal{g}(t) - \mathcal{g}\_{k(\varepsilon)}(t) \right| &< \frac{\varepsilon}{3\mu(S)} \qquad \forall \ t \in S, \quad \text{and} \\\left| \,^{(RL)} \int\_S \mathcal{g}\_{k(\varepsilon)} d\mu - \mu \right| &< \frac{\varepsilon}{3}. \end{aligned}$$

For every countable partition *P* := (*An*)*<sup>n</sup>* finer than *Pε*/3,*k*(*ε*) (the one that verifies Definition 4 for *gk*(*ε*)) and for every *tn* <sup>∈</sup> *An* we have that <sup>∑</sup>+<sup>∞</sup> *<sup>n</sup>*=<sup>0</sup> *g*(*tn*)*μ*(*An*) converges, since *μ* is of bounded variation.

In fact *<sup>g</sup>*(*tn*)*μ*(*An*) <sup>≤</sup> *<sup>L</sup>μ*(*An*) for every *<sup>n</sup>* <sup>∈</sup> <sup>N</sup> and, for every *<sup>k</sup>* <sup>∈</sup> <sup>N</sup>, it is 0 <sup>≤</sup> <sup>∑</sup>*<sup>k</sup> <sup>n</sup>*=<sup>0</sup> *μ*(*An*) ≤ *μ*(*S*). Moreover

$$\begin{split} \left| \left| \sum\_{n=0}^{+\infty} \mathcal{g}(t\_{\mathbb{R}}) \mu(A\_{\mathbb{R}}) - u \right| \right| &\leq \left| \sum\_{n=0}^{+\infty} \mathcal{g}(t\_{\mathbb{R}}) \mu(A\_{\mathbb{R}}) - \sum\_{n=0}^{+\infty} \mathcal{g}\_{k(\varepsilon)}(t\_{\mathbb{R}}) \mu(A\_{\mathbb{R}}) \right| + \left| \sum\_{n=0}^{+\infty} \mathcal{g}\_{k(\varepsilon)}(t\_{\mathbb{R}}) \mu(A\_{\mathbb{R}}) - \sum\_{n=0}^{+\infty} \mathcal{g}\_{k(\varepsilon)}(t\_{\mathbb{R}}) \mu(A\_{\mathbb{R}}) \right| \\ &+ \left| \left| \sum\_{n=0}^{+\infty} \mathcal{g}\_{k(\varepsilon)} \mu(A\_{\mathbb{R}}) \mu(A\_{\mathbb{R}}) \right| \leq \varepsilon. \end{split}$$

**Remark 3.** *We can extend Proposition 1 to the bounded sequences* (*gn*)*<sup>n</sup> that converge μ-almost uniformly on S (namely to the sequences* (*gn*)*<sup>n</sup> such that for every ε* > 0 *there exists B*(*ε*) ∈ A *with μ*(*B*(*ε*)) ≤ *ε and gn converges uniformly to g on S* \ *B*(*ε*)*), if we assume that even g is bounded.*

*We can proceed in fact in the same way, as in the previous proof, taking P*∗ *<sup>ε</sup>* := *Pε*/3,*k*(*ε*) ∧ {*S* \ *B*(*ε*), *B*(*ε*)} *and, for every countable partition P* := (*An*)*<sup>n</sup> finer than P*<sup>∗</sup> *<sup>ε</sup> , dividing* <sup>∑</sup>+<sup>∞</sup> *<sup>n</sup>*=<sup>0</sup> *g*(*tn*)*μ*(*An*) *in two parts: the one relative to S* \ *B*(*ε*)*, where the uniform convergence is assumed, and the remining part.*

*Convergence results in Gould integrability of functions with respect to a submeasure of finite variation are established for instance in [50].*

Given two submeasures *<sup>μ</sup>*1, *<sup>μ</sup>*<sup>2</sup> : A → <sup>R</sup><sup>+</sup> <sup>0</sup> with *μ*1(*A*) ≤ *μ*2(*A*) for every *A* ∈ A let *M* : A → *ck*(R<sup>+</sup> <sup>0</sup> ) defined by

$$M(A) = [\mu\_1(A), \mu\_2(A)].\tag{1}$$

*M* is called an *interval submeasure*. For results in this subject see for example [3,43].

Let *<sup>M</sup>* : A → *ck*(R<sup>+</sup> <sup>0</sup> ). We say that *M* is an *interval valued multisubmeasure* if


In literature the multimeasures that satisfy the first two statements are also called set valued fuzzy measures (see for example [4] (Definition 1), [3,11,42–44] and the references therein).

A very interesting case of interval-valued multisubmeasure was given, for the first time, in [6,8] where Dempster and Shefer proposed a mathematical theory of evidence using non additive measures: Belief and Plausibility in such a way for every set *A* the *Belief interval* of the set is [*Bel*(*A*), *Pl*(*A*)]. This theory is capable of deriving probabilities for a collection of hypotheses and it allows the system inferencing with the imprecision and uncertainty. If the target space is *ck*([0, 1]) it is used for example in decision theory.

We say that *M* is an *additive multimeasure* if *M*(*A* ∪ *B*) = *M*(*A*) + *M*(*B*) for every disjoint sets *A*, *B* ∈ A.

If a multimeasure *M* is countably additive in the Hausdorff metric *dH*, then it is called a *dH*-*multimeasure*. In this case we have that lim*n*→<sup>∞</sup> *dH* (∑*<sup>n</sup> <sup>k</sup>*=<sup>1</sup> *M*(*Ak*), *M*(*A*)) = 0, for every sequence of pairwice disjoint sets (*An*)*<sup>n</sup>* ⊂ A such that ∪*nAn* = *A*.

**Remark 4.** *By Ref. [43] (Remark 3.6) M*(*A*)=[*μ*1(*A*), *μ*2(*A*)] *is a multisubmeasure with respect to* ! *if and only if μ*1, *μ*<sup>2</sup> *are submeasures in the sense of Definition 2 (iii). Moreover M is monotone, finitely additive, order-continuous, exhaustive respectively if and only if the set functions μ*<sup>1</sup> *and μ*<sup>2</sup> *are the same (see [40] (Proposition 2.5, Remark 3.3)).*

**Definition 5.** *Let M* : A → *ck*(R<sup>+</sup> <sup>0</sup> )*. The variation of M is the set function M* : P(*S*) → [0, +∞] *defined by*

$$\overline{M}(E) = \sup \{ \sum\_{i=1}^{n} ||M(A\_i)||\_\prime \: \mid \{A\_i\}\_{i=1}^n \subset \mathcal{A}\_\prime \: \: A\_i \subseteq E, A\_i \cap A\_j = \mathcal{O}, \, i \neq j \} \dots$$

*M is said to be of finite variation if M*(*S*) < ∞*.*

**Remark 5.** *We can observe that if E* ∈ A*, then in the definition of M one may consider the supremum over all finite partitions* {*Ai*}*<sup>n</sup> <sup>i</sup>*=1∈ P*E. If M is finitely additive, then M*(*A*) = *M*(*A*), *for every A* ∈ A.

*If M is subadditive (countably subadditive, respectively) of finite variation, then M is finitely additive (countably additive, respectively). Finally, if M*(*A*)=[*μ*1(*A*), *μ*2(*A*)]*, for every A* ∈ A*, then M* = *μ*2*.*

#### **3.** *RL* **Interval Valued Integral and Its Properties**

In this section, we introduce and study Riemnn-Lebesgue integrability of interval-valued multifunctions with respect to interval-valued set multifunctions, pointing out various properties of this integral. For this, unless stated otherwise, in what follows suppose *S* is a nonempty set, with card *S* ≥ ℵ<sup>0</sup> (card *S* is the cardinality of *S*), A is a *σ*-algebra of subsets of *S*.

The multisubmeasure *M* here considered is an interval-valued one and satisfies (1).

Given *<sup>g</sup>*1, *<sup>g</sup>*<sup>2</sup> : *<sup>S</sup>* <sup>→</sup> <sup>R</sup><sup>+</sup> <sup>0</sup> with *<sup>g</sup>*1(*s*) <sup>≤</sup> *<sup>g</sup>*2(*s*) for all *<sup>s</sup>* <sup>∈</sup> *<sup>S</sup>*, let *<sup>G</sup>* : *<sup>S</sup>* <sup>→</sup> *ck*(R<sup>+</sup> <sup>0</sup> ) be the interval-valued multifunction defined by *G*(*s*)=[*g*1(*s*), *g*2(*s*)] for every *s* ∈ *S*. For every countable tagged partition Π := {(*Bn*,*sn*), *n* ∈ N} of *S* we denote by

$$\begin{aligned} \sigma\_{G,M}(\Pi) &:= \sum\_{n=1}^{\infty} G(s\_n) \cdot M(B\_n) = \sum\_{n=1}^{\infty} [\lg\_1(s\_n)\mu\_1(B\_n), \lg\_2(s\_n)\mu\_2(B\_n)] = 1, \\\ &= \{ \sum\_{n=1}^{\infty} y\_n, \ y\_n \in [\lg\_1(s\_n)\mu\_1(B\_n), \lg\_2(s\_n)\mu\_2(B\_n)], \ n \in \mathbb{N} \}. \end{aligned}$$

By [20] (Lemma 2.2) the set *σG*,*M*(Π) is closed and convex in R<sup>+</sup> <sup>0</sup> , so it is an interval [*u*(Π) *<sup>G</sup>*,*M*, *v* (Π) *<sup>G</sup>*,*M*].

**Definition 6.** *A multifunction <sup>G</sup>* : *<sup>S</sup>* <sup>→</sup> *ck*(R<sup>+</sup> <sup>0</sup> ) *is called Riemann-Lebesgue RL integrable with respect to M (on S) if there exists* [*a*, *<sup>b</sup>*] <sup>∈</sup> *ck*(R<sup>+</sup> <sup>0</sup> ) *such that for every ε* > 0*, there exists a countable partition P<sup>ε</sup> of S, so that for every tagged partition P* = {(*An*, *tn*)}*n*∈<sup>N</sup> *of S with P* ≥ *<sup>P</sup>ε, the series <sup>σ</sup>G*,*M*(*P*) *is convergent and*

$$d\_H(\sigma\_{G,M}(P), [a, b]) < \varepsilon. \tag{2}$$

[*a*, *b*] *is called the Riemann-Lebesgue integral of G with respect to M and it is denoted*

$$\left[a,b\right] = \left(\mathbb{R}\mathbb{L}\right) \int\_{S} G \, dM.$$

*Obviously, if it exists, is unique.*

**Example 1.** *Suppose <sup>S</sup>* <sup>=</sup> {*sn*|*<sup>n</sup>* <sup>∈</sup> <sup>N</sup>} *is countable,* {*sn*}∈A*, for every <sup>n</sup>* <sup>∈</sup> <sup>N</sup>*, and let <sup>G</sup>* : *<sup>S</sup>* <sup>→</sup> *ck*(R<sup>+</sup> <sup>0</sup> ) *be such that the series* ∞ ∑ *n*=0 *gi*(*sn*)*μi*({*tn*}), *i* = 1, 2 *are convergent. Then G is RL integrable with respect to M and*

$$\mathcal{L}(RL)\int\_{S} G \, dM = \left[\sum\_{n=0}^{\infty} \mathcal{g}\_1(s\_n)\mu\_1(\{s\_n\}), \sum\_{n=0}^{\infty} \mathcal{g}\_2(s\_n)\mu\_2(\{s\_n\})\right].$$

*Observe moreover that, in this case, the RL-integrability of such G with respect to M implies that the product G · G, as defined in* **i***), is integrable in the same sense. In particular, if such G is a discrete or countable interval-valued signal, the* (*RL*) *S G · G dM represents the energy of the signal.*

If *<sup>M</sup>* is of bounded variation and *<sup>G</sup>* : *<sup>S</sup>* <sup>→</sup> *ck*(R<sup>+</sup> <sup>0</sup> ) is bounded and such that *G* = {0} *M*-a.e., then, by [34] (Theorem 3.4), *G* is *M*-integrable and (*RL*) *GdM* = {0}.

*S* From now on we suppose that *G* is bounded and *μ*<sup>2</sup> is of finite variation.

**Proposition 2.** *An interval multifunction G* = [*g*1, *g*2] *is RL integrable with respect to M on S if and only if gi are RL integrable with respect to μi, i* = 1, 2 *and*

$$\int\_{S} \mathbb{G} dM = \left[ (\mathbb{R} \mathbb{L}) \int\_{S} \mathcal{g}\_1 d\mu\_{1\prime} \left( \mathbb{R} \mathbb{L} \right) \int\_{S} \mathcal{g}\_2 d\mu\_2 \right]. \tag{3}$$

**Proof.** Suppose that *G* = [*g*1, *g*2] is *RL* integrable with respect to *M* = [*μ*1, *μ*2], that means there exists [*a*, *<sup>b</sup>*] <sup>∈</sup> *ck*(R<sup>+</sup> <sup>0</sup> ) such that for every *ε* > 0, there exists a countable partition *P<sup>ε</sup>* of *S*, so that for every tagged partition *<sup>P</sup>* = {(*An*, *tn*)}*n*∈<sup>N</sup> of *<sup>S</sup>* with *<sup>P</sup>* ≥ *<sup>P</sup>ε*, the series *<sup>σ</sup>G*,*M*(*P*) is convergent and

$$d\_H([u\_{G,M'}^{(P)}, v\_{G,M}^{(P)}]\_\prime, [a, b]) := \max\{|u\_{G,M}^{(P)} - a|\_\prime, |v\_{G,M}^{(P)} - b|\} < \varepsilon.$$

By this inequality it follows that

$$\max \{ |\sum\_{n=1}^{\infty} g\_1(t\_n)\mu\_1(A\_{\boldsymbol{n}}) - a|\_{\boldsymbol{\prime}} \mid \sum\_{n=1}^{\infty} g\_2(t\_n)\mu\_2(A\_{\boldsymbol{n}}) - b|\_{\boldsymbol{\cdot}} \} \le \varepsilon, \qquad \forall \boldsymbol{n} \in \mathbb{N}\_{\boldsymbol{\cdot}}$$

for every tagged partition *<sup>P</sup>* = {(*An*, *tn*)}*n*∈<sup>N</sup> of *<sup>S</sup>* with *<sup>P</sup>* ≥ *<sup>P</sup><sup>ε</sup>* and then *gi* are *RL* integrable with respect to *μi*, *i* = 1, 2. Formula (3) follows from the convexity of the *RL* integral.

For the converse, for every *ε* > 0, let *Pε*,*gi* , *i* = 1, 2 two countable partitions that verify the definition of *RL* integrability for *gi*, *i* = 1, 2. Let *P<sup>ε</sup>* be a countable partition of *S* with *P<sup>ε</sup>* ≥ *Pε*,*g*<sup>1</sup> ∧ *Pε*,*g*<sup>2</sup> . Then, for every *P* := {*Bn*, *n* ∈ N} ≥ *P<sup>ε</sup>* and for every *tn* ∈ *Bn* it is

$$\left| \sum\_{n=0}^{+\infty} \lg\_i(t\_n) \mu\_i(B\_n) - (\text{RL}) \int\_{\mathcal{S}} \lg\_i d\mu\_i \right| < \varepsilon, \quad i = 1, 2.$$

Since *gi*, *i* = 1, 2 are selections of *G* this means that

$$d\_H\left(\left[\mu\_{G,M'}^{(P)}\sigma\_{G,M}^{(P)}\right]\_\prime \left[\left(\mathbb{R}\mathbb{L}\right)\int\_S \mathcal{G}\_1 d\mu\_{1\prime\prime}\left(\mathbb{R}\mathbb{L}\right)\int\_S \mathcal{G}\_2 d\mu\_2\right]\right) \le \varepsilon^2$$

and then the assertion follows.

**Remark 6.** *By Definition 6 and Proposition 2 we obtain the following definitions for the following cases:*

• *If <sup>M</sup>* <sup>=</sup> {*μ*} : A → <sup>R</sup><sup>+</sup> <sup>0</sup> *is an arbitrary set function and G* = [*g*1, *g*2] *with g*1(*s*) ≤ *g*2(*s*) *for every s* ∈ *S then*

$$\int\_{S} \mathcal{G}dM = \left[ (\mathbb{R}L) \int\_{S} \mathcal{G}\_1 d\mu\_{\prime\prime}(\mathbb{R}L) \int\_{S} \mathcal{G}\_2 d\mu \right].$$

• *If M* = [*μ*1, *<sup>μ</sup>*2] *as in (1) and G* <sup>=</sup> {*g*} : *<sup>S</sup>* <sup>→</sup> <sup>R</sup><sup>+</sup> <sup>0</sup> *then*

$$\int\_{S} \mathcal{G}dM = \left[ (\mathbb{R} \mathbb{L}) \int\_{S} \mathcal{G}d\mu\_{1\prime\prime}(\mathbb{R}\mathbb{L}) \int\_{S} \mathcal{G}d\mu\_{2} \right].$$

**Proposition 3.** *Let G be an interval valued multifuncion. The RL integrability with respect to M is hereditary on subsets A* ∈ A *. Moreover G is RL integrable with respect to M on A if and only if Gχ<sup>A</sup> (where χ<sup>A</sup> is the characteristic function of the set A) is RL integrable with respect to M on S. In this case, for every A* ∈ A*,*

$$\mathcal{S}^{(\mathbb{R}L)} \int\_{A} G \, dM = \left. \begin{pmatrix} \mathbb{R}L \\ \end{pmatrix} \right|\_{S} G \chi\_{A} \, dM.$$

**Proof.** Assume that *G* is *RL* integrable in *S* with respect to *M*. Let *A* ∈ A and denote by [*a*, *b*] the integral of *G*; then, for every *ε* > 0, there exists a countable partition *P<sup>ε</sup>* of *S*, such that, for every finer countable partition *<sup>P</sup>* := {*An*}*n*∈<sup>N</sup> and for every *tn* ∈ *An* it is

$$d\_H\left(\sigma\_{G,M}(P')\_\prime, [a, b]\right) \le \varepsilon.$$

Let *P*<sup>0</sup> be a partition such that *P*<sup>0</sup> ≥ *P<sup>ε</sup>* ∧ {*A*, *T* \ *A*}, and we denote by *PA* ⊂ *P*<sup>0</sup> the corresponding partition of the set *A*. Let Π*<sup>A</sup>* be a partition of *A* finer than *PA*, and extend it with a common partition of *S* \ *A* in such a way the new partition is finer than *Pε*.

It is possible to prove that *σG*,*M*(Π*A*) satisfy a Cauchy principle in *ck*(R<sup>+</sup> <sup>0</sup> ), and so the first claim follows by the completeness of the space. The equality follows from [34] (Theorem 3.2) and Proposition 2.

**Remark 7.** *It is easy to see that, if G is RL integrable with respect to M, for every α* ≥ 0 *it is:*

*(a) αG is RL integrable with respect to M and* (*RL*) *S αG dM* = *α*(*RL*) *S G dM*. *(b) G is RL integrable with respect to αM and* (*RL*) *S G d*(*αM*) = *α*(*RL*) *S G dM*.

**Theorem 1.** *If <sup>G</sup> is an interval valued RL integrable with respect to <sup>M</sup> multifunction, then IG* : A → *ck*(R<sup>+</sup> 0 ) *defined by*

$$I\_G(A) := \left(\mathbb{R}\mathbb{L}\right) \int\_A G \, dM$$

*is a finitely additive multimeasure.*

**Proof.** By Proposition <sup>3</sup> we have that *IG*(*A*) <sup>∈</sup> *ck*(R<sup>+</sup> <sup>0</sup> ) for every *A* ∈ A. In order to prove the additivity we can observe that, for every *A*, *B* ∈ A with *A* ∩ *B* = ∅

$$I\_G(A \cup B) \quad = \quad (\mathbb{R}\mathbb{L}) \int\_S \mathbb{G} \chi\_{A \cup B} \, dM = (\mathbb{R}\mathbb{L}) \int\_S (\mathbb{G} \chi\_A + \mathbb{G} \chi\_B) \, dM. \tag{4}$$

If we prove that for every pair of interval valued *RL* integrable with respect to *M* multifunctions *G*1, *G*<sup>2</sup> we have that

$$\mathcal{L}(\mathbb{R}\mathcal{L}) \int\_{\mathcal{S}} (\mathcal{G}\_1 + \mathcal{G}\_2) \, dM = (\mathbb{R}\mathcal{L}) \int\_{\mathcal{S}} \mathcal{G}\_1 \, dM + (\mathbb{R}\mathcal{L}) \int\_{\mathcal{S}} \mathcal{G}\_2 \, dM \tag{5}$$

the assertion follows. In order to prove formula (5) let *ε* > 0 be fixed. Since *G*1, *G*<sup>2</sup> are *RL* integrable with respect to *M*, for every *ε* > 0 there exists a countable partition *P<sup>ε</sup>* ∈ P such that for every *<sup>P</sup>* = {*An*}*n*∈<sup>N</sup> ≥ *<sup>P</sup><sup>ε</sup>* and every *tn* ∈ *An*, *<sup>n</sup>* ∈ N, the series *<sup>σ</sup>Gi*,*M*(*P*), *<sup>i</sup>* = 1, 2 are convergent and

$$d\_H\left(\sigma\_{G\_i,M}(P)\_{\prime}(\mathbb{R}L)\int\_S G\_i \,dM\right) < \frac{\varepsilon}{2}, \qquad i = 1,2.$$

Then *σG*1+*G*2,*M*(*P*) is convergent and, by [48] (Proposition 1.17),

$$d\_H \left( \sigma\_{G\_1 + G\_2, M}(P)\_{\prime \prime}(\mathbb{R}L) \int\_S G\_1 \, dM + (\mathbb{R}L) \int\_S G\_2 \, dM \right) < \varepsilon. $$

So *G*<sup>1</sup> + *G*<sup>2</sup> is *RL* integrable with respect to *M* and formula (5) is satisfied.

Now applying formula (5) with *G*<sup>1</sup> = *GχA*, *G*<sup>2</sup> = *Gχ<sup>B</sup>* to formula (4) we obtain the additivity of *IG*.

The set-valued integral is monotone relative to the order relation "!" and the inclusion one, with respect to the interval-valued integrands.

**Proposition 4.** *If F*, *G are two RL integrable with respect to M interval valued multifunctions with F* ! *G then, for every A* ∈ A*, IF*(*A*) ! *IG*(*A*)*.*

**Proof.** We will prove for *A* = *S*. Let *F*(*s*) := [ *f*1(*s*), *f*2(*s*)], *G*(*s*)=[*g*1(*s*), *g*2(*s*)]. By the integrability of *F* and *G* we have, by Proposition 2

$$\begin{aligned} I\_F(S) &:= (\mathbb{R}\mathbb{L}) \int\_S F \, d\mathcal{M} = \left[ (\mathbb{R}\mathbb{L}) \int\_S f\_1 d\mu\_{1\prime\prime}(\mathbb{R}\mathbb{L}) \int\_S f\_2 d\mu\_2 \right], \\ I\_G(S) &:= (\mathbb{R}\mathbb{L}) \int\_S G \, d\mathcal{M} = \left[ (\mathbb{R}\mathbb{L}) \int\_S \mathcal{G}\_1 d\mu\_{1\prime\prime}(\mathbb{R}\mathbb{L}) \int\_S \mathcal{G}\_2 d\mu\_2 \right]. \end{aligned}$$

Since *fi*(*s*) ≤ *gi*(*s*) for all *s* ∈ *S* and *i* = 1, 2 by [34] (Theorem 3.10) we have that

$$(\mathbb{R}\!L)\int\_{S} f\_{1}d\mu\_{1} \leq (\mathbb{R}\!L)\int\_{S} \mathbb{g}\_{1}d\mu\_{1}, \qquad (\mathbb{R}\!L)\int\_{S} f\_{2}d\mu\_{2} \leq (\mathbb{R}\!L)\int\_{S} \mathbb{g}\_{2}d\mu\_{2}.$$

and so by the weak interval order, **iii**), we have that *IF*(*S*) ! *IG*(*S*).

**Corollary 1.** *If F*, *G*, *F* ∧ *G*, *F* ∨ *G are RL integrable with respect to an interval valued multisubmeasure M then, for every A* ∈ A*,*

$$I(a) \qquad (\mathbb{R}L) \\ \int \,\_{S} F \wedge G \, dM \preceq I\_{F}(A) \wedge I\_{G}(A);$$

$$\text{(b)}\qquad I\_F(A) \lor I\_G(A) \preceq \text{(RL)} \int\limits\_S F \lor G \, d\mathcal{M}.$$

**Proof.** Let *F*(*s*)=[ *f*1(*s*), *f*2(*s*)], *G*(*s*)=[*g*1(*s*), *g*2(*s*)], *h*∗(*s*) = min{ *f*1(*s*), *g*1(*s*)}, *h*∗(*s*) = min{ *f*2(*s*), *g*2(*s*)}. By [34] (Theorem 3.10) (*RL*) *S h*∗*dμ*<sup>1</sup> ≤ (*RL*) *S f*1*dμ*1, (*RL*) *S g*1*dμ*<sup>1</sup> and an analogous result holds for (*RL*) *S h*∗*dμ*2. So the result given in 1.a) follows from the definition of ! and ∧.

The second statement follows analogously.

**Proposition 5.** *Let <sup>F</sup>*, *<sup>G</sup>* : *<sup>S</sup>* <sup>→</sup> *ck*(R<sup>+</sup> <sup>0</sup> ) *be bounded so that F*, *G are RL integrable with respect to M. If F* ⊆ *G*, *then IF*(*A*) ⊆ *IG*(*A*) *for all A* ∈ A*.*

**Proof.** As before we will prove for *S*. Let *ε* > 0 be arbitrary. Since *F*, *G* are *RL* integrable with respect to *M*, there exists a countable partition Π*<sup>ε</sup>* of *S* so that for every other countable partition <sup>Π</sup> = {*Bn*}*n*∈<sup>N</sup> ∈ P, with <sup>Π</sup> ≥ <sup>Π</sup>*<sup>ε</sup>* and every choise of points *sn* ∈ *Bn*, *<sup>n</sup>* ∈ N, the series

$$\sum\_{n=0}^{\infty} F(s\_n) \cdot M(B\_n)\_{\prime} \quad \sum\_{n=0}^{\infty} G(s\_n) \cdot M(B\_n)\_{\prime}$$

are convergent and

$$d\_H\left(I\_F(S), \sum\_{n=0}^{\infty} F(s\_n) \cdot M(B\_{\mathbb{N}})\right) < \frac{\varepsilon}{3}; \quad d\_H\left(I\_G(S), \sum\_{n=0}^{\infty} G(s\_n) \cdot M(B\_{\mathbb{N}})\right) < \frac{\varepsilon}{3}.$$

Then, by the triangular property of the eccess *e*,

$$\begin{split} e(I\_{F}(S), I\_{G}(S)) &\leq & d\_{H} \left( I\_{F}(S), \sum\_{n=0}^{\infty} F(s\_{n}) \cdot M(B\_{n}) \right) + e(\sum\_{n=0}^{\infty} F(s\_{n}) \cdot M(B\_{n}), \sum\_{n=0}^{\infty} G(s\_{n}) \cdot M(B\_{n})) + e(\sum\_{n=0}^{\infty} F(s\_{n}) \cdot M(B\_{n}), \sum\_{n=0}^{\infty} G(s\_{n}) \cdot M(B\_{n})) \\ &+ & d\_{H} \left( \sum\_{n=0}^{\infty} G(s\_{n}) \cdot M(B\_{n}), I\_{G}(S) \right) < \frac{2\varepsilon}{3} + e(\sum\_{n=0}^{\infty} F(s\_{n}) \cdot M(B\_{n}), \sum\_{n=0}^{\infty} G(s\_{n}) \cdot M(B\_{n})). \end{split}$$

Since the series ∞ ∑ *n*=0 *F*(*sn*)*· M*(*Bn*) and ∞ ∑ *n*=0 *<sup>G</sup>*(*sn*)*· <sup>M</sup>*(*Bn*) are convergent in *ck*(R<sup>+</sup> <sup>0</sup> ), and, by hypothesis, ∞ ∑ *n*=0 *F*(*sn*) *· M*(*Bn*) ⊆ ∞ ∑ *n*=0 *G*(*sn*) *· M*(*Bn*), then *e*( ∞ ∑ *F*(*sn*) *· M*(*Bn*), ∞ ∑ *G*(*sn*) *· M*(*Bn*)) = 0.

*n*=0 Consequently, from the arbitrariety of *ε* > 0, *e*(*IF*(*S*), *IG*(*S*)) = 0, which implies *IF*(*S*) ⊆ *IG*(*S*).

We can observe moreover that

**Proposition 6.** *If G is bounded and RL integrable with respect to M, with M of bounded variation, then*

$$\begin{aligned} (a) \qquad ||I\_G(\mathcal{S})|| &= (\mathbb{R} \mathcal{L}) \int\_{\mathcal{S}} \mathcal{g}\_2 \, d\mu\_2 = (\mathbb{R} \mathcal{L}) \int\_{\mathcal{S}} ||G|| \, d||\mathcal{M}||.\\ (b) \end{aligned} $$

*n*=0

$$\begin{aligned} \mathsf{T\_G}(S) &= \sup \{ \sum\_{i=1}^n |I\_{\mathbb{G}}(A\_i)|, \{A\_i, i = 1, \dots, n\} \in \mathcal{P} \} = \\ &= \sup \{ \sum\_{i=1}^n (\mathbb{R}\mathbb{L}) \int\_{A\_i} \mathbb{S}\_2 \, d\mu\_{2'} \, \left\{ A\_{i\prime} \, i = 1, \dots, n \right\} \in \mathcal{P} \} = (\mathbb{R}\mathbb{L}) \int\_{\mathcal{S}} \mathbb{S}\_2 \, d\mu\_2. \end{aligned}$$

**Proof.** It is a consequence of the properties of *dH* and [34] (Proposition 3.3, Theorem 3.5).

**Proposition 7.** *Let <sup>G</sup>* : *<sup>S</sup>* <sup>→</sup> *ck*(R<sup>+</sup> <sup>0</sup> ) *be a bounded multifunction such that G is RL integrable with respect to M on every set A* ∈ A*.*


**Proof.** The statements easily follow by Proposition 6.

Moreover

**Theorem 2.** *Let <sup>G</sup>* : *<sup>S</sup>* <sup>→</sup> *ck*(R<sup>+</sup> <sup>0</sup> ) *be a multifunction such that G is RL integrable with respect to M on every set A* ∈ A*. The following statements hold:*


**Proof.** Let *A*, *B* ∈ A with *A* ⊆ *B*. By monotonicity *μi*(*A*) ≤ *μi*(*B*) for *i* = 1, 2. We divide *B* in *A*, *B* \ *A* and we apply [34] (Theorem 3.2, Corollary 3.6). The conclusion follows by (**iii**).

Since *M* is a *dH*-multimeasure, then *M* is countably additive too and o-continuous. Applying Proposition <sup>7</sup> *IG* is o-continuous too. Let (*An*)*n*∈<sup>N</sup> ⊂ A be an arbitrary sequence of pairwise

disjoint sets, with ! ∞ *n*=1 *An* = *<sup>A</sup>* ∈ A. We denote by *Bn* the set *Bn* := ! ∞ *k*=*n*+1 *Ak*. Since *Bn* ∅, then lim*n*→<sup>∞</sup> *IG*(*Bn*) = 0. Since *IG* is finitely additive, we have

$$\lim\_{n \to \infty} d\_H(I\_G(A)\_\prime \sum\_{k=1}^n I\_G(A\_k)) = \lim\_{n \to \infty} d\_H(\sum\_{k=1}^n I\_G(A\_k) + I\_G(B\_n)\_\prime \sum\_{k=1}^n I\_G(A\_k)) \le \lim\_{n \to \infty} \|I\_G(B\_n)\| = 0$$

which ensures that *IG* is a *dH*-multimeasure.

Proceeding as in to the proof of the formula (5) and applying [34] (Theorem 3.8) we obtain the following result:

**Proposition 8.** *Let be <sup>M</sup>*1, *<sup>M</sup>*<sup>2</sup> : A → *ck*(R<sup>+</sup> <sup>0</sup> ) *, with <sup>M</sup>*1(∅) = *<sup>M</sup>*2(∅) = {0} *and suppose <sup>G</sup>* : *<sup>S</sup>* <sup>→</sup> *ck*(R<sup>+</sup> 0 ) *is RL integrable with respect to both <sup>M</sup>*<sup>1</sup> *and <sup>M</sup>*2*. If <sup>M</sup>* : A → *ck*(R<sup>+</sup> <sup>0</sup> ) *is the interval-valued multisubmeasure defined by M*(*A*) = *M*1(*A*) + *M*2(*A*)*, for every A* ∈ A*, then G is RL integrable with respect to M and*

$$\mathcal{S}(\mathbb{R}\mathcal{L}) \int\_{\mathcal{S}} \mathcal{G} \, d(\mathcal{M}\_1 + \mathcal{M}\_2) = (\mathbb{R}\mathcal{L}) \int\_{\mathcal{S}} \mathcal{G} \, dM\_1 + (\mathbb{R}\mathcal{L}) \int\_{\mathcal{S}} \mathcal{G} \, dM\_2.$$

**Theorem 3.** *Let <sup>M</sup> be of bounded variation and <sup>F</sup>*, *<sup>G</sup>* : *<sup>T</sup>* <sup>→</sup> *ck*(R<sup>+</sup> <sup>0</sup> ) *be bounded interval-valued multifunctions. If F*, *G are RL integrable with respect to M, then*

$$d\_H\left( (\mathbb{R}\boldsymbol{L}) \int\_S \boldsymbol{F} d\boldsymbol{M}, (\mathbb{R}\boldsymbol{L}) \int\_S \boldsymbol{G} \, d\boldsymbol{M} \right) \le \sup\_{\boldsymbol{s} \in \mathcal{S}} d\_H(\boldsymbol{F}(\boldsymbol{s}), \boldsymbol{G}(\boldsymbol{s})) \cdot \overline{\boldsymbol{M}}(\boldsymbol{S}).$$

**Proof.** Since *F*, *G* are *M*-integrable then *f*1, *g*<sup>1</sup> are *μ*1-integrable and *f*2, *g*<sup>2</sup> are *μ*2-integrable functions. According to [34] (Theorem 3.9), we have for *i* = 1, 2,

$$\left| (\mathbb{R} \mathbb{L}) \int\_{S} f\_{\mathrm{i}} \, d\mu\_{\mathrm{i}} - (\mathbb{R} \mathbb{L}) \int\_{S} g\_{\mathrm{i}} \, d\mu\_{\mathrm{i}} \right| \leq \sup\_{s \in S} |f\_{\mathrm{i}}(s) - g\_{\mathrm{i}}(s)| \overline{\mu}\_{\mathrm{i}}(S). \tag{6}$$

Therefore, by (6) and Remark 5, it follows

$$\begin{split} \operatorname{Ad}\_{H}\left( (\operatorname{KL}) \int\_{S} \operatorname{F} \operatorname{dM}\_{\tau}(\mathbb{R}\mathbf{1}) \int\_{S} \mathbb{G} \operatorname{dM} \right) &= \max \left\{ \left| (\operatorname{KL}) \int\_{S} f\_{1} \operatorname{f} \mu\_{1} - (\operatorname{\mathbb{KL}}) \int\_{S} \operatorname{g} \operatorname{1} \mu\_{1} \right|, \left| (\operatorname{\mathbb{KL}}) \int\_{S} f\_{2} \operatorname{f} \mu\_{2} - (\operatorname{\mathbb{KL}}) \int\_{S} \operatorname{g} \operatorname{2} \mu\_{2} \right| \right\} \\ &\leq \max \left\{ \sup\_{s \in S} |f\_{1}(s) - \operatorname{g}\_{1}(s)| \overline{\mu\_{1}}(S), \sup\_{s \in S} |f\_{2}(s) - \operatorname{g}\_{2}(s)| \overline{\mu\_{2}}(S) \right\} \leq \\ &\leq \max \left\{ \sup\_{s \in S} |f\_{1}(s) - \operatorname{g}\_{1}(s)|, \sup\_{s \in S} |f\_{2}(s) - \operatorname{g}\_{2}(s)| \right\} \overline{\mu\_{2}}(S) \leq \\ &= \sup\_{s \in S} d\_{H}(F(s), G(s)) \overline{M}(S). \end{split}$$

**Theorem 4.** *Let <sup>M</sup>*1, *<sup>M</sup>*<sup>2</sup> : A → *ck*(R<sup>+</sup> <sup>0</sup> ) *and <sup>G</sup>* : *<sup>S</sup>* <sup>→</sup> *ck*(R<sup>+</sup> <sup>0</sup> ) *be RL integrable with respect to both M*<sup>1</sup> *and M*2*. Then*

$$\begin{aligned} (a) \qquad &\text{If } M\_1 \preceq M\_2, \text{ then } (\text{RL}) \int\_S \text{G } dM\_1 \preceq (\text{RL}) \int\_S \text{G } dM\_2. \\ (b) \qquad &\text{If } M\_1 \subseteq M\_2, \text{ then } (\text{RL}) \int\_S \text{G } dM\_1 \subseteq (\text{RL}) \int\_S \text{G } dM\_2. \end{aligned}$$

**Proof.** Let *M*<sup>1</sup> := [*μ*∗, *μ*∗] and *M*<sup>2</sup> := [*ν*∗, *ν*∗]. Both the results are consequences of Theorem 2 and [34] (Theorem 3.11). It is enough to observe that if *M*<sup>1</sup> ! *M*<sup>2</sup> then *μ*<sup>∗</sup> ≤ *ν*<sup>∗</sup> and *μ*<sup>∗</sup> ≤ *ν*∗, while if *M*<sup>1</sup> ⊆ *M*<sup>2</sup> then *ν*<sup>∗</sup> ≤ *μ*<sup>∗</sup> ≤ *μ*<sup>∗</sup> ≤ *ν*∗.

As a particular case of Theorem 4 and Corollary 1 we have that for every *G* which is *RL* integrable with respect to both positive submeasures *μ*<sup>1</sup> and *μ*<sup>2</sup> then

$$\delta(\mathbb{R}\mathcal{L}) \int\limits\_{S} G \, d(\mu\_1 \wedge \mu\_2) \preceq (\mathbb{R}\mathcal{L}) \int\limits\_{S} G \, d\mu\_1 \wedge (\mathbb{R}\mathcal{L}) \int\limits\_{S} G \, d\mu\_2 \dots$$

Moreover a convergence result can be obtained using Proposition 1.

**Theorem 5.** *Let Gn* = [*g* (*n*) <sup>1</sup> , *g* (*n*) <sup>2</sup> ] *be a sequence of bounded RL-integrable interval valued multifunction with respect to M* = [*μ*1, *μ*2] *such that Gn* ! *Gn*+<sup>1</sup> *for every n* ∈ N*. If M is of bounded variation and there exists a function G* = [*g*1, *g*2] *such that:*


*then G is RL-integrable with respect to M and*

$$\lim\_{n \to \infty} d\_H \left( (\mathbb{R} \boldsymbol{L}) \int \mathop{\rm G}\_{\mathcal{S}} \mathop{\rm d} \nolimits \nolimits \mathcal{M}\_{\boldsymbol{\nu}} (\mathbb{R} \boldsymbol{L}) \int \mathop{\rm G}\_{\mathcal{S}} \mathop{\rm d} \nolimits \mathcal{M} \right) = 0.$$

**Proof.** Since *Gn* ! *Gn*+<sup>1</sup> we have that *g* (*n*) *<sup>i</sup>* ↑ for *i* = 1, 2, this is a consequence of Proposition 4 and Definition 6. By *dH*(*Gn*, *G*) → 0 uniformly we have that max{|*g* (*n*) *<sup>i</sup>* − *gi*|, *i* = 1, 2} converges uniformly to zero. We can use now Proposition 1 and we obtain

$$\lim\_{n \to \infty} \, \_{(\mathbb{R}\boldsymbol{L})} \int \, \_{\mathrm{S}} \prescript{(n)}{}{\mathrm{d}} d\mu\_{i} = \, (\mathbb{R}\boldsymbol{L}) \int \, \_{\mathrm{S}} \prescript{}{}{\mathrm{d}} d\mu\_{i\prime} \quad i = 1, 2.$$

For every *ε* > 0 let *k*(*ε*) ∈ N be such that

$$d\_H(\mathcal{G}(t), \mathcal{G}\_{k(\varepsilon)}(t)) < \varepsilon \,\,\forall t \in \mathcal{S}\_{\prime} \quad \text{and} \quad \left| (\mathbb{R} \mathcal{L}) \int\_{\mathcal{S}} \mathcal{g}\_i^{(k(\varepsilon))} d\mu\_i - (\mathbb{R} \mathcal{L}) \int\_{\mathcal{S}} \mathcal{g}\_i d\mu\_i \right| < \varepsilon, \; i = 1, 2.$$

So,

$$\operatorname{ad}\_H \left( (\mathbb{R} \mathbf{L}) \int\_S \mathbf{G}\_{k(\varepsilon)} \, d\mathcal{M}\_\prime \left[ (\mathbb{R} \mathbf{L}) \int\_S \mathcal{g}\_1 d\mu\_{1\prime\prime}(\mathbb{R} \mathbf{L}) \int\_S \mathcal{g}\_2 d\mu\_2 \right] \right) \le \varepsilon.$$

Let *P<sup>ε</sup>* be the countable partition of *S* given by / *<sup>i</sup>*=1,2 *Pε*,*i*, (the ones that verify Definition 4 for *g k*(*ε*) *<sup>i</sup>* , *<sup>i</sup>* = 1, 2 respectively). Then, for every countable partition *<sup>P</sup>* = {*An*}*n*∈<sup>N</sup> of *<sup>S</sup>* with *<sup>P</sup>* ≥ *<sup>P</sup><sup>ε</sup>* and for every *tn* ∈ *An* the series *σG*,*M*(*P*) is convergent and

$$\begin{split} d\_{H}\left(\sigma\_{\operatorname{G},\operatorname{M}}(\operatorname{P})\_{\operatorname{r}}\left[ (\operatorname{\mathbb{R}L})\int\limits\_{\operatorname{\mathcal{S}}}\operatorname{\mathcal{S}}\_{1}d\mu\_{1\star}(\operatorname{\mathbb{R}L})\int\limits\_{\operatorname{\mathcal{S}}}\operatorname{\mathcal{S}}\_{2}d\mu\_{2} \right] \right) &< d\_{H}\left(\sigma\_{\operatorname{G},\operatorname{M}}(\operatorname{P}),\sigma\_{\operatorname{G}\_{\operatorname{\mathbb{R}}(\boldsymbol{\varepsilon})},\operatorname{M}}(\operatorname{P})\right) + \\ &+ d\_{H}\left(\sigma\_{\operatorname{G}\_{\operatorname{\mathbb{R}}(\boldsymbol{\varepsilon})},\operatorname{M}}(\operatorname{P}),\left(\operatorname{\mathbb{R}L}\right)\int\limits\_{\operatorname{\mathcal{S}}}\operatorname{\mathcal{G}\_{\operatorname{\mathbb{R}}(\boldsymbol{\varepsilon})}}\operatorname{\mathcal{S}}\_{1}d\mu\_{1}\right) + \\ &+ d\_{H}\left(\left(\operatorname{\mathbb{R}L}\right)\int\limits\_{\operatorname{\mathcal{S}}}\operatorname{\mathcal{G}\_{\operatorname{\mathbb{R}}(\boldsymbol{\varepsilon})}}\operatorname{\mathcal{M}}\_{\operatorname{r}}\left(\operatorname{\mathcal{S}}1d\mu\_{1},\left(\operatorname{\mathcal{S}}1\right)\int\limits\_{\operatorname{\mathcal{S}}}\operatorname{\mathcal{S}}\_{2}d\mu\_{2}\right)\right). \end{split}$$

From previous inequalities and by the arbitrariety of *ε* the *RL*-integrability of *G* follows.

**Remark 8.** *Since this research starts from the papers [34,43], this part ends with a comparison between the two types of integral considered: the RL integral with the Gould one given in [43] (Definition 4.7).*

*If the interval-valued multifunction F is bounded and μ*<sup>2</sup> *is of finite variation then, analogously to Proposition 2 it is, by [43] (Proposition 4.9),*

$$\mathcal{E}(\mathcal{G}) \int\_{S} \mathcal{F} dM = \left[ (\mathcal{G}) \int\_{S} f\_1 \, d\mu\_1 \, / \mathcal{G} \right] \int\_{S} f\_2 d\mu\_2 \right].$$

*So, the two kinds of integral coincide on bounded interval-valued multifunctions with values in ck*(R<sup>+</sup> <sup>0</sup> ) *when μi*, *i* = 1, 2 *are complete countably additive measures by [34] (Proposition 4.5) or μi*, *i* = 1, 2 *are monotone, countably -subadditive by [34] (Theorem 4.7).*

*Without countable additivity the equivalence does not hold; an example can be constructed using [34] (Example 4.6). In the general case only partial results can be obtained on atoms when μi*, *i* = 1, 2 *are monotone, null additive and satisfy property* (*σ*)*: the proof follows from [34] (Theorem 4.8).*

*Accordingly with the comparison between Gould and Birkoff integrals given in [28] we have that Birkhoff, Gould, RL integrals of the bounded single valued functions agree in the countably additive case, see [28] (Theorem 3.10), while in [43] (Remark 5.5) an analogous comparison is given with the Choquet integral.*

*A comparison between simple Birkhoff and RL integrabilities, introduced in [23,28], in this non additive setting can be obtained using [34] (Theorem 4.2).*

Finally we would like to observe that the Rådström's embedding tell us that (*ck*(*X*), *dH*, ⊆), when *X* finite dimensional, is a near vector space with 0 element and order unit *BX*. In this case, using [51] (Theorem 5.1), it is a near vector space (see [51] (Definition 2.1) for its definition) that could be embedded, for example, in -<sup>∞</sup> or in *C*(Ω) with Ω compact and Hausdorff in such a way the embedding is an isometric isomorphism which takes into account the ordering on the hyperspace.

If we consider instead (*ck*(R<sup>+</sup> <sup>0</sup> ), *dH*, !), since in general there is no relation between "!" and "⊆" the Rådström embedding provide only the integrability of the interval-valued functions and does not take the weak interval order into account. For this reason we preferred to give the the construction of the *RL* integral and the proofs, both related to !, independently of the Rådström's embedding.

#### *3.1. Applications of Interval Valued Multifunctions*

Now, in order to explain what could be the benefits of this approach we give an example of an application of interval valued multifunctions on interval valued multisubmeasure in image processing. In fact a signal can be modeled as an interval-valued multifunction as in [12]. In fact, when the value of the points can not be assigned with precision, it might be preferable to use a measure-based approach.

The advantage of using the notion of interval-valued multifunction in signal analysis is that this formalism allows to include in a unique framework possible uncertainty or the noise on the value of a point.

This situation usually occurs in signal and image processing when images are derived by a measure process, as happens for instance for biomedical images (in CT images, MR images, etc), and in several other applied sciences. In particular, we can apply this representation to a digital image in such a way:

**Example 2.** *To each pixel (or to a set of pixels) of the image is associated an interval which measures the* round-off error *which is that committed on the detection on the signal due by the tolerances and by the limits on computational accuracy of the measurements tools ([52]).*

*When we consider* subsets of pixels *we are taking into account the so-called* time-jitter error*, i.e., the error that occur in the measure of a given signal when the sampling values can not be matched exactly at the* theoretical node *but just in a neighborhood of it (see, e.g., [53]).*

*In this sense, if I* = (*mi*,*j*) *is the matrix associated to a n* × *m static, gray-scale image, we can consider the space* <sup>S</sup> := (0, *<sup>n</sup>*] <sup>×</sup> (0, *<sup>m</sup>*] <sup>⊂</sup> <sup>R</sup>2*, and hence the interval-valued multifunction UI* : S→K<sup>+</sup> *<sup>C</sup> corresponding to I, will be given by:*

$$\mathcal{U}\_{\mathsf{I}}(\mathsf{x}) := [\mathsf{u}\_{\mathsf{I}}(\mathsf{x}), \mathsf{u}\_{\mathsf{I}}(\mathsf{x})], \qquad \mathsf{x} \in \mathsf{S}.$$

*The model of a digital image by an interval-valued multifunction as UI, and obtained by a certain discretization (algorithm) of an analogue image, allows to control the round-off error in the sense that, the true value assumed by original signal at the pixel x belongs to the interval* [*u*1(*x*), *u*2(*x*)]*, in fact providing a lower and an upper bound on the possible oscillations of the sampled image.*

*For example, in fractal image coding, the functions u*<sup>1</sup> *and u*<sup>2</sup> *represent respectively the lower and upper contraction maps of an image, which take into account of the round-off error in the contraction procedure, and can be chosen as follows:*

$$\mu\_1(\mathbf{x}) := \mathfrak{a}\_1 \mathfrak{u}(\mathbf{x}) + \beta\_1(\mathbf{x}), \quad \mathfrak{u}\_2(\mathbf{x}) := \mathfrak{a}\_2 \mathfrak{u}(\mathbf{x}) + \beta\_2(\mathbf{x}), \quad \mathbf{x} \in \mathcal{S}\_{\epsilon}$$

*where αi, i* = 1, 2*, are suitable integer scaling parameters, β<sup>i</sup>* : S → N*, i* = 1, 2*, are suitable functions, and u* : S → N *is the continuous model associated to the starting image I. The functions u*<sup>1</sup> *and u*<sup>2</sup> *provide for each pixel the interval containing the true value of the compressed image.*

*In particular, in the algorithm considered in [15], the functions u*<sup>1</sup> *and u*<sup>2</sup> *are piecewise constant, and for a starting image of* 225 × 225 *pixel size, they have been defined as follows:*

$$\mathcal{U}I\_l(\mathbf{x}) = [\mu\_1(\mathbf{x}), \mu\_2(\mathbf{x})] = [\mu(\mathbf{x}) - \beta(\mathbf{x}), \mu(\mathbf{x}) + \beta(\mathbf{x})], \qquad \mathbf{x} \in \mathcal{S}, \tag{7}$$

*where:*

$$\mu(\mathbf{x}) := m\_{\hat{i}, \hat{j}}, \quad \mathbf{x} \in (\mathbf{i} - \mathbf{1}, \mathbf{i}] \times (\mathbf{j} - \mathbf{1}, \mathbf{j}], \quad \mathbf{i}, \mathbf{j} = 1, \dots, 225,$$

*and*

$$\beta(\mathbf{x}) := \begin{cases} \mathbf{0}, & \mathbf{x} \in (0, 115] \times (0, 115], \\ \mathbf{40}, & \mathbf{x} \in (115, 225] \times (115, 225], \\ \mathbf{20}, & \text{otherwise}. \end{cases} \tag{8}$$

*As an example we use the interval-valued multifunction (7) to operate with the well-known image of "Baboon" given in Figure 1 (left); the images generated by u*<sup>1</sup> *and u*<sup>2</sup> *using the function β defined in (8) are given in Figure 1 (center and right).*

**Figure 1.** Baboon (**left**); The images generated by *u*<sup>1</sup> (**center**) and *u*<sup>2</sup> (**right**) using the interval valued multifunction (7), with *β* defined in (8).

*Here, also numerical truncation have been taken into account, in order to maintain the values of the pixels in the (integer) gray scale* [0, 255]*.*

*For other examples of functions u*<sup>1</sup> *and u*2*, see, e.g., [13,54]. For instance, in [13] the image representation by multifunctions is used for the implementation of edge detection algorithms, and in this case the corresponding functions u*<sup>1</sup> *and u*<sup>2</sup> *are:*

$$u\_1(\mathbf{x}) := \max\left\{0, \min\_{\mathbf{x'} \in n(\mathbf{x})} \left\{ I(\mathbf{x'}) - 1 \right\} \right\}, \qquad u\_2(\mathbf{x}) := \min\left\{ 255, \max\_{\mathbf{x'} \in n(\mathbf{x})} \left\{ I(\mathbf{x'}) + 1 \right\} \right\},$$

*where I*(*x*) *represents the value of a pixel at a position x* ∈ S*, while n*(*x*) *denotes any set of* 3 × 3 *pixels centered at x. For more details, or other applications, see [13,18].*

This example was built with the aim to highlight a useful link between the abstract theory of the interval-valued multifunction and the concrete application to image processing. One of the crucial tool in the above set-valued theory is provided by the Hausdorff distance between sets. This special metric plays an important role in the context of digital image processing, where it is used, for example, in order to measure the accuracy of certain class of algorithms, such as those of edge detection, already mentioned in the previous list of possible applications. More precisely, if *A* is the region of interest (ROI) of a given image and *B* is the corresponding approximation of the ROI *A* detected by a suitable edge detection algorithm, the Hausdorff distance measure the displacement between *A* and *B*, in fact evaluating the accuracy (i.e., the approximation error) of the method. For instance, in [55] the Hausdorff distance has been used in order to evaluate the degree of accuracy of an algorithm for the detection of the pervious area of the aorta artery from CT images without contrast medium. This procedure is useful, for example, in the diagnosis of aneurysms of the abdominal aorta artery, especially for patients with severe kidneys pathology for which CT images with contrast medium can not be performed. A similar use of the Hausdorff distance could be done for the edge detection algorithms considered in [13,18].

#### **4. Conclusions**

A Riemann Lebesgue integral is defined for interval-valued multifunction with respect to interval-valued multisubmeasures. Properties of the integral are established showing in particular that the multimeasure generated is finitely additive. Sufficient conditions for the monotonicity, the order continuity, bounded variation and convergence results are also obtained. A comparison with other integrals is sketchced; an example of an applications in image processing is given highlighting that the advantage of using the notion of interval-valued multifunction in signal analysis is that this formalism allows to include in a unique framework possible uncertainty or the noise on the evaluation of an image at any given pixel. In a future research we will generalize these results in the setting of Banach lattices and we will compare this method with other DIP (digital image processing) algorithms.

**Author Contributions:** Writing, review and editing, D.C., A.C., A.G., A.I. and A.R.S. All authors have read and agreed to the published version of the manuscript.

**Funding:** The first and the last authors are members of the working group "Teoria dell'Approssimazione e Applicazioni" of the Italian Mathematical Union (U.M.I.) and they were partially supported by: Grant "Analisi reale, teoria della misura ed approssimazione per la ricostruzione di immagini" (2020) of GNAMPA – INDAM (Italy) and University of Perugia—Fondo Ricerca di Base 2019.

**Conflicts of Interest:** The authors declare no conflict of interest.

**Copyright of the Pictures**: The Figure 1 (baboon) are adapted from the original image contained in the repository http://calvados.di.unipi.it/dokuwiki/lib/exe/detail.php?id=ffnamespace:apps&media=ffnamespace: baboon.jpg under the following license: CC Attribution-Share Alike 3.0 Unported.

#### **References**


**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
