Research Article  Open Access
A ThreeDimensional Wireless Indoor Localization System
Abstract
Indoor localization, an emerging technology in location based service (LBS), is now playing a more and more important role both in commercial and in civilian industry. Global position system (GPS) is the most popular solution in outdoor localization field, and the accuracy is around 10 meter error in positioning. However, with complex obstacles in buildings, problems rise in the “last mile” of localization field, which encourage a momentum of indoor localization. The traditional indoor localization system is either rangebased or fingerprintingbased, which requires a lot of time and efforts to do the predeployment. In this paper, we present a 3dimensional ondemand indoor localization system (3DODIL), which can be fingerprintfree and deployed rapidly in a multistorey building. The 3DODIL consists of two phases, vertical localization and horizontal localization. On vertical direction, we propose multistorey differential (MSD) algorithm and implement it to fulfill the vertical localization, which can greatly reduce the number of anchors deployed. We use enhanced field division (EFD) algorithm to conduct the horizontal localization. EFD algorithm is a rangefree algorithm, the main idea of which is to dynamically divide the field within different signature area and position the target. The accuracy and performance have been validated through our extensive analysis and systematic experiments.
1. Introduction
The next evolutionary in smart environments develops fast in building, utilities, industrial, home, shipboard, and transportation system automation. One of the crucial technologies and scientific researches of smart environments is based on locationaware systems. Locationaware systems have improved the quality of life and provide people with convenience in working, outing, shopping, dining, and many other activities.
In recent years, with the opensourced Android operating system and the continuous costreduction in hardware, smart phone is getting more and more popular. Android/IOS based smart phones can be seen everywhere, and we can easily get the location of our friends with the widely used locationaware systems, like global position system (GPS). GPS is the most popular solution in outdoor localization field, and the accuracy is around 10 meter error in positioning. It provides good userexperience in outdoor scenarios. However, when we try to use GPS in indoor environment, like in an office or at home, we can barely receive the signals of satellites, and the error of accuracy will be out of tolerance. The root cause is that with complex obstacles in buildings the satellite signals cannot transmit the walls and floors.
Thus, we need an indoor localization system to do the “last mile” localization. The “last mile” or “last kilometer” is a phrase used by the telecommunications, cable television, and Internet industries to refer to the final leg of the telecommunications networks delivering communications connectivity to retail customers. In localization field, we also have the “last mile” problems, and we need the indoor localization systems that can make up the weakness of GPS system.
Over the past decade, indoorlocalizationbased technology and algorithms are flourishing, and almost all indoor localization systems are based on signal analysis of wireless sensor networks (WSNs). Generally speaking, indoor location mechanisms can be separated into two parts: rangebased method and rangefree method. Rangebased methods are the mainstream ways to do indoor location in the early years, such as TOA [1], TDOA [2], AOA [3], and RSSI [4]. Since the rangebased method requires expensive hardware and precise statistics, the rangefree algorithms develop fast in recent years, such as RSD, EZ, BT, and UnLoc. However, some limitations restrain these algorithms from generalized indoor location, for fingerprinting technology needs lots of precollection of data.
Thus, we propose the idea of 3dimensional ondemand indoor localization system (3DODIL), which can be fingerprintfree and deployed effectively in a multistorey building. The contributions are as follows.
In this paper, we propose the idea of 3dimensional ondemand indoor localization system (3DODIL), which can be fingerprintfree and deployed rapidly in a multistorey building. The 3DDDIL consists of two parts, horizontal localization and vertical localization.
On horizontal direction, we use enhanced field division (EFD) algorithm to conduct the horizontal localization. EFD algorithm is a rangefree algorithm, which implements the idea that the whole map is dynamically divided into fields that with signature sequences mark and locate the target device. The key idea of EFD provides a new developing direction for indoor localization. Traditional ways focus on the actual value of the received signal strength (RSS), but the value of which is fluctuant and greatly impacted due to the complex environments with obstacles. The EFD algorithm lays emphasis on the possible position on the map and utilizes the highlow relationship between the RSSs from each anchor to decide which areas the receiver should lies in. This allows for the minimization of the RF fading effect. What is more, taking advantage of field division with optimal anchor placements, we avoid the issues of intensive fingerprint sampling. Fingerprinting technology is a highly precise way to do indoor localization, but the weakness is also obvious: it need more time to find optimal anchor placements.
On vertical direction, we implement multistorey differential (MSD) algorithm to fulfill the vertical localization, which can greatly reduce the number of anchors deployed. Vertical localization, or called multistorey positioning, is an underestimated aspect to be discussed in indoor localization field. With anchors predeployed in the building in typical system, we can equip anchors with the same allocation plan on each floor, and it is simple to determine the floor: the floor of the anchor which has the strongest signal. However, in ondemand situation, we need to do rapid deployment, as few anchors as possible. We need to divide more floors with limited anchors. The MSD algorithm provides an efficient way to do vertical localization with few anchors.
Finally, we build a localization system platform, as shown in Figure 7, from offtheshelf commercialized products. The system includes three major components: the mobile anchor platform, a mobile localization unit, and a backend server. The mobile anchor platform can be a portable LTE/GSM miniradio station or a traditional WiFi router, which act as a radio signal generator. The mobile localization unit is an Androidbased smart phone, which is a commercial device and we can use it to receive signals and communicate with the server. The backend server plays the role of data storage and computing unit.
The rest of the paper is organized as follows. Multistorey differential algorithm is presented in Section 2. Enhanced field division (EFD) algorithm is discussed in Section 3. Analysis data of MSD algorithm is illustrated in Section 4. The performance of the prototype system in practical environments is evaluated in Section 5. Section 6 briefly discusses related work, and Section 7 concludes the paper.
2. Multistorey Differential Algorithm
Threedimensional localization system should provide both horizontal and vertical localization functions as a comprehensive system. In contrast to horizontal location, vertical localization has quite different model. The insertion loss of the reinforced wall (about 20–40 dB) is larger comparing to the brick wall (about 10 dB to 20 dB) [5]. Since the separation of reinforced concrete walls between all floors, it is inaccurate to use typical trilateral RSSI localization system to do vertical localization. Wireless indoor localization without site survey (WILL) [6] provides a fingerprinting system to do 3dimensional localization. Unfortunately, there is no mature and general vertical localization which can be rapidly deployed ondemand. Thus, we designed multistorey differential (MSD) algorithm.
2.1. Fundamental Theory of Vertical Localization
Figure 1 shows the basic building model we use. Assume that the height of the building is , and the height of each floor is . We set up two anchors on the first floor () and top floor (), which are at the same horizontal position. Since we focus on the vertical localization, we can simplify the 3dimensional system of coordinate to 2dimensional system of coordinate, of which axis is the vector from to and axis is the distance between our test point and axis . The position of the test point () is .
Definition 1 (FSPL). FSPL is the abbreviation of freespace path loss, which is the difference between transmitter power and received signal power, and also takes the gain of two antennas into consideration: where is the received signal power, is the base station (BS) transmitter power, is the transmitting antenna gain, and is the receiving antenna gain.
Definition 2 (MSPL). MSPL is the abbreviation of multistorey path loss, which is the extension of FSPL with the insertion loss of reinforced concrete wall: where is the number of walls inserted in the path and is the insertion loss of reinforced concrete wall (about 20–40 dm) [5].
Theorem 3. Consider
Proof of Theorem 3. Corresponding values of path losses were computed from the measured received signal power using the expression [7]
where EIRP is the effective isotropic power of the BS transmitting antenna. For the freespace model, the expression for the received signal power, with directional base station antenna, is [8, 9]
where is the wavelength and is the distance between the BS antenna and mobile station (MS) receiver antenna. If and are expressed in dBm and the gains are in dB, (2) becomes [10]
Consider the FSPL, which is defined in (1). When substituting , (3) and (4) become
Also, we can get the MSPL equation
For typical radio applications, it is common to find measured in units of MHz and in meter. Assuming that we use WiFi routers as the anchors ( MHz), the MSPL equation becomes where is a constant value, and we can calculate the exact value as 40.05 in our assumption.
As we can see, equals the floor integer of , which is . At test point, we can receive the signals from both and . The MSPL for is
And the for is
Considering the variant of and ,
In our scenario, our target is to conduct localization in multistorey building, and we should combine our math equation with the given conditions. Thus, we focus on the exact floor rather than the precise height.
Theorem 4. When we are on different floors of a multistorey building, we can distinguish the value regardless of the horizontal point we stay at.
Proof of Theorem 4. Since we are aiming to distinguish the floor, we assume that the vertical coordinates are evenly spaced, and the absolute height is , where is the floor number. To prove our theorem, we receive the at two different points, which are noted as and :
Consider the relationship between and . Obviously, if we can distinguish different points of the floor and floor , we can easily distinguish each floor of the building. Thus, we need to make research on the extreme condition that , and we generalize as . Thus, (13) becomes
If we want to use this theorem to localize the floor, we need to decrease the impact of . We separate it into two parts:
.
We name from . Equation (14) becomes
Take as variable; we can get the derivative of the function:
When , the derivative , which means that is a monotonic decreasing function with variable when . When , (15) becomes
When , (15) becomes
.
We assume that . If we need to evaluate the impact of variables , , we can analyze the extreme condition: , . Equation (14) becomes
Equations (17), (18), and (19) demonstrate the different condition of , the first item of them has minor impact (about 0–10 dB), and the second item has a major impact (about 40–50 dB).
With the deduction above, we can use as the reference value to distinguish different floors.
2.2. Modeling with Multiple Anchors on Each Floor
We can separately do the vertical localization with the value of theoretically, but, in real conditions, with the complex indoor environment and the fluctuation feature of wireless signal, we need to design a stable and multianchor model to enhance the accuracy of vertical localization.
The fundamental theory of vertical localization analyzes the condition of 1pairanchor deployment, and we want to generalize it to multianchor model.
Consider the actual deployment scenario. We need to deploy wireless anchors in the whole building. For floor , we put anchors on it.
As we name our system ondemand, we should minimize the number of anchors deployed with less impact on the accuracy of the whole system. With the fundamental theory of vertical localization, we can distribute the anchors on some specific floors, which have an interval of floors: The of each floor is decided by the anchors needed for horizontal localization system.
If a building has more storeys than , we can degrade this building to several subbuildings with storeys and filter the receiving wireless signal only from floor (the first floor of subbuilding) and floor .
Now we have a storey subbuilding model, and we need to decide the interval and figure out the anchor algorithm.
The interval is decided by the value that how many anchors on th floor can be detected on the th floor. We note the value as . Obviously, must be a nonnegative number, and is the minimum of . We name the threshold of the subbuilding model, which equals .
To do the vertical localization, we need to guarantee that . Referring to (3) and (8) and substituting : Consider the normal usage conditions: dBm, dBm, and dBm. The smart phone can detect the wireless signals that is more than −99 dBm. Thus, If m and dB, we can get the .
When we receive wireless signals on the th floor, we need to take advantage of all received signals to enhance the accuracy and the success rate of vertical localization.
We assume that the value of on each floor is in a narrow interval, for we can separate the value from different floors. Let us define the weight assigned to different anchors as ; is the weight for anchor . As such, the value of multianchor mode will be determined as
Now we should consider the determination of : where is defined as the average of the signal, is the standard deviation of , and defines a relative interference coefficient that describes the loss in signal strength due to environmental interference. The will be higher for a signal unfriendly environment such as a concrete wall building than for an outdoor environment known for less interface. Therefore, to improve the accuracy, a signature must be selected. Let us define the weight assigned to different anchors as ; is the weight for . As such, the optimization problem will be to determine so that the total utility is maximized; that is,
Lemma 5. The optimal solution of will be
Proof of (26). Under our definition, the total utility can be written as in which, is the degradation correlation between anchors and , and is the number of anchors. The numerical solution of can be solved by Lagrangian derivation. Let us form where is a Lagrange multiplier. Then, taking derivatives of the Lagrangian with respect to vector , we get As a result, this proves (26).
The optimal can then be applied to our measurement of as to derive the optimal field signature .
3. The Algorithm of Enhanced Field Division
3.1. Field Division
During our previous research, we find that the relative signal strength rather than the absolute one can be employed to mark different unique regions in an environment. We then introduce the concept of field division with each divided subarea owning a unique signature (which will be described later). Without loss of generality, we start the discussion within a twodimensional paradigm although our algorithm can be applied into a threedimensional paradigm. We focus on how to use relative signal strength to determine subarea. For convenience, we first take a look at one region without interference. (Note that the practical factors will be considered in detail later on.) Ignoring the environmental noise, the received signal strength (RSS) will demonstrate a monotonic effect on the geographic distance. Though using the relative signal strength, we consider that a RSS model can facilitate our discussion. A wellstudied RSS distribution model in db (decibel) as in (31) [11] meets the purpose: in which, , , and are the constants parameters for anchor and is the distance between the current location and anchor . Letting be the number of anchor nodes in a map, we define a highdimensional location signature function as follows: where is signature of the target location and , , is the target’s RSSs relative to the anchor node . The function is modeled by a descending sorting function. Therefore, we have for all for every signature . Given the modified RSS transmission model, we denote the subarea in Definition 6 for field division.
Definition 6 (subarea). Subarea is the point set whose members own the same signature in the field segmentation.
And the boundary is defined as a curve that is separating two adjacent subareas. Given any paired anchors, , , the dedicated part of the boundary attributed to them consists of positions, where . Taking a region with two anchors, and , as an example, the positions on the boundary can be derived as where and represent the relative distance to anchors and . Solving this equation, we can get a relationship of and as , where .
3.2. Model Adjustment
Note that it is not the absolute accuracy of the RSSs, but rather the relative strength that EFD algorithm relies on in localization. Based on our measurements, we discovered the RSS degradation due to the impact of interference toward different anchors to be highly correlated. To accommodate RF fading and the multipath effects in environments, we can make an adjustment to the RSS model in (31) for the sake of convenience. Therefore, we model the effects by a factor to so that absorbs the environmental interference pattern. The ideal RSS model will be adjusted as where is the interference factor of anchors and is the adjustment error. Because factors , are highly correlated, the huge interference does not actually make that much difference in the division signatures after adjusting the model. As a result, the problem turns into estimating to minimize the estimation . If we denote and assume that the follows a normal distribution, , and there are prior beliefs about the and , and represent the best guess for and the uncertainty of the guess, which comes from prior experiments or specification from service providers; we can modify the inputs to (34) based on the in situ estimation of and [12]. A standard estimator of expectations is the sample mean : where is the number of available sampling series. However, the sample mean is a highly inefficient estimator as the sampling estimation varies widely when different sampling series are fed into the estimation process [12]. One way to cope with this issue is to use a more efficient balance estimator: where is our best guess, , and is the balance factor. The purpose is then to minimize balance under any given . We can formalize our model as where is covariance matrix between and and is the vector .
The first solution to that formalization is
Thus, the optimal balance factors given by the answer are used to derive the , as the optimal representation of in (34).
3.3. Grid Tracking Strategy
As mentioned previously, one of the biggest challenges that the EFD system faces is the status quo issue. In this section, we will discuss another key technique used in EFD, the grid tracking strategy, in which ordinal gridding based on in situ information such as the magnetic field is used to assist localization within a subarea. Denote by the estimated velocity of the target at time and as the moving direction of the target. These two properties, importantly describing the motion characteristics, can be estimated precisely by approaches such as the geomagnetic field analysis by the EFD location platform.
3.3.1. Magnetic Theory
Geomagnetic field is relatively stable and easy to use in giving a direction, so we estimate the directions via the geomagnetic field. To preprocess the magnetic field data, we adapt a recursive method to track the dynamic measurement by assuming that it is approximately stationary. Denote by the measured signal at time slot , which corresponds to , , or component of the magnetic field data and may suffer from the measurement noise or other uncertainties. The recursive method intends to track the noisefree data by updating where denotes the estimated signal at time slot and is a weight factor that balances between earlier data prior to time and the current measurement . Clearly, the larger is, the more dependent is on earlier measurements. By this way, it is possible to remove the noise effects by relying on the historical data. However, there is a clear tradeoff on the value of , since a larger would introduce difficulties in capturing sudden changes in . To this end, the smoothed magnetic fields with various choices of are plotted in Figure 2 to illustrate the tradeoff in selecting an appropriate value for . Lastly, by initializing , the recursive update in (39) is extremely simple to implement since it involves only linear updates. Interestingly, this preprocessing step can effectively smooth out the perturbation present in the magnetic field data, as shown in Figure 2 with .
(a)
(b)
(c)
(d)
3.3.2. Details of the Grid Tracking
While the target is moving, fictitious edges for four twodimensional directions, S, N, W, and E, are dynamically created to gauge the motion range of a target at each cycle . By default, the edges matrix is set to the boundary () of the subarea : where is the coordinates of the subarea when the target enters initially. Then, the at time is updated by where is a matrix; for example, and limit to the north.
If the edge exceeds the boundary of subarea , this will be replaced with the boundary. Note that the process continues until a detection of a signature change, which indicates that the target is crossing the boundary of the subarea. To move the , we define a velocity matrix as
If a target moves within one subarea, . However, once it crosses the boundary, the limit will make a specific change.
3.4. Correction with Boundary Passing
As the target crosses the boundary separating two subareas and , the corresponding signatures from EFD’s mobile localization platform will change. This actually offers a great opportunity for EFD to calibrate the location estimation. Denote by the direction indicator as in which , and . As the signatures change, the correction algorithm will be carried out immediately as described in Algorithm 1.
4. System Analysis
To illustrate the performance of the MSD system, we executed extensive computer analysis and conducted localization tests based on our softwarehardware 3DODDIL system platform. In this section, the analysis results are demonstrated.
4.1. MSD with One Anchor Pair
First, we design the analytical environment: we plot 1 pair of anchors on the 1st and th floor and calculate the values under this circumstance. The anchor pair is at the same point vertically, and we evaluate the value when the horizontal distance varies from 0 to 20 meter. We repeat the computing with different values to see whether our algorithm is universal. And the results are shown in Figure 3. We define the floor height meters and the insertion loss of reinforced concrete wall dB. The max distance from test point to the anchors horizontally is 20 meter.
(a) Interval
(b) Interval
(c) Interval
(d) Interval
4.2. MSD with Multipair Anchors
In this section, we analyze the condition that there are multianchors in MSD algorithm. The test environment is an expanded condition of MSD with 1 pair of anchors. Figure 4 illustrates different conditions with 3 anchor pairs. Figure 5 illustrates different conditions with 5 anchor pairs. We define as the horizontal distance for anchor pair . We set interval . The figures illustrate the value curve when multipair anchors are detected.
(a) d1 = 0, d2 = 10, and d3 = 20
(b) d1 = 5, d2 = 10, and d3 = 20
(c) d1 = 10, d2 = 15, and d3 = 20
(d) d1 = 12, d2 = 15, and d3 = 18
(a) d1 = 0, d2 = 5, d3 = 10, d4 = 15, and d5 = 20
(b) d1 = 6, d2 = 8, d3 = 16, d4 = 18, and d5 = 18
4.3. Effects of Multipair Anchors
As we can see in the results, the analytical points on each floor fluctuate in narrow space, and the line of each floor is almost parallel to each other. When more anchors are detected, the result curve is smoother and improves to distinguish the floors. We can conclude that, under ideal circumstance, the MSD algorithm can fulfill the task of floor localization.
5. Experiment and Analysis
To quantify the performance of the proposed MSD algorithm, we carried out extensive vertical location tests by using our developed localization system platform. In this section, we categorize the experiments we have done into two sets in different environments with increasing complexity.
5.1. Localization System Platform
As a generic localization system, the 3DODDL systems can be suitable to any types of anchors with different radio signal frequencies. For the convenience deployment of experiments, we build our platform with WiFi routers as the anchors and Androidbased smart phones/tablets as the signal receivers. We choose Tenda W3000R routers as our anchors and use Google Nexus 7 as the signal receiver. The whole platform is shown in Figure 6.
(a)
(b)
5.2. Experiments in a Dormitory Building
5.2.1. Received Signals in the Building
Figure 8 is the signal received map when we jog on the first floor. We set one anchor on each floor from floor 1 to floor 3. All the anchors are on the same horizontal position. We can see that the signals of each anchor fluctuate as time flies.
(a)
(b)
(c)
(d)
(e)
(f)
5.2.2. MSD Experiments with Two Anchor Pairs
In this part, we set two anchor pairs on floor 1 and floor 3 of the building and receive signals on floors 1, 2, and 3. The distance between each anchor pair on the same floor is 10 meters. The sampling points on horizontal direction are 2 m, 4 m, 6 m, 8 m, and 10 m. We record 20 times for all sampling points with every 5 seconds. Figure 8 shows the test data, with the coordinates: axis sampling time in sec, axis distance on horizontal direction in meter, and axis received signals in dB.
Then, we use our weighting model to calculate the weightedaverage of received signals on each floor: Finally, we use the MSD algorithm to calculate the values of . We can see the key feature in Table 1.

