US20100190511A1 - Method for ambiguity resolution in location determination - Google Patents
Method for ambiguity resolution in location determination Download PDFInfo
- Publication number
- US20100190511A1 US20100190511A1 US12/754,348 US75434810A US2010190511A1 US 20100190511 A1 US20100190511 A1 US 20100190511A1 US 75434810 A US75434810 A US 75434810A US 2010190511 A1 US2010190511 A1 US 2010190511A1
- Authority
- US
- United States
- Prior art keywords
- location determination
- ranging signals
- ranking
- ambiguous
- solution
- 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.)
- Abandoned
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/0257—Hybrid positioning
- G01S5/0268—Hybrid positioning by deriving positions from different combinations of signals or of estimated positions in a single positioning system
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/0252—Radio frequency fingerprinting
- G01S5/02521—Radio frequency fingerprinting using a radio-map
- G01S5/02524—Creating or updating the radio-map
- G01S5/02527—Detecting or resolving anomalies in the radio frequency fingerprints of the radio-map
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/0257—Hybrid positioning
- G01S5/0258—Hybrid positioning by combining or switching between measurements derived from different systems
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01L—SEMICONDUCTOR DEVICES NOT COVERED BY CLASS H10
- H01L2224/00—Indexing scheme for arrangements for connecting or disconnecting semiconductor or solid-state bodies and methods related thereto as covered by H01L24/00
- H01L2224/01—Means for bonding being attached to, or being formed on, the surface to be connected, e.g. chip-to-package, die-attach, "first-level" interconnects; Manufacturing methods related thereto
- H01L2224/42—Wire connectors; Manufacturing methods related thereto
- H01L2224/47—Structure, shape, material or disposition of the wire connectors after the connecting process
- H01L2224/48—Structure, shape, material or disposition of the wire connectors after the connecting process of an individual wire connector
- H01L2224/4805—Shape
- H01L2224/4809—Loop shape
- H01L2224/48091—Arched
Definitions
- the disclosed subject matter relates generally to methods for location determination, and more particularly to methods using various heuristics techniques for resolving ambiguity in location determination in environments with or without noise.
- range-based location determination systems measurements of ranging signals from a plurality of sources are converted to distance information associated with the source of each ranging signal. Distances to different sources with known locations are combined to solve for the unknown user location via geometric techniques known, for example as trilateration (a.k.a. triangulation). If delay of ranging signals cannot be known reliably (e.g. in asynchronous systems where the user clock is not synchronized to the network), location determination algorithms may treat user clock temporal bias as another unknown, to be solved for by the trilateration process, using an additional measurement.
- trilateration a.k.a. triangulation
- a location determination system is ambiguous if more than one location determination solution set of user coordinates and clock temporal bias is consistent with a set of distance measurements.
- Location determination systems can produce ambiguous location determination solutions in three distinct ways: first, ambiguity can be caused by insufficient measurements; second, ambiguity may be introduced by the properties of the algorithm employed in location determination; and third, ambiguity may be introduced by the presence of noisy measurements.
- a system has an insufficient number of measurements when the number of unknowns is greater than or equal to the number of independent measurements.
- the unknowns are the two-dimensional user spatial coordinates and user clock temporal bias.
- FIG. 1 There are three unknowns, namely the mobile station latitude, longitude and clock temporal bias.
- Given three independent distance measurements there are two possible location determination solutions, depicted at the intersection of each set of circles.
- the nature of the algorithm used for locating a user can also be a source of ambiguity.
- a well known algorithm that is susceptible to ambiguity is described in the U.S. Pat. No. 6,289,280.
- This algorithm solves for unknowns using a closed form system of equations. Because it solves for the user location algebraically, this algorithm runs efficiently, making it suitable to applications and devices with time or resource constraints.
- the solution uses linear algebra manipulations to combine the measurements into a system of quadratic equations where the number of equations equals the number of unknowns. Two solutions are produced associated with the two roots of the quadratic equations. The two solutions form an ambiguous set of solutions which needs to be resolved by additional means.
- the unknowns are the two-dimensional user spatial coordinates and user clock temporal bias.
- the system can be said to have a sufficient number of measurements to unambiguously solve for the user location.
- the algebraic method is used, the four measurements are combined into three “average” measurements and two solutions corresponding to these averages are identified, as shown in FIG. 2 .
- noisy measurements can lead to error in the determination of user location.
- a method for location determination in noisy environments by assuming the noise to be a discrete variable with known or computable statistical parameters.
- a set of adjusted measurements and corresponding solution are generated for each assumed noise level.
- Such location determination system is ambiguous, thus also warranting ambiguity resolution techniques.
- FIG. 3 consider three noise levels, each 100 meters apart, associated with the measurement from base station BS 2 .
- a set of circles is plotted as before, with a radius corresponding to the sum of distance measurement (in the case of BS 2 , this measurement is adjusted by the assumed noise level) and the clock temporal bias computed.
- a method for resolving ambiguity in location determination with N ambiguous location determination solutions using clock temporal bias value includes the following steps: generating a probability distribution function (PDF) model for clock temporal bias, obtaining N clock temporal bias values for the N ambiguous location determination solutions, inserting each of the N clock temporal bias values into the PDF model, evaluating the PDF model to get N PDF values, setting N goodness metrics to the N PDF values, comparing the N goodness metrics, defining a maximum goodness metric as the largest of the N goodness metrics, and selecting a final location determination solution having the maximum goodness metric.
- PDF probability distribution function
- a method for resolving ambiguity in location determination with a plurality of ambiguous location determination solutions using the order of receipt (a.k.a. time of arrival) of a plurality of ranging signals includes the following steps: ranking the plurality of ranging signals based on the order of actual time of arrival (TOA) from earliest to latest, ranking the plurality of ranging signals based on an expected order of time of arrival corresponding to each of the plurality of ambiguous location determination solutions, and comparing the ranking of the ranging signals based on the order of the actual time of arrival (TOA) and based on the expected order of time of arrival corresponding to each of the plurality of ambiguous location determination solutions.
- TOA actual time of arrival
- TOA actual time of arrival
- a method for resolving ambiguity in location determination with a plurality of ambiguous location determination solutions using the order of received power levels of a plurality of ranging signals includes the following steps: ranking the plurality of ranging signals based on the order of received power levels from strongest to weakest, ranking the plurality of ranging signals based on an expected order of received power levels corresponding to each of the plurality of ambiguous location determination solutions, and comparing the ranking of the ranging signals based on the order of received power levels and based on the expected order of received power levels corresponding to each of the plurality of ambiguous location determination solutions.
- a method for resolving ambiguity in location determination with a plurality of ambiguous location determination solutions includes the following steps: deriving a plurality of distances D 1 corresponding to each of the plurality of ambiguous location determination solutions for each of M sources, obtaining a plurality of original distances D 2 , comparing each of the plurality of distances D 1 with each of its corresponding plurality of original distances D 2 and computing a plurality of error measurements, and selecting a final location determination solution as having the lowest value of the plurality of error measurements.
- a method for resolving ambiguity in location determination in a noisy environment includes the following steps: selecting L noise levels iteratively within a range [A, B] with an increment I for a first of Q ranging signals, repeating the selecting step Q minus 1 times for each of the rest of the Q ranging signals to produce a plurality of L noise levels, creating a plurality of ambiguous location determination solutions for each of the plurality of L noise levels based on a geometric technique, selecting the final location determination solution from the plurality of ambiguous location determination solutions based on a discriminator function.
- FIG. 1 graphically illustrates an example of three measurements with unknown user spatial coordinates and clock temporal bias value with two ambiguous location determination solutions.
- FIG. 2 graphically illustrates an example of four measurements with three unknowns with two ambiguous location determination solutions.
- FIG. 3 graphically illustrates an example of noisy measurements with three ambiguous location determination solutions.
- FIG. 4 graphically illustrates the relationship between the clock temporal bias value and the horizontal error for each of the ambiguous location determination solutions for the Manhattan dataset.
- FIG. 5 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution #1 and solution #2) and the chosen final location determination solution by resolving ambiguity through clock temporal bias value constraints and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Manhattan dataset.
- CDF cumulative distribution function
- FIG. 6 graphically illustrates the histogram of the clock temporal bias produced by each of the two ambiguous location determination solutions (solution #1 and solution #2) for the Manhattan dataset.
- FIG. 7 graphically illustrates the relationship between the clock temporal bias value and the horizontal error for each of the ambiguous location determination solutions for the Campbell dataset.
- FIG. 8 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution #1 and solution #2) and the chosen final location determination solution by resolving ambiguity through clock temporal bias value constraints and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Campbell dataset.
- CDF cumulative distribution function
- FIG. 9 graphically illustrates the relationship between the clock temporal bias value and the horizontal error for each of the ambiguous location determination solutions for the Japan dataset.
- FIG. 10 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution #1 and solution #2) and the chosen final location determination solution by resolving ambiguity through clock temporal bias value constraints and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Japan dataset.
- CDF cumulative distribution function
- FIG. 11 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution #1 and solution #2) and the chosen final location determination solution by resolving ambiguity through consistency in order of ranging signal arrival and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Manhattan dataset.
- CDF cumulative distribution function
- FIG. 12 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution #1 and solution #2) and the chosen final location determination solution by resolving ambiguity through consistency in order of ranging signal arrival and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Japan dataset.
- CDF cumulative distribution function
- FIG. 13 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution #1 and solution #2) and the chosen final location determination solution by resolving ambiguity through consistency of distance values and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Manhattan dataset.
- CDF cumulative distribution function
- FIG. 14 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution #1 and solution #2) and the chosen final location determination solution by resolving ambiguity through consistency of distance values and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Japan dataset.
- CDF cumulative distribution function
- FIG. 15 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution #1 and solution #2) with the best ambiguous location determination solutions with and without noise removal and are plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Manhattan dataset.
- CDF cumulative distribution function
- Various heuristic approaches for ambiguity resolution are disclosed that exploit known or expected behavior of one or more unknowns or combinations thereof, consistency of various quantifiable parameters, and/or redundancy from measurements or different parameters within a measurement.
- a main advantage of the heuristic approaches disclosed is that redundant measurements are not required for the purpose of location determination solution disambiguation.
- the Manhattan database consists of 1316 fixes from about 25 locations across the city, with about 20 to 100 position fixes per location.
- the Japan dataset consists of about 1912 fixes from all over Japan.
- the Campbell dataset consists of a stationary location with about 2000 fixes.
- the method for ambiguity resolution uses physical constraints on the behavior of some of the unknowns and combinations thereof.
- asynchronous systems can benefit from a set of assumptions or information concerning clock temporal bias value, which can be exploited to disambiguate the user location.
- user clock temporal bias value may include the time the reference timing signal takes to propagate from the source of the reference timing signal to the user, some multipath and receiver processing overhead and clock slewing error.
- PDF probability distribution function
- User clock temporal bias values may be characterized by known statistical distributions, with distribution parameters dictated by various factors, such as the receiver design, the user clock design characteristics, multipath characteristics, and, in a terrestrial system, by network deployment issues such as base station density.
- the distribution parameters may be set statically or learned dynamically for a given location system.
- the goodness metric of a final location determination solution can be derived from the value of the computed clock temporal bias value and the probability distribution function (PDF) of clock temporal bias value for that scenario.
- PDF probability distribution function
- the following steps are needed: First, generate a probability distribution function (PDF) model for the N clock temporal bias values.
- the PDF model is parameterized by its mean and standard deviation.
- Fourth, for each of the ambiguous location determination solutions evaluate the PDF model at the inserted clock temporal bias value to get a PDF value. The result should be N PDF values for the N clock temporal bias values.
- N goodness metrics corresponding to the N PDF values which correspond to the N clock temporal bias values.
- Sixth compare the N goodness metrics and defined a maximum goodness metric as the largest of the N goodness metrics.
- Seventh select the final location determination solution from the ambiguous location determination solution with the maximum goodness metric.
- the quantity N is defined as an integer greater than one.
- clock temporal bias value is modeled to follow a Gaussian probability distribution function (PDF) with a mean and standard deviation of minus 100 meters and 100 meters respectively. If the computed clock temporal bias value (multiplied by the speed of light) is minus 200 meters, then goodness metric can be set to the value of the Gaussian probability distribution function one standard deviation away from the mean.
- PDFs other than a Gaussian PDF (for example, a finite set of statistical parameters) may be used.
- clock temporal bias value is assumed to have a static, simple uniform distribution, with permissible clock temporal bias value between some value A and another value B.
- a and B are set to minus 1000 meter and 1000 meters respectively.
- a and B are set to minus 2000 meters and 100 meters to better suit the Japan dataset as observed empirically.
- FIG. 4 the relationship between clock temporal bias value and the horizontal error of each ambiguous location determination solution is shown in FIG. 4 .
- FIG. 5 For large errors, it is observed that there is a strong, generally linear dependence of error on clock temporal bias value.
- the disambiguation algorithm works well in that the final location determination solution statistically outperforms either of the two ambiguous location determination solutions across the entire dataset. It is noted that the ambiguous location determination solution labeled “Solution 2 ” generally outperforms the ambiguous location determination solution labeled as “Solution 1 ”. Such discrepancy is apparent and explainable by the behavior of clock temporal bias value as produced by each ambiguous location determination solution, and shown in FIG. 6 . One can see how the second ambiguous location determination solution has a “tight” distribution while the first ambiguous location determination solution produced highly-varying results.
- FIG. 7 graphically illustrates the relationship between the clock temporal bias value and the horizontal error for each of the ambiguous location determination solutions for the Campbell dataset.
- FIG. 8 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution #1 and solution #2) and the chosen final location determination solution plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Campbell dataset.
- CDF cumulative distribution function
- FIG. 9 graphically illustrates the relationship between the clock temporal bias value and the horizontal error for each of the ambiguous location determination solutions
- FIG. 10 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution #1 and solution #2) and the chosen final location determination solution plotted as a cumulative distribution function (CDF) of the horizontal errors in meters.
- CDF cumulative distribution function
- ranging signal order There are various ways for deriving ranging signal order, such as the time of arrival (TOA) of ranging signals and/or the power level of ranging signals.
- TOA time of arrival
- the method for ambiguity resolution uses the order of receipt (a.k.a. time of arrival) of ranging signals from a plurality of sources to resolve ambiguity.
- the expected order corresponding to each solution is derived and compared to the actual order in which the ranging signals are received.
- the chosen final location determination solution is the one with the closest match in ordering, or a combination of a number of ambiguous location determination solutions with the closest match in ordering.
- the following steps are needed: First, rank the ranging signals based on the order of actual time of arrival (TOA) from earliest to latest. Second, for each of the ambiguous location determination solution, rank the ranging signals based on the expected order of time of arrival corresponding to each ambiguous location determination solution. Third, compare the ranking of the ranging signals in the order of the actual time of arrival (TOA) from earliest to latest to the ranking of the ranging signals based on the expected order of arrival corresponding to each ambiguous location determination solution. Fourth, select the final location determination solution as the one with the closest match in ordering, or a combination of a number of ambiguous location determination solutions with the closest match in ordering.
- the a-posteriori ranging signal order is compared to the a-priori ranging signal order derived from the time-of-arrival (TOA) information.
- TOA time-of-arrival
- For each ranging signal its rankings in the a-priori and a-posteriori lists are compared and labeled a “1” in case of a match and a “0” in case of a mismatch.
- the solution with the most matches (i.e., the most “1”s) across all ranging signals is picked as the location determination solution. In case of a tie, the average is used.
- the results for this disambiguation technique applied to the Manhattan dataset are shown in FIG. 11 and for the Japan dataset are shown in FIG. 12 .
- the final location determination solution fares similarly or statistically outperforms either of the two ambiguous location determination solutions.
- the method for ambiguity resolution uses the received power level of ranging signals to discriminate among ambiguous location determination solutions.
- To determine the final location determination solution from a set of ambiguous location determination solutions using the received power levels of the ranging signals the following steps are needed: First, rank the ranging signals based on the order of the received power level received, from strongest to weakest. Second, for each of the ambiguous location determination solution, rank the ranging signals based on the expected order of the received power levels corresponding to each ambiguous location determination solution. The expected order of the received power levels is assumed to be the same as the expected order of the time of arrival corresponding to each ambiguous location determination solution.
- the method for ambiguity resolution uses the numerical values of the distances to M sources and combinations thereof to resolve ambiguity.
- the quantity of sources is defined as M whereby M is an integer greater than one.
- the ambiguous location determination solution with the closest match (as defined by the root-mean-square of D 1 minus D 2 or by the mean square of D 1 minus D 2 ), or some combination of a number of values with the closest matches is chosen as the final location determination solution.
- the following steps are needed: First, for each of the M sources, derive the distance D 1 corresponding to each ambiguous location determination solution. Second, obtain the original distance D 2 (which is the original distance adjusted by clock temporal bias value). Third, compare the distance D 1 with the original distance D 2 and compute their error measurement. In one embodiment, the error measurement is the root-mean-square of D 1 ⁇ D 2 . In another embodiment, the error measurement is the root-sum-square of D 1 ⁇ D 2 . And, in yet another embodiment, the error measurement is the root-mean-square of D 1 ⁇ D 2 normalized by the size of the vector D 1 (or D 2 ). Fourth, select the final location determination solution as having the lowest value of error measurement.
- noisy measurements can also lead to errors in location determination.
- a subset of Q signals is selected by some method so that the algorithm may run efficiently.
- the selection method may select the Q signals randomly or via a more systematic method, such as choosing those signals with the lowest power level (that are more likely to incur noise overhead).
- L noise levels are postulated.
- the location determination solutions correspond to each of L times Q combinations of noise levels.
- a method for noise disambiguation involves the following steps. First, for the first of the Q ranging signals, select L noise levels iteratively within a range [A, B] with an increment I.
- One skilled in the art would know the range [A, B] and the increment I (and hence, the quantity of L) based on the chosen location determination system.
- Second repeat the above step Q minus 1 times for each of the rest of the Q minus 1 ranging signals to produce a plurality of L noise levels.
- a discriminator function known to one skilled in the art, is used to select the final location determination solution from the plurality of ambiguous location determination solutions produced for each of the plurality of noise levels.
- the discriminator function can use one or a combination of goodness metrics such as power ranking, distance values, clock temporal bias constraints, order of receipt of the ranging signals. In one embodiment, the discriminator function jointly maximizes the consistency of power ranking and distance values, while meeting the clock temporal bias constraint.
- Other discriminator functions are known to one skilled in the art and may be used without departing from the spirit of the present invention.
- N, A, B are set to 2, 0, and 500 meters respectively.
- the increment I values is investigated in the set of ⁇ 100 meters, 50 meters ⁇ .
- Disambiguation among solutions is done by means of a combination of power ranking, range consistency and clock temporal bias value.
- the power ranking algorithm works as follows: Using the power level of incoming ranging signals, the a-priori order of ranging signals is computed. Then, for each solution corresponding to a noise setting, the a-posteriori order of ranging signals is computed using geographical distances. The a-priori and a-posteriori orders are compared for a given ranging signal.
- the power ranking is computed as the sum of matches across all ranging signals, weighted for each ranging signal in the following way: if the ranging signal is the first ranging signal in the a-priori list (that is, for the base station which is assumed closest), then the weight is 1 ⁇ 2. Otherwise, the weight is 0.5 times (the number of ranging signals- 1 ). Note the weights are normalized, as they add up to 1. Also, the bias in the weighting allows the strongest ranging signal to be more heavily weighted which has been shown to give improved results.
- the root-mean-square differences in the a-priori and a-posteriori ranges are also computed for each solution. Range consistency is calculated as 1—RRMS/max(RRMS).
- the clock temporal bias value enforces a [ ⁇ 1000 meters to 1000 meters] range on the accepted solutions.
- the final location determination solution is the one that meets the clock temporal bias value constraint while maximizing the product of the power ranking and range consistency. The results are shown in FIG. 15 . It is noted that an improvement is observed over the best algebraic method.
Abstract
A method using various heuristics techniques for resolving ambiguity in location determination in environments with or without noise. A final location determination solution may be determined from a set of ambiguous location determination solutions by using clock temporal bias value, by using consistency information of ranging signal order such as the time of arrival and/or the received power level of the ranging signals, by using the distances to the sources, and/or by using other discriminator functions to select the final location determination solution from a plurality of ambiguous location determination solutions. The main advantage of the heuristic approaches is that redundant measurements are not required for location determination solution disambiguation.
Description
- This application is a divisional application of U.S. patent application Ser. No. 11/034,705, filed on Jan. 12, 2005, which claims the benefit of priority from U.S. Provisional Patent Application No. 60/629,011, filed Nov. 17, 2004, each of which is incorporated by reference for all purposes.
- 1. Field
- The disclosed subject matter relates generally to methods for location determination, and more particularly to methods using various heuristics techniques for resolving ambiguity in location determination in environments with or without noise.
- 2. Background
- In range-based location determination systems, measurements of ranging signals from a plurality of sources are converted to distance information associated with the source of each ranging signal. Distances to different sources with known locations are combined to solve for the unknown user location via geometric techniques known, for example as trilateration (a.k.a. triangulation). If delay of ranging signals cannot be known reliably (e.g. in asynchronous systems where the user clock is not synchronized to the network), location determination algorithms may treat user clock temporal bias as another unknown, to be solved for by the trilateration process, using an additional measurement.
- However, a location determination system is ambiguous if more than one location determination solution set of user coordinates and clock temporal bias is consistent with a set of distance measurements. Location determination systems can produce ambiguous location determination solutions in three distinct ways: first, ambiguity can be caused by insufficient measurements; second, ambiguity may be introduced by the properties of the algorithm employed in location determination; and third, ambiguity may be introduced by the presence of noisy measurements.
- First, a system has an insufficient number of measurements when the number of unknowns is greater than or equal to the number of independent measurements. For example, consider the case where the unknowns are the two-dimensional user spatial coordinates and user clock temporal bias. Consider the case depicted in
FIG. 1 . There are three unknowns, namely the mobile station latitude, longitude and clock temporal bias. There are three base stations, namely BS1, BS2 and BS3, and three associated distance measurements. Circles are plotted centered at a particular base station, with radii given by the sum of the distance between the mobile station and the base station as measured at the mobile station, and the computed clock temporal bias corresponding to a fitting solution. Given three independent distance measurements, there are two possible location determination solutions, depicted at the intersection of each set of circles. - Second, the nature of the algorithm used for locating a user can also be a source of ambiguity. A well known algorithm that is susceptible to ambiguity is described in the U.S. Pat. No. 6,289,280. This algorithm solves for unknowns using a closed form system of equations. Because it solves for the user location algebraically, this algorithm runs efficiently, making it suitable to applications and devices with time or resource constraints. The solution uses linear algebra manipulations to combine the measurements into a system of quadratic equations where the number of equations equals the number of unknowns. Two solutions are produced associated with the two roots of the quadratic equations. The two solutions form an ambiguous set of solutions which needs to be resolved by additional means.
- For example, consider the case where the unknowns are the two-dimensional user spatial coordinates and user clock temporal bias. With four measurements, the system can be said to have a sufficient number of measurements to unambiguously solve for the user location. Yet, when the algebraic method is used, the four measurements are combined into three “average” measurements and two solutions corresponding to these averages are identified, as shown in
FIG. 2 . - Third, noisy measurements can lead to error in the determination of user location. Consider a method for location determination in noisy environments by assuming the noise to be a discrete variable with known or computable statistical parameters. A set of adjusted measurements and corresponding solution are generated for each assumed noise level. Such location determination system is ambiguous, thus also warranting ambiguity resolution techniques. For example, in
FIG. 3 , consider three noise levels, each 100 meters apart, associated with the measurement from base station BS2. For each noise level, a set of circles is plotted as before, with a radius corresponding to the sum of distance measurement (in the case of BS2, this measurement is adjusted by the assumed noise level) and the clock temporal bias computed. There are three ambiguous solutions, associated with each noise level, shown by the intersection of circles. - Accordingly, since more than one possible solution is presented by these prior art algorithms, it would be desirable to provide a method for selecting the correct (a.k.a. final) location determination solution from a set of ambiguous location determination solutions.
- Disclosed are methods for selecting the correct (a.k.a. final) location determination solution from a set of ambiguous location determination solutions using various heuristics and/or noise removal.
- According to one aspect, a method for resolving ambiguity in location determination with N ambiguous location determination solutions using clock temporal bias value includes the following steps: generating a probability distribution function (PDF) model for clock temporal bias, obtaining N clock temporal bias values for the N ambiguous location determination solutions, inserting each of the N clock temporal bias values into the PDF model, evaluating the PDF model to get N PDF values, setting N goodness metrics to the N PDF values, comparing the N goodness metrics, defining a maximum goodness metric as the largest of the N goodness metrics, and selecting a final location determination solution having the maximum goodness metric.
- According to another aspect, a method for resolving ambiguity in location determination with a plurality of ambiguous location determination solutions using the order of receipt (a.k.a. time of arrival) of a plurality of ranging signals includes the following steps: ranking the plurality of ranging signals based on the order of actual time of arrival (TOA) from earliest to latest, ranking the plurality of ranging signals based on an expected order of time of arrival corresponding to each of the plurality of ambiguous location determination solutions, and comparing the ranking of the ranging signals based on the order of the actual time of arrival (TOA) and based on the expected order of time of arrival corresponding to each of the plurality of ambiguous location determination solutions.
- According to yet another one aspect, a method for resolving ambiguity in location determination with a plurality of ambiguous location determination solutions using the order of received power levels of a plurality of ranging signals includes the following steps: ranking the plurality of ranging signals based on the order of received power levels from strongest to weakest, ranking the plurality of ranging signals based on an expected order of received power levels corresponding to each of the plurality of ambiguous location determination solutions, and comparing the ranking of the ranging signals based on the order of received power levels and based on the expected order of received power levels corresponding to each of the plurality of ambiguous location determination solutions.
- According to still another aspect, a method for resolving ambiguity in location determination with a plurality of ambiguous location determination solutions includes the following steps: deriving a plurality of distances D1 corresponding to each of the plurality of ambiguous location determination solutions for each of M sources, obtaining a plurality of original distances D2, comparing each of the plurality of distances D1 with each of its corresponding plurality of original distances D2 and computing a plurality of error measurements, and selecting a final location determination solution as having the lowest value of the plurality of error measurements.
- According to yet another aspect, a method for resolving ambiguity in location determination in a noisy environment includes the following steps: selecting L noise levels iteratively within a range [A, B] with an increment I for a first of Q ranging signals, repeating the selecting step Q minus 1 times for each of the rest of the Q ranging signals to produce a plurality of L noise levels, creating a plurality of ambiguous location determination solutions for each of the plurality of L noise levels based on a geometric technique, selecting the final location determination solution from the plurality of ambiguous location determination solutions based on a discriminator function.
- It is understood that other embodiments will become readily apparent to those skilled in the art from the following detailed description, wherein it is shown and described various embodiments by way of illustration. The drawings and detailed description are to be regarded as illustrative in nature and not as restrictive.
-
FIG. 1 graphically illustrates an example of three measurements with unknown user spatial coordinates and clock temporal bias value with two ambiguous location determination solutions. -
FIG. 2 graphically illustrates an example of four measurements with three unknowns with two ambiguous location determination solutions. -
FIG. 3 graphically illustrates an example of noisy measurements with three ambiguous location determination solutions. -
FIG. 4 graphically illustrates the relationship between the clock temporal bias value and the horizontal error for each of the ambiguous location determination solutions for the Manhattan dataset. -
FIG. 5 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution #2) and the chosen final location determination solution by resolving ambiguity through clock temporal bias value constraints and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Manhattan dataset. -
FIG. 6 graphically illustrates the histogram of the clock temporal bias produced by each of the two ambiguous location determination solutions (solution # 1 and solution #2) for the Manhattan dataset. -
FIG. 7 graphically illustrates the relationship between the clock temporal bias value and the horizontal error for each of the ambiguous location determination solutions for the Campbell dataset. -
FIG. 8 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution #2) and the chosen final location determination solution by resolving ambiguity through clock temporal bias value constraints and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Campbell dataset. -
FIG. 9 graphically illustrates the relationship between the clock temporal bias value and the horizontal error for each of the ambiguous location determination solutions for the Japan dataset. -
FIG. 10 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution #2) and the chosen final location determination solution by resolving ambiguity through clock temporal bias value constraints and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Japan dataset. -
FIG. 11 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution #2) and the chosen final location determination solution by resolving ambiguity through consistency in order of ranging signal arrival and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Manhattan dataset. -
FIG. 12 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution #2) and the chosen final location determination solution by resolving ambiguity through consistency in order of ranging signal arrival and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Japan dataset. -
FIG. 13 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution #2) and the chosen final location determination solution by resolving ambiguity through consistency of distance values and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Manhattan dataset. -
FIG. 14 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution #2) and the chosen final location determination solution by resolving ambiguity through consistency of distance values and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Japan dataset. -
FIG. 15 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution #2) with the best ambiguous location determination solutions with and without noise removal and are plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Manhattan dataset. - The description set forth below in connection with the appended drawings is intended as a description of various embodiments of the invention and is not intended to represent the only embodiments in which the invention may be practiced. Each embodiment is provided merely as an example or illustration, and should not necessarily be construed as preferred or advantageous over other embodiments. Specific details are used to provide an understanding of the invention. However, it will be apparent to those skilled in the art that the invention may be practiced without these specific details. Acronyms and other descriptive terminology may be used merely for convenience and clarity and are not intended to limit the scope of the invention.
- Various heuristic approaches for ambiguity resolution are disclosed that exploit known or expected behavior of one or more unknowns or combinations thereof, consistency of various quantifiable parameters, and/or redundancy from measurements or different parameters within a measurement. A main advantage of the heuristic approaches disclosed is that redundant measurements are not required for the purpose of location determination solution disambiguation.
- For the purpose of illustration only, simulations using real user data from Manhattan, Japan, and Campbell are presented. The Manhattan database consists of 1316 fixes from about 25 locations across the city, with about 20 to 100 position fixes per location. The Japan dataset consists of about 1912 fixes from all over Japan. The Campbell dataset consists of a stationary location with about 2000 fixes.
- In one embodiment, the method for ambiguity resolution uses physical constraints on the behavior of some of the unknowns and combinations thereof. In particular, asynchronous systems can benefit from a set of assumptions or information concerning clock temporal bias value, which can be exploited to disambiguate the user location. In a system where the user clock timing is learned from some reference timing signal, user clock temporal bias value may include the time the reference timing signal takes to propagate from the source of the reference timing signal to the user, some multipath and receiver processing overhead and clock slewing error. A probability distribution function (PDF) of clock temporal bias value can be derived for that network and used to disambiguate the user location.
- User clock temporal bias values may be characterized by known statistical distributions, with distribution parameters dictated by various factors, such as the receiver design, the user clock design characteristics, multipath characteristics, and, in a terrestrial system, by network deployment issues such as base station density. The distribution parameters may be set statically or learned dynamically for a given location system. The goodness metric of a final location determination solution can be derived from the value of the computed clock temporal bias value and the probability distribution function (PDF) of clock temporal bias value for that scenario.
- To determine the final location determination solution from a set of N ambiguous location determination solutions using N clock temporal bias values, the following steps are needed: First, generate a probability distribution function (PDF) model for the N clock temporal bias values. In one embodiment, the PDF model is parameterized by its mean and standard deviation. Second, for each of the N ambiguous location determination solutions, obtain a clock temporal bias value, resulting in N clock temporal bias value. Third, for each of the ambiguous location determination solutions, insert its corresponding clock temporal bias value into the PDF model. Fourth, for each of the ambiguous location determination solutions, evaluate the PDF model at the inserted clock temporal bias value to get a PDF value. The result should be N PDF values for the N clock temporal bias values. Fifth, for each of the ambiguous location determination solutions, set the goodness metric of the ambiguous location determination solution to the PDF value. Again, there should be N goodness metrics corresponding to the N PDF values which correspond to the N clock temporal bias values. Sixth, compare the N goodness metrics and defined a maximum goodness metric as the largest of the N goodness metrics. Seventh, select the final location determination solution from the ambiguous location determination solution with the maximum goodness metric. The quantity N is defined as an integer greater than one.
- In one embodiment, clock temporal bias value is modeled to follow a Gaussian probability distribution function (PDF) with a mean and standard deviation of minus 100 meters and 100 meters respectively. If the computed clock temporal bias value (multiplied by the speed of light) is minus 200 meters, then goodness metric can be set to the value of the Gaussian probability distribution function one standard deviation away from the mean. One skilled in the art would know that other PDFs other than a Gaussian PDF (for example, a finite set of statistical parameters) may be used.
- In an example, clock temporal bias value is assumed to have a static, simple uniform distribution, with permissible clock temporal bias value between some value A and another value B. For Manhattan and Campbell datasets, A and B are set to minus 1000 meter and 1000 meters respectively. In the case of Japan, A and B are set to minus 2000 meters and 100 meters to better suit the Japan dataset as observed empirically. For the Manhattan dataset, the relationship between clock temporal bias value and the horizontal error of each ambiguous location determination solution is shown in
FIG. 4 . For large errors, it is observed that there is a strong, generally linear dependence of error on clock temporal bias value. For the Manhattan dataset, the performance results are shown inFIG. 5 . It is observed that the disambiguation algorithm works well in that the final location determination solution statistically outperforms either of the two ambiguous location determination solutions across the entire dataset. It is noted that the ambiguous location determination solution labeled “Solution 2” generally outperforms the ambiguous location determination solution labeled as “Solution 1”. Such discrepancy is apparent and explainable by the behavior of clock temporal bias value as produced by each ambiguous location determination solution, and shown inFIG. 6 . One can see how the second ambiguous location determination solution has a “tight” distribution while the first ambiguous location determination solution produced highly-varying results. -
FIG. 7 graphically illustrates the relationship between the clock temporal bias value and the horizontal error for each of the ambiguous location determination solutions for the Campbell dataset.FIG. 8 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution #2) and the chosen final location determination solution plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Campbell dataset. For the Japan dataset,FIG. 9 graphically illustrates the relationship between the clock temporal bias value and the horizontal error for each of the ambiguous location determination solutions whileFIG. 10 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution #2) and the chosen final location determination solution plotted as a cumulative distribution function (CDF) of the horizontal errors in meters. - Alternative methods for ambiguity resolution use consistency of ranging signal order. There are various ways for deriving ranging signal order, such as the time of arrival (TOA) of ranging signals and/or the power level of ranging signals.
- In one embodiment, the method for ambiguity resolution uses the order of receipt (a.k.a. time of arrival) of ranging signals from a plurality of sources to resolve ambiguity. The expected order corresponding to each solution is derived and compared to the actual order in which the ranging signals are received. The chosen final location determination solution is the one with the closest match in ordering, or a combination of a number of ambiguous location determination solutions with the closest match in ordering.
- To determine the final location determination solution from a set of ambiguous location determination solutions using the order of receipt (a.k.a. time of arrival) of ranging signals, the following steps are needed: First, rank the ranging signals based on the order of actual time of arrival (TOA) from earliest to latest. Second, for each of the ambiguous location determination solution, rank the ranging signals based on the expected order of time of arrival corresponding to each ambiguous location determination solution. Third, compare the ranking of the ranging signals in the order of the actual time of arrival (TOA) from earliest to latest to the ranking of the ranging signals based on the expected order of arrival corresponding to each ambiguous location determination solution. Fourth, select the final location determination solution as the one with the closest match in ordering, or a combination of a number of ambiguous location determination solutions with the closest match in ordering.
- In an example, the a-posteriori ranging signal order is compared to the a-priori ranging signal order derived from the time-of-arrival (TOA) information. For each ranging signal, its rankings in the a-priori and a-posteriori lists are compared and labeled a “1” in case of a match and a “0” in case of a mismatch. The solution with the most matches (i.e., the most “1”s) across all ranging signals is picked as the location determination solution. In case of a tie, the average is used. The results for this disambiguation technique applied to the Manhattan dataset are shown in
FIG. 11 and for the Japan dataset are shown inFIG. 12 . In general, the final location determination solution fares similarly or statistically outperforms either of the two ambiguous location determination solutions. - In another embodiment, the method for ambiguity resolution uses the received power level of ranging signals to discriminate among ambiguous location determination solutions. One can use received power level to discriminate among location determination solutions. To determine the final location determination solution from a set of ambiguous location determination solutions using the received power levels of the ranging signals, the following steps are needed: First, rank the ranging signals based on the order of the received power level received, from strongest to weakest. Second, for each of the ambiguous location determination solution, rank the ranging signals based on the expected order of the received power levels corresponding to each ambiguous location determination solution. The expected order of the received power levels is assumed to be the same as the expected order of the time of arrival corresponding to each ambiguous location determination solution. Third, compare the ranking of the ranging signals in the order of the received power level received from strongest to weakest to the ranking of the ranging signals based on the expected order of the received power levels corresponding to each ambiguous location determination solution. Fourth, select the final location determination solution as the one with the closest match in ordering, or a combination of a number of ambiguous location determination solutions with the closest match in ordering.
- In another embodiment, the method for ambiguity resolution uses the numerical values of the distances to M sources and combinations thereof to resolve ambiguity. One can derive each of the distance D1 to each of the M sources corresponding to each ambiguous location determination solution and compare them to the original distance adjusted by the clock temporal bias value (labeled as original distance D2). The quantity of sources is defined as M whereby M is an integer greater than one. The ambiguous location determination solution with the closest match (as defined by the root-mean-square of D1 minus D2 or by the mean square of D1 minus D2), or some combination of a number of values with the closest matches is chosen as the final location determination solution. It will be understood by one skilled in the art that the present invention is not confined to the root-mean-square or mean-square of D1 minus D2 and that other error measures such as, but not limited to, root-sum-square (RSS) or sum-square (SS) may also be used without violating the spirit of the invention. The results for the Manhattan and Japan dataset are shown in
FIG. 13 andFIG. 14 respectively. In both cases, the final location determination solution statistically outperforms either of the two ambiguous location determination solutions across the entire dataset. - To determine the final location determination solution from a set of ambiguous location determination solutions using numerical values of the distances to the M sources, the following steps are needed: First, for each of the M sources, derive the distance D1 corresponding to each ambiguous location determination solution. Second, obtain the original distance D2 (which is the original distance adjusted by clock temporal bias value). Third, compare the distance D1 with the original distance D2 and compute their error measurement. In one embodiment, the error measurement is the root-mean-square of D1−D2. In another embodiment, the error measurement is the root-sum-square of D1−D2. And, in yet another embodiment, the error measurement is the root-mean-square of D1−D2 normalized by the size of the vector D1 (or D2). Fourth, select the final location determination solution as having the lowest value of error measurement.
- Noisy measurements can also lead to errors in location determination. Out of a set of P signals, a subset of Q signals is selected by some method so that the algorithm may run efficiently. The selection method may select the Q signals randomly or via a more systematic method, such as choosing those signals with the lowest power level (that are more likely to incur noise overhead). For each of the Q signals selected, L noise levels are postulated. The location determination solutions correspond to each of L times Q combinations of noise levels. A method for noise disambiguation involves the following steps. First, for the first of the Q ranging signals, select L noise levels iteratively within a range [A, B] with an increment I. One skilled in the art would know the range [A, B] and the increment I (and hence, the quantity of L) based on the chosen location determination system. Second, repeat the above step Q minus 1 times for each of the rest of the Q minus 1 ranging signals to produce a plurality of L noise levels. Third, create a plurality of ambiguous location determination solutions for each of the plurality of noise levels based on known geometric techniques such as trilateration. Fourth, a discriminator function, known to one skilled in the art, is used to select the final location determination solution from the plurality of ambiguous location determination solutions produced for each of the plurality of noise levels. The discriminator function can use one or a combination of goodness metrics such as power ranking, distance values, clock temporal bias constraints, order of receipt of the ranging signals. In one embodiment, the discriminator function jointly maximizes the consistency of power ranking and distance values, while meeting the clock temporal bias constraint. Other discriminator functions are known to one skilled in the art and may be used without departing from the spirit of the present invention.
- In an example, the noise removal techniques for fixes in the Manhattan dataset. N, A, B are set to 2, 0, and 500 meters respectively. The increment I values is investigated in the set of {100 meters, 50 meters}. Disambiguation among solutions is done by means of a combination of power ranking, range consistency and clock temporal bias value. The power ranking algorithm works as follows: Using the power level of incoming ranging signals, the a-priori order of ranging signals is computed. Then, for each solution corresponding to a noise setting, the a-posteriori order of ranging signals is computed using geographical distances. The a-priori and a-posteriori orders are compared for a given ranging signal. Its ranking in the a-priori and a-posteriori lists are compared and a match is labeled a “1” while a mismatch returns a “0.” The power ranking is computed as the sum of matches across all ranging signals, weighted for each ranging signal in the following way: if the ranging signal is the first ranging signal in the a-priori list (that is, for the base station which is assumed closest), then the weight is ½. Otherwise, the weight is 0.5 times (the number of ranging signals-1). Note the weights are normalized, as they add up to 1. Also, the bias in the weighting allows the strongest ranging signal to be more heavily weighted which has been shown to give improved results. The root-mean-square differences in the a-priori and a-posteriori ranges (RRMS) are also computed for each solution. Range consistency is calculated as 1—RRMS/max(RRMS). The clock temporal bias value enforces a [−1000 meters to 1000 meters] range on the accepted solutions. The final location determination solution is the one that meets the clock temporal bias value constraint while maximizing the product of the power ranking and range consistency. The results are shown in
FIG. 15 . It is noted that an improvement is observed over the best algebraic method. - The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention.
Claims (17)
1. A method for resolving ambiguity in location determination with a plurality of ambiguous location determination solutions for a mobile station using time of arrival of a plurality of ranging signals from base stations in a cellular communication network, the method comprising:
receiving the plurality of ranging signals at the mobile station;
forming a first ranking of the plurality of ranging signals based on an order of actual time of arrival (TOA) of the ranging signals from earliest to latest;
forming second rankings of the plurality of ranging signals based on an expected order of time of arrival corresponding to each of the plurality of ambiguous location determination solutions;
comparing the first and second rankings of the ranging signals; and
determining a final location determination solution for the mobile station from the plurality of ambiguous location determination solutions based on a result of comparing the first and second rankings of the ranging signals.
2. The method of claim 1 wherein determining the final location determination solution comprises determining the ambiguous location determination solution for which the second ranking is a closest match to the first ranking.
3. The method of claim 1 wherein determining the final location determination solution comprises determining a combination of the plurality of ambiguous location determination solutions with a closest match in ordering.
4. The method of claim 1 further comprising assigning a score to each of the ambiguous location determination solutions based on comparing the first and second rankings of the ranging signals, wherein determining the final location determination solution comprises determining the ambiguous location determination solution having a highest score.
5. The method of claim 4 wherein assigning the score comprises determining a number of matches between the ranging signals of the first ranking and the ranging signals of each of the second rankings.
6. The method of claim 5 further comprising:
weighting each of the matches according to an estimated distance of the base station from which the corresponding ranging signal is received.
7. A method of determining a final location determination solution for a mobile station from a plurality of candidate location determination solutions, the method comprising:
receiving a plurality of ranging signals at the mobile station;
determining an actual time of arrival for each ranging signal;
forming a first ranking of the ranging signals based on an order of the actual time of arrival of the ranging signals;
forming second rankings of the ranging signals based on an expected order of time of arrival corresponding to each of the plurality of candidate location determination solutions;
assigning a score to each of the second rankings; and
determining the final location determination solution for the mobile station based on the scores of the second rankings.
8. The method of claim 7 wherein determining the final location determination solution comprises determining the candidate location determination solution for which the corresponding second ranking has a highest score.
9. The method of claim 7 wherein assigning a score to each of the second rankings comprises determining a number of matches between the first ranking of the ranging signals and each of the second rankings.
10. The method of claim 7 wherein forming the first ranking comprises ranking the ranging signals from earliest to latest based on the actual time of arrival of each of the ranging signals.
11. A computer-readable medium having stored thereon a processor-readable program of instructions configured to cause a processor to resolve ambiguity in location determination with a plurality of ambiguous location determination solutions for a mobile station using a time of arrival of each of a plurality of ranging signals, the program of instructions being configured to cause the processor to:
form a first ranking of the plurality of ranging signals based on an order of actual time of arrival (TOA) from earliest to latest;
form second rankings of the plurality of ranging signals based on an expected order of time of arrival corresponding to each of the plurality of ambiguous location determination solutions;
compare the first and second rankings of the ranging signals; and
determine a final location determination solution for the mobile station from the plurality of ambiguous location determination solutions based on a result of comparing the first and second rankings of the ranging signals.
12. The computer readable medium of claim 11 wherein the instructions to configured to cause the processor determine the final location determination solution comprise instructions configured to cause the processor to determine an ambiguous location determination solution for which the second ranking is a closest match to the first ranking.
13. The computer readable medium of claim 11 wherein the instructions configured to cause the processor to determine the final location determination solution comprise instructions configured to cause the processor to determine a combination of the plurality of ambiguous location determination solutions with a closest match in ordering.
14. The computer readable medium of claim 11 wherein the program further comprises instructions configured to cause the processor to assign a score to each of the ambiguous location determination solutions based on comparing the first and second rankings of the ranging signals, and wherein the instructions configured to cause the processor to determine the final location determination solution comprise instructions configured to cause the processor to determine the ambiguous location determination solution having a highest score.
15. The computer readable medium of claim 14 wherein the instructions configured to cause the processor to assign the score comprise instructions configured to cause the processor to determine a number of matches between the ranging signals of the first ranking and the ranging signals of each of the second rankings.
16. The computer readable medium of claim 15 wherein the instructions configured to cause the processor to determine the final location determination solution comprise instructions configured to cause the processor to determine the ambiguous location determination solution whose corresponding second ranking has a largest total number of matches with the first ranking.
17. The computer readable medium of claim 15 wherein the program further comprises instructions configured to cause the processor to weight each match according to an estimated distance of a base station from which the corresponding ranging signal is received.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/754,348 US20100190511A1 (en) | 2004-11-17 | 2010-04-05 | Method for ambiguity resolution in location determination |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US62901104P | 2004-11-17 | 2004-11-17 | |
US11/034,705 US7881905B2 (en) | 2004-11-17 | 2005-01-12 | Method for ambiguity resolution in location determination |
US12/754,348 US20100190511A1 (en) | 2004-11-17 | 2010-04-05 | Method for ambiguity resolution in location determination |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/034,705 Division US7881905B2 (en) | 2004-11-17 | 2005-01-12 | Method for ambiguity resolution in location determination |
Publications (1)
Publication Number | Publication Date |
---|---|
US20100190511A1 true US20100190511A1 (en) | 2010-07-29 |
Family
ID=35892584
Family Applications (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/034,705 Expired - Fee Related US7881905B2 (en) | 2004-11-17 | 2005-01-12 | Method for ambiguity resolution in location determination |
US12/754,348 Abandoned US20100190511A1 (en) | 2004-11-17 | 2010-04-05 | Method for ambiguity resolution in location determination |
US12/754,360 Abandoned US20100190512A1 (en) | 2004-11-17 | 2010-04-05 | Method for ambiguity resolution in location determination |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/034,705 Expired - Fee Related US7881905B2 (en) | 2004-11-17 | 2005-01-12 | Method for ambiguity resolution in location determination |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US12/754,360 Abandoned US20100190512A1 (en) | 2004-11-17 | 2010-04-05 | Method for ambiguity resolution in location determination |
Country Status (6)
Country | Link |
---|---|
US (3) | US7881905B2 (en) |
EP (3) | EP2211197A3 (en) |
JP (2) | JP5065040B2 (en) |
CA (1) | CA2587649A1 (en) |
RU (1) | RU2384859C2 (en) |
WO (1) | WO2006055790A2 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100190512A1 (en) * | 2004-11-17 | 2010-07-29 | Qualcomm Incoporated | Method for ambiguity resolution in location determination |
CN102368234A (en) * | 2011-09-27 | 2012-03-07 | 邵阳神风动力制造有限责任公司 | Method for converting PDF (portable document format) documents into DWG (drawing) documents |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8626503B2 (en) | 2005-07-14 | 2014-01-07 | Erik Gosuinus Petrus Schuijers | Audio encoding and decoding |
US8089399B2 (en) * | 2008-06-06 | 2012-01-03 | Skyhook Wireless, Inc. | System and method for refining a WLAN-PS estimated location using satellite measurements in a hybrid positioning system |
KR102051965B1 (en) * | 2012-06-28 | 2019-12-04 | 한국전자통신연구원 | Apparatus and method for measuring position |
CA2903460A1 (en) | 2013-03-15 | 2014-09-25 | Nextnav, Llc | Techniques to improve the performance of a fixed, timing-based radio positioning network using external assistance information |
CA3001490C (en) * | 2015-11-09 | 2024-03-19 | Wiser Systems, Inc. | Methods for synchronizing multiple devices and determining location based on the synchronized devices |
CN110139303B (en) * | 2019-04-23 | 2021-06-11 | 四川九洲电器集团有限责任公司 | Rapid simulation method and device for equivalent signal level TOA measurement |
JP7067541B2 (en) | 2019-10-31 | 2022-05-16 | Jfeスチール株式会社 | Rolling mill control method and control device |
Citations (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5293642A (en) * | 1990-12-19 | 1994-03-08 | Northern Telecom Limited | Method of locating a mobile station |
US5623414A (en) * | 1995-01-24 | 1997-04-22 | Massachusetts Inst Technology | Clock-aided satellite navigation receiver system for enhanced position estimation and integrity monitoring |
US5740048A (en) * | 1992-08-07 | 1998-04-14 | Abel; Jonathan S. | Method and apparatus for GPS positioning, filtering and integration |
US5878033A (en) * | 1994-03-22 | 1999-03-02 | Nortel Marta Cellular | Process and equipment for broadcasting messages to mobile radio communication stations |
US6114992A (en) * | 1997-05-22 | 2000-09-05 | Conexant Systems, Inc. | Satellite acquisition and measurement system and process |
US6249252B1 (en) * | 1996-09-09 | 2001-06-19 | Tracbeam Llc | Wireless location using multiple location estimators |
US6263208B1 (en) * | 1999-05-28 | 2001-07-17 | Lucent Technologies Inc. | Geolocation estimation method for CDMA terminals based on pilot strength measurements |
US6266537B1 (en) * | 1998-03-27 | 2001-07-24 | Nec Corporation | Radio communication system |
US6289280B1 (en) * | 1999-12-10 | 2001-09-11 | Qualcomm Incorporated | Method and apparatus for determining an algebraic solution to GPS terrestrial hybrid location system equations |
US20010022558A1 (en) * | 1996-09-09 | 2001-09-20 | Tracbeam Llc | Wireless location using signal fingerprinting |
US20020135462A1 (en) * | 2001-03-26 | 2002-09-26 | Nels Michael J. | Method for proactive personalization and alerts to mobile devices |
US6496489B1 (en) * | 1993-11-01 | 2002-12-17 | Telefonaktiebolaget Lm Ericsson | Method and apparatus for distinguishing between a digital control channel and a digital traffic channel in a radiocommunication system |
US6496701B1 (en) * | 1998-08-25 | 2002-12-17 | Lucent Technologies Inc. | Pattern-recognition-based geolocation |
US20030040300A1 (en) * | 2001-08-27 | 2003-02-27 | Alcatel | System of interoperability between MMS messages and SMS/EMS messages and an associated exchange method |
US6542121B1 (en) * | 1999-03-03 | 2003-04-01 | Qinetiq Limited | Sensor system |
US6564065B1 (en) * | 1999-04-20 | 2003-05-13 | Lucent Technologies Inc. | Bayesian-update based location prediction method for CDMA systems |
US20030119496A1 (en) * | 2001-10-29 | 2003-06-26 | Peter Gaal | Base station time calibration using position measurement data sent by mobile stations during regular position location sessions |
US20030146871A1 (en) * | 1998-11-24 | 2003-08-07 | Tracbeam Llc | Wireless location using signal direction and time difference of arrival |
US6889053B1 (en) * | 1999-07-26 | 2005-05-03 | Lucent Technologies Inc. | Likelihood-based geolocation prediction algorithms for CDMA systems using pilot strength measurements |
US20100190512A1 (en) * | 2004-11-17 | 2010-07-29 | Qualcomm Incoporated | Method for ambiguity resolution in location determination |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1309747A1 (en) | 1985-05-06 | 1999-09-10 | Г.Н. Громов | SYSTEM OF DETERMINATION OF THE LOCATION OF THE AIRCRAFT |
US4963889A (en) | 1989-09-26 | 1990-10-16 | Magnavox Government And Industrial Electronics Company | Method and apparatus for precision attitude determination and kinematic positioning |
RU2060513C1 (en) | 1992-10-06 | 1996-05-20 | Олег Евгениевич Антонов | Method for identifying object location |
RU2098297C1 (en) | 1996-01-05 | 1997-12-10 | Герман Федорович Иоилев | Method of and system for detection and identification of object |
CN1208630C (en) * | 1999-01-08 | 2005-06-29 | 真实定位公司 | Bandwidth synthesis for wireless location system |
CA2296812A1 (en) * | 1999-02-17 | 2000-08-17 | Lucent Technologies Inc. | Method for combining multiple measurements to determine the position of a mobile transceiver |
JP2000244967A (en) * | 1999-02-24 | 2000-09-08 | Mitsubishi Electric Corp | Mobile communication system, mobile unit and base station configuring the system and method for detecting position of the mobile unit in the system |
RU2152050C1 (en) | 1999-12-01 | 2000-06-27 | Закрытое акционерное общество "НПО Космического Приборостроения" | Satellite navigation system for detection of object position |
JP2003534531A (en) * | 2000-01-24 | 2003-11-18 | インテグリノーティックス コーポレイション | Multi-frequency pseudolite for carrier-based differential position determination |
RU2161317C1 (en) | 2000-05-11 | 2000-12-27 | Закрытое акционерное общество "НПО Космического Приборостроения" | System for high-accuracy localization of objects- navigational information users by navigational radio signals with authorized access in regime of differential corrections |
JP4396801B2 (en) * | 2001-02-27 | 2010-01-13 | パイオニア株式会社 | POSITIONING DEVICE, POSITIONING METHOD, AND POSITIONING SYSTEM |
JP2003075524A (en) * | 2001-08-31 | 2003-03-12 | Denso Corp | Mobile communication terminal device, program and position server therefor |
BRPI0215377B1 (en) | 2001-12-27 | 2016-05-17 | Qualcomm Inc | creating and using base station almanac information in a wireless communication system having a position locating capability |
FI113092B (en) * | 2002-05-31 | 2004-02-27 | Ekahau Oy | Measures of position differences and applications |
US7257411B2 (en) * | 2002-12-27 | 2007-08-14 | Ntt Docomo, Inc. | Selective fusion location estimation (SELFLOC) for wireless access technologies |
EP1445970B1 (en) | 2003-02-05 | 2009-04-01 | Cambridge Positioning Systems Limited | A method and system for locating a mobile radio receiver in a radio system with multiple tranmitters |
EP1683285A2 (en) * | 2003-05-23 | 2006-07-26 | Symbol Technologies, Inc. | Self calibration of signal strength location system |
-
2005
- 2005-01-12 US US11/034,705 patent/US7881905B2/en not_active Expired - Fee Related
- 2005-11-17 EP EP10162058A patent/EP2211197A3/en not_active Withdrawn
- 2005-11-17 EP EP10162057A patent/EP2211196A3/en not_active Withdrawn
- 2005-11-17 RU RU2007122504/09A patent/RU2384859C2/en not_active IP Right Cessation
- 2005-11-17 CA CA002587649A patent/CA2587649A1/en not_active Abandoned
- 2005-11-17 EP EP05825049A patent/EP1812806A2/en not_active Withdrawn
- 2005-11-17 JP JP2007543275A patent/JP5065040B2/en not_active Expired - Fee Related
- 2005-11-17 WO PCT/US2005/041847 patent/WO2006055790A2/en active Application Filing
-
2010
- 2010-04-05 US US12/754,348 patent/US20100190511A1/en not_active Abandoned
- 2010-04-05 US US12/754,360 patent/US20100190512A1/en not_active Abandoned
-
2011
- 2011-02-03 JP JP2011021813A patent/JP5129356B2/en not_active Expired - Fee Related
Patent Citations (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5293642A (en) * | 1990-12-19 | 1994-03-08 | Northern Telecom Limited | Method of locating a mobile station |
US5740048A (en) * | 1992-08-07 | 1998-04-14 | Abel; Jonathan S. | Method and apparatus for GPS positioning, filtering and integration |
US6496489B1 (en) * | 1993-11-01 | 2002-12-17 | Telefonaktiebolaget Lm Ericsson | Method and apparatus for distinguishing between a digital control channel and a digital traffic channel in a radiocommunication system |
US5878033A (en) * | 1994-03-22 | 1999-03-02 | Nortel Marta Cellular | Process and equipment for broadcasting messages to mobile radio communication stations |
US5623414A (en) * | 1995-01-24 | 1997-04-22 | Massachusetts Inst Technology | Clock-aided satellite navigation receiver system for enhanced position estimation and integrity monitoring |
US7525484B2 (en) * | 1996-09-09 | 2009-04-28 | Tracbeam Llc | Gateway and hybrid solutions for wireless location |
US20010022558A1 (en) * | 1996-09-09 | 2001-09-20 | Tracbeam Llc | Wireless location using signal fingerprinting |
US6249252B1 (en) * | 1996-09-09 | 2001-06-19 | Tracbeam Llc | Wireless location using multiple location estimators |
US6114992A (en) * | 1997-05-22 | 2000-09-05 | Conexant Systems, Inc. | Satellite acquisition and measurement system and process |
US6266537B1 (en) * | 1998-03-27 | 2001-07-24 | Nec Corporation | Radio communication system |
US6496701B1 (en) * | 1998-08-25 | 2002-12-17 | Lucent Technologies Inc. | Pattern-recognition-based geolocation |
US20030146871A1 (en) * | 1998-11-24 | 2003-08-07 | Tracbeam Llc | Wireless location using signal direction and time difference of arrival |
US6542121B1 (en) * | 1999-03-03 | 2003-04-01 | Qinetiq Limited | Sensor system |
US6564065B1 (en) * | 1999-04-20 | 2003-05-13 | Lucent Technologies Inc. | Bayesian-update based location prediction method for CDMA systems |
US6263208B1 (en) * | 1999-05-28 | 2001-07-17 | Lucent Technologies Inc. | Geolocation estimation method for CDMA terminals based on pilot strength measurements |
US6889053B1 (en) * | 1999-07-26 | 2005-05-03 | Lucent Technologies Inc. | Likelihood-based geolocation prediction algorithms for CDMA systems using pilot strength measurements |
US6289280B1 (en) * | 1999-12-10 | 2001-09-11 | Qualcomm Incorporated | Method and apparatus for determining an algebraic solution to GPS terrestrial hybrid location system equations |
US20020135462A1 (en) * | 2001-03-26 | 2002-09-26 | Nels Michael J. | Method for proactive personalization and alerts to mobile devices |
US20030040300A1 (en) * | 2001-08-27 | 2003-02-27 | Alcatel | System of interoperability between MMS messages and SMS/EMS messages and an associated exchange method |
US20030119496A1 (en) * | 2001-10-29 | 2003-06-26 | Peter Gaal | Base station time calibration using position measurement data sent by mobile stations during regular position location sessions |
US7006834B2 (en) * | 2001-10-29 | 2006-02-28 | Qualcomm Incorporated | Base station time calibration using position measurement data sent by mobile stations during regular position location sessions |
US20100190512A1 (en) * | 2004-11-17 | 2010-07-29 | Qualcomm Incoporated | Method for ambiguity resolution in location determination |
US7881905B2 (en) * | 2004-11-17 | 2011-02-01 | Qualcomm Incorporated | Method for ambiguity resolution in location determination |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100190512A1 (en) * | 2004-11-17 | 2010-07-29 | Qualcomm Incoporated | Method for ambiguity resolution in location determination |
CN102368234A (en) * | 2011-09-27 | 2012-03-07 | 邵阳神风动力制造有限责任公司 | Method for converting PDF (portable document format) documents into DWG (drawing) documents |
Also Published As
Publication number | Publication date |
---|---|
EP2211197A2 (en) | 2010-07-28 |
CA2587649A1 (en) | 2006-05-26 |
US7881905B2 (en) | 2011-02-01 |
RU2384859C2 (en) | 2010-03-20 |
EP2211197A3 (en) | 2010-12-15 |
EP2211196A3 (en) | 2010-12-15 |
EP2211196A2 (en) | 2010-07-28 |
US20100190512A1 (en) | 2010-07-29 |
JP2008521009A (en) | 2008-06-19 |
US20060106573A1 (en) | 2006-05-18 |
JP5065040B2 (en) | 2012-10-31 |
JP5129356B2 (en) | 2013-01-30 |
WO2006055790A3 (en) | 2007-01-04 |
RU2007122504A (en) | 2008-12-27 |
WO2006055790A2 (en) | 2006-05-26 |
JP2011107158A (en) | 2011-06-02 |
EP1812806A2 (en) | 2007-08-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20100190511A1 (en) | Method for ambiguity resolution in location determination | |
Shah et al. | 3D weighted centroid algorithm & RSSI ranging model strategy for node localization in WSN based on smart devices | |
McGuire et al. | Data fusion of power and time measurements for mobile terminal location | |
DK1514130T3 (en) | Probabilistic model for a positioning technique | |
CN101536591B (en) | Extended clustering for improved positioning | |
Tayebi et al. | The application of ray-tracing to mobile localization using the direction of arrival and received signal strength in multipath indoor environments | |
Kul et al. | IEEE 802.11 WLAN based real time indoor positioning: Literature survey and experimental investigations | |
CN102695269B (en) | A kind of positioning correction method, relevant apparatus and system | |
WO2010077819A1 (en) | Method for position estimation using generalized error distributions | |
EP3039446A1 (en) | Location detection system and method | |
CN108882149B (en) | NLOS compensation positioning method of distance correlation probability | |
CN102325370A (en) | High-precision three-dimensional positioner for wireless sensor network node | |
KR100714120B1 (en) | Method and apparatus for locating a remote unit | |
Li et al. | Cramer-rao lower bound analysis of data fusion for fingerprinting localization in non-line-of-sight environments | |
KR101258832B1 (en) | Method for determining location using the access point and apparatus thereof | |
Pradhan et al. | Mathematical analysis of line intersection algorithm for TOA trilateration method | |
Zhao et al. | Signal path reckoning localization method in multipath environment | |
Filípek et al. | Indoor localization based on beacons and calculated by particle filter | |
Go et al. | An efficient non-line-of-sight error mitigation method for TOA measurement in indoor environments | |
Zhang et al. | Improved KNN algorithm with historical information fusion for indoor positioning | |
Başak et al. | Accurate indoor localization with optimized fingerprinting algorithm | |
Wu et al. | A Regional Positioning Algorithm for Non-line-of-sight (NOLS) Propagation | |
Shixun et al. | Single base station hybrid TOA/AOD/AOA localization algorithms with the synchronization error in dense multipath environment | |
CN103149551A (en) | Convex optimization range-based localization method based on linear combination | |
Khalaf-Allah | Nonparametric Bayesian filtering for location estimation, position tracking, and global localization of mobile terminals in outdoor wireless environments |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: QUALCOMM INCORPORATED, CALIFORNIA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:SEIBERT, CRISTINA A;REEL/FRAME:024236/0677 Effective date: 20100414 |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO PAY ISSUE FEE |