1 Introduction
The sampling Kantorovich operators have been introduced to approximate and reconstruct not necessarily continuous signals. In [4], the authors introduced this operators starting from the wellknown generalized sampling operators (see e.g. [11, 13, 10, 23, 5]) and replacing, in their definition, the sample values with . Clearly, this is the most natural mathematical modification to obtain operators which can be welldefined also for general measurable, locally integrable functions, not necessarily continuous. Moreover, this situation very often occur in Signal Processing, when one cannot match exactly the sample at the point : this represents the socalled "timejitter” error. The theory of sampling Kantorovich operators allow us to reduces the timejitter error, calculating the information in a neighborhood of rather that exactly in the node . These operators, as the generalized sampling operators, represent an approximate version of classical sampling series, based on the WhittakerKotelnikovShannon sampling theorem (see e.g. [2]).
Subsequently, the sampling Kantorovich operators have been studied in various settings. In [16, 17] the multivariate version of these operators were introduced. Results concerning the order of approximation are shown in [18]. Extensions to more general contexts are presented in [24, 25, 26, 7].
The multivariate sampling Kantorovich operators considered in this paper are of the form:
where is a locally integrable function such that the above series is convergent for every . The symbol
denotes vectors where each
, is a certain strictly increasing sequence of real numbers with . Note that, the sequences are not necessary equally spaced. We denote by the sets:and , . Moreover, the function is a kernel satisfying suitable assumptions.
For the operators in (I) we recall some convergence results. We have that the family converges pointwise to , when is continuous and bounded and converges uniformly to , when is uniformly continuous and bounded. Moreover, to cover the case of not necessarily continuous signal, we study our operators in the general setting of Orlicz spaces . For functions belonging to and generated by the convex function , the family of sampling Kantorovich operators is "modularly” convergent to , being the latter the natural concept of convergence in this setting.
The latter result, allow us to apply the theory of the sampling Kantorovich operators to approximate and reconstruct images. In fact, static gray scale images are characterized by jumps of gray levels mainly concentrated in their contours or edges and this can be translated, from a mathematical point of view, by discontinuities (see e.g. [19]).
Here, we introduce and analyze in detail some practical applications of the sampling Kantorovich algorithm to thermographic images, very useful for the analysis of buildings in seismic engineering. The thermography is a remote sensing technique, performed by the image acquisition in the infrared. Thermographic images are widely used to make noninvasive investigations of structures, to analyze the story of the building wall, to make diagnosis and monitoring buildings, and to make structural measurements. A further important use, is the application of the texture algorithm for the separation between the bricks and the mortar in masonries images. Through this procedure the civil engineers becomes able to determine the mechanical parameters of the structure under investigation. Unfortunately, the direct application of the texture algorithm to the thermographic images, can produce errors, as an incorrect separation between the bricks and the mortar.
Then, we use the sampling Kantorovich operators to process the thermographic images before to apply the texture algorithm. In this way, the result produced by the texture becomes more refined and therefore we can apply structural analysis after the calculation of the various parameters involved. In order to show the feasibility of our applications, we present in detail a realworld casestudy.
2 Preliminaries
In this section we recall some preliminaries, notations and definitions.
We denote by (resp. ) the space of all uniformly continuous and bounded (resp. continuous and bounded) functions endowed with the usual supnorm , , and by the subspace of the elements having compact support. Moreover, will denote the linear space of all (Lebesgue) measurable real functions defined on .
The function is said to be a function if it satisfies the following assumptions: , and for every ; is continuous and non decreasing on ; .
The functional (see e.g. [20, 8]) defined by
is a modular in . The Orlicz space generated by is given by
The space is a vector space and an important subspace is given by
is called the space of all finite elements of . It is easy to see that the following inclusions hold: Clearly, functions belonging to and are not necessarily continuous. A norm on , called Luxemburg norm, can be defined by
We will say that a family of functions is norm convergent to a function , i.e., for , if and only if , for every . Moreover, an additional concept of convergence can be studied in Orlicz spaces: the "modular convergence". The latter induces a topology (modular topology) on the space ([20, 8]).
We will say that a family of functions is modularly convergent to a function if , for some . Obviously, norm convergence implies modular convergence, while the converse implication does not hold in general. The modular and norm convergence are equivalent if and only if the function satisfies the condition, see e.g., [20, 8]. Finally, as last basic property of Orlicz spaces, we recall the following.
Lemma 2.1 ([6]).
The space is dense in with respect to the modular topology, i.e., for every and for every there exists a constant and a function such that .
3 The sampling Kantorovich operators
In this section we recall the definition of the operators with which we will work. We will denote by a vector where each element , is a sequence of real numbers with , , for every , and such that there exists , for which , for every . Note that, the elements of are not necessary equally spaced. In what follows, we will identify with the symbol the sequence .
A function will be called a kernel if it satisfies the following properties:

