# Fingerprinting Indoor Positioning Method Based on Kernel Ridge Regression with Feature Reduction

Authors: Yanfen Le; Shijialuo Jin; Hena Zhang; Weibin Shi; Heng Yao
Journal: Wireless Communications and Mobile Computing
 (2021)

---

## Abstract

An important goal of indoor positioning systems is to improve positioning accuracy as well as reduce power consumption. In this paper, we propose an indoor positioning method based on the received signal strength (RSS) fingerprint. The proposed method used a certain criterion to select fixed access points (FPs) in an offline phase instead of an online phase for location estimation. Principal component analysis (PCA) was applied to reduce the features of the RSS measurements but retain the most information possible for establishing the positioning model. Then, a kernel-based ridge regression method was used to obtain the nonlinear relationship between the principal components of the RSS measures and the position of the target. We thoroughly investigated the performance of the proposed method in realistic wireless local area network (WLAN) and wireless sensor network (WSN) indoor environments and made comparisons with recently developed methods. The experimental results indicated that the proposed method was less dependent on the density of the reference points and had higher positioning accuracy than the commonly used positioning methods, and it adapts to different application environments.


---

## Body

## 1. Introduction

Recent advances in information science have made it practical and accessible to provide indoor positioning services (IPS), for applications such as indoor personal navigation, healthcare, and environmental monitoring [1–3]. The key task for such systems is to determine the position of the user or portable device. Accurate positioning in complicated indoor environments is quite challenging and has received wide attention in the past decade [4–6].A variety of wireless technologies have been used for indoor positioning. Some of them achieve high accuracy in the order of tens of centimeters, which usually requires additional specialized hardware [7–10]. Techniques based on radio frequency (RF) signals are considered as the most common cost-effective solutions for indoor positioning. Among them, the fingerprinting technique, which makes use of only received signal strength (RSS) values, has gained significant interest [11–14]. Generally, a predeployment site survey process or offline phase is required during which a radio map is constructed by collecting RSS samples (fingerprint) from location fixed access points (FP), which are access points in a WLAN or anchor nodes in a WSN at different reference points (RP) over the whole positioning area. Then, the target is located during the online phase by matching the online RSS values with the prestored fingerprints. In general, such methods require neither the location of the FPs nor a wireless channel propagation model. The fingerprinting method was first proposed in [15], where the Euclidean distance between online and offline RSS values was used to select the matched RPs and then weights the positions of the RPs to obtain the estimated target position; this technique is also known as weighted K-nearest neighbor (WKNN) algorithm. The concept of the fingerprinting method is based on the fact that the RSS value at any given point is determined primarily by the surrounding environment and the location of the FPs and is, thus, unique. However, some shortcomings of RSS measurement, such as its vulnerability to environmental changes and random fluctuations, prevent the positioning accuracy from reaching a satisfactory level [16].Various pattern recognition algorithms have been applied for better positioning performance [17, 18]. One aim of these algorithms is to improve the positioning accuracy even if there are anomalous or missing RSS measurements. Classification algorithms such as support vector machine (SVM) [1, 19, 20], decision tree [21], or even neural networks classifier [22] have been used for indoor localization. In general, the number of classifier models increases with the number of RPs when indoor localization is considered as a multiclass classification problem, which incurs high time and memory consumption. Some models use regression algorithms with features extracted from the RSS values to build a robust and adaptive relationship between the positions and the RSS measures. The sparse recovery algorithm, the least absolute shrinkage and selection operator (LASSO) algorithm, kernel ridge regression (KRR), and Elastic-Net algorithm have been used to model linear or nonlinear relations for better performance [23–26].Feature extraction is a commonly used method for data reduction and noise elimination [27, 28]. The new method we proposed uses feature reduction before a relationship is established between the locations and the RSS values. The proposed technique uses principal component analysis (PCA) to find a new presentation of data in terms of least square [29, 30]. Thus, a new set of variables called principal components (PCs) instead of RSS values is obtained. The advantages of using PCA to reduce features are as follows: first, PCs are compact descriptions of the measured RSS values. Choosing a proper number of PCs eliminates redundancy but retains maximal information, which is more conducive to establishing the positioning model. The significance of PCA also lies in the reduction of the computational complexity. In general, the dimension of RSS features, or the number of FPs, is an important factor in determining computational complexity during the entire positioning process. Instead of choosing a subset of FPs as reported in previous works [23, 24, 31], PCA selects a subset of PCs but retains information from all the FPs with low dimensionality. Then, a kernel-based model is used to explore a proper linear combination of kernels with selected PCs to represent the position. In the experiments, the localization system is evaluated thoroughly using the collected RSS data from two actual indoor environments, including a corridor with WLAN signals and a hall with WSN signals.The rest of this paper is organized as follows: the details of the proposed positioning algorithm are described in Section2. The experimental setup and a discussion of the results are introduced in Section 3. Finally, the conclusion is presented in Section 4.


## 2. Proposed Algorithm Based on PCA and KRR

This section describes the proposed fingerprinting method as shown in Figure1, which consists of two phases: the offline phase and the online phase. In the offline phase, a site survey is carried out to collect the RP’s information and the corresponding RSS values, then spatial filter based on FP coverage and FP selection according to certain criterion is used to select the proper set of FPs to build a radio map, and then the model training with the proposed PCA and KRR algorithm is performed. In the online phase, the localization process using the PCA reduction and KRR is performed to get the target position.Figure 1
Flow diagram of the proposed indoor positioning method.


