# Patent

## Title
System and method for three-dimensional geolocation of emitters based on energy measurements

## Abstract

According to an embodiment of the present invention, a three-dimensional (3-D) energy-based emitter geolocation technique determines the geolocation of a radio frequency (RF) emitter based on energy or received signal strength (RSS) of transmitted signals. The technique may be employed with small unmanned air vehicles (UAV), and obtains reliable geolocation estimates of radio frequency (RF) emitters of interest.

## Background
BACKGROUND
1. Technical Field
The present invention embodiments pertain to determining locations of emitters. In particular, the present invention embodiments pertain to determining locations of radio frequency (RF) emitters in a three-dimensional space based on energy or received signal strength (RSS) measurements of the emitters at various locations.
2. Discussion of Related Art
Conventional techniques for geolocation of a radio frequency (RF) emitter are commonly based on measurements of a received signal strength (RSS) of signals transmitted from the emitter. Accordingly, this RSS-based technique may also be referred to as an energy-based technique. The received signal strength (RSS) is usually integrated over the duration of the transmitted signal in order to obtain a signal energy measurement and enhance signal to noise ratio. Since the transmitted radio frequency (RF) signal attenuates during propagation through space, the received signal strength (RSS) of the signal may be used to estimate the distance from the receiver to the emitter. However, this technique may not be very accurate due to multipath fading and shadowing effects that may significantly distort the expected received signal strength (RSS).
In order to improve the accuracy, the received signal strength (RSS) or energy-based geolocation technique may employ a radio frequency (RF) propagation map of the environment. The map is basically a database with information created from known terrain data and learned from observed energy measurements at different combinations of emitter and receiver locations. The improved geolocation technique determines the best fit in the energy space to potential emitter locations based on the radio frequency (RF) propagation map. However, this improved technique requires a large number of real-time measurements and/or terrain modeling. Thus, this technique can only be used in applications in which the radio frequency (RF) propagation map is available, and sufficient computing capacity exists to process the large amount of data.
SUMMARY
An embodiment of the present invention pertains to a three-dimensional (3-D) energy-based emitter geolocation technique that determines the geolocation of a radio frequency (RF) emitter based on energy or received signal strength (RSS) of transmitted signals. The technique may be employed with small unmanned air vehicles (UAV), and obtains reliable geolocation estimates of radio frequency (RF) emitters of interest.
Present invention embodiments provide several advantages. For example, the technique of present invention embodiments provides the simplicity and the performance robustness required by a low-cost, compact system. The use of a small unmanned air vehicle (UAV) provides a cost-effective manner to reliably measure received signal strength (RSS) data generated from the radio frequency (RF) emitter of interest. Further, a Line-of-Sight (LOS) condition offered by the unmanned air vehicle (UAV) minimizes errors from path loss modeling and the effects from signal fading and shadowing. Each RSS measurement that feeds the geolocation algorithm is the measurement with the maximum signal to noise ratio selected from a block of consecutive RSS data (which is referred to herein as the maximum signal to noise ratio (MSNR) rule). The block size can be determined using the assessment of the spaced-frequency spaced-time correlation function of the propagation channel. The MSNR rule further enhances the quality of the RSS measurements used by the geolocation algorithm. In addition, the combination of the technique with the use of an unmanned air vehicle (UAV) enables an overall system to be small, compact, flexible, reliable, and of low-cost.
The above and still further features and advantages of present invention embodiments will become apparent upon consideration of the following detailed description of example embodiments thereof, particularly when taken in conjunction with the accompanying drawings wherein like reference numerals in the various figures are utilized to designate like components.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a diagrammatic illustration of an example environment for determining geolocation of a radio frequency (RF) emitter according to an embodiment of the present invention.
FIG. 2 is a block diagram of a system for determining geolocation of a radio frequency (RF) emitter according to an embodiment of the present invention.
FIG. 3 is a procedural flow chart illustrating a manner in which to determine geolocation of a radio frequency (RF) emitter according to an embodiment of the present invention.
FIG. 4 is a graphical representation of simulation results for an embodiment of the present invention illustrating the relationship between geolocation error and the quantity of locations for received signal strength (RSS) measurements.
FIG. 5 illustrates representations of simulation results for an embodiment of the present invention showing determined geolocations for a radio frequency (RF) emitter under various conditions.



