Efficient Radio Map Construction Based on Low-Rank Approximation for Indoor Positioning

Fingerprint-based positioning in a wireless local area network (WLAN) environment has received much attention recently. One key issue for the positioning method is the radio map construction, which generally requires significant effort to collect enough measurements of received signal strength (RSS). Based on the observation that RSSs have high spatial correlation, we propose an efficient radio map construction method based on low-rank approximation. Different from the conventional interpolation methods, the proposed method represents the distribution of RSSs as a low-rank matrix and constructs the dense radio map from relative sparse measurements by a revised low-rank matrix completion method. To evaluate the proposed method, both simulation tests and field experiments have been conducted. The experimental results indicate that the proposed method can reduce the RSS measurements evidently. Moreover, using the constructed radio maps for positioning, the positioning accuracy is also improved.


Introduction
Recently, with the rapidly increasing location-based services (LBS), such as positioning, tracking, navigation, and location-based security, the positioning issue has been extensively studied.In an indoor environment, the wireless local area network (WLAN) is considered to be a basic infrastructure and can be easily implemented because of its ubiquitous coverage.For a WiFi-integrated mobile device, its location can be determined using the received signal strength (RSS) by many positioning methods [1,2], in which the fingerprintbased method is regarded as a potential method with acceptable positioning accuracy.
Generally, the fingerprint-based positioning method is implemented in two phases: the off-line training phase and the online positioning phase.In the off-line training phase, at the selected reference points, the RSSs from different access points (APs) are measured by a mobile device.The RSS measurements and the correspondent locations are generally formulated as the radio map, which infers the relation between the RSS distribution and the spatial locations.In the online positioning phase, an observed RSS measurement is matched to the radio map, and the location can be estimated by many proposed methods, such as -nearest neighbor algorithm (KNN) [1], kernel-based algorithm [3], or Bayesian estimation method [4].
The radio map is vital for the fingerprint-based positioning method, but the construction of a radio map is a timeconsuming work.To get complete RSS samples in the target area, the measurements should be recorded at dense reference points.However, in a complicated indoor environment, it is impractical to get RSSs anywhere because of the limitations of building structure and no admittance to some rooms.So some researchers attempt to reduce the measurements of radio map construction and estimate the unknown measurements using radio propagation model or interpolation methods.Compared with the intensive researches on the positioning methods, less attention has been paid to the fundamental work of radio map construction.How to construct the radio map efficiently is still a challenging issue for the fingerprint-based positioning.
In this paper, we propose an efficient radio map construction method based on low-rank approximation.Different from the existing interpolation methods, the proposed method represents the RSS distribution as a low-rank matrix based on two observations.For one thing the RSS data has high spatial correlation which implies low-rank property; for another, as an extension of compressive sensing theory, the low-rank approximation is a proven method in data reconstruction with high efficiency and good performance [5,6].So based on the low-rank matrix complement model [7], we construct the RSSs radio map from relatively sparse measurements.Additionally, to eliminate the interference of noises and obtain practical results, the basic low-rank complement model is revised by combining it with the signal spatial consistency, namely, a smoothing low-rank (SLR) model.The proposed method is evaluated on both simulation and field test data compared with the previous methods.The results show that the measurements of radio map construction can be significantly reduced and the positioning accuracy is also improved using the constructed radio map for positioning.The novelty of this paper is that we represent the RSS radio map construction problem as a revised low-rank matrix complement model.Additionally, the proposed optimization solution to the SLR model is another contribution of this paper.
The remainder of this paper is organized as follows.In Section 2, the related work of RSS radio map construction methods is briefly described.In Section 3, the low-rank complement theory and the details of the proposed SLR model are presented.In Section 4, the proposed radio map construction method is evaluated and the positioning experiments with the constructed radio map are implemented.The conclusions are drawn in Section 5.

