US20070139269A1  System and method for modelfree position estimation and tracking  Google Patents
System and method for modelfree position estimation and tracking Download PDFInfo
 Publication number
 US20070139269A1 US20070139269A1 US11303140 US30314005A US2007139269A1 US 20070139269 A1 US20070139269 A1 US 20070139269A1 US 11303140 US11303140 US 11303140 US 30314005 A US30314005 A US 30314005A US 2007139269 A1 US2007139269 A1 US 2007139269A1
 Authority
 US
 Grant status
 Application
 Patent type
 Prior art keywords
 position
 recited
 map
 corresponding
 method
 Prior art date
 Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
 Granted
Links
Images
Classifications

 G—PHYSICS
 G01—MEASURING; TESTING
 G01S—RADIO DIRECTIONFINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCEDETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
 G01S5/00—Positionfixing by coordinating two or more direction or position line determinations; Positionfixing by coordinating two or more distance determinations
 G01S5/02—Positionfixing by coordinating two or more direction or position line determinations; Positionfixing by coordinating two or more distance determinations using radio waves
 G01S5/0252—Positionfixing by coordinating two or more direction or position line determinations; Positionfixing by coordinating two or more distance determinations using radio waves by comparing measured values with prestored measured or simulated values
Abstract
Description
 The present invention is directed, in general, to position tracking and, more particularly, to a system and method for position tracking, particularly advantageous with respect to indoor premises, that is not required to depend upon either a propagation model for the region at which the tracking takes place or a motion model for the object being tracked.
 Dynamically estimating the position of a moving person or object, or “tracking,” has been the subject of intense study for decades. Tracking people and objects indoors from signal strength measurements has applications as diverse as security monitoring, selfguided museum tours and personalization of communications services. Accurate dynamic tracking in realtime has been elusive, though, because signal propagation in buildings and the paths that people follow are complex.
 The widespread deployment of wireless networks in buildings provides an opportunity to estimate the position of people and objects in realtime for emergency response, protection of corporate assets, and personalized, positionbased communications. Consequently, many approaches to static position and dynamic tracking from wireless signal strength measurements have been developed, at least in theory. Typically, tracking involves updating static position estimates using a motion model.
 Relating received signal strength (RSS) to position is key to estimating the position of either a static or moving object. Several parametric propagation models that describe path loss across a region have been devised, but none has been unequivocally accepted (see, e.g., Molkdar, “Review on Radio Propagation Into and Within Buildings,” IEEE ProceedingsH, vol. 138, pp. 6173, 1991; and Hashemi, “The Indoor Radio Propagation Channel,” Proceedings of the IEEE, vol. 81, pp. 943968, 1993). Many require restrictive assumptions.
 The most common models are based on inverse exponent laws: 10 log_{10}L(d)=γ log_{10}d+L_{ref}, where d is distance, L(d) is actual path loss and L_{ref }is path loss at a reference distance. Values for γ range from about 1.5 to 5 depending on channel conditions, such as whether the target position is lineofsight down a corridor or whether the transmitter and receiver are on adjacent floors (see, e.g., Rappaport, et al., “UHF Fading in Factories,” IEEE Journal Selected Area of Communications, vol. 7, pp. 4048, 1989; Rappaport, “Indoor Radio Communications for Factories of the Subsequent,” IEEE Communications Magazine, pp. 1524, May 1989; Saleh, et al., “A Statistical Model for Indoor Multipath Propagation,” IEEE Journal Selected Area of Communications, vol. 5, pp. 128137, 1987; Hashemi, supra; and references therein). There can also be additional lognormal shadowing loss with a standard deviation ranging from 3 decibels (dB) to 15 dB that affects RSS readings (id.). All these sources of variation need to be modeled for position estimation and tracking based on a propagation model to be effective.
 More ambitious propagation models do not rely solely on distance but also include the loss incurred from walls, floors, partitions and other obstacles that may lie between the transmitter and the receiver. Path loss can then be determined by adding up attenuation estimates for the various obstacles (see, e.g., Seidel, et al., “900 MHz Path Loss and Prediction Techniques for In Building Communications System Designs,” Proceedings 41st Vehicular Technology Conference, St. Louis Mo., VTC 91, 1991). In practice, however, these extended path loss models offer little improvement over exponent models.
 Including detailed information about the composition of the walls, floors and other structures in the building can help, but this information is difficult to obtain. Moreover, even sophisticated raytracing models that take detailed building maps and materials into account, such as the wellknown WISE software tool (see, e.g., Fortune, et al., “WISE Design of Indoor Wireless Systems: Practical Computation and Optimization,” IEEE Computational Science and Engineering, vol. 2, pp. 5868), have typical errors of 6 dB, which is far too high for accurate tracking. These tools may however be used in planning out the position and number of access points (APs) as well as the distribution of training positions.
 Nonetheless, some position estimation algorithms have been based on path loss models (see, e.g., Chen, et al., “Signal Strength Based Indoor Geoposition,” Proceedings IEEE Conference on Communications, New York, 2002). The parameters of the path loss model may be estimated from a set of training data with known positions using a regression model that includes wall attenuation and a path loss exponent. Position estimates are then based on the fitted model.
 Several approaches to indoor tracking have avoided path loss models and instead drawn an empirical mean RSS map for each of a set of N_{AP }APs by interpolating the mean RSSs obtained at a set of known training positions. RSS measurements at an unknown position are then compared to the RSS maps to estimate the receiver's position. Because there is no path loss model, the positions of the APs need not be known, and indeed are not used in any way.
 Some systems require training measurements to be constantly generated by transceivers, called “emitters,” at fixed positions (see, e.g., Krishnan, et al. “A System for LEASE: System Position Estimation Assisted by Stationary Emitters for Indoor RF Wireless Networks,” In Proceedings of IEEE Infocom, 2004). The emitters allow the current propagation environment, which depends on factors that change over time, such as building occupancy, to be reestimated routinely. A signal strength mean map is estimated by fitting a nonparametric smooth function to the training data using the coordinates of their positions as the covariates in the fitted model. To estimate a position, the J strongest RSS readings may be compared to the corresponding J smoothed RSS maps, ignoring the N_{AP}J weaker signals. The position of the nearest neighbor using Euclidean distance is declared to be the estimated position.
 Perhaps the most wellknown algorithm for tracking is the Kalman filter, which has been used for indoor tracking in, e.g., Guvenc, et al., “Enhancements to RSS Based Indoor Tracking Systems Using Kalman Filters,” International Signal Processing Conference and Global Signal Processing Expo, Dallas, Tex., 2003, and Fox, et al., “Bayesian Techniques for Position Estimation,” IEEE Pervasive Computing, vol. 2, no. 3, 2003, for example. Kalman filters are based on linear motion. Unfortunately, people often do not take linear paths indoors because doors, walls and corners frequently interfere with them. Particle filters have sometimes been used, but their computations are burdensome and not suited to online tracking.
 Voronoi filters that restrict movement to a graph that by design prohibits wall crossing have also been used (see, e.g., Liao, et al., “Voronoi Tracking: Position Estimation Using Sparse and Noisy Sensor Data,” Proceedings of the International Conference on Intelligent Robots and Systems,” IEEE/RSJ, 2003). This procedure can be more accurate, but it also requires significantly more computation.
 Roundtrip delay (see, e.g., Low, et al., “Pulse Detection Algorithm for LineofSight (LOS) UWB Ranging Applications,” IEEE Antennas and Propagation Letters, vol. 4, pp. 6367, 2004) has been used to estimate position and is potentially far more accurate than techniques based on signal strength, but it requires high bandwidth (ultra wideband) signals and specialized hardware. For example, ranging involves detecting the return pulse reflected from a passive device on the target and measuring the round trip delay. For lineofsight positions, this method can achieve centimeter accuracy over distances of about 20 m (id.).
 All of the aforementioned approaches suffer from extensive training, undue computational complexity, inadequate positional accuracy or a combination of these. What is needed in the art is a simple, system and method for performing dynamic tracking that performs well with limited training data. More specifically, what is needed in the art is a system and method that is not required to depend upon either a propagation model or a motion model of the region in which the tracking takes place.
 To address the abovediscussed deficiencies of the prior art, the present invention provides, in one aspect, a system for position estimation of an object relative to a region. In one embodiment, the system includes: (1) a map generator configured to create, without requiring a propagation model, a mean map for the region and (2) a measurement matching engine associated with the map generator and configured to match received measurements regarding the object to a number (e.g., k) of nearest neighboring grid points in the mean map.
 In another aspect, the present invention provides a method of position tracking an object relative to a region. In one embodiment, the method includes: (1) creating, without requiring a propagation model, a mean map for the region and (2) matching received measurements regarding the object to a number (e.g., k) of nearest neighboring grid points in the mean map.
 The foregoing has outlined preferred and alternative features so that those skilled in the pertinent art may better understand the detailed description that follows. Additional features will be described hereinafter that form the subject of the claims. Those skilled in the pertinent art should appreciate that they can readily use the disclosed conception and specific embodiment as a basis for designing or modifying other structures for carrying out the same purposes. Those skilled in the pertinent art should also realize that such equivalent constructions do not depart from the spirit and scope of the present invention.
 For a more complete understanding of the invention, reference is now made to the following descriptions taken in conjunction with the accompanying drawing, in which:

