Next Article in Journal
Application Research of Multi-Mode Relay in Future Heterogeneous Networks
Next Article in Special Issue
Improving Accuracy and Reliability of Bluetooth Low-Energy-Based Localization Systems Using Proximity Sensors
Previous Article in Journal
Preparation and Characterization of Chitosan–Alginate Polyelectrolyte Complexes Loaded with Antibacterial Thyme Oil Nanoemulsions
Previous Article in Special Issue
A Novel Method of Adaptive Kalman Filter for Heading Estimation Based on an Autoregressive Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Indoor Localization Based on Wi-Fi Received Signal Strength Indicators: Feature Extraction, Mobile Fingerprinting, and Trajectory Learning

1
Department of Electrical, Electronic and Control Engineering, Hankyong National University, Anseoung 17579, Korea
2
Department of Military Digital Convergence, Ajou Univeristy, Suwon 16499, Korea
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(18), 3930; https://doi.org/10.3390/app9183930
Submission received: 6 August 2019 / Revised: 11 September 2019 / Accepted: 17 September 2019 / Published: 19 September 2019
(This article belongs to the Special Issue Recent Advances in Indoor Localization Systems and Technologies)

Abstract

:
This paper studies the indoor localization based on Wi-Fi received signal strength indicator (RSSI). In addition to position estimation, this study examines the expansion of applications using Wi-Fi RSSI data sets in three areas: (i) feature extraction, (ii) mobile fingerprinting, and (iii) mapless localization. First, the features of Wi-Fi RSSI observations are extracted with respect to different floor levels and designated landmarks. Second, the mobile fingerprinting method is proposed to allow a trainer to collect training data efficiently, which is faster and more efficient than the conventional static fingerprinting method. Third, in the case of the unknown-map situation, the trajectory learning method is suggested to learn map information using crowdsourced data. All of these parts are interconnected from the feature extraction and mobile fingerprinting to the map learning and the estimation. Based on the experimental results, we observed (i) clearly classified data points by the feature extraction method as regards the floors and landmarks, (ii) efficient mobile fingerprinting compared to conventional static fingerprinting, and (iii) improvement of the positioning accuracy owing to the trajectory learning.

Wi-Fi RSSI based indoor localization [1,2,3,4] is one of the standard approaches for indoor localization. It is able to utilize the RSSI measurements received from a large number of access points (APs) that are already built in construction. However, the Wi-Fi RSSI as a function of distance between a receiver (smartphone) and a transmitter (wireless AP) is nonlinear and varying due to interference of the indoor environments such as the other radio signals, walls, and obstacles. To address this problem, many machine learning based localization methods [5,6,7,8] have been developed, which learn a pattern of the RSSI measurements corresponding locations across the interested positioning area. In addition, due to its unbiased estimation capability, it is likely to be combined with other kinds of localization, such as pedestrian localization using inertial measurement unit (IMU) [9,10], visual localization [11,12], and magnetic sensor-based localization [13,14].
In particular, semisupervised learning algorithms have been recently suggested for efficient indoor localization, which reduce the human effort necessary for collecting training data [15,16,17,18,19,20]. For example, for indoor localization, a large amount of unlabeled data can be easily collected by recording only Wi-Fi RSSI measurements without requiring position labels, which can save resources for collection and calibration. By contrast, labeled training data have to be created manually. Adding a large amount of unlabeled data in the semisupervised learning framework can prevent the decrement in the estimation accuracy that occurs when using only a small amount of labeled data.
Given the advantage of semisupervised learning, this study describes (i) feature extraction, (ii) mobile fingerprinting, and (iii) mapless localization for efficient Wi-Fi RSSI-based localization. Figure 1 describes the interconnections between the research parts and the flow of our localization system.
First, the feature extraction considered in this paper aims to find a pattern from raw Wi-Fi RSSI and to reduce dimensionality. It allows not only recognizing different floor levels and different landmarks (e.g., toilet, room, and elevator), but also boosting calculation time. This study implements the multistory estimation including the classification of the floors and landmarks, by using the single Wi-Fi RSSI measurement set.
Second, most approaches to indoor localization have used the conventional static fingerprinting method in the training phase, where a trainer has to collect labeled training data stationarily at every position (or grid) while measuring and labeling Wi-Fi RSSI measurements. For more rapid and efficient collection, this study proposes the mobile fingerprinting method that allows a trainer to continue walking during the collection. Instead of the trainer’s necessity to label the positions corresponding the measurements, the proposed algorithm automatically pseudolabels (or estimates position of) the unlabeled data with a small amount of labeled data. For accurate pseudolabeling, we design a semisupervised regression by considering both the spatial and temporal relationship of the Wi-Fi RSSI sets.
Third, in indoor localization, it is common to apply a filtering method, such as a particle filter, to estimate the position or to boost accuracy [21,22,23]. However, this implies the assumption that a localization service provider can give the floor plan of the area of interest, which causes a large volume of data being transmitted to users. In this paper, the trajectory learning algorithm is integrated with the floor and landmark classification, mobile fingerprinting, and positioning for the expanded multistory building experiment.
Last, a position estimation algorithm is created by combining a particle filter and Gaussian process (GP) [24], to exploit the learned trajectories as a prior function and to use probabilistic GP likelihood by modeling the relationship between Wi-Fi RSSI and position.
We evaluate the proposed algorithms in a multifloor building through the experience of a few users. The experimental results show that (i) clearly classified data points with respect to different floors and landmarks by the feature extraction, (ii) efficient mobile fingerprinting compared to the conventional fingerprinting method, and (iii) improvement of positioning accuracy up to the average 1.2 m in comparison to standard approach thanks to the trajectory learning.
The paper is organized as follows. Section 1 surveys some studies relevant to indoor localizations. In Section 2, details about the mobile fingerprinting are presented. Section 3 examines characteristics of Wi-Fi RSSI and describes a semisupervised discriminant analysis for extracting features from raw Wi-Fi RSSI observations, with a description of floor classification and landmark detection. Section 4 introduces the mapless localization and Section 5 summarizes the experimental results. Finally, Section 6 presents the conclusion.

1. Related Studies

This section provides a survey of studies relevant to indoor localization. Different semisupervised learning techniques are reviewed for feature extraction in Section 1.1 and for mobile fingerprinting in Section 1.2. Mapless localization, in Section 1.3, describes trajectory learning and position estimation.

1.1. Semisupervised Feature Extraction

A Wi-Fi RSSI dataset consists of RSSI values corresponding to a user’s position, obtained from Wi-Fi access points (APs). The dimension of a raw Wi-Fi RSSI dataset is defined as the number of APs scanned in the entire area of interest in a building. Typically, the dimensionality is so high that it is impractical to perform real-time operations. Additionally, many of the elements in a raw RSSI database are usually empty because APs cannot cover a wide area. Therefore, reconstructing the data of a raw database is of paramount importance. Deep learning approaches [25,26,27,28] have been applied to the RSSI indoor localization. The high accuracy reported is due to its capability for feature extraction by the deep neural network structure, which learns complex pattern of the RSSI observations. The feature extraction in the standard deep learning except the autoencoder method [29] is computed among the hidden neural network layers, and then the hidden layers are connected to the latest localization layer. Thus, they cannot be used for our purpose to obtain the low-dimensional feature vector separately.
This paper combines two feature extraction methods: Fisher discriminant analysis (FDA) and principal component analysis (PCA). FDA is an a supervised dimensionality reduction method, and PCA is an unsupervised dimensionality reduction method. Supervised FDA tends to find embedded spaces overfitted to labeled samples. Therefore, it is effective to add unlabeled data, which can be collected easily and in large volumes. By using both labeled and unlabeled data, the combination of FDA and PCA can produce better accuracy compared to that achieved by solely using supervised FDA or unsupervised PCA. In this study, the role of semisupervised discriminant analysis is two-fold: dimension reduction and detection of floor level and landmarks.

