System and methods for non-parametric technique based geolocation and cognitive sensor activation
09804253 · 2017-10-31
Assignee
Inventors
Cpc classification
G01S5/28
PHYSICS
International classification
Abstract
The present invention relates to a geolocation system and method for a multi-path environment. The geolocation system comprises one or more emitters (201a . . . 201n), one or more sensors (202a . . . 202n) comprising at least one processor. A first processor (204) estimates angle of arrival (AOA) and time of arrival (TOA) from the signals received from said one or more emitters (201a . . . 201n). A second processor (205) determines clusters based on the (AOA) and (TOA) data. The system also comprises a central node (207) in communication with at least one sensor (202a . . . 202n) and configured to estimate geolocation of one or more emitters (201a . . . 201n) wherein, said second processor (205) clusters data for the one or more emitters (201a . . . 201n) by executing a non-parametric Bayesian technique and said central node (207) utilizes hybrid angle of arrival-time difference of arrival (AOA-TDOA) technique to determine geolocation of each of the emitters (201a . . . 201n).
Claims
1. A geolocation system for a multi-path environment comprising: one or more emitters; one or more sensors comprising at least one processor; a first processor configured to estimate an angle of arrival (AOA) and a time of arrival (TOA) from the signals received from the one or more emitters; a second processor configured to determine clusters based on the AOA and TOA data received from said one or more emitters; a central node in communication with at least one sensor and configured to estimate geolocation and a next time instant position of the one or more emitters based on the determined clusters; and a cognitive sensor activation unit in communication with the central node, the cognitive sensor activation unit configured to selectively activate the one or more sensors to detect the corresponding one or more emitters based on the estimated geolocation and next time instant position of the one or more emitters; wherein said second processor clusters data for the one or more emitters in said multi-path environment by executing a non-parametric Bayesian technique; and said central node utilizes hybrid angle of arrival-time difference of arrival (AOA-TDOA) technique to determine the geolocation of each of the one or more emitters; wherein said second processor utilizes a Dirichlet process mixture model (DPMM) technique to determine an actual number of emitters comprising the one or more emitters.
2. The geolocation system as claimed in claim 1, wherein the first processor includes a digital signal processor (DSP), and the second processor includes a general purpose processor (GPP).
3. The geolocation system as claimed in claim 1, wherein each of said sensors comprises an antenna array configured to receive signals in the multi-path environment.
4. The geolocation system as claimed in claim 1, wherein each of said sensor comprises a transceiver configured to transmit the clustered data to the central node.
5. The geolocation system as claimed in claim 1, wherein said central node and said cognitive sensor activation unit are connected to a graphical user interface.
6. A method for providing geolocation information in a multi-path environment comprising the steps of: receiving a mixture of signals from one or more emitters at one or more sensors in said multi-path environment; determining angle of arrival (AOA) and time of arrival (TOA) data from each of the received signals in said multi-path environment; clustering the mixture of received signals based on the determined AOA and TOA data into one or more clusters arising from one or more emitters utilizing a non-parametric Bayesian technique; transmitting by the one or more sensors the one or more clusters to a central node for further processing; estimating by the central node the time difference of arrival (TDOA) profile for each of the one or more emitters based on the one or more clusters and extracting the AOA data for said emitters from within the distinct one or more clusters; processing by the central node the AOA data and TDOA profile data for each of the one or more emitters to perform a hybrid AOA-TDOA technique for determining the geolocation of the one or more emitters; wherein the central node further utilizes a Dirichlet process mixture model (DPMM) technique to determine an actual number of emitters comprising the one or more emitters; and selectively activating the one or more sensors to detect the corresponding one or more emitters based on the determined geolocation of the one or more emitters.
7. The method as claimed in claim 6, comprising displaying geolocation data on a graphical user interface.
8. A method for selectively activating one or more sensors for tracking one or more emitters in a multi-path environment, comprising the steps of: receiving signals from one or more emitters at one or more sensors in said multi-path environment; determining angle of arrival (AOA) and time of arrival (TOA) from each of the received signals in said multi-path environment; clustering into clusters the mixture of received signals arising from one or more emitters utilizing a non-parametric Bayesian technique; transmitting by the one or more sensors said clusters to a central node for further processing; estimating by the central node the time difference of arrival (TDOA) profile for each of the one or more emitters and extracting the AOA data for said emitters from within the distinct clusters; processing by the central node the TDOA profile and AOA data information for each of the one or more emitters to perform a hybrid AOA-TDOA technique for determining the geolocation of the one or more emitters, wherein the central node further utilizes a Dirichlet process mixture model (DPMM) technique to determine an actual number of emitters comprising the one or more emitters; tracking the position and a next time instant position of each of the emitters based on said geolocation of the one or more emitters; and selectively activating the one or more sensors in proximity to one or more emitters based on the determined geolocation and next time instant position of the one or more emitters.
9. The method as claimed in claim 8, further comprising discretizing the multi-path environment into Voronoi regions based on each sensor being located at a center of a corresponding Voronoi region; wherein selectively activating the one or more sensors includes activating the one or more sensors based on the Voronoi region discretization.
10. The method as claimed in claim 8, wherein selectively activating the one or more sensors in proximity to a targeted emitter reduces the sensor transmission power and noise contributions.
Description
BRIEF DESCRIPTION OF DRAWINGS
(1)
(2)
(3)
(4)
(5)
(6)
(7)
(8)
DETAILED DESCRIPTION OF THE INVENTION
(9)
(10) Referring to
(11) The sensor also has a field programmable gate array or FPGA board which may embody a digital signal processor (204) which performs the angle of arrival-time difference of arrival (AOA-TDOA) estimation. The sensor may further comprise a general purpose power PC or a general purpose processor (GPP) (205) for implementing a non-parametric Bayesian clustering technique according to an embodiment of the present invention using a hardware description language (HDL) programming. In an alternate embodiment, the clustering technique may be performed by the digital signal processor (204) which is also responsible for correctly estimating the number of emitters based on the AOA-TOA data. For example, assuming that z is a set of multi-path signals from one or more emitters received at the antenna array (203), if z consists of sequence of observations z.sub.1 to z.sub.n, then processor (204) estimates the TOA and AOA data from the down-converted signal received from the transceiver (206).
(12) The sensor (202) is operably connected to a central node (207). In an alternate embodiment, the central node (207) may be located at a remote location and may be in wireless communication with the plurality of sensors (202a . . . 202n). The central node (207) is further connected to a cognitive sensor activation unit (208) which may be configured to execute geolocation techniques for identifying emitter locations from the clusters of observations and TDOA profiles. For example, and as described in reference to
(13)
(14) Next, in Step 303, the sensors (202a, . . . , 202n) perform clustering individually to discriminate distinct emitter contributions within the AOA-TOA data obtained from processor (204). In the exemplary embodiment, in the general purpose processor (205), a non-parametric Bayesian clustering technique called Dirichlet Process Mixture Model (DPMM) based clustering of AOA-TOA data is employed for obtaining the actual number of emitters in a environment. The number of emitters and their received signal distribution parameters are unknown a priori, and are to be inferred directly from the data. Using a Dirichlet Process (DP) permits the assumption of a priori a countably infinite number of clusters, and to directly utilize the raw multi-path data in order to make an accurate estimate of the number of emitters without relying on any initial assumptions.
(15) An exemplary DPMM technique for identifying the cluster of observations in accordance with an exemplary embodiment is illustrated hereinafter.
(16) If z is a set of AOA-TOA components derived from multi-path components or observation set received at antenna array, AOA-TOA multi-path data at each sensor is processed by processor (204) to obtain sequence of data Z.sub.1, . . . , Z.sub.n. Thus, for each i=1, . . . , n, z.sub.i={θ.sub.i, φ.sub.i, τ.sub.i}, where θ.sub.i is the azimuth angle data associated with a particular multi-path component of an emitter (multi-path component LOS or NLOS) contribution and φ.sub.i is the corresponding elevation angle information and τ.sub.i is the TOA of the multi-path component.
(17) In the exemplary embodiment, it is assumed that z.sub.i has a multivariate Gaussian distribution with mean vector μ.sub.i and covariance matrix σ.sub.i. If ψ.sub.i={μ.sub.i, σ.sub.i} is the parameter of interest for data z.sub.i, the goal will be to find the posterior distribution of (ψ.sub.1, . . . , ψ.sub.n) given the data (z.sub.1, . . . , z.sub.n). If we make sequence of observations (z.sub.1, . . . , z.sub.n), where for each i=1, . . . , n, z.sub.i˜F(•|ψ.sub.i), and ψ.sub.iεΨ is a parameter describing the observation distribution, employing a Dirichlet process of clustering the data using Gibbs sampling technique, this prior is chosen to be a stochastic process which leads to a model with very rich features. Specifically, the Dirichlet Process is a distribution over the space of all probability measures on Ψ. A random distribution G on Ψ is then drawn from this distribution, and given G, the parameters Ψ.sub.i, i=1, . . . , n, are independent and identically distributed according to G. To define the Dirichlet Process, we first let G.sub.0 be a probability distribution over Ψ, which represents our prior belief about a parameter, and α be a positive number that serves as a weight between our prior belief and information inferred from observed data. We say that G is distributed as a Dirichlet Process, denoted as G˜DP(α, G.sub.0), if for any finite measurable partition A.sub.1, . . . , A.sub.r of Ψ, we have
(G(A.sub.1), . . . ,G(A.sub.r))˜Dir(αG.sub.0(A.sub.1), . . . ,αG.sub.0(A.sub.r)),
where Dir(•) is the Dirichlet distribution. From this definition, we see that the DP is a stochastic process. Thus, the DPMM has the following representation,
G˜DP(α,G.sub.0),
ψ.sub.i|G˜G,
z.sub.i|Ψ.sub.i˜F(•|ψ.sub.i).
(18) Let ψ.sub.−i=(ψ.sub.1, . . . , ψ.sub.i−1, ψ.sub.i+1, . . . , ψ.sub.n) be the vector of parameters excluding ψ.sub.i. Assuming that all distributions have a density with respect to some dominating σ-finite measure and using the same symbols to denote the distribution as well as the density, the posterior distribution of ψ.sub.i, conditioned on the data z and ψ.sub.−i is then given by
p(ψ.sub.i|ψ.sub.−i,z.sub.i)∞F(z.sub.i|ψ.sub.i)p(ψ.sub.i|ψ.sub.−i), (2)
since given ψ.sub.−i, ψ.sub.i depends only on z.sub.i. From the Blackwell-MacQueen Polya-Urn scheme, the conditional distribution of ψ.sub.i given ψ.sub.−i is
(19)
where δ.sub.ψ is the Dirac delta function at ψ. Thus the posterior distribution is given by
(20)
where β=1/(αq.sub.0+Σ.sub.j≠iF(z.sub.i|ψ.sub.j)) is a normalizing constant, and
q.sub.0=∫G.sub.0(ψ)F(z.sub.i|ψ)dψ (5)
is the marginal density of z.sub.i at its realization. In order to evaluate the integral (5), we choose G.sub.0 to be a conjugate prior to the Gaussian distribution F(z.sub.i|ψ.sub.i). In this work, the Normal-Wishart distribution for G.sub.0 is used. A Gibbs sampler can now be designed to obtain the posterior distribution of ψ.sub.i given all the data. Let p.sub.z(ψ.sub.i|ψ.sub.−i) be the conditional distribution of ψ.sub.i given all the data z. We sample ψ.sub.i according to
(21)
(22) We initialize the Gibbs sampler by considering each data z.sub.i as being in its own set, with ψ.sub.i.sup.(0)=z.sub.i. Subsequently the Gibbs sampling for the k th step is done in the following way. Sample ψ.sub.1.sup.k from p.sub.z(•|ψ.sub.2=ψ.sub.2.sup.(k-1), ψ.sub.3=ψ.sub.3.sup.(k-1), . . . , ψ.sub.n=ψ.sub.n.sup.(k-1)) Sample ψ.sub.2.sup.k from p.sub.z(•|ψ.sub.1=ψ.sub.1.sup.(k), ψ.sub.3=ψ.sub.3.sup.(k-1), . . . , ψ.sub.n=ψ.sub.n.sup.(k-1)) Sample ψ.sub.n.sup.k from P.sub.z=(•|ψ.sub.1=ψ.sub.1.sup.(k), ψ.sub.2=ψ.sub.2.sup.(k), . . . , ψ.sub.n-1=ψ.sub.n-1.sup.(k))
(23) The DPMM technique of the present invention detects the distinct components or clusters in the raw data without making any a priori assumptions about the number of emitters or any candidate locations.
(24) Next, in Step 304, the sensors collaborate and evaluate the TDOA profiles for each emitter signals based on the clustered TOA in the previous step. Once a set of optimum AOA-TOA clusters is determined and TDOA profile is obtained at each node, the clusters and TDOA profile are then sent to the central node (207) for data fusion and further processing (Step 305). In Step 306, the central node (207) then establishes each emitter position based on the AOA-TDOA information received. This is done by, for example, evaluating a least square solution based on the TDOA data. This coarse TDOA estimate is then refined in Step 307 by a hybrid AOA-TDOA technique.
(25) In an exemplary embodiment, the central node (207) executes the hybrid (AOA)-(TDOA) technique and Levenberg-Marquardt technique and utilizes the data θ.sub.mean, φ.sub.mean, τ.sub.mean, which are obtained from cluster of discriminated emitters determined using non-parametric Bayesian technique, preferably a Dirichlet Process Mixture Model (DPMM) technique, to achieve the geolocation of each of the emitters in a multipath environment. The inferred set of parameters Ψ={ψ.sub.1.sup.k, ψ.sub.2.sup.k, . . . , ψ.sub.n.sup.k} is utilized to form hybrid geolocation estimate. As shown in
(26) Based on the above hybrid angle of arrival (AOA) and time difference of arrival (TDOA), the central node predicts the position of each of the emitter and in the next time instant by utilizing tracking technique like extended Kalman filtering, and then works out a sensor activation process in order to select the best sensors in the proximity of the emitter to be selected for sensing based on the Voronoi region discretization (Step 308). The Voronoi region discretization of the embodiment identifies the nearest sensor to the emitter node; for example if the emitter location falls within a Voronoi cell/region, then the nearest sensor is the sensor at the center of the Voronoi region. This property of placing the sensors at the centers of Voronoi cells helps to identify the nearest set of sensors around the target emitter.
(27) Thus, in this exemplary approach we track the location of the emitter before trying to detect its presence. This track-before-detect approach conserves crucial sensor power and also contributes in reducing false alarms and noise contributions given by redundant and far away sensors. In an alternative approach, based on this position prediction the central node (207) activates the most optimum set of sensors and the whole process of discrimination and prediction is repeated. In an advantageous embodiment, feedback from the individual sensors (202a, . . . , 202n) to the central node (207) enables the central node to perceive the sensor node network and this cognition would help the sensor node unit to design adaptive beamforming, power and spectrum allocation strategies, to enhance the geolocation performance and mitigate interference and jamming.
(28) In yet another advantageous embodiment, the cognitive sensor activation framework according to present invention enables to conserve the crucial sensor node transmission power and reduces false alarm and noise contributions if redundant sensors are activated in the next time instant.
(29)
(30) The DPMM technique of the present invention has several advantages as compared to known systems and methods. It makes the geolocation system adaptable by providing a unique machine learning approach to make inferences directly from the observed data, without having to resort to any intermediate parameters for estimation of the mixing components. Moreover it adapts its estimation mechanism on the fly in accordance with the observed data (by updating the posterior distribution with each observed data instance during Gibbs Sampling as shown previously).
(31) Advantageously, the DPMM technique of the present invention is also very flexible as it presents a very diffused model which can accommodate higher dimensional data and provides a framework to perform unsupervised mixture component analysis without the need of any a priori information on the number of emitters or the mixing proportion of their AOA-TOA data. Processing data utilizing the DPMM technique reveals the underlying latent mixing components through an iterative procedure taking into account each data instance at a time through update of the posterior distribution. The clustering approach of the DPMM technique is also flexible in terms of accommodating higher dimensional data, for example as stated previously, emitter discrimination can be based on multivariate data (elevation angle, azimuth angle and TOA).
(32) The DPMM technique described hereinabove can also accommodate further data parameters like received signal strength, frequency of arrival to facilitate a better clustering of mixing components. Unlike other approaches which take into account or need sensor-scene-specific information (scatterer layout, topology, maps), DPMM provides a very generic model which can be used to learn about any sensor node network under consideration.
(33) Advantageously, the DPMM technique of the present invention enhances the emitter discrimination capability resulting in improved geolocation estimate for each emitter. For example, in the exemplary embodiment, an accuracy of approximately 20 m to 25 m in 3D could be achieved with TOA spread of 1 μsec and AOA spread of 15° in the multi-path components of different emitters in an urban topology of 2 km×2 km×2 km. The DPMM technique is robust and relies upon coarse TDOA estimate to aid the hybrid estimation and thus can handle worse AOA spreads (more than 20°) without affecting the geolocation accuracy.
(34)
(35)
(36)
(37)
(38) From the above it is evident that the cognitive sensor activation framework would activate only those sensor nodes which are in proximity to detected emitters thereby improving the working of cognitive sensor activation framework.
(39) While preferred embodiments of the present subject matter have been described, it is to be understood that the embodiments described are illustrative only and that the scope of the invention is to be defined solely by the appended claims when accorded a full range of equivalence, many variations and modifications naturally occurring to those of skill in the art from a perusal hereof.