### 2.1. Radio Map Construction with FP Selection

Consider an indoor environment where RF signals from FPs can be received throughout the area. A number of RPs, denoted asP1,P2,⋯,PN, are set, and their coordinates pPi=xi,yi are recorded. The RF signals are scanned at each RP, and RSS values collected from the same FP are averaged and then stored as a fingerprint into radio map in a vector Fi=fi1,⋯,fij,⋯,fiK, where fij,i=1,⋯,N,j=1,⋯,K, holds the time average RSS values from FPj at location Pi. Let χ=FP1,⋯,FPK be the set of FPs. Because not all FPs detected at the site are available at each RP, a certain value, -100 dBm for WSN-based anchors and 0 dBm for WLAN-based APs in our experiments, is experimentally set to imply an FP’s unavailability.It is also worth mentioning that in an indoor environment, especially in a WLAN-based positioning scenario, up to dozens of APs can usually be detected throughout the site. For example, in our experiments, 198 APs in total were detectable on one floor, which is far more than the number of FPs required for positioning. The signal from each FP could provide some information about the location, while too many FPs used for positioning increases the cost of storage and the computation of the algorithm. In addition, FPs having a large variance of RSS values are not suitable for positioning and could contribute to large errors. Thus, a certain criterion should be used to select a particular set of FPs that present the characteristics of the signal distribution effectively.In our work, we considered three FP section methods to find a balance between the number of FPs and the positioning accuracy they can achieve. Because the FP selection is applied in the offline phase in our method, the radio map constructed with the selected FPs is to be used for further positioning model setup. Therefore, it is necessary that FPs detected in the online phase be as consistent as possible with those selected in the offline phase. Therefore, a spatial filter based on the coverage is first applied to find the FPs with the largest coverage. We define the FP coverage as(1)Cj=∑i=1Ncij,j=1,2,⋯,K,where cij=1 if the RSS value from FPj is available at RPi and is 0 otherwise. The FPj with Cj≥C0 is selected as the basis of the subsequent FP section, where C0 is a threshold varying with the number of RPs in different scenarios. This value is set to 9 in one of our implementations, which means at least 9 RPs can receive RSS data from each filtered FP. After the spatial filter, a group of ∣H∣, H⊂χ, FPs are selected. FP selection is then carried out within this group.Three FP selection criteria, denoted as Strongest FPs, Least-variance FPs, and Combined criterion, are evaluated in our work.(1)
Strongest FPs. This criterion assigns a score to an FP according to the strength of its signal throughout the site, which is defined as:(2)ξj1=∑i=1Nfij,j∈Η,where fij is the time-average RSS value. The idea behind this criterion is that bigger RSS values can reveal more features of the site. The FPs are sorted in decreasing order according to their scores, and then a subset of ∣S1∣, S1⊆H, FPs with the largest ξ1 are selected.
(2)
Least-variance FPs. In this criterion, scores are also assigned to FPs but with focus on the variance of the RSS signals throughout the collection time. The score of FPj is calculated as follows:(3)ξj2=∑i=1NΔij,j∈Η,where Δij is the unbiased estimated variance of the RSS readings from FPj at location Pi. An FP with a low variance indicates that it can provide stable signals over time, which is conducive to a correct match between the online RSS and the fingerprints. It is important to note that some RPs cannot receive the signal from FPj as mentioned earlier; the variance of these RPs needs to be set to a large value (e.g., 100 in our implementation). The FPs are sorted in increasing order according to their scores, and then a subset of ∣S2∣, S2⊆H, FPs with the lowest scores are selected.
(3)
Combined criterion selected FPs. The criterion accounts for both the spatial distribution of RSS across all RPs and the RSS variance over the collecting time. According to this criterion, each FP is assigned a score calculated as follows:(4)ξj3=∑i=1Nfij−f¯j2∑i=1NΔij,j∈Η,where
(5)f¯j=1N∑i=1Nfij,j∈Η.A higher score means better stability of RSS over time as well as greater discriminability of FPs across the site. A subset of∣S3∣, S3⊆H, FPs with the highest ξ3 are selected.Now assume thatK′ FPs are retained after selection to build the radio map Ψ, which is denoted by:
(6)Ψ=Pi,Fi=p1F1p2F2⋯⋯pNFN=x1y1f1,1⋯f1,K′x2y2f2,1⋯f2,K′⋯⋯⋯⋯⋯xNyNfN,1⋯fN,K′,The total fingerprintF=F1,F2,⋯,FN′T can be considered prior knowledge of the signal characteristics of the surveyed site. In the online phase, the target gets a new Ft to determine its coordinate p^t, which is the key task of the indoor positioning. In this way, a radio map Ψ is crucially important for the proper training of the regression model.


### 2.2. Feature Reduction with PCA