1.2. Semisupervised Learning for Mobile Fingerprinting

In this section, the semisupervised learning methods for regression are presented, which are used to build the relationship between the Wi-Fi RSSI data sets and the locations. We then derive the contribution of the proposed mobile fingerprinting algorithm.
Utilization of the unlabeled data has been studied in the semisupervised learning to improve the efficiency and accuracy of the indoor localization. In those works, semisupervised deep learning methods [18,19,30] have been recently developed. The mobile fingerprinting requires the light computation and should be implemented as fast as the mobile device or the server deals with huge amount of the unlabeled data. Therefore, rather than the heavily computational deep learning methods, the support vector machine (SVM)-based [31] semisupervised learning algorithms are applied in this paper, which solves a convex optimization problem.
Semisupervised least square (SSL) [32] adds manifold regularization of unlabeled data into the standard Least Square SVR framework [33]. Because of its linearized setup, this algorithm is fast and useful for real-time application. Semisupervised colocalization (SSC) [34] builds an optimization framework consisting of a singular value decomposition, a manifold regularization, and a loss function. Because SSC estimates the locations of the APs as well as a target’s location, it requires the known locations of the AP, whereas SVR-based semisupervised algorithms do not. Moreover, the large number of tuning parameters and the heavy computation may be a burden. Both SSL and SSC apply the unlabeled data only for making a manifold regularization, through the graph Laplacian. Further progressed utilization of unlabeled data occurs in transductive support vector machine (TSVM) [35] and Laplacian Embedded Regression Least Squares (LapERLS) [36]. TSVM attempts to find the labels of unlabeled data and obtain the decision function. Because finding a solution requires searching all candidate-labels of the unlabeled data, TSVM is not feasible in most applications.
LapERLS introduces an intermediate variable as a substitute for the original labeled data. During the optimization process via Karush–Kuhn–Tucker (KKT) conditions, pseudolabels and a transformed kernel matrix are generated. Then, the pseudolabels and transformed kernel matrix are used as a substitution for the original labeled data and the original kernel matrix, respectively. This algorithm is useful for large-scale problems, due to the light computation necessary for obtaining the transformed kernel matrix, because the standard kernel matrix and Laplacian matrix are decoupled. However, LapERLS becomes inaccurate when only a small number of labeled training data points are used.
We adopt the idea of pseudolabeling from LapERLS because the pseudolabels can compensate for the lack of labeled data. To improve the accuracy of the pseudolabels, we propose adding a temporal relation to unlabeled training data that are collected as time-series. A study [37] employing a Hodrick–Prescott (H–P) filter that captures a smoothed-curve representation for a time series from training data is helpful for assigning time-series pseudolabels. Consequently, our pseudolabeling is able to consider both the spatial and temporal aspect of the training data sets. Note that this pseudolabeling technique based on this semisupervised regression is used for the mobile fingerprinting, not for estimating the user’s position.

1.3. Mapless Localization

Accommodating a mapless situation might be valuable for secure localization operations by keeping information private. The pedestrian localization based on IMU and camera sensors easily produces the integral error so the user should carefully hold the receiver such as smartphone stationary and should not rotate it. This limitation restricts its practical usage. This study eliminates the need to restrict the users’ behavior and the need to assume on the accurate signal propagation model.
Crowd-sourcing has been a useful tool for indoor localization. Because the crowdsourced data are collected from a huge number of different users conveying various mobile devices, it has potential to help solving challenging problems such as heterogeneous hardware [38] and security issues [39]. In this study, we apply crowdsourced data to learn the hidden trajectory and extract the floor plan, to compensate for the absence of true map information. The trajectory learning algorithm originates in demonstration learning [40]. For the purpose of this study, trajectory learning is combined with semisupervised feature extraction in Section 1.1.
Finally, as the position estimator, the particle filter is employed for two reasons: First, this filter can use the learned trajectories as a prior distribution. Second, in the particle filter framework, the likelihood function can be defined as the function referring to the relationship between the RSSI measurement sets and the positions. In this study, the likelihood is defined as the probabilistic model by the Gaussian process.

2. Mobile Fingerprinting Based on Semisupervised Learning

Wi-Fi fingerprint localization estimates a location by matching the currently received Wi-Fi RSSI measurements to those in a training database. For creating this database, the conventional fingerprinting method requires a trainer to manually labels all the Wi-Fi RSSI measurements at every point of the grid. Instead of the time-consuming conventional method, this section introduces a new mobile fingerprinting data collection algorithm that allows the trainer to continue walking without the stationary calibration. In the training phase for data collection, it is common sense that data collecters recognize which floor they are located on. Thus, the proposed mobile fingerprint method aims to obtain 2D position of the unlabeled data.
Suppose y = [ y 1 , , y d ] T R d is the set of the Wi-Fi RSSI measurements received from d Wi-Fi APs. In 2-D space, the user’s location is defined as ( x X , x Y ). The l number of the labeled training data points are defined as the set { ( y i , [ x X i , x Y i ] T ) } i = 1 l with y i X R d . The u number of unlabeled data set { y i } i = l + 1 l + u comprises only the RSSI measurements.
It is desired to find the separate mappings f X : X R and f Y : X R , which denote the relationships between Wi-Fi signal strength and location of the smartphone, using the labeled training data { ( y i , x X i ) } i = 1 l and { ( y i , x Y i ) } i = 1 l , and the unlabeled data { y i } i = l + 1 l + u . Because the models f X and f Y are learned independently, we omit the subscripts of f X , f Y , and of x X , x Y , for simplification.
In the SVM-based semisupervised learning framework, the optimization formulation is as follows,
f * = arg min f H k c i V ( y i , x i , f ) + γ A f A 2 + γ I f I 2 ,
where V is a loss function; f A 2 is the norm of the function in the Reproducing Kernel Hilbert Space H k ; f I 2 is the norm of the function in the low dimensional manifold; and c, γ A , and γ I are the regularization weight parameters. To represent the manifold f I 2 , the method uses graph Laplacian, the so-called graph-based semisupervised learning, and it is also called a semisupervised support vector machine when we use an ϵ -insensitive loss function. Then, the solution is achieved by iterative quadratic programming. More details about the semisupervised algorithm can be found in [37].

2.1. Hodric–Prescott Filter

