Research Article  Open Access
Construction of TimeStamped Mobility Map for Path Tracking via SmithWaterman Measurement Matching
Abstract
Path tracking in wireless and mobile environments is a fundamental technology for ubiquitous locationbased services (LBSs). In particular, it is very challenging to develop highly accurate and costefficient tracking systems applied to the anonymous areas where the floor plans are not available for security and privacy reasons. This paper proposes a novel path tracking approach for large WiFi areas based on the timestamped unlabeled mobility map which is constructed from SmithWaterman received signal strength (RSS) measurement matching. Instead of conventional location fingerprinting, we construct mobility map with the technique of dimension reduction from the raw measurement space into a lowdimensional embedded manifold. The feasibility of our proposed approach is verified by the realworld experiments in the HKUST campus WiFi networks, sMobileNet. The experimental results prove that our approach is adaptive and capable of achieving an adequate precision level in path tracking.
1. Introduction
The recent decade has witnessed a growing interest in the locationbased applications and services for both indoor and outdoor environments [1–4]. Since the WiFi networks are now widely available, the possibility of tracking people’s motion paths by using the WiFi received signal strength (RSS) allows the ubiquitous contextawareness and several potential innovative services [5]. For instance, if the shoppers’ paths are tracked by the retailers in a store, the sales information and the related advertisements could be pushed in based on the shoppers’ realtime locations [6]. As another example, the hospitals can utilize the patients’ path information to identify whether they are in an emergency situation and also assign the closest doctors or nurses to see the patients, if necessary [7].
A variety of wireless network techniques have been considered for location tracking in indoor or outdoor environments. Although the popular and widely used GPS can provide accurate information for outdoor localization and path navigation services, the positioning signals are generally blocked in the indoor or underground scenarios [8, 9]. To solve this problem, the WiFi network is chosen as the favorite technique to achieve indoor localization and tracking due to the popularity in public hotspots and low cost for the deployment in practice [10–13]. In the most recent WiFi localization and tracking approaches, the sitesurvey measurement on RSS fingerprints is required in the offline phase to construct the RSS radio map associated with the target area [14–18]. However, the adaptation degradation problem occurs due to the time consuming and labor intensive work on fingerprint recording [19]. To solve this problem, Wang et al. in [20] introduced a new idea of mapping the people’s motion paths into a mobility map in which the location points (LPs) are connected by transition relations. As discussed in [20], there are three categories of LPs involved in people’s motion paths: (i) personal common locations (PCLs) which many people have spent a lot of time in, (ii) crucial locations (CLs) where multiple adjacent paths intersect, and (iii) ordinary locations which are used to describe the transition relations between neighboring LPs. Each LP is formed by merging the similar measurements which are recorded from assisted GPS (AGPS), WiFi, and cellular networks. However, the constructed mobility map in [20] fails to consider the timestamp relations of measurements. In our previous work [21], we found that, for the calculation of measurement similarities, the timestamps and signal strengths are two sides of a coin. With this idea, we performed spectral clustering on RSS shotgun reads based on the combination of timestamps and signal strengths and also refer to KullbackLeibler divergence of RSS distributions in different LPs to conduct mobility map construction [21]. The most significant problem to limit the practical use for the mobility map in [21] is the low precision and ambiguity in PCL identification, which means that the PCLs cannot be precisely and uniquely identified from the mobility map.
To overcome the disadvantages of the conventional approaches, we propose the tracking solution based on the timestamped mobility map constructed from WiFi RSS measurement matching in this paper. This solution complies with three basic prerequisites: (i) it can be applied to the large anonymous WiFi areas by using physically unlabeled highdimensional measurements; (ii) mobility map is constructed from significant LPs which are involved in many people’s motion paths; and (iii) people’s motion paths are tracked in an adequate precision level. To meet these goals, we divide our approach into the following four main steps: (i) measurement quantization in a lowdimensional manifold which lies in the raw RSS space, (ii) LP identification by SmithWaterman measurement matching [22], (iii) LP assembling into the mobility map in a temporal logic manner [23], and (iv) people’s motion path tracking in mobility map.
The rest of this paper is organized as follows. Section 2 gives an overview of some relevant tracking approaches which have been used so far. In Section 3, we describe the detailed steps involved in our proposed approach. Section 4 presents the experimental results and analysis. Finally, we conclude this paper and provide some future directions in Section 5.
2. Related Work
As the WiFi technique becomes prevalent wireless solution in public hotspots, there are a largely increasing number of different approaches used to track people’s paths by using WiFi technique [24, 25]. In general, these approaches fall into five main categories: proximity sensing, location fingerprinting, pattern matching, time trilateration, and angle triangulation.
2.1. Proximity Sensing
The proximity sensing is recognized as the simplest way to track people’s locations in a realtime manner [26, 27]. The location calculation is done based on the density of access points (APs) and granularity of divided cells in target area. In most cases, the target is located at the closest cell which it most probably belongs to. In [26], the authors divided the target area into several disjoint cells and fitted the Gaussian RSS distributions for the hearable base stations from the recorded RSSs in each cell. Then, when a localization request arrives, the Bayesian probabilistic method is employed to locate the target into the cell which has the highest confidence probability. Finally, the Markov chain is used for path tracking. As another example of the proximity sensingbased location tracking, the Herecast in [27] conducted the WiFi localization by using a database consisting of the APs’ service set identifiers (SSIDs) and the signal coverage range of each AP. For any location request, the area corresponding to the coverage of the AP which has been detected as the strongest AP, namely, the AP associated with the largest RSS, is referred to as the receiver’s estimated location. Based on this approach, it is extremely difficult to perform a finer path tracking due to the imprecise localization results.
2.2. Location Fingerprinting
The location fingerprinting has been most widely used in current location tracking systems in WiFi environments [28–30]. This approach requires the constructed radio map of fingerprints. Each fingerprint is a vector of RSS associated with its physical locations which are calibrated in the offline phase. In the online phase, the target or the location server retrieves the radio map to estimate the location which has the most similar fingerprint to each newly recorded RSS measurement. The first representative RADAR system [28] was designed based on the assumption that the physically adjacent locations have the same fingerprints as in signal space. The operation of RADAR system consists of two phases. The radio map is first constructed in the offline phase to be afterwards used for location estimation. In the online phase, the target’s locations are tracked by using the nearest neighbor(s) in signal space (or nearest neighbor(s) (KNN) algorithm). The Horus [30] and Nibble [29] are another two prominent fingerprintbased location tracking systems. Both the Horus and Nibble systems work based on the Bayesian inference approach, while the major difference between them is about the way to depict the RSS distributions at reference points (RPs). In Horus system, a Gaussian distribution curve for each hearable AP is fitted from the recorded RSSs at each RP, while the Nibble system uses a histogram to record the frequencies of recorded RSSs at each RP. Moreover, from the study of the problems about RSS correlation, variations of RSSs with respect to the environmental changes, and relations of RSSs and spatial characteristics, the Horus system is featured with high accuracy and low computation cost compared to the Nibble system.
2.3. Pattern Matching
Reference [31] proposed a new location tracking system, LENSR, adopting the nearest neighborhood vector mappingaided topological counter propagation network. Fang and Lin in [32] studied the discriminantadaptive neural network (DANN) for location tracking in WiFi environment. Different from the conventional pattern matching approaches, DANN extracts the lowdimensional discriminative components for neural network training. Other similar works on pattern matchingbased location tracking can be found in [33, 34]. The pattern matching approach addressed in [33] relies on the multilayer perceptron architecture by onestep secant training. In [34], the pattern matching approach with well training process is proved to perform better localization accuracy than the conventional nearest neighbor(s) and Bayesian inference approaches. However, the major drawback of pattern trainingbased location tracking system is that it should be conducted by sufficient training before it works.
2.4. Trilateration and Triangulation
The basic idea of trilateration and triangulation approaches comes from the time of arrival (TOA) and angle of arrival (AOA) measurements. To enable the localization in 2dimensional areas, the signal measurements from at least three and two APs should be made for the TOA and AOA systems, respectively [35, 36]. In TOA systems [35], the trilateration approach is conducted on the distances between the APs and tracking target which are calculated by the measured propagation time between them. Moreover, the exact time synchronization is also required for the measurement of propagation time. The main advantages of AOA systems [36] are that there are as few as two APs for the purpose of 2dimensional localization; meanwhile, the time synchronization between the APs and tracking target is not required. However, the location precision could degrade when the signal is blocked by the walls and infrastructures or the target is located far away from the APs.
Since the TOA and AOA location systems involve significant changes on hardware devices and infrastructures which make these two systems difficult to be widely applied in practice, the RSSbased trilateration approach is more preferred by current work [37]. Different from the TOA systems, the distances between the APs and tracking target are calculated by the RSS propagation models. In [37], Narzullaev compared three representative models used for WiFi RSSbased trilateration approach: (i) logdistance loss model which assumes that the mean of RSSs approximately decreases logarithmically with the propagation distance, (ii) multislope loss model which achieves a larger granularity of the predicted locations and requires a shorter sample collection time, and (iii) multiwall loss model which carefully takes the path loss caused by the walls and floors into account.
In all, applying the aforementioned location tracking approaches into the large WiFi environments could be a challenging work by the reasons of the inaccurate localization results in proximity sensing, laboring cost for fingerprint calibration and training process in location fingerprinting and pattern matching, respectively, and extra devices and infrastructures required by trilateration and triangulation approaches. The main contribution of this paper is to develop a better solution to track people’s motion path in large WiFi environments by using RSSbased timestamped mobility map without any fingerprint.
3. System Description
3.1. System Overview
Our proposed system consists of two phases: offline training phase and online tracking phase, as shown in Figure 1. The offline training phase is conducted on the network side with a large amount of computation resource, while the online tracking phase is conducted on the sourceweak client side.
In the offline training phase, we first record RSS measurements to conduct measurement quantization. The quantized RSS measurements are then used to identify the raw LPs by performing the SmithWaterman measurement matching. Each LP corresponds to a significant location which is involved in many people’s motion paths. Finally, we do the LP assembling to construct the mobility map corresponding to the target area. In the online tracking phase, we first quantize each new RSS measurement into a discrete level. Then, the matching LP with respect to each new online fragment can be determined based on the fine LP matching. Finally, people’s motion paths are tracked by connecting every two consecutive matching LPs along the shortest path in mobility map. For the sake of convenience, a list of notations used in this paper is given in Notation.
3.2. RSS Measurement Recording
In our system, the RSS measurements are sporadically recorded by our planned volunteers equipped with WiFi mobile receivers following their routine activities in target area. A measurement is a vector of RSS which consists of the RSS values from all the hearable APs [38]. Each string of consecutive measurements is called a fragment. As discussed in [28, 39], the measurements could be similar if they are recorded at nearby locations. We define the two fragments containing the common similar measurements as a growing fragment pair in which each overlapped piece of common similar measurements forms a raw LP to be afterwards used for LP merging and splitting to construct the timestamped mobility map.
We set as the th fragment where and stand for the numbers of fragments and measurement in , respectively, and is the th measurement. If there are hearable APs, we can obtain , where is the RSS value from AP . In each growing fragment pair (i.e., ), the th LP is denoted as . After all the LPs are obtained, the mobility map we seek to construct is recognized as a graph in which and stand for the sets of LPs and timestamped transition relations between neighboring LPs, as previously discussed in [21].
3.3. Measurement Quantization
For the sake of applying SmithWaterman measurement matching technique to construct mobility map, we need to quantize the RSS measurements into different discrete levels based on the similarities of RSS measurements. Specifically, we use Laplacian embeddingbased spectral clustering to quantize the RSS measurements which have been merged into the same cluster in the same quantization level. Thus, the number of clusters by spectral clustering equals the number of quantization levels. The detailed steps of measurement quantization process are provided as follows.
Step 1. We calculate the similarity of and as where .
Step 2. Considering the problem of mapping the raw measurements into a dimensional space, we can represent the mapped measurements as a matrix , where the superscript “” denotes the transpose operation anddenotes the mapped vector of . Based on the Laplacian embedding [40], we can obtain the optimal objective function as where “” denotes the trace operation, , and . When and , we have ; otherwise, we have . As discussed in [21], the solution to the optimization problem in (1) can be given by the eigenvectors associated with the smallest eigenvalues of the eigenvalue problem in (2):
Step 3. We perform means clustering on the mapped dimensional vectors to obtain the clusters, , where denotes the th cluster. Then, we quantize the RSS measurements corresponding to the mapped vectors in the same cluster into the same quantization level.
3.4. SmithWaterman Measurement Matching
The objective of SmithWaterman measurement matching is to identify the raw LPs for the construction of mobility map associated with the target area. To meet this goal, we adopt the SmithWaterman alignment approach to find the winning paths in the scoring space for each growing fragment pair and then perform measurement matching to identify the raw LPs. The steps of the raw LP identification are as follows.
Step 1. In growing fragment pair , , when and are in the same quantization level, we set a positive matching score, , for the measurement pair ; otherwise, we set a negative mismatching score, , for . The negative missing scores, and , are set when there is no measurement in and matched with and , respectively. We have the relations of “matching score missing score mismatching score.” Then, we can obtain the scoring space, , for the growing fragment pair , as shown in the following:
Step 2. We select the measurement pair, , which has the highest score in scoring space as the first point on the winning path, such that . We require that the score of the first point should be higher than (i.e., ), where is the threshold for measurement matching. In our experiments, we set .
Step 3. We compare the scores of three previous measurement pairs, , , and , and select the pair which has the highest score among them as the second point on the winning path. We repeat this process until the selected pair has the score zero. At this point, the selected pair with the score zero is defined as the last point on the winning path.
Step 4. After the winning path in scoring space is obtained, we identify the corresponding raw LP by merging the matched measurement pairs. Based on the SmithWaterman alignment, the three measurement matching criteria are provided as follows.(i)Criterion 1: measurements and are matched when there is a diagonal jump from to in scoring space.(ii)Criterion 2: measurement is not matched with any measurement in fragment when there is a topdown jump from to in scoring space.(iii)Criterion 3: measurement is not matched with any measurement in fragment when there is a leftright jump from to in scoring space.
To identify the other raw LPs from the scoring space, we continue to select the measurement pair which has the highest score in the remaining measurement pairs which are not involved in the previous winning paths as the first point of a new winning path. We follow Steps 3 and 4 until this new winning path arrives at a measurement pair which has the score zero or is involved in the previous winning paths. We name this measurement pair as the last point on this new winning path. For simplicity, we only focus on the situation that only one raw LP exists in a scoring space (i.e., for the growing fragment pair ) since the situation of multiple raw LPs can be avoided by manually chopping each longlength fragment into several shorter ones. The length of a fragment is defined as the number of measurements contained in this fragment.
3.5. Mobility Map Construction
After all the raw LPs have been identified, the next work is to assemble the raw LPs into the mobility map in a temporal logic manner. As discussed before, since the measurements in each raw LP are labeled by timestamps, we can approximately represent each raw LP as a time interval which starts at the last point and ends at the first point on its corresponding winning path. Then, the raw LP assembling process can be converted into a temporal reasoning problem, as introduced in [23]. The detailed steps are described below.
Step 1. Based on Allen’s interval algebra (i.e., 13 temporal logic relations: , , , , , , , , , , , , and ) in [23], we can capture the temporal logic relations between the raw LPs in each growing fragment pair. Specifically, when and are two raw LPs for the growing fragment pair , we obtain the following:
In (4), when the timestamp of the last point in a raw LP is larger than the timestamp of the first point of another raw LP, we define “the last point is after the first point”; otherwise, we define “the last point is before the first point.”
Step 2. The temporal logic relations between the raw LPs (i.e., ) and their belonging growing fragment pair (i.e., ) are given as follows: where the start and end points in (or ) are defined as the measurements which have the smallest and largest timestamps in (or ), respectively.
Step 3. Since the mobility map we seek to construct is a connected graph, the temporal logic relations of any two raw LPs can be obtained by Allen’s interval algebra based on the timestamped transitions between the LPs and fragments. To illustrate this result clearer, we use the transitivity table in [23] to show the temporal logic relations between the different raw LPs. Table 1 gives the possible temporal logic relations between any two LPs (i.e., and ) belonging to the two different growing fragment pairs (i.e., and ).
We take the relations of and , for instance. Based on the transitivity table, there are three possible temporal logic relations between and (i.e., , , and ), such that
Finally, the block diagram for the LP assembling into a mobility map is shown in Figure 2. We also take the relations of and , for instance. Based on (6) and Figure 2, (i) if the first points in and are the same (i.e., ), we merge into to form a new LP consisting of all the measurement pairs in and ; (ii) if the first point in is after the first point in (i.e., ), we merge all the overlapped measurement pairs in into and then delete all the overlapped measurement pairs in ; and (iii) if the first point in is after the first point in (i.e., ), we merge all the overlapped measurement pairs in into and then delete all the overlapped measurement pairs in .
 