6. Related Work
Localization in WSN has become an active research topic again recently, and a continuum of algorithms has been proposed [13–25]. Sensors and beacons are utilized to detect the movements for rangebased localization in [14, 16, 22, 23]. To undermine the effect from physical channel, rangefree algorithms based on a relative distance using trilateration and gridbased methods based on landmarks in school and office contexts emerged in RSD [26], RND [27], UnLoc [28], and EZ [29]. Both rangefree and gridbased methods make great improvement in roubustness and reliability, and each of them owns unique merits in indoor localization. However, they understate the complexity of implementing 3dimensional localization, albeit they depict reliable methods for 2dimensional localization. In this paper, we focus on the 3dimensional localization, analyze the characteristics of indoor vertical distribution, and propose a distinctive vertical localization and a horizontal localization in a hybrid idea of rangefree and gridbased methods.
7. Conclusion
This paper introduces a new mechanism of localization called 3dimensional ondemand indoor localization system (3DODIL), which can increase the accuracy and stability of localization of multistorey buildings. On horizontal direction, we use enhanced field division (EFD) algorithm to conduct the horizontal localization. On vertical direction, we implement multistorey differential (MSD) algorithm to fulfill the vertical localization, which can greatly reduce the number of anchors deployed. To test the mechanism of the MSD, we conduct a series of experiments. And we also build a localization system platform to conduct realenvironmental experiments. The result proves that the MSD is more reliable than other approaches for vertical localization. More importantly, the lowcost design and ondemand deployment allow the 3DODIL for largescale applications and wide utilization.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research was supported in part by the National Key Basic Research Program of China (2013CB329603), the National Natural Science Foundation of China (nos. 61271220 and 61170164), and NSF Grant CNS1217791.
References
 X. Wang, Z. Wang, and B. O'Dea, “A TOAbased location algorithm reducing the errors due to nonlineofsight (NLOS) propagation,” IEEE Transactions on Vehicular Technology, vol. 52, no. 1, pp. 112–116, 2003. View at: Publisher Site  Google Scholar
 H. Ni, G. Ren, and Y. Chang, “A TDOA location scheme in OFDM based WMANs,” IEEE Transactions on Consumer Electronics, vol. 54, no. 3, pp. 1017–1021, 2008. View at: Publisher Site  Google Scholar
 D. Niculescu and B. Nath, “Ad hoc positioning system (APS) using AOA,” in Proceedings of the 22nd Annual Joint Conference on the IEEE Computer and Communications Societies (INFOCOM '03), vol. 3, pp. 1734–1743, San Francisco, Calif, USA, April 2003. View at: Publisher Site  Google Scholar
 F. Viani, L. Lizzi, P. Rocca, M. Benedetti, M. Donelli, and A. Massa, “Object tracking through RSSI measurements in wireless sensor networks,” Electronics Letters, vol. 44, no. 10, pp. 653–654, 2008. View at: Publisher Site  Google Scholar
 T. B. Gibson and D. C. Jenn, “Prediction and measurement of wall insertion loss,” IEEE Transactions on Antennas and Propagation, vol. 47, no. 1, pp. 55–57, 1999. View at: Publisher Site  Google Scholar
 C. Wu, Z. Yang, Y. Liu, and W. Xi, “WILL: wireless indoor localization without site survey,” IEEE Transactions on Parallel and Distributed Systems, vol. 24, no. 4, pp. 839–848, 2013. View at: Publisher Site  Google Scholar
 J. Walfisch and H. L. Bertoni, “A theoretical model of UHF propagation in urban environments,” IEEE Transactions on Antennas and Propagation, vol. 36, no. 12, pp. 1788–1796, 1988. View at: Publisher Site  Google Scholar
 V. Garg and J. Wilkes, Principles and Applications of GSM, Prentice Hall, Upper Saddle River, NJ, USA, 1999.
 M. Yacoub, Propagation Considerations for Mobile Radio System, Brazil Publication, University of Campinas, Campinas, Brazil, 1996.
 K. Feher, Wireless Digital Communication: Modulation and Spread Spectrum Applications, Prentice Hall, Old Tappan, NJ, USA, 1995.
 C. Papamanthou, F. P. Preparata, and R. Tamassia, “Algorithms for location estimation based on RSSI sampling,” in Algorithmic Aspects of Wireless Sensor Networks, vol. 5389 of Lecture Notes in Computer Science, pp. 72–86, Springer, Berlin, Germany, 2008. View at: Google Scholar
 Q. Zhang, G. Sobelman, and T. He, “Gradientdriven target acquisition in mobile wireless sensor networks,” in Mobile AdHoc and Sensor Networks, vol. 4325 of Lecture Notes in Computer Science, pp. 365–376, Springer, Berlin, Germany, 2006. View at: Publisher Site  Google Scholar
 C. Zhang, Y. Zhang, and Y. Fang, “Localized algorithms for coverage boundary detection in wireless sensor networks,” Wireless Networks, vol. 15, no. 1, pp. 3–20, 2009. View at: Publisher Site  Google Scholar
 C. Wang, J. Chen, Y. Sun, and X. Shen, “Wireless sensor networks localization with isomap,” in Proceedings of the IEEE International Conference on Communications (ICC '09), pp. 1–5, June 2009. View at: Publisher Site  Google Scholar
 I. Stojmenovic and X. Lin, “Poweraware localized routing in wireless networks,” IEEE Transactions on Parallel and Distributed Systems, vol. 12, no. 11, pp. 1122–1133, 2001. View at: Publisher Site  Google Scholar
 M. T. Isik and O. B. Akan, “A three dimensional localization algorithm for underwater acoustic sensor networks,” IEEE Transactions on Wireless Communications, vol. 8, no. 9, pp. 4457–4463, 2009. View at: Publisher Site  Google Scholar
 S. Basagni, M. Nati, and C. Petrioli, “Localization errorresilient geographic routing for wireless sensor networks,” in Proceedings of the IEEE Global Telecommunications Conference (GLOBECOM '08), pp. 870–875, December 2008. View at: Publisher Site  Google Scholar
 S. P. Tarzia, P. A. Dinda, R. P. Dick, and G. Memik, “Indoor localization without infrastructure using the acoustic background spectrum,” in Proceedings of the 9th International Conference on Mobile Systems, Applications, and Services (MobiSys '11), pp. 155–168, July 2011. View at: Publisher Site  Google Scholar
 E. Martin, O. Vinyals, G. Friedland, and R. Bajcsy, “Precise indoor localization using smart phones,” in Proceedings of the 18th ACM International Conference on Multimedia ACM Multimedia (MM '10), pp. 787–790, October 2010. View at: Publisher Site  Google Scholar
 H.L. Chang, B. Tian Jr., T.T. Lai, H.H. Chu, and P. Huang, “Spinning beacons for precise indoor localization,” in Proceedings of the 6th ACM Conference on Embedded Networked Sensor Systems (SenSys '08), pp. 127–140, ACM, November 2008. View at: Publisher Site  Google Scholar
 P. Barsocchi, S. Lenzi, S. Chessa, and G. Giunta, “Virtual calibration for RSSIbased indoor localization with IEEE 802.15.4,” in Proceedings of the IEEE International Conference on Communications (ICC '09), pp. 1–5, Dresden, Germany, June 2009. View at: Publisher Site  Google Scholar
 H.S. Ahn and W. Yu, “Environmentaladaptive RSSIbased indoor localization,” IEEE Transactions on Automation Science and Engineering, vol. 6, no. 4, pp. 626–633, 2009. View at: Publisher Site  Google Scholar
 H.S. Kim and J.S. Choi, “Advanced indoor localization using ultrasonic sensor and digital compass,” in Proceedings of the International Conference on Control, Automation and Systems (ICCAS '08), pp. 223–226, Seoul, Republic of Korea, October 2008. View at: Publisher Site  Google Scholar
 A. Taok, N. Kandil, and S. Affes, “Neural networks for fingerprintingbased indoor localization using ultrawideband,” Journal of Communications, vol. 4, no. 4, pp. 267–275, 2009. View at: Google Scholar
 L. Ding, X. Gao, W. Wu, W. Lee, X. Zhu, and D. Du, “Distributed construction of connected dominating sets with minimum routing cost in wireless networks,” in Proceedings of the 30th IEEE International Conference on Distributed Computing Systems (ICDCS '10), pp. 448–457, Genoa, Italy, June 2010. View at: Publisher Site  Google Scholar
 Z. Zhong and T. He, “Achieving rangefree localization beyond connectivity,” in Proceedings of the 7th ACM Conference on Embedded Networked Sensor Systems (SenSys '09), pp. 281–294, November 2009. View at: Publisher Site  Google Scholar
 G. Wu, S. Wang, B. Wang, Y. Dong, and S. Yan, “A novel rangefree localization based on regulated neighborhood distance for wireless ad hoc and sensor networks,” Computer Networks, vol. 56, no. 16, pp. 3581–3593, 2012. View at: Publisher Site  Google Scholar
 H. Wang, S. Sen, A. Elgohary, M. Farid, M. Youssef, and R. R. Choudhury, Unsupervised Indoor Localization, Duke University, 2012.
 K. Chintalapudi, A. P. Iyer, and V. N. Padmanabhan, “Indoor localization without the pain,” in Proceedings of the 16th Annual International Conference on Mobile Computing and Networking (MobiCom '10), pp. 173–184, Chicago, Ill, USA, September 2010. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Ping Yi 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.