FIG. 1 illustrates a block diagram of one embodiment of a position estimation or tracking system constructed according to the principles of the present invention; 
FIG. 2 illustrates an exemplary graphical representation of the square root of the mean squared error (√{square root over (MSE)}) as a function of distance traveled between two consecutive measurements when α=0.75 and σ=3 dB; 
FIG. 3 illustrates a flow diagram of one embodiment of a method of position estimation or tracking constructed according to the principles of the present invention; 
FIGS. 4A and 4B together illustrate exemplary graphical representations of contours of the estimated mean signal strength and estimated standard deviation of signal strength for one AP for a particular set of sample data; 
FIG. 5 illustrates an exemplary graphical representation of errors in position estimates for an experiment with the mean and standard deviation maps ofFIGS. 4A and 4B ; 
FIG. 6 illustrates an exemplary graphical representation of a simulated distribution of median position error over a region with a training set of 100 RSS measurements from each of 30 positions and each of the 88 positions estimated from one randomly chosen RSS measurement from the position; 
FIG. 7A illustrates an exemplary graphical representation of a true path, position estimates based on the moving average, and smoothed, tracking estimates; and 
FIG. 7B illustrates an exemplary graphical representation of an estimated path using a moving average with a span of eight followed by loess smoothing with a span of 40.  The position tracking system disclosed herein is appropriate for use with any region or premises: outdoor, indoor or a combination of outdoor and indoor. However, as has been described, indoor premises have more complex radio frequency propagation characteristics, making position and tracking more difficult. Therefore, the position tracking system described herein finds particularly advantageous use in indoor premises. The region may be of any size or shape and may have any plural number of associated APs. If the region has more than one floor, the APs may be on any one or combination of floors.
 Various embodiments of a straightforward algorithm contained within a system or carried out in a method for position tracking will be described. The described embodiments require neither a propagation model to determine position nor a motion model to track movement. The algorithm is also relatively simple to compute; many embodiments require only conventional tools for performing Delaunay triangulation, linear interpolation, moving averaging and local regression. Experiments with real data and simulations based on real data will also be set forth below. The experiments will show that the embodiment of the algorithm that is the subject of the experiment is effective.
 Referring initially to
FIG. 1 , illustrated is a block diagram of one embodiment of a position tracking system constructed according to the principles of the present invention. The illustrated embodiment of the indoor position tracking system is generally designated 100 and includes a map generator 110, a mean (and standard deviation) maps database 120, a signal strength measurement matching engine 130 and a path smoother 140.  Before proceeding, it should be noted that the embodiments described in this Detailed Description happen to be based on signal strength. While signal strength is certainly an basis on which position estimates can be made, the present invention is not limited to estimating position or tracking based on signal strength. To the contrary, the present invention fully encompasses round trip delay or any other information a receiver, such as an AP, can provide that can form a basis upon which the position or movement of an object can be determined.
 The map generator 110 is configured to create, without requiring a propagation model, a signal strength mean map for the particular region. In one specific embodiment, a regular grid of points is used to represent the floor of the region. Using the average of the signal strength values obtained through training or in existence as described above, the mean signal strength at each grid point is determined by interpolation for each AP in turn. The result is a signal strength mean map in which a vector is associated with every grid point, each vector having a number of components that equals the number of APs.
 The term “AP” is broadly defined to include any radio transceiver capable of providing signal strength measurements and includes, for example, base stations in cellular telecommunication networks.
 The map generator 110 may also be configured to create, without requiring a propagation model, a signal strength standard deviation map for the region. In the specific embodiment described immediately above, the signal strength standard deviation map is constructed like the signal strength mean map, except that standard deviations are used instead of means.
 The map generator 110 may create the maps through a training process or may simply compile preexisting training data into maps. The mean (and, optionally, standard deviation) maps are stored in the mean (and standard deviation) maps database 120.
 The signal strength measurement matching engine 130 is associated with the map generator 110 and is configured to match signal strength measurements regarding the object to a number (e.g., k) of nearest neighboring grid points in the mean map.
 “Nearest,” as the term is used herein, is based on the distance between (1) a vector composed of signal strength measurements regarding a given object made by multiple APs and (2) vectors corresponding to grid points in the signal strength mean map. Thus, “nearest” is in terms of signal strength, not physical distance. Typically, the distance is a Euclidean distance, but other measures of distance fall within the broad scope of the present invention. Also, those skilled in the art should understand that “nearest” may be expressed in terms of round trip delay, if such is the basis upon which maps are created.
 The means of the knearest neighboring grid points taken from the signal strength mean map are then averaged to yield a static position estimate for the object.
 In one embodiment, the signal strength measurement matching engine 130 is also configured to employ a signal strength standard deviation map to produce the static position estimate. In such case, the average of the means of the knearest neighboring grid points taken from the signal strength mean map is weighted based on standard deviations associated with the knearest neighboring grid points taken from the signal strength standard deviation map to yield the static position estimate for the object. As
FIG. 1 shows, the signal strength measurement matching engine 130 produces at least static position estimates, and perhaps also rough path estimates, assuming that one wants to track the movement of an object and not just determine its static position.  The path smoother 140 is associated with the signal strength measurement matching engine 130 and is configured to make smooth, without requiring a motion model, the static position estimates regarding the object to form smooth path estimates.
 The specific embodiment of the map generator 110 operates as follows. For each of N_{AP }APs, the mean received signal strength (RSS) is mapped across the region by linearly interpolating sample means over a triangular partition of the space, ignoring height. In a more specific embodiment, the map generator 110 further operates to map the standard deviation of RSS by linearly interpolating an estimate of standard deviation over the same triangles. The resulting signal strength mean map and, optionally, the resulting signal strength standard deviation map are stored in the mean (and standard deviation) maps database 120. It should be noted, however, that the present invention can use preexisting training data and is not restricted to triangular partitioning or linear interpolation. Other partitioning and interpolation techniques fall within its broad scope.
 The specific embodiment of the signal strength measurement matching engine 130 operates as follows. Signal strength is estimated at the current position by taking a moving average of the m most recent RSS measurements. Then, the knearest grid point neighbors of this average are found using the N_{AP}dimensional signal strength mean map. Then, a weighted average is taken of their corresponding k positions, the weights being based on the signal strength standard deviation map. With Gaussian weights, this estimate approximates a Bayes estimate of position for a uniform prior density on position, Gaussian RSS measurements and squared error loss. By this process, the illustrated embodiment of the signal strength measurement matching engine 130 produces a stream of static position and rough path estimates. Rough path estimates are paths made up by joining static positions with line segments; they are typically jagged and not smooth.
 The specific embodiment of the path smoother 140 operates as follows. The moving average position estimates are made smooth using locally linear regression. Locally linear regression is known (see, e.g., Cleveland, et al., “Locally Weighted Regression: An approach to Regression Analysis by Local Fitting,” J. Amer. Stat. Assoc., vol. 12, pp. 596610, 1988, incorporated herein by reference). Loess smoothing is one type of locally linear regression and is employed in some embodiments of the present invention. Making the moving average position estimates smooth using locally linear regression increases the stability of the estimated path, and as a result decreases wall crossings. Path smoothing is not present in all embodiments of the present invention and may be performed by a variety of conventional techniques known to those skilled in the pertinent art.
 In the illustrated embodiment, mapping is done offline, and position estimation and tracking are done online whenever RSS measurements for the tracked object become available.