“Noinfo” means that all the temporal logic relations are applied. 
3.6. Path Tracking in Mobility Map
There are two main steps involved in path tracking: (i) coarse RSS quantization and (ii) fine LP matching. The path tracking in mobility map is conducted as follows.
Step 1 (coarse RSS quantization). As discussed in Section 3.2, after the offline RSS measurement quantization, we can obtain clusters associated with the quantization levels. Then, for each new measurement, , in the online fragment, , where is the number of new measurements in , we calculate the Euclidean distance between and the average measurement in each cluster (i.e., ), , and then quantize the new measurement in a discrete level of cluster , such that
Step 2 (fine LP matching). We select the fragment which has the longest length in each LP as the labeling fragment. The labeling fragment in is denoted as , where is the number of measurements contained in . When the labeling fragment satisfies the relation of
we set as the matching LP. After that, all the new measurements before (or the new measurements with IDs not larger than ) are deleted to form a new online fragment (i.e., ) and then continue to search for the next matching LP. We repeat this process until there is no new measurement which remains in the online fragment or the score for the newly formed fragment is lower than the threshold, . In our experiments, we set . After all the matching LPs are obtained, we track the people’s motion paths by connecting every two consecutive matching LPs along the shortest path in mobility map. The shortest path is defined as the path which passes by the smallest number of LPs.
Some of the raw RSS fragments recorded may be very long. When the number of RSS measurements in a fragment is too long, the computation problem may arise for the process of LP assembling into a mobility map, while the main computation cost is involved in the offline training phase. In the online tracking phase, when the user sends a location query with its new RSS fragment, our system retrieves the cluster centers and returns the quantization level as well as the highest score in scoring space. The LP corresponding to the highest score is selected as the matching LP. At this point, the calculation complexity when our system tracks hundreds and thousands of people walking around in the target area forms an interesting work in future. To clearly show the computation cost required in offline training phase, we take the spectral clustering, for example. By using the MATLAB 7.10.0 (R2010a) under WINDOWS XP system, Table 2 compares the time cost for spectral clustering in different numbers of measurements conditions. All the computations are run on a PC with Intel Core i32120 CPU. In Table 2, we can find that as the number of measurements increases, the time cost for spectral clustering will also increase.