PCA is a fast and efficient technique of dimensionality reduction used widely [32]. Instead of directly using all the selected FPs, our approach replaces the RSS values with a subset of PCs, which is obtained by a transformation based on the PCA. Given a transformation between vector F and ψ as ψ=FW, where F∈R1∗K′, ψ∈R1∗K, and k<K′. F should be reconstructed from ψWT and labeled as F^. PCA seeks to solve the problem of minimizing the mean square reconstruction error, which can be expressed by:
(7)argminW∑i=1NF^i−Fi22=argminW∑i=1NψiWT−Fi22,where W is the transformation matrix and N is the number of training samples. This optimization problem has been solved well, and W can be obtained by the following steps. First, calculate the covariance matrix of all RSS samples defined by:
(8)Cσ=∑i=1NFi−F¯Fi−F¯T,where F¯=1/N∑i=1NFi is the mean of the samples, Cσ is a positive-semidefinite symmetric matrix, and its eigenvalues are easy to compute and sort in descending order as λ1,λ2,⋯,λK′. The corresponding normalized eigenvectors w1,w2,⋯,wK′ are geometrically orthonormal and statistically uncorrelated. Then, the K′×k transformation matrix W has the form:
(9)W=w1⋯wkT.Herek is the number of PCs selected for the best transmission from RSS values to PCs with reduced dimension for maximal extraction of signal features. Then, the proper value of k is determined by cross-validation with the sampled fingerprints in the radio map during the offline phase. When the matrix W is available, the new PC-based fingerprinting Φ transformed from the RSS-based F is given by:
(10)Φ=FW=ψ1ψ2⋮ψN=φ11φ21⋮φN1φ12φ22⋮φN2⋯⋯⋮⋯φ1kφ2k⋮φNk.A regression model can be built between the newΦ and the locations.During the positioning stage, the online PCs can be directly extracted from the online RSS measurements using the trained matrixW:
(11)ψt=FtW=φt1,φt2,⋯,φtk.With the extractedΨt, the regression model can be used to determine the position of the target.


### 2.3. KRR Approach Based on PCs

A brief review is given for a general understanding of the theory of following adopted kernel ridge regression. Given a training setψ1,p1,⋯,ψN,pN, where N is the number of samples, each ψi∈ℜk is a row vector in Φ denoting an input sample with a corresponding out pi∈ℜ. The ridge regression algorithm entails solving the following optimization problem [33].
(12)ω^=argminωP−ωTP−Φω+λω2,where ω∈Rk is the weight vector, λ≥0 is a regularization parameter tuned to control the compromise between the training error and the complexity of the solution. Then, a regression model is found to describe the linear relationship between input vector ψ and output p, such that
(13)p=ωTψ,where ω=ΦTΦΦT+λIk−1P, and Ik is a k×k identity matrix.Generally, in nonlinear cases, a kernel-based method is introduced to solve the problems, which map the samples into a higher dimensional feature space where the problems become linear separable, such that(14)ϕ:ψ∈ℝk⟶ϕψ∈ℝΛ,where superscript Λ stands for a higher-dimensional space. The mapping ϕ which is chosen to convert nonlinear relation between the output and the independent input variables into linear relation is not necessary to know. Then, the regression can be constructed in the feature space, and the solution to the regression problem only depends on the dot product in the feature space. The kernel function satisfying Mercer’s condition is introduced as a format of a dot product, Kψi,ψj=ϕψi·ϕψj. Now, the weight vector can be rewritten as:
(15)ω=ΦTK+λIN−1P=ΦTα=∑i=1Nαiψi,where K is an N×N kernel matrix with the element of κψi,ψj, P is the coordinate sets of the RPs, and α=K+λIk−1P can be determined by the PC-based fingerprinting Φ.In the online phase, the target’s position can be directly estimated as(16)p^=ωTψt=∑i=1NαiψiTψt=∑i=1Nαiκψt,ψi,In our method, a Gaussian kernel is used:(17)κψi,ψj=exp−ψi−ψj22σ2,where σ is the bandwidth of the Gaussian kernel. Considering that our output is two coordinates, x^ and y^, two parameters, σx and σy, should be set and determined by the training RSS data.


## 2.1. Radio Map Construction with FP Selection