and is bounded in a neighborhood of ;

for every , ;

for some ,
where denotes the usual Euclidean norm.
We now recall the definition of the linear multivariate sampling Kantorovich operators introduced in [16]. Define:
(1) 
where is a locally integrable function such that the above series is convergent for every , where
and , . The operators in (1) have been introduced in [4] in the univariate setting.
Remark 3.1.
(a) Under conditions and , the following properties for the kernel can be proved:
and, for every
(2) 
uniformly with respect to , see [16]. (b) By (a), we obtain that with are welldefined. Indeed,
for every , i.e. .
Remark 3.2.
Note that, in the onedimensional setting, choosing , for every , condition is equivalent to
where ,
, denotes the Fourier transform of
; see [12, 4, 16, 18].4 Convergence results
In this section, we show the main approximation results for the multivariate sampling Kantorovich operators. In [16], the following approximation theorem for our operators has been proved.
Theorem 4.1.
Let . Then, for every ,
In particular, if , then
Proof.
Here we highlight the main points of the proof.
Let and be fixed. By the continuity of we have that for every fixed there exists such that for every , . Then, by we obtain:
For and we have for sufficiently large, then by the continuity of we obtain (see Remark 3.1 (a)). Moreover, by the boundedness of and (2) we obtain for sufficiently large, then the first part of the theorem follows since is arbitrary. The second part of the theorem follows similarly replacing with the parameter of the uniform continuity of . ∎
In order to obtain a modular convergence result, the following normconvergence theorem for the sampling Kantorovich operators (see [16]) can be formulated.
Theorem 4.2.
Let be a convex function. For every we have
Now, we recall the following modular continuity property for , useful to prove the modular convergence for the above operators in Orlicz spaces.
Theorem 4.3.
Let be a convex function. For every there holds
for some . In particular, maps in .
Now, the main result of this section follows (see [16]).
Theorem 4.4.
Let be a convex function. For every , there exists such that
Proof.
The setting of Orlicz spaces allows us to give a unitary approach for the reconstruction since we may obtain convergence results for particular cases of Orlicz spaces. For instance, choosing , , we have that and , where is the usual norm. Then, from Theorem 4.3 and Theorem 4.4 we obtain the following corollary.
Corollary 4.5.
For every , , the following inequality holds:
Moreover, we have:
The corollary above, allows us to reconstruct signals (in sense), therefore signals/images not necessarily continuous. Other examples of Orlicz spaces for which the above theory can be applied can be found e.g., in [21, 20, 8, 4, 16]. The theory of sampling Kantorovich operators in the general setting of Orlicz spaces allows us to obtain, by means of a unified treatment, several applications in many different contexts.
5 Examples of special kernels
One important fact in our theory is the choice of the kernels, which influence the order of approximation that can be achieved by our operators (see e.g. [18] in onedimensional setting).
For instance, one can take into consideration radial kernels, i.e., functions for which the value depends on the Euclidean norm of the argument only. Example of such a kernel can be given, for example, by the BochnerRiesz kernel, defined as follows , for , where , is the Bessel function of order and is the Euler function. For more details about this matter, see e.g. [11].
To construct, in general, kernels satisfying all the assumptions , is not very easy.
For this reason, here we show a procedure useful to construct examples using product of univariate kernels, see e.g. [11, 16, 17]. In this case, we consider the case of uniform sampling scheme, i.e., .
Denote by , the univariate functions , , satisfying the following assumptions:
(3) 
for some , is bounded in a neighborhood of the origin and
(4) 
for every , for . Now, setting , , it is easy to prove that is a multivariate kernel for the operators satisfying all the assumptions of our theory, see e.g., [11, 16].
As a first example, consider the univariate Fejér’s kernel defined by , , where the sinc function is given by
Clearly, is bounded, belongs to
and satisfies the moment conditions
for , as shown in [12, 4, 16]. Furthermore, taking into account that the Fourier transform of is given by (see [12])we obtain by Remark 3.2 that condition is fulfilled. Then, we can define , the multivariate Fejér’s kernel, satisfying the condition upon a multivariate kernel. The Fejér’s kernel and the bivariate Fejér’s kernel are plotted in Figure 1.
The Fejér’s kernel is an example of kernel with unbounded support, then to evaluate our sampling Kantorovich series at any given , we need of an infinite number of mean values . However, if the function has compact support, this problem does not arise. In case of function having unbounded support the infinite sampling series must be truncated to a finite one, which leads to the socalled truncation error. In order to avoid the truncation error, one can take kernels with bounded support. Remarkable examples of kernels with compact support, can be constructed using the wellknown central Bspline of order , defined by
where the function denotes the positive part of (see [4, 24, 16]). The central Bspline and the bivariate Bspline kernel are plotted in Figure 2.
We have that, the Fourier transform of is given by , , and then, if we consider the case of the uniform spaced sampling scheme, condition (4) is satisfied by Remark 3.2. Clearly, are bounded on , with compact support , and hence , for all . Moreover, it easy to deduce that condition is fulfilled for every , see [4]. Hence , , is the multivariate Bspline kernel of order .
Finally, other important examples of univariate kernels are given by the Jacksontype kernels, defined by , , with , , where the normalization coefficients are given by . Since are bandlimited to , i.e., their Fourier transform vanishes outside this interval, condition is satisfied by Remark 3.2 and is satisfied since , as , . In similar manner, by the previous procedure we can construct a multivariate version of Jackson kernels. For more details about Jacksontype kernels, and for others useful examples of kernels see e.g. [12, 8, 3, 4].
6 The sampling Kantorovich algorithm for image reconstruction
In this section, we show applications of the multivariate sampling Kantorovich operators to image reconstruction. First of all, we recall that every bidimensional gray scale image is represented by a suitable matrix and can be modeled as a step function , with compact support, belonging to , . The definition of arise naturally as follows:
where ,
, are the characteristics functions of the sets
(i.e. , for and otherwise).The above function is defined in such a way that, to every pixel it is associated the corresponding gray level .
We can now consider approximation of the original image by the bivariate sampling Kantorovich operators based upon some kernel .
Then, in order to obtain a new image (matrix) that approximates in sense the original one, it is sufficient to sample (for some ) with a fixed sampling rate. In particular, we can reconstruct the approximating images (matrices) taking into consideration different sampling rates and this is possible since we know the analytic expression of .
If the sampling rate is chosen higher than the original sampling rate, one can get a new image that has a better resolution than the original one’s. The above procedure has been implemented by using MATLAB and tools of matrix computation in order to obtain an optimized algorithm based on the multivariate sampling Kantorovich theory.
In the next sections, examples of thermographic images reconstructed by the sampling Kantorovich operators will be given to show the main applications of the theory to civil engineering. In particular, a real world casestudy is analyzed in term of modal analysis to obtain a model useful to study the response of the building under seismic action.
7 An application of thermography to civil engineering
In the present application, thermographic images will be used to locate the resisting elements and to define their dimensions, and moreover to investigate the actual texture of the masonry wall, i.e., the arrangement of blocks (bricks and/or stones) and mortar joints. In general, thermographic images have a resolution too low to accurately analyze the texture of the masonries, therefore it has to be increased by means of suitable tools.
In order to obtain a consistent separation of the phases, that is a correct identification of the pixel which belong to the blocks and those who belong to mortar joints, the image is converted from grayscale representation to blackandwhite (binary) representation by means of an image texture algorithm, which employs techniques belonging to the field of digital image processing. The image texture algorithm, described in details in [14], leads to areas of white pixels identified as blocks and areas of black pixels identified as mortar joints. However, the direct application of the image texture algorithm to the thermographic images (see Figure 3), can produce errors, as an incorrect separation between the bricks and the mortar (see Figure 3). Therefore, we can use the sampling Kantorovich operators to process the thermographic images. In particular, here we used the operators based upon the bivariate Jacksontype kernel with (see Section 5) and the parameter (see Figure 3). The application of the image texture algorithm produces a consistent separation of the phases (see Figure 3).
In order to perform structural analysis, the mechanical characteristics of an homogeneous material equivalent to the original heterogeneous material are sought. The equivalence is in the sense that, when subjected to the same boundary conditions (b.c.), the overall responses in terms of mean values of stresses and deformations are the same. The equivalent elastic properties are estimated by means of the “testwindow” method
[15].8 Casestudy
The image texture algorithm described previously has been used to analyze a realworld casestudy: a building consisting of two levels and an attic, with a very simple architectural plan, a rectangle with sides 11 m and 11.4 m (see Figure 4). The vertical structural elements consist of masonry walls. At first level the masonry walls have thickness of 40 cm with blocks made of stones, while at second level the masonry is made of squared tuff block and it has thickness of about 35 cm. Both surfaces of the walls are plastered. The slab can be assumed to be rigid in the horizontal plane and the building is placed in a medium risk seismic area. According to a preliminary visual survey, the structure does not have strong asymmetries between principal directions and the distribution of the masses are quite uniform.
The response of the building is evaluated by means of three different models. In the first model, only information that can be gathered by the visual survey have been used, and the characteristics of material has been assume according to Italian Building Code. In the second model information concerning the actual geometry of vertical structural elements acquired by means of thermographic survey have been used; for example, the survey showed that the openings at ground level in the main facade were once larger and have been subsequently partially filled with nonstructural masonry, and thus the dimensions of all the masonry walls in the facade have to be reduced. Eventually, in the third model information about the actual texture of masonries of ground level (chaotic masonry, Fig. 5) and of second level (periodic masonry, Fig. 5) have been used. In particular, the actual textures where established using the reconstructed thermographic images.
The main characteristics of the models are reported in Table 1. For the mechanical characteristic the Young’s modulus and the shear modulus are shown. In the third model, the following mechanical characteristics (Young’s modulus and Poisson’s ratio ) of the constituent phases have been used: for the stones of the ground level, , for the bricks of the second level , for the mortar of both levels .
ground level  second level  
geometry  
Model #1  visual survey  3346  1115  1620  540 
Model #2  thermogr. survey  3346  1115  1620  540 
Model #3  thermogr. survey  7050  2957  1996  833 
The behavior under seismic actions is estimated by means of modal analysis, using a commercial code based on the Finite Element Method. The periods of the first three modes and the corresponding mass participating ratio for the two principal direction of seismic action are reported in Tab. 2. For all the models, the first mode is along the axis, the second is along the axis, and the third is mainly torsional. Anyway, the second mode is not a pure translational one since the asymmetry in the walls distribution produces also a torsional component.
Periods  Mass partecipating ratio  

