# Deployment of Sensors According to Quasi-Random and Well Distributed Sequences for Nonparametric Estimation of Spatial Means of Random Fields

Ewa Skubalska-Rafajlowicz, Ewaryst Rafajlowicz

<figure><img src="image_1.png" /></figure>

To cite this version:

Ewa Skubalska-Rafajlowicz, Ewaryst Rafajlowicz. Deployment of Sensors According to Quasi-Random and Well Distributed Sequences for Nonparametric Estimation of Spatial Means of Random Fields. 26th Conference on System Modeling and Optimization (CSMO), Sep 2013, Klagenfurt, Austria. pp.306-316, 10.1007/978-3-662-45504-3_30 . hal-01286440

HAL Id: hal-01286440

https://inria.hal.science/hal-01286440

Submitted on 10 Mar 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Deployment of sensors according to
quasi-random and well distributed sequences for
nonparametric estimation of spatial means of
random fields

Ewa Skubalska-Rafajłowicz and Ewaryst Rafajłowicz

Institute of Computer Engineering control and Robotics
{ewa.rafajlowicz,ewaryst.rafajlowicz}@pwr.wroc.pl

**Abstract.** Our aim is to discuss advantages of quasi-random points (also known as uniformly distributed (UD) points [8]) and their sub-class recently proposed by the authors [17] that are well-distributed (WD) as sensors' positions in estimating the spatial mean. UD and WD sequences have many interesting properties that are useful both for wireless sensors networks (coverage an d connectivity) and for large area networks such as radiological or environment pollution monitoring stations. In opposite to most popular parameter estimation approaches, we consider a nonparametric estimator of the spatial mean. We shall prove the estimator convergence in the integrated mean square-error sense.

**Keywords:** sensor networks, quasi-random sequences, space-filling curves, nonparametric estimation

# 1 Introduction

Spatial sampling is a crucial issue for proper estimation of parameters in
spatio-temporal dynamical models [9], [18] and for estimation of spatial
fields (see [14]). For wireless sensor network (WSN) at least three require-
ments are crucial, namely, efficient energy usage, coverage and connec-
tivity. Here, we concentrate on the last two of them from the view-point
of sensors' deployment. Coverage (or information coverage) is the ability
of WSN to cover the whole area, assuming that a single sensor has the
ability of collecting information from its (usually circular) neighborhood.
The connectivity requires that wireless senors can transmit information
from one to another and finally to a sink. Both connectivity and coverage
require sensors to be evenly placed in the area.

Our aim is to discuss advantages of equidistributed (EQD) (also known
as uniformly distributed (UD) or quasi-random points [8]) and their
sub-class recently proposed by the authors [17] that are well-distributed
(WD). Furthermore, we shall prove that it is possible to construct a non-
parametric estimator of the mean of (possibly correlated) random field
that is based on observations from such points.

We refer the reader to [1] and [3] for surveys on WSN and to [5], [10],
[14], [11], [19] for an excerpt of approaches recently proposed for sensors
placement.
## 2 Equidistributed sequences – good candidates for sensors’ sites

Equidistributed sequences, also called uniformly distributed, quasi-random or quasi Monte-Carlo sequences, are well known in the theory of numerical integration (see, e.g., [8], [7]). Here, we summarize some of their basic properties, putting emphasis on those, which indicate that they are good candidates for sensors’ positions in WSN.

### 2.1 Definition

Define $I_d = [0, 1]^d$ as $d$-dimensional unit cube, which is our space for sensors’ deployment. Clearly, most WSN are considered for $d = 2$, but it can also be of interest in some applications to consider WSN in 3D space, e.g., for air pollution.

A deterministic sequence $(x_i)_{i=1}^n$ is called EQD sequence in $I_d$ iff

$$ \lim_{n \to \infty} n^{-1} \sum_{i=1}^{n} g(x_i) = \int_{I_d} g(x) dx \quad (1) $$