Let us describe a scenario of mobile fingerprinting where a trainer collects the training data during walking. The observations are naturally recorded in time-series. This section introduces capturing the temporal property from the mobile fingerprint data by exploiting the H–P filter [41]. By using the H–P filter, the optimization problem can be formulated as follows,
min f i = 1 K f ( y i ) x i 2 + γ T i = 3 K f ( y i ) + f ( y i 2 ) 2 f ( y i 1 ) 2
where { ( y i , x i ) } i = 1 K is the time-series labeled training data over discrete time horizon K. The second term renders the sequential functional values f ( y i ) , f ( y i 1 ) , f ( y i 2 ) on a line in the embedded space. The solution of Equation (2) in matrix form is
f = I + γ T D D T 1 x ,
where x = [ x 1 , , x K ] T and
D = 0 0 0 0 0 0 0 1 2 1 0 0 0 1 2 1 0 0 0 0 1 2 1 K × K .
In the following section, this H–P filter-based optimization is combined with the semisupervised learning framework to assign the temporal aspect to the unlabeled data.

2.2. Semisupervised Pseudolabeling

Here, the semisupervised optimization to create accurate pseudolabels for the unlabeled data by considering both the temporal and spatial aspect is presented.
Given the labeled and unlabeled training data { ( y i , x i ) } i = 1 K arranged in chronological order, the optimization for generating pseudolabels X ˜ R K is
min X ˜ R n c 2 E T E + 1 2 μ 1 X ˜ T L X ˜ + 1 2 μ 2 X ˜ T D D T X ˜ subject to : Λ X ˜ X = E ,
where Λ is a diagonal matrix of trade-off parameters with Λ i i = 1 if y i is a labeled data point and Λ i i = 0 if y i is an unlabeled data point. μ 1 and μ 2 represent a trade-off relationship between spatial and temporal correlation, matrix D is defined in Equation (3), and X R K is given by
X = x i if y i labeled 0 else if y i unlabeled ,
and graph Laplacian L is defined as L = B 1 / 2 ( B C ) B 1 / 2 , with the adjacency matrix C and the diagonal matrix B given by B i i = j = 1 K C i j . In general, the edge weights C i j are defined as a Gaussian function of Euclidean distance, given by
C i j = exp x i x j 2 / σ w 2 ,
where σ w 2 is the kernel width parameter.
With the introduction of a multiplier α P , the Lagrangian of Equation (4) is given by
L = c 2 E T E + 1 2 μ 1 X ˜ T L X ˜ + 1 2 μ 2 X ˜ T D D T X ˜ + α P Λ X ˜ X E
The derivatives of Equation (5) with respect to the variables X ˜ , α P , and e set to zero are
L X ˜ = μ 1 L X ˜ + μ 2 D D T X ˜ + Λ α P = 0
L e = c E α P = 0
L α P = Λ X ˜ X E = 0
Substituting Equation (7) into Equation (8) gives
Λ X ˜ 1 c α P = X .
Then, the linear algebraic equations satisfying the KKT conditions are defined as follows,
Φ A = F ,
where
Φ = μ 1 L + μ 2 D D T Λ Λ 1 c
and,
A = X ˜ α P , F = 0 X
Therefore, the optimal pseudolabels are obtained by inverse calculation of the matrix Φ . The advantages of this optimization compared to traditional semisupervised learning are as follows.
  • A closed-form solution can be obtained from the linear algebraic Equation (10), which is faster than for other semisupervised algorithms and requires an iterative quadratic programming (QP) solver.
  • The solution for X ˜ also inherits the sparsity characteristics of a support vector machine [42]. This is beneficial when we manage training data storage by saving reliable training data only, that is, large values of α P i.e., the support vectors.
  • In addition to spatial representation, a time-series representation is considered, where the relevant term, μ 2 D D T , is inserted independently into the optimization problem. Without a substantial increase in computational time, it improves the accuracy of the pseudolabels. Its performance is analyzed in the next section.
In a summary, we address the role of the mobile fingerprinting and the connection to the next sections. The mobile fingerprinting is a method for database construction. Especially, it is useful when amount of the labeled data is not enough by its capability of accurately labelling the positions of the unlabeled data. As shown in Figure 1, the calibrated data from the mobile fingerprinting are used for the following feature extraction in Section 3 and the localization in Section 4.

3. Feature Extraction and Application to Floor Classification and Landmark Detection

Section 3.1 overviews the characteristics of Wi-Fi RSSIs and the need for feature extraction. Section 3.2 describes a semisupervised discriminant algorithm for the floor classification and landmark detection.

3.1. Wi-Fi RSSI Characteristics

This section examines characteristics of Wi-Fi RSSIs and discovers the necessity of the feature extraction for Wi-Fi RSSI-based indoor localization.

3.1.1. Nonlinearity and Uncertainty

According to a propagation model of Wi-Fi RSSIs [43,44,45], the RSSI value decreases exponentially when the distance between a transmitter and a receiver increases linearly. In reality, however, the interruption of the multipath generates uncertain Wi-Fi RSSIs because of the existence of a number of walls and the interference from other radio signals.
For validation, we record Wi-Fi RSSI values on an office floor shown in Figure 2a. Figure 2b shows the RSSI values according to the distance from Wi-Fi AP2; one graph shows the RSSI curve when a user moves along the corridor and the other is when the path is interrupted by walls. Because this study does not assume the that a map is available, it is impossible to predict the change in RSSI patterns with respect to the distance.

3.1.2. Sparsity

The other important characteristic is the sparsity of the raw Wi-Fi RSSI measurements. In a typical building, many Wi-Fi APs are installed, such as commercial APs, private internet-connected devices, and WLAN printers. For example, 193 APs were scanned in our experimental floors, and 531 APs are used in the the UJIIndoorLoc open dataset [46].
Not all the APs can be scanned at one point position because a single AP does not cover the entire positioning area. Thus, the RSSIs of APs located far from a user are recorded empty, as illustrated in Figure 3. The typical method to deal with the empty elements is to replace them by the possibly minimum RSSI value, e.g., −100 dBm. However, the sparsity still exists because the minimum RSSI values replace most of the elements, which may deteriorate localization performance. Therefore, we need to extract a feature. Further, substantially reduced dimensionality is helpful to reduce computational load.

3.2. Semisupervised Discriminant Analysis

Both principal component analysis (PCA) and Fisher discriminant analysis (FDA) have been widely used for feature extraction. These methods play the role of reducing dimensionality by highlighting meaningful elements from the original data vector. FDA is a supervised learning as it uses only labeled data, and PCA is an unsupervised learning because it uses only unlabeled data. The combination of FDA and PCA is categorized into the semisupervised learning.

3.2.1. Generalized Eigenvalue Problem

Remember that y = [ y 1 , , y d ] T R d is the set of the raw Wi-Fi RSSIs. The feature set is given z R r ( 1 r d ) , which is the resultant low-dimensional set by a feature extraction implementation. Let T R d × r be a transformation matrix; then, it has the following relationship,
z = T T y ,
where · T represents the transpose operation. Both FDA and PCA solve the following optimization problem,
T opt = argmax T tr T T B T ( T T C T ) 1 ,
where B is a quantity we want to increase, C is a quantity to be hopefully decreased, and tr ( · ) represents the trace of a matrix.
Further, suppose that { λ i } i = 1 d is the set of generalized eigenvalues and { φ i } i = 1 d is the set of the corresponding generalized eigenvectors. Then, the generalized eigenvalue problem can be defined as
B φ = λ C φ .
The generalized eigenvectors are orthogonal to matrix C:
φ i T C φ j = 0 ,
for i j , and the eigenvectors are normalized:
φ i T C φ i = 1 for i = 1 , 2 , , d .
Additionally, assume that the eigenvalues are organized in descending order:
λ 1 λ 2 λ d .
Finally, the transformation matrix can be summarized as
T opt = λ 1 φ 1 | | λ r φ r | | λ d φ d .
From Equation (19), it can be seen that the influence of the transformation becomes deemphasized as the dimension order increases, because the eigenvalues and eigenvectors are decreasing.