Consider an indoor environment where RF signals from FPs can be received throughout the area. A number of RPs, denoted asP1,P2,⋯,PN, are set, and their coordinates pPi=xi,yi are recorded. The RF signals are scanned at each RP, and RSS values collected from the same FP are averaged and then stored as a fingerprint into radio map in a vector Fi=fi1,⋯,fij,⋯,fiK, where fij,i=1,⋯,N,j=1,⋯,K, holds the time average RSS values from FPj at location Pi. Let χ=FP1,⋯,FPK be the set of FPs. Because not all FPs detected at the site are available at each RP, a certain value, -100 dBm for WSN-based anchors and 0 dBm for WLAN-based APs in our experiments, is experimentally set to imply an FP’s unavailability.It is also worth mentioning that in an indoor environment, especially in a WLAN-based positioning scenario, up to dozens of APs can usually be detected throughout the site. For example, in our experiments, 198 APs in total were detectable on one floor, which is far more than the number of FPs required for positioning. The signal from each FP could provide some information about the location, while too many FPs used for positioning increases the cost of storage and the computation of the algorithm. In addition, FPs having a large variance of RSS values are not suitable for positioning and could contribute to large errors. Thus, a certain criterion should be used to select a particular set of FPs that present the characteristics of the signal distribution effectively.In our work, we considered three FP section methods to find a balance between the number of FPs and the positioning accuracy they can achieve. Because the FP selection is applied in the offline phase in our method, the radio map constructed with the selected FPs is to be used for further positioning model setup. Therefore, it is necessary that FPs detected in the online phase be as consistent as possible with those selected in the offline phase. Therefore, a spatial filter based on the coverage is first applied to find the FPs with the largest coverage. We define the FP coverage as(1)Cj=∑i=1Ncij,j=1,2,⋯,K,where cij=1 if the RSS value from FPj is available at RPi and is 0 otherwise. The FPj with Cj≥C0 is selected as the basis of the subsequent FP section, where C0 is a threshold varying with the number of RPs in different scenarios. This value is set to 9 in one of our implementations, which means at least 9 RPs can receive RSS data from each filtered FP. After the spatial filter, a group of ∣H∣, H⊂χ, FPs are selected. FP selection is then carried out within this group.Three FP selection criteria, denoted as Strongest FPs, Least-variance FPs, and Combined criterion, are evaluated in our work.(1)
Strongest FPs. This criterion assigns a score to an FP according to the strength of its signal throughout the site, which is defined as:(2)ξj1=∑i=1Nfij,j∈Η,where fij is the time-average RSS value. The idea behind this criterion is that bigger RSS values can reveal more features of the site. The FPs are sorted in decreasing order according to their scores, and then a subset of ∣S1∣, S1⊆H, FPs with the largest ξ1 are selected.
(2)
Least-variance FPs. In this criterion, scores are also assigned to FPs but with focus on the variance of the RSS signals throughout the collection time. The score of FPj is calculated as follows:(3)ξj2=∑i=1NΔij,j∈Η,where Δij is the unbiased estimated variance of the RSS readings from FPj at location Pi. An FP with a low variance indicates that it can provide stable signals over time, which is conducive to a correct match between the online RSS and the fingerprints. It is important to note that some RPs cannot receive the signal from FPj as mentioned earlier; the variance of these RPs needs to be set to a large value (e.g., 100 in our implementation). The FPs are sorted in increasing order according to their scores, and then a subset of ∣S2∣, S2⊆H, FPs with the lowest scores are selected.
(3)
Combined criterion selected FPs. The criterion accounts for both the spatial distribution of RSS across all RPs and the RSS variance over the collecting time. According to this criterion, each FP is assigned a score calculated as follows:(4)ξj3=∑i=1Nfij−f¯j2∑i=1NΔij,j∈Η,where
(5)f¯j=1N∑i=1Nfij,j∈Η.A higher score means better stability of RSS over time as well as greater discriminability of FPs across the site. A subset of∣S3∣, S3⊆H, FPs with the highest ξ3 are selected.Now assume thatK′ FPs are retained after selection to build the radio map Ψ, which is denoted by:
(6)Ψ=Pi,Fi=p1F1p2F2⋯⋯pNFN=x1y1f1,1⋯f1,K′x2y2f2,1⋯f2,K′⋯⋯⋯⋯⋯xNyNfN,1⋯fN,K′,The total fingerprintF=F1,F2,⋯,FN′T can be considered prior knowledge of the signal characteristics of the surveyed site. In the online phase, the target gets a new Ft to determine its coordinate p^t, which is the key task of the indoor positioning. In this way, a radio map Ψ is crucially important for the proper training of the regression model.


## 2.2. Feature Reduction with PCA

PCA is a fast and efficient technique of dimensionality reduction used widely [32]. Instead of directly using all the selected FPs, our approach replaces the RSS values with a subset of PCs, which is obtained by a transformation based on the PCA. Given a transformation between vector F and ψ as ψ=FW, where F∈R1∗K′, ψ∈R1∗K, and k<K′. F should be reconstructed from ψWT and labeled as F^. PCA seeks to solve the problem of minimizing the mean square reconstruction error, which can be expressed by:
(7)argminW∑i=1NF^i−Fi22=argminW∑i=1NψiWT−Fi22,where W is the transformation matrix and N is the number of training samples. This optimization problem has been solved well, and W can be obtained by the following steps. First, calculate the covariance matrix of all RSS samples defined by:
(8)Cσ=∑i=1NFi−F¯Fi−F¯T,where F¯=1/N∑i=1NFi is the mean of the samples, Cσ is a positive-semidefinite symmetric matrix, and its eigenvalues are easy to compute and sort in descending order as λ1,λ2,⋯,λK′. The corresponding normalized eigenvectors w1,w2,⋯,wK′ are geometrically orthonormal and statistically uncorrelated. Then, the K′×k transformation matrix W has the form:
(9)W=w1⋯wkT.Herek is the number of PCs selected for the best transmission from RSS values to PCs with reduced dimension for maximal extraction of signal features. Then, the proper value of k is determined by cross-validation with the sampled fingerprints in the radio map during the offline phase. When the matrix W is available, the new PC-based fingerprinting Φ transformed from the RSS-based F is given by:
(10)Φ=FW=ψ1ψ2⋮ψN=φ11φ21⋮φN1φ12φ22⋮φN2⋯⋯⋮⋯φ1kφ2k⋮φNk.A regression model can be built between the newΦ and the locations.During the positioning stage, the online PCs can be directly extracted from the online RSS measurements using the trained matrixW:
(11)ψt=FtW=φt1,φt2,⋯,φtk.With the extractedΨt, the regression model can be used to determine the position of the target.