holds for every $g$ continuous on $I_d$.

Thus, formally, EQD sequence behave as uniformly distributed random sequences, since (1) mimics a law of large numbers. As we shall see later, EQD sequences are in some sense "more uniform" than uniformly distributed random sequences.

The well known EQD sequences include Corput, Halton, Hammersley, Korobov, Zaremba and Sobol (see, e.g., [6], [7] [8]). For our purposes, the following example is important.

**Important example** – the Weyl sequence is defined as follows:

$$ t_i = \text{fractional part}(i\theta), \quad i = 1, 2, \dots, $$

where $\theta$ is selected irrational number. E.g., quadratic irrational – $\theta = (\sqrt{5}-1)/2$ behaves quite well in practice. In general, as it is more difficult to approximate $\theta$ by rational numbers, then it is better candidate for generating the Weyl sequence.

### 2.2 Discrepancy as indicator of WSN coverage and connectivity

A discrepancy $D_n$ of EQD sequence is well known measure of its uniformity (see, e.g., [6], [7]). $D_n$ discrepancy of $(x_i)_{i=1}^n$ is defined as follows:

$$ D_n = \sup_{A \subset I_d} \left| \mu_d(A) - \frac{N_n(A)}{n} \right|, \quad (2) $$

where $A$ is any parallelepiped and supremum is taken with respect to all such $A \subset I_d$, $\mu_d(A)$ is $d$-dimensional Lebesgue measure (the area or volume) of $A$, while $N_n(A)$ is number of $x_i$'s in $A$.
As in the classical Monte-Carlo method, one can expect that for fairly spaced sensors $\mu_d(A)$ is close to the fraction of sensors in A. For interpreting (8) assume that the supremum is attained for a certain set $\hat{A}_n \subset I_d$. Then, $\hat{A}_n$ is the set that is most unevenly covered by sensors, i.e., it contains too small or too large number of sensors in comparison to other areas.

For this reason, in our opinion, $D_n$ is a good measure for connectivity and coverage of a sensors' net. It is, however, not easy to find $\hat{A}_n$ numerically and for this reason in the theory of numerical integration a simplified, but still useful, version of discrepancy, called $D_n^*$ discrepancy, is more frequently used.

Let $\Pi(x)$ be the parallelepiped in $I_d$ with vertices in $(0, \dots, 0)$ and $x$. $D_n^*$ discrepancy of $(x_i)_{i=1}^n$ is defined as follows:

$$D_n^* = \sup_{x \in I_d} \left| \mu_d(\Pi(x)) - \frac{N_n(x)}{n} \right|, \quad (3)$$

where $N_n(x)$ is the number of $x_i$'s in $\Pi(x)$.

Notice that $D_n^*$ is an analog to Kolmogorov-Smirnoff statistics for testing the uniformity of a distribution. $D_n^*$ can be efficiently calculated for a given set of points in $1D$, $2D$. Furthermore, it can be proved (see [7]) that

$$D_n^* \le D_n \le 2^d D_n^*. \quad (4)$$

For "good" known EQD sequences:

$$D_n^* = \frac{\log^d(n)}{n}.$$

This is much better than for "usual" uniformly distributed random variables, for which

$$D_n^* \sim \frac{1}{\sqrt{n}}.$$

The same order $1/\sqrt{n}$ is obtained when equidistant grid is considered in $\mathbb{R}^2$.

## 3 Basic properties of space-filling curves

Our idea is to transform a sequence of one dimensional EQD sequence by a space-filling curve (SFC) in order to obtain multidimensional EQD sequence with good properties. Note that a similar construction has been already proposed and used in the theory of numerical integration, but the transformed sequence was equidistant, i.e., $i/n$, $i = 1, 2, \dots, n$.

Below, we summarize known properties of space filling curves that are either directly used in the rest of the paper or indicate why we can expect that our construction leads to sensors' site with good coverage and connectivity.