3.2.2. PCA

Let { y } i = 1 u be a set of u unlabeled RSSI observations and S p be the total scatter matrix, given by
S p = i = 1 u ( y i μ ) ( y i μ ) T ,
where μ is the mean of all the observations:
μ = 1 u i = 1 u y i .
Equation (20) can be also expressed in a pairwise form as
S p = i = 1 u y i y i T n μ μ T = 1 2 i , j = 1 u y i y i T 1 2 i , j = 1 u y i y j T = 1 2 i , j = 1 u W i j p ( y i y j ) ( y i y j ) T ,
with W i j p = 1 / u . PCA finds the transformation matrix T p c a where
T pca = argmax T tr T T S p T ( T T T ) 1 .
Note that PCA is an unsupervised method that does not require label information such as floor level. Thus, PCA alone cannot be used for a supervised problem such as floor classification and landmark detection in this paper.

3.2.3. FDA

Let { ( y i , w i ) } i = 1 l be a set of l labeled training data points, where w i { 1 , 2 , , C } is the label of the RSSI vector y i . For example, the label can be either floor level or a user-defined landmark. In FDA, the between-class covariance matrix S b and the within-class covariance matrix S w are defined as
S b = c = 1 C l c ( μ c μ ) ( μ c μ ) T ,
S w = c = 1 C i : w i = c ( Y i μ c ) ( Y i μ c ) T ,
where l c is the number of the labeled data points in class c { 1 , 2 , , C } with l = c = 1 C l c , i : w i = c indicates summation over i such that w i = c , and μ c is the mean vector of the data points in class c, given by
μ c = 1 l c i : w i = c Y i .
The transformation matrix for FDA is given by
T fda = argmax T tr T T S b T ( T T S w T ) 1 .
Equations (24) and (25) can be expressed by the following weight matrices.
W i j b = 1 l 1 l c if w i = w j 1 l if w i w j
W i j w = 1 l c if w i = w j 0 if w i w j

3.2.4. Semisupervised Combination of FDA and PCA

To extend to a semisupervised version, it modifies the original FDA in a way that utilizes the unlabeled data. Let { ( y i , w i ) } i = 1 l + u be the set of both the labeled and unlabeled data points, where l is the number of the labeled data points and u is the number of the unlabeled data points. If { ( y i , w i ) } is the labeled data, w i { 1 , 2 , , C } denotes the class labels; for the unlabeled data, w j = 0 . Analogously to traditional semisupervised discriminant analysis, we define the new weight matrices, modified from Equations (27) and (28):
W i j b S = 1 l 1 l c if w i = w j = c 1 l if w i w j 1 u if w i = 0 or w j = 0
W i j w S = 1 l c if w i = w j = c 0 if w i w j 0 if w i = 0 or w j = 0
The corresponding matrices S b S and S w S are
S b S = 1 2 i , j = 1 l + u W i j b S ( y i y j ) ( y i y j ) T ,
S w S = 1 2 i , j = 1 l + u W i j w S ( y i y j ) ( y i y j ) T .
The generalized eigenvalue problem for the semisupervised version is
T semi = argmax T tr T T S b s T ( T T S w s T ) 1 ,
S b s = α b a l S b S + ( 1 α b a l ) S p ,
S w s = α b a l S w S + ( 1 α b a l ) I ,
where I is the identity matrix and α b a l is a parameter to adjust the balance between PCA and FDA.

3.3. Application to Floor Classification and Landmark Detection

3.3.1. Floor Classification

To apply the proposed semisupervised discriminant analysis algorithm in Section 3.2.4 to the floor classification, the label of the training data should be defined as the floor level. As a result, the transformation matrix plays a role as dividing the classes with respect to each floor level. In the test phase, when a test RSSI measurements set is arrived, it is first processed by the transformation matrix in Equation (13). Second, the floor level is estimated by k-NN method, which selects the averagely nearest class point in the training data to the test data.

3.3.2. Landmark Detection

Landmark detection intends to recognize whether a user is located at preliminarily defined points or not. Landmarks in an indoor environment can be toilets, elevators, and rooms, which create a particular pattern of Wi-Fi RSSI measurements such that their distinct features can be distinctly extracted. In this paper, the landmark detection is used for the trajectory learning, which will be introduced in Section 4.2.
The landmark detection implementation is as follows. First, similar to the floor classification, the label of the training data should be defined as the landmark index for the semisupervised discriminant analysis algorithm in the training phase, which is described in Section 3.2.4. In the test phase, when a test RSSI measurements set has come, the distance on the signal space is calculated between each landmark’s feature set and the current RSSI feature vector. Let D c be the distance given by
D c = μ ¯ c z * , c { landmark 1 , landmark 2 , } ,
where μ ¯ c = T semi T μ c is the center of the training data points belong to the same landmark and z * = T semi T y * is the test RSSI feature vector.
In sum, the semisupervised feature extraction is proposed to deal with the nonlinearity and sparsity of the raw Wi-Fi RSSI measurements. Two independent feature extraction models according to the different label type are applied to the floor classification and the landmark detection. In particular, the landmark detection can trim a trajectory sample by detecting two landmarks as the start and end points of a path, which is used for the trajectory learning in the next section.

4. Mapless Localization

In this section, we achieve localization that does not require true map information. Section 4.1 formalizes the position estimation based on a particle filter and Gaussian process. Section 4.2 introduces learning trajectories collected from a crowd for creating map information.

4.1. 2-D Position Estimator Based on Particle Filter and Gaussian Process

A particle filter involves obtaining a recursive estimate of the posterior distribution x R 2 at current time k, given all the observations z 1 : k R r . When we define { x k i , w k i } i = 1 N p as the set of N p particles and corresponding weights, the posterior density function is
p x k | z 1 : k i = 1 N p w k i δ ( x k x k i ) .
In Equation (37), δ ( · ) is the Dirac delta function. The weights are normalized so that i = 1 N p w k i = 1 . The estimate of the state x k is given by
x ^ k = E [ x k | z 1 : k ] i = 1 N p w k i x k i ,
and the weights are updated using the likelihood p ( z k | x k i ) :
w k i = w k 1 i p z k | x k i .
In this study, the likelihood is defined as a Gaussian process to achieve a nonlinear relationship between positions and RSSI observations, and is given by
p z k | x k = j = 1 r p ( z k j | x k ) = j = 1 r N ( μ x k j , σ x k j ) ,
where N ( μ x k j ( · ) , σ x k j ( · ) ) is a Gaussian distribution whose mean and variance are as follows,
μ x k j = k * T ( K ( X ˜ , X ˜ ) + σ G P 2 I ) 1 z * j σ x k j 2 = k ( x k , x k ) k * T ( K ( X ˜ , X ˜ ) + σ G P 2 I ) 1 k * .
In Equation (40), training input X ˜ R n × 2 is defined as the pseudolabels of the x-y positions obtained in Section 2.2, and training input z * j R n is the Wi-Fi observation set corresponding to the j-th pseudolabels in X ˜ . Further, the kernel function and matrix, k ( · , · ) and K, respectively, are defined by Gaussian kernel, and k * is the n × 1 vector of covariances between x k and X ˜ . More details about the derivation of the Gaussian process and parameter selection can be found in the work by the authors of [47].
Here, z j , for j = 1 , , r ( r = 10 in this paper) are the PCA-driven observations from Section 3.2.2, that is, z = T pca T y , where y is the raw Wi-Fi RSSIs and T pca is the transformation matrix in Equation (23). As r = 10 , ten different Gaussian process models as in Equation (39) are used.

