# Patent

## Title
Method and system for radio map filtering via adaptive clustering

## Abstract

A method for estimating a location of a wireless device in a wireless local network is provided. The method includes forming a first set comprised of the signal strength received from access points that the wireless device received a signal from and an indicator of no signal strength measured for access points that the wireless device did not receive a signal from. Next, a scan subset can be formed comprised of access points in the first set that has associated signal strength. Next, a cluster comprised of the calibration points can be formed based on the scan subset. A distance between the first set and each of the calibration point in the cluster can be calculated. Then, the smallest distance can be selected as the location estimate.

## Background
TECHNICAL FIELD OF THE INVENTION
This invention relates to the field of real time location systems and, more particularly, to a method and system for location estimation in wireless local area networks via adaptive clustering.
BACKGROUND OF THE INVENTION
Object location based on the received signal strength (RSSI) of a signal transmitted from reference objects has applications in many different fields of endeavor. For example, objects having small transceiver affixed to them can be tracked throughout a warehouse. Also, user of wireless devices can determine their locations in an environment using the RSSI received from reference devices.
One prior art method for locating an object RSSI calculates the signal space distance from calibration points within a wireless location system. In one particular prior art method, the location can be determined by calculating a distance in signal space between the current sample point and calibration points.
As an example, an area has a total of four reference-devices. At a given point in time, a wireless device will scan the wireless channels to determine the RSSI received from each reference device. The result can be written as a set of n RSSI measurements. For example, if the set of the reference device RSSIs at a given location can be called the sample, S, then the set of RSSI measured from each reference point will give S: {AP1, AP2, AP3, AP4}. The sample set can then be compared to the RSSIs of various calibration points. Each of the calibration points will have their own set of received signal strengths measurements from each reference device. The comparison can be done by calculating the Euclidean distance from each calibration point to the sample using the measured signal strengths and choosing the calibration point closest (having the smallest separation) to the sample set. The sample can be estimated to be located near that calibration point. This algorithm is discussed in “RADAR: An in-building RF based user location and tracking system”, by Parumuir Bahl and Venkata N. Padmanabhan, and published in Preceding of INFOCOM, 2000.
This method has several drawbacks. First, it assumes that the receiver can receive signal from all transmitters at all the time, thus failing to handle the situation where only a subset of the transmitters can be heard by the receiver, which is the typical case in large scale network. Second, since it fails to filter calibration points, all calibration points must be checked against the sample. As the number of reference devices and calibration points increase, the computational load and power consumption increases as well. Also, it is possible that the location algorithm can be thrown off by certain calibration points if all calibration points are used. What is need is a method and system for location estimation on wireless local networks via adaptive clustering.
SUMMARY OF THE INVENTION
A method for estimating a location of a wireless device in a wireless network and calibration points is provided in accordance with one exemplary embodiment of the present invention. The method includes forming a first set comprised of the signal strength received from reference devices that the wireless device received a signal from and an indicator of no signal strength measured for reference devices that the wireless device did not receive a signal from. Next, a scan subset can be formed comprised of reference devices in the first set that has associated signal strength. A cluster comprised of the calibration points can be formed based on the scan subset. A distance between the first set and each of the calibration point in the cluster can be calculated. Then, the smallest distance can be selected as the location estimate.
In one aspect of the present invention, the step of forming a cluster further comprises forming a cluster comprising calibration points that contain exactly the same reference devices as the scan subset.
In another exemplary embodiment of the present invention, a system for estimating the location of a wireless device in a wireless network is disclosed. The system includes a plurality of reference devices distributed throughout the wireless network. The system also includes a plurality of calibration points distributed throughout the wireless network. A radio map can be formed from the signal strength received from the reference devices at each of the calibration points. The wireless device can be located in the network and can be operable to form a scan set comprising the signal strength received from each of the plurality of reference devices, and form a scan subset of the scan set comprising reference devices in the scan set that have an associated signal strength. The wireless device can be further operable to form a cluster set comprising cluster points selected from the radio map based on the scan subset. The location of the wireless device can be estimated by calculating the distance between scan subset and each of the calibration points in the cluster subset and choosing the calibration point associated with the smallest distance as the location estimate.



BRIEF DESCRIPTION OF THE DRAWINGS
The present invention will hereinafter be described in conjunction with the following drawing figures, wherein like numerals denote like elements, and:
FIG. 1 illustrates multiple reference devices and calibration points in a wireless network in accordance with one exemplary embodiment of the present invention;
FIG. 2 is a flowchart of a clustering method in accordance with an exemplary embodiment of the present invention; and
FIG. 3 is a flowchart of an adaptive clustering method in accordance with an exemplary embodiment of the present invention.