## 2.3. KRR Approach Based on PCs

A brief review is given for a general understanding of the theory of following adopted kernel ridge regression. Given a training setψ1,p1,⋯,ψN,pN, where N is the number of samples, each ψi∈ℜk is a row vector in Φ denoting an input sample with a corresponding out pi∈ℜ. The ridge regression algorithm entails solving the following optimization problem [33].
(12)ω^=argminωP−ωTP−Φω+λω2,where ω∈Rk is the weight vector, λ≥0 is a regularization parameter tuned to control the compromise between the training error and the complexity of the solution. Then, a regression model is found to describe the linear relationship between input vector ψ and output p, such that
(13)p=ωTψ,where ω=ΦTΦΦT+λIk−1P, and Ik is a k×k identity matrix.Generally, in nonlinear cases, a kernel-based method is introduced to solve the problems, which map the samples into a higher dimensional feature space where the problems become linear separable, such that(14)ϕ:ψ∈ℝk⟶ϕψ∈ℝΛ,where superscript Λ stands for a higher-dimensional space. The mapping ϕ which is chosen to convert nonlinear relation between the output and the independent input variables into linear relation is not necessary to know. Then, the regression can be constructed in the feature space, and the solution to the regression problem only depends on the dot product in the feature space. The kernel function satisfying Mercer’s condition is introduced as a format of a dot product, Kψi,ψj=ϕψi·ϕψj. Now, the weight vector can be rewritten as:
(15)ω=ΦTK+λIN−1P=ΦTα=∑i=1Nαiψi,where K is an N×N kernel matrix with the element of κψi,ψj, P is the coordinate sets of the RPs, and α=K+λIk−1P can be determined by the PC-based fingerprinting Φ.In the online phase, the target’s position can be directly estimated as(16)p^=ωTψt=∑i=1NαiψiTψt=∑i=1Nαiκψt,ψi,In our method, a Gaussian kernel is used:(17)κψi,ψj=exp−ψi−ψj22σ2,where σ is the bandwidth of the Gaussian kernel. Considering that our output is two coordinates, x^ and y^, two parameters, σx and σy, should be set and determined by the training RSS data.


## 3. Experimental Results and Analysis

In this section, we introduce the experimental setups and evaluate the performance of the proposed algorithm by comparing it with other algorithms. The performance of the proposed method is measured by the average positioning error (AE) of all test points and the empirical cumulative distribution function (CDF) of errors. The former is intended to calculate the Euclidean distance between the estimated and the actual location of the test point, while the latter indicates the maximum and minimum errors.


### 3.1. Experimental Setups

The experiments were conducted at two office sites of the Optical Engineering and Technology building at the University of Shanghai for Science and Technology, China, as shown in Figure2. The first site labeled as Map 1 includes four hallways, where 92 locations are arranged as PRs with a grid spacing of 1.8 m. The BSSID (MAC address) and RSS data of the available WLAN APs were collected with a scanning interval of 3 seconds using a HUAWEI Android cell phone within 27 seconds at each RP. A total of 198 APs, all in the rooms or from the next floors, were detected throughout the site. The RSS data of the desired APs were averaged to construct the radio map. The variances of the RSS data were also obtained. Online RSS measures were collected on different days using the same device at 42 test points. The second site labeled as Map 2 includes two hallways and one lobby, which is equipped with WSN nodes based on TI CC2430 including 90 RPs spaced 1.8 m apart and 18 anchors working as FPs. The RSS measurements were obtained at each RP by moving node with sniffer software. During the offline phase, the RSS data were recorded for a period of 120 seconds (2 samples per second) over each RP. The online observations were collected in the same way with 36 selected test points.Figure 2
Maps of experimental sites. The red dots indicate the RPs, and the blue triangles represent the test points.
(a)
Map 1(b)
Map 2Figure3 shows the number of received FPs at each RP of the two sites in the offline phase. It can be seen from Figure 2 that most of the RPs receive 10~30 FPs on Map 1, while on Map 2, this number is 10~15. No RP can receive all the FPs’ signals on both maps. In these two typical indoor environments, especially that of Map 1, the dense distribution of rooms and walls greatly attenuated or obstructed the radio signals. Although a large number of FPs could be detected in the whole area, the FPs available at each RP were greatly reduced.Figure 3
The number of received FPs at each RP in the offline phase.
(a)
Map 1(b)
Map 2


### 3.2. Experimental Results with Different FP Selections