DETAILED DESCRIPTION OF EXAMPLE EMBODIMENTS
Embodiments of the present invention pertain to a three-dimensional (3-D) energy-based geolocation technique that obtains reliable geolocation estimates of a radio frequency (RF) emitter based on energy or received signal strength (RSS) of transmitted signals. The geolocation of a radio frequency (RF) emitter is a critical need for many applications. The technique of present invention embodiments may be employed with unmanned air vehicles (UAV) that are usually small, utilized for low altitudes, and employ typical guidance technologies for operation (e.g., following pre-planned or manually provided paths or waypoints). These types of vehicles are well suited for enabling three-dimensional (3-D) geolocation of radio frequency (RF) emitters of interest.
An example environment for determining the geolocation of a radio frequency (RF) emitter in a three-dimensional space is illustrated in FIG. 1. Specifically, the environment includes a radio frequency (RF) emitter 120, and a mobile sensor 100 (e.g., an unmanned air vehicle (UAV) or other platform with a radio frequency (RF) sensor, etc.). The mobile sensor travels along a pre-planned path 110 (e.g., a pre-planned flight path in the case of an unmanned air vehicle (UAV)). Mobile sensor 100 includes an antenna 130 that receives signals from radio frequency (RF) emitter 120 in order to measure the strength of those signals as described below. The radio frequency (RF) emitter and mobile sensor are located within a three-dimensional space of the environment (e.g., defined by X, Y, and Z axes as illustrated in FIG. 1). Locations within the three-dimensional space may be represented by coordinates that indicate a position along each of the respective X, Y, and Z axes. By way of example, radio frequency (RF) emitter 120 is positioned at an unknown location (x, y, z) within the three-dimensional space, while mobile sensor 100 receives signals transmitted from the radio frequency (RF) emitter at known locations along path 110 within the three-dimensional space (e.g., locations (x0, y0, z0), (x1, y1, z1), and (x2, y2, z2) as viewed in FIG. 1). The Z axis represents the height or altitude, and indicates the offset between the mobile sensor and pre-planned path 110 (e.g., distances z0, z1, z2 as viewed in FIG. 1).
Mobile sensor 100 measures at selected locations (e.g., (x0, y0, z0), (x1, y1, z1), and (x2, y2, z2) as viewed in FIG. 1) the received signal strength (RSS) (e.g., p0, p1, p2 as viewed in FIG. 1) of radio frequency (RF) signals emitted by emitter 120. The received signal strength (RSS) at each location is proportional to the distance (e.g., r0, r1, r2 as viewed in FIG. 1) between that location and radio frequency (RF) emitter 120. The received signal strength (RSS) measurement can be viewed as a special case of the signal energy in which only a single signal sample is used for the measurement at each location. Each RSS measurement that feeds the geolocation algorithm is the measurement with the maximum signal to noise ratio selected from a block of consecutive RSS data (which is referred to herein as the maximum signal to noise ratio (MSNR) rule). The block size can be determined using the assessment of the spaced-frequency spaced-time correlation function of the propagation channel.
Once mobile sensor 100 collects the received signal strength (RSS) measurements, the geolocation estimate of radio frequency (RF) emitter 120 is determined based on those measurements as described below. The received signal strength (RSS) measurements may be collected by using an unmanned air vehicle (UAV) or other platform along a flight or other pre-planned path, or by using plural unmanned air vehicles (UAV) or other platforms each collecting a measurement at one or more locations along that path. In other words, measurements from plural locations may be ascertained via a single platform traveling to different locations, or via plural platforms each positioned at different locations and networking or otherwise sharing the collected data for the geolocation determination. Since measurement errors exist due to path loss modeling, signal fading, shadowing effects, noise/interference, antenna pattern effects, time-varying channel and transmit power effects, and implementation errors, a Least Mean Square (LMS) technique is preferably employed to determine the location of radio frequency (RF) emitter 120 as described below. Although FIG. 1, by way of example only, indicates measurements at certain locations (e.g., (x0, y0, z0), (x1, y1, z1), and (x2, y2, z2) as viewed in FIG. 1), any quantity of received signal strength (RSS) measurements (e.g., pi, where i=0 to N) may be collected at any corresponding locations ((xi, yi, zi), where i=0 to N) within the three-dimensional space.
Present invention embodiments resolve the location of radio frequency (RF) emitter 120 by estimating the energy or received signal strength (RSS) of signals emitted from emitter 120 via the received signal strength (RSS) measurements ascertained from plural locations (e.g., p0, p1, p2 measured at locations (x0, y0, z0), (x1, y1, z1), and (x2, y2, z2) as viewed in FIG. 1) along path 110. The received signal strength (RSS) measurements are each proportional to the distance between the location of that measurement and radio frequency (RF) emitter 120 (e.g., r0, r1, r2 as viewed in FIG. 1) as described above. The measurements are utilized in a set of simultaneous equations to determine the location of the radio frequency (RF) emitter within the three-dimensional space as described below.
An example system 200 for determining the geolocation of a radio frequency (RF) emitter according to an embodiment of the present invention is illustrated in FIG. 2. Initially, system 200 preferably resides on mobile sensor 100 (FIG. 1) to measure the received signal strength (RSS) and determine the geolocation of the radio frequency (RF) emitter. However, the processing and one or more other portions of system 200 may be remote from the mobile sensor and receive the received signal strength (RSS) measurements for the geolocation determination. In particular, system 200 includes antenna 130, a receiver 210, and a processing device 230. Antenna 130 is preferably implemented by an omni-directional antenna, and directs received signals into receiver 210. The antenna may be implemented by any conventional or other antenna configurable to receive the signals emitted from radio frequency (RF) emitter 120.
Receiver 210 includes a radiometer or energy detector 220 that provides an energy measure (e.g., received signal strength (RSS)) of the signals received from antenna 130. The receiver may be implemented by any conventional or other receiving device capable of receiving the emitted radio frequency (RF) signals, while the radiometer may be implemented by any conventional or other device to measure the energy or received signal strength (RSS) of a received signal. Based on the MSNR rule, the selected, received signal strength (RSS) measurements are provided to processing device 230 to determine the geolocation of radio frequency (RF) emitter 120 as described below.
Processing device 230 may include a processor 250, a memory 260, and an interface unit 270. Processor 250 determines the geolocation of radio frequency (RF) emitter 120 based on the measurements received from receiver 210 and provides corresponding geolocation data 240. The processor includes one or more location modules to determine the location of radio frequency (RF) emitter 120 from a set of simultaneous equations incorporating a Least Mean Square (LMS) technique as described below. The processor may be implemented by any conventional or other computer or processing unit (e.g., a microprocessor, a microcontroller, systems on a chip (SOCs), fixed or programmable logic, etc.), where the one or more location modules may be implemented by any combination of any quantity of software and/or hardware modules or units. Memory 260 may be included within or external of processor 250, and may be implemented by any conventional or other memory unit with any type of memory (e.g., random access memory (RAM), read only memory (ROM), etc.). The memory may store the one or more location modules for execution by processor 250, and data for performing the geolocation technique of present invention embodiments. Interface unit 270 enables communication between system 200 and other devices or systems, and may be implemented by any conventional or other communications device (e.g., wireless communications device, etc.).
The manner in which processor 250 (e.g., via one or more location modules) determines the geolocation of a radio frequency (RF) emitter based on received signal strength (RSS) at various locations is illustrated in FIGS. 1 and 3. Initially, one or more mobile sensors 100 measure received signal strength (RSS) of signals emitted from radio frequency (RF) emitter 120 at one or more locations (e.g., a quantity of locations from 0 through N as described below) along path 110 at step 300. A set of simultaneous equations to determine the geolocation of the radio frequency (RF) emitter based on the received signal strength (RSS) measurements are determined, and converted into matrix form. In particular, the location of radio frequency (RF) emitter 120 within the three-dimensional space may be represented by the coordinates (x, y, z), while the position of mobile sensor 100 ascertaining a measurement at an ith location along path 110 may be represented by the coordinates (xi, yi, zi). The distance (e.g., ri, as viewed in FIG. 1) in the three-dimensional space between the location of the radio (RF) frequency emitter (e.g., (x, y, z) as viewed in FIG. 1) and the ith measuring location (e.g., (xi, yi, zi)), may be expressed as the following:

ri2=(x−xi)2+(y−yi)2+(z−zi)2; for i=0 to N.  (Equation 1)

The distance (e.g., di, for i=0 to N) between a reference origin in the three-dimensional space (e.g., (0, 0, 0)) and a location of mobile sensor 100 (e.g., (xi, yi, zi)) may be expressed as the following:

di2=xi2+yi2+zi2; for i=0 to N.  (Equation 2)

The difference of the square of the distances (e.g., ri2−r02) for the ith measuring location (e.g., (xi, yi, zi)) and an arbitrary reference location of mobile sensor 100 (e.g., (x0, y0, z0)) may be expressed (based on Equations 1 and 2) as the following:

ri2−r02=di2−d02−2x(xi−x0)−2y(yi−y0)−2z(zi−z0), for i=1 to N,  (Equation 3)

where this equation (Equation 3) may be equivalently expressed as the following equation:





  
    
      
        
          
            
              [
              
                
                  
                    r
                    i
                    2
                  
                  
                    r
                    0
                    2
                  
                
                -
                1
              
              ]
            
            ⁢
            
              r
              0
              2
            
          
          +
          
            2
            ⁢
            
              x
              ⁡
              
                (
                
                  
                    x
                    i
                  
                  -
                  
                    x
                    0
                  
                
                )
              
            
          
          +
          
            2
            ⁢
            
              y
              ⁡
              
                (
                
                  
                    y
                    i
                  
                  -
                  
                    y
                    0
                  
                
                )
              
            
          
          +
          
            2
            ⁢
            
              z
              ⁡
              
                (
                
                  
                    z
                    i
                  
                  -
                  
                    z
                    0
                  
                
                )
              
            
          
        
        =
        
          
            d
            i
            2
          
          -
          
            
              d
              0
              2
            
            .
          
        
      
    
    
      
        (
        
          Equation
          ⁢
          
            
          
          ⁢
          4
        
        )
      
    
  




The above equation (Equation 4) may be simplified by employing a parameter, βi, which corresponds to the ith measuring location, and may be expressed as follows:





  
    
      
        
          
            β
            i
          
          =
          
            [
            
              
                
                  r
                  i
                  2
                
                
                  r
                  0
                  2
                
              
              -
              1
            
            ]
          
        
        ,
        
          
            for
            ⁢
            
              
            
            ⁢
            i
          
          =
          
            1
            ⁢
            
              
            
            ⁢
            to
            ⁢
            
              
            
            ⁢
            
              N
              .
            
          
        
      
    
    
      
        (
        
          Equation
          ⁢
          
            
          
          ⁢
          5
        
        )
      
    
  




In addition, the terms of the above equation (Equation 4) may be converted to matrix form and employ the parameter, βi (from Equation 5). The equation terms may be expressed by matrices P (e.g., representing terms on the left side of the equal sign in Equation 4) and R (e.g., representing terms on the right side of the equal sign in Equation 4) as follows:





  
    P
    =
    
      [
      
        
          
            
              β
              1
            
          
          
            
              2
              ⁢
              
                (
                
                  
                    x
                    1
                  
                  -
                  
                    x
                    0
                  
                
                )
              
            
          
          
            
              2
              ⁢
              
                (
                
                  
                    y
                    1
                  
                  -
                  
                    y
                    0
                  
                
                )
              
            
          
          
            
              2
              ⁢
              
                (
                
                  
                    z
                    1
                  
                  -
                  
                    z
                    0
                  
                
                )
              
            
          
        
        
          
            
              β
              2
            
          
          
            
              2
              ⁢
              
                (
                
                  
                    x
                    2
                  
                  -
                  
                    x
                    0
                  
                
                )
              
            
          
          
            
              2
              ⁢
              
                (
                
                  
                    y
                    2
                  
                  -
                  
                    y
                    0
                  
                
                )
              
            
          
          
            
              2
              ⁢
              
                (
                
                  
                    z
                    2
                  
                  -
                  
                    z
                    0
                  
                
                )
              
            
          
        
        
          
            ...
          
          
            ...
          
          
            ...
          
          
            ...
          
        
        
          
            ...
          
          
            ...
          
          
            ...
          
          
            ...
          
        
        
          
            
              β
              N
            
          
          
            
              2
              ⁢
              
                (
                
                  
                    x
                    N
                  
                  -
                  
                    x
                    0
                  
                
                )
              
            
          
          
            
              2
              ⁢
              
                (
                
                  
                    y
                    N
                  
                  -
                  
                    y
                    0
                  
                
                )
              
            
          
          
            
              2
              ⁢
              
                (
                
                  
                    z
                    N
                  
                  -
                  
                    z
                    0
                  
                
                )
              
            
          
        
      
      ]
    
  
  ,
  
    

  
  ⁢
  
    R
    =
    
      [
      
        
          
            
              
                d
                1
                2
              
              -
              
                d
                0
                2
              
            
          
        
        
          
            
              
                d
                2
                2
              
              -
              
                d
                0
                2
              
            
          
        
        
          
            ...
          
        
        
          
            ...
          
        
        
          
            
              
                d
                N
                2
              
              -
              
                d
                0
                2
              
            
          
        
      
      ]
    
  