FIG. 2 , described below, particularly illustrates this division.  In certain embodiments described herein, tracking involves only Delaunay triangulation of the region, linear interpolation over triangles, determination of knearest neighboring grid points, weighted averaging, moving averaging and locally linear smoothing. Wellknown, fast algorithms exist for all of these procedures (see, e.g., Fortune, “Voronoi Diagrams and Delaunay Triangulations,” CRC Handbook of Discrete and Computational Geometry, pp. 377388, 1997, for Delaunay triangulation; Clarkson, “Nearest Neighbor Queries in Metric Spaces,” Discrete Comput. Geom., vol. 22, pp. 6393, 1999, and Fortune, et al., “A Note on Rabin's Nearest Neighbor Algorithm,” Information Processing Letters, vol. 8, pp. 2023, 1979 for knearest neighboring grid points; and Cleveland, et al., “LocallyWeighted Regression: An Approach to Regression Analysis by Local Fitting,” Journal of the American Statistical Association, vol. 12, pp. 596610, 1988, and Cleveland, et al., “Computational Methods for Local Fitting,” Statistics and Computing, vol. 12, pp. 4762, 1991, for locally linear smoothing, all of which being incorporated herein by reference).
 Rough guidelines will be provided below for choosing the number of nearest neighboring grid points and the window of the moving average. It will be shown that static position estimation may be carried out to great advantage by taking many neighbors, as many as 100, and weighting them probabilistically.
 The illustrated embodiment of the map generator 110 creates a purely empirical signal strength mean map by partitioning the region (“floor space,” ignoring height) with Delaunay triangulation and then applying trianglebased linear interpolation to the sample mean RSS at a set of training positions.
 The illustrated embodiment of the signal strength measurement matching engine 130 estimates a position by finding its nearest neighboring grid points in the set of empirical maps for substantially all APs, and not just the maps corresponding to the J strongest signals at the position. Further, the average of k nearest neighboring grid points is used to estimate position. Again, k may be many dozens, if not more, closest neighbors. A heuristic for choosing k will be described below.
 The illustrated embodiment of the signal strength measurement matching engine 130 estimates a position by using Gaussian weights with means and standard deviations taken from the maps of the means and standard deviations of signal strengths. This probabilistic weighting gives a Bayes estimate of position for squared error loss when the prior distribution is uniform across the map and RSS measurements have normal distributions.
 Map Generation by One Embodiment of the Map Generator 110.
 The sample mean {circumflex over (μ)}(x,y)=
R (x,y) of RSS measurements from an AP are computed at each position (x,y) in a training set. To take into account changing variability over the region, the mean and standard deviation (MAD) estimate of the standard deviation {circumflex over (σ)}(x,y)=1.4826·medianR−median(R) at each training position is also computed. {circumflex over (σ)} may alternatively be estimated from RSS measurements taken during tracking.  The factor 1.4826 is chosen so that the MAD is an unbiased estimate of the standard deviation for Gaussian measurements. The MAD estimate is used because it is known to be robust against outliers and nonnormality. Empirical mean and standard deviation maps for each AP are obtained by linearly interpolating the sample means and MAD estimates over Delaunay triangles (see, e.g., Fortune, supra).
 Static Position Estimation by One Embodiment of the Signal Strength Measurement Matching Engine 130.
 The number of nearest neighboring grid points, k, is based on an uncertainty in a static position estimate. An unknown position is estimated by finding the knearest neighboring grid points in the empirical mean maps {{circumflex over (μ)}_{i}(x,y)} for i=1, . . . ,N_{AP }to the observed vector of RSS R=(R_{1}, . . . ,R_{N} _{ AP }). If the jnearest neighbor (x_{j},y_{j}) has an estimated mean ({circumflex over (μ)}_{j1}, . . . , {circumflex over (μ)}_{jN} _{ AP }) and an estimated standard deviation ({circumflex over (σ)}_{j1}, . . . , {circumflex over (σ)}_{jN} _{ AP }), it is assigned a weight, e.g.:
${w}_{j}=\prod _{i=1}^{{N}_{\mathrm{AP}}}{\left(\sqrt{2\pi}{\hat{\sigma}}_{\mathrm{ji}}\right)}^{1}\mathrm{exp}\left({\left(\mathrm{Ri}{\hat{\mu}}_{\mathrm{ji}}\right)}^{2}{\left(2{\hat{\sigma}}_{\mathrm{ji}}^{2}\right)}^{1}\right).$  The particular weight given in this equation is the Gaussian probability density function evaluated at R. Given n test measurements at the unknown position, the average measurement is used for finding the knearest neighboring grid points. Each w_{j }is the product of the n Gaussian probability densities evaluated at the n test measurements respectively. The estimated position (x,y) is defined to be the weighted average:
$\begin{array}{cc}\left(\frac{\sum _{j=1}^{k}{w}_{j}{x}_{j}}{\sum _{j=1}^{k}{w}_{j}},\frac{\sum _{j=1}^{k}{w}_{j}{y}_{j}}{\sum _{j=1}^{k}{w}_{j}}\right).& \text{\hspace{1em}}\end{array}$  Conventionally, only a few nearest neighboring grid points, e.g., 1, 2 or 3, and uniform weights have been used, but using a much larger k with Gaussian weights can yield significantly better results. For example, using Gaussian weights over the entire mean map gives an approximate posterior mean for the unknown position when the RSS measurements are normally distributed with a common variance and the prior mean is uniform (see, e.g., Roos, et al., “A Statistical Modeling Approach to Position Estimation,” IEEE Transactions on Mobile Computing, vol. 1, no. 1, JanuaryMarch 2002). The RSS measurements do not have a common variance, so the portion of the map used to estimate the position should be large enough to capture a large fraction of the posterior distribution but not so large as to include many positions with very different distributions of RSS.
 Roughly, if the unknown position is estimated from the sample mean of n test measurements from the unknown position, its estimated standard error is about {circumflex over (σ)}/√{square root over (n)}, where {circumflex over (σ)} is an average, possibly weighted, of the standard deviation estimates across positions. The uncertainty interval should contain at least one standard error centered at the position corresponding to the sample mean at the test position. Thus, k should be approximately the number of grid points that lie in a circle of radius α{circumflex over (σ)}/√{square root over (n)}, where α is the error in distance caused by a 1 dB shift in RSS. When the grid points for a mean map are a distance d apart, it is implied that:
$\begin{array}{cc}k\approx \frac{\pi \text{\hspace{1em}}{a}^{2}{\hat{\sigma}}^{2}}{{\mathrm{nd}}^{2}}& \left(1\right)\end{array}$  In the experiments analyzed in Section 5, which have four APs, α is less than 1. In fact, α=0.75. Furthermore, a threshold, w_{thresh}=(√{square root over (2π)}{circumflex over (σ)})^{−nN} ^{ AP }exp(−0.5nN_{AP}), is used to exclude neighbors with weights less than w.
 Modelfree Online Tracking by One Embodiment of the Signal Strength Measurement Matching Engine 130.
 Linear motion models work well when the trajectory is simple, but not when people stop, turn, reverse direction, and change speeds abruptly, as they do indoors. On the other hand, most paths are at least locally smooth. Thus, some embodiments of the present invention assume with good reason that the path itself is locally linear, with the linear coefficients changing over time. It is important to note that the RSS measurements along a locally linear path need not have a locally linear mean, because indoor propagation is complex.
 Thus, at every RSS reading, the path is updated in two steps. It is assumed that the mean and standard deviation of RSS readings are roughly constant over short paths, so the current position is first estimated from an average, possibly weighted, of RSS measurements over the recent past using the weighted knearest neighbor method described above. It is then assumed that the motion path is locally linear over longer time spans, which allows a smooth, locally linear model to be fitted to a longer segment of the estimated positions in the path obtained so far. Note that the order of operations affects the outcome; the moving average controls for the variability in RSS measurements and the smoothing imposes regularity on the path. These two steps will now be described.
 Constant Position Estimation by One Embodiment of the Signal Strength Measurement Matching Engine 130.
 If the target moves slowly relative to the rate at which the RSS measurements are taken, the mean signal strength on its path should be roughly constant for short time spans. Thus, recent RSS measurements, not just the current RSS measurement alone, provide reasonable estimates of the current position. Since averaging RSS measurements reduces their variance, computing a position estimate for an average RSS, possibly weighted, over a short time span should give more stable position estimates without incurring much bias. Therefore, to estimate the position at a time t at which an RSS measurement is taken:
 (1) the average R of the m most recent RSS measurements is computed, and
 (2) the weighted average of the positions of the knearest neighboring grid points of R in the empirical signal strength mean map is found as described above.
 The position estimate therefore corresponds to a moving average of RSS measurements within a short window or span.
 The bias introduced by using past RSS measurements in the moving average when the target is not stationary may be reduced by delayed tracking. That is, to estimate the position at time t, one may consider not only m past and current measurements, but also n subsequent measurements. The delay in tracking that results from using n subsequent measurements may be acceptable if n is small and RSS readings are obtained frequently.
 The optimal choice of m minimizes the mean squared error (MSE) in estimated position, which depends on the variance σ^{2 }in RSS measurements and the distance traveled over m measurements. If the moving average uses only the current and past measurements (no subsequent measurements) and the distance ν between two adjacent RSS measurements is constant, the average distance of the sites of the measurements from the current position is mν/2. The MSE in the estimate of the current position is then:
MSE=variance+bias^{2}=α^{2}σ^{2}/m+m^{2}ν^{2}/4,
where α is the average change in distance that results from a ldB shift in RSS, as described above. The above equation shows that (1) the variance of the error decreases linearly with 1/m and (2) the bias term increases linearly with m. If σ and ν are roughly constant over the m measurements, the root MSE (RMSE) is minimized at 1.09·(α^{2}σ^{2}ν)^{1/3 }by taking$m\approx {\left(\frac{\sqrt{2}a\text{\hspace{1em}}\sigma}{v}\right)}^{2/3}.$  For example, suppose α=0.75, the RSS measurements are received every 0.3 seconds, σ is about 3 dB along the recent path, and the velocity is about 1 meters per second (m/sec) (i.e., ν=0.3 m).
 Turning now to
FIG. 2 , illustrated is an exemplary graphical representation of RMSE as a function of distance traveled between two consecutive measurements when α=0.75 and σ=3 dB.FIG. 2 shows how the minimal RMSE and the optimal m depend on velocity. The solid line shows a cuberootlaw fit to the data; the numbers beside the solid points show the optimal m for that velocity. If at least one of the m RSS measurements is a subsequent measurement (a measurement taken ahead of the current position), the bias is reduced and the total number of measurements in the moving average can be larger than Equation (2) implies.  Locally Linear Motion
 Computing a position estimate from a moving average of RSS measurements amounts to assuming that the target does not move over short time intervals. These estimates can be used to track a moving target, but they are likely to be noisy because they are based on only a few measurements. (Otherwise, the assumption that the target is not moving might be inadequate.) To make use of older past information, a motion model may be assumed. In the illustrated embodiment, it will be assumed that the path is locally linear, allowing the wellknown loess regression model (see, e.g., Cleveland, et al., “LocallyWeighted Regression . . . ,” supra; and, Cleveland, et al., “Computational Methods . . . ,” supra) to be used, e.g., to make the path smooth.
 As those skilled in the pertinent art understand, the general form of the loess regression model is:
x_{1}=g(t_{i})+ε_{i}, i=1, . . . ,n,
where g is an unknown, locally linear regression curve and the ε_{i }are independent, constant variance Gaussian errors. The assumption of constant variance is approximately true over time spans that are not too large relative to the rate at which the target is moving. The loess model is fit by iteratively reweighted least squares over a span (e.g., 40) of nearest neighboring grid points, where the weight on (t_{i},x_{i}) depends on the distance of the fitted value (t_{i},{circumflex over (x)}_{i}) from the observed (see, Cleveland, et) al., “LocallyWeighted Regression . . . ,” supra, for a detailed explanation of loess smoothing.) Because only a short delay in prediction is desirable in tracking, many more times t_{i }are typically to the left of t than to the right of t in the span of t. A fast algorithm for fitting a loess model can be found in Cleveland, et al., “Computational Methods . . . ,” supra, and is available in R, a wellknown, publicly available, opensource system for statistical computations and graphics (see, e.g., http://cran.rproject.org). Not only is loess smoothing capable of being carried out without an exact motion model, it is also robust to outliers and is not adversely affected by longtailed, nonGaussian distributions.  Finally, isolated measurements may occur at the limit of detection (about −90 dB in the data set forth below). While the previous and subsequent measurements are not close to the limit of detection, the current measurement is. Several explanations may exist for these outliers, such as passing obstacles, but in any case spurious “missing” readings may be corrected by replacing each measurement at the limit of detection by the average of the previous two measurements, and then proceeding as if the replacement were the true reading. Replacement by the mean is not optimal, but is relatively easy to implement. Replacement with a random draw from a normal distribution would probably produce better results.
 Turning now to
FIG. 3 , illustrated is a flow diagram of one embodiment of a method of position estimation or tracking constructed according to the principles of the present invention. The method begins in a start step 310 wherein it is desired to estimate the position of, or track the movement of, objects relative to a particular region.  In a step 320, a training set of RSS measurements is taken from known positions in the region. In a step 330, RSS mean maps for the region are created. In a step 340, RSS standard deviation maps for the region are also created. As
FIG. 3 indicates, the steps 310, 320, 330, 340 are advantageously performed offline and a priori. Creating may simply involve loading a database with preexisting training data pertaining to the premises.  At this point, the method is advantageously performed online or in real time. In a step 350, RSS measurements are matched to RSS mean maps and perhaps RSS standard deviation maps for the knearest neighboring grid points. This yields one or more position estimations. In a step 360, a moving window of static position estimates is locally linearly smoothed to provide tracking. The method ends in an end step 370. A more specific embodiment of a dynamic tracking algorithm will now be described.
 For a nearest neighboring grid points span k, a moving averaging span m, (including n forward RSS measurements), loess spans of s and r past and subsequent estimated positions, the specific embodiment of the dynamic tracking algorithm has the following steps.
 1. For a set of training data with N_{AP }APs,
 (a) For each AP, create an empirical propagation model by triangle based linear interpolation of sample means of RSS.
 (b) For each AP, create an empirical map of the standard deviation of RSS by trianglebased linear interpolation of MAD estimates.
 2. For t=1,2, . . . , given an N_{AP}dimensional vector of RSS measurements R_{t},
 (a) If R_{t }is at the limit of detection, replace R_{t }with the average of the two previous measurements.
 (b) Compute the moving average maR(t) of R_{t−m+1}, . . . , R_{t+n}.
 (c) Estimate the position ({circumflex over (x)}_{t},ŷ_{t}) of maR(t) by finding the Gaussianweighed average of the positions of its knearest neighboring grid points in the empirical signal strength mean map.
 (d) Fit a loess model to ({circumflex over (x)}_{t−s+1}, . . . ,{circumflex over (x)}_{t+r}) to estimate {circumflex over (X)}_{t}.
 (e) Fit a loess model to (ŷ_{t−s+1}, . . . ,ŷ_{t+r}) to estimate Ŷ_{t}.
 (f) Return ({circumflex over (X)}_{t},Ŷ_{t}) as the estimated position at time t.
 The first major step constitutes training and can be computed offline and a priori; the second major step defines tracking and can be computed online and in real time. The final positions {({circumflex over (X)}_{t},Ŷ_{t})} are smooth with respect to t and serve as the estimated moving path. It would be preferable, but more complicated, to make the bivariate (x_{t},y_{t}) smooth in one step, taking into account the correlation in the two dimensions. This may further reduce wallcrossings.
 Turning now to
FIGS. 4A and 4B , together are illustrated exemplary graphical representations of contours of the estimated mean signal strength and estimated standard deviation of signal strength for one AP 410. Solid lines in the background indicate a building map 420; a closed circle indicates the AP 410, which is located in the northwest stairwell asFIGS. 4A and 4B are oriented. The dotted lines show the Delaunay triangulation based on 30 training positions. As a starting point, sample wifi path data set forth in Evennou, et al., “MapAided Indoor Mobile Positioning System Using Particle Filter,” IEEE Wireless Communications and Networking Conference, New Orleans, Louisiana, March 2005 (incorporated herein by reference), will be used.  Evennou, et al., obtained about 100 RSS measurements from four APs at each of 88 positions in a space of about 35 m×40 m that includes corridors and offices around a walled central courtyard 430 with stairwells in the corners. APs, including the AP 410, were located in each of the four corners (including a northeast corner 440 and a southwest corner 450). RSS measurements from 30 fixed positions that are roughly uniformly spaced over the region of interest (ignoring the courtyard 430, for example) will be used for training. It has been found that using more training positions is costly and has little effect on the quality of position estimation.
 Again,
FIGS. 4A and 4B show the estimated signal strength MAD maps based on all the measurements for the 30 training positions for the AP 410 and an interpolation grid size of 0.25 m. (The mean and standard deviation maps for the other APs are similar, with obvious rotations that depend on the position of the AP.) As would be expected, the mean signal strength diminishes from about −50 dB to about −90 dB (apparently, the limit of detection) with distance from the AP. Surprisingly, the mean signal strength falls off with the square root of distance in both the corridors and offices, with an additional constant penalty for positions in offices. This assumption does not hold with respect to all APs “behind the courtyard wall,” however, so it will not be used.  The most prominent features of the standard deviation maps are the hills in each corner. For the AP 410, the standard deviation is highest in the northeast corner 440 and lowest in the southwest corner 450, probably because the signal in the southwest corner 450 is weak and bounded below by −90 dB. Over the entire space, the estimated standard deviation ranges from 1 dB to 5 dB, with a mean and median standard deviation over the region of about 2.8 dB and 3.0 dB respectively.
 First, the mean and standard deviation maps of
FIGS. 4A and 4B are used to estimate all 88 positions using an arbitrary subset of ten RSS measurements from each AP for each position. Thus, the standard error at any point in the empirical signal strength mean map is about 2.8/√{square root over (10)}˜0.95. Because the grid size is 0.25 m, equation (1) gives k˜22. Thus, position is estimated with a Gaussian weighted average of the 25 nearest neighboring grid points.  Turning now to
FIG. 5 , illustrated is an exemplary graphical representation of errors in position estimates for an experiment with the mean and standard deviation maps shown inFIGS. 4A and 4B . The maps were based on 100 measurements at each of the 30 training positions; positions were estimated from ten randomly chosen RSS measurements from each position. The circles at the end of the line show the true position; the other end of the line denotes the estimated position. The errors are generally small, except outside the convex hull of the training positions. Moreover, the errors are generally larger along the northsouth corridors (not referenced) than they are along the eastwest corridors (also not referenced). The errors also tend to be larger in the corners, i.e. , very close to the APs. The 0.1, 0.5 and 0.9 quartiles of the errors are 0.35 m, 0.96 m and 2.86 m respectively. Similar maps are obtained for other trials under the same conditions (ten RSS measurements from each position to be estimated).  Turning now to
FIG. 6 , illustrated is an exemplary graphical representation of a simulated distribution of median position error over the region with a training set of 100 RSS measurements from each of 30 positions and each of the 88 positions estimated from one randomly chosen RSS measurement from the position. The white boxplots are for uniform weights; the gray boxplots are for Gaussian weights with means and standard deviations taken from the empirical mean and standard deviation maps. Each pair of boxplots represents a different number k of nearest neighboring grid points. The effect of using Gaussian weights rather than uniform weights can be seen inFIG. 6 . There, static position estimates for the 88 positions are based on one randomly selected RSS measurement from each position, using the training maps shown inFIGS. 4A and 4B . The optimal number of nearest neighboring grid points k from equation (1) is now 220, or roughly 200. InFIG. 6 , the boxplots are based on computing the median error over the map for 1000 replications of this experiment. InFIG. 6 , the boxes show the 0.25, 0.5 and 0.75 quantiles, the open circles show outliers, and the dashed lines show the tails of the data without the outliers. As would be expected, Gaussian weights are not much better than uniform weights if k is small; but their importance increases with k. There is about a 5% reduction in the average median error for k=200 and Gaussian weights. Note that the difference between using k=1 and k=5 is negligible, with or without weights.  Now, the performance of an exemplary embodiment of the dynamic tracking algorithm will be analyzed for the sample wifi data (RSS measurements and positions) using the training maps of
FIGS. 4A and 4B . However, the issue of whether the training maps are appropriate for the tracking data should first be evaluated.  The tracking data were collected days after the training data, and the signal strength map(s) might have changed in the interim. In fact, comparing the mean observed RSS in the tracking experiments with the mean observed RSS in the training data, shows that adjusting the RSS measurements for APs 1, 2, 3 and 4 by 1 dB, 1.5 dB, 1.5 dB and 1 dB respectively corrects for the average drift in the mean. More generally, drift in calibration could be evaluated using data routinely collected from known positions, as proposed in, e.g., Krishnan, et al., supra.

FIG. 5 shows the results of applying our tracking algorithm to the sample wifi data with k=100 nearest neighboring grid points and a moving average of m=8 past and current measurements. The 0.1, 0.5, and 0.9 error quartiles over the map are 0.3 m, 1.3 m and 2.5 m respectively for the moving average estimates; the values after smoothing are 0.5 m, 1.4 m and 2.6 m. Thus, smoothing in this case does not reduce errors, but the smooth path has fewer wall crossings and unrealistic jumps.  Finally, simulated paths are constructed by fixing the positions along an actual path and simulating RSS measurements by sampling from normal distributions with means and standard deviations determined by the empirical RSS maps. Accordingly, turning now to
FIGS. 7A and 7B , illustrated are exemplary graphical representations of true and estimated paths.FIG. 7A shows a true path (a light gray line 710), position estimates based on the moving average (a dotted black line 720), and smoothed, tracking estimates (a solid black line 730). The filled circles denote the start of the path; the open circles denote the end of the path.FIG. 7B shows an estimated path using a moving average with a span of 8 (a dotted black line 740) followed by loess with a span of 40 (a solid black line 750) and the true path (a light gray line 760) shadowed by the estimated path, which goes into one room 770 off the corridor.  The same training data used in
FIGS. 4A and 4B were used inFIGS. 7A and 7B . Here the 0.1, 0.5 and 0.9 quartiles of position error are 0.27 m, 0.9 m and 1.8 m respectively after the moving average step. After smoothing the path, the 0.1, 0.5 and 0.9 quartile errors are 0.4 m, 0.8 m and 1.8 m, and the paths are again more realistic.  Table 1, below is a table of averages (over 100 replications) of median error using simulated RSS measurements along the true path in
FIG. 7A as a function of the number of nearest neighboring grid points k, the number m of past and current RSS measurements in the moving average, and whether zero or one subsequent RSS measurements were used in the moving average. Table 1 shows the effect of using k nearest neighboring grid points and a moving average over m−1 past RSS measurements, the current measurement, and either no or one subsequent measurement. Table 1 reports median error over the path, averaged over 100 simulated paths. As would be expected, using one subsequent (t=n+1) RSS measurement to estimate position at time n reduces the median error along the path noticeably (by about 10%). Moreover, with one subsequent RSS measurement, the optimal choices appear to be obtained with an average window of m=8 and k=25 nearest neighboring grid points. Equation (1) suggests k=π0.75^{2}{circumflex over (σ)}^{2}/ 8/16)≈28 in this case, which is quite close to simulated optimum for both the nosubsequentmeasurement and onesubsequentmeasurement cases in the table.TABLE 1 Averages of Median Error Using Simulated RSS Measurements k m 1 5 10 25 50 100 200 500 No Subsequent Measurement 1 2.02 1.99 1.97 1.94 1.91 1.87 1.83 1.82 2 1.45 1.43 1.41 1.39 1.36 1.35 1.34 1.33 4 1.06 1.04 1.03 1.01 1.00 0.99 0.99 0.99 6 0.93 0.91 0.90 0.89 0.88 0.88 0.88 0.88 8 0.91 0.89 0.88 0.87 0.87 0.87 0.87 0.86 10 0.95 0.93 0.92 0.92 0.91 0.91 0.91 0.91 12 1.04 1.02 1.01 1.00 1.00 1.00 1.00 1.00 One Subsequent Measurement 1 2.02 1.99 1.97 1.94 1.91 1.87 1.83 1.82 2 1.46 1.43 1.41 1.39 1.36 1.35 1.34 1.33 4 1.05 1.02 1.01 0.99 0.98 0.97 0.97 0.97 6 0.89 0.87 0.86 0.85 0.84 0.83 0.83 0.83 8 0.85 0.82 0.81 0.80 0.80 0.79 0.79 0.79 10 0.86 0.84 0.83 0.82 0.82 0.82 0.82 0.82 12 0.92 0.90 0.89 0.89 0.88 0.88 0.88 0.88  Although the present invention has been described in detail, those skilled in the pertinent art should understand that they can make various changes, substitutions and alterations herein without departing from the spirit and scope of the invention in its broadest form.
Claims (24)
 1. A system for position estimation of an object relative to a region, comprising:a map generator configured to create, without requiring a propagation model, a mean map for said region; anda measurement matching engine associated with said map generator and configured to configured to match received measurements regarding said object to a number of nearest neighboring grid points in said mean map.
 2. The system as recited in