4.2. Trajectory Learning from a Crowd

In the particle filter framework, the sampling relies on the prior probability p ( x t | x t 1 ) . Under the unknown-map situation, the learned trajectory compensates for the absence of the true map. The prior function is defined as follows,
p ( x t | x t 1 ) = P H F · P T L ,
where P T L is a learned map and P H F is for capturing a smoothed-curve representation of a time-series trajectory using the H–P filter introduced in Section 2.1:
P H F = N x t x t 2 2 x t 1 2 ; 0 , σ v ,
In Equation (42), it does not require the estimation of velocity.
Now, we describe how to build P T L . In indoor spaces, people trace similar trajectories to save their travel distance and time. The underlying idea for trajectory learning is that people tend to follow similar trajectories when they have the same departure point and destination. The departure and destination points are automatically obtained by the landmark detection algorithm described in Section 3.3.2. Suppose that we sample the trajectories X j k obtained from M different people k = 0 , , M 1 and that the trajectories may have different trajectory lengths j = 0 , , N ( k ) 1 . Then, there might exist a hidden intended trajectory h 1 : k , representative of all X j k . For example, h 1 : k can be the average trajectory of X j k .
The goal is to learn the hidden intended trajectory h t . Dynamic time warping with Kalman smoothing [40] is applied to this problem. The hidden trajectory h t has length O at t = 0 , , O 1 . The length O is initialized to an ample size, namely, O = 2 / M k = 1 M N ( k ) . The trajectory learning method treats the samples X j k as observations of h t such that
h t + 1 = f ( h t ) + w t h , w h N ( 0 , Σ h ) X j k = h τ j k + w j X , w X N ( 0 , Σ X ) .
Here, the covariance Σ h and Σ X of the Gaussian noises w h and w X should be estimated. The subscript τ j k is the time index of h corresponding to X j k .
Log-likelihood maximization is used to estimate the hidden trajectory h t and the time indices τ j k , as follows,
max τ , Σ ( · ) log p ( h , τ ; Σ z , Σ X ) ,
where Σ ( · ) refers to Σ h and Σ X . In the work by the authors of [40], an iterative algorithm for solving Equation (43) is introduced. First, while keeping τ constant, it updates the covariance matrix Σ ( · ) by separate E and M steps. In the E step, a Kalman smoother is applied to obtain the pairwise marginals over the latent variables h 1 : t ; the M step updates the covariance matrix Σ ( · ) . Then, dynamic time warping is applied to calculate τ ^ with the given Σ ( · ) through the following optimization.
τ ^ = arg max τ log p ( h , τ ; Σ ( · ) ) = arg max τ k = 0 M 1 j = 0 N k 1 [ log p ( X j k | h τ j k | T , τ j k ) + log p ( τ j k | τ j 1 k ) ] .
The details of the dynamic programming for solving Equation (44) are described in the work by the authors of [40].
Now, we describe how to collect observations, that is, M trajectories X j k . By the landmark detection introduced in Section 3.3.2, we can estimate the edge of any pair of trajectories. Additionally, the elements of the trajectories are filled with the estimates obtained from the particle filter from Section 4.1.
To generalize map learning, suppose that we sample a set of n learned trajectories h 1 : n = { h ( 1 ) , h ( 2 ) , , h ( n ) } , where each h ( · ) might have different start and end points and each h 1 : n is exploited to obtain P T L in Equation (41). We assume that P T L follows the Gaussian distribution given by
P T L N x k ; μ T L , Σ T L
and
μ T L = arg h min h 1 : n x k ,
where μ T L is the trajectory among h 1 : n closest to x k . The variance Σ T L is set to the estimated covariance Σ h of the learn trajectory, which is obtained from Kalman smoothing in Equation (43). The covariance Σ X , which is also estimated in Equation (43), indicates how far the samples’ trajectory is apart from the learned trajectory. Therefore, by inspecting Σ X , we might detect an outlier that might arise when someone does not follow the common trajectory. In this paper, the outliers are filtered out by 95% confidence interval of the trajectory samples.

5. Experiment

The experimental field is a three-story building, where the area of each floor is 47 m × 36 m. The total number of scanned Wi-Fi APs is 193 and ten people are employed for training and testing.
Training data are divided into labeled data, composed of RSSI values and the corresponding labels and unlabeled data, which only consists of RSSI measurements. The labels of the labeled data are of three different categories, namely, floor level, landmark, and 2D position. The algorithm can be seen as a hierarchical structure; the floor level is first determined, and then the position is estimated. For the trajectory learning, the trajectory samples are generated whenever two landmarks indicating the start and end points (of the trajectory) are detected. Then, they are sent to a server to learn the hidden trajectories. Subsequently, the newly learned map information is used to update the particle filter algorithm.

5.1. Mobile Fingerprint

In this experiment, we compare all the benchmarks, that is, SSL, SSC, LapERLS, and the additional supervised k-NN algorithm. The parameters of the proposed algorithm are defined as μ 1 = 5 , μ 2 = 2 , c = 1 in Equation (5), and the other parameters of the compared methods are selected for by best performance. For the experimental study, we vary the number of the used labeled training data points out of the fixed total of the training data. When 283 number of additional test data points are used, and the results using 283 and 93 training data points are shown in Figure 4a,b, respectively. From the results shown in both parts of this figure, we observe that our algorithm outperforms the compared methods. In the cases of 100% and 75% labeled data in Figure 4b, SSC provides a slightly smaller error than the other methods. However, considering the advantage given to SSC, that is, as it knows the locations of the Wi-Fi APs, this result is not noteworthy. The major contribution of our algorithm can be seen in the case where a very small number of labeled data points are used. From Figure 4a,b, our algorithm shows a slightly increasing error as fewer labeled data points are used, whereas the others exhibit a substantially increasing error when comparing the case with 25% labeled data to that with 10% case.
Finally, Figure 4c shows the CPU running time of the compared algorithms with respect to the percentage of labeled data. SSC requires more computational time than the other methods. The proposed algorithm needs slightly more time than LapERLS and SSL (an additional 0.2 s at most), due to the additional time-series term to be solved in the optimization process.

5.2. Floor Classification and Landmark Detection