The overall equation (Equation 4) may be represented by the following matrix equation:




  
    
      
        
          P
          ·
          
            [
            
              
                
                  
                    r
                    0
                    2
                  
                
              
              
                
                  x
                
              
              
                
                  y
                
              
              
                
                  z
                
              
            
            ]
          
        
        =
        R
      
    
    
      
        (
        
          Equation
          ⁢
          
            
          
          ⁢
          6
        
        )
      
    
  




The terms xi, yi, zi, (for i=0 to N) within matrix P represent the known positions or coordinates in the three-dimensional space where mobile sensor 100 ascertains the received signal strength (RSS) measurements, while the terms r02, x, y, and z in the solution matrix are unknown and to be solved by the above equation (Equation 6). The determined values for x, y, and z represent the coordinates (or location) of radio frequency (RF) emitter 120 within the three-dimensional space, while the determined value for r02 represents the square of the distance between radio frequency (RF) emitter 120 and the known reference location (e.g., at coordinates x0, y0, and z0 within the three-dimensional space) of mobile sensor 100.
The values for the unknown variables (e.g., r02, x, y, and z) indicating the location of radio frequency (RF) emitter 120 may be determined by solving for these variables in Equation 6, thereby providing the following expression:




  
    
      
        
          
            [
            
              
                
                  
                    r
                    0
                    2
                  
                
              
              
                
                  x
                
              
              
                
                  y
                
              
              
                
                  z
                
              
            
            ]
          
          =
          
            
              
                (
                
                  
                    P
                    T
                  
                  ⁢
                  P
                
                )
              
              
                -
                1
              
            
            ⁢
            
              P
              T
            
            ⁢
            R
          
        
        ,
      
    
    
      
        (
        
          Equation
          ⁢
          
            
          
          ⁢
          7
        
        )
      
    
  




where PT represents the transpose of matrix P, and (PTP)−1 represents the inverse of the product of matrix P and the transpose of matrix P.

In order to determine the unknown variables (e.g., r02, x, y, and z) indicating the location of radio frequency (RF) emitter 120 in the above equation (Equation 7), the parameter, βi, of matrix P may be estimated based on the measurements of received signal strength (RSS) obtained by mobile sensor 100. Considering the line of sight (LOS) propagation loss between mobile sensor 100 (e.g., unmanned air vehicle (UAV)) and radio frequency (RF) emitter 120, the received signal power, pi at the ith location along path 110 is inversely proportional to the square law of the distance, ri, between the mobile sensor (e.g., unmanned air vehicle (UAV)) and the radio frequency (RF) emitter. Assuming the power of radio frequency (RF) emitter 120 remains constant during the measurements of received signal strength (RSS) along path 110, the parameter, βi, may be estimated based on the received signal strength (RSS) or power measurements as follows:




  
    
      
        
          
            β
            i
          
          =
          
            
              [
              
                
                  
                    r
                    i
                    2
                  
                  
                    r
                    0
                    2
                  
                
                -
                1
              
              ]
            
            ≅
            
              [
              
                
                  
                    p
                    0
                  
                  
                    p
                    i
                  
                
                -
                1
              
              ]
            
          
        
        ,
        
          
            for
            ⁢
            
              
            
            ⁢
            i
          
          =
          
            1
            ⁢
            
              
            
            ⁢
            to
            ⁢
            
              
            
            ⁢
            
              N
              .
            
          
        
      
    
    
      
        (
        
          Equation
          ⁢
          
            
          
          ⁢
          8
        
        )
      
    
  




At least four independent equations (or at least four rows of matrices P and R) are required to determine the four unknown variables (e.g., r02, x, y, and z) and, hence, the location of radio frequency (RF) emitter 120. However, measurements from at least five locations are required to provide estimates for the parameter, βi (e.g., a reference measurement for p0, and a measurement for each pi, for i=1 to 4).
The estimates for the parameter, βi (for i=1 to N), and the various terms that can be derived from the known measuring locations of mobile sensor 100 (e.g., xi, yi, zi (for i=0 to N); di2 (for i=0 to N), etc.) are applied to matrices P and R at step 302. The applied values within matrices P and R are utilized in Equation 7 to determine the values for the unknown variables (e.g., r02, x, y, and z) in the solution matrix at step 304. Since there are path loss model errors, signal fading and/or shadowing effects, noise, interference, and implementation errors that impact the measurement, the above determination (Equations 1-7) is formulated to provide a Least Mean Square (LMS) solution for the variables in the solution matrix.
The determined Least Mean Square (LMS) values for x, y, and z within the solution matrix (derived from Equation 7) represent the coordinates of radio frequency (RF) emitter 120 within the three-dimensional space, and are utilized to provide the Least Mean Square (LMS) location of the radio frequency (RF) emitter within that space at step 306. The determined location of the radio frequency (RF) emitter may be used for various applications at step 308. For example, the location information may be processed by processor 230 or forwarded to another system via interface unit 270. The location information may be processed to direct or control a vehicle or other platform to an emitter at a location of interest (e.g., to provide assistance at that location, etc.). Further, the location information may be utilized to generate an image of the area and indicate the emitter locations.
In order to assess the accuracy of the Least Mean Square (LMS) solution, an error vector, ξ, representing the error for the Least Mean Square (LMS) solution may be expressed by the following relationship:




  
    
      
        ξ
        =
        
          
            P
            ·
            
              [
              
                
                  
                    
                      r
                      0
                      2
                    
                  
                
                
                  
                    x
                  
                
                
                  
                    y
                  
                
                
                  
                    z
                  
                
              
              ]
            
          
          -
          R
        
      
    
    
      
        (
        
          Equation
          ⁢
          
            
          
          ⁢
          9
        
        )
      
    
  