Space-filling curve is a mapping $\Phi: [0, 1] \xrightarrow{\text{onto}} I_d$ such that

* $\Phi(t)$ a continuous function in $I_1 = [0, 1]$

* maps $I_1 = [0, 1]$ onto $I_d - d$-dim. cube.

The Hilbert, Peano and Sierpinski are well known examples of SFCs
[15]. All these curves: – preserve areas and neighbors (see SFC2, SFC 3 below),

– fill the space uniformly and as densely as desired.

– can be generated recursively and efficiently.

For the Peano, Hilbert and Sierpiński curves the following properties hold:

SFC 1) $\forall g : I_d \to R$, g - continuous

$$ \int_{I_d} g(x) dx = \int_{0}^{1} g(\Phi(t)) dt, \qquad (5) $$

where $x = [x^{(1)}, x^{(2)}, \dots, x^{(d)}]$

SFC 2) Hölder continuity:

$$ \exists_{C_{\Phi}>0} ||\Phi(t) - \Phi(t')|| \leq C_{\Phi}|t - t'|^{1/d}, \qquad (6) $$

where $||.||$ is the Euclidean norm in $R^d$.

It is well known that $\Phi$ does not have the inverse. However, for any Borel set $A \subset I_d$ we can define $\Phi^{-1}(A)$ as preimage of $A$.

SFC 3) $\Phi$ preserves the Lebesgue measure in the sense that

$$ \forall A \subset I_d, \mu_d(A) = \mu_1(\Phi^{-1}(A)), \qquad (7) $$

where $\mu_1, \mu_d$ are the Lebesgue measures in $R_1, R_d$, respectively. In other words, volume of Borel set $A$ in $I^d$ is equal to the length of $\Phi^{-1}(A)$ in $I_1$.

We stress that the Peano, Hilbert and Sierpiński curves can be approximated efficiently, using $O(\lceil \frac{d}{\epsilon} \rceil)$ arithmetic operations, where $\epsilon > 0$ is the approximation accuracy (see [2], [16]). Notice that for our purpose we shall calculate $\Phi(t)$ once for $t = t_1, t_2, \dots, t_n$. Approximation of the Hilbert SFC is shown in Fig. 1.

# 4 Proposed EQD sequences

The following algorithm generates EQD sequences that can be used as sensors' positions.

## Algorithm 1

Step 1) Generate EQD sequence in [0, 1] by the Weyl method: $t_i = \text{frac}(i\theta)$, $i = 1, 2, \dots, n$, where $\theta$ is irrational number, e.g., $\theta = (\sqrt{5} - 1)/2$.

Step 2*) Sort $t_i$'s and get $t_{(1)} < t_{(2)} < \dots < t_{(n)}$. This step is optional, it serves mainly for theoretical purposes, but it can also be used to determine ordering of senors along SFC, which – in turn – can be used to indicate the ordering of information transmission between sensors.

Step 3) Select SFC and generate $x_i$'s as follows: $x_i = \Phi(t_{(i)})$, $i = 1, 2, \dots, n$.

Below we state properties of the above generated sequences that justify their usefulness as sensors' positions.

Algorithm 1 generates sequences that are *extendable*, i.e., one can add points without recalculating positions of earlier sites. This property is
<figure><img src="image_3.png" /><figcaption>Fig. 1. Approximation of the Hilbert SFC using 1024 points.</figcaption></figure>

not shared by many other methods of generating multidimensional EQD sequences.

A deterministic sequence $(x_i)_{i=1}^n$ is called **well distributed** (WD) in $I_d$ iff

$$ \lim_{n \to \infty} n^{-1} \sum_{i=p}^{p+n} g(x_i) = \int_{I_d} g(x) dx $$

holds uniformly in $p$, for every $g \in C(I_d)$, i.e., continuous on $I_d$.