claim 1 wherein said map generator is further configured to create, without requiring a propagation model, a standard deviation map for said region and said measurement matching engine is further configured to employ said standard deviation map to weight an average of said number of nearest neighboring grid points.  3. The system as recited in
claim 1 wherein said mean map is a signal strength mean map.  4. The system as recited in
claim 2 wherein said standard deviation map is a signal strength standard deviation map.  5. The system as recited in
claim 1 wherein said nearest neighboring grid points are determined by generating vectors from moving averages of test measurements and determining distances between said vectors.  6. The system as recited in
claim 1 further comprising a path smoother associated with said measurement matching engine and configured to make smooth static position estimates regarding said object without requiring a motion model.  7. The system as recited in
claim 1 wherein said map generator is configured to create said mean map by interpolating sample means over a triangular partition of a given area, ignoring height.  8. The system as recited in
claim 7 wherein said map generator is configured to create said mean map by linearly interpolating said sample means.  9. The system as recited in
claim 2 wherein said map generator is configured to create said standard deviation map by interpolating an estimate of standard deviation over a partition of a given area, ignoring height.  10. The system as recited in
claim 9 wherein said map generator is configured to create said standard deviation map by linearly interpolating said estimate of said standard deviation.  11. The system as recited in
claim 6 wherein said path smoother is further configured to locally linearly smooth a moving window of said static position estimates.  12. A method of determining a position of an object in a region having a plurality of access points, comprising:measuring a vector of components representing the position, each component being a characteristic of a corresponding one of the access points, the characteristic being either a transmission signal intensity of the corresponding access point at the position or a time delay for signals sent between the corresponding access point and the position; andestimating the position of the object based on the measured vector, a set of reference vectors corresponding to grid points in the region, and error vectors corresponding to the reference vectors; andwherein each component of each reference vector is an estimate of the value of the characteristic for one of the access points at the corresponding grid point; andwherein each error vector provides estimates of errors in the corresponding reference vector.
 13. The method as recited in