The values for the solved variables (e.g., r02, x, y, and z) are utilized to determine the error vector. If a vector, Xo, provided Least Mean Square (LMS) values for the unknown variables (e.g., r02, x, y, and z), the Least Mean Square (LMS) error vector, ξo, may be expressed as follows:

ξo=P·Xo−R  (Equation 10)

The transpose of the P matrix must be orthogonal to the Least Mean Square (LMS) error vector, ξo, thereby providing the following expression:

PTξo=0  (Equation 11)

Accordingly, the Least Mean Square (LMS) solution vector, Xo, must satisfy the orthogonal property, which may be expressed as follows based on the above equations (Equations 10 and 11):

PTξo=PTP·Xo−PTR=0  (Equation 12)

Therefore, the Least Mean Square (LMS) solution vector, Xo, may expressed as follows:

Xo=(PTP)−1PTR  (Equation 13)

Substituting the Least Mean Square (LMS) solution vector, Xo, into the error equation (Equation 10) gives the following expression for the error vector (e.g., indicating the error of the Least Mean Square (LMS) solution (e.g., values for variables r02, x, y, and z from Equation 7) for each of the measuring locations of mobile sensor 100):

ξo=[P(PTP)−1PT−I]R,  (Equation 14)

where I is the identity matrix, (PTP)−1PT is a projection matrix that projects the matrix R onto the column space of matrix P, and [(PTP)−1PTI] is a projection matrix that projects the matrix R onto the orthogonal complement of the column space of matrix P. Geometrically, the projection of the matrix R onto the orthogonal complement of the column space of matrix P assures the orthogonally between the error vector and the column space of matrix P that provides the Least Mean Square (LMS) solution.

The geolocation technique of a present invention embodiment employing small unmanned air vehicles (UAV) has been modeled and simulated using Matlab tools available from The Mathworks, Inc. of Natick, Mass. A graphical illustration of the simulation results providing the relationship between geolocation error and the quantity of locations for received signal strength (RSS) measurements is illustrated in FIG. 4. In the simulation, the locations of an unmanned air vehicle (UAV) were randomly generated, and the following conditions were assumed: the signal to noise ratio (SNR) was 10 dB; the emitter power remained constant during the measurements; the variance of Lognormal shadow varied from 0 to 5 dB; and the path loss follows the line of sight (LOS) law.
As viewed in FIG. 4, a Root Mean Square (RMS) error of the geolocation estimates (e.g., derived from the error vector of Equation 14) for a radio frequency (RF) emitter converges to a robust level with eight or more measuring locations. The geolocation estimates perform well for the variance of Lognormal shadowing less than 2 dB, which is achievable with a line of sight (LOS) path loss condition offered by the unmanned air vehicle (UAV).
Two samples of results of a Monte Carlo simulation conducted for various operating scenarios and randomized parameter variations are illustrated in FIG. 5. The simulations provide the geolocation characteristics for eight measuring locations (e.g., indicators 10 as viewed in FIG. 5) and one emitter location (e.g., indicator 20 as viewed in FIG. 5). The simulation was performed for 100 randomized paths with a signal to noise ratio (SNR) of 10 dB, and a variance of Lognormal Shadow of 0 or 2 dB. The simulation shows that the geolocation technique of present invention embodiments provides an unbiased three-dimensional geolocation estimate (e.g., where the geolocation estimates are shown by indicators 30, and the averaged estimate is shown by indicator 40 as viewed in FIG. 5). The simulation results indicate that the present invention geolocation technique is compatible with small unmanned air vehicles (UAV), and provides geolocation estimates of radio frequency (RF) emitters with enhanced reliability.
It will be appreciated that the embodiments described above and illustrated in the drawings represent only a few of the many ways of implementing a system and method for three-dimensional geolocation of emitters based on energy or received signal strength (RSS) measurements.
The environment of the present invention embodiments may include any quantity of mobile sensors, and emitters. The emitters may be implemented by any quantity of any conventional or other devices emitting radio frequency (RF) or any other suitable energy signals (e.g., energy signals in any suitable bands (e.g., infrared, microwave, optical, etc.)). The emitters may be located at any quantity of any desired locations within the three-dimensional space of the environment. The mobile sensors may be implemented by any quantity of any conventional or other mobile or stationary vehicle or platform (e.g., unmanned air vehicle (UAV), air vehicle, ground vehicle, platform or structure mounted at a location or on a vehicle, etc.), and may include any quantity of any conventional or other sensing device (e.g., RF or other sensor, etc.). The mobile sensors may each measure any desired characteristics of emitted signals at any one or more locations within the environment.
The pre-planned path may traverse any desired locations within the environment, where any quantity of measurements may be obtained during traversal of the path. Further, measurements may be obtained at any locations residing within a specified offset or range from the pre-planned path. Alternatively, the path may be determined in random fashion.
The antenna may be implemented by any conventional or other antenna (e.g., omni-directional, directional, etc.) configurable to receive the signals emitted from the one or more emitters. The receiver may be implemented by any conventional or other receiving device capable of receiving the emitted radio frequency (RF) or other energy signals. The radiometer may be implemented by any conventional or other device to measure the energy or received signal strength (RSS) or other characteristics of a received signal. The radiometer may be included within or separate from the receiver.
The processor may be implemented by any quantity of any conventional or other computer systems or processing units (e.g., a microprocessor, a microcontroller, systems on a chip (SOCs), fixed or programmable logic, etc.), and may include any commercially available or custom software (e.g., communications software, location modules, etc.).
It is to be understood that the software (e.g., location modules, etc.) for the processor of the present invention embodiments may be implemented in any desired computer language and could be developed by one of ordinary skill in the computer arts based on the functional descriptions contained in the specification and flow charts illustrated in the drawings. Further, any references herein of software performing various functions generally refer to computer systems or processors performing those functions under software control. The processor of the present invention embodiments may alternatively be implemented by any type of hardware and/or other processing circuitry. The various functions of the processor may be distributed in any manner among any quantity of software modules or units, processing or computer systems and/or circuitry, where the computer or processing systems may be disposed locally or remotely of each other and communicate via any suitable communications medium (e.g., LAN, WAN, Intranet, Internet, hardwire, modem connection, wireless, etc.). For example, the functions of the present invention embodiments may be distributed in any manner among the processor, receiver, and/or external devices. The software and/or algorithms described above and illustrated in the flow charts may be modified in any manner that accomplishes the functions described herein. In addition, the functions in the flow charts or description may be performed in any order that accomplishes a desired operation.
The software of the present invention embodiments (e.g., location modules, etc.) may be available on a program product apparatus or device including a recordable or computer usable medium (e.g., magnetic or optical mediums, magneto-optic mediums, floppy diskettes, CD-ROM, DVD, memory devices, etc.) for use on stand-alone systems or systems connected by a network or other communications medium, and/or may be downloaded (e.g., in the form of carrier waves, packets, etc.) to systems via a network or other communications medium. Further, the tangible recordable or computer usable medium may be encoded with instructions or logic to perform the functions described herein (e.g., embedded logic such as an application specific integrated circuit (ASIC), digital signal processor (DSP) instructions, software that is executed by a processor, etc.).
The memory may be included within or external of the processor, and may be implemented by any conventional or other memory unit with any suitable storage capacity and any type of memory (e.g., random access memory (RAM), read only memory (ROM), etc.). The memory may store any desired information for performing the geolocation technique of present invention embodiments (e.g., location modules, data, etc.). The interface unit may be implemented by any quantity of any conventional or other communications device (e.g., wireless communications device, wired communication device, etc.), and may be configured for communication over any desired network (e.g., wireless, cellular, LAN, WAN, Internet, Intranet, VPN, etc.).
Present invention embodiments may employ any quantity of variables or equations to determine the estimated location of one or more emitters, provided that the quantity of equations is greater than or equal to the quantity of unknown variables. Further, any conventional or other techniques may be employed to produce the location estimate with minimal error (e.g., Least Mean Square (LMS), etc.). The equations may be represented in any desired form (e.g., matrix form, vectors, scalars, etc.), and be solved in any desired fashion to enable determination of the emitter location. The location estimate may be produced and/or converted to any desired form, and may be provided with respect to any desired reference (e.g., coordinates within the space, longitude and latitude indications, GPS coordinates, etc.).
The resulting location estimate may be utilized for any suitable applications (e.g., generation of a map image of the area, vehicle or other platform guidance systems to direct the vehicle or platform toward or away from areas, radar or other detection systems, etc.).
The various indices (e.g., i, N, etc.) are preferably integers, but may be any types of numbers with any suitable numeric ranges.
It is to be understood that the terms “top”, “bottom”, “front”, “rear”, “side”, “height”, “length”, “width”, “upper”, “lower”, “vertical” and the like are used herein merely to describe points of reference and do not limit the present invention to any particular orientation or configuration.
From the foregoing description, it will be appreciated that the invention makes available a novel system and method for three-dimensional geolocation of emitters based on energy or received signal strength (RSS) measurements, wherein locations of radio frequency (RF) emitters in a three-dimensional space are determined based on energy or received signal strength (RSS) measurements of the emitters at various locations.
Having described example embodiments of a new and improved system and method for three-dimensional geolocation of emitters based on energy or received signal strength (RSS) measurements, it is believed that other modifications, variations and changes will be suggested to those skilled in the art in view of the teachings set forth herein. It is therefore to be understood that all such variations, modifications and changes are believed to fall within the scope of the present invention as defined by the appended claims.