For the feature extraction, the dimensionality of the Wi-Fi RSSI set is reduced from 193 to 10, i.e., d = 193 and r = 10 . All the RSSI values are scaled into the range [0, 1]. Figure 5a compares FDA (supervised) to the semisupervised discriminant analysis for floor level estimation. Note that because PCA is an unsupervised method, the PCA alone cannot be used for the supervised floor and landmark detections. To estimate the floor level, the k-NN method with k = 1 is employed. In Figure 5a, the ratio of the used labeled data varies from 10% to 100%. The developed algorithm provides better accuracy than FDA. The most noticeable result appears when the ratio of the labeled data is small. Although 10% of the data is labeled, semisupervised analysis achieves 95% accuracy whereas supervised analysis results in 60% accuracy.
Selecting a tuning the balancing parameter in Equation (33) depends on a cross validation. The optimal value is selected by minimizing the training error. In Figure 5b, the floor estimation accuracy is shown with respect to the variation of parameter values α b a l introduced in Equation (33) and the ratio of the labeled training data. Setting a relatively large α b a l value intends to focus more on FDA than on PCA in the semisupervised learning analysis. In this study, α b a l = 0 . 2 is used in the rest of all experiments.
Landmarks in an indoor environment, such as toilets, elevators, and rooms, create a particular pattern of Wi-Fi RSSI measurement sets so that their distinct features can be distinctly extracted. Figure 6a clarifies the extracted features at some landmarks using the semisupervised discriminant analysis algorithm. Figure 6b shows the distance on the signal space according to a user’s path. The proposed algorithm detects the moment when the user arrives at each pre-defined landmark site with the threshold 0.3 in Figure 6b. Consequently, we can apply the landmark detection algorithm to calibrate any user’s trajectory, by determining the start and end positions. Note that these trimmed trajectories are used for the trajectory learning in Section 4.

5.3. Trajectory Learning

Given the trajectory samples obtained from the landmark detection, Figure 7 shows the learned map on a floor with respect to the number of the participants. In Figure 7, seven different trajectories having different start and end landmarks are drawn with the different colors. We can observe that some learned trajectories are still inaccurate in Figure 7a, as they do not match the area on the true map when only two participants join. In Figure 7b, the more participants join in, the closer the learned map approaches the ground truth.

5.4. Localization

To confirm the positioning error, two localization cases are compared with 50 number of particles setup for executing the particle filter. Figure 8a,b contains the localization results when using the learned trajectories as the prior for the particle filter. Figure 8c,d contains the localization results without using the learned trajectories. Because the user is moving, it is unable to measure the ground truths of the user’s path all the time. Instead, some waypoints are designated to alarm the user to record its location with time stamps. Based on these waypoints, the average error is calculated. Ten positioning experiments by ten different users for each case are implemented. The average error of the developed algorithm is 2.2 m. By contrast, the average error for the localization without the learned trajectories is 3.4 m. In the worst case, the proposed algorithm has 2.8 m error and the benchmark one has 5.1 m. In the best case, each has 1.3 m and 2.5 m error, respectively. Thus, the map learning improves the localization accuracy by 1.2 m on average. It is noted that in Figure 8c the positioning is remarkably inaccurate because of the indoor environment composing majority of open space, whereas the localization on the aisles are relatively more precise. On the other hand, as shown in Figure 8a, the accurately learned map information supports to overcome this environmental restriction by enhancing the accuracy. The last remark of the experimental result analysis between Figure 8b,d is the positioning regarding the south room. Because the true map information is not given in this paper, it is unable for the localization in Figure 8d to recognize the wall between the room and aisle. On the other hand, due to aid of the trajectory learning, the proposed method recognizes the isolation by the wall, in which Figure 8b shows the clear–separate position estimations between the room and the aisle.

6. Conclusions

In this study, we investigated indoor localization performing simultaneous floor classification, landmark detection, positioning, and map learning. The study was divided into three topics: (i) feature extraction from Wi-Fi RSSIs for floor classification and landmark detection, (ii) mobile fingerprinting and pseudolabeling for positioning, and (iii) mapless localization.
In the first part, characteristics of Wi-Fi RSSIs were determined by pattern recognition, using a semisupervised discriminant analysis. The proposed algorithm extracted the features from the noisy Wi-Fi signal data by reducing them from a high to a low dimension. During this process, worthless elements were removed to obtain clustered data points according to the different floors and different landmarks. At the same time, by investigating the distance between a test point and the training data on the reduced signal space, we successfully detected the floor and landmark changes.
The second part addressed the efficiency of fingerprinting. Compared to conventional static fingerprinting, our algorithm improved the efficiency for collecting the training data because a user could be mobile during the collection without manually labelling the position at every grid point in the area of interest. The proposed pseudolabeling algorithm based on the semisupervised regression aimed to obtain accurate pseudolabeled positions for the unlabeled training data. Considering both the spatial and temporal aspect, we formalized the optimization based on a graph Laplacian and H–P filter. Further, the optimization provided a closed-form solution, which enabled fast computation.
In the last part, we considered the situation where the true map information is not available. The key idea was crowdsourcing, from which we can obtain trajectory samples. From the experimental results, as more participants joined in, the learned map was updated more accurately. The experiments conducted in this study involved floors in a multifloor office building. Many people participated in validating our algorithm. The integration of all the parts led to accurate localization.

Author Contributions

J.Y. designed and developed the learning based localization algorithm, and implemented and validated the experiments; J.P. reviewed and edited the paper.