DETAILED DESCRIPTION OF THE DRAWINGS
The following detailed description is merely exemplary in nature and is not intended to limit the invention or the application and uses of the invention. Furthermore, there is no intention to be bound by any expressed or implied theory presented in the preceding technical field, background, brief summary or the following detailed description
A method for determining the location of an object in a wireless network using a clustering method in which any calibration point not in the same cluster is not used in the calculation to determine location. Thresholding techniques can be used to further improve the calculation.
An exemplary environment for the use of the present invention is illustrated in FIG. 1. FIG. 1 illustrates an area 102, which can be any area 102, enclosed or otherwise, where a number of wireless reference devices 104 are deployed, typically as part of a wireless network. Also provided are a number of calibration points 106. A wireless device 108, which can be any mobile device capable of wireless communications, is within the area 102.
The reference devices 104 transmit signals to aid in the location of objects within the area covered by the reference devices. In one embodiment, reference device 104 can be an access point in a wireless local area network.
Wireless device 108 can be any device capable of wireless communications. For example, wireless device 108 can be a radio frequency ID tag, a wireless handled computer device such as a personal digital assistant with integrated wireless communication capabilities, a laptop computer with wireless communication ability, a wireless data terminal and the like. In one embodiment, wireless device 108 includes a receiver or transceiver for receiving signals from reference devices, a memory for saving data needed to determine a location and a processor to perform such calculations.
Prior to determining the location of the wireless device 108, a radio map is formed. This can be done by taking a survey at each calibration point 106. To take a survey, a wireless device can be used at each calibration point 106 to listen for transmissions from the reference devices 104. The signal strength from each reference device 104 heard at that calibration point can be recorded. This process can be repeated several times and the signal strength values are typically averaged. After the survey is completed, the RSSI values for each of the reference devices 104 as measured at each calibration point 106 and the coordinate of each calibration point 106 are stored at any location in the network. For example, the radio map can be stored at the wireless device 108, at the reference device 104, at a server computer coupled to the wireless network (not pictured) or another location.
Once the survey has been completed at every calibration points the radio map is compete. Then, a wireless device 108 that is to be located can perform a scan. A scan is like a survey that is done by the wireless device 108 but with an unknown location. In the scan, the wireless device records a signal strength for each reference device 104 in the network. If the wireless device 108 did not receive a signal from a reference device 104, the strength can be set as not available (NA) or no available number (NAN). The signal strength for each reference device 104 as recorded by a wireless device 108 can be arranged as a set of signal strength. Once this information is obtained, the signal strengths found in the scan can be compared with the signal strengths of the calibration point 106 to find an estimated location of the wireless device 108. The calculation can be done at the wireless device 108, the reference device 104, at a server computer or another location.
Prior to discussion of the location methods of the present, two lemmas need to be established. The first lemma is that the variance of received signal strength when expressed in decibels is a constant. The second lemma is that the higher the received signal strength from a specific access point, the higher the probability that the data (typically in terms of packets) can be received from this access point.
To establish that the variance of RSSI (in decibels) is a constant number, it is first noted that in indoor mobile radio channels, the well known Rayleigh distribution can be used to describe the statistical time varying nature of the received envelope of a flat fading signal. Based on this assumption, we can calculate the variance of RSSI by use of the probability density function of the Rayleigh distribution.
The probability density function (pdf) of received envelope v is given by:




  
    
      
        
          
            f
            ⁡
            
              (
              v
              )
            
          
          =
          
            
              1
              
                σ
                2
              
            
            ⁢
            
              exp
              ⁡
              
                (
                
                  -
                  
                    
                      v
                      2
                    
                    
                      2
                      ⁢
                      
                        σ
                        2
                      
                    
                  
                
                )
              
            
          
        
        ,
        
          0
          ≤
          v
          ≤
          ∞
        
      
    
    
      
        (
        1
        )
      
    
  




Then the pdf of received power, P, where P=v2, is the following exponential:




  
    
      
        
          
            
              
                
                  f
                  ⁡
                  
                    (
                    P
                    )
                  
                
                =
                
                  
                    1
                    
                      σ
                      2
                    
                  
                  ⁢
                  
                    exp
                    ⁡
                    
                      (
                      
                        -
                        
                          P
                          
                            2
                            ⁢
                            
                              σ
                              2
                            
                          
                        
                      
                      )
                    
                  
                
              
              ,
              
                0
                ≤
                P
                ≤
                ∞
              
            
          
        
        
          
            
              =
              
                
                  1
                  μ
                
                ⁢
                
                  exp
                  ⁡
                  
                    (
                    
                      -
                      
                        P
                        μ
                      
                    
                    )
                  
                
              
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