Related Work
Constructing a fingerprint radio map is an essential step of positioning and has great influence on the positioning accuracy.Generally, high positioning accuracy can be obtained by a dense radio map.To achieve dense fingerprint measurements and reduce the calibration efforts, inserting some virtual data from the real measurements is a direct way.Many researchers utilize the radio propagation model to get dense measurements.For example, Xiang et al. [8] and Kuo and Tseng [9] use the radio propagation model to construct the radio map from a small number of on-site signal measurements.Due to the diversity of indoor environments, the common radio propagation model cannot model the RSS multipath fading and interference.So the interpolation methods, such as radial basis function (RBF) [10] and inverse distance weighting (IDW) [9], are used to construct a dense radio map.These methods generally estimate the missing RSS data by the interpolation of the measurements at local adjacent reference points.The measurements at local neighborhood usually have correlative relation.Therefore, the positioning accuracy of the extended radio map by interpolation usually outperforms the result of the original radio map.
Although the performance of the interpolation methods is acceptable, it still demands too much labeled data for obtaining accurate positioning results.So some new explorations are made to realize positioning with less RSS measurements.For example, considering the signal sparsity, the compressive sensing (CS) theory is used for indoor positioning systems [11].Inspired by this work, we utilize the low-rank approximation method to reduce the effort of constructing the fingerprint radio map.As an extension of the CS theory, low-rank matrix representation is considered more suitable for modeling the RSS spatial distribution because the signal sparse property is described in the inherent 2D space, while CS only describes the signal sparsity in the transformation coefficient domain.However, the low-rank matrix completion problem and its solution are intensively discussed recently [7,[12][13][14].Through proper sampling, the low-rank optimization can exactly recover the whole matrix merely from a small number of entries.So the low-rank completion method is widely used in many applications.For example, in WSNs, Zhang et al. [15] present a novel approach using the low-rank approximation to estimate the missing values in traffic matrices based on nuclear norm minimization.Cheng et al. [6] propose an efficient data collection approach based on a low-rank completion technique.Nikitaki et al. [16] present an efficient fingerprint positioning training method using matrix completion, which is the most related work to our method and can be regarded as the basic low-rank (BLR) model in our paper.However, in their method, the spatial consistence of RSSs is not considered and the solution of the matrix complement model is not discussed.Additionally, the validation of their method is only verified by simulation data.So we will continue to explore the application of low-rank completion method in RSS radio map construction and propose a revised low-rank complement model.

Basic Low-Rank Matrix Completion Model.
If the positioning area is represented as a rectangle, the regular RSS measurements at  1 ×  2 reference points will generate an  1 ×  2 matrix for each AP.If the sampling rate of RSS measurement is dense enough, the positioning will get high accuracy.However, it is impractical to get complete measurements at every reference point in the dense grid.Usually the RSS value at a small subset of the reference point can be collected, so the RSS measurement matrix is often incomplete.Here we adopt a mask operator A(⋅) to represent the subset sampling procedure: where  is the complete measurement matrix at all reference points and the element (, ) represents the RSS value of an AP at the reference point (, ). is an incomplete matrix with only a sparse set of elements having valid values, which means the RSS values at these reference points are available.For the sake of clarity, the operator A(⋅) can be specified as an  1 × 2 matrix , such that By this formulation, A() =  *  = , where * represents an element-wise product; that is, (, ) = (, )(, ).
As the RSS signal distribution has highly spatial and temporary correlation, the measurement matrix  collected at the target area is considered as having low-rank property approximately.Recently, Candes and Recht have proved that it is highly possible to recover a low-rank matrix from a subset of its entries [7].Thus, we can recover the measurement matrix by solving the following optimization problem: min rank () , This is the basic low-rank matrix completion model for fingerprint radio map construction.

Smoothing Low-Rank Matrix Approximation Model.
Although the low-rank matrix completion method is widely used for data reconstruction, the data properties, such as the distribution, consistence, and smoothness, are not well specified in practice.Using the basic low-rank model, the reconstruction error is usually little at the known reference points, but the reconstruction matrix is not consistent.However, the actual RSS signal spatial propagation is normally consistent at local area except for some noisy or singular measurements.To obtain a practical and local continuous reconstructed matrix, we propose a smoothing low-rank matrix approximation model by adding the smoothness constraint term in the basic low-rank model.So the low-rank model in (3) is revised as the following form: min rank () +  () , where () is the smooth term; that is, small value of () infers a good consistence of . is the weight to balance the low-rank term and the smooth term, which is usually determined by experiments.We define () via the diversity of matrix horizontal and vertical difference; that is, where  (6)

Model Solution Based on Singular Value Decomposition.
However, the minimization of matrix rank is not practical because it is an NP-hard problem [7].To work out this optimization problem, the nuclear norm ‖ ⋅ ‖ * is proposed to replace the matrix rank approximately [13].Yet the time complexity is still high for a large dimension matrix.In this paper, a solution based on singular value decomposition (SVD) is used to solve the low-rank matrix issue.Firstly, the measurement matrix  is decomposed by SVD technique into three matrices: where  is an  1 × Considering that the measurements received by the mobile devices are often not accurate and the matrix in our scenario does not completely have the low-rank property, we relax the condition of A(  ) =  and transform the constrained minimization model to an unconstrained model as the following form: where ‖A(  ) − ‖ 2  represents the reconstruction error at the sampling subset  with an adjusting weight .This is the final smoothing low-rank matrix completion model for fingerprint radio map construction.
To solve the minimum problem in (10), we adopt an alternant iteration algorithm for  and .Firstly,  and  are initialized randomly.Then, we fix  and optimize  by a linear least square method.After  is updated by the optimization, we fix  and let  be the optimization variable.Repeat the above alternant iteration procedure until the target

Input:
The sparse measurement matrix ; The estimated rank of the measurement matrix ; The model weights  and ; The error threshold for the iteration ; Output: The constructed measurement matrix X; (1) Initialize ,  randomly, denoted by  0 ,  0 .

Experiments Setup.
To evaluate the proposed method, we implement the fingerprint radio map construction experiment both on simulated and real data, and the constructed results are compared with IDW, RBF, and the basic lowrank (BLR) method.The positioning performance of the constructed radio map is also validated in this section.
In the simulation experiment, we assume that the RSS values at 5000 reference points are collected in a 100 m × 50 m square field with stepsize of 1m in horizontal and vertical directions.We randomly deploy 50 APs in this area.To simulate the RSS signal spatial distribution, we use the following radio propagation path loss model to simulate the signal attenuation: where   () is the receiving power at distance  from the AP,   is the transmitting power, and ( 0 ) is the average path loss value at the reference distance  0 (generally set to 1 m). is the path loss exponent, and   represents Gaussian noise with the distribution (0, ).Therefore, the RSSs can be calculated if the path loss exponent is known.We assume that transmission distance is 30 meters; in other words, if the distance between the AP and the reference point is more than 30 meters, the RSS value will be set to −100 dB.We set the path loss exponent  = 4.The field experiment is performed in a laboratory area with the dimension of 53 m × 15 m (shown in Figure 1) in the third floor of the information building in Beijing University of Technology.In this area, about 20 APs can be detected.In the off-line RSSs collection, a person carrying a mobile device moves through the area, and the RSS value and its location are recorded.Totally at 337 reference points (the black points in Figure 1), the RSS measurements are collected.In order to avoid systematic errors and obtain a reliable radio map, at each reference point we collect 10 samples for each AP.The average of the 10 measurements is recorded as the final RSS value for this reference point.To recover the whole measurements in the test area by the proposed SLR model, we randomly select a part of the 337 reference points as the known measurements, that is,  in (10).The remainder measurements at the 337 reference points are used for evaluating the constructed results.

Radio Map Construction Results.
To intuitively display the radio map construction results, the radio map is rendered as an image with the color values representing the RSS measurements.In the simulation experiment, the original simulated measurements of an AP and the constructed results From the results, it can be found that SLR and IDW yield good performance.As the noise has heavy effects on fingerprint, the radio map construction, and the positioning performance, the experiments are also implemented with different variances of noise.The result (Figure 3(b)) indicates that SLR and IDW behave robustly, while other methods provide bad results when the noise variance increases.It is denoted that the construction error is defined as the mean difference between the original and the constructed values at the reference points which are not included in the samples used for radio map construction.
In the above experiments, it is important to select suitable parameters  and  for the minimum optimization model.
In practice, the optimal parameters can be figured out from the statistics of the errors between the original and the constructed results.Concretely, given a subset of the measurements (e.g., a random subset of 20% measurements), we can construct the whole measurement matrix by the SLR method.Thus, the construction error can be computed at all measurement points.For different parameters, the experiments are implemented and the average construction errors can be computed to get optimal parameters.The statistics of the construction errors with different  in 20% sampling are shown in Figure 5.The construction errors of all points and the points in , which are used for SLR complement, are both computed.It is shown that if there is no smooth constraint, that is,  = 0, the construction errors of the points in  are minor, but the errors of all points are large.In other words, the BLR method meets perfectly the A() =  condition, but the whole construction error is not acceptable, as the local consistence of RSSs is invalid.So the function of the smoothing parameter  is to balance the two factors, and the optimal  in our experiments is set to 0.008.The optimal parameter  can be obtained in a similar  procedure.It is denoted that if the experimental environment is changed, the parameters should be adjusted to get ideal constructed results.
In the real instance, the constructed fingerprint radio map of an AP using 50% measurements of 337 reference points is shown in Figure 4.It is shown that the distribution of the RSSs in real environments is very complex compared with the simulation.From the constructed images, it can be found that the result of SLR method is continuous while the results of other methods are discontinuous.Although there are no complete 15 × 53 measurements, the constructed error can be computed by the average difference of the rest known measurements of the 337 reference points.The construction errors with different numbers of reference points are shown in Figure 6.Here, the reference points are randomly selected from the 337 points for radio map construction.It is shown that the construction error of the BLR method is much larger than others, so it is not included in Figure 6.From the results, it is concluded that the proposed SLR method provides the best performance, though the errors are bigger than those in simulated situation.To get more practical results, the construction experiments are repeated 5 times for a fixed number of reference points and the mean error of the 5 experiments is regarded as the final result for the number.For convenience, the parameter  is also set to 0.08 and the parameter  is set to 1.0 in these experiments.As the actual noise of the real measurement is unknown, the noise variation experiments are not implemented on the real data.

Positioning Results.
To further evaluate the efficiency of the constructed radio maps, these constructed results are used for positioning.Here we adopt the widely used positioning method, the kernel-based algorithm [3], in the online positioning phase.In our experiments, the estimated location is the mean location of 20 nearest neighbors in the sense of RSS value similarity.In the simulated situation, parts of measurements of the reference points are used for the radio map construction while the rest are served as the test.The positioning error is defined as the average distance between the estimated location and its real location.Figure 7  the positioning results are consistent with the radio map construction results.The positioning performance of the SLR method has robust behavior.In the real environment, the constructed fingerprint radio maps with different numbers of reference points are also used for positioning.The positioning errors are shown in Figure 8.Compared with the other methods, the result of the SLR method has the best positioning accuracy, and for the same positioning accuracy, the SLR method needs less RSS measurements.For example, the SLR method provides 1.53 m precision using 150 reference points while others require 200 reference points at least.So the SLR method yields a more efficient radio map construction solution.It is denoted that the positioning errors are also the average values of 5 times positioning results.

Conclusion
In this paper, an efficient radio map construction method is proposed based on low-rank approximation.Observing that the signal has high spatial correlation, the RSS distribution is modeled as low-rank matrix.By the low-rank completion theory, given relative sparse entries, the whole measurement matrix is recovered using an optimization method.As the basic low-rank completion model does not preserve the consistence of the radio map, a revised low-rank model is proposed to get practical results, in which the consistence is represented as the smooth term in the model.The solution of the proposed model is also discussed in detail.The proposed method is validated through both simulation and real experiments.The results of the experiments show that the proposed method significantly reduces the number of measurements for radio map construction and behaves robustly to noise interference.The constructed radio maps are also used for positioning.It is shown that the positioning accuracy is improved compared with the state-of-the-art methods.In future works, we plan to explore the radio map updating mechanisms and apply the proposed method in complicated environment.

15 mFigure 1 :
Figure 1: The experimental area in the third floor of the information building in Beijing University of Technology.The black points are reference points.

4
and the average of the path loss value ( 0 ) = −35 dB.The variance of the noise is confined into the interval [0, 16].Under this experiment setup, all of the RSS values between reference points and APs can be obtained.Totally, the radio map has 50 × 100 × 50 measurements and each AP has 100×50 measurements, which forms an instance of the measurement matrix .In the radio map construction experiments, a subset of the measurements at the 100 × 50 reference points (e.g., 20%) is randomly selected for low-rank completion.

Figure 2 :
Figure 2: The original radio map and the constructed results by different methods in the simulated experiments.
Different variances of noise

Figure 3 :Figure 4 :
Figure 3: Radio map construction errors on the simulated data.

Figure 5 :Figure 6 :
Figure 5: The construction errors of different smoothing parameters.

Figure 7 :
Figure 7: Positioning errors on the simulated data.

Figure 8 :
Figure 8: Positioning errors on the real data.
1 unitary matrix and  is an  2 ×  2 unitary matrix.Σ is an  1 × 2 diagonal matrix containing the singular values   , which are arranged in a decreasing order.Then we can factorize the matrix  as  = Σ  =   , where  = Σ 1/2 and  = Σ 1/2 .So the model in (6) is replaced by the following minimization model:min rank (  ) +  (      D  (  )     By the low-rank property of , the dimensions of  and  can be reduced to  1 ×  and  2 × , respectively, by omitting minimal singular values, where  is a preestimated value which is determined by the rank of .From the lemma in [17], if the restricted isometry property holds on A() and rank () < rank (  ), then (8) is equivalent to the following model: