Patent Application
20100121797
This disclosure relates generally to chemical detection and more specifically to the detection of a chemical through signature extraction and classification.
The threat of attack on military and civilian targets employing chemical warfare agents and toxic industrial chemicals is of growing concern. Various technologies to detect and identify such chemicals are currently under development. Standoff chemical detectors are one example of a technology that can identify these chemicals. Standoff chemical detectors allow for real time, on the move detection for contamination avoidance and reconnaissance operations.
This disclosure provides a system and method for standoff chemical detection, including the detection of nitric acid.
In one embodiment, a method is disclosed that includes obtaining at least one measurement in a spectral domain of a sample and computing one or more qualities of the measurements in the spectral domain. The computing of the one or more qualities results in at least one peak within the spectral domain. This method also includes comparing the one or more computed qualities against a chemical signature of nitric acid. In addition, this method includes determining if the chemical signature is present in the sample.
In another embodiment, a system is disclosed for detecting a chemical signature that includes a processor for preprocessing an interferogram from received scene spectral information. The processor is configured to extract one or more features from the preprocessed interferogram corresponding to one or more predefined feature templates representative of one or more chemical vapor clouds and to process the features to determine if a chemical signature is present.
In yet another embodiment, a computer readable medium having instructions for causing a processor to perform a method detecting chemical is provided. The method includes receiving an interferogram which captures scene radiance information, performing apodization on the interferogram, and performing a chirp Fast Fourier Transform on the apodized interferogram. This method also includes applying a calibration curve and matching the signatures derived from the corrected spectrum to signatures derived from target chemical templates.
For a more complete understanding of this disclosure, reference is now made to the following description, taken in conjunction with the accompanying drawings, in which:
Nitric acid (HNO3) is a powerful oxidizing agent that is a highly corrosive and toxic strong acid. Nitric acid is used in a number of industrial and military applications. These military applications include explosives such as nitroglycerin and cyclotrimethylenetrinitramine.
The detection of nitric acid is difficult because of the reactivity of nitric acid, the false positives created by the presence of other compounds in a sample, and the interference caused by other compounds in a sample. These and other technical problems are overcome using feature signatures of nitric acid that has been obtained through measurements in the infrared (IR) spectral domain (hereinafter “spectral domain”). These feature signatures are created from the extracted features of the nitric acid as matched to other near extracted features of typical interferences. While the spectral domain described herein references the IR spectral domain, it is explicitly understood that any domain, including any domain within the electromagnetic spectrum, may be used consistent with this disclosure.
A feature signature refers to a set of one or more feature measurements, for example, amplitude and mse. Feature measurements quantify a chemical signature, which refers to the unique spectral characteristics that a particular compound will have at a given concentration within a frequency band in the spectral domain. A spectral characteristic can represent an emission peak or an absorption trough in the spectral band. Sometimes a spectral characteristic represents a partial peak or trough in a narrower spectral band. This is necessary when part of the peak or trough is affected by system artifacts or common background chemicals. The presence or partial presence of a peak or a trough in a spectral region can be referred to as the presence of a chemical signature. A chemical signature of a target chemical can also represent a spectral band where no peak or trough is present. This is referred to as absence of a chemical signature.
One of the complexities in detecting a chemical is that its spectrum changes based on the characteristics of the chemical cloud and its environment. That is under a specific condition (e.g., a cloud at long distance) the spectrum, thus the chemical, may be represented by only 1 peak (chemical signature). The same chemical, under a different condition, e.g., a cloud at short distance will have a different spectrum, e.g., two representative peaks (chemical signatures). As another example, a spectral peak may become broader and saturated when the concentration of the chemical cloud becomes very high. Hence, a chemical can have multiple sets of chemical signatures. The present disclosure detects a target chemical with multiple sets of chemical signatures.
Salient features are extracted from a chemical signature. One set of chemical signatures is quantified by a set of features. Consequently, a chemical can be represented by multiple sets of feature signatures, each of which corresponds to a set of salient features.
In order to detect the target chemical and the target chemical in the presence of other interfering chemicals while rejecting other interfering chemicals, the target chemical feature signature may be augmented with salient features of selected interfering chemicals. This creates an augmented feature signature of the target chemical. Thus, a target chemical is represented by multiple sets of augmented feature signatures. Each augmented feature signature is classified by the neural net. The target chemical is detected when any one of the augmented feature signatures is classified positively.
The present disclosure uses nitric acid as an example of a material that may be detected using the presently disclosed systems and methods. However, it is explicitly understood that any number of different compounds or elements may be detected using similar methods. Therefore, the present disclosure should not be limited to the detection of nitric acid.
A chemical detection system 100 for use in detecting chemicals is shown in
One type of chemical detection system utilized employs passive sensing of infrared (IR) emissions. The IR emissions along with background emissions are received through a window 132 mounted in the enclosure 110 and focused by a second lens system 136 onto a beam splitter element 140. Some of the IR is transmitted by a first stationary mirror 144 mounted behind the beam splitter element 140. The rest of the IR is reflected by splitter element 140 onto a moving mirror 146. The reflected beams from the stationary mirror 144 and moving mirror 146 combine to create an interference pattern, which is detected by an IR detector 148. An output of the IR detector 148 is sampled in one of two modes to create an interferogram, which is processed at a processor 160 to provide an output 170 such as a decision regarding whether or not the signatures of the chemical exists.
Therefore, when the target chemical is detected at the low resolution mode, the mode is switched (block 220) to a confirmation mode. The confirmation mode identifies the chemical at high resolution followed by sequential decision making (block 230). Sequential decision making determines if the feature signature is present by verifying it is repeatable. If the presence of the feature signature is confirmed, the chemical compound relating to the feature signature is mapped to provide an indication of the location of the chemical compound (block 240). Algorithms are utilized to detect chemicals as shown in
It is understood that the presently disclosed systems and methods can use the reduced resolution (16 wavenumber) “search mode” and then, upon the detection of part of the chemical signature, switch to a 4-wavenumber resolution “confirm mode”. The time to scan the entire field of range at 4-wavenumber resolution would exceed time constraints in many applications and would not provide enough time to take protective measures or to take evasive action for contamination avoidance. The time to acquire radiometrically equivalent 16-wavenumber resolution data is 16 times less than that for 4-wavenumber resolution data. The 16-wavenumber resolution data does not provide as much detail as the 4-wavenumber resolution data, hence the chemical differentiation and false alarm performances of the 16-wavenumber resolution mode can be poorer than that of the 4-wavenumber resolution mode. Therefore, a dual “search” and “confirmation” mode approach is used in one embodiment, in which the 16- and 4-wavenumber resolution modes are used in concert to meet timing and detection requirements. Of course, given faster measurement systems and processors, a single high-resolution mode approach will be feasible, or a single mode of suitable resolution may be used. The present disclosure is not limited to 4- or 16-wavenumber resolution.
Effectively, the search mode operation detects all regions of interest (ROI) that potentially have chemicals with the known chemical signatures. It may need to do this with a reasonably low rate of false triggers but with the same sensitivity as the confirmation mode because to miss a compound in search-mode would result in the failure to detect the compound. A rule is defined such that the search mode can be switched immediately to confirmation mode without scanning the entire field of range. This happens in the mode switch block 220 when the search mode result reaches a high confidence decision that a chemical cloud is present. Thus, the processing can detect the chemical in the shortest time. The confirmation mode applies a step and stare operation, in which high-resolution (4 cm-1) data is collected and analyzed to confirm the presence of and classify the types of compounds in the field of view. Any false triggers from the search-mode are rejected.
A further challenge is that the algorithms may need to detect down to very low signature strengths that approach the noise level of the system with a very low false alarm rate. The small signal detection capabilities are dictated by the concentration and size of the cloud 125, cloud distance and cloud-to-background temperature difference. Furthermore, the small chemical signal may need to be detected under many variations, which could be due to system-to-system differences or changes in operational environments. For example, the frequency of a laser diode that provides the data sampling reference in the sensor varies slightly from one laser to the next. As a result, the spectral resolution may vary from system to system. As another example, the detector response is affected by temperature, and consequently the spectral characteristics will be affected. Extracting the consistent feature signature amid the noise and signal variations may be critical to the success in the chemical detection.
The confirm mode utilizes a sequential decision process whereby a final detection decision is based on N-out-of-M detections from a sequence of confirm mode scans in the same field of view. When a sequential decision is invoked, the final decision at any instance of time can be one of three: “chemical detected,” “no chemical detected,” or “no final decision yet.” A final “chemical detected” is made only when strong evidence of the chemical is accumulated, such as a majority of the single decisions is consistent. On the other hand, a final decision on “no chemical detected” is made based on very weak or no evidence of chemical presence. Thus, any spurious, single scan, false detection will be rejected. In such cases, the detection cycle returns back to the previous stage. No final decision is made when the number of cumulative detected chemicals does not support nor deny the presence of a chemical. If no final decision is made, additional sequential scans are incorporated until the target chemical or no chemical decision is made. The process rules include an upper bound to the value of ‘M’ as a time constraint. Hence, sequential decision making reduces the false alarm rate and increases the confidence that a chemical is present when the final “chemical detected” decision is made.
Once sequential decision making confirms the presence of a chemical signature, the detection cycle switches to the chemical cloud mapping process (block 240). The chemical cloud mapping process locates the extents of the chemical cloud based on a search pattern.
The search and confirmation modes both process interferograms to make a decision on the presence and class of the chemical, if any. Both modes may utilize the same algorithm as shown in
Preprocessing 310, the first stage of the detection algorithm, transforms measured interferograms 410 into spectra as illustrated in
The preprocessing stage also includes the following functions as shown in
A chirp-Fast Fourier Transform (FFT) converts spectra to an identical frequency comb of 4-wavenumbers for all systems at 430. Each sensor may have a different sampling reference. The chirp-FFT allows sampling of data at selected frequencies and interpolates to a selected frequency comb to calibrate between the sensors.
Frequency dependent gain/offset correction is applied at 440 to provide a spectrum comprising amplitude for each wavenumber at 450. The gain/offset correction is derived from a calibration process. In
The feature extraction stage is shown in
In one embodiment, the spectrum output 520 first undergoes a normalization process, wherein the spectrum output 520 is divided by a Planck's function whose temperature is estimated from several points of the spectrum output 520. The output is a normalized spectrum, which has peaks and valleys around nominal values of one.
The feature extractor is designed to be sensitive to peaks and valleys in the spectrum. When the system is aimed at a blackbody scene, all elements in the resulting feature vector are zero except for noise. Warm chemical clouds relative to the scene produce emission peaks in the spectrum and corresponding positive amplitude feature vector elements. Cool chemical clouds produce absorption valleys and negative amplitude feature elements. Each feature measurement in the feature signature is extracted from the spectrum by a match filter that has been tailored to a particular peak in the absorption coefficient curve of a chemical or common interferent. Thus, each feature in the feature signature characterizes the peak or valley in the scene spectrum at a frequency band with a shape that corresponds to a known chemical/interferent absorption phenomenon.
The match filters utilized by the feature extractor 530 are selected using a heuristic approach with the objective of maximizing detection sensitivity and discriminating capability. An initial set of potential match filters is derived from the target chemical and interferent absorption coefficient curves in the frequency range of interest. This initial set consists of several hundred potential match filters. The most prominent match filters from each target chemical are chosen since these provide the greatest detection sensitivity relative to the noise in the system. Some prominent interferent match filters are also chosen because these can sometimes provide discriminating capability. Typically, a set of 20 to 40 match filters are selected and packaged into a feature matrix 510 that can be loaded into the system via an interface 515.
The feature extractor 530 produces a discriminating feature signature 540. For the subset of scenes that are relatively simple—a target chemical or interferent cloud against a relatively benign background, a chemical decision can be made based on a threshold on the feature signature. For more complex scenes, with multiple target chemicals and/or interferents and feature-rich backgrounds, a further classifier stage described below is utilized.
Template 645 illustrates matching of template 620 to detected spectra. Curves 640, 641, 642, and 643 are shown superimposed on the graph with spectral band from the normalized spectra 650. A least squares fit algorithm is applied to determine matches. The fit algorithm computes amplitude, slope, offset, mean square error of fit (mse), and skew of the fit between the template and the spectral region. Given a shape template, S, whose first and second moment are zero, and the corresponding spectral region, Y, (both Y and S are vectors of length n), the amplitude, slope, offset, and mse are computed as follow:
amplitude=Y.S′
slope=Y′.L
offset=mean(Y)
mse=square_root(((P(i)−Y(i))2)/n)
where:
L=(L0−mean(L0))/norm(L0−mean(L0);
P=offset*U+slope*L+amplitude*S/S(i)2; and
The third stage of the detection algorithm is the classification algorithm 700, as shown in
The feature signatures is represented at 705 in
The first process of the classifier is a preconditioning step, where the classifier performs a normalization step process 730 . . . 760 to be able to detect or classify a wider range of chemical signatures. In one embodiment, the normalization step is an option determined by a parameter in the classifier parameters 710. Also involved is a noise threshold test 720 . . . 755, which measures and removes very weak signals. The measure is a weighed sum of features that are predefined for each target chemical classifier, and is compared with a threshold. This threshold is adaptively set according to minimum detection requirements for each target chemical, the false alarm requirements and the SNR for the system in operation. When the measure does not exceed the threshold, a weak signal and no-chemical detection for that target chemical classifier is declared without exercising that feature signature neural network.
The classification algorithm may be implemented through a neural network bank 740 . . . 770, in which each of the neural networks is trained to detect a particular feature signature corresponding to a target chemical and reject other non-similar target chemicals, different interferents and background signatures. The neural network is based on the backpropagation architecture with one hidden layer. The size of the hidden layer was carefully chosen in order to classify the target chemical under different scenarios and not over generalize the detection scheme. An output threshold 780 is associated with each neural network that is tuned based on detection performance and false alarm rate. Since there are usually multiple templates per target chemical deriving the key discriminating features for that target chemical, not all of the feature measurements computed in the feature extraction process need be run through the neural network for it to arrive at the target chemical decision. The selected feature indices for each classifier are stored in the classifier parameters 710. While a neural network is shown in
It is understood that
While this disclosure has described certain embodiments and generally associated methods, alterations and permutations of these embodiments and methods will be apparent to those skilled in the art. Accordingly, the above description of example embodiments does not define or constrain this disclosure. Other changes, substitutions, and alterations are also possible without departing from the spirit and scope of this disclosure, as defined by the following claims.
The invention described herein was made with U.S. Government support under subcontract number LS97-00001 under Prime Contract Number DAAM01-97-C-0030 awarded by U.S. Army, SBCCOM, Edgewood, Md. The United States Government has certain rights in the invention.
