INFORMATICAInformatica1822-88440868-49520868-4952Vilnius UniversityINFOR39510.15388/20-INFOR395Research ArticleIoT Devices Signals Processing Based on Shepard Local Approximation Operators Defined in Riesz MV-AlgebrasNojeDandan.noje@primus-tech.ro12
D. Noje received his BSc (eq.MSc) in mathematics and informatics (1996) and MSc in mathematics (1997) from University of Oradea. He was a lecturer at University of Oradea between 1996 and 2016. Since 2016 he is the head of Research Department of Primus Technologies SRL. His current research focus is on IoT, artificial inteligence, innovation and technology transfer. He has published 3 books and 25 scientific papers. He is a member of the EUSDR PA8 Working Group Innovation and Technology Transfer.
DzitacIoanioan@dzitac.ro34
I. Dzitac received his BSc (eq.MSc) in mathematics (1977) and PhD in information sciences (2002) from Babes-Bolyai University of Cluj-Napoca. In 2019 he has defended his Habilitation Thesis “Soft Computing for Decision Making” at Alexandru Ioan Cuza University of Iasi. He is a mathematics and information sciences professor at Agora University of Oradea (since 2017) and at Aurel Vlaicu University of Arad (since 2009). Also he is a PhD supervisor at University of Craiova. He was an adjunct professor at University of Chinese Academy of Sciences, Beijing, China (2013–2016) and since 2016 he is an advisory board member at Graduate School of Management of Technology, Hoseo University, South Korea. His current research interests include different aspects of artificial intelligence, soft computing and applications of fuzzy logic in economy (decision making, quantitative management, qualitative management, e-learning platforms management, etc.). He is a senior member of IEEE (since 2011), member of Publication Committee of International Academy of Information Technology and Quantitative Management (IAITQM, since 2014), etc. He was an invited/tutorial speaker and/or invited special sessions’ organizer and chairman in China (2013: Beijing, Suzhou and Chengdu, 2015: Dalian, 2016: Beijing), India (2014: Madurai, 2017: Delhi), Russia (2014: Moscow), Brazil (2015: Rio), Lithuania (2015: Druskininkai), South Korea (2016: Asan), USA (2018: Nebraska). He has published 3 books, 12 courses and materials for students, 6 proceedings and more than 80 scientific papers in journals and conferences’ proceedings (over 40 in ISI Web of Science indexed journals with over 200 citations). In Google Scholar he has over 510 citations.
PopNicolaenicpop@gmail.com5
N. Pop received his BSc (eq.MSc) in mathematics (1973) from West University of Timisoara, Romania and PhD in mathematics (1997) from Babes-Balyai University of Cluj-Napoca, Romania. In 2016 he has defended his habilitation thesis at School of Advanced Studies of the Romanian Academy, Bucharest. He is a professor emeritus Dr. at North University Center, Baia Mare, in Technical University of Cluj-Napoca, Romania, and a senior researcher at the Institute of Solid Mechanics of the Romanian Academy, Bucharest, Department of Robotics & Mechatronics. He is a specialist in variational analysis and numerical methods for frictional contact problems in elasticity and in solving partial differential equations; in methods control in dynamic systems and piecewise-smooth dynamical systems with friction, with applications in robotics; in stability, balance, optimal control and tracking of mobile robot trajectories; in advanced robotics intelligent control research and solid mechanics integrated into the versatile, intelligent, and portable robot platform. He has published many research papers and books that have an international reputation.
TarcaRadurtarca@uoradea.ro1
R. Tarca received his BSc (eq.MSc) “cum laudae” in mechanical engineering (1992) from Technical University of Cluj Napoca, and Bachelor in Economics from University of Oradea. He received his PhD title “cum laudae” in robotics (2001) from University Politehnica Timisoara. He received his professor title in robotics in 2004. From 2004 to 2016, he was the head of Mechatronics Department, and from 2008 he is a Co-Chairholder of the UNESCO chair in Information Technologies at the University of Oradea. He is the director of the “PRODUCTICA IMT” Research Centre since 2007. Currently he is the director of the council for doctoral studies at the University of Oradea. He has published over 150 papers in international journals and conferences’ proceedings and he has coordinated seven national and international grants as a project manager. His current research interests include aspects related to robotics, sensors and artificial intelligence.
The Industry 4.0 and smart city solutions are impossible to be implemented without using IoT devices. There can be several problems in acquiring data from these IoT devices, problems that can lead to missing values. Without a complete set of data, the automation of processes is not possible or is not satisfying enough. The aim of this paper is to introduce a new algorithm that can be used to fill in the missing values of signals sent by IoT devices. In order to do that, we introduce Shepard local approximation operators in Riesz MV-algebras for one variable function and we structure the set of possible values of the IoT devices signals as Riesz MV-algebra. Based on these local approximation operators we define a new algorithm and we test it to prove that it can be used to fill in the missing values of signals sent by IoT devices.
IoT devicessignal processingShepard local approximation operatorslocal approximation operatorsapproximation algorithmsRiesz MV-algebrasvectorial MV-algebrasUEFISCDI198PED /2017This research was funded by UEFISCDI Romania, grant number 198PED /2017, grant code: PN-III-P2-2.1-PED-2016-1955. With this grant, a new fiber optic Bragg grating sensor system designed to monitor the ethanol fermentation during the bioethanol and wine production was developed. Introduction
As part of the new industry revolution, the so-called Industry 4.0 www (2016), the automation of processes takes a more and more crucial role (Wollschlaeger et al., 2017). The process automation is based on near real time data collected by IoT devices. There are several challenges in developing applications to automate processes based on IoT devices (Heinis et al., 2017; Kamienski et al., 2017), one of them is to ensure sets of complete and validated data. Missing data can be caused by many things, but most of the times it is due to a malfunction of an IoT device or a communication problem between the IoT device and the processing applications. There is a continuous focus on finding new methods to fill in the missing data using various mathematical methods (Zhao and Zheng, 2017; Ruan et al., 2017; Leturiondo et al., 2017; Xu et al., 2017), methods that can be used to develop software modules that act as input validators for industrial automated control systems. In reality, the signal collected by IoT devices creates a discrete-time signal from a continuous process, called sample (Rajeshwari and Rao, 2008). The method we propose in this paper can be applied on both signals and samples. If signals are considered, in order to fill in some missing data, a sample can be considered in the temporal vicinity of the missing value and the proposed method can be applied on it. Considering this, further in the paper we will refer to signals.
In this paper, new Shepard local approximation operators are introduced in Riesz MV-algebras (Bede and Di Nola, 2004; Di Nola et al., 2003), and based on the Riesz MV-algebra structure of IoT devices signals, a new algorithm that can fill in the missing data is defined and tested to prove that it is suitable for the role it was designed. Since several kernels can be used by the Shepard local approximation operators, the most known ones will be used in numerical experiments, considering several parametrizations, in order to determine which are suitable for real applications. In order to have a comprehensive view about the performance of the proposed method, in the numerical experiments the signal-to-noise ratio (SNR) was also determined.
In Noje et al. (2003) it was proved that RGB model has the structure of vectorial MV-algebras. The same algebraic structure is used in this paper to model IoT devices signals. This leads us to the idea of further applications of the new proposed method in image processing, like image zooming or reconstruction of missing parts of images.
Materials and Methods
In 1958, multivalued algebras, shortly named MV-algebras, were introduced by Chang (1958; 1959) as the algebraic structures corresponding to the ∞-valued Lukasiewicz logic.
An MV-algebra is a structure A=A,⊕,¬,0A if and only if the following axioms are fulfilled: A,⊕,¬,0Ais an abelian monoid,¬¬x=x,x⊕¬0A=¬0A,¬¬x⊕y⊕y=¬¬y⊕x⊕x.
In a MV-algebra A, the constant 1A and the binary operations ⊙ and ⊖ can be defined as follows: 1A=¬0A,x⊙y=¬(¬x⊕¬y),x⊖y=x⊙¬y.
Also we can define a distance function d:A×A→A as follows:
d(x,y)=(x⊖y)⊕(y⊖x).
This distance, as it is defined, is a metric and plays a very important role in image and signal processing.
By introducing an additional external operation, in 2003, the concept of Vectorial MV-algebras (Noje and Bede, 2003), shortly named VMV-algebras, was defined. It is an algebraic structure that is used in image processing (Noje and Bede, 2001; Noje et al., 2003; Noje, 2002). Let consider an MV-algebra A and an external operation defined as follows:
∙:R+×A→A.
The MV algebra A is an VMV-algebra if and only if the following axioms are fulfilled: 1∙x=x,∀x∈A,(a+b)∙x=a∙x⊕b∙x,∀x∈Aand∀a,b∈R+,a∙(b∙x)⩽(a·b)∙x,∀x∈Aand∀a,b∈R+,d(a∙x,a∙y)⩽a∙d(x,y),∀x,y∈Aand∀a∈R+.
VMV-algebras inspired new algebraic structures, MV-modules and Riesz MV-algebras, structures that were introduced and studied in Bede and Di Nola (2004), Di Nola et al. (2003).
An MV-algebra A is a truncated MV-module over the unital latticeal ring (R,v) if an external operation ∙:R+×A→A is defined, such that the following properties are fulfilled for ∀α,β∈R+ and ∀x,y∈A. (α+β)∙x=α∙x⊕β∙x,α∙(x⊕y)=α∙x⊕α∙y,ifx⩽¬y,α∙(β∙x)=(α·β)∙x,ifα,β∈[0,v].
If property
v∙x=x,
is also fulfilled, then A is an unital MV-module over unital ring (R,v).
If an MV-algebra is a truncated unital module over (R,1), it is a Riesz MV-algebra.
In Di Nola et al. (2003), it was proved that in any Riesz MV-algebra the following properties are fulfilled: α∙(β∙x)⩽(α·β)∙x,d(α∙x,α∙y)⩽α∙d(x,y),0∙x=0,α∙0=0,x⩽y⇒α∙x⩽α∙y,α⩽β⇒α∙x⩽β∙x,α∙(x⊕y)⩽α∙x⊕α∙y, for any x,y∈A and α,β∈R+.
It was also proved that any Riesz MV-algebra is an VMV-algebra, but the reciprocal statement is not true.
If we consider a Boolean algebra A and we define α∙x=x,∀x∈A and α∈R+, it is easy to prove that all axioms of VMV-algebras are fulfilled, but since 0∙x=x the third axiom of Riesz MV-algebras is not fulfilled, thus A is not a Riesz MV-algebra. There are examples of VMV algebras that are also Riezs MV-algebras, like the RGB model (Noje and Bede, 2001; Noje et al., 2003; Noje, 2002) and the fuzzy sets.
Local approximation operators (Bittner, 2002; Lazzaro and Montefusco, 2002; Renka, 1988a; Zuppa, 2004) are used in data processing (Renka, 1988b, 1988c). Two variable Shepard local approximation operators, operators with application in image processing and similar structure data, were introduced (Shepard, 1968).
In Bede and Di Nola (2004) it was proved that Riesz MV-algebras are algebraic and topological structures for data processing, because any method developed in the classical numerical analysis is applicable in Riesz MV-algebras if the Riesz MV-algebras operations are used.
Based on this statement, we introduce Shepard local approximation operators on Riesz MV-algebras. Let us consider a Riesz MV-algebra A, and a function f:[0,n]→A and a Shepard kernel (Shepard, 1968), which is a strictly decreasing function K:[0,1]→R+. Also, we consider the set
B(x,r)={y∈[0,n]∣x−y⩽r}.
A Shepard local approximation operator is a function S:[0,n]→A defined as follows:
S(f,x)=⊕xi∈B(x,r)K(|x−xi|r)∑xi∈B(x,r)K(|x−xi|r)∙f(xi),
where ⊕ and ∙ are the Riesz MV-algebra operations.
Considering the statement that any method developed in the classical numerical analysis is applicable in Riesz MV-algebras if the Riesz MV-algebras operations are used, it is easy to see that all properties of local Shepard approximation operators hold.
Results
For industrial applications, the signals received from IoT devices are processed using computers. If we consider the numerical data types used to store information in computer memory, if they are stored using t bits, it means that the possible values for data are in the interval [0,2t−1]. Some of the operations are partial operations, for instance + is a partial operation. This leads us to consider that a suitable algebraic structure for these data types is Riesz MV-algebra.
Shepard Local Approximation Operators for IoT Device Signal Processing
It was proved that the structure ([0,2t−1],⊕,¬,0) is an MV-algebra (Noje and Bede, 2001), if the following definitions are used: x⊕y=defmin(2t−1,x+y),¬x=def2t−1−x,∀x,y∈[0,2t−1].
We consider the external operation ∙:R+×[0,2t−1]→[0,2t−1], defined as follows:
a∙x=defmin(2t−1,a·x),∀a∈R+ and ∀x∈[0,2t−1].
It was proved that the structure ([0,2t−1],⊕,¬,0,∙) is a vectorial MV-algebra and it is easy to see that it is also a Riesz MV-algebra (Noje, 2002).
If we use the above definition of ⊕ and ∙ operations and the formula of the Shepard local approximation operator from Definition 5, we can define an algorithm that can be used to fill in the missing data of signals received from IoT devices.
In this paper we consider the most known kernels: K(u)=1u2λ,Shepard kernel,K(u)=e−λu2,Exponential kernel,K(u)=sinqπusinπu2λ,Shepard–Jackson kernel, where λ is a parameter that can influence the performance of obtained results, and q is the degree of the Shepard–Jackson kernel. Several other types of kernels can be considered (Jun-Bao et al., 2014; Xiuyuan et al., 2016; Xiaodan and Bohu, 2001), but testing the method using these three is enough to prove its efficiency.
The Missing Data Fill-in Algorithm and Testing Results
We consider that the data sent by an IoT device on regular time intervals is a time-based function
f:[0,n]→[0,2t−1].
We also consider that 0 is the moment when the first signal was sent, and that 1 is the length of the time interval when a new signal is transmitted by an IoT device. In the formula of Definition 5, xi are the time intervals when new data was received and x are the time intervals when no data was received and that has to be approximated.
The fill-in algorithm has the following steps:
A kernel has to be selected;
The parameter λ is set;
If Shepard–Jackson kernel if used, the degree q of the kernel is set;
The radius r, that influence how many received values are considered in the approximation of missing values, is set;
All missing values are approximated.
In order to be able to determine the dependency of the algorithm accuracy depending on the kernel used and on parametrization, several experiments were performed. In each experiment all kernels and several random parametrizations were considered. We selected f(x)=xsin(πx20) as test function and a set of 30 signals that should be received. Several values were supposed to be missing and were approximated. To reduce the code complexity, we considered a situation where each z number of values is missing. A continuous function was considered, since one of the purposes of this algorithm is to be used to approximate missing values collected by IoT sensor system designed to monitor the ethanol fermentation during the bioethanol and wine production.
The aim of this paper is to determine which kernel of the three considered is producing better results. The performance comparison of the new method introduced in this paper using the three considered kernels and several parametrizations has been done using the formula
P=∑x∈[0,n]|f(x)−S(f,x)|,
because we are interested in the overall cumulated error.
We can raise the question: what results does this new method produce compared to other existing methods? For this reason we decided to also calculate the Signal-to-noise ratio (SNR) (Johnson, 2019; González and Woods, 2008), but we will not go deeper in its analysis, this being the target of a later work.
SNR=∑x∈[0,n]S(f,x)2∑x∈[0,n]f(x)−S(f,x)2.
Numerical Results
In this experiment we consider two situations: one is when each second value is missing, and the other is when each third value is received. Also, several parametrizations are considered. After running the tests, we get the approximation errors listed as follows (Tables 1, 2, 3, 4).
Approximation errors when each second value is missing using formula (P).
Parametrization
Shepard kernel
Exponential kernel
Shepard–Jackson kernel
r=5, λ=2, q=15
4.33074
17.7671
90.8879
r=3, λ=10, q=15
3.94205
3.9453
33.3609
r=5, λ=10, q=15
3.94205
4.87258
90.8879
r=10, λ=10, q=15
3.94205
15.685
93.3031
r=10, λ=10, q=60
3.94205
15.685
233.883
Approximation errors when each third value is received using using formula (P).
Parametrization
Shepard kernel
Exponential kernel
Shepard–Jackson kernel
r=5, λ=2, q=15
33.5505
26.5892
181.714
r=3, λ=10, q=15
38.5835
35.4093
69.8386
r=5, λ=10, q=15
38.5835
17.8535
182.038
r=10, λ=10, q=15
38.5835
23.5176
227.683
r=10, λ=10, q=60
38.5835
23.5176
296.172
Approximation errors when each second value is missing using formula (SNR).
Parametrization
Shepard kernel
Exponential kernel
Shepard–Jackson kernel
r=5, λ=2, q=15
861.081277
89.34913291
2.630481509
r=3, λ=10, q=15
958.1572323
957.2236846
18.79983289
r=5, λ=10, q=15
958.1572323
739.1964225
2.630481509
r=10, λ=10, q=15
958.1572323
113.0449305
3.490112323
r=10, λ=10, q=60
958.1572323
113.0449305
0.683421307
Approximation errors when each third value is received using using formula (SNR).
Parametrization
Shepard kernel
Exponential kernel
Shepard–Jackson kernel
r=5, λ=2, q=15
47.08291553
43.50653079
1.39908321
r=3, λ=10, q=15
37.82563498
43.90402941
11.7046936
r=5, λ=10, q=15
37.82563498
106.9724316
1.393220187
r=10, λ=10, q=15
37.82563498
53.8759769
0.769694034
r=10, λ=10, q=60
37.82563498
53.8759769
0.574528505
As we can see, by using Shepard and exponential kernels we get the best results. This result leads us to further consider only the usage of Shepard and exponential kernels as suitable for applications. In what follows, we present the pattern of approximated values printed over the original function for two different situations.
The Pattern of Approximated Values in the Case when Each Second Value is Missing. Parametrization: <inline-formula id="j_infor395_ineq_095">
<alternatives><mml:math><mml:mi mathvariant="italic">r</mml:mi><mml:mo>=</mml:mo><mml:mn>3</mml:mn></mml:math><tex-math><![CDATA[$
r=3$]]></tex-math></alternatives></inline-formula>, <inline-formula id="j_infor395_ineq_096">
<alternatives><mml:math><mml:mi mathvariant="italic">λ</mml:mi><mml:mo>=</mml:mo><mml:mn>10</mml:mn></mml:math><tex-math><![CDATA[$
\lambda =10$]]></tex-math></alternatives></inline-formula> (Fig. <xref rid="j_infor395_fig_001">1</xref>)
In the considered example, the approximated values generated using the Shepard and exponential kernels deliver very similar results considering both the error of approximation and shapes of the original and the approximated functions. It has to be mentioned that we get a very small advantage by using the exponential kernel.
The pattern of approximated values in the case when each second value is missing. For r=3, λ=10 using Shepard/exponential kernel.
The Pattern of Approximated Values in the Case when Each Third Value is Received. Parametrization: <inline-formula id="j_infor395_ineq_099">
<alternatives><mml:math><mml:mi mathvariant="italic">r</mml:mi><mml:mo>=</mml:mo><mml:mn>5</mml:mn></mml:math><tex-math><![CDATA[$
r=5$]]></tex-math></alternatives></inline-formula>, <inline-formula id="j_infor395_ineq_100">
<alternatives><mml:math><mml:mi mathvariant="italic">λ</mml:mi><mml:mo>=</mml:mo><mml:mn>10</mml:mn></mml:math><tex-math><![CDATA[$
\lambda =10$]]></tex-math></alternatives></inline-formula> (Fig. <xref rid="j_infor395_fig_002">2</xref>)
In this second example, the approximated values generated using the Shepard and exponential kernels deliver very different results both considering the error of approximation and the shapes of the original and the approximated functions. The shape of the function obtained using the exponential kernel fits much better to the shape of the original function than the shape of the function obtained using the Shepard kernel. Considering this, we can state that by using the exponential kernel we get a very clear advantage, especially when a large amount of data is missing.
The pattern of approximated values in the case when each third value is received. For r=5, λ=10 using Shepard/exponential kernel.
Conclusion
As mentioned above, we further consider only the usage of Shepard and exponential kernels for industrial applications. The results are influenced the most by value of λ when Shepard kernel is used, but there does not exist a clear dependence of approximation error on parametrization when exponential kernel is used. Thus, methods should be further developed to determine the proper set of parameters for each of the kernels. This parametrization may depend also on the shape of the function that has to be approximated and on the volume of missing data. In this optimization process, other error measuring methods can be considered, depending on the real industrial process that is to be modelled.
Another research direction is to introduce Shepard local approximation operators to approximate two-dimension functions, and a more general case for multi-dimension functions, since in most of the cases, a value of a production system is influenced by several parameters, not only by one (Noje et al., 2019).
The structure ([0,2t−1],⊕,¬,0,∙) that models the model of IoT devices signals is also the structure that models the RGB model (Noje et al., 2003). This leads us to the idea of further applications of the new proposed method in image processing, like image zooming or reconstruction of missing parts of images.
ReferencesBede, B., Di Nola, A. (2004). Elementary calculus in Riesz MV-algebras. International Journal of Approximate Reasoning, 36, 129–149.Bittner, K. (2002). Direct and inverse approximation theorems for local trigonometric bases. Journal of Approximation Theory, 117, 74–102.Chang, C.C. (1958). Algebraic analysis of many valued logics. Transactions of the American Mathematical Society, 88, 467–490.Chang, C.C. (1959). A new proof of the completeness of the Lukasiewicz axioms. Transactions of the American Mathematical Society, 93, 74–80.Di Nola, A., Flondor, P., Leustean, I. (2003). MV-modules. Journal of Algebra, 261, 21–40.González, R.C., Woods, R.E. (2008). Digital Image Processing. Prentice Hall.Heinis, T., Martinho, C.G., Meboldt, M. (2017). Fundamental challenges in developing Internet of things applications for engineers and product designers. In: Conference: 21st International Conference on Engineering Design, ICED17, pp. 279–288.Johnson, D.H. (2019). Signal-to-noise ratio. https://doi.org/10.4249/scholarpedia.2088. Available online: http://www.scholarpedia.org/article/Signal-to-noise_ratio (accessed on 1 July 2019).Jun-Bao, L., Shu-Chuan, C., Jeng-Shyang, P. (2014). Kernel Learning Algorithms for Face Recognition. Springer, New York.Kamienski, C., Jentsch, M., Eisenhauer, M., Kiljander, J., Ferrera, E., Rosengrene, P., Thestrup, J., Souto, E., Andrade, W.S., Sadok, D. (2017). Application development for the Internet of things: a context-aware mixed criticality systems development platform. Computer Communications, 104, 1–16.Lazzaro, D., Montefusco, L.B. (2002). Radial basis functions for the multi-variate interpolation of large scattered data sets. Journal of Computational and Applied Mathematics, 140, 521–536.Leturiondo, U., Salgado, O., Cianic, L., Galarb, D., Catelanic, M. (2017). Architecture for hybrid modelling and its application to diagnosis and prognosis with missing data. Measurement, 108, 152–162.Noje, D. (2002). Using Bernstein polynomials for image zooming. In: Proceedings of the Symposium Zilele Academice Clujene, Computer Science Section, pp. 99–102.Noje, D., Bede, B. (2001). The MV-algebra structure of RGB model. Studia Universitatis Babes-Bolyai, Informatica, 56(1), 77–86.Noje, D., Bede, B. (2003). Vectorial MV-algebras. Soft Computing, 7(4), 258–262.Noje, D., Bede, B., Kos, V. (2003). Image contrast modifiers using vectorial MV-algebras. In: Proceedings of the 11th Conference on Applied and Industrial Mathematics Vol. 2, pp. 32–35.Noje, D., Tarca, R., Dzitac, I., Pop, N. (2019). IoT devices signals processing based on multi-dimensional shepard local approximation operators in Riesz MV-algebras. International Journal of Computers Communications & Control, 14(1), 56–62.Rajeshwari, R., Rao, B.V. (2008). Signals and Systems. PHI Learning Pvt. Ltd.Renka, R.J. (1988a). Multivariate interpolation of large sets of scattered data. ACM Transactions on Mathematical Software, 14(2), 139–148.Renka, R.J. (1988b). Algorithm 660 QSHEP2D: quadratic shepard method for bivariate interpolation of scattered data. ACM Transactions on Mathematical Software, 14(2), 149–150.Renka, R.J. (1988c). Algorithm 661 QSHEP3D: quadratic shepard method for trivariate interpolation of scattered data. ACM Transactions on Mathematical Software, 14(2), 151–152.Ruan, W., Xu, P., Sheng, Q.Z., Falkner, N.J.G., Li, X., Zhang, W.E. (2017). Recovering Missing Values from Corrupted Spatio-Temporal Sensory Data via Robust Low-Rank Tensor Completion, Database Systems for Advanced Applications, DASFAA, 2017. Lecture Notes in Computer Science, Vol. 10177. Springer, Cham.Shepard, D.D. (1968). A two dimensional interpolation function for irregularly spaced data. In: Proceedings of 23rd Natiobal Conference ACM, pp. 517–524.Wollschlaeger, M., Sauter, T., Jasperneite, J. (2017). The future of industrial communication: automation networks in the era of the Internet of things and industry 4.0. IEEE Industrial Electronics Magazine, 11(1), 17–27.Xiaodan, X., Bohu, L. (2001). A unified framework of multiple kernels learning for hyperspectral remote sensing big data. Journal of Information Hiding and Multimedia Signal Processing, 7(2), 296–303.Xiuyuan, C., Xiyuan, P., Jun-Bao, L., Yu, P. (2016). Overview of deep kernel learning based techniques and applications. Journal of Network Intelligence, 1(3), 82–97.Xu, P., Ruan, W., Sheng, Q.Z., Gu, T., Yao, L. (2017). Interpolating the missing values for multi-dimensional spatial-temporal sensor data: a tensor SVD approach. In: Proceedings of the 14th EAI International Conference on Mobile and Ubiquitous Systems: Computing, Networking and Services, Melbourne, VIC, Australia, November 7–10, 2017(MobiQuitous 2017).Zhao, L., Zheng, F. (2017). Missing data reconstruction using adaptively updated dictionary in wireless sensor networks. In: Proceedings of Science CENet, 040.Zuppa, C. (2004). Error estimates for modified local Shepard’s interpolation formula. Applied Numerical Mathematics, 49, 245–259.www: Plattform Industrie 4.0 (2016). Aspects of the research roadmap in application scenarios. Federal Ministry for Economic Affairs and Energy, Berlin, Germany. Available online: https://www.plattform-I40.de/I40/Redaktion/eN/Downloads/Publikation/aspects-of-the-research-roadmap.html (accessed on 1 October 2018).