The received signal strength indication (RSSI) is measured in on a decibel scale, where PdB=10 log 10(P/10−3). The mean and variance of PdB is as following (where E(P)=μ and Var(P)=μ2):




  
    
      
        
          
            
              
                E
                ⁡
                
                  (
                  
                    P
                    dB
                  
                  )
                
              
              =
              
              ⁢
              
                
                  ∫
                  0
                  ∞
                
                ⁢
                
                  10
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      log
                      10
                    
                    ⁡
                    
                      (
                      
                        1000
                        ⁢
                        P
                      
                      )
                    
                  
                  ⁢
                  
                    f
                    ⁡
                    
                      (
                      P
                      )
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    d
                    P
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  10
                  ⁢
                  
                    
                      ∫
                      0
                      ∞
                    
                    ⁢
                    
                      
                        
                          log
                          10
                        
                        ⁡
                        
                          (
                          P
                          )
                        
                      
                      ⁢
                      
                        f
                        ⁡
                        
                          (
                          P
                          )
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        d
                        P
                      
                    
                  
                
                +
                
                  30
                  ⁢
                  
                    
                      ∫
                      0
                      ∞
                    
                    ⁢
                    
                      
                        f
                        ⁡
                        
                          (
                          P
                          )
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        d
                        P
                      
                    
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  10
                  ⁢
                  
                    
                      ∫
                      0
                      ∞
                    
                    ⁢
                    
                      
                        log
                        10
                      
                      ⁢
                      
                        P
                        ·
                        
                          f
                          ⁡
                          
                            (
                            P
                            )
                          
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        d
                        P
                      
                    
                  
                
                +
                30
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  10
                  ⁢
                  
                    
                      ∫
                      0
                      ∞
                    
                    ⁢
                    
                      
                        log
                        10
                      
                      ⁢
                      
                        P
                        ·
                        
                          1
                          μ
                        
                      
                      ⁢
                      
                        exp
                        ⁡
                        
                          (
                          
                            -
                            
                              P
                              μ
                            
                          
                          )
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        d
                        P
                      
                    
                  
                
                +
                30
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  10
                  ⁢
                  
                    
                      ∫
                      0
                      ∞
                    
                    ⁢
                    
                      
                        
                          
                            log
                            10
                          
                          ⁡
                          
                            (
                            
                              
                                P
                                μ
                              
                              ⁢
                              μ
                            
                            )
                          
                        
                        ·
                        
                          exp
                          ⁡
                          
                            (
                            
                              -
                              
                                P
                                μ
                              
                            
                            )
                          
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        d
                        
                          P
                          μ
                        
                      
                    
                  
                
                +
                30
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  10
                  ⁢
                  
                    
                      ∫
                      0
                      ∞
                    
                    ⁢
                    
                      
                        
                          
                            log
                            10
                          
                          ⁡
                          
                            (
                            
                              y
                              ⁢
                              
                                
                              
                              ⁢
                              μ
                            
                            )
                          
                        
                        ·
                        
                          exp
                          ⁡
                          
                            (
                            
                              -
                              y
                            
                            )
                          
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        d
                        y
                      
                    
                  
                
                +
                30
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  10
                  ⁢
                  
                    
                      ∫
                      0
                      ∞
                    
                    ⁢
                    
                      
                        
                          
                            log
                            10
                          
                          ⁡
                          
                            (
                            y
                            )
                          
                        
                        ·
                        
                          exp
                          ⁡
                          
                            (
                            
                              -
                              y
                            
                            )
                          
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        d
                        y
                      
                    
                  
                
                +
                
                  10
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      log
                      10
                    
                    ⁡
                    
                      (
                      μ
                      )
                    
                  
                
                +
                30
              
            
          
        
      
    
    
      
        (
        3
        )
      
    
  
  
    
      
        
          
            
              
                E
                ⁡
                
                  (
                  
                    P
                    dB
                    2
                  
                  )
                
              
              =
              
              ⁢
              
                
                  ∫
                  0
                  ∞
                
                ⁢
                
                  100
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      log
                      10
                      2
                    
                    ⁡
                    
                      (
                      
                        1000
                        ⁢
                        P
                      
                      )
                    
                  
                  ⁢
                  
                    f
                    ⁡
                    
                      (
                      P
                      )
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    d
                    P
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                100
                ⁢
                
                  
                    ∫
                    0
                    ∞
                  
                  ⁢
                  
                    
                      (
                      
                        9
                        +
                        
                          6
                          ⁢
                          
                            
                          
                          ⁢
                          
                            log
                            10
                          
                          ⁢
                          P
                        
                        +
                        
                          
                            log
                            10
                            2
                          
                          ⁢
                          P
                        
                      
                      )
                    
                    ⁢
                    
                      f
                      ⁡
                      
                        (
                        P
                        )
                      
                    
                    ⁢
                    
                      
                    
                    ⁢
                    
                      d
                      P
                    
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                900
                +
                
                  600
                  ⁢
                  
                    
                      ∫
                      0
                      ∞
                    
                    ⁢
                    
                      
                        log
                        10
                      
                      ⁢
                      
                        Pf
                        ⁡
                        
                          (
                          P
                          )
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        d
                        P
                      
                    
                  
                
                +
                
                  100
                  ⁢
                  
                    
                      ∫
                      0
                      ∞
                    
                    ⁢
                    
                      
                        (
                        
                          
                            log
                            10
                            2
                          
                          ⁢
                          P
                        
                        )
                      
                      ⁢
                      
                        f
                        ⁡
                        
                          (
                          P
                          )
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        d
                        P
                      
                    
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                900
                +
                
                  600
                  ⁢
                  
                    
                      ∫
                      0
                      ∞
                    
                    ⁢
                    
                      
                        
                          log
                          10
                        
                        ⁡
                        
                          (
                          
                            
                              P
                              μ
                            
                            ·
                            μ
                          
                          )
                        
                      
                      ⁢
                      
                        exp
                        ⁡
                        
                          (
                          
                            -
                            
                              P
                              μ
                            
                          
                          )
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        d
                        
                          P
                          μ
                        
                      
                    
                  
                
                +
                100
              
            
          
        
        
          
            
              
              ⁢
              
                
                  ∫
                  0
                  ∞
                
                ⁢
                
                  
                    
                      (
                      
                        
                          
                            log
                            10
                          
                          ⁢
                          
                            P
                            μ
                          
                        
                        +
                        
                          
                            log
                            10
                          
                          ⁢
                          μ
                        
                      
                      )
                    
                    2
                  
                  ⁢
                  
                    f
                    ⁡
                    
                      (
                      P
                      )
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    d
                    P
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                900
                +
                
                  600
                  ⁡
                  
                    [
                    
                      
                        
                          ∫
                          0
                          ∞
                        
                        ⁢
                        
                          
                            (
                            
                              
                                log
                                10
                              
                              ⁢
                              y
                            
                            )
                          
                          ⁢
                          
                            exp
                            ⁡
                            
                              (
                              
                                -
                                y
                              
                              )
                            
                          
                          ⁢
                          
                            
                          
                          ⁢
                          
                            d
                            y
                          
                        
                      
                      +
                      
                        
                          log
                          10
                        
                        ⁢
                        μ
                      
                    
                    ]
                  
                
                +
              
            
          
        
        
          
            
              
              ⁢
              
                100
                ⁡
                
                  [
                  
                    
                      ∫
                      0
                      ∞
                    
                    ⁢
                    
                      
                        
                          (
                          
                            
                              
                                log
                                10
                              
                              ⁢
                              
                                P
                                μ
                              
                            
                            +
                            
                              
                                log
                                10
                              
                              ⁢
                              μ
                            
                          
                          )
                        
                        2
                      
                      ⁢
                      
                        exp
                        ⁡
                        
                          (
                          
                            -
                            
                              P
                              μ
                            
                          
                          )
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        d
                        
                          P
                          μ
                        
                      
                    
                  
                  ]
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                900
                +
                
                  600
                  ⁡
                  
                    [
                    
                      
                        
                          ∫
                          0
                          ∞
                        
                        ⁢
                        
                          
                            (
                            
                              
                                log
                                10
                              
                              ⁢
                              y
                            
                            )
                          
                          ⁢
                          
                            exp
                            ⁡
                            
                              (
                              
                                -
                                y
                              
                              )
                            
                          
                          ⁢
                          
                            
                          
                          ⁢
                          
                            d
                            y
                          
                        
                      
                      +
                      
                        
                          log
                          10
                        
                        ⁢
                        μ
                      
                    
                    ]
                  
                
                +
              
            
          
        
        
          
            
              
              ⁢
              
                100
                ⁡
                
                  [
                  
                    
                      ∫
                      0
                      ∞
                    
                    ⁢
                    
                      
                        
                          (
                          
                            
                              
                                log
                                10
                              
                              ⁢
                              y
                            
                            +
                            
                              
                                log
                                10
                              
                              ⁢
                              μ
                            
                          
                          )
                        
                        2
                      
                      ⁢
                      
                        exp
                        ⁡
                        
                          (
                          
                            -
                            y
                          
                          )
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        d
                        y
                      
                    
                  
                  ]
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                900
                +
                
                  600
                  ⁡
                  
                    [
                    
                      
                        
                          ∫
                          0
                          ∞
                        
                        ⁢
                        
                          
                            (
                            
                              
                                log
                                10
                              
                              ⁢
                              y
                            
                            )
                          
                          ⁢
                          
                            exp
                            ⁡
                            
                              (
                              
                                -
                                y
                              
                              )
                            
                          
                          ⁢
                          
                            
                          
                          ⁢
                          
                            d
                            y
                          
                        
                      
                      +
                      
                        
                          log
                          10
                        
                        ⁢
                        μ
                      
                    
                    ]
                  
                
                +
              
            
          
        
        
          
            
              
              ⁢
              
                100
                [
                
                  
                    
                      ∫
                      0
                      ∞
                    
                    ⁢
                    
                      
                        
                          (
                          
                            
                              log
                              10
                            
                            ⁢
                            y
                          
                          )
                        
                        2
                      
                      ⁢
                      
                        exp
                        ⁡
                        
                          (
                          
                            -
                            y
                          
                          )
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        d
                        y
                      
                    
                  
                  +
                  
                    
                      (
                      
                        
                          log
                          10
                        
                        ⁢
                        μ
                      
                      )
                    
                    2
                  
                  +
                
              
            
          
        
        
          
            
              
              ⁢
              
                2
                ⁢
                
                  (
                  
                    
                      log
                      10
                    
                    ⁢
                    μ
                  
                  )
                
                ⁢
                
                  
                    ∫
                    0
                    ∞
                  
                  ⁢
                  
                    
                      (
                      
                        
                          log
                          10
                        
                        ⁢
                        y
                      
                      )
                    
                    ⁢
                    
                      exp
                      ⁡
                      
                        (
                        
                          -
                          y
                        
                        )
                      
                    
                    ⁢
                    
                      
                    
                    ⁢
                    
                      d
                      y
                    
                  
                
              
              ]
            
          
        
      
    
    
      
        (
        4
        )
      
    
  
  
    
      
        
          
            
              
                Var
                ⁡
                
                  (
                  PdB
                  )
                
              
              =
              
              ⁢
              
                
                  E
                  ⁡
                  
                    (
                    PdB2
                    )
                  
                
                -
                
                  
                    [
                    
                      E
                      ⁡
                      
                        (
                        PdB
                        )
                      
                    
                    ]
                  
                  ⁢
                  2
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                900
                +
                
                  600
                  ⁡
                  
                    [
                    
                      
                        
                          ∫
                          0
                          ∞
                        
                        ⁢
                        
                          
                            (
                            
                              
                                log
                                10
                              
                              ⁢
                              y
                            
                            )
                          
                          ⁢
                          
                            exp
                            ⁡
                            
                              (
                              
                                -
                                y
                              
                              )
                            
                          
                          ⁢
                          
                            
                          
                          ⁢
                          
                            d
                            y
                          
                        
                      
                      +
                      
                        
                          log
                          10
                        
                        ⁢
                        μ
                      
                    
                    ]
                  
                
                +
              
            
          
        
        
          
            
              
              ⁢
              
                100
                [
                
                  
                    
                      ∫
                      0
                      ∞
                    
                    ⁢
                    
                      
                        
                          (
                          
                            
                              log
                              10
                            
                            ⁢
                            y
                          
                          )
                        
                        2
                      
                      ⁢
                      
                        exp
                        ⁡
                        
                          (
                          
                            -
                            y
                          
                          )
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        d
                        y
                      
                    
                  
                  +
                  
                    
                      (
                      
                        
                          log
                          10
                        
                        ⁢
                        μ
                      
                      )
                    
                    2
                  
                  +
                
              
            
          
        
        
          
            
              
                
                ⁢
                
                  2
                  ⁢
                  
                    (
                    
                      
                        log
                        10
                      
                      ⁢
                      μ
                    
                    )
                  
                  ⁢
                  
                    
                      ∫
                      0
                      ∞
                    
                    ⁢
                    
                      
                        (
                        
                          
                            log
                            10
                          
                          ⁢
                          y
                        
                        )
                      
                      ⁢
                      
                        exp
                        ⁡
                        
                          (
                          
                            -
                            y
                          
                          )
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        d
                        y
                      
                    
                  
                
                ]
              
              -
            
          
        
        
          
            
              
              ⁢
              
                
                  [
                  
                    
                      10
                      ⁢
                      
                        
                          ∫
                          0
                          ∞
                        
                        ⁢
                        
                          
                            
                              
                                log
                                10
                              
                              ⁡
                              
                                (
                                y
                                )
                              
                            
                            ·
                            
                              exp
                              ⁡
                              
                                (
                                
                                  -
                                  y
                                
                                )
                              
                            
                          
                          ⁢
                          
                            
                          
                          ⁢
                          
                            d
                            y
                          
                        
                      
                    
                    +
                    
                      10
                      ⁢
                      
                        
                      
                      ⁢
                      
                        
                          log
                          10
                        
                        ⁡
                        
                          (
                          μ
                          )
                        
                      
                    
                    +
                    30
                  
                  ]
                
                2
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  100
                  ⁢
                  
                    
                      ∫
                      0
                      ∞
                    
                    ⁢
                    
                      
                        
                          (
                          
                            
                              log
                              10
                            
                            ⁢
                            y
                          
                          )
                        
                        2
                      
                      ⁢
                      
                        exp
                        ⁡
                        
                          (
                          
                            -
                            y
                          
                          )
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        d
                        y
                      
                    
                  
                
                -
              
            
          
        
        
          
            
              
              ⁢
              
                
                  100
                  ⁡
                  
                    [
                    
                      
                        ∫
                        0
                        ∞
                      
                      ⁢
                      
                        
                          
                            
                              log
                              10
                            
                            ⁡
                            
                              (
                              y
                              )
                            
                          
                          ·
                          
                            exp
                            ⁡
                            
                              (
                              
                                -
                                y
                              
                              )
                            
                          
                        
                        ⁢
                        
                          
                        
                        ⁢
                        
                          d
                          y
                        
                      
                    
                    ]
                  
                
                2
              
            
          
        
      
    
    
      
        (
        5
        )
      
    
  