4. Experimental Results and Analysis
In this section, we will evaluate the performance of mobility map construction and motion path tracking based on the actual RSS fragments (of dimensions 650) recorded on five representative paths in HKUST campus. The five fragments are recorded on path 1, , which is from North Bus Stop to Library and with the length of 191; path 2, , which is from Lab 2149 to Library and with the length of 195; path 3, , which is from Lab 2149 to Coffee Shop and with the length of 85; path 4, , which is from Lab 2149 to Office 2514 and with the length of 106; and path 5, , which is from Lab 2149 to LTJ theater and with the length of 80 [41]. Each path consists of several physically adjacent traces. We take the fragment recorded on path 1 (i.e., ), for instance. contains two consecutive segments, and , which are recorded on trace 1 (between North Bus Stop and Atrium) and trace 2 (between Atrium and Library), respectively. The traces labeled by superscript “*” (i.e., trace 2, trace 6, and trace 7) are the LPs to be assembled in our testing. The notation “” indicates the path direction. A summary of these fragments and the corresponding traces is shown in Figure 3.
First of all, based on the raw measurement space consisting of 657 measurements with dimensions of 650 in Figure 4, we can calculate the similarity of any two measurements in Figure 5. The large similarity values (in range of ) represent that the corresponding measurement pairs are extremely similar. Moreover, the largest similarity value (or value 1) can be achieved by the similarity between any measurement and itself as expected.
As discussed in Section 3.2, the optimization problem in (1) can be converted into the generalized eigenvalue problem in (2). Figure 6 shows the three eigenvectors (of dimensions 657) associated with the smallest eigenvalues. Since the first eigenvector is a constant vector with the eigenvalue zero, we only use the second and third eigenvectors associated with the eigenvalues 0.91 and 0.94, respectively, as the basis of the mapped twodimensional space. Then, we obtain the 20 clusters (or quantization levels) in the mapped twodimensional space in Figure 7. For the SmithWaterman measurement matching, the 20 quantization levels can be recognized as the 20 amino acids. We denote the 20 amino acids as follows: Cys (C), Ser (S), Thr (T), Pro (P), Ala (A), Gly (G), Asn (N), Asp (D), Glu (E), Gln (Q), His (H), Arg (R), Lys (K), Met (M), Ile (I), Leu (L), Val (V), Phe (F), Tyr (Y), and Trp (W) [42]. After that, each RSS fragment can be represented by a string of amino acids, as shown in Table 3.