It is known that Weyl seq. $t_i = frac(i\theta)$ is WD. As far as we know, multidimensional, extendable WD sequences are not known. The only exception is the above proposed sequence.

**Theorem:** If $\theta$ irrational and SFC has the property SFC 1), then $x_i$'s generated by Algorithm 1 are not only EQD but also well distributed.

**Proof – EQD property.** $\forall g \in C(I_d)$ we have:

$$ n^{-1} \sum_{i=1}^{n} g(x_i) = n^{-1} \sum_{i=1}^{n} g(\Phi(t_i)) \to \\ \to \int_{0}^{1} g(\Phi(t))dt = \int_{I_d} g(x) dx, \qquad (8) $$

since $\{t_i\}_{i=1}^n$ are EQD, $g(\Phi(\cdot))$ is also continuous, while last equality follows from SFC 1). The proof of WD property uses the Weyl criterion, generalized to WD that is too complicated to be presented here (see [17] for details).

**Corollary 1.**

Under the same assumptions as in Thm. 1, for points generated by Algorithm 1 we have:

$$ D_n \to 0 \quad \text{and} \quad D_n^* \to 0 \quad \text{as} \quad n \to \infty \qquad (9) $$
The first statement follows from (8) by selecting $g$ as indicator functions of parallelepipeds. The second one is a direct consequence of (4).

Points generated by Algorithm 1, using the Hilbert SFC and $\theta = (\sqrt{5} - 1)/2$ are shown in Fig. 2. We do not have a closed form formula for

<figure><img src="image_4.png" /><figcaption>**Fig. 2.** $n = 256$ points generated by Algorithm 1, using the Hilbert SFC, and linked by thin lines, indicating their ordering on SFC that can be used as ordering for transmitting information between sensors.</figcaption></figure>

$D_n^*$ of sequences generated by Algorithm 1. Instead, we have performed extensive simulations for $n$ ranging from 50 to 30 000. Then, values of $D_n^*$ were calculated and the least-squares method was used to fit the dependence of $D_n^*(d)$ on $n$ for $d = 2$ and $d = 3$. The results are the following:

$$0.2 \log^2(n)/n \quad \text{for } d=2, \qquad (10)$$

$$0.06 \log^3(n)/n \quad \text{for } d=3. \qquad (11)$$

They follow a general pattern of "good" EQD sequences

$$D_n^*(d) = O\left(\frac{\log^d(n)}{n}\right). \qquad (12)$$

# 5 Nonparametric estimation of spatial means of random fields

Our aim in this section is to show how to estimate the spatial mean of a random field using WD points generated by Algorithm 1. The proposed
estimator is similar to those proposed earlier in [12] and [13], but it differs
in the following two respects.

- In [12] and [13] points were generated using the Halton-Hammersley sequences, which are EQD, but it is not known whether they are WD or not, while our sensors' positions are WD.

- Here, we allow correlated observations.

Let us assume that observations $y_i$ of a scalar random field with unknown mean $f(x)$ are collected at spatial points $x_i$'s, generated by Algorithm 1. More precisely,

$$y_i = f(x_i) + \epsilon_i, \quad i = 1, 2, \dots, n, \tag{13}$$

where $\epsilon_i$ are random variables (r.v.'s) for which the following conditions hold:

ERR1) $\epsilon_i$'s have zero mean, finite variance $\sigma^2$,

ERR2) for each $i = 1, 2, \dots, n$ the covariance $E(\epsilon_i \epsilon_j)$ is not equal to zero only for a finite number, $0 \le G(n) < n$ of indices $j = 1, 2, \dots, n$.

This assumption means that we allow for correlations between observations from $i$-th sensor and observations from only a finite number $G(n)$ of its nearest neighborhood sensors with indices in a set denoted by $O(i)$. As we shall see later, we may allow $G(n)$ to grow with $n$, but rather slowly. It is a reasonable assumption since the statistical dependence between sensors fades out rapidly as a function of a distance between them. Clearly, if errors are uncorrelated, we have $O(i) = \{i\}$.