## Claims

1. A system for locating an emitter within an area comprising:
an unmanned air vehicle that traverses a pre-planned path through said area to obtain signal strength measurements at a plurality of locations within said area with a line of sight condition that reduces effects from signal fading and shadowing, wherein said unmanned air vehicle includes a receiver to receive signals transmitted by said emitter and obtain a series of measurements of signal strength for said received signals at each of said plurality of different locations within said area; and
a processor to process said signal strength measurements to locate said emitter within said area, wherein said processor includes:
a location module to process said signal strength measurements and select for each location the signal strength measurement with a maximum signal to noise ratio from the series of said signal strength measurements for that location, and determine a three-dimensional location of said emitter within said area based on relationships of distances between said emitter and each of said plurality of locations, wherein said distances are unknown, said selected signal strength measurements are proportional to said unknown distances, and said location module includes:
a variable module to determine said three-dimensional location by solving a set of simultaneous equations relating to said distances, wherein said set of simultaneous equations includes a set of variables representing said unknown distances between said emitter and said plurality of locations, and unknown variables representing coordinates of said three-dimensional location of said emitter within said area, and wherein said set of variables representing said unknown distances within said set of simultaneous equations includes values that are approximated by corresponding selected signal strength measurements to solve for said unknown variables and determine said three-dimensional location.

2. The system of claim 1, wherein said variable module solves said set of simultaneous equations using a Least Mean Square technique.

3. The system of claim 1, wherein said location module includes:
a parameter module to determine a parameter for said set of simultaneous equations based on said selected signal strength measurements.

4. The system of claim 3, wherein said parameter is based on a ratio of said selected signal strength measurements.