Therefore, although the mean and variance of the received power are correlated, the variance of RSSI is constant, assuming no noise.
Additionally, the higher the RSSI from a specific reference device 104, the higher the probability that the packet from this reference device 104 can be received. The relationship between the RSSI value and the possibility of receiving a packet is known as the wireless channel information. Since the variance of RSSI due to the Rayleigh fading is a constant number as shown above, and because the receiver cannot resolve the received signal below a certain threshold due to the limitation of the sensitivity of the receiver, it is obvious that the lower the RSSI, the higher the probability that the RSSI will below the receiver's sensitivity level.
The receiver's ability to receive a given signal depends on several factors such as the sensitivity of the receiver, the gain of the automatic gain control in the receiver and the like. Different receiver designs can affect the ability of a receiver.
Using the above information, location methods can be developed in accordance with the teachings of the present invention. The following definitions are used in the discussions that follow.
1) Calibration Sets: C={c1, c2 . . . cW}, calibration points within the area. Here W is the number of calibration points.
2) Current scan: S: {pr(1), . . . pr(M)}
3) Reference Devices Set: AP={AP1, AP2 . . . APM}, total M reference devices.
4) RSSI set: Pr(n)={pr(n,1), pr(n,2) . . . pr(n, M)}T. M is the number of reference devices. pr(n, m) is RSSI from the mth reference device.
5) Signal level threshold: Th1 or Th2 value below which the signal will be ignored.
6) RSSIs of calibration points: PC={PrC(1), PrC(2) . . . PrC(W)}, M×W matrix.
7) Clustering: The clustering hereafter refers to the grouping of calibration points that share some features. For example, c1:{11,nan,10,10} and c2={12,nan,10,11} have the same pattern F1, where F1: {Can see AP1, AP3, AP4}. The pattern is that signals from AP1, AP3 and AP4 can be received. Thus, c1 and c2 can be grouped into a cluster G1. The calibration point c3={13,nan,nan,nan} belongs to another cluster where only signals from AP1 can be received: F1: {Only AP1 can be seen}.
A first exemplary location method in accordance with the teachings of the present invention can be described with reference to the flowchart of FIG. 2. In a first step, a scan is completed (step 202). The scan produces a set S, where S: {pr(1), pr(2) . . . pr(M)}, and where M is the total number of access points. The wireless device may not receive signals from some reference devices and instead of reporting a signal strength, will report a not available (na) or no available number (nan).
Next, a subset of S, {APV}, can be formed comprising the reference devices from which a signal was received. That is, the reference devices that were seen in the scan are in set {APV} (step 204). Then for each of the calibration points, a subset, {APC(w)}, can be formed comprising the calibration points that report an RSSI from a reference device (step 206). Next, calibration points that contain exactly the same reference devices that were seen in the current scan and represented by the set {APV} are grouped together to form a group Gn (step 208). That is, Gn contains calibration points where {APC(w)}={APV}.
If Gn is a null set, then the comparison calculation will be done as before, with the scan set S: {pr(1) . . . pr(M)} and each of the calibration points being used to calculate a Euclidean distance between the scan set and the calibration points. The calibration point closest to the scan point using the Euclidean distance formula is declared as the estimated location (step 210).
If the set Gn is not empty, then the scan set S: {pr(1), . . . pr(n)} can be compared to the calibration points within the set Gn (step 212). Again, the calibration point that is closest in distance can be chosen as the estimate for the location of the mobile unit (step 214).
The above method is an improvement over prior methods because by performing calculations within a subset of the calibration points the computational load is greatly reduced. However, there is the possibility that in this method an error can occur due to searching in the wrong cluster. Searching in the wrong cluster can be caused by noise, network traffic and newly installed reference devices. The problem of searching the wrong cluster can have three main causes. First, the current scan may include a reference device that was not seen by one or more of the calibration points during the survey. Second, the current scan may have missed a reference device it should have received a signal from because of noise or because the reference device was busy. Third, new reference devices can be installed after the survey.
As mentioned above, it can be possible that one or more calibration points may have missed a reference device during a survey because the reference device was interfered with by noise or heavy network traffic. However, as noted earlier when discussing the second lemma, the greater the RSSI from a reference device the greater the probability of a packet being received. Conversely, the RSSI of reference devices that were missed should not be very large because if the reference device had a high RSSI, the reference device should have been found during the survey.
For example, if an average RSSI of 8 for a reference device results in loss of signal 10% of the time, then the probability that the reference device will not be seen in N tries is 10−N. If the average RSSI of 5 results in loss of signal 90% of the time, then the probability that the reference device will not be seen in N tries is 0.9−N. For five scans in a survey, the probability that the reference device will not be seen by a calibration point if the average RSSI is 8 is 10−5 or 0.001%. For five (5) scans, the probability that the reference device will not be seen if the RSSI is five (5) is 0.95 or 59%.
Therefore, by setting a threshold as a filter criteria and selecting reference devices that exceed a certain RSSI, the chance that there are reference devices missed in the survey that are seen in the scan can be greatly decreased. The threshold value to choose is a matter of choice and can be based on the number of scans taken during the survey and the relationship between RSSI and the probability of receiving a packet, the sensitivity or other parameters of the wireless receiver and the like.
Another problem encountered in the first method discussed in conjunction with FIG. 2 is that the scan may not see a particular reference device seen by the calibration points because of either heavy network traffic or strong noise at the time of the scan. Because the first method as discussed in conjunction with FIG. 2 searched calibration points that see exactly the same set of reference devices as in the scan set, there can be a chance that the search for the calibration point will be conducted in the wrong cluster. To avoid this, the search “radius” (number of possible calibration points to search) can be increased by searching not the calibration points that have the same set of reference devices that received signals, but calibration points having at least the same reference devices. To avoid expanding the search radius too much, the number of additional reference devices above the amount in the scan set can be set to a fixed number, K, known also as the portion. In one embodiment, K can be set to be 10% of the total number of reference devices to limit the extra number of access points, although any value of K can be chosen. The larger the K value chosen, the larger the potential calculation overhead due to the larger size of the cluster in which to search. The selection of the portion K can, in one embodiment, be partially determined by the wireless channel information and the properties of the particular receiver.
To solve the problem of a scan reading signals from new reference devices that were not part of the survey, any RSSI value from an unknown reference device can be ignored.
FIG. 3 is a flowchart illustrating a second exemplary method for location in accordance with the teachings of the present invention that implements using thresholds and increasing the search radius. In a first step, a scan is performed and a set S is formed, where S: {pr(1), pr(2) . . . pr(M)} (step 302). If the scan contains any information from an unknown reference device that does not belong to AP, the unknown reference device can be considered as a newly installed reference device and the RSSI from the newly found reference device will not be used to calculate Euclidean distance. Next, the reference devices in set S that have a RSSI exceeding a fixed, predetermined threshold, Th1, are extracted to form a new set {APV} (step 304).
Next, for each of the calibration points a set can be formed comprising reference devices having a RSSI that exceed a certain threshold, Th2 (step 306). This set can be labeled {APC1(w)}. A second set can be formed for each calibration points comprising reference devices for which an RSSI is available (step 308). This set can be denoted {APC2(w)}.
Once those sets are formed, calibration points where {APV}⊂{APC2(w)} and where length(APC1(w)−APC1(w)∩APV)≦K, will form a cluster, Gn. Note that length is a function that determines the number of elements in a set. Thus, length(APC1(w)−APC1(w)∩APV) will result in the number of elements that belong to {APC1(w)} and do not belong to {APV}.
Therefore, to form cluster Gn, calibration points in set {APC1(w)} that includes at least the same reference devices in the set {APV} and the calibration points in set {APC1(w)} that have no more than K elements than the number of elements in {APC1(w)∩APV} are chosen.
If the cluster, Gn, is empty, then calibration points are used to calculate a series of Euclidean distances between the calibration points and the scan set. If there are calibration points within cluster Gn, then the Euclidean distance between the calibration point in cluster Gn and the scan set can be used (step 308). The calibration point with the minimum distance from the scan set can be the location estimation. The Euclidean distance is calculated as:




  
    d
    w
  
  =
  
    
      
        ∑
        
          m
          =
          1
        
        M
      
      ⁢
      
        (
        
          
            
              p
              r
            
            ⁡
            
              (
              
                n
                ,
                m
              
              )
            
          
          -
          
            
              {
              
                
                  P
                  rC
                
                ⁡
                
                  (
                  m
                  )
                
              
              )
            
            2
          
        
      
    
  




Alternatively, to reduce computational load, the square of the distance can be calculated and the smallest value chosen as the location estimation.




  
    d
    w
    2
  
  =
  
    
      ∑
      
        m
        =
        1
      
      M
    
    ⁢
    
      (
      
        
          
            p
            r
          
          ⁡
          
            (
            
              n
              ,
              m
            
            )
          
        
        -
        
          
            {
            
              
                P
                rC
              
              ⁡
              
                (
                m
                )
              
            
            )
          
          2
        
      
    
  




In an example of the second method, assume that there are four reference devices and four calibration points distributed throughout an area. The radio map for the calibration points is shown in Table I.














Calibration Points
Pr (w, 1)
Pr (w, 2)
Pr (w, 3)
Pr (w, 4)



















w = 1 (X1, Y1, Z1)
10
NAN
NAN
NAN


w = 2 (X2, Y2, Z2)
9
21
7
5


w = 3 (X3, Y3, Z3)
12
22
7
21


w = 4 (X4, Y4, Z4)
13
21
8
5









Also, assume that the current scan is S: {12, 23, 5, nan} and Th1=10, Th2=20, and K=0.
First, a subset, {APV}, of the current scan is formed using the reference devices which have an RSSI exceeding Th1, which in this example is 10. Thus {APV}={AP1, AP2}.
Next, a subset, {APC1(w)}, for each calibration point is formed using reference devices that have an RSSI exceeding a threshold. In this step the threshold is Th2, which in this example is 20.
{APC1(1)}={ }; {APC1(2)}={AP2}; {APC1(3)}={AP2, AP4}; {APC1(4)}={AP2}
Next, another subset, {APC2(w)}, for each calibration point is formed using reference devices have available signal strength for that calibration point.
{APC2(1)}={AP1}; {APC2(2)}={AP1, AP2, AP3, AP4}; {APC2(3)}={AP1, AP2, AP3, AP4}; {APC2(4)}={AP1, AP2, AP3, AP4}
Next, the calibration points from the subsets {APC1(w)} and {APC2(w)} are used to form a cluster. The cluster will contain calibration points where: {APV}⊂{APC2(w)} and length ({APC1(W)}−{APC1(W)∩{APV})≦K.
Considering {APV⊂APC2(w)}, APC2(1) is eliminated from consideration since APC2(1)} only has AP1 as part of the set.
Next, length(APC1(w)−{APC1(w)}∩{APV}) is determined:

    
    
        For w=2, length ({APC1(2)}−{APC1(2)}∩{APV})=length ({ })=0
        For w=3, length ({APC1(3)}−{APC1(3)}∩{APV})=length ({AP4})=1
        For w=4, length ({APC1(4)}−{APC1(4)}∩{APV}=length ({ })=0
    
    


Since K=0, the third calibration point is eliminated from contention. Therefore, Gn={calibration point 2, calibration point 4}. Now, the Euclidean distance between the set {APV} and the calibration points in Gn will be calculated using:




  
    d
    w
  
  =
  
    
      
        ∑
        
          m
          =
          1
        
        M
      
      ⁢
      
        (
        
          
            
              p
              r
            
            ⁡
            
              (
              
                n
                ,
                m
              
              )
            
          
          -
          
            
              {
              
                
                  P
                  rC
                
                ⁡
                
                  (
                  m
                  )
                
              
              )
            
            2
          
        
      
    
  




For w=2:

d2=√{square root over ((12−9)2+(23−21)2+(5−7)2)}{square root over ((12−9)2+(23−21)2+(5−7)2)}{square root over ((12−9)2+(23−21)2+(5−7)2)}

d2=4.12

And for w=4

d4=√{square root over ((12−12)2+(23−22)2+(5−7)2)}{square root over ((12−12)2+(23−22)2+(5−7)2)}{square root over ((12−12)2+(23−22)2+(5−7)2)}

d4=3.74

The smallest distance is d4; therefore d4 is the location estimate. Alternatively, the square of the distance could be calculated.
While at least one exemplary embodiment has been presented in the foregoing detailed description, it should be appreciated that a vast number of variations exist. It should also be appreciated that the exemplary embodiment or exemplary embodiments are only examples, and are not intended to limit the scope, applicability, or configuration of the invention in any way. Rather, the foregoing detailed description will provide those skilled in the art with a convenient road map for implementing the exemplary embodiment or exemplary embodiments. It should be understood that various changes can be made in the function and arrangement of elements without departing from the scope of the invention as set forth in the appended claims and the legal equivalents thereof.

## Claims

1. A method for determining a location estimation of a wireless device in a wireless location system having a plurality of reference devices, the method comprising:
receiving signals from one or more of the plurality of reference devices at the wireless device;
forming a first set comprising a plurality of elements, each element of the plurality of elements comprising a signal strength measurements of signals received at the wireless device from each of the plurality of reference devices from which the wireless device receives signals and an indicator of no signal strength measured for each of the plurality of reference devices from which the wireless device does not receive signals;
forming a scan subset comprising each of the plurality of reference devices in the first set that has an associated signal strength measurement;
forming a cluster set comprising calibration points selected from a plurality of calibration points based on the scan subset, wherein the plurality of calibration points comprise a set of received signal strength measurements from signals received at the calibration points from each of the plurality of reference devices; and
calculating a calculated distance between the first set and each of the calibration points in the cluster set, when the cluster set has at least one element;
calculating the calculated distance between the first set and each of the plurality of calibration points, when the cluster set has no elements; and
selecting the calibration point having the smallest calculated distances as the location estimation.

2. The method of claim 1 wherein the step of forming a cluster set further comprises forming a cluster set comprising calibration points selected from the plurality of calibration points that contain exactly the same reference devices as the scan subset.

3. The method of claim 1 wherein the step of forming a first set further comprises detecting newly installed reference devices by identifying signals from reference devices that did not appear in any calibration point within the plurality of calibration points.

4. The methods of claim 1 wherein the step of forming a scan subset further comprises forming a scan subset comprising each of the plurality of reference devices in the first set that has an associated signal strength measurement above a fixed threshold.

5. A system for estimating a location of a wireless device in a wireless network comprising:
a plurality of reference devices distributed throughout the wireless network, each of the plurality of reference devices operable to send a signal;
a plurality of calibration points distributed throughout the wireless network, the calibration points comprising a set of received signal strength measurements from signals received from each of the plurality of reference devices; and
the wireless device located in the network, the wireless device configured to:
form a scan set comprising a received signal strength indicator (RSSI) value from each of the plurality of reference devices from which the wireless device receives signals from and an indicator of no signal received from each of the plurality of reference devices that the wireless device does not receive signals as elements of the scan set;
form a scan subset of the scan set comprising reference devices in the scan set that have associated signal strength;
form a cluster set comprising calibration points selected from the plurality of calibration points based on the scan subset;
determine an estimated location by selecting the smallest distance between the scan subset and each of the calibration points in the cluster set, when the cluster set has at least one element; and
determine an estimated location by selecting the smallest distance between the scan subset and each of the plurality of calibration points, when the cluster set has no elements.

6. The system of claim 5 wherein the cluster set comprises calibration points selected from the plurality of calibration points that contain the exact same reference devices as the scan subset.

7. The system of claim 6 wherein the scan set contains reference devices having an RSSI above a fixed threshold.

8. The system of claim 5 wherein any reference devices in the scan set that are not in the plurality of calibration points are considered newly installed reference devices.

9. The system of claim 5 wherein the wireless device is a radio frequency identification tag.

10. The system of claim 5 wherein the wireless device is a handheld computer unit.

11. The system of claim 5 wherein the wireless device is a computer with a wireless network transceiver.

12. The system of claim 5 wherein the distance is calculated as the square of the distance between the scan subset and each of the calibration points in the system.

13. The system of claim 5 wherein the reference device comprises an access point of a wireless network.

14. A wireless device for use in a wireless network, the wireless device comprising:
a memory configured to store a plurality of calibration points, each of the plurality of calibration points including a set of received signal strength indicator (RSSI) values for signals received from each of a plurality of reference devices;
a receiver configured to receive signals from one or more of the plurality of reference devices; and
a processor coupled to the receiver and the memory, the processor operable to:
determine RSSI values for each of the signals received from the one or more of the plurality of reference devices;
form a scan set comprising the RSSI value from each of the plurality of reference devices that the receiver received signals from and further comprising an indicator of no signal received from each of the plurality of reference devices that the receiver did not receive signals;
form a scan subset comprising reference devices in the scan set that have an associated RSSI value;
form a cluster set comprising calibration points selected from the plurality of calibration points based on the scan set;
determine an estimated location by selecting the smallest distance between the scan subset and each of the calibration points of the cluster set, when the cluster set has at least one element; and
determine an estimated location by selecting the smallest distance between the scan subset and each of the plurality of calibration points, when the cluster set has no elements.

15. The device of claim 14 wherein the wireless device is part of a radio frequency identification tag.

16. The device of claim 14 wherein the wireless device is a wireless data entry terminal.

17. A method for estimating a location of a wireless device in a wireless location system having a plurality of reference devices and a plurality of calibration points, each of the calibration points associated with a set of received signal strength measurements received at the calibration point from each of the plurality of reference devices from which a signal can be receive, the method comprising:
forming a first set comprising a plurality of elements, each element of the plurality of elements corresponding to signal strength measurements of signals received at the wireless device from each of the plurality of reference devices from which the wireless device receives signals and an indicator of no signal strength measured for any of the plurality of reference devices from which the wireless device does not receive signals:
forming a scan subset comprising each of the plurality of reference devices in the first set that has an associated signal strength measurement above a first fixed threshold;
forming a cluster set by:
forming one or more first calibration subset for each calibration point comprising reference devices that have an associated signal strength above a second fixed threshold;
forming one or more second calibration subset for each calibration point comprising reference devices having an available received signal strength indicator;
selecting calibration points that include reference devices that are in the scan subset and contained in the one or more second calibration subset and where a difference between the number of reference devices that belong to the one or more first calibration subset that do not belong to the intersection of the reference devices in the scan subset and the reference devices in the one or more first calibration subset is equal to or does not exceed a fixed portion of the total reference devices; and

calculating a calculated distance between the first set and each of the calibration points in the cluster set when the cluster set has at least one element;
calculating the calculated distance between the first set and each of the plurality of calibration points when the cluster set has no elements; and
selecting the smallest of the calculated distances as a location estimate.

18. A wireless device for use in a wireless network, the wireless device comprising:
a memory configured to store a plurality of calibration points, each of the plurality of calibration points including a set of received signal strength indicator (RSSI) values for signals received from each of a plurality of reference devices;
a receiver configured to receive signals from one or more of the plurality of reference devices; and
a processor coupled to the receiver and the memory, the processor operable to:
determine a received signal strength indicator (RSSI) value for each of the signals received from the one or more of the plurality of reference devices;
form a scan set comprising the RSSI value from each of the plurality of reference devices that the receiver received a signal from and further comprising an indicator of no signal received from each of the plurality of reference devices that the receiver did not receive a signal;
form a scan subset comprising reference devices in the scan set that have an associated RSSI value above a first fixed threshold;
form one or more first calibration subset for each calibration point comprising reference devices that has an associated signal strength above a second fixed threshold; and
form one or more second calibration subset for each calibration subset comprising reference devices having an available received signal strength indicator (RSSI); and
form a cluster set by selecting calibration points that include reference devices that are in the scan subset and contained in the one or more second calibration subset and where a difference between the number of reference devices that belong to the one or more first calibration subset that do not belong to the intersection of the reference devices in the scan subset and the reference devices in the one or more first calibration subset is less than or equal to a fixed portion of the total reference devices; and
determine an estimated location by selecting the smallest distance between the scan subset and each of the calibration points of the cluster set when the cluster set has at least one element; and
determine an estimated location by selecting the smallest distance between the scan subset and each of the plurality of calibration points when the cluster set has no elements.