US6292592B1 - Efficient multi-resolution space-time adaptive processor - Google Patents
Efficient multi-resolution space-time adaptive processor Download PDFInfo
- Publication number
- US6292592B1 US6292592B1 US09/175,222 US17522298A US6292592B1 US 6292592 B1 US6292592 B1 US 6292592B1 US 17522298 A US17522298 A US 17522298A US 6292592 B1 US6292592 B1 US 6292592B1
- Authority
- US
- United States
- Prior art keywords
- matrix
- providing
- tap weights
- coefficients
- response
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/295—Means for transforming co-ordinates or for evaluating data, e.g. using computers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/78—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using electromagnetic waves other than radio waves
- G01S3/782—Systems for determining direction or deviation from predetermined direction
- G01S3/785—Systems for determining direction or deviation from predetermined direction using adjustment of orientation of directivity characteristics of a detector or detector system to give a desired condition of signal derived from that detector or detector system
- G01S3/786—Systems for determining direction or deviation from predetermined direction using adjustment of orientation of directivity characteristics of a detector or detector system to give a desired condition of signal derived from that detector or detector system the desired condition being maintained automatically
- G01S3/7864—T.V. type tracking systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/52—Discriminating between fixed and moving objects or between objects moving at different speeds
- G01S13/522—Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves
- G01S13/524—Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves based upon the phase or frequency shift resulting from movement of objects, with reference to the transmitted signals, e.g. coherent MTi
- G01S13/5246—Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves based upon the phase or frequency shift resulting from movement of objects, with reference to the transmitted signals, e.g. coherent MTi post processors for coherent MTI discriminators, e.g. residue cancellers, CFAR after Doppler filters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/66—Radar-tracking systems; Analogous systems
- G01S13/72—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar
- G01S13/723—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar by using numerical data
- G01S13/726—Multiple target tracking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
- G01S17/89—Lidar systems specially adapted for specific applications for mapping or imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/483—Details of pulse systems
- G01S7/486—Receivers
- G01S7/487—Extracting wanted echo signals, e.g. pulse detection
Definitions
- the present invention relates to radar and electro-optic image data processing systems and techniques. More specifically, the present invention relates to space time adaptive array processing systems and techniques for use in radar and electro-optic image data processing applications.
- Active and passive imaging technologies are employed in diverse applications where there is a need to track an object as it moves through space.
- radar is often used to track targets and other aircraft.
- Electro-optic technologies including laser based systems are also used for such applications.
- Clutter is often due to the detection of objects other than a desired target and may result from natural as well as artificial objects. Further, the clutter may vary in size and number and may be static or dynamic. Interference may be intentional and, if so, it may be designed to overpower the sensing technology or it may be designed to cause a malfunction or misread of a true target location.
- STAP space time adaptive array processing
- a first composite input signal is provided based on plurality of data values output from a sensor in response to a scene including a target and clutter.
- the composite signal is processed to provide a plurality of tap weights.
- the tap weights are generated by the matrix of data values which is first filtered by a wavelet transform to provide a set of coefficients.
- the coefficients are sparsened to provide a sparse matrix.
- the sparse matrix is then inverse wavelet transformed to provide the tap weights.
- the tap weights are applied to the composite signal to yield a clutter reduced output signal.
- the matrix is a covariance matrix.
- a method for implementing the teachings of the invention in the data domain is also disclosed.
- the sparsed matrix is inverted and a set of steering vectors is applied to create the tap weights.
- the invention affords an enhanced Signal-to-Interference+Noise Ratio (SINR) because (i) wavelets provide better bases for nonstationary processes and therefore offer improved sample support performance and (ii) coefficient thresholding in wavelet domain removes noisy data that is difficult to estimate.
- SINR Signal-to-Interference+Noise Ratio
- FIG. 1 shows an Electronic Counter Countermeasure (ECCM) scenario that may occur in advanced air-to-air encounters.
- ECCM Electronic Counter Countermeasure
- FIG. 2 a shows how a notch can be placed at the location of the jammer in the main antenna beam and indicates how the effects of the jammer can be canceled.
- FIG. 2 b illustrates how aperture-limited resolution provided by conventional nonadaptive processing is enhanced to provide subbeamwidth resolution that permits the separation of the two closely spaced targets.
- FIG. 3 illustrates range doppler spread in a terrain scattered interference jamming scenario.
- FIG. 4 is a block diagram illustrative of the efficient multi-resolution space time adaptive processing system and method of the present invention.
- FIG. 5 illustrates a pyramidal algorithm
- FIG. 6 illustrates an extension of a one-dimensional pyramidal algorithm to two dimensions.
- FIG. 7 represents a wavelet coefficient mapping of a decomposed matrix and clearly shows the multiple resolutions.
- FIG. 8 is a diagram illustrating a sparsening of a covariance matrix using a discrete wavelet transform in accordance with the method of the present invention.
- FIG. 9 a shows a sample covariance matrix
- FIG. 9 b show a discrete wavelet transform of the sample covariance matrix of FIG. 9 a.
- FIG. 10 depicts SNR as function of degree of sparsening of covariance matrix.
- FIG. 11 depicts a 2-D histogram of significant wavelet coefficients.
- FIGS. 12 a and 12 b depict SNR vs. sparsening using a coefficient template.
- FIG. 1 shows an Electronic Counter Countermeasure (ECCM) scenario 10 that may occur in advanced air-to-air encounters.
- ECCM Electronic Counter Countermeasure
- Two particular jamming possibilities are illustrated: (i) a Stand Off Jammer (SOJ) 12 or 14 that places a jamming signal in the mainlobe 16 of an antenna beam from a missile 20 ; and (ii) two closely spaced targets 22 and 24 residing within the mainlobe 16 .
- SOJ Stand Off Jammer
- FIG. 2 a shows how a notch can be placed at the location of the jammer in the main antenna beam and indicates how the effects of the jammer can be canceled.
- Case (ii) which is an example of superresolution, e.g., the MUSIC algorithm, is shown in FIG. 2 b and illustrates how the aperture-limited resolution provided by conventional nonadaptive processing is enhanced to provide subbeamwidth resolution that permits the separation of the two closely spaced targets.
- FIG. 3 illustrates range doppler spread in a terrain scattered interference jamming scenario.
- the Terrain Scattered Interference (TSI) has a significant spread in range and doppler.
- STAP Space-Time-Doppler Adaptive Processing
- STAP Space-Time-Doppler Adaptive Processing
- the faithful representation of the covariance matrix requires a large number of degrees of freedom. It may also be possible to handle the highly nonstationary clutter environment by using fewer degrees of freedom, but processing at a much faster update rate. Either of the scenarios require multi-GFLOP processing.
- the purpose of adaptive processing is to remove the interference to extract obscured targets.
- STAP works by estimating the covariance matrices from sample data that does not contain the target.
- FIG. 4 is a block diagram illustrative of the efficient multi-resolution space time adaptive processing system and method of the present invention.
- the system 30 includes a sensor 32 and a STAP processor 40 .
- the sensor 32 is typically electromagnetic, e.g., radar, or electro-optic, e.g., laser.
- the STAP processor 40 is implemented in accordance with the present teachings preferably in an application specific integrated circuit (ASIC) or a field programmable gate array (FPGA).
- ASIC application specific integrated circuit
- FPGA field programmable gate array
- the STAP processor 40 is implemented within a digital processing system (not shown) provided between the sensor 32 and a system (not shown) for providing track processing, control or display as will be appreciated by those skilled in the art.
- the sensor 32 outputs an electrical data signal that represents a target and interference.
- the sensor output is input to a filter bank 33 and a covariance estimator 34 .
- a set of weight vectors, W n is output by the estimator 34 from a sample covariance matrix, ⁇ , and a set of steering vectors, S n .
- the steering vectors S n is supplied by a cueing system 36 as is common in the art.
- the tap weights are used by the filter bank to filter the clutter from the sensor output in a conventional manner.
- the system 40 solves, in real time, the matrix equation:
- this is a numerically intensive process that drives the STAP throughput requirements.
- the present invention reduces the computational requirements of STAP processing.
- the present invention treats the matrices that arise in RF STAP—Range Doppler Maps, Data Matrices and Covariance Matrices—as images and then exploit techniques that have been developed in multiresolution image analysis.
- Specific properties of the wavelet transform that are useful include: sparsening, denoising, edge detection and efficient bases for nonstationary processes.
- a brief introduction to wavelets and wavelet packets is provided below to facilitate an understanding of the present invention. The reader is referred to a number of excellent papers for more detailed information. (See for example, “Ten Lectures on Wavelets,” by I. Daubeshies, SIAM, Philadelphia, Pa., 1992 and “A theory for Multiresolution Signal Decomposition: The Wavelet Representation,” by S. Mallat, IEEE Trans., Pattern Anal. Mach. Intel., Vol. 11, pp. 674-693, 1989.)
- the wavelet transform is characterized by a dilation parameter and a translation parameter.
- Continuous wavelet transforms have the properties that they possess a decomposition of the identity and are overcomplete. For this reason, they are not efficient computationally and thus have not found widespread engineering applications.
- FIG. 5 illustrates the pyramidal algorithm introduced by Mallat.
- a Quadrature Mirror Filter pair, g(n) and h(n) is first constructed.
- the data vector V n 0 is then processed into a smoothed portion, V n 1 , and a details portion, W n 1 .
- the first octave smoothed vector is again decomposed into the second octave smoothed and detail data vector.
- the process is repeated for each smoothed vector; the result is an (n) orthogonal subspace decomposition of the data vector with no redundancy.
- the decomposition of the vector into a series of orthogonal vector spaces allows each octave to be processed independently. This so-called multi-resolution decomposition is the basis for many proposed applications.
- FIG. 6 illustrates an extension of a one-dimensional pyramidal algorithm to two dimensions.
- the 2-D DWT is used to decompose and reconstruct a two dimensional data vector, i.e., a matrix.
- FIG. 7 represents a wavelet coefficient mapping of a decomposed matrix and clearly shows the multiple resolutions. This multi-resolution property may provide a useful approach to target discrimination based on feature extraction since many target characteristics appear differently at each resolution level.
- the discrete wavelet packet is a generalization of the DWT in that both the details and smoothed portions of the data vector are filtered at each octave.
- the DWP tree is then pruned and a ‘best basis’ is selected based on entropy or similar criteria.
- FIG. 8 is a diagram illustrating a sparsening of a covariance matrix using a discrete wavelet transform in accordance with the method of the present invention.
- the present invention uses the wavelet transform to solve Eq. [1] in a computationally efficient manner.
- the sample covariance matrix ⁇ is first transformed using a Discrete Wavelet Transform (DWT) into ⁇ W .
- DWT Discrete Wavelet Transform
- the transformed matrix ⁇ W is sparsed.
- small coefficients are thresholded to zero, resulting in a sparsened matrix ⁇ circumflex over ( ⁇ ) ⁇ W .
- the sparsened matrix ⁇ circumflex over ( ⁇ ) ⁇ W is then inverted to yield ⁇ circumflex over ( ⁇ ) ⁇ W ⁇ 1 and transformed using an Inverse DWT (IDWT) to yield ⁇ circumflex over ( ⁇ ) ⁇ ⁇ 1 .
- This value ⁇ circumflex over ( ⁇ ) ⁇ ⁇ 1 is then inserted into equation [1] above to calculate the tap weights.
- thresholding techniques including global thresholding, zero trees, and vector quantization See “Image Coding Using Wavelet Transform,” by Antonini, M. Barlaud, P. Mathieu, and I. Daubechies, IEEE, Trans. Image Processing, Vol. 1, pp. 205-220, 1992 and “Image Coding Using Vector Quantization: A Review,” by Nasrabadi and R. King, IEEE, Trans. Commun., Vol. 36, August 1988.
- FIG. 9 a shows a sample covariance matrix.
- FIG. 9 b show a discrete wavelet transform of the sample covariance matrix of FIG. 9 a .
- the physical environment is representative of TSI and was generated from a high fidelity simulation.
- a banding structure, indicative of a sparse matrix, is clearly present in the transformed sample covariance matrix. This banding structure is known to be a rather general property of many types of scattering and does not depend on the specific details of the TSI model.
- FIG. 10 depicts SNR as function of degree of sparsening of covariance matrix.
- An appropriate metric for evaluating the wavelet transform is the SNR as a function of sparsening of the sample covariance matrix. From FIG. 10, it can be seen that a sparsening of approximately 90% can be achieved without materially degrading the SNR. Interestingly, a slight increase in SNR is evident for some sparsening; there are several interpretations of this: de-noising, improved training data, etc.
- the ability to achieve significant sparsening is an indication that the DWT is in some sense a close approximation to a data independent Karhunen-Loeve Transform.
- FIG. 11 depicts a 2-D histogram of significant wavelet coefficients. What is plotted here is a histogram of significant wavelet coefficients across 128 Pulse Repetition Intervals (PRIs). The majority of the histogram values were either 0 or 128, indicating that the coefficients were either significant or insignificant for all PRIs. This was further borne out (c.f., FIG. 12) by calculating the SNR as a function of sparsening by forming a template on the first PRI and using that for all remaining PRIs.
- PRIs Pulse Repetition Intervals
- FIGS. 12 a and 12 b depict SNR vs. sparsening using a coefficient template.
- A is an N ⁇ M data matrix
- R is an M ⁇ M upper triangular matrix
- Q * ⁇ ⁇ ⁇ ⁇ Q [ Q 0 * 0 0 0 0 Q 0 * 0 0 0 0 Q 0 * 0 0 0 Q 0 * ] ⁇ [ ⁇ 0 ⁇ - 1 ⁇ - 2 ⁇ - 3 ⁇ 1 0 ⁇ - 1 ⁇ - 2 ⁇ 2 ⁇ 1 ⁇ 0 ⁇ - 1 ⁇ 3 ⁇ 2 ⁇ 1 ⁇ 0 ]
Abstract
An image processing system and method. In accordance with the inventive method, adapted for use in an illustrative image processing application, a first composite input signal is provided based on plurality of data values output from a sensor in response to a scene including a target and clutter. The composite signal is processed to provide a plurality of tap weights. The tap weights are generated by the matrix of data values which is first filtered by a wavelet transform to provide a set of coefficients. The coefficients are sparsened to provide a sparse matrix. The sparse matrix is then inverse wavelet transformed to provide the tap weights. Finally, the tap weights are applied to the composite signal to yield a clutter reduced output signal. In the illustrative implementation, the matrix is a covariance matrix. However, a method for implementing the teachings of the invention in the data domain is also disclosed. In the illustrative implementation, the sparsed matrix is inverted and a set of steering vectors is applied to create the tap weights. The invention affords an enhanced Signal-to-Interference+Noise Ratio (SINR) because (i) wavelets provide better bases for nonstationary processes and therefore offer improved sample support performance and (ii) coefficient thresholding in wavelet domain removes noisy data that is difficult to estimate.
Description
1. Field of the Invention
The present invention relates to radar and electro-optic image data processing systems and techniques. More specifically, the present invention relates to space time adaptive array processing systems and techniques for use in radar and electro-optic image data processing applications.
2. Description of the Related Art
Active and passive imaging technologies are employed in diverse applications where there is a need to track an object as it moves through space. In military and commercial applications, for example, radar is often used to track targets and other aircraft. Electro-optic technologies including laser based systems are also used for such applications.
Unfortunately, as is well known in the art, tracking of targets including aircraft and spacecraft using radar and electro-optic techniques and the tracking of vessels using sonar may be problematic due to the presence of clutter and other sources of interference. Clutter is often due to the detection of objects other than a desired target and may result from natural as well as artificial objects. Further, the clutter may vary in size and number and may be static or dynamic. Interference may be intentional and, if so, it may be designed to overpower the sensing technology or it may be designed to cause a malfunction or misread of a true target location.
Accordingly, the elimination of clutter and other sources of interference has received considerable attention from designers of passive and active tracking systems. One technique which offers the promise of being of use for such applications is called space time adaptive array processing or “STAP”. Considered for use with systems that employ an array of sensing elements, STAP would involve the creation of a covariance matrix in the vicinity of the target. The covariance matrix would be used to provide an estimate of the clutter. The estimate would then be used to remove the clutter in a gate around the target.
STAP attempts to suppress spatio-temporal interference, hence covariance matrix estimates must be updated in real time to handle rapidly changing interference statistics. Unfortunately, clutter is often nonstationary due to movement or jamming. As a result, the steps of estimating the covariance matrix and canceling the clutter using STAP would be computationally intensive. Hence, the data processing requirements for current applications are often considerable, e.g., on the order of ten billion floating point operations per second (10 GFLOPS).
Accordingly, there is a need in the art for an efficient system and method for processing sensor outputs to eliminate clutter and interference for current and future military, commercial and industrial applications.
The need in the art is addressed by the image processing system and method of the present invention. In accordance with the inventive method, adapted for use in an illustrative image processing application, a first composite input signal is provided based on plurality of data values output from a sensor in response to a scene including a target and clutter. The composite signal is processed to provide a plurality of tap weights. The tap weights are generated by the matrix of data values which is first filtered by a wavelet transform to provide a set of coefficients. The coefficients are sparsened to provide a sparse matrix. The sparse matrix is then inverse wavelet transformed to provide the tap weights. Finally, the tap weights are applied to the composite signal to yield a clutter reduced output signal.
In the illustrative implementation, the matrix is a covariance matrix. However, a method for implementing the teachings of the invention in the data domain is also disclosed. In the illustrative implementation, the sparsed matrix is inverted and a set of steering vectors is applied to create the tap weights.
The invention affords an enhanced Signal-to-Interference+Noise Ratio (SINR) because (i) wavelets provide better bases for nonstationary processes and therefore offer improved sample support performance and (ii) coefficient thresholding in wavelet domain removes noisy data that is difficult to estimate.
FIG. 1 shows an Electronic Counter Countermeasure (ECCM) scenario that may occur in advanced air-to-air encounters.
FIG. 2a shows how a notch can be placed at the location of the jammer in the main antenna beam and indicates how the effects of the jammer can be canceled.
FIG. 2b illustrates how aperture-limited resolution provided by conventional nonadaptive processing is enhanced to provide subbeamwidth resolution that permits the separation of the two closely spaced targets.
FIG. 3 illustrates range doppler spread in a terrain scattered interference jamming scenario.
FIG. 4 is a block diagram illustrative of the efficient multi-resolution space time adaptive processing system and method of the present invention.
FIG. 5 illustrates a pyramidal algorithm.
FIG. 6 illustrates an extension of a one-dimensional pyramidal algorithm to two dimensions.
FIG. 7 represents a wavelet coefficient mapping of a decomposed matrix and clearly shows the multiple resolutions.
FIG. 8 is a diagram illustrating a sparsening of a covariance matrix using a discrete wavelet transform in accordance with the method of the present invention.
FIG. 9a shows a sample covariance matrix.
FIG. 9b show a discrete wavelet transform of the sample covariance matrix of FIG. 9a.
FIG. 10 depicts SNR as function of degree of sparsening of covariance matrix.
FIG. 11 depicts a 2-D histogram of significant wavelet coefficients.
FIGS. 12a and 12 b depict SNR vs. sparsening using a coefficient template.
Illustrative embodiments and exemplary applications will now be described with reference to the accompanying drawings to disclose the advantageous teachings of the present invention.
While the present invention is described herein with reference to illustrative embodiments for particular applications, it should be understood that the invention is not limited thereto. Those having ordinary skill in the art and access to the teachings provided herein will recognize additional modifications, applications, and embodiments within the scope thereof and additional fields in which the present invention would be of significant utility.
FIG. 1 shows an Electronic Counter Countermeasure (ECCM) scenario 10 that may occur in advanced air-to-air encounters. Two particular jamming possibilities are illustrated: (i) a Stand Off Jammer (SOJ) 12 or 14 that places a jamming signal in the mainlobe 16 of an antenna beam from a missile 20; and (ii) two closely spaced targets 22 and 24 residing within the mainlobe 16. Both of these threats can be addressed with adaptive processing techniques. Case (i) is illustrated in FIG. 2a which shows how a notch can be placed at the location of the jammer in the main antenna beam and indicates how the effects of the jammer can be canceled. Case (ii), which is an example of superresolution, e.g., the MUSIC algorithm, is shown in FIG. 2b and illustrates how the aperture-limited resolution provided by conventional nonadaptive processing is enhanced to provide subbeamwidth resolution that permits the separation of the two closely spaced targets.
FIG. 3 illustrates range doppler spread in a terrain scattered interference jamming scenario. In the typical jamming scenario illustrated in FIG. 3, the Terrain Scattered Interference (TSI) has a significant spread in range and doppler. The extraction of targets from nonstationary clutter may thus require the use of Space-Time-Doppler Adaptive Processing (STAP). In this case, the faithful representation of the covariance matrix requires a large number of degrees of freedom. It may also be possible to handle the highly nonstationary clutter environment by using fewer degrees of freedom, but processing at a much faster update rate. Either of the scenarios require multi-GFLOP processing.
The purpose of adaptive processing is to remove the interference to extract obscured targets. STAP works by estimating the covariance matrices from sample data that does not contain the target.
FIG. 4 is a block diagram illustrative of the efficient multi-resolution space time adaptive processing system and method of the present invention. The system 30 includes a sensor 32 and a STAP processor 40. The sensor 32 is typically electromagnetic, e.g., radar, or electro-optic, e.g., laser. The STAP processor 40 is implemented in accordance with the present teachings preferably in an application specific integrated circuit (ASIC) or a field programmable gate array (FPGA). The STAP processor 40 is implemented within a digital processing system (not shown) provided between the sensor 32 and a system (not shown) for providing track processing, control or display as will be appreciated by those skilled in the art. The sensor 32 outputs an electrical data signal that represents a target and interference. The sensor output is input to a filter bank 33 and a covariance estimator 34. A set of weight vectors, Wn, is output by the estimator 34 from a sample covariance matrix, Λ, and a set of steering vectors, Sn. The steering vectors Sn is supplied by a cueing system 36 as is common in the art. The tap weights are used by the filter bank to filter the clutter from the sensor output in a conventional manner. The system 40 solves, in real time, the matrix equation:
In accordance with conventional teachings, this is a numerically intensive process that drives the STAP throughput requirements. However, as set forth below, the present invention reduces the computational requirements of STAP processing.
STAP Processing Using Wavelet Transforms
The present invention treats the matrices that arise in RF STAP—Range Doppler Maps, Data Matrices and Covariance Matrices—as images and then exploit techniques that have been developed in multiresolution image analysis. Specific properties of the wavelet transform that are useful include: sparsening, denoising, edge detection and efficient bases for nonstationary processes. A brief introduction to wavelets and wavelet packets is provided below to facilitate an understanding of the present invention. The reader is referred to a number of excellent papers for more detailed information. (See for example, “Ten Lectures on Wavelets,” by I. Daubeshies, SIAM, Philadelphia, Pa., 1992 and “A theory for Multiresolution Signal Decomposition: The Wavelet Representation,” by S. Mallat, IEEE Trans., Pattern Anal. Mach. Intel., Vol. 11, pp. 674-693, 1989.)
Introduction to Wavelets and Wavelet Packets
The wavelet transform is characterized by a dilation parameter and a translation parameter. Continuous wavelet transforms have the properties that they possess a decomposition of the identity and are overcomplete. For this reason, they are not efficient computationally and thus have not found widespread engineering applications.
In 1989, Mallat introduced a pyramidal algorithm for the discrete wavelet transform (DWT) that, like the FFT, is a fast, linear operation that operates on a data vector.
FIG. 5 illustrates the pyramidal algorithm introduced by Mallat. A Quadrature Mirror Filter pair, g(n) and h(n), is first constructed. The data vector Vn 0 is then processed into a smoothed portion, Vn 1, and a details portion, Wn 1. The first octave smoothed vector is again decomposed into the second octave smoothed and detail data vector. The process is repeated for each smoothed vector; the result is an (n) orthogonal subspace decomposition of the data vector with no redundancy. The decomposition of the vector into a series of orthogonal vector spaces allows each octave to be processed independently. This so-called multi-resolution decomposition is the basis for many proposed applications.
FIG. 6 illustrates an extension of a one-dimensional pyramidal algorithm to two dimensions. The 2-D DWT shown here is separable; it is the tensor product of two 1-D DWTs: T(x,y)=T(x){circle around (x)}T(y). The 2-D DWT is used to decompose and reconstruct a two dimensional data vector, i.e., a matrix.
FIG. 7 represents a wavelet coefficient mapping of a decomposed matrix and clearly shows the multiple resolutions. This multi-resolution property may provide a useful approach to target discrimination based on feature extraction since many target characteristics appear differently at each resolution level.
The discrete wavelet packet (DWP) is a generalization of the DWT in that both the details and smoothed portions of the data vector are filtered at each octave. The DWP tree is then pruned and a ‘best basis’ is selected based on entropy or similar criteria.
While both the FFT and DWT map the data vector from one vector space into another, their basis functions are significantly different. The basis function for the FFT are complex exponentials, which are localized in frequency but not in time; the basis functions for the DWT are much more general and are localized in both frequency and time. It is this dual localization that renders large classes of functions sparse in the wavelet domain.
STAP Implementation using the DWT
FIG. 8 is a diagram illustrating a sparsening of a covariance matrix using a discrete wavelet transform in accordance with the method of the present invention. The present invention uses the wavelet transform to solve Eq. [1] in a computationally efficient manner. As shown schematically in FIG. 8, the sample covariance matrix Λ is first transformed using a Discrete Wavelet Transform (DWT) into ΛW. Next, the transformed matrix ΛW is sparsed. During the sparsing operation, small coefficients are thresholded to zero, resulting in a sparsened matrix {circumflex over (Λ)}W. The sparsened matrix {circumflex over (Λ)}W is then inverted to yield {circumflex over (Λ)}W −1 and transformed using an Inverse DWT (IDWT) to yield {circumflex over (Λ)}−1. This value {circumflex over (Λ)}−1 is then inserted into equation [1] above to calculate the tap weights. Those skilled in the art will appreciate that the invention is not limited to the thresholding technique employed in the sparsening operation. A number of thresholding techniques can be employed, including global thresholding, zero trees, and vector quantization See “Image Coding Using Wavelet Transform,” by Antonini, M. Barlaud, P. Mathieu, and I. Daubechies, IEEE, Trans. Image Processing, Vol. 1, pp. 205-220, 1992 and “Image Coding Using Vector Quantization: A Review,” by Nasrabadi and R. King, IEEE, Trans. Commun., Vol. 36, August 1988.
FIG. 9a shows a sample covariance matrix. FIG. 9b show a discrete wavelet transform of the sample covariance matrix of FIG. 9a. The physical environment is representative of TSI and was generated from a high fidelity simulation. A banding structure, indicative of a sparse matrix, is clearly present in the transformed sample covariance matrix. This banding structure is known to be a rather general property of many types of scattering and does not depend on the specific details of the TSI model.
While significant sparsening can clearly be achieved, there are two major concerns with working in the DWT domain: performance degradation and transform overhead. These issues are addressed below.
Performance Evaluation
FIG. 10 depicts SNR as function of degree of sparsening of covariance matrix. An appropriate metric for evaluating the wavelet transform is the SNR as a function of sparsening of the sample covariance matrix. From FIG. 10, it can be seen that a sparsening of approximately 90% can be achieved without materially degrading the SNR. Interestingly, a slight increase in SNR is evident for some sparsening; there are several interpretations of this: de-noising, improved training data, etc. The ability to achieve significant sparsening is an indication that the DWT is in some sense a close approximation to a data independent Karhunen-Loeve Transform.
Temporal Stability of the DWT Coefficients and Data Domain Processing
In this section, the difficult problem of reducing the overhead required when constructing and implementing the DWT is considered. Two techniques are examined separately: temporal stability and data domain processing.
Temporal Stability
This section discusses the investigation into temporal stability of the wavelet coefficients. Those skilled in the art will appreciate that although the magnitude of the significant wavelet coefficients could change rapidly in time, their location might not. It turns out that this is indeed the case, as illustrated in FIG. 11.
FIG. 11 depicts a 2-D histogram of significant wavelet coefficients. What is plotted here is a histogram of significant wavelet coefficients across 128 Pulse Repetition Intervals (PRIs). The majority of the histogram values were either 0 or 128, indicating that the coefficients were either significant or insignificant for all PRIs. This was further borne out (c.f., FIG. 12) by calculating the SNR as a function of sparsening by forming a template on the first PRI and using that for all remaining PRIs.
FIGS. 12a and 12 b depict SNR vs. sparsening using a coefficient template.
Data Domain Processing
In practice, STAP calculations may be performed in the data (voltage) domain using the M×M sample covariance matrix, Λ, computed as
where A is an N×M data matrix.
Using QR-decomposition:
where R is an M×M upper triangular matrix.
Similarly, one can factor the covariance matrix using Cholesky decomposition; to within a scale factor:
Eq. [1] can now be written as:
Since L is lower triangular, the solution is found by two back-substitutions. The wavelet transform can be applied to the data matrix A and coefficient thresholding applied. Using this approach, only the significant elements of Λ need be computed. Similarly, we know from the previous section that only certain covariance matrix elements are significant in the wavelet domain. These can be related to the data domain through Eq.2:
DWT and Block Toeplitz Symmetry
In this section, the inventive method by which the wavelet transformation can be constructed to preserve a block Toeplitz structure is disclosed. Consider the example of Λ being a 32×32 matrix of 8×8 Toeplitz subblocks. We first define an orthogonal 8×8 wavelet transform, Q0 and then form the tensor product Q:
Thus, a novel scheme for the efficient implementation of multiresolution STAP is disclosed. In addition to potential throughput reductions, the multiresolution formulation affords much broader and more flexible processing capabilities.
The present invention has been described herein with reference to a particular embodiment for a particular application. Those having ordinary skill in the art and access to the present teachings will recognize additional modifications applications and embodiments within the scope thereof.
It is therefore intended by the appended claims to cover any and all such applications, modifications and embodiments within the scope of the present invention.
Accordingly,
Claims (12)
1. A system for processing a composite signal from a sensor and providing a plurality of tap weights in response thereto, said system including:
means for providing a matrix of data values;
means for performing a wavelet transformation of said matrix and providing a set of coefficients in response thereto, wherein the wavelet transformation is constructed to preserve a block Toeplitz structure;
means for sparsing said coefficients to provide a sparsed matrix; and
means for performing an inverse wavelet transformation of said sparsed matrix to provide said tap weights.
2. An image processing system comprising:
sensing means for providing a plurality of data values corresponding to a scene including a target and clutter and providing a first composite input signal in response thereto;
estimation means for processing said composite signal and providing a plurality of tap weights in response thereto, said estimation means including:
means for providing a matrix of said data values,
means for performing a wavelet transformation of said matrix and providing a set of coefficients in response thereto,
means for sparsing said coefficients to provide a sparsed matrix, and
means for performing an inverse wavelet transformation of said sparsed matrix to provide said tap weights; and
means for applying said tap weights to said composite signal to provide a clutter reduced output signal.
3. The invention of claim 2 wherein said sensor means includes an array of sensor elements.
4. The invention of claim 2 wherein said sensor means includes a radar sensor.
5. The invention of claim 2 wherein said sensor means includes an electro-optic sensor.
6. The invention of claim 2 wherein said matrix is a covariance matrix.
7. The invention of claim 2 wherein said matrix is a data matrix.
8. The invention of claim 2 further including means for applying a set of steering vectors to said estimation means.
9. The invention of claim 2 further including the step of inverting said sparsed matrix.
10. The invention of claim 2 wherein the wavelet transformation is constructed to preserve a block Toeplitz structure.
11. A method for processing a composite signal from a sensor and providing a plurality of tap weights in response thereto, said method including the steps of:
providing a matrix of data values;
performing a wavelet transformation of said matrix and providing a set of coefficients in response thereto, wherein the wavelet transformation is constructed to preserve a block Toeplitz structure;
sparsing said coefficients to provide a sparsed matrix, and
performing an inverse wavelet transformation of said sparsed matrix to provide said tap weights.
12. An image processing method including the steps of:
providing a plurality of data values corresponding to a scene including a target and clutter and providing a first composite input signal in response thereto;
processing said composite signal and providing a plurality of tap weights in response thereto, said processing step further includes the steps of:
providing a matrix of said data values,
performing a wavelet transformation of said matrix and providing a set of coefficients in response thereto,
sparsing said coefficients to provide a sparsed matrix, and
performing an inverse wavelet transformation of said sparsed matrix to provide said tap weights; and
applying said tap weights to said composite signal to provide a clutter reduced output signal.
Priority Applications (8)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/175,222 US6292592B1 (en) | 1998-10-19 | 1998-10-19 | Efficient multi-resolution space-time adaptive processor |
PCT/US1999/024265 WO2000023818A1 (en) | 1998-10-19 | 1999-10-15 | Efficient multi-resolution space-time adaptive processor |
EP99953209A EP1121607B1 (en) | 1998-10-19 | 1999-10-15 | Efficient multi-resolution space-time adaptive processor |
DE69940672T DE69940672D1 (en) | 1998-10-19 | 1999-10-15 | EFFICIENT MULTI-RESOLUTION SPACE-TIME-ADAPTIVE PROCESSOR |
ES99953209T ES2325234T3 (en) | 1998-10-19 | 1999-10-15 | ADAPTIVE PROCESSOR SPACE-TEMPORARY MULTIRESOLUTION EFFICIENT. |
IL14271699A IL142716A0 (en) | 1998-10-19 | 1999-10-15 | Efficient multi-resolution space-time adaptive processor |
AT99953209T ATE427505T1 (en) | 1998-10-19 | 1999-10-15 | EFFICIENT MULTIPLE RESOLUTION SPACE-TIME ADAPTIVE PROCESSOR |
IL142716A IL142716A (en) | 1998-10-19 | 2001-04-19 | Efficient multi-resolution space-time adaptive processor |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/175,222 US6292592B1 (en) | 1998-10-19 | 1998-10-19 | Efficient multi-resolution space-time adaptive processor |
Publications (1)
Publication Number | Publication Date |
---|---|
US6292592B1 true US6292592B1 (en) | 2001-09-18 |
Family
ID=22639445
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US09/175,222 Expired - Lifetime US6292592B1 (en) | 1998-10-19 | 1998-10-19 | Efficient multi-resolution space-time adaptive processor |
Country Status (7)
Country | Link |
---|---|
US (1) | US6292592B1 (en) |
EP (1) | EP1121607B1 (en) |
AT (1) | ATE427505T1 (en) |
DE (1) | DE69940672D1 (en) |
ES (1) | ES2325234T3 (en) |
IL (2) | IL142716A0 (en) |
WO (1) | WO2000023818A1 (en) |
Cited By (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6393137B1 (en) * | 1999-06-17 | 2002-05-21 | Raytheon Company | Multi-resolution object classification method employing kinematic features and system therefor |
US6480139B1 (en) * | 2000-10-05 | 2002-11-12 | Lockheed Martin Corporation | Jammer counting using cholesky factor diagonal elements |
US6771793B1 (en) * | 1999-02-17 | 2004-08-03 | Fuji Photo Film Co., Ltd. | Image processing method and apparatus |
US6873287B2 (en) * | 2000-11-02 | 2005-03-29 | Telefonaktiebolaget Lm Ericsson | Signal processing arrangement |
US20060056722A1 (en) * | 2004-09-14 | 2006-03-16 | Nathan Moroney | Edge preserving method and apparatus for image processing |
US7197095B1 (en) | 2001-09-26 | 2007-03-27 | Interstate Electronics Corporation | Inverse fast fourier transform (IFFT) with overlap and add |
US7259714B1 (en) * | 2005-05-04 | 2007-08-21 | Cataldo Thomas J | Unique space time adaptive system (USS) |
US7630344B1 (en) * | 2001-03-30 | 2009-12-08 | Science Applications International Corporation | Multistage reception of code division multiple access transmissions |
US20090324117A1 (en) * | 2008-06-27 | 2009-12-31 | Microsoft Corporation | Image denoising techniques |
US20100134342A1 (en) * | 2007-04-17 | 2010-06-03 | Thales | Method for Cleaning Signals for Centralized Antijamming |
US7733465B2 (en) | 2004-05-26 | 2010-06-08 | Bae Systems Information And Electronic Systems Integration Inc. | System and method for transitioning from a missile warning system to a fine tracking system in a directional infrared countermeasures system |
US7778367B1 (en) * | 2005-09-23 | 2010-08-17 | Rockwell Collins, Inc. | Method and apparatus for receiving a geo-location signal |
CN101819269A (en) * | 2010-03-19 | 2010-09-01 | 清华大学 | Space-time adaptive processing method under non-homogeneous clutter environment |
US20130261981A1 (en) * | 2012-04-03 | 2013-10-03 | Westerngeco L.L.C. | Covariance estimation using sparse wavelet representation |
CN103364769A (en) * | 2013-07-05 | 2013-10-23 | 南京邮电大学 | Adaptive estimation method for spatially non-uniform sea clutter covariance matrix |
CN103675768A (en) * | 2013-12-17 | 2014-03-26 | 中国人民解放军国防科学技术大学 | Suppression method for satellite-borne SAR radio-frequency interference |
US20150139560A1 (en) * | 2013-11-21 | 2015-05-21 | Bae Systems Information And Electronic Systems Integration Inc. | Coded image system and method thereof |
CN107167782A (en) * | 2017-06-27 | 2017-09-15 | 西安电子科技大学 | Based on the sparse reconstructing method of radar three-dimensional isomery battle array that letter miscellaneous noise ratio is maximum |
CN107728115A (en) * | 2017-09-11 | 2018-02-23 | 电子科技大学 | Ambient interferences suppressing method based on SVM after a kind of radar target imaging |
US10101445B2 (en) * | 2014-04-29 | 2018-10-16 | Research Foundation Of The City University Of New York | Power centroid radar |
US10145961B1 (en) * | 2016-04-27 | 2018-12-04 | Rockwell Collins, Inc. | Efficient transform-based quantization processing |
US10228449B2 (en) * | 2012-03-09 | 2019-03-12 | The United States Of America As Represented By The Secretary Of The Army | Method and system for jointly separating noise from signals |
CN109490840A (en) * | 2018-11-22 | 2019-03-19 | 中国人民解放军海军航空大学 | Based on the noise reduction and reconstructing method for improving the sparse radar target HRRP from encoding model |
CN109669182A (en) * | 2018-12-05 | 2019-04-23 | 南京邮电大学 | Dynamic/quiet target joint sparse imaging method of passive double-base SAR |
US20190228233A1 (en) * | 2008-05-09 | 2019-07-25 | Intuvision Inc. | Video tracking systems and methods employing cognitive vision |
CN110196412A (en) * | 2019-04-25 | 2019-09-03 | 电子科技大学 | A kind of STAP method of joint sparse |
CN110244271A (en) * | 2019-05-17 | 2019-09-17 | 中国人民解放军战略支援部队信息工程大学 | Radar emitter sorting recognition methods and device based on multiple simultaneous compressed transform |
CN111624571A (en) * | 2020-05-19 | 2020-09-04 | 哈尔滨工业大学 | Non-uniform Weibull background statistical distribution parameter estimation method based on self-adaptive tight frame |
JPWO2021161551A1 (en) * | 2020-02-12 | 2021-08-19 | ||
CN113406628A (en) * | 2021-05-12 | 2021-09-17 | 中国科学院国家空间科学中心 | Data resampling method for interference imaging altimeter |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102012500A (en) * | 2010-10-18 | 2011-04-13 | 南京鹏力科技有限公司 | Wave-measuring radar universal video signal acquisition processing circuit |
RU2510863C1 (en) * | 2012-09-04 | 2014-04-10 | Открытое акционерное общество "Федеральный научно-производственный центр "Нижегородский научно-исследовательский институт радиотехники" | Inter-scanning device for mapping passive jamming using chirp signals |
CN103605121B (en) * | 2013-11-18 | 2016-08-31 | 南京理工大学 | Wideband radar data fusion method based on rapid sparse Bayesian learning algorithm |
CN104459658B (en) * | 2014-06-19 | 2017-01-25 | 西安电子科技大学 | Unite-domain positioning space-time two-dimensional self-adaptive processing method based on data fitting |
RU2565237C1 (en) * | 2014-08-07 | 2015-10-20 | Федеральное государственное бюджетное учреждение науки Тихоокеанский океанологический институт им. В.И. Ильичева Дальневосточного отделения Российской академии наук (ТОИ ДВО РАН) | Evaluation of complex signal carrier frequency doppler shift |
CN106844284A (en) * | 2016-11-24 | 2017-06-13 | 天津津航计算技术研究所 | For the reconfigurable system of Radar Signal Processing |
RU2661941C1 (en) * | 2017-12-14 | 2018-07-23 | Российская Федерация, от имени которой выступает Государственная корпорация по космической деятельности "РОСКОСМОС" | Method for radar imaging using the radar station with the synthesized aperture |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5602760A (en) | 1994-02-02 | 1997-02-11 | Hughes Electronics | Image-based detection and tracking system and processing method employing clutter measurements and signal-to-clutter ratios |
US5612700A (en) | 1995-05-17 | 1997-03-18 | Fastman, Inc. | System for extracting targets from radar signatures |
US5619998A (en) | 1994-09-23 | 1997-04-15 | General Electric Company | Enhanced method for reducing ultrasound speckle noise using wavelet transform |
US5748143A (en) * | 1996-12-09 | 1998-05-05 | The United States Of America As Represented By The Secretary Of The Air Force | Adaptive post-doppler sequential beam processor |
US5848193A (en) * | 1997-04-07 | 1998-12-08 | The United States Of America As Represented By The Secretary Of The Navy | Wavelet projection transform features applied to real time pattern recognition |
US5907302A (en) * | 1997-12-19 | 1999-05-25 | The United States Of America As Represented By The Secretary Of The Air Force | Adaptive elevational scan processor statement of government interest |
US5933546A (en) * | 1996-05-06 | 1999-08-03 | Nec Research Institute, Inc. | Method and apparatus for multi-resolution image searching |
US5946038A (en) * | 1996-02-27 | 1999-08-31 | U.S. Philips Corporation | Method and arrangement for coding and decoding signals |
US6014897A (en) * | 1998-09-02 | 2000-01-18 | Mo; Larry Y. L. | Method and apparatus for improving sidelobe performance of sparse array using harmonic imaging |
US6038197A (en) * | 1998-07-14 | 2000-03-14 | Western Atlas International, Inc. | Efficient inversion of near singular geophysical signals |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
IT1199170B (en) * | 1984-07-27 | 1988-12-30 | Selenia Ind Elettroniche | DIGITAL RADAR SIGNAL PROCESSOR ABLE TO PERFORM THE ADAPTIVE CANCELLATION OF THE CLUTTER BY PARAMETRIC ESTIMATOR |
-
1998
- 1998-10-19 US US09/175,222 patent/US6292592B1/en not_active Expired - Lifetime
-
1999
- 1999-10-15 AT AT99953209T patent/ATE427505T1/en active
- 1999-10-15 ES ES99953209T patent/ES2325234T3/en not_active Expired - Lifetime
- 1999-10-15 DE DE69940672T patent/DE69940672D1/en not_active Expired - Lifetime
- 1999-10-15 WO PCT/US1999/024265 patent/WO2000023818A1/en active Application Filing
- 1999-10-15 IL IL14271699A patent/IL142716A0/en active IP Right Grant
- 1999-10-15 EP EP99953209A patent/EP1121607B1/en not_active Expired - Lifetime
-
2001
- 2001-04-19 IL IL142716A patent/IL142716A/en not_active IP Right Cessation
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5602760A (en) | 1994-02-02 | 1997-02-11 | Hughes Electronics | Image-based detection and tracking system and processing method employing clutter measurements and signal-to-clutter ratios |
US5619998A (en) | 1994-09-23 | 1997-04-15 | General Electric Company | Enhanced method for reducing ultrasound speckle noise using wavelet transform |
US5612700A (en) | 1995-05-17 | 1997-03-18 | Fastman, Inc. | System for extracting targets from radar signatures |
US5946038A (en) * | 1996-02-27 | 1999-08-31 | U.S. Philips Corporation | Method and arrangement for coding and decoding signals |
US5933546A (en) * | 1996-05-06 | 1999-08-03 | Nec Research Institute, Inc. | Method and apparatus for multi-resolution image searching |
US5748143A (en) * | 1996-12-09 | 1998-05-05 | The United States Of America As Represented By The Secretary Of The Air Force | Adaptive post-doppler sequential beam processor |
US5848193A (en) * | 1997-04-07 | 1998-12-08 | The United States Of America As Represented By The Secretary Of The Navy | Wavelet projection transform features applied to real time pattern recognition |
US5907302A (en) * | 1997-12-19 | 1999-05-25 | The United States Of America As Represented By The Secretary Of The Air Force | Adaptive elevational scan processor statement of government interest |
US6038197A (en) * | 1998-07-14 | 2000-03-14 | Western Atlas International, Inc. | Efficient inversion of near singular geophysical signals |
US6014897A (en) * | 1998-09-02 | 2000-01-18 | Mo; Larry Y. L. | Method and apparatus for improving sidelobe performance of sparse array using harmonic imaging |
Non-Patent Citations (4)
Title |
---|
"A Theory for Multiresolution Signal Decomposition: The Wavelet Representation," S.G. Mallat, IEEE transactions on Pattern Analysis and Machine Intelligence, vol. 11, No.7, Jul., 1989, pp. 674-693.* |
"Image Coding Using Vector Quantization: A Review," N.M. Nasrabadi et al., IEEE Transactions on Communications, vol. 36, No. 8, Aug., 1988, pp. 957-971.* |
"Image Coding Using Wavelet Transform,"M. Antonini et al., IEEE Transactions in Image Processing, vol. 1, No.2, Apr., 1992, pp. 205-220.* |
"Target Extraction From Clutter Images Using Wavelet Packet Analysis," H.J. Kim et al., Proceedings of the 1998 IEEE Radar Conference, May, 1998, 195-200. * |
Cited By (51)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6771793B1 (en) * | 1999-02-17 | 2004-08-03 | Fuji Photo Film Co., Ltd. | Image processing method and apparatus |
US6393137B1 (en) * | 1999-06-17 | 2002-05-21 | Raytheon Company | Multi-resolution object classification method employing kinematic features and system therefor |
US6480139B1 (en) * | 2000-10-05 | 2002-11-12 | Lockheed Martin Corporation | Jammer counting using cholesky factor diagonal elements |
US6873287B2 (en) * | 2000-11-02 | 2005-03-29 | Telefonaktiebolaget Lm Ericsson | Signal processing arrangement |
US7630344B1 (en) * | 2001-03-30 | 2009-12-08 | Science Applications International Corporation | Multistage reception of code division multiple access transmissions |
US20080025446A1 (en) * | 2001-09-26 | 2008-01-31 | L-3 Communications Corporation | Efficient space-time adaptive processing (stap) filter for global positioning system (gps) receivers |
US7197095B1 (en) | 2001-09-26 | 2007-03-27 | Interstate Electronics Corporation | Inverse fast fourier transform (IFFT) with overlap and add |
US7471744B2 (en) | 2001-09-26 | 2008-12-30 | L-3 Communications Corporation | Efficient space-time adaptive processing (STAP) filter for global positioning system (GPS) receivers |
US7292663B1 (en) * | 2001-09-26 | 2007-11-06 | L-3 Communications Corporation | Efficient space-time adaptive processing (STAP) filter for global positioning system (GPS) receivers |
US20080018533A1 (en) * | 2001-09-26 | 2008-01-24 | L-3 Communications Corporation | Efficient space-time adaptive processing (stap) filter for global positioning system (gps) receivers |
US20110142285A1 (en) * | 2004-05-26 | 2011-06-16 | Bae Systems Information And Electronic Systems Integration Inc. | System and method for transitioning from a missile warning system to a fine tracking system in a directional infrared countermeasures system |
US20110170087A1 (en) * | 2004-05-26 | 2011-07-14 | Bae Systems Information And Electronic Systems Integration Inc. | System and method for transitioning from a missile warning system to a fine tracking system in a directional infrared countermeasures system |
US7733465B2 (en) | 2004-05-26 | 2010-06-08 | Bae Systems Information And Electronic Systems Integration Inc. | System and method for transitioning from a missile warning system to a fine tracking system in a directional infrared countermeasures system |
US8184272B2 (en) | 2004-05-26 | 2012-05-22 | Bae Systems Information And Electronic Systems Integration Inc. | System and method for transitioning from a missile warning system to a fine tracking system in a directional infrared countermeasures system |
US8023107B2 (en) | 2004-05-26 | 2011-09-20 | Bae Systems Information And Electronic Systems Integration Inc. | System and method for transitioning from a missile warning system to a fine tracking system in a directional infrared countermeasures system |
US20060056722A1 (en) * | 2004-09-14 | 2006-03-16 | Nathan Moroney | Edge preserving method and apparatus for image processing |
US7259714B1 (en) * | 2005-05-04 | 2007-08-21 | Cataldo Thomas J | Unique space time adaptive system (USS) |
US20070200750A1 (en) * | 2005-05-04 | 2007-08-30 | Cataldo Thomas J | Unique space time adaptive system (uss) |
US7778367B1 (en) * | 2005-09-23 | 2010-08-17 | Rockwell Collins, Inc. | Method and apparatus for receiving a geo-location signal |
US20100134342A1 (en) * | 2007-04-17 | 2010-06-03 | Thales | Method for Cleaning Signals for Centralized Antijamming |
US8299955B2 (en) * | 2007-04-17 | 2012-10-30 | Thales | Method for cleaning signals for centralized antijamming |
US20190228233A1 (en) * | 2008-05-09 | 2019-07-25 | Intuvision Inc. | Video tracking systems and methods employing cognitive vision |
US8280185B2 (en) * | 2008-06-27 | 2012-10-02 | Microsoft Corporation | Image denoising techniques |
US20090324117A1 (en) * | 2008-06-27 | 2009-12-31 | Microsoft Corporation | Image denoising techniques |
CN101819269A (en) * | 2010-03-19 | 2010-09-01 | 清华大学 | Space-time adaptive processing method under non-homogeneous clutter environment |
US10228449B2 (en) * | 2012-03-09 | 2019-03-12 | The United States Of America As Represented By The Secretary Of The Army | Method and system for jointly separating noise from signals |
US20130261981A1 (en) * | 2012-04-03 | 2013-10-03 | Westerngeco L.L.C. | Covariance estimation using sparse wavelet representation |
CN103364769A (en) * | 2013-07-05 | 2013-10-23 | 南京邮电大学 | Adaptive estimation method for spatially non-uniform sea clutter covariance matrix |
CN103364769B (en) * | 2013-07-05 | 2015-06-03 | 南京邮电大学 | Adaptive estimation method for spatially non-uniform sea clutter covariance matrix |
US9445115B2 (en) * | 2013-11-21 | 2016-09-13 | Bae Systems Information And Electronic Systems Integration Inc. | Coded image system and method thereof |
US20150139560A1 (en) * | 2013-11-21 | 2015-05-21 | Bae Systems Information And Electronic Systems Integration Inc. | Coded image system and method thereof |
CN103675768A (en) * | 2013-12-17 | 2014-03-26 | 中国人民解放军国防科学技术大学 | Suppression method for satellite-borne SAR radio-frequency interference |
US10101445B2 (en) * | 2014-04-29 | 2018-10-16 | Research Foundation Of The City University Of New York | Power centroid radar |
US10145961B1 (en) * | 2016-04-27 | 2018-12-04 | Rockwell Collins, Inc. | Efficient transform-based quantization processing |
CN107167782A (en) * | 2017-06-27 | 2017-09-15 | 西安电子科技大学 | Based on the sparse reconstructing method of radar three-dimensional isomery battle array that letter miscellaneous noise ratio is maximum |
CN107167782B (en) * | 2017-06-27 | 2020-04-10 | 西安电子科技大学 | Radar three-dimensional heterogeneous array sparse reconstruction method based on signal-to-noise-ratio maximization |
CN107728115A (en) * | 2017-09-11 | 2018-02-23 | 电子科技大学 | Ambient interferences suppressing method based on SVM after a kind of radar target imaging |
CN107728115B (en) * | 2017-09-11 | 2020-08-11 | 电子科技大学 | SVM-based background interference suppression method after radar target imaging |
CN109490840A (en) * | 2018-11-22 | 2019-03-19 | 中国人民解放军海军航空大学 | Based on the noise reduction and reconstructing method for improving the sparse radar target HRRP from encoding model |
CN109490840B (en) * | 2018-11-22 | 2022-05-31 | 中国人民解放军海军航空大学 | Radar target HRRP (high-resolution Radar response) denoising and reconstructing method based on improved sparse self-coding model |
CN109669182A (en) * | 2018-12-05 | 2019-04-23 | 南京邮电大学 | Dynamic/quiet target joint sparse imaging method of passive double-base SAR |
CN109669182B (en) * | 2018-12-05 | 2022-05-17 | 南京邮电大学 | Passive bistatic SAR moving/static target joint sparse imaging method |
CN110196412A (en) * | 2019-04-25 | 2019-09-03 | 电子科技大学 | A kind of STAP method of joint sparse |
CN110196412B (en) * | 2019-04-25 | 2023-04-25 | 电子科技大学 | STAP method combining sparsity |
CN110244271A (en) * | 2019-05-17 | 2019-09-17 | 中国人民解放军战略支援部队信息工程大学 | Radar emitter sorting recognition methods and device based on multiple simultaneous compressed transform |
CN110244271B (en) * | 2019-05-17 | 2022-06-10 | 中国人民解放军战略支援部队信息工程大学 | Radar radiation source sorting and identifying method and device based on multiple synchronous compression transformation |
WO2021161551A1 (en) * | 2020-02-12 | 2021-08-19 | 株式会社デンソー | Radar device |
JPWO2021161551A1 (en) * | 2020-02-12 | 2021-08-19 | ||
CN111624571A (en) * | 2020-05-19 | 2020-09-04 | 哈尔滨工业大学 | Non-uniform Weibull background statistical distribution parameter estimation method based on self-adaptive tight frame |
CN113406628A (en) * | 2021-05-12 | 2021-09-17 | 中国科学院国家空间科学中心 | Data resampling method for interference imaging altimeter |
CN113406628B (en) * | 2021-05-12 | 2023-03-17 | 中国科学院国家空间科学中心 | Data resampling method for interference imaging altimeter |
Also Published As
Publication number | Publication date |
---|---|
DE69940672D1 (en) | 2009-05-14 |
WO2000023818A1 (en) | 2000-04-27 |
EP1121607B1 (en) | 2009-04-01 |
IL142716A0 (en) | 2002-03-10 |
IL142716A (en) | 2006-08-20 |
EP1121607A1 (en) | 2001-08-08 |
ES2325234T3 (en) | 2009-08-28 |
ATE427505T1 (en) | 2009-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6292592B1 (en) | Efficient multi-resolution space-time adaptive processor | |
Stankwitz et al. | Nonlinear apodization for sidelobe control in SAR imagery | |
Bergin et al. | Improved clutter mitigation performance using knowledge-aided space-time adaptive processing | |
Melvin et al. | Knowledge-aided signal processing: a new paradigm for radar and other advanced sensors | |
Çetin et al. | Feature-enhanced synthetic aperture radar image formation based on nonquadratic regularization | |
Wicks et al. | Space-time adaptive processing: A knowledge-based perspective for airborne radar | |
US9075129B2 (en) | Method and system for forming images by comparing subsets of image data | |
US20150061926A1 (en) | Target detection utilizing image array comparison | |
CN107480619B (en) | Noise-reduction method and system based on EEMD and the Ground Penetrating Radar B-scan image for arranging entropy | |
Huang et al. | A novel tensor technique for simultaneous narrowband and wideband interference suppression on single-channel SAR system | |
Rangaswamy et al. | A parametric multichannel detection algorithm for correlated non-Gaussian random processes | |
Lawton | A new polar Fourier transform for computer-aided tomography and spotlight synthetic aperture radar | |
Soumekh | Signal subspace fusion of uncalibrated sensors with application in SAR and diagnostic medicine | |
Bergin et al. | STAP with knowledge-aided data pre-whitening | |
Penna et al. | (Non-) homomorphic approaches to denoise intensity SAR images with non-local means and stochastic distances | |
Mao et al. | Angular superresolution of real aperture radar with high-dimensional data: Normalized projection array model and adaptive reconstruction | |
Dong et al. | Segmentation of radar imagery using the Gaussian Markov random field model | |
Zweig | Super-resolution Fourier transforms by optimisation, and ISAR imaging | |
Khosravi et al. | A tutorial and performance analysis on ENVI tools for SAR image despeckling | |
Ashour et al. | Adaptive window bandwidth selection for direction of arrival estimation of uniform velocity moving targets based relative intersection confidence interval technique | |
Jagadesh et al. | Modeling Target Detection and Performance Analysis of Electronic Countermeasures for Phased Radar. | |
Zenghui et al. | Local degrees of freedom of airborne array radar clutter for STAP | |
Yun et al. | Distance estimation scheme exploiting IR-UWB radar with clutter suppressing algorithm in indoor environments | |
Berizzi et al. | A new fast method for the reconstruction of 2-D microwave images of rotating objects | |
Leung et al. | A spatial temporal dynamical model for multipath scattering from the sea |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: RAYTHEON COMPANY, CALIFORNIA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:BRAUNREITER, DENNIS C.;SCHMITT, HARRY A.;CHEN, HAI-WEN;REEL/FRAME:010176/0156;SIGNING DATES FROM 19990616 TO 19990628 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
FEPP | Fee payment procedure |
Free format text: PAYOR NUMBER ASSIGNED (ORIGINAL EVENT CODE: ASPN); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
FPAY | Fee payment |
Year of fee payment: 4 |
|
FPAY | Fee payment |
Year of fee payment: 8 |
|
FPAY | Fee payment |
Year of fee payment: 12 |