Localization performance is always related to the number of FPs (K′) used for positioning. In our method, FPs with certain coverage areas were first selected using a spatial filter. The threshold C0 varied with the layout of the map, as well as the distributed RPs and FPs. Figure 4 shows the CDF of the positioning error when different thresholds were applied to filter the FPs. On Map 1, although there are up to 198 FPs, only 52 FPs can be received by more than 15 RPs out of a total of 92. The number K′ increases to 89 when the threshold drops to 9. With the FP selection criterion, 80 FPs were used for positioning, and almost no performance difference was observed as shown in Figure 4(a). Figure 4(b) shows the experimental results for Map 2. When C0≤50, all 18 FPs were received by all RPs. This may be attributed to the regular distribution of WSN anchor nodes and an indoor environment without many barriers. Considering this, no spatial filter was used in the following experiments for Map 2.Figure 4
Positioning error of the proposed method with different thresholds for the spatial filter.
(a)
Map 1(b)
Map 2Figure5 depicts the average location errors under the three different FP selection criteria, namely, the Strongest FPs, Least-variance FPs, and Combined criterion, when a proper subset of PCs is used for KRR-based localization.Figure 5
Average positioning error of the proposed method using three different FP selection criteria.
(a)
Map 1(b)
Map 2It can be seen that no matter which FP selection criterion is adopted, the positioning method with more selected FPs obtains higher positioning accuracy at both sites. As shown in Figure5(a), when 80 FPs are adopted, that is, 80 FPs are selected from the spatial filtered 98 FPs according to the FP selection criterion, the AE under the three criteria is all approximately 2.4 m. Almost the same result can be observed in Map 2. It is worth mentioning that even if we choose 80 FPs in Map 1 according to a certain FP selection criterion, such as the strongest FPs, only 18 of the RPs throughout the whole site can receive signals from more than 20 FPs. The situation is even worse when only 40 FPs are used; only one of the 92 RPs can receive signals from more than 20 FPs, which may explain the relationship between the positioning error and the number of FPs in Figure 5(a). For the positioning method that needs to establish the model between positions and signals in the offline phase, it is necessary to collect enough FP RSS measurements to extract features, especially in indoor environments with rooms and walls that greatly attenuate or completely block signal propagation. This is a completely different approach from those that apply FP selection in the online phase, as proposed in [23, 24, 31].


### 3.3. Analysis of the Number of PCs

The proper choice ofk (the number of PCs) is important to the proposed method. A larger value of k seems to retain more information after PCA transformation, but it may also retain redundant information and noise, which can be seen from the experimental results. Figure 6 compares the average positioning errors for the different number of FPs versus the number of PCs at both sites.Figure 6
Average positioning error versus the number of PCs. (a) FPs were selected with Least-variance FP selection. (b) FPs were selected with Combined criterion FP selection.
(a)
Map 1(b)
Map 2As illustrated in Figure6, no matter how many FPs are used in the experiments, the positioning accuracy varies with the number of PCs, and there is always an optimal value kp for the highest positioning accuracy, but this value is certainly not the maximal one. For example, in the case of the 80 FPs adopted in Map 1, the optimal kp is 23. Further research found that this value kp ensured that the selected PCs retained more than 85% of the total information if eigenvalues were used to quantify the information contribution of each PC, which provides a relatively simple way to determine this optimal value. It is also worth mentioning that no matter how many PCs are extracted, the premise is to have sufficient FPs’ information for high positioning accuracy. From Figure 6(a), the best performance in the case of 40 FPs is 2.77 m, while with 80 FPs, the worst performance is 2.56 m. The same is true for Map 2, where the corresponding values for the best and worst performance are 2.28 m and 1.88 m. It can also be inferred from Figure 5 that the number of FPs has a greater impact on performance compared with the number of PCs, which means considering enough FPs before choosing the right number of PCs.


### 3.4. Reduction in the Number of RPs

Figure7 shows the comparison of positioning error among Euclidian-WKNN, LASSO, KRR, and the proposed PCA-KRR method with respect to the number of RPs. For Map 1, a total number of 80 FPs selected with the Combined criterion were used with the optimal value of kp. The number of FPs is 18 for Map 2. From the results displayed in Figure 7, we can see that the accuracy of the positioning gradually drops for both maps as the number of RPs used for positioning decreases, while the proposed PCA-KRR method achieves significantly better accuracy for both maps than the other methods regardless of the number of RPs. For instance, when the number of PRs used was 92 and a quarter of RPs (23) in Map 1, the AE of the proposed PCA-KRR method is, respectively, 2.34 m and 2.88 m, which are lower than the 2.71 m and 3.93 m of WKNN and the 2.78 m and 3.20 m of KRR. It is worth noting that the LASSO algorithm shows a different positioning performance for each of the two maps. For Map 1, the reduction in the number of RPs does not affect the positioning accuracy. While for Map 2, the result is different. The reduction of RPs introduces a large positioning error, indicating that LASSO is not robust enough for various environments. Furthermore, from these experimental results, it can be deduced that the density of RP does affect the positioning accuracy of our method, but the influence is limited. For example, when only a quarter of the RPs are used for positioning, i.e., the space between the RPs is as large as 7.2 m, the AE of the PCA-KRR method still reaches 2.88 m and 2.13 m, respectively, in Map 1 and Map 2 when enough FPs are guaranteed. Therefore, it is suggested that the number of FPs should be as high as possible while the density of RPs can be appropriately reduced to ensure the positioning accuracy while reducing the cost of labor.Figure 7
Average positioning error versus the number of RPs of WKNN, LASSO, KRR, and the proposed method. 1/2, 1/3, and 1/4 RPs mean, respectively, a half, a third, and a quarter of the total RPs in the maps.
(a)
Map 1(b)
Map 2


### 3.5. Tuning the Parameters