claim 12 wherein the estimating includes evaluating a weighted average over a proper subset of the grid points, the proper subset being those grid points whose corresponding reference vectors are close to the measured vector.  14. The method as recited in
claim 12 wherein the weighted average depends on components of the reference vectors and components of the error estimates.  15. The method as recited in
claim 12 further comprising identifying those of the reference vectors that are at a smaller Euclidean distance from the measured vector in a space of said vectors than others of the reference vectors to determine the proper subset of the grid points.  16. The method as recited in
claim 15 wherein the averaging including weighting each grid point by a corresponding weight that depends on the reference vector corresponding to the same grid point.  17. The method as recited in
claim 16 wherein each weight depends on the error vector corresponding to the same grid point.  18. A method of tracking an object in a region having a plurality of access points, comprising:measuring vectors corresponding to a sequence of positions of the object, each component of one of the vectors being a characteristic of a corresponding of one of the access points, the characteristic being either a transmission signal intensity of the corresponding access point at the position corresponding to the one of the vectors or a time delay for transmitted a signal between the corresponding access point and the position corresponding to the one of the vectors;at each position of the sequence, estimating a position of the object from the measured vectors corresponding to a number of nearby ones of the positions; andwherein the estimating steps are performed such that the number of nearby ones of the positions varies with the position being estimated.
 19. The method as recited in