5. The system of claim 1, wherein said set of simultaneous equations includes:



  
    
      
        [
        
          
            
              
                r
                0
                2
              
            
          
          
            
              x
            
          
          
            
              y
            
          
          
            
              z
            
          
        
        ]
      
      =
      
        
          
            (
            
              
                P
                T
              
              ⁢
              P
            
            )
          
          
            -
            1
          
        
        ⁢
        
          P
          T
        
        ⁢
        R
      
    
    ,
    where
  




  
    
      P
      =
      
        [
        
          
            
              
                β
                1
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      x
                      1
                    
                    -
                    
                      x
                      0
                    
                  
                  )
                
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      y
                      1
                    
                    -
                    
                      y
                      0
                    
                  
                  )
                
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      z
                      1
                    
                    -
                    
                      z
                      0
                    
                  
                  )
                
              
            
          
          
            
              
                β
                2
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      x
                      2
                    
                    -
                    
                      x
                      0
                    
                  
                  )
                
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      y
                      2
                    
                    -
                    
                      y
                      0
                    
                  
                  )
                
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      z
                      2
                    
                    -
                    
                      z
                      0
                    
                  
                  )
                
              
            
          
          
            
              ...
            
            
              ...
            
            
              ...
            
            
              ...
            
          
          
            
              ...
            
            
              ...
            
            
              ...
            
            
              ...
            
          
          
            
              
                β
                N
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      x
                      N
                    
                    -
                    
                      x
                      0
                    
                  
                  )
                
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      y
                      N
                    
                    -
                    
                      y
                      0
                    
                  
                  )
                
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      z
                      N
                    
                    -
                    
                      z
                      0
                    
                  
                  )
                
              
            
          
        
        ]
      
    
    ,
    
      

    
    ⁢
    
      R
      =
      
        [
        
          
            
              
                
                  d
                  1
                  2
                
                -
                
                  d
                  0
                  2
                
              
            
          
          
            
              
                
                  d
                  2
                  2
                
                -
                
                  d
                  0
                  2
                
              
            
          
          
            
              ...
            
          
          
            
              ...
            
          
          
            
              
                
                  d
                  N
                  2
                
                -
                
                  d
                  0
                  2
                
              
            
          
        
        ]
      
    
    ,
    
      

    
    ⁢
    
      
        β
        i
      
      =
      
        [
        
          
            
              p
              0
            
            
              p
              i
            
          
          -
          1
        
        ]
      
    
    ,
  


and wherein x, y, z represent said coordinates of said three-dimensional location of said emitter within said area, r02 represents a square of a distance between said emitter and a reference location, xi, yi, zi, (for i=0 to N) represent coordinates of said locations for obtaining said signal strength measurements, di2 (for i=0 to N) represents a square of a distance between a reference origin in said area and said locations for obtaining said signal strength measurements, pi (for i=0 to N) represents a selected signal strength measurement at said locations, and N represents a quantity of said locations for obtaining said signal strength measurements.

6. The system of claim 1, wherein said signals transmitted by said emitter include radio frequency signals.

7. A method for locating an emitter within an area comprising:
(a) receiving signals transmitted by said emitter via a receiver and obtaining a series of measurements of signal strength for said received signals at each of a plurality of different locations within said area, wherein said receiver is mounted on an unmanned air vehicle, and step (a) further includes:
(a.1) traversing a pre-planned path through said area with said unmanned air vehicle to obtain said signal strength measurements at said plurality of locations with a line of sight condition that reduces effects from signal fading and shadowing; and

(b) processing said signal strength measurements, via a processor, and selecting for each location the signal strength measurement with a maximum signal to noise ratio from the series of said signal strength measurements for that location, and determining a three-dimensional location of said emitter within said area based on relationships of distances between said emitter and each of said plurality of locations, wherein said distances are unknown, said selected signal strength measurements are proportional to said unknown distances, and step (b) further includes:
(b.1) determining said three-dimensional location by solving a set of simultaneous equations relating to said distances, wherein said set of simultaneous equations includes a set of variables representing said unknown distances between said emitter and said plurality of locations, and unknown variables representing coordinates of said three-dimensional location of said emitter within said area, and wherein said set of variables representing said unknown distances within said set of simultaneous equations includes values that are approximated by corresponding selected signal strength measurements to solve for said unknown variables and determine said three-dimensional location.

8. The method of claim 7, wherein step (b.1) further includes:
(b.1.1) solving said set of simultaneous equations using a Least Mean Square technique.

9. The method of claim 7, wherein step (b) further includes:
determining a parameter for said set of simultaneous equations based on said selected signal strength measurements.

10. The method of claim 9, wherein said parameter is based on a ratio of said selected signal strength measurements.

11. The method of claim 7, wherein said set of simultaneous equations includes:



  
    
      
        [
        
          
            
              
                r
                0
                2
              
            
          
          
            
              x
            
          
          
            
              y
            
          
          
            
              z
            
          
        
        ]
      
      =
      
        
          
            (
            
              
                P
                T
              
              ⁢
              P
            
            )
          
          
            -
            1
          
        
        ⁢
        
          P
          T
        
        ⁢
        R
      
    
    ,
    where
  




  
    
      P
      =
      
        [
        
          
            
              
                β
                1
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      x
                      1
                    
                    -
                    
                      x
                      0
                    
                  
                  )
                
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      y
                      1
                    
                    -
                    
                      y
                      0
                    
                  
                  )
                
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      z
                      1
                    
                    -
                    
                      z
                      0
                    
                  
                  )
                
              
            
          
          
            
              
                β
                2
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      x
                      2
                    
                    -
                    
                      x
                      0
                    
                  
                  )
                
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      y
                      2
                    
                    -
                    
                      y
                      0
                    
                  
                  )
                
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      z
                      2
                    
                    -
                    
                      z
                      0
                    
                  
                  )
                
              
            
          
          
            
              ...
            
            
              ...
            
            
              ...
            
            
              ...
            
          
          
            
              ...
            
            
              ...
            
            
              ...
            
            
              ...
            
          
          
            
              
                β
                N
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      x
                      N
                    
                    -
                    
                      x
                      0
                    
                  
                  )
                
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      y
                      N
                    
                    -
                    
                      y
                      0
                    
                  
                  )
                
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      z
                      N
                    
                    -
                    
                      z
                      0
                    
                  
                  )
                
              
            
          
        
        ]
      
    
    ,
    
      

    
    ⁢
    
      R
      =
      
        [
        
          
            
              
                
                  d
                  1
                  2
                
                -
                
                  d
                  0
                  2
                
              
            
          
          
            
              
                
                  d
                  2
                  2
                
                -
                
                  d
                  0
                  2
                
              
            
          
          
            
              ...
            
          
          
            
              ...
            
          
          
            
              
                
                  d
                  N
                  2
                
                -
                
                  d
                  0
                  2
                
              
            
          
        
        ]
      
    
    ,
    
      

    
    ⁢
    
      
        β
        i
      
      =
      
        [
        
          
            
              p
              0
            
            
              p
              i
            
          
          -
          1
        
        ]
      
    
    ,
  