Second, after the raw measurements are quantized into 20 discrete levels, we will calculate the winning paths from the scoring space for each growing fragment pair. By setting , we can obtain four winning paths from the four growing fragment pairs , , , and in Figure 8. The matching score, missing score, and mismatching score are determined by PAM250 mutation matrix [42]. In Figure 9, the mobility map associated with the target area is constructed by LP assembling. From Figures 3 and 9, the mobility map constructed by LP assembling is extremely similar to the trajectory of the people’s actual motion. Then, the LP assembling is proved to perform well in layout description for the anonymous areas where the floor plans are not available.
(a) Fragment pair
(b) Fragment pair
(c) Fragment pair
(d) Fragment pair
Finally, to verify the efficiency of the mobility map constructed by LP assembling further, we use three paths for the testing of path tracking: (i) TP 1 Trace 1Trace 2 which is the same as path 1, (ii) TP 2 Trace 2Trace 3Trace 6Trace 7 which is in the opposite direction of path 2, and (iii) TP 3 Trace 5Trace 7Trace 7Trace 4 which is a combined path from paths 4 and 5. The previous two traces on TP 3 are in the opposite direction of path 5, while the following two traces are the same as path 4. In coarse RSS quantization step, we quantize the new fragments recorded on TP 1, TP 2, and TP 3 (with the lengths of , 195, and 186) into discrete levels, as shown in Figure 10.
(a) New measurement space
(b) Quantization results
Then, after the new fragment recorded on each path has been quantized, we conduct the fine LP matching to identify the matching LPs in mobility map. Without loss of generality, all the traces in mobility map (i.e., Traces 1–7) are considered for LP matching in our experiments. Then, the traces with the corresponding labeling fragments which satisfy the relations in (8) are selected as the matching LPs. The results of LP matching for the new online fragments recorded on TP 1, TP 2, and TP 3 are shown in Figures 11, 12, and 13, respectively.
(a) Scoring spaces for the raw new fragment and labeling fragments on different traces
(b) Scoring spaces for the 1st newly formed fragment and labeling fragments on different traces
(a) Scoring spaces for the raw new fragment and labeling fragments on the four similar traces
(b) Scoring spaces for the 1st newly formed fragment and labeling fragments on the four similar traces
(a) Scoring spaces for the raw new fragment and labeling fragments on the four similar traces
(b) Scoring spaces for the 1st newly formed fragment and labeling fragments on the four similar traces
(c) Scoring spaces for the 2nd newly formed fragment and labeling fragments on the four similar traces
In Figure 11(a), we find that the highest score, 61, is captured at the position (74, 104) in scoring space for the raw new fragment and the labeling fragment on Trace 1. Thus, we delete all the new measurements with the IDs not larger than from the raw new fragment to obtain the 1st newly formed fragment (with the measurement IDs from 105 to 191). As can be seen from Figure 11(b), Trace 2 with the highest score, 45.33, to the 1st newly formed fragment is selected as the second matching LP. Since the scores for any traces to the 2nd newly formed fragment are lower than , we construct the people’s motion path in mobility map as Trace 1Trace 2 which is the same as the actual motion path, TP 1.
Figures 12 and 13 show the scoring spaces for the new online fragments recorded on TP 2 and TP 3, respectively, and labeling fragments on the four most similar traces. The four traces which have the highest scores to the new online fragments are defined as the four most similar traces. The scores for the four most similar traces with respect to the new fragments are shown in Table 4.