direction  direction  
1st mode  2nd mode  3rd mode  1st mode  2nd mode  3rd mode  1st mode  2nd mode  3rd mode  
Model #1  0.22  0.15  0.13  0.00  0.64  0.13  0.73  0.00  0.00 
Model #2  0.22  0.14  0.13  0.00  0.51  0.27  0.73  0.01  0.00 
Model #3  0.18  0.13  0.11  0.00  0.46  0.26  0.68  0.01  0.00 
As can be noted, the first two models have the same period for the fundamental mode in direction; nevertheless, in direction there is a slight difference due to the reduced dimensions of masonry walls width discovered by means of thermographic survey. In fact, the reduction of dimensions leads to a decrease of global stiffness while the total mass is almost the same. The periods of third model show a reduction of about 20% due to the greater value of equivalent elastic moduli estimated by means of homogenized texture. As already noted, for seismic action in direction the structure response is dominated by second mode, with also a significant contribution from the third mode (which is torsional); for seismic action in the response is dominated by first mode. It is also worth noting that Model #3 shows a reduced mass participating ratio for the fundamental mode in each direction, and therefore a greater number of modes should be considered in the evaluation of seismic response in order to achieve a suitable accuracy.
9 Concluding remarks
The sampling Kantorovich operators and the corresponding MATLAB algorithm for image reconstruction are very useful to enhance the quality of thermographic images of portions of masonry walls. In particular, after the processing by sampling Kantorovich algorithm, the thermographic image has higher definition with respect to the original one and therefore, it was possible to estimate the mechanical characteristics of homogeneous materials equivalent to actual masonries, taking into account the texture (i.e., the arrangement of blocks and mortar joints). These materials were used to model the behavior of a case study under seismic action. This model has been compared with others constructed by wellknow methods, using “naked eyes” survey and the mechanical parameters for materials taken from the Italian Building Code. Our method based on the processing of thermographic images by sampling Kantorovich operators enhances the quality of the model with respect to that based on visual survey only.
In particular the proposed approach allows to overcome some difficulties that arise when dealing with the vulnerability analysis of existing structures, which are: i) the knowledge of the actual geometry of the walls (in particular the identification of hidden doors and windows); ii) the identification of the actual texture of the masonry and the distribution of inclusions and mortar joints, and from this iii) the estimation of the elastic characteristics of the masonry. It is noteworthy that, for item i) the engineer has usually limited knowledge, due to the lack of documentation, while for items ii) and iii) he usually use tables proposed in technical manuals and standards which however give large bounds in order to encompass the generality of the real masonries. Instead, the use of reconstruction techniques on thermographic images coupled with homogenization permits to reduce these latter uncertainty on the estimation of the mechanical characteristics of the masonry.
References
 [1]
 [2] L. Angeloni and G. Vinti, Rate of approximation for nonlinear integral operators with applications to signal processing, Differential and Integral Equations, 18 (8), 855890, 2005.
 [3] C. Bardaro, P.L. Butzer, R.L. Stens and G. Vinti, Approximation of the Whittaker Sampling Series in terms of an Average Modulus of Smoothness covering Discontinuous Signals, J. Math. Anal. Appl., 316, 269306, 2006.
 [4] C. Bardaro, P.L. Butzer, R.L. Stens and G. Vinti, KantorovichType Generalized Sampling Series in the Setting of Orlicz Spaces, Sampling Theory in Signal and Image Processing, 6 (1), 2952, 2007.
 [5] C. Bardaro, P.L. Butzer, R.L. Stens and G. Vinti, Prediction by samples from the past with error estimates covering discontinuous signals, IEEE, Transaction on Information Theory, 56 (1), 614633, 2010.
 [6] C. Bardaro and I. Mantellini, Modular Approximation by Sequences of Nonlinear Integral Operators in MusielakOrlicz Spaces, Atti Sem. Mat. Fis. Univ. Modena, special issue dedicated to Professor Calogero Vinti, 46, 403425, 1998.
 [7] C. Bardaro and I. Mantellini, On convergence properties for a class of Kantorovich discrete operators, Numer. Funct. Anal. Optim.. 33 (4), 374396, 2012.
 [8] C. Bardaro, J. Musielak and G. Vinti, Nonlinear Integral Operators and Applications, De Gruyter Series in Nonlinear Analysis and Applications, New York, Berlin, 9, 2003.
 [9] C. Bardaro and G. Vinti, Some Inclusion Theorems for Orlicz and MusielakOrlicz Type Spaces, Annali di Matematica Pura e Applicata, 168, 189203, 1995.
 [10] C. Bardaro and G. Vinti, A general approach to the convergence theorems of generalized sampling series, Applicable Analysis, 64, 203217, 1997.
 [11] P.L. Butzer, A. Fisher and R.L. Stens, Generalized sampling approximation of multivariate signals: theory and applications, Note di Matematica, 10 (1), 173191, 1990.
 [12] P.L. Butzer and R.J. Nessel, Fourier Analysis and Approximation, I, Academic Press, New YorkLondon, 1971.

