The present invention relates generally to data fusion. It particularly relates to a data fusion method that provides an adaptive weighting technique using reliability functions to integrate data from a plurality of sensors.
Sensor systems incorporating a plurality of sensors (multi-sensor systems) are widely used for a variety of military applications including ocean surveillance, air-to-air and surface-to-air defense (e.g., self-guided munitions), battlefield intelligence, surveillance and target detection (classification), and strategic warning and defense. Also, multi-sensor systems are used for a plurality of civilian applications including condition-based maintenance, robotics, automotive safety, remote sensing, weather forecasting, medical diagnoses, and environmental monitoring (e.g., weather forecasting).
To obtain the full advantage of a multi-sensor system, an efficient data fusion method (or architecture) may be selected to optimally combine the received data from the multiple sensors. For military applications (especially target recognition), a sensor-level fusion process is widely used wherein data received by each individual sensor is fully processed at each sensor before being output to a system data fusion processor. The data (signal) processing performed at each sensor may include a plurality of processing techniques to obtain desired system outputs (target reporting data) such as feature extraction, and target classification, identification, and tracking. The processing techniques may include time-domain, frequency-domain, multi-image pixel image processing techniques, and/or other techniques to obtain the desired target reporting data.
An exemplary, prior art example of a multi-sensor, sensor-level fusion (process) system 100 for automatic target recognition (ATR) is shown in
Many multi-sensor systems (such as system 100 in
Currently, a data fusion method (strategy) that is widely used for feature-level systems is multiplicative fusion (e.g., Bayes or Dempster-Shafer methods). Commonly, the multiplicative fusion method multiplies a plurality of probability functions (generated from the received data from each individual sensor) to produce a single term (value). The generation of the single term makes it complex to weight contributions from the plurality of sensors (which may have different reliability values over different tracking time periods due to different sensor constraints, atmospheric conditions, or other factors) and thus may produce a less accurate data fusion output (decision output regarding target classification). Additionally, under certain conditions, a data fusion reliability output (using data from all sensors) may be worse than a single sensor reliability output.
Therefore, due to the disadvantages of the current multiplicative data fusion method, there is a need to provide a multi-sensor system that uses an additive data fusion method to produce multiple terms for weighting and determines a current, better performing (higher reliability) sensor to adaptively weight the contributions from the plurality of different sensors for improved reliability of target classification. Also, there is a need to provide a multi-sensor data fusion system that can adaptively weight multi-sensor reliability towards the better performing sensor (and away from a worse performing sensor) when predetermined conditions arise making the better single sensor reliability higher than the data fusion (combining all sensor data) reliability.
The method and system of the present invention overcome the previously mentioned problems by providing a multi-sensor data fusion system capable of adaptively weighting the contributions from a plurality of sensors in the system using an additive calculation of a sensor reliability function for each sensor. During a predetermined tracking period, data is received from each individual sensor in the system and a sensor reliability function is determined for each sensor based on the SNR (signal-to-noise ratio) for the received data from each sensor. Each sensor reliability function is individually weighted based on the SNR for each sensor and a comparison of predetermined sensor operation characteristics for each sensor and a better performing (higher reliability) single sensor. Additive calculations are performed on the reliability functions to produce both an absolute and a relative reliability function which provide a confidence level for the multi-sensor system relating to the correct classification (recognition) of targets and decoys.
a–3b show graphs of exemplary sensor performance results for a multi-sensor data fusion system in accordance with embodiments of the present invention.
a–4b show diagrams of exemplary sensor reliability results for a multi-sensor data fusion system in accordance with embodiments of the present invention.
a–5b show diagrams of exemplary sensor classification results for a multi-sensor system using multiplicative and additive data fusion in accordance with embodiments of the present invention.
a–6b show diagrams of exemplary sensor classification results for a multi-sensor system using additive data fusion in accordance with embodiments of the present invention.
a–7d show further diagrams of exemplary sensor classification for a multi-sensor system using additive fusion in accordance with embodiments of the present invention.
The plurality of sensors 205, using associated sensor processors, may each perform the well-known process of feature extraction to detect and pull out features which help discriminate the objects in each sensor's field of view and combine all the feature extractions (from each sensor) as a composite input to data fusion processor 208. Data fusion processor 208 may perform, as described in detail later, all levels of discrimination (detection, classification—recognition, identification, and tracking) of the object (target) using a predetermined data fusion algorithm to recognize the object of interest, differentiate the object from decoys (false targets), and produce at least one (system) weighted, reliability function that links the observed object to a predetermined target with some confidence level. The system reliability function may be used to generate a decision output 210 (target report) for target detection such as “validated target” or “no desired target encountered”. Also, alternatively, plurality of sensors 205 may feed-through (without processing or with minimal processing) received data to processor 208 for feature extraction and target discrimination processing.
The particular combination of sensors 205 for system 200 may include a number of different sensors selected to provide exemplary predetermined system attributes (parameters) including temporal and spatial diversity (fusion), sensitivity, bandwidth, noise, operating range, transmit power, spatial resolution, polarization, and other system attributes. These different sensors may include, but are not limited to, passive and/or active sensors operating in the RF (radio frequency) range such as MMW (millimeter-wave) sensors, IR (infrared) sensors (e.g., Indium/Antimony—InSb focal plane array), laser sensors, and other passive and/or active sensors useful in providing the exemplary predetermined system attributes.
During exemplary operation as described herein and in accordance with the flow process diagram shown in
For multi-sensor system 200, there may be variations in sensor reliability among the plurality of sensors 205 (e.g., based on variations in SNR and other factors) during the tracking period such that the data fusion processor 208 (when performing data fusion) may determine and assign a higher weight to a best performing sensor (with the highest SNR) than a (lower) weight assigned to a worse (or worst) performing sensor (e.g., with a lower SNR) such that a fused result (combined reliability function for the plurality of sensors) may be weighted more towards the best performing (highest reliability) sensor. The variations in sensor reliabilities for the plurality of sensors 205 may be caused by a number of factors including weather conditions, different sensor attributes such as better range accuracy of an RF sensor than an IR sensor at longer ranges, or other factors causing at least one sensor to perform better than another sensor during a predetermined tracking period.
Advantageously during operation as described herein, the SNR may be used by data fusion processor 208 as a measure of sensor reliability during a predetermined tracking period to help generate a sensor reliability function for each one of the plurality of sensors 205. Thereafter, data fusion processor 208 may execute (perform) a predetermined data fusion algorithm incorporating additive and/or multiplicative calculation (of each individual sensor reliability function) to generate at least one overall (combined) reliability function for the multi-sensor system (full plurality of sensors). As part of generating the overall reliability function (for the plurality of sensors) in accordance with the fusion algorithm (process), data fusion processor 208 may adaptively weight (for a predetermined number of frames) each sensor reliability function based on the SNR (a measure of individual sensor reliability or confidence level) for each sensor during the tracking period.
However, under certain conditions (e.g., conditions causing a false alarm rate above a predetermined threshold), the fused (combined) reliability result determined (generated) by fusion processor 208 for the (entire) plurality of sensors (during the tracking period) may not be better than the individual sensor reliability result calculated from the performance of a better single sensor (e.g., the higher reliability sensor having the higher SNR). Therefore, the fusion processor 208 may use at least one additional predetermined sensor parameter (attribute) to better determine individual sensor reliability (function) weighting based on whether or not a data fusion result (generated from each sensor contributing) provides a more reliable result than a reliability result from a (better performing) single sensor.
Relying on predetermined measurements and analysis (e.g., testing and/or computer simulation of sensor operation), data fusion processor 208 may use the comparative (relative) received operating characteristics (ROC) between each sensor as the additional sensor parameter to help determine reliability weighting for each one of the plurality of sensors 205 during a predetermined tracking period. The ROC performance (curve) for each one of the plurality of sensors 205 may be generated (determined) using likelihood functions to represent (characterize) sensor information (during target tracking) such as (target) detections, no detections, measured SNRs, and other sensor information obtained from sensor measurements, observations, or other sensor data outputs. Thereafter, the ROC likelihood function for each sensor may be combined to generate likelihood (probability) functions of correct classification (recognition) of target and decoy (false target) for system 200.
For multi-sensor system 200, generation of the likelihood (probability) functions for correct classification (Pcc) of target and decoy using ROC likelihood function generation may include predetermination of the likelihood function for individual sensor noise caused by temporal fusion (diversity) as each sensor (auto) correlates data from multiple time frames (e.g., 120 time frames) during a predetermined tracking period. The temporal noise measurements (errors) for each one of the plurality of sensors 205 may be represented as a random variable (RV) where the negative impact of RV may be reduced using a plurality of methods including spatial and temporal fusion methods (used to combine data from differently located sensors and/or a single sensor outputting a plurality of data frames) to increase the probability of correct classification for a target and/or decoy (Pcc, Pct). Spatial and temporal fusion methods may be used to generate a combined likelihood (pdf) function for differently located sensors and/or sensors having a plurality of data time frames.
ROC (received operating characteristics) performance curves may be generated using a plurality of methods including calculation of the combined probability density function (pdf or likelihood function) for a plurality of different fusion methods. The plurality of different fusion methods may include, but are not limited to, additive fusion, multiplicative fusion (e.g., Bayes and/or Dempster-Shafer), fuzzy logic fusion using minimum and/or maximum calculations, and other fusion methods (strategies) that help to minimize the errors associated with noise (represented by RV). Likelihood function (pdf) calculations for each fusion method that are combined using this ROC method are shown in Appendix A. Each fusion method may be based on a two-object (e.g., target—t, decoy—d), spatial fusion example (e.g., IR and RF sensor) where the likelihood functions (representing Pcc) may be expressed as p(t1), p(d1) for a first sensor (sensor1—IR), and by p(t2), p(d2) for a second sensor (sensor2—RF).
Alternatively, ROC curves may be generated using computer simulations (calculations) to generate a high number of random samples (e.g., 10,000) to represent RVs with different pdfs. Thereafter, the combined pdfs may be determined from the histograms of combined RVs based on the different fusion methods (shown in Appendix A). Exemplary diagrams of ROC performance curves (generated using the alternative method) representing the probability of correct classification (versus probability of false alarm) for the plurality of sensors 205 of system 200 are shown in
a, 7b show exemplary curves of the probability of correct classification of decoy (Pcd—
As shown in
As shown in
Also, as shown in
Individual sensor reliability functions for each one of the plurality of sensors 205 may be expressed using the calculations and conditions shown in Appendix B. An optimal value (e.g., 0.5) for linear coefficient a may be determined from testing/simulation of weather (flight) conditions (of the sensor's field of view) for multi-sensor system 200. SNR may be dynamically calculated (estimated) from time frame to time frame for each one of the plurality of sensors 205 by measuring (calculating) the feature signal intensity for each frame, and dividing that value by the measured noise standard deviation for each frame. Alternatively, SNR may be calculated from a summation of the ignorance sets from a Dempster-Shafer computation as a measure of the noise intensity as disclosed in the cross-referenced provisional application Ser. No. 60/367,282, filed Mar. 26, 2002.
Data fusion processor 208 may use a plurality of fusion methods (algorithms) to generate relative and absolute reliability functions (levels) for multi-sensor system 200. The plurality of fusion methods may include a fusion method based on SNR (shown in Appendix C), a fusion method based on F/S ratio (shown in Appendix D), and a fusion method based on SNR and F/S ratio (shown in Appendix E). The methods may be based on a two-object (e.g., target—t, decoy—d), spatial fusion example (e.g., IR and RF sensor) where the likelihood functions (representing Pcc) may be expressed as p(t1), p(d1) for a first sensor (sensor1—IR), and by p(t2), p(d2) for a second sensor (sensor2—RF), and where the reliability for sensor1 at a particular time frame may be defined as rel1 and the reliability for sensor2 (at the same particular time frame) may be defined as rel2.
In accordance with embodiments of the present invention, (simulated) results of the additive data fusion algorithm (as shown in Appendices C, D, E) performed by multi-sensor system 200 are shown in
As shown in
a, 4b show diagrams of exemplary IR and RF sensor reliability functions (e.g., rel(t) as shown in Appendix B) for 300 time frames of decoy and target performance data, respectively. As shown in
a, 6b show the resulting diagrams of combined additive fusion (with adaptive, reliability weighting as described herein) using equations (7) and (8) from Appendix C. As shown in
A plurality of advantages may be provided in accordance with embodiments of the present invention including a data fusion method (incorporating additive fusion) that adaptively weights the contributions from different sensors (within a multi-sensor system) to generate at least two system reliability functions (relative and absolute reliability) where SNR and relative ROC performance between sensors may be used as measures of reliability.
Although the invention is primarily described herein using particular embodiments, it will be appreciated by those skilled in the art that modifications and changes may be made without departing from the spirit and scope of the present invention. As such, the method disclosed herein is not limited to what has been particularly shown and described herein, but rather the scope of the present invention is defined only by the appended claims.
Additive Fusion
p(t)=p(t1)+p(t2),
and
p(d)=p(d1)+p(d2)
For two independent random variables (RVs), X and Y, the combined pdf of the summation of these two RVs (Z=X+Y) may be calculated as the convolution of the two individual pdfs:
fz(z)=∫fx(x)fY(z−x)dx
(from 0 to ∞).
For an additive fusion example,
fp(t)(p(t))=∫fp(t1)(p(t1))fp(t2)(p(t)−p(t1)) dp(t1)
(from 0 to ∞),
and
fp(d)(p(d))=∫fp(d1)(p(d1)) fp(d2)(p(d)−p(d1)) dp(d1)
(from 0 to ∞),
The fused classification performance of the ROC curves may be estimated from the combined probability density functions (pdfs) in the above equations where the results are shown in
Multiplication (the Bayes) Fusion
p(t)=p(t1)*p(t2), and p(d)=p(d1)*p(d2).
The two independent RVs, X and Y, the combined pdf of the multiplication of these two RVs (Z=X*Y) may be calculated as the nonlinear convolution of the two individual pdfs:
fz(z)=·(1/|x|)Fx(x)fY(z/x) dx (from 0 to ∞).
For a multiplication fusion example,
fp(t)(p(t))=·1/|p(t1)|fp(t1)(p(t1)) fp(t2)(p(t)/p(t1)) (from 0 to ∞),
and
fp(d)(p(d))=·1/|p(d1)|fp(d1)(p(t1) fp(d2)(p(d))/p(d1))dp(d1) (from 0 to ∞).
The Relationship Between Additive and Multiplication Fusions
If the logarithm on both sides of the above-identified multiplication fusion equations is performed, then
In [p(t)]=In[p(t1)+In [p(t2)], and In[p(d)]=In[p(d1)+In[p(d2).
The one multiplication term becomes two additive terms of logarithm functions in each of the equation. If two RVs have log-normal pdfs, the equations above indicate that the multiplicative fusion of two RVs with log-normal distributions may be equivalent to the additive fusion of two RVs with normal distributions.
MIN, MAX, and MINMAX Fusions
The conjunction (AND) and disjunction (OR) are two well-known functions used in Fuzzy Logics. For two independent RVs: X and Y, the combined pdf of the conjunction of these two RVs [Z=min (X, Y)] may be given as:
fz(z)=fx(z)[1−FY(z)]+fY(z)[1−Fx(z)],
where F(z) is the cumulative distribution function.
Similarly, for two independent RVs: X and Y, the combined pdf of the disjunction of these two RVs [Z=max(X, Y)] may be given as:
Fx(z)=fx(z)FY(z)+fY(z)Fx(z).
For a two-object example, the MIN (conjunction) fusion may be defined as:
p(t)=min[p(t1), p(t2)], and p(d)=min [p(d1), p(d2)]
The MAX (disjunction) fusion may be defined as:
p(t)=max[p(t1), p(t2], and p(d)=max[p(d1), p(d2)].
The MINMAX fusion stategy may enhance one class (e.g., the target) over the other (e.g., the decoy). In some situations (e.g., target is the incoming missile warhead and the system is attempting to intercept the missile), the Pcc for the target should be as high as possible although the false alarm rate will be similarly increased. At the detection stage, the CFAR (constant false-alarm ratio) may be modified to change the Pd. At higher feature levels of fusion, the MINMAX fusion method may be used for this purpose.
For enhancing the target,
For enhancing the decoy,
For 0≦ rel (t)≦1, a reliability function, “rel (t)”, may be defined as a linear function of signal-to-noise ratio (SNR):
If rel2>rel1, then the relative reliability (rrel) may be expressed as:
For an exemplary scenario, if rel1=0.6 and rel2=0.8, then
Fr rel2>rel1, a combination of additive and multiplicative fusion may be expressed as:
Alternatively, if rel1>rel2, then
Also, for additive fusion, if rel2>rel1, then
Alternatively (for additive fusion), if rel1>rel2, then
p(t)=rrel2*p(t2)+p(t1), (9)
p(d)=rrel2*p(d2)+p(d1). (10)
As shown in equations (2), (3), when the relative reliability function (rrel1) −1, then the equation dissolves into the traditional multiplicative fusion.
Also, absolute reliability may be expressed as:
For a combination of additive and multiplicative fusion,
p(t)=rel1*rel2*[p(t1)*p(t2)]+(1−rel1)*p(t2)+(1−rel2)*p(t1), (11)
p(d)=rel1*rel2*[p(d1)*p(d2)]+(1−rel1)*p(d2)+(1−rel2)*p(d1), (12)
And for additive fusion,
P(t)=rel1*p(t1)+rel2*p(t2), (13)
P(d)=rel1* p(d1)+rel2*p(d2) (14)
Reliability Function Using F/S Ratio
For 0≦rl(t)≦1, a reliability function may be defined as a linear function of F/S ratio:
rl(t)={Rfs(fa,t), or 1 if Rfs(fa,t)>1 (15)
where t is the time frame number, fa is the false alarm rate, Rfs( ) is the F/S ratio, and
Rfs(fa,t)={Rfs(a,t), or 1 if fa<CFA (16)
Where CFA is the critical false alarm rate
For a combination of additive and multiplicative fusion at a specific time frame, and when comparing any two sensors of a plurality of sensors, if sensor2 is the better sensor (the sensor with better ROC performance), then
p(t)=rl*[p(t1)*p(t2)]+(1−rl)*p(t2), (17)
p(d)=rl*[p(d1)*p(d2)]+(1−rl)*p(d2). (18)
If sensor1 is the better sensor, then
p(t)=rl*[p(t1)*p(t2)](1−rl)*p(d2), (19)
p(d)=rl*[p(d1)*p(d2)]+(1−rl)*p(d1). (20)
For additive fusion, if sensor2 is the better sensor, then
p(t)=rl*p(t1)+p(t2), (21)
p(d)=rl*p(d1)+p(d2). (22)
If sensor1 is the better sensor, then
p(t)=p(t1)+rl*p(t2), (23)
p(d)=p(d1)+rl*p(d2). (24).
Reliability Function Using Both SNR and F/S Ratio
For additive fusion and using relative reliability for SNR, if rel2>rel1 and sensor1 is the better sensor, then
p(t)=rrel1*p(t1)+rl*p(t2), (25)
p(d)=rrel1*p(d1)+rl*p(d2), (26)
If sensor2 is the better sensor, then
p(t)=rl*rrel1*p(t1)+p(t2), (27)
p(d)=rl*rrel1*p(d1)+p(d2), (28)
If rel1>rel2 and sensor1 is the better sensor, then
p(t)=rl*p(t1)+rrel2*p(t2), (29)
p(d)=rl*p(d1)+rrel2*p(d2), (30)
For additive fusion and using absolute reliability for SNR, if sensor1 is the better sensor, then
p(t)=rel1*p(t1)+rl*rel2*p(t2), (31)
p(d)=rel1*p(d1)+rl*rel2*p(d2), (32)
If sensor2 is the better sensor, then
p(t)=rl*rel1*p(t1)+rel2*p(t2), (33)
p(d)=rl*rel1p(d1)+rel2*p(d2). (34)
This application claims the benefit of U.S. provisional application Ser. No. 60/367,282, filed Mar. 26, 2002.
Number | Name | Date | Kind |
---|---|---|---|
5341142 | Reis et al. | Aug 1994 | A |
6006163 | Lichtenwalner et al. | Dec 1999 | A |
6215983 | Dogan et al. | Apr 2001 | B1 |
6453272 | Slechta | Sep 2002 | B1 |
6487523 | Jarman et al. | Nov 2002 | B1 |
6512996 | Praskovsky et al. | Jan 2003 | B1 |
6559657 | McCarthy et al. | May 2003 | B1 |
6763322 | Potyrailo et al. | Jul 2004 | B1 |
6868346 | Larson et al. | Mar 2005 | B1 |
6909997 | Chen et al. | Jun 2005 | B1 |
6944566 | Chen et al. | Sep 2005 | B1 |
Number | Date | Country | |
---|---|---|---|
20030186663 A1 | Oct 2003 | US |
Number | Date | Country | |
---|---|---|---|
60367282 | Mar 2002 | US |