There are several ways to tune the regularization parameterλ and the bandwidth of the Gaussian kernel σ involved in the proposed method. One effective approach for setting these parameters is the well-known cross-validation (CV). Another approach is to use the training data to find the best values that make the positioning algorithm the most accurate, which has been applied in this paper. We use some of the collected RSS data at the test points for parameter training. The experimental results show that the best value of λ changes only with the positioning site and does not change with the number of RPs and PCs used in the method. This value is 0.001 and 0.0001 in Map 1 and Map 2, respectively. As mentioned earlier, there are also two kernel widths σx and σy for two coordinates. Table 1 shows the positioning error for the different σ values used for positioning in the two maps. When σx and σy take the same value, there is no significant impact on the positioning accuracy, which can be observed from the results for both maps. Whether the two parameters choose the same value or not affects the positioning accuracy within 2%. In our experimental results above, σx and σy share the same value for the simplification of parameter turning.Table 1
Positioning error with different bandwidth of the Gaussian kernelσ used in the method.
SiteNumber of RPsNNumber of PCskσxσyAE/mMap 1922211112.341072.30462512122.5013112.49Map 290814141.771581.7545915151.7811151.74


## 3.1. Experimental Setups

The experiments were conducted at two office sites of the Optical Engineering and Technology building at the University of Shanghai for Science and Technology, China, as shown in Figure2. The first site labeled as Map 1 includes four hallways, where 92 locations are arranged as PRs with a grid spacing of 1.8 m. The BSSID (MAC address) and RSS data of the available WLAN APs were collected with a scanning interval of 3 seconds using a HUAWEI Android cell phone within 27 seconds at each RP. A total of 198 APs, all in the rooms or from the next floors, were detected throughout the site. The RSS data of the desired APs were averaged to construct the radio map. The variances of the RSS data were also obtained. Online RSS measures were collected on different days using the same device at 42 test points. The second site labeled as Map 2 includes two hallways and one lobby, which is equipped with WSN nodes based on TI CC2430 including 90 RPs spaced 1.8 m apart and 18 anchors working as FPs. The RSS measurements were obtained at each RP by moving node with sniffer software. During the offline phase, the RSS data were recorded for a period of 120 seconds (2 samples per second) over each RP. The online observations were collected in the same way with 36 selected test points.Figure 2
Maps of experimental sites. The red dots indicate the RPs, and the blue triangles represent the test points.
(a)
Map 1(b)
Map 2Figure3 shows the number of received FPs at each RP of the two sites in the offline phase. It can be seen from Figure 2 that most of the RPs receive 10~30 FPs on Map 1, while on Map 2, this number is 10~15. No RP can receive all the FPs’ signals on both maps. In these two typical indoor environments, especially that of Map 1, the dense distribution of rooms and walls greatly attenuated or obstructed the radio signals. Although a large number of FPs could be detected in the whole area, the FPs available at each RP were greatly reduced.Figure 3
The number of received FPs at each RP in the offline phase.
(a)
Map 1(b)
Map 2


## 3.2. Experimental Results with Different FP Selections

Localization performance is always related to the number of FPs (K′) used for positioning. In our method, FPs with certain coverage areas were first selected using a spatial filter. The threshold C0 varied with the layout of the map, as well as the distributed RPs and FPs. Figure 4 shows the CDF of the positioning error when different thresholds were applied to filter the FPs. On Map 1, although there are up to 198 FPs, only 52 FPs can be received by more than 15 RPs out of a total of 92. The number K′ increases to 89 when the threshold drops to 9. With the FP selection criterion, 80 FPs were used for positioning, and almost no performance difference was observed as shown in Figure 4(a). Figure 4(b) shows the experimental results for Map 2. When C0≤50, all 18 FPs were received by all RPs. This may be attributed to the regular distribution of WSN anchor nodes and an indoor environment without many barriers. Considering this, no spatial filter was used in the following experiments for Map 2.Figure 4
Positioning error of the proposed method with different thresholds for the spatial filter.
(a)
Map 1(b)
Map 2Figure5 depicts the average location errors under the three different FP selection criteria, namely, the Strongest FPs, Least-variance FPs, and Combined criterion, when a proper subset of PCs is used for KRR-based localization.Figure 5
Average positioning error of the proposed method using three different FP selection criteria.
(a)
Map 1(b)
Map 2It can be seen that no matter which FP selection criterion is adopted, the positioning method with more selected FPs obtains higher positioning accuracy at both sites. As shown in Figure5(a), when 80 FPs are adopted, that is, 80 FPs are selected from the spatial filtered 98 FPs according to the FP selection criterion, the AE under the three criteria is all approximately 2.4 m. Almost the same result can be observed in Map 2. It is worth mentioning that even if we choose 80 FPs in Map 1 according to a certain FP selection criterion, such as the strongest FPs, only 18 of the RPs throughout the whole site can receive signals from more than 20 FPs. The situation is even worse when only 40 FPs are used; only one of the 92 RPs can receive signals from more than 20 FPs, which may explain the relationship between the positioning error and the number of FPs in Figure 5(a). For the positioning method that needs to establish the model between positions and signals in the offline phase, it is necessary to collect enough FP RSS measurements to extract features, especially in indoor environments with rooms and walls that greatly attenuate or completely block signal propagation. This is a completely different approach from those that apply FP selection in the online phase, as proposed in [23, 24, 31].


## 3.3. Analysis of the Number of PCs