claim 18 wherein the estimating steps are performed such that the number of nearby ones of the positions depends on the speed of the object.  20. The method as recited in
claim 18 wherein the estimating steps are performed such that the number of nearby ones of the positions decreases as the speed of the object increases.  21. The method as recited in
claim 18 wherein the estimating steps are performed such that the number of nearby ones of the positions depends on estimated errors of the measured vectors at the position being estimated.  22. The method as recited in
claim 18 wherein the estimating steps are performed such that the number of nearby ones of the positions increases as the estimated errors of the measured vectors at the positions being estimated increase.  23. The method as recited in
claim 18 further comprising determining a track of the object from the estimated positions by a method, the method including weighting the estimated positions to determine the track.  24. The method as recited in
claim 23 wherein the method includes performing a loess smoothing of the sequence formed by the estimated positions.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

US11303140 US8089407B2 (en)  20051216  20051216  System and method for modelfree position estimation and tracking 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

US11303140 US8089407B2 (en)  20051216  20051216  System and method for modelfree position estimation and tracking 
Publications (2)
Publication Number  Publication Date 

US20070139269A1 true true US20070139269A1 (en)  20070621 
US8089407B2 US8089407B2 (en)  20120103 
Family
ID=38172813
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

US11303140 Active 20281220 US8089407B2 (en)  20051216  20051216  System and method for modelfree position estimation and tracking 
Country Status (1)
Country  Link 

US (1)  US8089407B2 (en) 
Cited By (39)
Publication number  Priority date  Publication date  Assignee  Title 

US20070258420A1 (en) *  20060508  20071108  Farshid AlizadehShabdiz  Estimation of speed of travel using the dynamic signal strength variation of multiple WLAN access points 
US20070259624A1 (en) *  20060508  20071108  Farshid AlizadehShabdiz  Estimation of speed and direction of travel in a WLAN positioning system 
US20080037689A1 (en) *  20060809  20080214  Tolga Kurt  Adaptive Kalman filtering for fast fading removal 
US20080248808A1 (en) *  20070405  20081009  Farshid AlizadehShabdiz  Estimation of position, speed and bearing using time difference of arrival and received signal strength in a wlan positioning system 
US20080248741A1 (en) *  20070405  20081009  Farshid AlizadehShabdiz  Time difference of arrival based estimation of direction of travel in a wlan positioning system 
US20090105939A1 (en) *  20071022  20090423  Toyota Motor Engineering & Manufacturing North America, Inc.  Vehicle navigation system with obstacle avoidance 
US20090109095A1 (en) *  20071025  20090430  Shohei Hido  Location estimation system, method and program 
US20090252138A1 (en) *  20060508  20091008  Skyhook Wireless, Inc.  Estimation of speed and direction of travel in a wlan positioning system using multiple position estimations 
KR100938806B1 (en)  20080430  20100128  대구대학교 산학협력단  Wireless Sensor Node Position Tracking Method Using Probability Filtering of RF Received Signal Strength 
US20100067388A1 (en) *  20061201  20100318  Jun Whan Kim  Method for reserving propagation time by estimating channel environment in wireless personal area network 
US20100127933A1 (en) *  20081127  20100527  Industrial Technology Research Institute  Algorithm of collecting and constructing training location data in a positioning system and the positioning method therefor 
US20100322077A1 (en) *  20061107  20101223  Electronics And Telecommunications Research Institute  Environment analysis system and method for indoor wireless location 
US20110043373A1 (en) *  20090820  20110224  Trimble Navigation Limited  Reverse RFID location system 
WO2011071199A1 (en) *  20091210  20110616  (주)한울옵틱스  System and method for positiontracking 
US20110148711A1 (en) *  20091221  20110623  Industrial Technology Research Institute  Positioning method and communication system using the same 
US20110260922A1 (en) *  20100421  20111027  Motorola, Inc.  Mapping locations based on received signal strengths 
US8049668B1 (en) *  20090312  20111101  Agilent Technologies, Inc.  Method and system for locating signal emitters using iterative residual weighting 
US8077090B1 (en) *  20100615  20111213  Microsoft Corp.  Simultaneous localization and RF modeling 
US8098200B1 (en)  20090312  20120117  Agilent Technologies, Inc.  Method and system for locating signal emitters using residual values 
US20120169530A1 (en) *  20101230  20120705  Honeywell International Inc.  Portable housings for generation of building maps 
US20120225665A1 (en) *  20090619  20120906  Paul Dean Alexander  Characterisation of a wireless communications link 
EP2557433A1 (en) *  20110812  20130213  Sony Mobile Communications AB  Indoor positioning with rake receivers 
EP2637035A1 (en) *  20120309  20130911  BlackBerry Limited  Location estimation using radio scene signatures 
US20130237278A1 (en) *  20120309  20130912  Research In Motion Limited  Signature based connectivity prediction for portable electronic devices 
US8750617B2 (en)  20120309  20140610  Blackberry Limited  Signature representation of data having high dimensionality 
WO2014093000A1 (en) *  20121214  20140619  Apple Inc.  Location fingerprinting 
US8909245B2 (en)  20061107  20141209  Skyhook Wireless, Inc.  System and method for estimating positioning error within a WLANbased positioning system 
US9042907B2 (en)  20120309  20150526  Blackberry Limited  Location identification of a portable electronic device based on characteristics of an operating environment of the portable electronic device 
US9069634B2 (en)  20120309  20150630  Blackberry Limited  Signature representation of data with aliasing across synonyms 
US9103900B2 (en)  20060707  20150811  Skyhook Wireless, Inc.  System and method of gathering WLAN packet samples to improve position estimates of WLAN positioning device 
US9179248B2 (en) *  20120316  20151103  MCube Inc.  Methods and devices for position determination using regional minimal residual estimation 
US9194934B2 (en)  20120309  20151124  Blackberry Limited  Location estimation using radio scene signatures 
WO2016206849A1 (en) *  20150625  20161229  Bayerische Motoren Werke Aktiengesellschaft  Method, computerreadable medium, and system for generating reference data for determining a position of objects relative to a vehicle 
US9588217B2 (en)  20120327  20170307  Microsoft Technology Licensing, Llc  Locating a mobile device 
US9612121B2 (en)  20121206  20170404  Microsoft Technology Licensing, Llc  Locating position within enclosure 
WO2017112118A1 (en) *  20151222  20170629  Mastercard International Incorporated  Systems and methods for use in determining detailed locations for certain entities 
US20170227363A1 (en) *  20150108  20170810  Uti Limited Partnership  Method and apparatus for enhanced pedestrian navigation based on wlan and mems sensors 
US9955358B2 (en)  20060508  20180424  Skyhook Wireless, Inc.  Determining quality metrics utilized in building a reference database 
US9989988B2 (en)  20120203  20180605  Mcube, Inc.  Distributed MEMS devices time synchronization methods and system 
Families Citing this family (7)
Publication number  Priority date  Publication date  Assignee  Title 