and wherein x, y, z represent said coordinates of said three-dimensional location of said emitter within said area, r02 represents a square of a distance between said emitter and a reference location, xi, yi, zi, (for i=0 to N) represent coordinates of said locations for obtaining said signal strength measurements, di2 (for i=0 to N) represents a square of a distance between a reference origin in said area and said locations for obtaining said signal strength measurements, pi (for i=0 to N) represents a selected signal strength measurement at said locations, and N represents a quantity of said locations for obtaining said signal strength measurements.

12. The method of claim 7, wherein said signals transmitted by said emitter include radio frequency signals.

13. A program product apparatus comprising a computer readable memory device with computer program logic recorded thereon for locating an emitter within an area, said program product apparatus further comprising:
a location module to process signal strength measurements of signals transmitted by said emitter and obtained at a plurality of different locations within said area and select for each location the signal strength measurement with a maximum signal to noise ratio from a series of said signal strength measurements obtained for that location, and to determine a three-dimensional location of said emitter within said area based on relationships of distances between said emitter and each of said plurality of locations, wherein said signal strength measurements are obtained by an unmanned air vehicle that traverses a pre-planned path through said area to obtain said signal strength measurements at said plurality of locations with a line of sight condition that reduces effects from signal fading and shadowing, and wherein said distances are unknown, said selected signal strength measurements are proportional to said unknown distances, and said location module includes:
a variable module to determine said three-dimensional location by solving a set of simultaneous equations relating to said distances, wherein said set of simultaneous equations includes a set of variables representing said unknown distances between said emitter and said plurality of locations, and unknown variables representing coordinates of said three-dimensional location of said emitter within said area, and wherein said set of variables representing said unknown distances within said set of simultaneous equations includes values that are approximated by corresponding selected signal strength measurements to solve for said unknown variables and determine said three-dimensional location.

14. The apparatus of claim 13, wherein said variable module solves said set of simultaneous equations using a Least Mean Square technique.

15. The apparatus of claim 13, wherein said location module includes:
a parameter module to determine a parameter for said set of simultaneous equations based on said selected signal strength measurements.

16. The apparatus of claim 15, wherein said parameter is based on a ratio of said selected signal strength measurements.

17. The apparatus of claim 13, wherein said set of simultaneous equations includes:



  
    
      
        [
        
          
            
              
                r
                0
                2
              
            
          
          
            
              x
            
          
          
            
              y
            
          
          
            
              z
            
          
        
        ]
      
      =
      
        
          
            (
            
              
                P
                T
              
              ⁢
              P
            
            )
          
          
            -
            1
          
        
        ⁢
        
          P
          T
        
        ⁢
        R
      
    
    ,
    where
  




  
    
      P
      =
      
        [
        
          
            
              
                β
                1
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      x
                      1
                    
                    -
                    
                      x
                      0
                    
                  
                  )
                
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      y
                      1
                    
                    -
                    
                      y
                      0
                    
                  
                  )
                
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      z
                      1
                    
                    -
                    
                      z
                      0
                    
                  
                  )
                
              
            
          
          
            
              
                β
                2
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      x
                      2
                    
                    -
                    
                      x
                      0
                    
                  
                  )
                
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      y
                      2
                    
                    -
                    
                      y
                      0
                    
                  
                  )
                
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      z
                      2
                    
                    -
                    
                      z
                      0
                    
                  
                  )
                
              
            
          
          
            
              ...
            
            
              ...
            
            
              ...
            
            
              ...
            
          
          
            
              ...
            
            
              ...
            
            
              ...
            
            
              ...
            
          
          
            
              
                β
                N
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      x
                      N
                    
                    -
                    
                      x
                      0
                    
                  
                  )
                
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      y
                      N
                    
                    -
                    
                      y
                      0
                    
                  
                  )
                
              
            
            
              
                2
                ⁢
                
                  (
                  
                    
                      z
                      N
                    
                    -
                    
                      z
                      0
                    
                  
                  )
                
              
            
          
        
        ]
      
    
    ,
    
      

    
    ⁢
    
      R
      =
      
        [
        
          
            
              
                
                  d
                  1
                  2
                
                -
                
                  d
                  0
                  2
                
              
            
          
          
            
              
                
                  d
                  2
                  2
                
                -
                
                  d
                  0
                  2
                
              
            
          
          
            
              ...
            
          
          
            
              ...
            
          
          
            
              
                
                  d
                  N
                  2
                
                -
                
                  d
                  0
                  2
                
              
            
          
        
        ]
      
    
    ,
    
      

    
    ⁢
    
      
        β
        i
      
      =
      
        [
        
          
            
              p
              0
            
            
              p
              i
            
          
          -
          1
        
        ]
      
    
    ,
  


and wherein x, y, z represent said coordinates of said three-dimensional location of said emitter within said area, r02 represents a square of a distance between said emitter and a reference location, xi, yi, zi (for i=0 to N) represent coordinates of said locations for obtaining said signal strength measurements, di2 (for i=0 to N) represents a square of a distance between a reference origin in said area and said locations for obtaining said signal strength measurements, pi (for i=0 to N) represents a selected signal strength measurement at said locations, and N represents a quantity of said locations for obtaining said measurements.

18. The apparatus of claim 13, wherein said signals transmitted by said emitter include radio frequency signals.