Table 5 gives the new online fragments and the corresponding tracking paths which are connected by the matching LPs in mobility map. Therefore, our constructed timestamped mobility map performs well in people’s motion path tracking, especially for the motion paths which are in the same direction of the raw paths used for mobility map construction. For instance, on TPs 1 and 3, Trace 1Trace 2 and Trace 7Trace 4 which are in the same direction of the raw paths 1 and 4 are precisely tracked in mobility map in Figures 11 and 13.

5. Conclusions
In this paper, we introduced a novel WiFi RSS measurement matching approach to construct a finegranularity mobility map for layout description in anonymous areas without location fingerprinting. In our system framework, we first adopt SmithWaterman alignment approach to obtain the significant LPs from measurement matching in scoring space. Then, Allen’s interval algebra is used to conduct LP assembling into a timestamped mobility map in temporal logic manner. Our extensive mathematical analysis and experimental results show that our approach can deal with the problems of mobility map construction from sporadically recorded highdimensional WiFi RSS measurements and people’s motion path tracking in arbitrary directions. We will continue to investigate the integration of timestamped RSS measurements and motion sensors for the sake of constructing a highly precise mobility map. For instance, the relations of LPs can be described by not only the timestamps, but also the angle and speed measured by the offtheshelf smartphones.
Notation
:  th fragment 
:  th measurement in 
:  RSS value from the th AP in 
:  Growing fragment pair, and 
:  th LP in 
:  Mobility map 
:  Set of LPs 
:  Set of timestamped transition relations between neighboring LPs 
:  Similarity of and 
:  Normalized Euclidean distance of RSS between and 
:  Matrix of the mapped measurements 
:  Mapped vector of 
:  th mapped value in 
:  th eigenvalue value 
:  th cluster 
:  Positive matching score between and 
:  Negative mismatching score between and 
:  Negative missing score when there is no measurement in to be matched with 
:  Negative missing score when there is no measurement in to be matched with 
:  Scoring space with respect to 
:  Score between and 
:  New online fragment 
:  th new measurement in 
:  Average measurement in 
:  Euclidean distance between and 
:  Labeling fragment in 
:  Number of fragments 
:  Number of measurements in 
:  Number of hearable APs 
:  Number of LPs with respect to 
:  Number of dimensions of the mapping space 
:  Number of clusters/number of quantization levels 
:  Number of LPs after SmithWaterman measurement matching 
:  Number of new measurements in 
:  Number of measurements contained in 
:  Threshold for measurement matching 
:  Threshold for the score of the newly formed fragment 
:  Temporal logic relation between the LPs and 
, , , , , , , , , , , , and :  13 temporal logic relations in Allen’s interval algebra 
Superscript “”:  Transpose operation 
“”:  Trace operation. 
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors wish to thank the reviewers and editors for their valuable comments and useful corrections to improve this work. This work was supported in part by the Program for Changjiang Scholars and Innovative Research Team in University (IRT1299), the National Natural Science Foundation of China (61301126 and 61304197), the National Science and Technology Major Project (2012ZX03006002(3)), the Special Fund of Chongqing Key Laboratory (CSTC), the Fundamental and Frontier Research Project of Chongqing (cstc2013jcyjA40041, cstc2013jcyjA40032, and cstc2013jcyjA40034), the Scientific and Technological Research Program of Chongqing Municipal Education Commission (KJ130528), the Startup Foundation for Doctors of CQUPT (A201233), and the Science Foundation for Young Scientists of CQUPT (A201277).
References
 Z. Chen, J. Xia, and B. Irawan, “Development of fuzzy logic forecast models for locationbased parking finding services,” Mathematical Problems in Engineering, vol. 2013, Article ID 473471, 6 pages, 2013. View at: Publisher Site  Google Scholar
 Z. Tian, X. Liu, M. Zhou, and K. Xu, “Mobility tracking by fingerprintbased KNN/PF approach in cellular networks,” in Proceedings of the IEEE Wireless Communications and Networking Conference (WCNC '13), pp. 4570–4575, Shanghai, China, April 2013. View at: Publisher Site  Google Scholar
 J. Hightower and G. Borriello, “Location systems for ubiquitous computing,” IEEE Computer, vol. 34, no. 8, pp. 57–66, 2001. View at: Publisher Site  Google Scholar
 H. Wu, T. Wang, and Z. Zhong, “Location updating schemes for highspeed railway cellular communication systems,” Mathematical Problems in Engineering, vol. 2012, Article ID 802152, 15 pages, 2012. View at: Publisher Site  Google Scholar
 L. Zhang, F. Zhao, Z. Li, H. Luo, and J. Xu, “Design and implementation of WiFi network based mobile location service system,” in Proceedings of the 2nd International Conference on Business Computing and Global Informatization (BCGIN '12), pp. 641–644, Shanghai, China, October 2012. View at: Publisher Site  Google Scholar
 C. Evans, “Intelligent retail business: location based services for mobile customers,” in Proceedings of the 2nd International Conference on Pervasive Computing and Applications (ICPCA '07), pp. 354–359, Birmingham, UK, July 2007. View at: Publisher Site  Google Scholar
 T. R. Hansen, J. E. Bardram, and M. Soegaard, “Moving out of the lab: deploying pervasive technologies in a hospital,” IEEE Pervasive Computing, vol. 5, no. 3, pp. 24–31, 2006. View at: Publisher Site  Google Scholar
 F. van Diggelen, “Indoor GPS theory & implementation,” in Proceedings of the IEEE Position Location and Navigation Symposium, pp. 240–247, Palms Springs, Calif, USA, April 2002. View at: Publisher Site  Google Scholar
 M. Zhou, Q. Guo, and Z. Wang, “A novel stable clustering design method for hierarchical satellite network,” Chinese Journal of Aeronautics, vol. 23, no. 1, pp. 91–102, 2010. View at: Publisher Site  Google Scholar
 M. Zhou, Z. Tian, K. Xu, X. Yu, X. Hong, and H. Wu, “SCaNME: location tracking system in largescale campus WiFi environment using unlabeled mobility map,” Expert Systems with Applications, vol. 41, no. 7, pp. 3429–3443, 2014. View at: Publisher Site  Google Scholar
 M. Zhao, H. Yang, J. Liu, Y. Chen, and J. Zhou, “Directional WiFi based indoor location system for emergency,” in Proceedings of the 7th International Conference on Ubiquitous Intelligence & Computing/Autonomic & Trusted Computing, pp. 501–502, Xian, China, October 2010. View at: Publisher Site  Google Scholar
 M. Zhou, Z. Tian, K. Xu, X. Yu, and H. Wu, “Theoretical entropy assessment of fingerprintbased WiFi localization accuracy,” Expert Systems with Applications, vol. 40, no. 15, pp. 6136–6149, 2013. View at: Publisher Site  Google Scholar
 M. Zhou, Z. Tian, K. Xu, X. Hong, and X. Tang, “Towards statistical linear accuracy of neighbor matching localization in WiFi environment,” Journal of Computational Information Systems, vol. 9, no. 22, pp. 9177–9184, 2013. View at: Google Scholar
 K. Kaemarungsi and P. Krishnamurthy, “Properties of indoor received signal strength for WLAN location fingerprinting,” in Proceedings of the 1st Annual International Conference on Mobile and Ubiquitous Systems: Networking and Services (MOBIQUITOUS '04), pp. 14–23, Boston, Mass, USA, August 2004. View at: Publisher Site  Google Scholar
 S.H. Fang and T.N. Lin, “A dynamic system approach for radio location fingerprinting in wireless local area networks,” IEEE Transactions on Communications, vol. 58, no. 4, pp. 1020–1025, 2010. View at: Publisher Site  Google Scholar
 M. Zhou, Y. Xu, L. Ma, and S. Tian, “On the statistical errors of RADAR location sensor networks with builtin WiFi Gaussian linear fingerprints,” Sensors, vol. 12, no. 3, pp. 3605–3626, 2012. View at: Publisher Site  Google Scholar
 Z. Deng, Y. Xu, and L. Ma, “Indoor positioning via nonlinear discriminative feature extraction in wireless local area network,” Computer Communications, vol. 35, no. 6, pp. 738–747, 2012. View at: Publisher Site  Google Scholar
 M. Zhou, Y. Xu, and L. Tang, “Multilayer ANN indoor location system with area division in WLAN environment,” Journal of Systems Engineering and Electronics, vol. 21, no. 5, pp. 914–926, 2010. View at: Publisher Site  Google Scholar
 C. Koweerawong, K. Wipusitwarakun, and K. Kaemarungsi, “Indoor localization improvement via adaptive RSS fingerprinting database,” in Proceedings of the International Conference on Information Networking (ICOIN '13), pp. 412–416, Bangkok, Thailand, January 2013. View at: Publisher Site  Google Scholar
 X. Wang, A. K. Wong, and Y. Kong, “Mobility tracking using GPS, WiFi and cell ID,” in Proceedings of the the International Conference on Information Networking (ICOIN '12), pp. 171–176, Phuket, Thailand, February 2012. View at: Publisher Site  Google Scholar
 M. Zhou, A. K. Wong, Z. Tian, V. Y. Zhang, X. Yu, and X. Luo, “Adaptive mobility mapping for people tracking using unlabelled WiFi shotgun reads,” IEEE Communications Letters, vol. 17, no. 1, pp. 87–90, 2013. View at: Google Scholar
 T. F. Smith and M. S. Waterman, “Identification of common molecular subsequences,” Journal of Molecular Biology, vol. 147, no. 1, pp. 195–197, 1981. View at: Google Scholar
 J. F. Allen, “Maintaining knowledge about temporal intervals,” Communications of the ACM, vol. 26, no. 11, pp. 832–843, 1983. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. Zhou, Y. Xu, and L. Ma, “Radiomap establishment based on fuzzy clustering for WLAN hybrid KNN/ANN indoor positioning,” China Communications, vol. 7, no. 3, pp. 64–80, 2010. View at: Google Scholar
 Z. Tian, X. Tang, M. Zhou, and Z. Tan, “Fingerprint indoor positioning algorithm based on affinity propagation clustering,” EURASIP Journal on Wireless Communications and Networking, vol. 2013, article 272, 2013. View at: Publisher Site  Google Scholar
 E. Flannery, A. Haeberlen, L. Kavraki, A. Ladd, A. Rudys, and D. Wallach, “System and method for localization over a wireless network,” Patent WO 2006017266 A2, 2006. View at: Google Scholar
 M. Paciga and H. Lutfiyya, “Herecast: an open infrastructure for locationbased services using WiFi,” in Proceedings of the IEEE International Conference on Wireless and Mobile Computing, Networking and Communications (WiMob '05), vol. 4, pp. 21–28, Montreal, Canada, August 2005. View at: Publisher Site  Google Scholar
 P. Bahl and V. N. Padmanabhan, “RADAR: an inbuilding RFbased user location and tracking system,” in Proceedings of the 19th IEEE Annual Joint Conference of the IEEE Computer and Communications Societies (IEEE INFOCOM '00), pp. 775–784, TelAviv, Israel, March 2000. View at: Google Scholar
 P. Castro, P. Chiu, T. Kremenek, and R. Muntz, “A probabilistic room location service for wireless networked environments,” in Proceedings of the 3rd International Conference on Ubiquitous Computing (UbiComp '01), pp. 18–34, Atlanta, Ga, USA, March 2001. View at: Google Scholar
 M. Youssef and A. Agrawala, “The Horus location determination system,” Wireless Networks, vol. 14, no. 3, pp. 357–374, 2008. View at: Publisher Site  Google Scholar
 D. Kurt and M. Milos, “Wireless based object tracking based on neural networks,” in Proceedings of the 3rd IEEE Conference on Industrial Electronics and Applications (ICIEA '08), pp. 308–313, Singapore, June 2008. View at: Publisher Site  Google Scholar
 S.H. Fang and T.N. Lin, “Indoor location system based on discriminantadaptive neural network in IEEE 802.11 environments,” IEEE Transactions on Neural Networks, vol. 19, no. 11, pp. 1973–1978, 2008. View at: Publisher Site  Google Scholar
 R. Battiti, T. L. Nhat, and A. Villani, “Locationaware computing: a neural network model for determining location in wireless LANs,” Tech. Rep. DIT020083, Department of Information and Communication Technology, University of Trento, Trento, Italy, 2002. View at: Google Scholar
 S. Saha, K. Chaudhuri, D. Sanghi, and P. Bhagwat, “Location determination of a mobile device using IEEE 802.11b access point signals,” in Proceedings of the IEEE Wireless Communications and Networking (WCNC '03), vol. 3, pp. 1987–1992, New Orleans, La, USA, March 2003. View at: Publisher Site  Google Scholar
 S. A. Golden and S. S. Bateman, “Sensor measurements for WiFi location with emphasis on timeofarrival ranging,” IEEE Transactions on Mobile Computing, vol. 6, no. 10, pp. 1185–1198, 2007. View at: Publisher Site  Google Scholar
 A. V. Alejos, M. G. Sanchez, D. R. Iglesias, and I. Cuinas, “Incidence angle estimation algorithm for a lowcost AoA BASED INDOOR location system,” in Proceedings of the IEEE International Symposium on Antennas and Propagation and USNC/URSI National Radio Science Meeting (APSURSI '09), pp. 1–4, Charleston, SC, USA, June 2009. View at: Publisher Site  Google Scholar
 A. Narzullaev, Y. W. Park, and H. Y. Jung, “Accurate signal strength prediction based positioning for indoor WLAN systems,” in Proceedings of the IEEE/ION Position, Location and Navigation Symposium (PLANS '28), pp. 685–688, Monterey, Calif, USA, May 2008. View at: Publisher Site  Google Scholar
 M. Zhou, Z. Tian, K. Xu, X. Yu, and H. Wu, “Error analysis for RADAR neighbor matching localization in linear logarithmic strengthvarying WiFi environment,” The Scientific World Journal, vol. 2014, Article ID 647370, 15 pages, 2014. View at: Publisher Site  Google Scholar
 H. Alasti, K. Xu, and Z. Dang, “Efficient experimental path loss exponent measurement for uniformly attenuated indoor radio channels,” in Proceedings of the 10th IEEE Southeast Conference (SOUTHEASTCON '09), pp. 255–260, Atlanta, Ga, USA, March 2009. View at: Publisher Site  Google Scholar
 M. Belkin and P. Niyogi, “Laplacian eigenmaps for dimensionality reduction and data representation,” Neural Computation, vol. 15, no. 6, pp. 1373–1396, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 http://pathadvisor.ust.hk.
 M. O. Dayhoff, R. M. Schwartz, and B. C. Orcutt, “A model of evolutionary change in proteins,” Atlas of Protein Sequence and Structure, vol. 5, supplement 3, pp. 345–351, 1978. View at: Google Scholar
Copyright
Copyright © 2014 Mu Zhou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.