US8532962B2 (en) *  20091223  20130910  Honeywell International Inc.  Approach for planning, designing and observing building systems 
US8990049B2 (en)  20100503  20150324  Honeywell International Inc.  Building structure discovery and display from various data artifacts at scene 
US8538687B2 (en)  20100504  20130917  Honeywell International Inc.  System for guidance and navigation in a building 
US9342928B2 (en)  20110629  20160517  Honeywell International Inc.  Systems and methods for presenting building information 
CN103344941B (en) *  20130613  20150812  北京空间飞行器总体设计部  Realtime target detection method based wireless sensor network 
CN103744053B (en) *  20131230  20160203  电子科技大学  One kind of indoor positioning method using a multichannel match vector 
US9216508B2 (en) *  20140114  20151222  Qualcomm Incorporated  Connectivity maintenance using a quality of servicebased robot path planning algorithm 
Citations (7)
Publication number  Priority date  Publication date  Assignee  Title 

US6463290B1 (en) *  19990108  20021008  Trueposition, Inc.  Mobileassisted network based techniques for improving accuracy of wireless location system 
US6556832B1 (en) *  20000204  20030429  Qualcomm Incorporated  Method and apparatus for evaluation of position location performance 
US20040180671A1 (en) *  20030224  20040916  Spain David Stevenson  Location estimation of wireless terminals based on combinations of signal strength measurements and geometryofarrival measurements 
US20050020278A1 (en) *  20030722  20050127  Krumm John C.  Methods for determining the approximate location of a device from ambient signals 
US20050136845A1 (en) *  20030922  20050623  Fujitsu Limited  Method and apparatus for location determination using minibeacons 
US20050257748A1 (en) *  20020802  20051124  Kriesel Marshall S  Apparatus and methods for the volumetric and dimensional measurement of livestock 
US20060025142A1 (en) *  20040802  20060202  Bruce Hamilton  System and method for wireless coverage detection 
Patent Citations (7)
Publication number  Priority date  Publication date  Assignee  Title 

US6463290B1 (en) *  19990108  20021008  Trueposition, Inc.  Mobileassisted network based techniques for improving accuracy of wireless location system 
US6556832B1 (en) *  20000204  20030429  Qualcomm Incorporated  Method and apparatus for evaluation of position location performance 
US20050257748A1 (en) *  20020802  20051124  Kriesel Marshall S  Apparatus and methods for the volumetric and dimensional measurement of livestock 
US20040180671A1 (en) *  20030224  20040916  Spain David Stevenson  Location estimation of wireless terminals based on combinations of signal strength measurements and geometryofarrival measurements 
US20050020278A1 (en) *  20030722  20050127  Krumm John C.  Methods for determining the approximate location of a device from ambient signals 
US20050136845A1 (en) *  20030922  20050623  Fujitsu Limited  Method and apparatus for location determination using minibeacons 
US20060025142A1 (en) *  20040802  20060202  Bruce Hamilton  System and method for wireless coverage detection 
Cited By (60)
Publication number  Priority date  Publication date  Assignee  Title 

US7835754B2 (en) *  20060508  20101116  Skyhook Wireless, Inc.  Estimation of speed and direction of travel in a WLAN positioning system 
US20070259624A1 (en) *  20060508  20071108  Farshid AlizadehShabdiz  Estimation of speed and direction of travel in a WLAN positioning system 
US8103288B2 (en)  20060508  20120124  Skyhook Wireless, Inc.  Estimation of speed and direction of travel in a WLAN positioning system using multiple position estimations 
US8526967B2 (en)  20060508  20130903  Skyhook Wireless, Inc.  Estimation of speed and direction of travel in a WLAN positioning system 
US8090386B2 (en)  20060508  20120103  Skyhook Wireless, Inc.  Estimation of speed and direction of travel in a WLAN positioning system 
US9008690B2 (en)  20060508  20150414  Skyhook Wireless, Inc.  Estimation of speed and direction of travel in a WLAN positioning system 
US8014788B2 (en)  20060508  20110906  Skyhook Wireless, Inc.  Estimation of speed of travel using the dynamic signal strength variation of multiple WLAN access points 
US20090252138A1 (en) *  20060508  20091008  Skyhook Wireless, Inc.  Estimation of speed and direction of travel in a wlan positioning system using multiple position estimations 
US20110058495A1 (en) *  20060508  20110310  Skyhook Wireless, Inc.  Estimation of Speed and Direction of Travel in a WLAN Positioning System 
US9955358B2 (en)  20060508  20180424  Skyhook Wireless, Inc.  Determining quality metrics utilized in building a reference database 
US20070258420A1 (en) *  20060508  20071108  Farshid AlizadehShabdiz  Estimation of speed of travel using the dynamic signal strength variation of multiple WLAN access points 
US9103900B2 (en)  20060707  20150811  Skyhook Wireless, Inc.  System and method of gathering WLAN packet samples to improve position estimates of WLAN positioning device 
US9279877B2 (en)  20060707  20160308  Skyhook Wireless, Inc.  Technique for using cached information with a WLAN positioning system to obtain an estimate of a position of a mobile device 
US20080037689A1 (en) *  20060809  20080214  Tolga Kurt  Adaptive Kalman filtering for fast fading removal 
US9426613B2 (en)  20061107  20160823  Skyhook Wireless, Inc.  System and method for estimating positioning error within a WLANbased positioning system 
US20100322077A1 (en) *  20061107  20101223  Electronics And Telecommunications Research Institute  Environment analysis system and method for indoor wireless location 
US8909245B2 (en)  20061107  20141209  Skyhook Wireless, Inc.  System and method for estimating positioning error within a WLANbased positioning system 
US20100067388A1 (en) *  20061201  20100318  Jun Whan Kim  Method for reserving propagation time by estimating channel environment in wireless personal area network 
US8169975B2 (en) *  20061201  20120501  Electronics And Telecommunications Research Institute  Method for reserving propagation time by estimating channel environment in wireless personal area network 
US20080248741A1 (en) *  20070405  20081009  Farshid AlizadehShabdiz  Time difference of arrival based estimation of direction of travel in a wlan positioning system 
US20080248808A1 (en) *  20070405  20081009  Farshid AlizadehShabdiz  Estimation of position, speed and bearing using time difference of arrival and received signal strength in a wlan positioning system 
US20090105939A1 (en) *  20071022  20090423  Toyota Motor Engineering & Manufacturing North America, Inc.  Vehicle navigation system with obstacle avoidance 
US20090109095A1 (en) *  20071025  20090430  Shohei Hido  Location estimation system, method and program 
US8405551B2 (en)  20071025  20130326  International Business Machines Corporation  Location estimation system, method and program 
US8138974B2 (en) *  20071025  20120320  International Busines Machines Corporation  Location estimation system, method and program 
KR100938806B1 (en)  20080430  20100128  대구대학교 산학협력단  Wireless Sensor Node Position Tracking Method Using Probability Filtering of RF Received Signal Strength 
US20100127933A1 (en) *  20081127  20100527  Industrial Technology Research Institute  Algorithm of collecting and constructing training location data in a positioning system and the positioning method therefor 
US8102315B2 (en) *  20081127  20120124  Industrial Technology Research Institute  Algorithm of collecting and constructing training location data in a positioning system and the positioning method therefor 
US8049668B1 (en) *  20090312  20111101  Agilent Technologies, Inc.  Method and system for locating signal emitters using iterative residual weighting 
US8098200B1 (en)  20090312  20120117  Agilent Technologies, Inc.  Method and system for locating signal emitters using residual values 
US20120225665A1 (en) *  20090619  20120906  Paul Dean Alexander  Characterisation of a wireless communications link 
US9319157B2 (en) *  20090619  20160419  Cohda Wireless Pty. Ltd.  Characterisation of a wireless communications link 
US8797141B2 (en) *  20090820  20140805  Trimble Navigation Limited  Reverse RFID location system 
US20110043373A1 (en) *  20090820  20110224  Trimble Navigation Limited  Reverse RFID location system 
WO2011071199A1 (en) *  20091210  20110616  (주)한울옵틱스  System and method for positiontracking 
US20110148711A1 (en) *  20091221  20110623  Industrial Technology Research Institute  Positioning method and communication system using the same 
US20110260922A1 (en) *  20100421  20111027  Motorola, Inc.  Mapping locations based on received signal strengths 
US8077090B1 (en) *  20100615  20111213  Microsoft Corp.  Simultaneous localization and RF modeling 
US8773946B2 (en) *  20101230  20140708  Honeywell International Inc.  Portable housings for generation of building maps 
US20120169530A1 (en) *  20101230  20120705  Honeywell International Inc.  Portable housings for generation of building maps 
EP2557433A1 (en) *  20110812  20130213  Sony Mobile Communications AB  Indoor positioning with rake receivers 
US8879607B2 (en)  20110812  20141104  Sony Corporation  Indoor positioning with rake receivers 
US9989988B2 (en)  20120203  20180605  Mcube, Inc.  Distributed MEMS devices time synchronization methods and system 
US8948809B2 (en) *  20120309  20150203  Blackberry Limited  Signature based connectivity prediction for portable electronic devices 
US8908978B2 (en)  20120309  20141209  Blackberry Limited  Signature representation of data having high dimensionality 
US8750617B2 (en)  20120309  20140610  Blackberry Limited  Signature representation of data having high dimensionality 
US20130237278A1 (en) *  20120309  20130912  Research In Motion Limited  Signature based connectivity prediction for portable electronic devices 
EP2637035A1 (en) *  20120309  20130911  BlackBerry Limited  Location estimation using radio scene signatures 
US9194934B2 (en)  20120309  20151124  Blackberry Limited  Location estimation using radio scene signatures 
US9042907B2 (en)  20120309  20150526  Blackberry Limited  Location identification of a portable electronic device based on characteristics of an operating environment of the portable electronic device 
US9069634B2 (en)  20120309  20150630  Blackberry Limited  Signature representation of data with aliasing across synonyms 
US9179248B2 (en) *  20120316  20151103  MCube Inc.  Methods and devices for position determination using regional minimal residual estimation 
US9869748B2 (en)  20120327  20180116  Microsoft Technology Licensing, Llc  Locating a mobile device 
US9588217B2 (en)  20120327  20170307  Microsoft Technology Licensing, Llc  Locating a mobile device 
US9612121B2 (en)  20121206  20170404  Microsoft Technology Licensing, Llc  Locating position within enclosure 
US8977298B2 (en)  20121214  20150310  Apple Inc.  Location fingerprinting 
WO2014093000A1 (en) *  20121214  20140619  Apple Inc.  Location fingerprinting 
US20170227363A1 (en) *  20150108  20170810  Uti Limited Partnership  Method and apparatus for enhanced pedestrian navigation based on wlan and mems sensors 
WO2016206849A1 (en) *  20150625  20161229  Bayerische Motoren Werke Aktiengesellschaft  Method, computerreadable medium, and system for generating reference data for determining a position of objects relative to a vehicle 
WO2017112118A1 (en) *  20151222  20170629  Mastercard International Incorporated  Systems and methods for use in determining detailed locations for certain entities 
Also Published As
Publication number  Publication date  Type 