[13]
P.L. Butzer and R.L. Stens, Linear prediction by
samples from the past,
Advanced Topics in Shannon Sampling and Interpolation Theory
, (editor R.J. Marks II), SpringerVerlag, New York, 1993.  [14] N. Cavalagli, F. Cluni and V. Gusella, Evaluation of a Statistically Equivalent Periodic Unit Cell for a quasiperiodic masonry, International Journal of Solids and Structures, 50, 42264240, 2013.
 [15] F. Cluni and V. Gusella, Homogenization of nonperiodic masonry structures, International Journal of Solids and Structures, 41, 19111923, 2004.
 [16] D. Costarelli and G. Vinti, Approximation by Multivariate Generalized Sampling Kantorovich Operators in the Setting of Orlicz Spaces, Bollettino U.M.I., Special volume dedicated to Prof. Giovanni Prodi, 9 (IV), 445468, 2011.
 [17] D. Costarelli and G. Vinti, Approximation by Nonlinear Multivariate Sampling Kantorovich Type Operators and Applications to Image Processing, Numerical Functional Analysis and Optimization, 34 (8), 819844, 2013.
 [18] D. Costarelli and G. Vinti, Order of approximation for sampling Kantorovich operators, in print in: Journal of Integral Equations and Applications, 2014.
 [19] R. Gonzales and R. Woods, Digital Image Processing, PrenticeHall NJ  USA, 2002.
 [20] J. Musielak, Orlicz Spaces and Modular Spaces, SpringerVerlag, Lecture Notes in Math. 1034, 1983.
 [21] J. Musielak and W. Orlicz, On modular spaces, Studia Math., 28, 4965, 1959.
 [22] M.M. Rao and Z.D. Ren, Theory of Orlicz Spaces, Pure and Appl. Math., Marcel Dekker Inc. New YorkBaselHong Kong, 1991.
 [23] G. Vinti, A general approximation result for nonlinear integral operators and applications to signal processing, Applicable Analysis, 79, 217238, 2001.
 [24] G. Vinti and L. Zampogni, Approximation by means of nonlinear Kantorovich sampling type operators in Orlicz spaces, J. Approx. Theory, 161, 511528, 2009.
 [25] G. Vinti and L. Zampogni, A Unifying Approach to Convergence of Linear Sampling Type Operators in Orlicz Spaces, Advances in Differential Equations, 16 (56), 573600, 2011.
 [26] G. Vinti, L. Zampogni, Approximation results for a general class of Kantorovich type operators, in print in Advanced Nonlinear Studies (2014).
 [27]