Our aim is to estimate $f$ from $(x_i, y_i)$, $i = 1, 2, \dots, n$ in a nonparametric way, i.e, without imposing any finite parametrization of $f$. Instead, we impose some smoothness assumptions on $f$ and our aim is to construct estimator $\hat{f}_n(x)$ such that the integrated mean-square error (IMSE) $I(f, f_n) \stackrel{\text{def}}{=} \int_{I_d} E(f(x) - \hat{f}_n(x))^2 dx \to \infty$ as $n \to \infty$.

We select the following estimator $\hat{f}_n$ of $f$.

Algorithm 2

$$\hat{f}_n(x) = \sum_{k=1}^{N} \hat{a}_{kn} v_k(x), \tag{14}$$

$$\hat{a}_{kn} = \frac{1}{n} \sum_{i=1}^{n} y_i v_k(x_i), \quad k = 1, 2, \dots, N, \tag{15}$$

where $v_1, v_2, \dots$ is a complete sequence of orthonormal functions<fnref n="1" /> in the space $L_2(I_d)$ of square integrable functions on $I_d$. Denote by $a_k = \int_{I_d} f(x) v_k(x) dx$, $k = 1, 2, \dots$, the coefficients of the Fourier series of $f$ in the basis $(v_k)$, i.e.,

$$f(x) \sim \sum_{k=1}^{\infty} a_k v_k(x), \tag{16}$$

where convergence is understood in $L_2$ norm. In (14) one can recognize the truncated version of (16). In practice, the truncation point $N$ is also estimated. In asymptotical considerations below, we admit a slowly

<footnote n="1">¹ Orthonormality means that $\int v_k(x)^2 dx = 1$ and $\int v_k(x)v_\ell(x) dx = 0$, for all $k \ne \ell$</footnote>
growing dependence of $N$ on the number of sensors $n$ and we shall write $N(n)$.

**Remark 1.**

Notice that $\hat{a}_{kn}$ is the estimator of $a_k$, which is asymptotically unbiased. Indeed, equidistribution of $x_i$'s implies

$$E(\hat{a}_{kn}) = n^{-1} \sum_{i=1}^{n} f(x_i) v_k(x_i) \rightarrow \int_{I_d} f v_k = a_k. \quad (17)$$

But, due to WD property of $x_i$'s, we can say more, namely

$$E(\hat{a}_{kn}) = n^{-1} \sum_{i=p}^{n+p} f(x_i) v_k(x_i) \rightarrow a_k, \quad (18)$$

uniformly w.r.t $p$. This property seems to be important for WSN, because it frequently happens that a group of sensors fails (due to, e.g., a battery or communication faults) and can be replaced by differently located sensors. From (18) it follows that WSN based on sensors placed at WD points still can provide estimators with a small bias, provided that $n$, i.e., the number of active sensors, is sufficiently large.

**Remark 2.** Notice that (15) has the form that is well suited for collecting data along a network, because coefficients can be calculated recursively. Indeed, $v_k(x_i)$'s, ($k = 1, 2, \dots, N$) can be pre-computed and stored in $i$-th sensor. Sensors can be ordered along SFC (see Fig. 2) and values of partial sums, after adding $y_i v_k(x_i)$ to previous ones can be passed from one sensor to another. The role of the sink is to calculate (14).