US8089407B2 (en)  20120103  grant 
Similar Documents
Publication  Publication Date  Title 

Whitehouse et al.  A practical evaluation of radio signal strength for rangingbased localization  
Rizos et al.  Indoor positioning techniques based on wireless LAN  
Bekkali et al.  RFID indoor positioning based on probabilistic RFID map and Kalman filtering  
Venkatesh et al.  A linear programming approach to NLOS error mitigation in sensor networks  
US6148211A (en)  Method and system for estimating a subscriber's location in a cluttered area  
Bergamo et al.  Localization in sensor networks with fading and mobility  
Jourdan et al.  Position error bound for UWB localization in dense cluttered environments  
Musicki et al.  Mobile emitter geolocation and tracking using TDOA and FDOA measurements  
US6564065B1 (en)  Bayesianupdate based location prediction method for CDMA systems  
Shen et al.  Performance comparison of TOA and TDOA based location estimation algorithms in LOS environment  
Bose et al.  A practical path loss model for indoor WiFi positioning enhancement  
Rong et al.  Angle of arrival localization for wireless sensor networks  
Yang et al.  Indoor localization using improved rssbased lateration methods  
US20050243936A1 (en)  Method and system for determining user location in a wireless communication network  
US7317419B2 (en)  Selfcalibrated path loss position estimation and zone of impact determination  
Gopakumar et al.  Localization in wireless sensor networks using particle swarm optimization  
US8077090B1 (en)  Simultaneous localization and RF modeling  
Venkatraman et al.  Hybrid TOA/AOA techniques for mobile location in nonlineofsight environments  
US20110125323A1 (en)  Localization by learning of wavesignal distributions  
US6889053B1 (en)  Likelihoodbased geolocation prediction algorithms for CDMA systems using pilot strength measurements  
US6941144B2 (en)  Method and apparatus for detecting excess delay in a communication signal  
US7751829B2 (en)  Method and apparatus for location determination using minibeacons  
US6263208B1 (en)  Geolocation estimation method for CDMA terminals based on pilot strength measurements  
Goldoni et al.  Experimental analysis of RSSIbased indoor localization with IEEE 802.15. 4  
US20040127228A1 (en)  Method for correcting NLOS error in wireless positioning system 
Legal Events
Date  Code  Title  Description 

AS  Assignment 
Owner name: LUCENT TECHNOLOGIES INC., NEW JERSEY Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:CHEN, AIYOU;HARKO, CRISTINA LAURA;LAMBERT, DIANE;AND OTHERS;REEL/FRAME:017750/0605;SIGNING DATES FROM 20051214 TO 20060330 Owner name: LUCENT TECHNOLOGIES INC., NEW JERSEY Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:CHEN, AIYOU;HARKO, CRISTINA LAURA;LAMBERT, DIANE;AND OTHERS;SIGNING DATES FROM 20051214 TO 20060330;REEL/FRAME:017750/0605 

AS  Assignment 
Owner name: LUCENT TECHNOLOGIES INC., NEW JERSEY Free format text: CORRECTIVE ASSIGNMENT TO CORRECT THE CORRECTIVE ASSIGNMENT TO RERECORD ASSIGNMENT PREVIOUSLY RECORDED ON REEL 017750 FRAME 0605;ASSIGNORS:CHEN, AIYOU;HARKO, CRISTINA LAURA;LAMBERT, DIANE;AND OTHERS;REEL/FRAME:017808/0806;SIGNING DATES FROM 20051214 TO 20060330 Owner name: LUCENT TECHNOLOGIES INC., NEW JERSEY Free format text: CORRECTIVE ASSIGNMENT TO CORRECT THE CORRECTIVE ASSIGNMENT TO RERECORD ASSIGNMENT PREVIOUSLY RECORDED ON REEL 017750 FRAME 0605. ASSIGNOR(S) HEREBY CONFIRMS THE CORRECT EXECUTION DATE OF INVENTORS AIYOU CHEN AND PHILIP A. WHITING;ASSIGNORS:CHEN, AIYOU;HARKO, CRISTINA LAURA;LAMBERT, DIANE;AND OTHERS;SIGNING DATES FROM 20051214 TO 20060330;REEL/FRAME:017808/0806 

AS  Assignment 
Owner name: ALCATELLUCENT USA INC., NEW JERSEY Free format text: MERGER;ASSIGNOR:LUCENT TECHNOLOGIES INC.;REEL/FRAME:027047/0930 Effective date: 20081101 

AS  Assignment 
Owner name: ALCATEL LUCENT, FRANCE Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:ALCATELLUCENT USA INC.;REEL/FRAME:027069/0868 Effective date: 20111013 

AS  Assignment 
Owner name: CREDIT SUISSE AG, NEW YORK Free format text: SECURITY AGREEMENT;ASSIGNOR:LUCENT, ALCATEL;REEL/FRAME:029821/0001 Effective date: 20130130 Owner name: CREDIT SUISSE AG, NEW YORK Free format text: SECURITY AGREEMENT;ASSIGNOR:ALCATEL LUCENT;REEL/FRAME:029821/0001 Effective date: 20130130 

AS  Assignment 
Owner name: SOUND VIEW INNOVATIONS, LLC, NEW JERSEY Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:ALCATEL LUCENT;REEL/FRAME:033416/0763 Effective date: 20140630 

AS  Assignment 
Owner name: ALCATEL LUCENT, FRANCE Free format text: RELEASE BY SECURED PARTY;ASSIGNOR:CREDIT SUISSE AG;REEL/FRAME:033868/0001 Effective date: 20140819 

FPAY  Fee payment 
Year of fee payment: 4 