Funding

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean Government (MSIT) (No. 2019R1F1A1057516).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Zhang, L.; Valaee, S.; Xu, Y.; Ma, L.; Vedadi, F. Graph-based semisupervised learning for indoor localization using crowdsourced data. Appl. Sci. 2017, 7, 467. [Google Scholar] [CrossRef]
  2. Zhang, Z.; Tian, Z.; Zhou, M.; Li, Z.; Wu, Z.; Jin, Y. WIPP: Wi-Fi compass for indoor passive positioning with decimeter accuracy. Appl. Sci. 2016, 6, 108. [Google Scholar] [CrossRef]
  3. Hernández, N.; Ocaña, M.; Alonso, J.; Kim, E. Continuous space estimation: Increasing WiFi-based indoor localization resolution without increasing the site-survey effort. Sensors 2017, 17, 147. [Google Scholar] [CrossRef] [PubMed]
  4. Zheng, L.; Hu, B.; Chen, H. A high accuracy time-reversal based WiFi indoor localization approach with a single antenna. Sensors 2018, 18, 3437. [Google Scholar] [CrossRef] [PubMed]
  5. Tran, H.Q.; Ha, C. Improved Visible Light-Based Indoor Positioning System Using Machine Learning Classification and Regression. Appl. Sci. 2019, 9, 1048. [Google Scholar] [CrossRef]
  6. Wang, X.; Gao, L.; Mao, S. CSI phase fingerprinting for indoor localization with a deep learning approach. IEEE Internet Things J. 2016, 3, 1113–1123. [Google Scholar] [CrossRef]
  7. Yan, J.; Zhao, L.; Tang, J.; Chen, Y.; Chen, R.; Chen, L. Hybrid kernel based machine learning using received signal strength measurements for indoor localization. IEEE Trans. Veh. Technol. 2017, 67, 2824–2829. [Google Scholar] [CrossRef]
  8. Chen, Z.; Zou, H.; Yang, J.; Jiang, H.; Xie, L. WiFi Fingerprinting Indoor Localization Using Local Feature-Based Deep LSTM. IEEE Syst. J. 2019. [Google Scholar] [CrossRef]
  9. Ruiz, A.R.J.; Granja, F.S.; Honorato, J.C.P.; Rosas, J.I.G. Accurate pedestrian indoor navigation by tightly coupling foot-mounted IMU and RFID measurements. IEEE Trans. Instrum. Meas. 2011, 61, 178–189. [Google Scholar] [CrossRef]
  10. Zhao, Y.; Liang, J.; Sha, X.; Yu, J.; Duan, H.; Shi, G.; Li, W.J. Estimation of pedestrian altitude inside a multi-story building using an integrated micro-IMU and barometer device. IEEE Access 2019, 7, 84680–84689. [Google Scholar] [CrossRef]
  11. Brachmann, E.; Rother, C. Learning less is more-6d camera localization via 3d surface regression. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Salt Lake City, UT, USA, 18–22 June 2018; pp. 4654–4662. [Google Scholar]
  12. Elloumi, W.; Latoui, A.; Canals, R.; Chetouani, A.; Treuillet, S. Indoor pedestrian localization with a smartphone: A comparison of inertial and vision-based methods. IEEE Sens. J. 2016, 16, 5376–5388. [Google Scholar] [CrossRef]
  13. Wang, G.; Wang, X.; Nie, J.; Lin, L. Magnetic-based Indoor Localization using Smartphone via a Fusion Algorithm. IEEE Sens. J. 2019, 19, 6477–6485. [Google Scholar] [CrossRef]
  14. Shu, Y.; Bo, C.; Shen, G.; Zhao, C.; Li, L.; Zhao, F. Magicol: Indoor localization using pervasive magnetic field and opportunistic WiFi sensing. IEEE J. Sel. Areas Commun. 2015, 33, 1443–1457. [Google Scholar] [CrossRef]
  15. Zhou, M.; Tang, Y.; Nie, W.; Xie, L.; Yang, X. GrassMA: Graph-based semisupervised manifold alignment for indoor WLAN localization. IEEE Sens. J. 2017, 17, 7086–7095. [Google Scholar] [CrossRef]
  16. Zhou, M.; Tang, Y.; Tian, Z.; Xie, L.; Nie, W. Robust neighborhood graphing for semisupervised indoor localization with light-loaded location fingerprinting. IEEE Internet Things J. 2017, 5, 3378–3387. [Google Scholar] [CrossRef]
  17. Yoo, J.; Kim, H. Target localization in wireless sensor networks using online semisupervised support vector regression. Sensors 2015, 15, 12539–12559. [Google Scholar] [CrossRef] [PubMed]
  18. Mohammadi, M.; Al-Fuqaha, A.; Guizani, M.; Oh, J.S. Semisupervised deep reinforcement learning in support of IoT and smart city services. IEEE Internet Things J. 2017, 5, 624–635. [Google Scholar] [CrossRef]
  19. Gu, Y.; Chen, Y.; Liu, J.; Jiang, X. Semi-supervised deep extreme learning machine for Wi-Fi based localization. Neurocomputing 2015, 166, 282–293. [Google Scholar] [CrossRef]
  20. Xia, Y.; Ma, L.; Zhang, Z.; Wang, Y. Semisupervised Positioning Algorithm in Indoor WLAN Environment. In Proceedings of the 2015 IEEE 81st Vehicular Technology Conference (VTC Spring), Glasgow, UK, 11–14 May 2015. [Google Scholar]
  21. Wu, Z.; Jedari, E.; Muscedere, R.; Rashidzadeh, R. Improved particle filter based on WLAN RSSI fingerprinting and smart sensors for indoor localization. Comput. Commun. 2016, 83, 64–71. [Google Scholar] [CrossRef]
  22. Xie, H.; Gu, T.; Tao, X.; Ye, H.; Lu, J. A reliability-augmented particle filter for magnetic fingerprinting based indoor localization on smartphone. IEEE Trans. Mob. Comput. 2015, 15, 1877–1892. [Google Scholar] [CrossRef]
  23. Pak, J.M.; Ahn, C.K.; Shmaliy, Y.S.; Lim, M.T. Improving reliability of particle filter-based localization in wireless sensor networks via hybrid particle/FIR filtering. IEEE Trans. Ind. Inf. 2015, 11, 1089–1098. [Google Scholar] [CrossRef]
  24. Ko, J.; Fox, D. GP-BayesFilters: Bayesian filtering using Gaussian process prediction and observation models. Auton. Robots 2009, 27, 75–90. [Google Scholar] [CrossRef] [Green Version]
  25. Belmonte-Hernández, A.; Hernández-Peñaloza, G.; Gutiérrez, D.M.; Álvarez, F. SWiBluX: Multi-Sensor Deep Learning Fingerprint for precise real-time indoor tracking. IEEE Sens. J. 2019, 19, 3473–3486. [Google Scholar] [CrossRef]
  26. Adege, A.; Lin, H.P.; Tarekegn, G.; Jeng, S.S. Applying deep neural network (DNN) for robust indoor localization in multi-building environment. Appl. Sci. 2018, 8, 1062. [Google Scholar] [CrossRef]
  27. Wang, X.; Gao, L.; Mao, S.; Pandey, S. DeepFi: Deep learning for indoor fingerprinting using channel state information. In Proceedings of the 2015 IEEE Wireless Communications and Networking Conference (WCNC), New Orleans, LA, USA, 9–12 March 2015; pp. 1666–1671. [Google Scholar]
  28. Wang, X.; Gao, L.; Mao, S.; Pandey, S. A Deep Learning based Approach for Indoor Localization. Technology 2017, 66, 763–776. [Google Scholar]
  29. Khatab, Z.E.; Hajihoseini, A.; Ghorashi, S.A. A fingerprint method for indoor localization using autoencoder based deep extreme learning machine. IEEE Sens. Lett. 2017, 2, 1–4. [Google Scholar] [CrossRef]
  30. Jiang, X.; Chen, Y.; Liu, J.; Gu, Y.; Hu, L. FSELM: Fusion semisupervised extreme learning machine for indoor localization with Wi-Fi and Bluetooth fingerprints. Soft Comput. 2018, 22, 3621–3635. [Google Scholar] [CrossRef]
  31. Burges, C.J. A tutorial on support vector machines for pattern recognition. Data Min. Knowl. Discov. 1998, 2, 121–167. [Google Scholar] [CrossRef]
  32. Yoo, J.; Kim, H.J. Online Estimation using Semisupervised Least Square SVR. In Proceedings of the 2014 IEEE International Conference on Systems, Man, and Cybernetics (SMC), San Diego, CA, USA, 5–8 October 2014. [Google Scholar]
  33. Suykens, J.A.; Vandewalle, J. Least squares support vector machine classifiers. Neural Process. Lett. 1999, 9, 293–300. [Google Scholar] [CrossRef]
  34. Pan, J.J.; Pan, S.J.; Yin, J.; Ni, L.M.; Yang, Q. Tracking mobile users in wireless networks via semisupervised colocalization. IEEE Trans. Pattern Anal. Mach. Intell. 2012, 34, 587–600. [Google Scholar] [CrossRef] [PubMed]
  35. Chapelle, O.; Vapnik, V.; Weston, J. Transductive Inference for Estimating Values of Functions; NIPS: Denver, CO, USA, 1999; Volume 12, pp. 421–427. [Google Scholar]
  36. Chen, L.; Tsang, I.W.; Xu, D. Laplacian embedded regression for scalable manifold regularization. IEEE Trans. Neural Networks Learn. Syst. 2012, 23, 902–915. [Google Scholar] [CrossRef] [PubMed]
  37. Yoo, J.; Johansson, K.H. Semi-supervised learning for mobile robot localization using wireless signal strengths. In Proceedings of the 2017 International Conference on Indoor Positioning and Indoor Navigation (IPIN), Sapporo, Japan, 18–21 September 2017. [Google Scholar]
  38. Fang, S.H.; Wang, C.H. A novel fused positioning feature for handling heterogeneous hardware problem. IEEE Trans. Commun. 2015, 63, 2713–2723. [Google Scholar] [CrossRef]
  39. Fang, S.H.; Chuang, C.C.; Wang, C. Attack-resistant wireless localization using an inclusive disjunction model. IEEE Trans. Commun. 2012, 60, 1209–1214. [Google Scholar] [CrossRef]
  40. Abbeel, P.; Coates, A.; Ng, A.Y. Autonomous helicopter aerobatics through apprenticeship learning. Int. J. Robot. Res. 2010. [Google Scholar] [CrossRef]
  41. Ravn, M.O.; Uhlig, H. On adjusting the Hodrick–Prescott filter for the frequency of observations. Rev. Econ. Stat. 2002, 84, 371–376. [Google Scholar] [CrossRef]
  42. Suykens, J.A.; Lukas, L.; Vandewalle, J. Sparse approximation using least squares support vector machines. In Proceedings of the 2000 IEEE International Symposium on Circuits and Systems (ISCAS), Geneva, Switzerland, 28–31 May 2000; Volume 2, pp. 757–760. [Google Scholar]
  43. Zayets, A.; Steinbach, E. Robust WiFi-based indoor localization using multipath component analysis. In Proceedings of the 2017 International Conference on Indoor Positioning and Indoor Navigation (IPIN), Sapporo, Japan, 18–21 September 2017. [Google Scholar]
  44. Wu, D.; Zhang, D.; Xu, C.; Wang, H.; Li, X. Device-free WiFi human sensing: From pattern-based to model-based approaches. IEEE Commun. Mag. 2017, 55, 91–97. [Google Scholar] [CrossRef]
  45. Zhuang, Y.; Li, Y.; Lan, H.; Syed, Z.; El-Sheimy, N. Smartphone-based WiFi access point localisation and propagation parameter estimation using crowdsourcing. Electron. Lett. 2015, 51, 1380–1382. [Google Scholar] [CrossRef]
  46. Torres-Sospedra, J.; Montoliu, R.; Martínez-Usó, A.; Avariento, J.P.; Arnau, T.J.; Benedito-Bordonau, M.; Huerta, J. UJIIndoorLoc: A new multi-building and multifloor database for WLAN fingerprint-based indoor localization problems. In Proceedings of the 2014 International Conference on Indoor Positioning and Indoor Navigation (IPIN), Busan, Korea, 27–30 October 2014; pp. 261–270. [Google Scholar]
  47. Rasmussen, C.E. Gaussian Processes for Machine Learning; MIT Press: Cambridge, MA, USA, 2006. [Google Scholar]