Now, our aim is to sketch the proof of IMSE consistency of $\hat{f}_n$. For simplicity of formulas, we assume that $d = 2$ and ONS) $v_k(x^{(1)}, x^{(2)})$ are ordered and normalized trigonometric functions of the form: $1, \sin(x^{(1)}), \sin(x^{(2)}), \sin(x^{(1)}) \cos(x^{(2)}), \cos(x^{(1)}), \sin(x^{(2)}), \cos(x^{(1)}), \cos(x^{(2)}), \dots$. We have $N = N(n)$ such products in common and they are commonly bounded by $H > 0$, say, while their derivatives are commonly bounded by $N(n)H$. Notice that $N(n)$ implicitly depends on $d$, so for $d = 2$, $N(n) = M_1(n)M_2(n)$, where $M_1(n)$ and $M_2(n)$ are the numbers of trigonometric functions in $x^{(1)}$ and $x^{(2)}$ variables, respectively. Below, when we write $N(n) \rightarrow \infty$, we require that $M_1(n) \rightarrow \infty$ and $M_2(n) \rightarrow \infty$

One can weaken this assumption by allowing $H$ be dependent on $k$, as it is, e. g., for the Legendre polynomials.

The orthonormality of $(v_k)$ implies

$$IMSE(\hat{f}_n, f) = W_n + B_n^2 + R(N, f) \quad (19)$$

where $R(N, f) \stackrel{def}{=} \sum_{k=N+1}^{\infty} a_k^2$,

$$W_n \stackrel{def}{=} \sum_{k=1}^{N} \mathrm{Var}(\hat{a}_{kn}), \quad (20)$$

$$B_n^2 \stackrel{def}{=} \sum_{k=1}^{N} (E(\hat{a}_{kn}) - a_k)^2, \quad (21)$$
To estimate the variance, let us note that

$$ \mathrm{Var}(\hat{a}_{kn}) = \frac{1}{n^2} E \left( \sum_{i=1}^{n} \varepsilon_i v_k(x_i) \right)^2 = \quad (22) $$

$$ = \frac{1}{n^2} \sum_{i=1}^{n} \sum_{j=1}^{n} \mathrm{cov}(\varepsilon_i, \varepsilon_j) v_k(x_i) v_k(x_j) = $$

$$ = \frac{1}{n^2} \sum_{i=1}^{n} \sum_{j \in O(i)} \mathrm{cov}(\varepsilon_i, \varepsilon_j) v_k(x_i) v_k(x_j) $$

Due to ONS) and ERR2) we have

$$ \mathrm{Var}(\hat{a}_{kn}) \le \frac{H G(n) \sigma^2}{n} \left[ \frac{1}{n} \sum_{i=1}^{n} |v_k(x_i)| \right] \le \frac{H^2 G(n) \sigma^2}{n}. \quad (23) $$

Thus, $W_n \le \frac{N(n) H^2 G(n) \sigma^2}{n}$.

Denote by $V(I_d)$ the space of functions having bounded variation and let $V(g)$ denote the total variation of $g \in V(I_d)$. Then, by the Koksma-Hlavka inequality (see, e.g., [6]) we obtain for $f \in V(X) \cap C(X)$

$$ (E(\hat{a}_{kn}) - a_k)^2 \le (V(f \cdot v_k) \cdot D_n^*)^2 \le c_1 N^2(n) (D_n^*)^2, \quad (24) $$

where $c_1 > 0$ is a constant independent of $n$, while the second inequality can be obtained in a way similar to the one that was used in [12]. Thus, $B_n^2 \le c_1 N^3(n) (D_n^*)^2$ and finally, we obtain

$$ \mathrm{IMSE}(\hat{f}_n, f) \le c_1 N^3(n) (D_n^*)^2 + c_2 N(n)/n + R(N, f). \quad (25) $$

Notice that $f \in V(I_d) \cap C(I_d) \subset L^2(I_d)$, which implies $R(N, f) \to 0$ as $N(n) \to \infty$.

**Proposition:** If $f \in V(I_d) \cap C(I_d)$ and sequence $N(n)$ is selected in such a way that for $n \to \infty$

$$ N(n) \to \infty, \quad N^3(n) (D_n^*)^2 \to 0, \quad N(n)/n \to 0, \quad (26) $$

then $\mathrm{IMSE}(\hat{f}_n, f) \to 0$.

In the case of bivariate random fields, a particular choice of $N(n) = M_1(n) M_2(n)$ depends on the rate of decay $D_n^*$, which is typically of order $O(\frac{1}{n^{1-\epsilon}})$, where $\epsilon > 0$ is arbitrarily small. Thus, selecting $M_1(n) = M_2(n) = c n^{\beta/2}$, $c > 0$, $0 < \beta < 2(1 - \epsilon)/3$, we can assure that for $N(n) = M_1(n) M_2(n)$ conditions (26) hold.

## References

1. Akyildiz I. F., Su W., Sankarasubramaniam Y., Cayirci E.: Wireless sensor networks: a survey. Computer Networks 38, 393–422 (2002)

2. Butz A.R.: Alternative algorithm for Hilbert’s space-filling curve. IEEE Trans. Comput. C-20, 424–426 (1971)
3. Chong C-Y., Kumar S.P.: Sensor networks: Evolution, opportuni- ties, and challenges. Proceedings of the IEEE 91, 1247–1256 (2003)

4. Drmota M. and R.F. Tichy: Sequences, Discrepancies and Applica- tions. Springer, Berlin (1997)

5. Griffith D. A.: Statistical efficiency of model-informed geographic sampling designs. 7th International Symposium on Spatial Accu- racy Assesssm ent in Natural Resources and Environmental Sciences. Edited by M. Caetano and M. Painho, 91–98 (2006)

6. Kuipers, L., Niederreiter, H.: Uniform Distribution of Sequences. Wiley, New York (1974) Reprint, Dover, Mineola (2006)

7. Lemieux Ch.: Monte Carlo and Quasi-Monte Carlo Sampling. Springer, Berlin (2009)

8. Niederreiter H.: Random Number Generation and QuasiMonte Carlo Methods. SIAM, Philadelphia (1992)

9. Patan M.: Optimal sensor network scheduling in identification of distributed parameter systems. Lecture Notes in Control and Infor- mation Sciences, Springer-Verlag, Berlin (2012)

10. Pilz J., Spöck G.: Spatial sampling design for prediction taking account of uncertain covariance structure. 7th International Sym- posium on Spatial Accuracy Assessment in Natural Resources and Environmental Sciences. Edited by M. Caetano and M. Painho 109– 119 (2006)

11. Rafajłowicz E., Rafajłowicz W.: Sensors allocation for estimating scalar fields by wireless sensor networks. Proceedings 2010 Fifth In- ternational Conference on Broadband and Biomedical Communica- tions IB2Com, Malaga (2010)

12. Rafajłowicz E., Schwabe R.: Halton and Hammersley sequences in multivariate nonparametric regression. Statistics and Probability Letters 76, 803–812 (2006)

13. Rafajłowicz E., Schwabe R.: Equidistributed designes in nonpara- metric regression. Statistica Sinica 13, 129–142 (2003)

14. Rasch D., Pilz J., Verdooren L.R., Gebhardt A.: Optimal Experi- mental Design with R. Francis and Taylor (2011)

15. Sagan H.: Space-filling Curves. Springer-Verlag, Berlin (1994)

16. Skubalska-Rafajłowicz E.: Pattern recognition algorithm based on space-filling curves and orthogonal expansion. IEEE Transection on Information Theory 47, 1915–1927 (2001)

17. Skubalska-Rafajłowicz E., Rafajłowicz E.: Sampling multidimen- sional signals by a new class of quasi-random sequences. Multidim. Syst. Sign. Process. 23, 237253 (2012)

18. Uciński D.: Optimal Measurement Methods for Distributed Param- eter System Identification. CRC Press, London, New York (2005)

19. Uciński D.: An algorithm to configure a large-scale monitoring net- work for parameter estimation of distributed systems. in Proceedings of the European Control Conference 2007, Kos, Greece (2007)