The proper choice ofk (the number of PCs) is important to the proposed method. A larger value of k seems to retain more information after PCA transformation, but it may also retain redundant information and noise, which can be seen from the experimental results. Figure 6 compares the average positioning errors for the different number of FPs versus the number of PCs at both sites.Figure 6
Average positioning error versus the number of PCs. (a) FPs were selected with Least-variance FP selection. (b) FPs were selected with Combined criterion FP selection.
(a)
Map 1(b)
Map 2As illustrated in Figure6, no matter how many FPs are used in the experiments, the positioning accuracy varies with the number of PCs, and there is always an optimal value kp for the highest positioning accuracy, but this value is certainly not the maximal one. For example, in the case of the 80 FPs adopted in Map 1, the optimal kp is 23. Further research found that this value kp ensured that the selected PCs retained more than 85% of the total information if eigenvalues were used to quantify the information contribution of each PC, which provides a relatively simple way to determine this optimal value. It is also worth mentioning that no matter how many PCs are extracted, the premise is to have sufficient FPs’ information for high positioning accuracy. From Figure 6(a), the best performance in the case of 40 FPs is 2.77 m, while with 80 FPs, the worst performance is 2.56 m. The same is true for Map 2, where the corresponding values for the best and worst performance are 2.28 m and 1.88 m. It can also be inferred from Figure 5 that the number of FPs has a greater impact on performance compared with the number of PCs, which means considering enough FPs before choosing the right number of PCs.


## 3.4. Reduction in the Number of RPs

Figure7 shows the comparison of positioning error among Euclidian-WKNN, LASSO, KRR, and the proposed PCA-KRR method with respect to the number of RPs. For Map 1, a total number of 80 FPs selected with the Combined criterion were used with the optimal value of kp. The number of FPs is 18 for Map 2. From the results displayed in Figure 7, we can see that the accuracy of the positioning gradually drops for both maps as the number of RPs used for positioning decreases, while the proposed PCA-KRR method achieves significantly better accuracy for both maps than the other methods regardless of the number of RPs. For instance, when the number of PRs used was 92 and a quarter of RPs (23) in Map 1, the AE of the proposed PCA-KRR method is, respectively, 2.34 m and 2.88 m, which are lower than the 2.71 m and 3.93 m of WKNN and the 2.78 m and 3.20 m of KRR. It is worth noting that the LASSO algorithm shows a different positioning performance for each of the two maps. For Map 1, the reduction in the number of RPs does not affect the positioning accuracy. While for Map 2, the result is different. The reduction of RPs introduces a large positioning error, indicating that LASSO is not robust enough for various environments. Furthermore, from these experimental results, it can be deduced that the density of RP does affect the positioning accuracy of our method, but the influence is limited. For example, when only a quarter of the RPs are used for positioning, i.e., the space between the RPs is as large as 7.2 m, the AE of the PCA-KRR method still reaches 2.88 m and 2.13 m, respectively, in Map 1 and Map 2 when enough FPs are guaranteed. Therefore, it is suggested that the number of FPs should be as high as possible while the density of RPs can be appropriately reduced to ensure the positioning accuracy while reducing the cost of labor.Figure 7
Average positioning error versus the number of RPs of WKNN, LASSO, KRR, and the proposed method. 1/2, 1/3, and 1/4 RPs mean, respectively, a half, a third, and a quarter of the total RPs in the maps.
(a)
Map 1(b)
Map 2


## 3.5. Tuning the Parameters

There are several ways to tune the regularization parameterλ and the bandwidth of the Gaussian kernel σ involved in the proposed method. One effective approach for setting these parameters is the well-known cross-validation (CV). Another approach is to use the training data to find the best values that make the positioning algorithm the most accurate, which has been applied in this paper. We use some of the collected RSS data at the test points for parameter training. The experimental results show that the best value of λ changes only with the positioning site and does not change with the number of RPs and PCs used in the method. This value is 0.001 and 0.0001 in Map 1 and Map 2, respectively. As mentioned earlier, there are also two kernel widths σx and σy for two coordinates. Table 1 shows the positioning error for the different σ values used for positioning in the two maps. When σx and σy take the same value, there is no significant impact on the positioning accuracy, which can be observed from the results for both maps. Whether the two parameters choose the same value or not affects the positioning accuracy within 2%. In our experimental results above, σx and σy share the same value for the simplification of parameter turning.Table 1
Positioning error with different bandwidth of the Gaussian kernelσ used in the method.
SiteNumber of RPsNNumber of PCskσxσyAE/mMap 1922211112.341072.30462512122.5013112.49Map 290814141.771581.7545915151.7811151.74


## 4. Conclusion

In this paper, we have proposed a new RSS fingerprinting-based positioning method for indoor localization. We showed that PCA-based feature reduction can efficiently extract the RSS features related to positions, and the nonlinear ridge regression can build a proper model between the features of RSS and positions. The proposed positioning method, PCA-KRR, provides high positioning accuracy in various indoor environments. Extensive experimental results have demonstrated that the performance of the proposed method is sensitive to the number of APs or anchor nodes, while the sparsity of the RPs did not reduce the accuracy of the method much. This is conducive to practical application for reducing the labor cost of fingerprints collection. We also have shown that the proposed method can efficiently and effectively locate a target through WLAN signals or WSN signals.



---

*Source: 6631585-2021-01-11.xml*