#### 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 range-based or fingerprinting-based, which requires a lot of time and efforts to do the predeployment. In this paper, we present a 3-dimensional on-demand indoor localization system (3D-ODIL), which can be fingerprint-free and deployed rapidly in a multistorey building. The 3D-ODIL 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 range-free 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 location-aware systems. Location-aware 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 open-sourced Android operating system and the continuous cost-reduction 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 location-aware 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 user-experience 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, indoor-localization-based 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: range-based method and range-free method. Range-based 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 range-based method requires expensive hardware and precise statistics, the range-free 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 3-dimensional on-demand indoor localization system (3D-ODIL), which can be fingerprint-free and deployed effectively in a multistorey building. The contributions are as follows.

In this paper, we propose the idea of 3-dimensional on-demand indoor localization system (3D-ODIL), which can be fingerprint-free and deployed rapidly in a multistorey building. The 3D-DDIL 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 range-free 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 high-low 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 on-demand 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 off-the-shelf commercialized products. The system includes three major components: the mobile anchor platform, a mobile localization unit, and a back-end server. The mobile anchor platform can be a portable LTE/GSM mini-radio station or a traditional Wi-Fi router, which act as a radio signal generator. The mobile localization unit is an Android-based smart phone, which is a commercial device and we can use it to receive signals and communicate with the server. The back-end 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

Three-dimensional 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 3-dimensional localization. Unfortunately, there is no mature and general vertical localization which can be rapidly deployed on-demand. 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 3-dimensional system of coordinate to 2-dimensional 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 free-space 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 free-space 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 Wi-Fi 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 1-pair-anchor 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 on-demand, 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 two-dimensional paradigm although our algorithm can be applied into a three-dimensional 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 well-studied 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 high-dimensional 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 noise-free 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 trade-off 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 trade-off 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 two-dimensional 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 software-hardware 3D-ODDIL 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 3D-ODDL 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 Wi-Fi routers as the anchors and Android-based 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 weighted-average 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 range-based localization in [14, 16, 22, 23]. To undermine the effect from physical channel, range-free algorithms based on a relative distance using trilateration and grid-based methods based on landmarks in school and office contexts emerged in RSD [26], RND [27], UnLoc [28], and EZ [29]. Both range-free and grid-based 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 3-dimensional localization, albeit they depict reliable methods for 2-dimensional localization. In this paper, we focus on the 3-dimensional localization, analyze the characteristics of indoor vertical distribution, and propose a distinctive vertical localization and a horizontal localization in a hybrid idea of range-free and grid-based methods.

#### 7. Conclusion

This paper introduces a new mechanism of localization called 3-dimensional on-demand indoor localization system (3D-ODIL), 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 real-environmental experiments. The result proves that the MSD is more reliable than other approaches for vertical localization. More importantly, the low-cost design and on-demand deployment allow the 3D-ODIL for large-scale 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 CNS-1217791.