Figure 1. Wi-Fi received signal strength indicator (RSSI)-based indoor localization framework. In the training phase, the mobile fingerprinting for data collection and the feature extraction are conducted. The feature extraction models are applied to the floor level and landmark detections in the test phase. Also, the trajectory learning is proposed for improving position estimation accuracy. The positioning composes the 2D position estimation, floor classification, and landmark detections by using the same RSSI measurements.
Figure 1. Wi-Fi received signal strength indicator (RSSI)-based indoor localization framework. In the training phase, the mobile fingerprinting for data collection and the feature extraction are conducted. The feature extraction models are applied to the floor level and landmark detections in the test phase. Also, the trajectory learning is proposed for improving position estimation accuracy. The positioning composes the 2D position estimation, floor classification, and landmark detections by using the same RSSI measurements.
Applsci 09 03930 g001
Figure 2. (a) Experimental floor plan and (b) RSSI value of Wi-Fi AP2 according to the distance.
Figure 2. (a) Experimental floor plan and (b) RSSI value of Wi-Fi AP2 according to the distance.
Applsci 09 03930 g002
Figure 3. Sparsity of a raw Wi-Fi RSSI set.
Figure 3. Sparsity of a raw Wi-Fi RSSI set.
Applsci 09 03930 g003
Figure 4. Localization results of the compared algorithms when we vary the percentage of the labeled training data: (a) case with 283 training data points, (b) case with 93 training data points, and (c) running time of the compared semisupervised algorithms for pseudolabeling with respect to the percentage of used training data.
Figure 4. Localization results of the compared algorithms when we vary the percentage of the labeled training data: (a) case with 283 training data points, (b) case with 93 training data points, and (c) running time of the compared semisupervised algorithms for pseudolabeling with respect to the percentage of used training data.
Applsci 09 03930 g004
Figure 5. (a) Comparison of the supervised discriminant analysis and the semisupervised discriminant analysis for floor level estimation and (b) floor estimation accuracy with respect to the variation of α b a l parameter values introduced in Equation (33) and the ratio of the labeled training data.
Figure 5. (a) Comparison of the supervised discriminant analysis and the semisupervised discriminant analysis for floor level estimation and (b) floor estimation accuracy with respect to the variation of α b a l parameter values introduced in Equation (33) and the ratio of the labeled training data.
Applsci 09 03930 g005
Figure 6. (a) Result of feature extraction from dataset on the landmarks. (b) Variation of the distance metric D c in Equation (36), corresponding to a user’s movement in the following order; room 1 elevator 1 toilet 1 elevator 2 toilet 2 elevator 1 room 1 .
Figure 6. (a) Result of feature extraction from dataset on the landmarks. (b) Variation of the distance metric D c in Equation (36), corresponding to a user’s movement in the following order; room 1 elevator 1 toilet 1 elevator 2 toilet 2 elevator 1 room 1 .
Applsci 09 03930 g006
Figure 7. Experimental result of trajectory learning.
Figure 7. Experimental result of trajectory learning.
Applsci 09 03930 g007
Figure 8. Positioning results with the trajectory learning (a,b) and without the trajectory learning (c,d).
Figure 8. Positioning results with the trajectory learning (a,b) and without the trajectory learning (c,d).
Applsci 09 03930 g008

Share and Cite

MDPI and ACS Style

Yoo, J.; Park, J. Indoor Localization Based on Wi-Fi Received Signal Strength Indicators: Feature Extraction, Mobile Fingerprinting, and Trajectory Learning. Appl. Sci. 2019, 9, 3930. https://doi.org/10.3390/app9183930

AMA Style

Yoo J, Park J. Indoor Localization Based on Wi-Fi Received Signal Strength Indicators: Feature Extraction, Mobile Fingerprinting, and Trajectory Learning. Applied Sciences. 2019; 9(18):3930. https://doi.org/10.3390/app9183930

Chicago/Turabian Style

Yoo, Jaehyun, and Jongho Park. 2019. "Indoor Localization Based on Wi-Fi Received Signal Strength Indicators: Feature Extraction, Mobile Fingerprinting, and Trajectory Learning" Applied Sciences 9, no. 18: 3930. https://doi.org/10.3390/app9183930

APA Style

Yoo, J., & Park, J. (2019). Indoor Localization Based on Wi-Fi Received Signal Strength Indicators: Feature Extraction, Mobile Fingerprinting, and Trajectory Learning. Applied Sciences, 9(18), 3930. https://doi.org/10.3390/app9183930

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop