EP0562004B1 - Multi-dimensional data processing and display - Google Patents
Multi-dimensional data processing and display Download PDFInfo
- Publication number
- EP0562004B1 EP0562004B1 EP92902673A EP92902673A EP0562004B1 EP 0562004 B1 EP0562004 B1 EP 0562004B1 EP 92902673 A EP92902673 A EP 92902673A EP 92902673 A EP92902673 A EP 92902673A EP 0562004 B1 EP0562004 B1 EP 0562004B1
- Authority
- EP
- European Patent Office
- Prior art keywords
- matrix
- dimensional data
- data
- dimensional
- singular
- 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar 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/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/56—Display arrangements
- G01S7/62—Cathode-ray tube displays
- G01S7/6245—Stereoscopic displays; Three-dimensional displays; Pseudo-three dimensional displays
-
- 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar 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
- 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/80—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 ultrasonic, sonic or infrasonic waves
- G01S3/8006—Multi-channel systems specially adapted for direction-finding, i.e. having a single aerial system capable of giving simultaneous indications of the directions of different signals
-
- 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/80—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 ultrasonic, sonic or infrasonic waves
- G01S3/802—Systems for determining direction or deviation from predetermined direction
-
- 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/80—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 ultrasonic, sonic or infrasonic waves
- G01S3/86—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 ultrasonic, sonic or infrasonic waves with means for eliminating undesired waves, e.g. disturbing noises
-
- 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/04—Display arrangements
- G01S7/06—Cathode-ray tube displays or other two dimensional or three-dimensional displays
- G01S7/20—Stereoscopic displays; Three-dimensional displays; Pseudo-three-dimensional displays
-
- 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/66—Sonar 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/87—Combinations of sonar 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8977—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using special techniques for image reconstruction, e.g. FFT, geometrical transformations, spatial deconvolution, time deconvolution
Definitions
- the invention relates in general to data processing systems for real-time, multi-dimensional analysis and display. More particularly, the invention is directed to a real-time processing system that processes multi-dimensional signals received from a passive acoustic surveillance system.
- Prior art monitoring systems typically receive data from a beamformer and an array of hydrophones which are dragged behind, or mounted on, a vessel.
- the data received from the beamformer is the amplitude of the acoustic energy, within a frequency band, at a specified angle (or "bearing"), and at a discrete time interval.
- LOFAR data is the amplitude of the data received at a particular bearing in a time by frequency format. These prior art systems therefore can analyze data from only one direction, or bearing, at a single time, and cannot correlate data across several bearings while maintaining the bearing data separately.
- the present invention is not limited to processing sonar signals.
- Digital images can be formed from a variety of input data signals, including seismics, radar, radio, video and film. Because digital images require a very large amount of data storage, it is difficult to process digital images on a real-time basis.
- Yet another object of the invention is to provide a multi-dimensional image processing system that compresses data for remote transmission, and then reconstitutes the received data to accurately represent the original image.
- US-A-4 736 199 a method of selectively displaying data is known. This method analyzes large quantities of data in real-time.
- data blocks are derived from acoustic waves received by a beam forming apparatus, the data blocks being in the form of Lofargrams.
- the method described comprises selecting an initial set of data blocks from the series of data blocks, wherein the initial set of data blocks comprises an initial data block and every following nth data block in the series and then displaying this initial set.
- the multi-dimensional processing and display system of the present invention is used with a sonar tracking system.
- the sonar tracking system provides input sonar data, usually the amplitude of the received signal, at sequential time intervals, to the processing system.
- the data is arranged in a three-dimensional matrix, the dimensions being bearing, frequency and time.
- each data element in the matrix represents the amplitude of the signal received at the specified bearing, frequency and time.
- the input data (or raw data) is scaled to accentuate or suppress certain bearings, times, and/or frequencies.
- the three-dimensional matrix is separated into a number of matrices of two-dimensional data which are concatenated together along the bearing dimension to form one large two-dimensional matrix.
- the invention creates and maintains a historical database which is also concatenated with the two-dimensional matrix. This database allows loud, interfering sources and noise to be diminished and other sources of interest to be enhanced.
- the data can be analyzed efficiently using Singular Value Decomposition.
- the two-dimensional concatenated matrix is decomposed to obtain a compressed form of the matrix, in the preferred embodiment, singular vectors and singular values are obtained.
- the singular vectors are partitioned into one or more groups on the basis of their corresponding singular values.
- Certain data elements in the two-dimensional matrix are enhanced or diminished by modifying the singular values within each of the groups of singular vectors to enhance certain objects that are associated with the data.
- An enhanced two-dimensional concatenated matrix is generated by multiplying the two-dimensional concatenated matrix by a diagonal matrix of the modified singular values and a matrix of the singular vectors.
- the enhanced two-dimensional matrix has enhanced or diminished data values associated with certain data objects.
- the two-dimensional enhanced matrix is partitioned into a series of two dimensional matrices which are then arranged to form an enhanced three-dimensional matrix. All or portions of the enhanced three-dimensional matrix can then be displayed.
- the preferred embodiment displays data objects as "pipecleaner-like" objects within a three-dimensional transparent cube.
- the multi-dimensional image processing system of the present invention operates in "real-time". It processes input data signals and produces an image having three or more dimensions. It also allows the suppression of certain objects and the enhancement of other objects in the data.
- the data can be compressed for remote transmission, and then reconstituted when received at the remote location to accurately represent the original data. While the invention is described in its use in one specific application, an ocean-going vessel tracking system, it will be seen that it can easily be used in other applications.
- Sonar-processing systems are frequently deployed on ocean-going vessels to detect and track other vessels and objects, both surface and sub-surface.
- Present sonar-processing systems receive data from a set of hydrophones which may, for example, be towed behind the vessel. Data can also be obtained from hydrophones mounted on the bows of ships.
- the signals from the hydrophones are processed by a beamformer to generate a set of channels which are sensitive to acoustic energy coming from particular directions.
- the time-series of beam-formed data are passed through digital filters, and can then be displayed.
- Fig. 1 shows a monitoring system 30 which monitors acoustic energy received from seven angles, from 0° ("forward endfire") to 180° ("aft endfire”) at 30° intervals.
- the acoustic energy produced by an object 32, at a location 60° from the monitoring system 30, would be picked up most clearly by the beam 34 pointing in that direction.
- the monitoring system determines the amplitude of acoustic energy within a frequency band (or frequency bands) over a discrete time interval.
- the data is stored and analyzed, being continually updated as new data is received.
- a passing vessel 40 is shown crossing the path of the monitoring system 30 of Fig. 1.
- the monitoring system will detect the passing vessel at four discrete time intervals, labeled t0, t1, t2, and t3.
- the passing vessel is at 45° angle, or bearing, from the monitoring system.
- a typical system monitors as many as 50 angles, continuously updating output displays approximately once every ten seconds.
- the acoustic energy is received from natural sources, such as wind, waves, and marine organisms, or may come from man-made sources, such as merchant shipping, or submarines. This results in a mass of data received at any one time interval. As much as 46,000 bytes of data can be generated by the receiving devices every time interval.
- the present system processes input data in real-time, and eliminates unwanted data, such as loud sources or background noise.
- unwanted data such as loud sources or background noise.
- an operator can focus in on only a few vessels, which may be very quiet and would otherwise be obscured by other signals.
- the present system processes and displays data in a three-dimensional format, allowing data at each frequency to be processed individually and allowing all frequencies to be viewed at one time.
- Fig. 3 shows the amplitude vs. bearing representation of one acoustic source as received at time interval t 0 and at frequency f 1 .
- the present embodiment converts the graphical representation of the amplitude and bearing of the received signal to a greyscale format, or a "bearing time recording" (BTR).
- a greyscale image is formed from a set of scan lines. Each element within the scan line corresponds to a particular bearing, and the amplitude of the energy received at that bearing is displayed as a shade of grey, from white (0) to black (15), or black (0) to white (15). A collection of such scan lines results in a greyscale image, or BTR.
- the greyscale format thus represents, at one frequency band and at one time interval, the amplitude of the received signal over the range of bearings being monitored.
- the signals for the different bearings are represented as a horizontal vector, and the amplitude is represented by the intensity of the vector at different points.
- the signals shown in Fig. 3 are equivalent to the greyscale format shown in Fig. 4.
- Fig. 5 shows a series of scan intervals, shown as multiple greyscale lines, or BTR scans, received at multiple bearings over several time intervals.
- a BTR scan can be built up as time increases, showing how a signal changes over time. The use of a BTR thus allows the representation of the bearing and amplitude of the acoustic energy received at a particular frequency band to build up over time out of a stack of individual BTR scans.
- the BTR representation of bearing, amplitude and time allows the acoustic data to be represented by a matrix. Matrix analysis can then be used to describe patterns found in the BTR.
- the BTR is treated as a graphic representation of a matrix, and Singular Value Decomposition is used decompose the BTR matrix into its temporal and bearing components.
- Matrix analysis using singular values and singular vectors is well known in the prior art.
- the following publications describe such matrix analysis in detail: Digital Spectral Analysis with Applications , S.L. Marple, 1987; Matrix Computations , G.H. Golub and C.F. Van Loan, 1989; "Singular Value Decomposition and Least Squares Solutions", Numerical Math , G.H. Golub and C. Reinsch, 1970; LINPAC User's Guide , J.J. Dongarra, et.al., 1979; and "An Improved Algorithm for Computing Singular Value Decomposition", T. F. Chan, Communications of the ACM , 1982. Matrix analysis through the use of eigenvectors and eigenvalues are also well known in the prior art.
- Fig. 6 shows the output display of a BTR image 70.
- a matrix 72 provides a numerical representation of the BTR image for display, as illustrated. Integers from 0 to 15 are used to represent the amplitude of the signal, which had been shown in the BTR image by the intensity of the line.
- the BTR image is represented in a matrix X containing elements arranged in a two-dimensional format (shown in the box marked 60).
- This matrix can be decomposed, as described in the above references, into right singular vectors, left singular vectors, and singular values.
- the right singular vectors are arranged in matrix Q t (shown in the box marked 64).
- the left singular vectors are arranged in a matrix P (shown in the box marked 66).
- the singular values (shown in box 62) are arranged along the principal diagonal in a matrix D .
- Information in the matrix X which contains input data, can be represented by its right singular vectors Q t and the singular values D .
- the raw data can thus be represented in a substantially compressed form.
- the matrix X (60) is shown containing, in the present embodiment, the input BTR data.
- the singular values D (62) and the right singular vectors Q t (64) are calculated by an algorithm, such as outlined by Golub and Reinsch in
- the singular vectors and/or singular values are used by the real-time multi-dimensional processing system of the present invention as a filter to enhance and/or suppress features within the BTR data.
- the singular values D (62) are displayed in a diagonal form and represent the weights used to adjust the singular vectors.
- the left singular vectors P (66), in the preferred embodiment, are not explicitly calculated.
- the input data is post-multiplied by the right singular vectors to enhance the input data.
- the right singular vectors Q t (64) and singular values D (62) are used to represent important features within the input data, but in a substantially compressed form. In this numeric example, data is compressed by 82%. In larger, real-world applications, the degree of compression would be substantially greater and could reach up to 98% compression. This allows the data to be enhanced, or otherwise processed, and displayed, without losing any necessary data, but saving substantial amounts of computing time.
- the data processing system of the present invention substantially reduces the amount of computations that are necessary to process data in real time.
- a towed array of hydrophones in which data is updated every twenty seconds would perform 890,000 floating point instructions per second. Without the data "compression" before processing and enhancement, a similar system would require approximately 85,670,000 floating point instructions per second, making real-time image processing difficult, if not impossible.
- the data processing system uses singular value decomposition to describe patterns, remove unwanted components, and isolate and analyze components of interest in the BTR image.
- eigenvector decomposition of the cross-product matrix of the BTR data may be used to decompose the BTR data.
- the cross-product of the data matrix X is either X'X or XX' .
- Eigenvector decomposition is also well known in the prior art.
- the present invention stores multiple BTR images for multiple frequency bands in a three-dimensional matrix.
- the BTR images 80, 82, 84, 86 and 88 are stacked side-by side, creating a data "cube" 89.
- the dimensions of the cube are bearing by time by frequency.
- the amplitude of a signal can be determined by reading the value stored in the cube at any given time, bearing, and frequency.
- FIG. 9 one aspect of the present invention is shown, which is to display the data which is stored in a three-dimensional matrix as a three-dimensional "transparent cube" 90.
- the paths of passing vessels appear as "pipecleaner-like" objects within the transparent cube.
- the three-dimensional data cube 90 of Fig. 9 is shown on its three sides.
- the three-dimensional data cube, arranged as a matrix of data elements, can be "sliced” in any of the three directions, resulting in three faces.
- Face "A” (102) is a BTR image showing time vs. bearing at a specified frequency.
- Face "B” (104) is a "FRAZ" face, which shows, for a specified time, the bearing vs. frequency of the signals.
- Face "C” (100) the "LOFARgram” face shows the time vs. frequency for a specified bearing.
- the data can be accessed and visualized from any one of the three different faces of the sonar data cube 90.
- the "BTR" face is particularly useful because it displays multiple sources of acoustic energy as a function of the locations of the energy sources over time.
- Fig. 11 shows a BTR 110 for a particular frequency band.
- Natural sources of background energy 112 including wind, waves and marine organisms appear as a random, stippled structure 111 throughout the BTR display.
- Acoustic energy produced by merchant shipping, 114 and 116 appear as vertical traces 115 and 117, corresponding to their general direction and track histories.
- Other quieter sources, such as submarines 118 may appear as faint traces 119, or may be obscured by loud, interfering sources, such as merchant ships 116 and 114.
- individual acoustic sources may display acoustic energy at a number of different frequency BTR's.
- a merchant ship 116 produces acoustic energy on frequency BTR's 120, 122 and 124.
- Another merchant ship 114 produces acoustic energy on frequency BTR's 110 and 122.
- An unknown submarine source 118 produces acoustic energy on frequency BTR's 110 and 120.
- the multi-dimensional data processing system associates similar track histories across multiple frequency bands.
- frequency BTR's 110 and 120 are associated with submarine 118, and the BTR's are -correlated and displayed as vertical trace 119.
- the track history data base may be used to enhance data signals of interest, which are usually weak sources, and removes both unwanted signals, usually loud acoustic sources, and noise.
- the system of the present invention implements an algorithm based on singular value decomposition, to aid in detecting and tracking objects. Correlating across time and frequency identifies those frequency bands containing similar track histories. . This association enables an operator to detect and isolate acoustic sources, to suppress loud sources and noise, and to enhance other, quieter sources. This selective "partitioning" of sources on the basis of the magnitude of their singular values, or correlation, across frequency bands is important for detection of weak sources in the presence of loud interference and noise.
- a flow chart for the main algorithm of the present invention is shown.
- the algorithm operates on BTR scans, for example BTR slices 110, 120, 122, and 124, for a set of frequency bands over a certain period of time.
- input BTR data for one time interval is obtained (130) and added to the three-dimensional matrix (132), and the BTR data for the oldest time interval is removed from the three-dimensional matrix (134).
- the BTR data in the form of a three-dimensional matrix Y, is input into a preprocessing function (136).
- the preprocessing function 136 weights and reformats the data into a two-dimensional matrix X.
- the data is reformatted by first separating the three-dimensional matrix into a number of two-dimensional matrices, and then concatenating the two-dimensional matrices along a common dimension.
- the two-dimensional matrices are concatenated along the bearing dimension.
- the data, reformatted into a concatenated two-dimensional matrix is then sent to the detection function 138 where the BTR data is correlated across frequencies, then analyzed, and partitioned into subspaces in order to emphasize certain signals.
- the data is sent next to the enhancement and tracking function 140 where it is filtered, enhanced, and arranged back to the form of a three-dimensional enhanced data matrix.
- the data is displayed by display function 142.
- Prior history data is a special set of bearing traces that represents the position and status of acoustic sources within the surrounding area.
- the bearing traces are in the form of singular vectors which are associated with weak sources Q ws 144, and with bearing track locations of selected contacts Q c 146.
- some of the weak sources Q ws have been weighted, to form weighted weak sources Q wws 145.
- Track history data of weak sources Q ws 144 is calculated by the Detection Function 138 and is passed back to the Preprocessing Function 136 to be added to new data. The purpose of this feedback is to build the energy from weak sources to the point where the weak sources may be distinguished from background noise. Weighted weak source Q wws 145 track history data is also passed forward to the enhancement and tracking function 140 to allow signals of interest to be enhanced, and other, unwanted signals, to be diminished.
- the prior track history data of selected contacts, Q c 146 is calculated by the enhancement and tracking function 140 and is fed back to the preprocessing function 136.
- the preprocessing function 136 weights the three-dimensional input data, separates the three-dimensional matrix into a series of two-dimensional matrices, concatenates the two-dimensional matrices together along the bearing dimension, and then concatenates the two-dimensional matrix with the prior bearing history data.
- the preprocessing function initially obtains the bearing-time input data (BTR) for multiple frequencies, in the form of a three-dimensional data matrix (150).
- BTR bearing-time input data
- the preprocessing function at 152, weights the BTR data to accentuate or suppress certain bearings, times, and/or frequencies.
- Each time has a time coefficient t i 156 which is used to appropriately scale the amplitude data at each time interval.
- Each bearing also has an associated bearing coefficient B j 154, which is used to appropriately scale the amplitude data at specified bearings.
- each frequency has an associated frequency coefficient F k 158, which is used to scale the amplitude value of the data at that frequency. Increasing the magnitude of data at selected bearings, times, and frequencies increases the value of the associated data in subsequent analysis.
- the bearing, time and frequency coefficients can either be preset or can be alterable by an operator.
- the data is then reformatted at 160 from its three-dimensional format into a concatenated, two-dimensional format.
- the three-dimensional matrix is reformatted by first separating it into a series of two-dimensional matrices and then concatenating the two-dimensional matrices together along a common dimension.
- the common dimension is bearing: B f1 : B f2 : B f3 : . . . : B fn (Where ":" indicates concatenation along the bearing dimension for frequencies f 1 , f 2 , f 3 , ..., f n . )
- the data, which had been stored as a three-dimensional data matrix is separated, or sliced into a number of two-dimensional data matrices. These two-dimensional data matrices are then concatentated along the bearing dimension, resulting in one large concatenated two-dimensional matrix. This results in the weighted BTR scans for a given set of frequencies being concatenated together to produce BTR scans that are joined along the bearing dimension, thus creating a new, elongated two-dimensional matrix.
- the matrix has bearing in one dimension and a combination of time and frequency in the other dimension.
- the track histories of weak sources Q ws 144 are continuously updated by the detection function 138 and are passed back from the detection function to the preprocessing function.
- the track histories of weak sources Q ws 144 are in the form of right singular vectors, which are derived through singular value decomposition in the detection function.
- the track histories of the weak sources, Q ws 144 are concatenated with the scaled BTR data in the two-dimensional data matrix, shown at 162.
- the track history database retains track histories of weak sources Q ws (144a to 144h) and track histories of selected contacts Q c (146a to 146h).
- the BTR data is combined with track history data every 30 seconds.
- the entire data set, containing current BTR data and historical data, is then analyzed using singular value decomposition. In the preferred embodiment, this process is repeated every 30 seconds, and the historical data base is continually updated.
- the track history data base is created and continually updated using the track history received from the detection function and the enhancement and tracking function.
- the track history database contains historical data that shows the status of the surroundings as determined for several previous time intervals. In the preferred embodiment, data is stored for the previous 30 seconds; one minute, 10 minutes, 30 minutes, one hour, two hours, three hours and four hours.
- the track history database is updated by storing the most recently received track history data of weak sources Q ws000 144a and of selected contacts Q c000 146a every 30 + n seconds, where n refers to the additional time the system takes to perform the singular value decomposition calculations. Every one minute, the track history Q ws000 144a is "passed back", or stored in the database, to represent the historical status one minute previously. This data becomes Q ws001 (shown as 144b). Similarly, in the preferred embodiment, historical data is "passed back" every 10 minutes, 30 minutes, one hour, two hours, and three hours.
- the track history of selected contacts Q c is entered into a historical database in a manner similar to the track history of quiet sources Q ws .
- the present invention thus provides for a historical database representing the data at these time intervals to be created and continuously updated.
- Historic track data is efficiently maintained in compressed form.
- Each new analysis includes the compressed, historic data, which is equivalent to a full analysis of the full (uncompressed) historic data, yet at a fraction of the computational cost.
- the contact track history data Q c 146 is continuously updated by the enhancement/ tracking function 140, and is passed back to be used by the preprocessing function. This data is also weighted and concatenated with the two-dimensional weighted BTR data, shown at 164.
- the Detection Function 138 correlates the BTR track history data across frequency and time, analyzes the data in terms of its dominant singular values and associated right singular vectors, and partitions the data into subspaces on the basis of their singular value structure.
- the weighted BTR data in the form of a two-dimensional matrix is obtained from the preprocessing function at 170.
- the BTR data is analyzed at 172 by Singular Value Decomposition (SVD), or in an alternative embodiment, through Eigenvector Decomposition (EVD) of the bearing-by-bearing correlation matrix X'X .
- Singular value decomposition results in an associated matrix of singular values and singular vectors which has been derived from the BTR data.
- EVD (X'X) ---> D, Q t (Using Eigenvector Decomposition)
- the matrix X (60) is defined in terms of its singular values D (62), its right singular vectors Q t (64), and its left singular vectors P (66).
- the singular value shown in the first row and the first column of the singular value matrix D (180), and having a value of 34.85, indicates the magnitude of the corresponding right singular vector shown in the first row of the right singular vector matrix Q t .
- the right singular vector in the first row and centered at the fifth bearing location has its largest value of .85 (shown at 182) centered at the fifth bearing location indicating the location of a contact.
- the associated left singular vector P 1 shown in the first column of the matrix P is shown increasing in magnitude (i.e.
- singular value D 2 (190), having a value of 3.99, indicates the magnitude of another, weaker source.
- the right singular vector Q t 2 the second row in matrix Q t , has its largest value, .96, located at the second bearing location indicating the location of this contact.
- the left singular vector P 2 the second column of matrix P, shows the signal growing weaker with time, indicating that this contact is moving away from the array.
- the Detection Function next, at 174, classifies the singular values into subspaces.
- the track history data corresponding to weak sources Q ws is computed by subspace selection.
- the singular values that have been derived from the BTR data are classified into subspaces based on the magnitudes of the singular values. In the preferred embodiment, singular values are classified as one of three subspace categories: loud sources, weak sources, and noise.
- Fig. 17A shows a typical plot of the logarithm of the singular values ordered according their magnitudes.
- the graph shows a log-linear relationship 181.
- the linear portion of the graph is characteristic of the singular value structure of noise.
- the left hand side of the graph shows larger values, corresponding to loud sources, and the right hand side corresponds to smaller values, typically corresponding to noise.
- the singular values are classified into one of three categories, loud sources 180, weak sources 184 and noise 182.
- the entire range of singular values need not be calculated. Instead, to economize on computer calculations, the singular values and corresponding singular vectors may be derived sequentially. Once the logarithm of the successive singular values appears to decrease in a linear manner (as shown at 181 in Fig. 17A and 17B), the singular values do not need to be calculated any longer. In alternate embodiments, calculation of the singular values can be stopped when the singular values reach a certain value, when the derivative of the singular values reach a constant level, when a certain number of singular values have been calculated, or using some other criteria. For example, in one embodiment only seven singular values are calculated for every matrix.
- the detection function modifies selected singular values based on the subspace selection performed at 174.
- the BTR data is filtered and enhanced using the modified singular values and right singular vectors.
- the detection function modifies the weights of the singular values to reduce the importance of the singular vectors associated with loud sources 180 and noise 182 and also to enhance weak signals 184. Enhanced singular vectors whose singular values are below a certain threshold 190 may then be eliminated. In this manner, a subset of singular vectors Q ws is selected based on the weighted singular values of weak sources 184.
- Figs. 18A and 18B The method of enhancing the data is shown in Figs. 18A and 18B.
- the BTR data that was used in the example in Figs. 7A and 7B is shown.
- the singular values D enh (200) of the data have been modified to enhance certain weak signals and to diminish stronger signals.
- the singular value associated with a loud acoustic source, D 1 (180) has been reduced from 34.85 to 1.00.
- the singular value associated with a weak acoustic source, D 2 (190) has been increased from 3.99 to 10.00.
- Noise signals have similarly been diminished, in this example, D 3 (192) was reduced from 2.55 to 1.00, and D 4 (194) was reduced from 1.13 to 1.00.
- the values for the singular vectors remain the same.
- the resulting BTR data Y enh is shown at 182, and the graphic display of this data is shown at 204. In this display, the weak source is now prominently displayed, while the strong source has been diminished.
- the detection function at 178 computes a set of track history data corresponding to weak sources Qws , which describe the bearing histories of acoustic sources across many frequencies of the signal subspace. Qws may be rotated to maximize its loadings towards a particular bearing.
- the right singular vectors associated with weak sources, Q ws may be weighted by adjusting the singular values, resulting in weighted weak sources Q wws (145).
- the detection function 138 passes track history data in the form of right singular vectors Q ws 144 back to the Preprocessing Function 136. The detection function 138 then passes the weighted right singular vectors Q wws 145 forward to the enhancement and tracking function 140.
- the enhancement and tracking function 140 is shown in greater detail.
- the concatenated two-dimensional matrix of BTR data Y and the weighted weak source right singular vectors Q wws (184) are passed in from the detection function.
- the enhancement function uses the weighted weak source right singular vectors Q wws , which correspond to track histories of acoustic sources, to enhance the BTR data.
- the enhanced BTR data 60 is generated at 212 by projecting the BTR data onto the weighted singular vectors.
- the adjusted right singular vectors Q wws T may be pre-multiplied by their transpose Q wws and used as a spatial filter, (Q wws Q wws T ), to enhance and/or suppress information across all BTR data.
- Contact frequency signatures S c are derived at 214 from the enhanced BTR data and are correlated with the original incoming data to determine and produce, at 216, contact track history Q c 146.
- the track histories of weak sources for certain frequency bands can be selected for BTR display.
- the enhanced track histories of weak sources may be re-scaled and/or normalized for display purposes.
- the signal subspace associated with Q ws is computed either by summing the squared left singular vectors (i.e. Y Q wws ), or by squaring and summing the enhanced data.
- a contact of interest may be identified by its contact frequency signature by examining the enhanced data along a particular bearing and time, or by examining the data through a set of scans, linked along track in the summary energy BTR, or through examination of the transparent data cube.
- the frequency signature, S c is derived through analysis of the linked data or through analysis of the transparent data cube.
- Frequency signatures may be chosen, for example, by finding the beam/time which is most characteristic of features associated with acoustic sources and read off and placed in a frequency vector.
- an acoustic contact is identified by a unique frequency signature 230.
- the frequency signature is read from the enhanced sonar data cube at a particular time and bearing, or averaged along a curved time/bearing "slice". Alternately, the frequency signature can be computed through Singular Value Decomposition of linked data in a time/bearing slice, or alternatively through Eigenvector Decomposition of its associated cross-product matrix.
- the frequency signature is used to track an acoustic source across time and bearing, by correlation with incoming data. Frequency signatures may change erratically and unpredictably when changes occur in the acoustic source itself or when the acoustic energy is modified while traveling through the constantly changing ocean environment.
- a convolution function 232 may be used to broaden the features of the frequency signature.
- the resulting signature S c (234) may also be thresholded.
- the resulting signature S c (234) is less sensitive to slight frequency variations, which could be caused by doppler shifts, frequency instabilities, or other uncertainties.
- contact track history Q c and weak source track history Q ws can be displayed.
- a summary contact track history can be computed by summing the associated squared vector Q c along the frequency dimension, or by squaring and summing the enhanced data, which can then be re-scaled and normalized for display purposes.
- the "Enhanced BTR" data is again reformatted. However this time, the reverse procedure is followed.
- the concatenated two-dimensional matrix is first partitioned along the common dimension into a series of two-dimensional matrices.
- the series of two-dimensional matrices are then arranged into a three-dimensional format.
- the enhanced BTR data is then passed as a sonar data cube for display while the contact track history Q c is sent back to the preprocessing function to be used in analyzing incoming data.
- enhanced sonar data is displayed in the form of a transparent data cube 240 on an operator interface screen.
- the operator can interactively rotate the three-dimensional transparent cube through various angles. Contacts appear as "pipecleaner-like" objects within the transparent cube 90.
- the three-dimensional cube can be rotated and displayed from different perspectives.
- the transparent cube contains enhanced and thresholded three-dimensional data, and therefore the transparent cube displays an object in a true space-time format. This type of display is not possible in prior art systems which do not eliminate noise, nor correlate enhanced data across three dimensions.
- the screen display includes time, bearing, and frequency cursors, which allow an operator to freely “travel" through the cube, displaying the data at any time, frequency and bearing. The operator can use these cursors to control the display of "slices", or planes, through the cube.
- Frequency by time 242 is shown in LOFAR format for a particular bearing.
- Frequency by bearing 244 is shown in FRAZ format for a particular time.
- Bearing by time 246, 248 is shown in BTR format for a particular frequency band or bands.
- BTR display 246 could be used to display selected contacts of interest Q c while BTR display 248 could be used to display weak sources Q ws .
- the various displays are arranged so that they share a common dimension. For example, the FRAZ display and the LOFAR display each share the frequency dimension, and the LOFAR and BTR displays share the time dimension.
- the transparent data cube can only display enhanced, thresholded data.
- Cursors in the displays are linked so that as the operator moves one cursor, cursors in the other displays will move accordingly.
- the cursor movement controls the positioning and slicing of the sonar data cube and the -displays. For example, moving the time cursor in the LOFAR display would cause the FRAZ display to change rapidly. Similarly, moving the bearing cursor in the BTR display or in the FRAZ display would cause the LOFAR display to change accordingly.
- the frequency signature 234 at the specified time and bearing is displayed in a separate window on the screen.
- Data can be scanned rapidly within the BTR, FRAZ, and LOFAR windows, giving a "cartoon" presentation on the operator screen.
- a summary of contact track history including summary energy BTR and contact summary BTR, can be displayed in multiple windows. Subbands, subintervals and bearing sectors can also be displayed.
- Subsets of the data can be displayed in a variety of ways.
- Subcubes in multiple windows can be used to display sector BTR's (i.e. selected bearings), sector summary BTR's, summary energy BTR (i.e. all contacts), and summary contact BTR (i.e. all contacts being tracked by Q c ).
- Subbands such as portions of LOFAR bands or FRAZ subbands, can also be displayed in multiple windows. Any part of the cube can be displayed as a subcube.
- Acoustic data is obtained from a two-dimensional spherical sonar array of hydrophones 260.
- the dimensions of the data that is obtained are bearing 262 and depression elevation angle (D/E) 264, with a third dimension being time 266.
- a typical spherical sonar system might have 400 data channels, covering 40 different bearings and 10 different D/E angles.
- a measure of the total (filtered) energy would be updated many times a second.
- multiple frequency bands could be monitored, in which case the data processing system of the present invention would process data in a four-dimensional format.
- the data obtained from the spherical array of hydrophones is processed similarly to data obtained by a towed array of hydrophones.
- the data is typically displayed by a bearing-by-time display, which would be presented to an operator for a selected D/E angle in a three-part display; the top portion 270 is updated at a fast rate (e.g. 6 times a second), a middle portion 272 is averaged and updated at a slower rate (e.g. once every 10 seconds), and a bottom portion 274 is averaged and updated at an even slower rate (e.g. once every 40 seconds).
- a fast rate e.g. 6 times a second
- a middle portion 272 is averaged and updated at a slower rate
- a bottom portion 274 is averaged and updated at an even slower rate (e.g. once every 40 seconds).
- each D/E by bearing "slice" would be input to the processing function that was previously described.
- Fig. 23 shows data obtained from a fixed sonar system consisting of multiple sensor arrays distributed over a large geographic area. Each of the sensor arrays provides three-dimensional data.
- the three-dimensional data cubes are arranged as a two-dimensional matrix of three-dimensional matrices.
- a new surface 280, 282, 284, 286, 288, and 290 is added to the top of each data cube.
- the new data is concatenated into one, two-dimensional matrix 292, and the matrix is processed as described earlier.
- each of the sonar data cubes 300, 304, 308, 312, 316 and 320 is processed separately into a two-dimensional matrix 302, 306, 310, 314, 318 and 322.
- the singular values and right singular vectors for each of the resulting two-dimensional matrices is then calculated.
- the singular vectors are combined into one two-dimensional matrix 324 and processed according to the invention as described earlier. In this manner, higher dimensional data may be recursively reduced and processed according to the invention described earlier.
- Acoustic data is obtained from an array of hydrophones.
- the dimensions of the data obtained are bearing (B°), time (t), and delayed time ( ⁇ t).
- the data cube thus contains input time-series data.
- the data may be used for high-resolution, multi-dimensional signal processing.
- One of the two-dimensional faces of the data cube 332 has dimensions time (t) by delayed time ( ⁇ t), and is a time-delayed replica of the data.
- This type of matrix is referred to as a Toeplitz matrix.
- the right singular values describe the auto-correlation structure of the data and may be used for filtering and high-resolution frequency estimation. (See, for example, "Eigenanalysis-Based Frequency Estimation", by S. L. Marple in Digital Spectral Analysis with Applications , Prentice-Hall, 1987.)
- a second face of the data cube 334 has dimensions bearing (B) by delayed time ( ⁇ t), and may be used for filtering and high-resolution estimation of the direction of arrival of multiple signals.
- B bearing
- ⁇ t delayed time
- the right singular vectors of this matrix describe the location of energy along the bearing dimension.
- the left-singular vectors describe the auto-correlation and cross-correlation of signals.
- this alternate embodiment of the invention combines these two high-resolution techniques into one, unified operation which may be processed in real time.
- the multi-dimensional data processing and display system of the present invention may also be used to compress data, remotely transmit the data, and then reconstitute the data at a remote location.
- the input data is decomposed using either singular value decomposition or eigenvector decomposition. Selected singular values and singular vectors are then transmitted. After the singular values and singular vectors are received, the input data is reconstituted from the singular values and singular vectors.
- the input data can be filtered or enhanced by changing the value of the singular values.
- the compressed data can be searched and accessed directly, without searching the entire database. This increases the efficiency of database search and retrieval operations.
- the real-time multi-dimensional image processing system of the present invention is not limited to passive sonar detecting systems, but can also be used with active sonar systems. It is also useful in a variety of applications, such as, for example, a medical imaging system showing an working organ.
- the signal processing system is used to show the organ from any side, and to "travel" through the organ, shcwing slices at any position, or to visualize the data as a transparent data cube.
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Acoustics & Sound (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Control Of Indicators Other Than Cathode Ray Tubes (AREA)
- Radar Systems Or Details Thereof (AREA)
- Radio Relay Systems (AREA)
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US628337 | 1990-12-14 | ||
US07/628,337 US5175710A (en) | 1990-12-14 | 1990-12-14 | Multi-dimensional data processing and display |
PCT/US1991/009342 WO1992010880A1 (en) | 1990-12-14 | 1991-12-13 | Multi-dimensional data processing and display |
Publications (3)
Publication Number | Publication Date |
---|---|
EP0562004A1 EP0562004A1 (en) | 1993-09-29 |
EP0562004A4 EP0562004A4 (ja) | 1994-08-03 |
EP0562004B1 true EP0562004B1 (en) | 1998-07-01 |
Family
ID=24518466
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP92902673A Expired - Lifetime EP0562004B1 (en) | 1990-12-14 | 1991-12-13 | Multi-dimensional data processing and display |
Country Status (10)
Country | Link |
---|---|
US (1) | US5175710A (ja) |
EP (1) | EP0562004B1 (ja) |
JP (1) | JP2813064B2 (ja) |
KR (1) | KR970003700B1 (ja) |
AT (1) | ATE167967T1 (ja) |
AU (1) | AU650826B2 (ja) |
CA (1) | CA2098284C (ja) |
DE (1) | DE69129699T2 (ja) |
IL (1) | IL100320A (ja) |
WO (1) | WO1992010880A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10161681C2 (de) * | 2001-12-14 | 2003-12-24 | Joachim Falkenhagen | Verfahren zur Erkennung von Strukturen in (n+1)-dimensionalen Räumen durch Verschiebung und n-dimensionale Aggregation |
Families Citing this family (70)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5490516A (en) * | 1990-12-14 | 1996-02-13 | Hutson; William H. | Method and system to enhance medical signals for real-time analysis and high-resolution display |
US5245587A (en) * | 1990-12-14 | 1993-09-14 | Hutson William H | Multi-dimensional signal processing and display |
US5559940A (en) * | 1990-12-14 | 1996-09-24 | Hutson; William H. | Method and system for real-time information analysis of textual material |
US5348020A (en) * | 1990-12-14 | 1994-09-20 | Hutson William H | Method and system for near real-time analysis and display of electrocardiographic signals |
US6384847B1 (en) * | 1992-03-20 | 2002-05-07 | International Business Machines Corporation | Interactive graphical method for analyzing many-dimensional data sets |
JP2836684B2 (ja) * | 1992-07-14 | 1998-12-14 | 三菱電機株式会社 | レーダ信号表示装置 |
US5465396A (en) | 1993-01-12 | 1995-11-07 | Usa Digital Radio Partners, L.P. | In-band on-channel digital broadcasting |
US5333129A (en) * | 1993-02-03 | 1994-07-26 | The Regents Of The University Of California | Acoustic imaging in the ocean using ambinet noise |
US7171016B1 (en) | 1993-11-18 | 2007-01-30 | Digimarc Corporation | Method for monitoring internet dissemination of image, video and/or audio files |
US5956624A (en) * | 1994-07-12 | 1999-09-21 | Usa Digital Radio Partners Lp | Method and system for simultaneously broadcasting and receiving digital and analog signals |
IL114471A0 (en) * | 1994-07-12 | 1996-01-31 | Usa Digital Radio Partners L P | Method and system for simultaneously broadcasting and analog signals |
US6988202B1 (en) * | 1995-05-08 | 2006-01-17 | Digimarc Corporation | Pre-filteriing to increase watermark signal-to-noise ratio |
CN1186647C (zh) * | 1996-04-12 | 2005-01-26 | 环球核心实验室有限公司 | 处理和探测地震信号的方法及装置 |
GB9614398D0 (en) * | 1996-07-09 | 2013-12-11 | Marconi Gec Ltd | Targeting system |
DE19758441B3 (de) * | 1996-07-09 | 2016-10-13 | Thomson Marconi Sonar Ltd. | Zielwertermittlungssystem |
GB9621221D0 (en) * | 1996-10-11 | 2002-08-14 | Thomson Marconi Sonar Ltd | Target discriminator |
FR2771532B1 (fr) * | 1997-10-06 | 2000-09-22 | France Etat | Procede de representation sur un ecran de visualisation d'un champ de donnees multidirectionnelles notamment de donnees fournies par une source d'ecoute passive, et systeme pour la mise en oeuvre d'un tel procede |
US6993186B1 (en) * | 1997-12-29 | 2006-01-31 | Glickman Jeff B | Energy minimization for classification, pattern recognition, sensor fusion, data compression, network reconstruction and signal processing |
FR2777088B1 (fr) * | 1998-04-01 | 2000-05-19 | France Etat | Procede de tri de signaux acoustiques emis par un mobile |
US6501705B1 (en) * | 2000-06-22 | 2002-12-31 | Lockheed Martin Corporation | Acoustic surveillance system with signature and bearing data display |
US20020091678A1 (en) * | 2001-01-05 | 2002-07-11 | Miller Nancy E. | Multi-query data visualization processes, data visualization apparatus, computer-readable media and computer data signals embodied in a transmission medium |
ITSV20010017A1 (it) * | 2001-05-28 | 2002-11-28 | Esaote Spa | Dispositivo per il rilevamento d'immagini in risonanza magnetica nucleare |
US7231061B2 (en) * | 2002-01-22 | 2007-06-12 | Digimarc Corporation | Adaptive prediction filtering for digital watermarking |
US7886151B2 (en) | 2002-01-22 | 2011-02-08 | Purdue Research Foundation | Temporal synchronization of video and audio signals |
CA2436400A1 (en) * | 2002-07-30 | 2004-01-30 | Abel G. Wolman | Geometrization for pattern recognition, data analysis, data merging, and multiple criteria decision making |
US20040185778A1 (en) * | 2003-03-20 | 2004-09-23 | Biundo Marc C. | Peripheral communication |
US7557805B2 (en) * | 2003-04-01 | 2009-07-07 | Battelle Memorial Institute | Dynamic visualization of data streams |
DE10329252A1 (de) * | 2003-06-25 | 2005-01-13 | Rolls-Royce Deutschland Ltd & Co Kg | Regelsystem für ein Flugtriebwerk |
US7292245B2 (en) * | 2004-01-20 | 2007-11-06 | Sensitron, Inc. | Method and apparatus for time series graph display |
US7288066B2 (en) * | 2004-11-01 | 2007-10-30 | Medtronic, Inc. | Data compression method for implantable medical devices |
US7545376B2 (en) * | 2004-12-06 | 2009-06-09 | George Mason Intellectual Properties, Inc. | Interactive closed-loop data entry with real-time graphical feedback |
US7710827B1 (en) * | 2006-08-01 | 2010-05-04 | University Of Alaska | Methods and systems for conducting near-field source tracking |
WO2008112445A1 (en) * | 2007-03-12 | 2008-09-18 | Raytheon Company | Systems and methods for detection and analysis of amplitude modulation of underwater sound |
US8107321B2 (en) * | 2007-06-01 | 2012-01-31 | Technische Universitat Graz And Forschungsholding Tu Graz Gmbh | Joint position-pitch estimation of acoustic sources for their tracking and separation |
US8094154B1 (en) | 2007-12-20 | 2012-01-10 | The United States Of America As Represented By The Secretary Of The Navy | Intelligent leading indicator display |
US8204637B1 (en) | 2007-12-20 | 2012-06-19 | The United States Of America As Represented By The Secretary Of The Navy | Aircraft approach to landing analysis method |
US7797130B1 (en) | 2007-12-20 | 2010-09-14 | The United States Of America As Represented By The Secretary Of The Navy | Baseline comparative leading indicator analysis |
ES2422756T3 (es) * | 2008-03-12 | 2013-09-13 | Raytheon Co | Sistema y método de sonar autónomo |
US8114156B2 (en) * | 2008-05-30 | 2012-02-14 | Edwin Burton Hatch | Flexibly compliant ceramic prosthetic meniscus for the replacement of damaged cartilage in orthopedic surgical repair or reconstruction of hip, knee, ankle, shoulder, elbow, wrist and other anatomical joints |
US8326542B2 (en) * | 2009-11-19 | 2012-12-04 | International Business Machines Corporation | Method and system for retrieving seismic data from a seismic section in bitmap format |
US8612157B2 (en) * | 2010-07-08 | 2013-12-17 | Westerngeco L.L.C. | Method to attenuate strong marine seismic noise |
JP5079934B2 (ja) * | 2011-01-18 | 2012-11-21 | パナソニック株式会社 | 車両方向特定装置、車両方向特定方法、及びそのプログラム |
DE102013000802B3 (de) * | 2013-01-18 | 2013-12-05 | Bundesrepublik Deutschland, vertreten durch das Bundesministerium der Verteidigung, dieses vertreten durch das Bundesamt für Ausrüstung, Informationstechnik und Nutzung der Bundeswehr | Sonar-Verfahren |
EP3005569A4 (en) * | 2013-05-31 | 2017-01-11 | Nokia Technologies Oy | Calibration data |
DE102014106351A1 (de) * | 2014-05-07 | 2015-11-12 | Atlas Elektronik Gmbh | Verfahren zum Orts- und/oder Zeitauflösen einer ersten Unterwasserschallquelle von einer zweiten Unterwasserschallquelle und zugehörige Vorrichtung sowie Wasserfahrzeug |
US9846228B2 (en) | 2016-04-07 | 2017-12-19 | Uhnder, Inc. | Software defined automotive radar systems |
US9689967B1 (en) | 2016-04-07 | 2017-06-27 | Uhnder, Inc. | Adaptive transmission and interference cancellation for MIMO radar |
US10261179B2 (en) | 2016-04-07 | 2019-04-16 | Uhnder, Inc. | Software defined automotive radar |
US9599702B1 (en) | 2016-04-25 | 2017-03-21 | Uhnder, Inc. | On-demand multi-scan micro doppler for vehicle |
EP3449272B1 (en) | 2016-04-25 | 2022-11-02 | Uhnder, Inc. | Vehicle radar system with a shared radar and communication system, and method for managing such a system in a vehicle |
WO2017187243A1 (en) | 2016-04-25 | 2017-11-02 | Uhnder, Inc. | Vehicular radar sensing system utilizing high rate true random number generator |
WO2017187304A2 (en) | 2016-04-25 | 2017-11-02 | Uhnder, Inc. | Digital frequency modulated continuous wave radar using handcrafted constant envelope modulation |
WO2017187299A2 (en) | 2016-04-25 | 2017-11-02 | Uhnder, Inc. | Successive signal interference mitigation |
US10573959B2 (en) | 2016-04-25 | 2020-02-25 | Uhnder, Inc. | Vehicle radar system using shaped antenna patterns |
EP3449275A4 (en) | 2016-04-25 | 2020-01-01 | Uhnder, Inc. | REDUCTION OF PMCW-PCMW INTERFERENCE |
US9791551B1 (en) | 2016-04-25 | 2017-10-17 | Uhnder, Inc. | Vehicular radar system with self-interference cancellation |
WO2017187306A1 (en) | 2016-04-25 | 2017-11-02 | Uhnder, Inc. | Adaptive filtering for fmcw interference mitigation in pmcw radar systems |
US9753121B1 (en) | 2016-06-20 | 2017-09-05 | Uhnder, Inc. | Power control for improved near-far performance of radar systems |
US9949020B1 (en) * | 2016-08-12 | 2018-04-17 | Ocean Acoustical Services and Instrumentation System | System and method for including soundscapes in online mapping utilities |
WO2018051288A1 (en) | 2016-09-16 | 2018-03-22 | Uhnder, Inc. | Virtual radar configuration for 2d array |
US11454697B2 (en) | 2017-02-10 | 2022-09-27 | Uhnder, Inc. | Increasing performance of a receive pipeline of a radar with memory optimization |
US10908272B2 (en) | 2017-02-10 | 2021-02-02 | Uhnder, Inc. | Reduced complexity FFT-based correlation for automotive radar |
US10866306B2 (en) | 2017-02-10 | 2020-12-15 | Uhnder, Inc. | Increasing performance of a receive pipeline of a radar with memory optimization |
US11105890B2 (en) | 2017-12-14 | 2021-08-31 | Uhnder, Inc. | Frequency modulated signal cancellation in variable power mode for radar applications |
US11474225B2 (en) | 2018-11-09 | 2022-10-18 | Uhnder, Inc. | Pulse digital mimo radar system |
US11681017B2 (en) | 2019-03-12 | 2023-06-20 | Uhnder, Inc. | Method and apparatus for mitigation of low frequency noise in radar systems |
CN110151175A (zh) * | 2019-04-10 | 2019-08-23 | 杭州电子科技大学 | 基于ceemd与改进小波阈值的表面肌电信号消噪方法 |
CN111260055B (zh) * | 2020-01-13 | 2023-09-01 | 腾讯科技(深圳)有限公司 | 基于三维图像识别的模型训练方法、存储介质和设备 |
US11899126B2 (en) | 2020-01-13 | 2024-02-13 | Uhnder, Inc. | Method and system for multi-chip operation of radar systems |
CN113180680B (zh) * | 2021-05-17 | 2022-06-21 | 复旦大学 | 一种改进的基于奇异谱分析的心电信号降噪方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3894348A (en) * | 1966-12-19 | 1975-07-15 | Us Navy | Relator simulator system |
US4754282A (en) * | 1970-03-25 | 1988-06-28 | The United States Of America As Represented By The Secretary Of The Navy | Improved data analysis system |
GB2165727B (en) * | 1984-10-10 | 1988-05-18 | Plessey Co Plc | A method of selectively displaying data |
US4703349A (en) * | 1985-03-19 | 1987-10-27 | Picturetel Corporation | Method and apparatus for multi-dimensional signal processing using a Short-Space Fourier transform |
US4995011A (en) * | 1989-09-20 | 1991-02-19 | Woods Hole Oceanographic Institute | Acoustic mapping system using tomographic reconstruction |
-
1990
- 1990-12-14 US US07/628,337 patent/US5175710A/en not_active Expired - Fee Related
-
1991
- 1991-12-11 IL IL10032091A patent/IL100320A/en not_active IP Right Cessation
- 1991-12-13 EP EP92902673A patent/EP0562004B1/en not_active Expired - Lifetime
- 1991-12-13 AT AT92902673T patent/ATE167967T1/de not_active IP Right Cessation
- 1991-12-13 WO PCT/US1991/009342 patent/WO1992010880A1/en active IP Right Grant
- 1991-12-13 CA CA002098284A patent/CA2098284C/en not_active Expired - Fee Related
- 1991-12-13 JP JP4502871A patent/JP2813064B2/ja not_active Expired - Lifetime
- 1991-12-13 DE DE69129699T patent/DE69129699T2/de not_active Expired - Fee Related
- 1991-12-13 AU AU91371/91A patent/AU650826B2/en not_active Ceased
- 1991-12-13 KR KR1019930701771A patent/KR970003700B1/ko not_active IP Right Cessation
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10161681C2 (de) * | 2001-12-14 | 2003-12-24 | Joachim Falkenhagen | Verfahren zur Erkennung von Strukturen in (n+1)-dimensionalen Räumen durch Verschiebung und n-dimensionale Aggregation |
Also Published As
Publication number | Publication date |
---|---|
AU650826B2 (en) | 1994-06-30 |
CA2098284A1 (en) | 1992-06-15 |
CA2098284C (en) | 1998-12-29 |
IL100320A (en) | 1996-01-19 |
JPH06503888A (ja) | 1994-04-28 |
DE69129699T2 (de) | 1999-04-08 |
WO1992010880A1 (en) | 1992-06-25 |
KR930703752A (ko) | 1993-11-30 |
IL100320A0 (en) | 1992-09-06 |
JP2813064B2 (ja) | 1998-10-22 |
KR970003700B1 (ko) | 1997-03-21 |
ATE167967T1 (de) | 1998-07-15 |
EP0562004A1 (en) | 1993-09-29 |
AU9137191A (en) | 1992-07-08 |
DE69129699D1 (de) | 1998-08-06 |
EP0562004A4 (ja) | 1994-08-03 |
US5175710A (en) | 1992-12-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP0562004B1 (en) | Multi-dimensional data processing and display | |
US5379268A (en) | Multi-dimensional signal processing and display | |
US7035166B2 (en) | 3-D forward looking sonar with fixed frame of reference for navigation | |
EP1121607B1 (en) | Efficient multi-resolution space-time adaptive processor | |
US5216640A (en) | Inverse beamforming sonar system and method | |
Cho et al. | Acoustic beam profile-based rapid underwater object detection for an imaging sonar | |
Stergiopoulos | Implementation of adaptive and synthetic-aperture processing schemes in integrated active-passive sonar systems | |
Zarnich | A fresh look at broadband passive sonar processing | |
Zhong | Three-dimensional higher-order compressive sensing for raypath separation in an acoustic waveguide | |
US6771561B2 (en) | Sonar system especially for shallow water application such as in littoral environments | |
Bell et al. | Ziv-Zakai lower bounds in bearing estimation | |
Stergiopoulos et al. | Theory and Implementation of Advanced Signal Processing for Active and Passive Sonar Systems | |
Bjornstad | Noisefield Estimation, Array Calibration and Buried Threat Detection in the | |
Cuschieri | 3-D imaging using an electronically scanned FLS | |
James et al. | Blind Estimation of Sonar Images from Diverse Speckle-Scene Models | |
LeBlanc et al. | Electronically steered and focused forward-looking scan sonar | |
Cao | Digital signal processing for a high-resolution three-dimensional sonar imaging system for autonomous underwater vehicles | |
Lawrence et al. | A comparison of spatial spectra from vertical line arrays of different aperture using maximum entropy method (MEM) beamforming | |
Rao | Optimal quadratic array detection using quadratic time-frequency and time-scale representations | |
Walker et al. | Modal decomposition for detection of plane waves with a line array receiver | |
Baggeroer et al. | Source parameter estimation in correlated spatial noise | |
Reynolds et al. | Using a multi‐dimensional matched filter algorithm to separate multipaths |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
17P | Request for examination filed |
Effective date: 19930713 |
|
AK | Designated contracting states |
Kind code of ref document: A1 Designated state(s): AT BE CH DE DK ES FR GB GR IT LI LU MC NL SE |
|
A4 | Supplementary search report drawn up and despatched | ||
AK | Designated contracting states |
Kind code of ref document: A4 Designated state(s): AT BE CH DE DK ES FR GB GR IT LI LU MC NL SE |
|
17Q | First examination report despatched |
Effective date: 19960321 |
|
GRAG | Despatch of communication of intention to grant |
Free format text: ORIGINAL CODE: EPIDOS AGRA |
|
GRAG | Despatch of communication of intention to grant |
Free format text: ORIGINAL CODE: EPIDOS AGRA |
|
GRAH | Despatch of communication of intention to grant a patent |
Free format text: ORIGINAL CODE: EPIDOS IGRA |
|
GRAH | Despatch of communication of intention to grant a patent |
Free format text: ORIGINAL CODE: EPIDOS IGRA |
|
GRAA | (expected) grant |
Free format text: ORIGINAL CODE: 0009210 |
|
AK | Designated contracting states |
Kind code of ref document: B1 Designated state(s): AT BE CH DE DK ES FR GB GR IT LI LU MC NL SE |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: NL Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 19980701 Ref country code: LI Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 19980701 Ref country code: IT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT;WARNING: LAPSES OF ITALIAN PATENTS WITH EFFECTIVE DATE BEFORE 2007 MAY HAVE OCCURRED AT ANY TIME BEFORE 2007. THE CORRECT EFFECTIVE DATE MAY BE DIFFERENT FROM THE ONE RECORDED. Effective date: 19980701 Ref country code: GR Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 19980701 Ref country code: FR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 19980701 Ref country code: ES Free format text: THE PATENT HAS BEEN ANNULLED BY A DECISION OF A NATIONAL AUTHORITY Effective date: 19980701 Ref country code: CH Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 19980701 Ref country code: BE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 19980701 Ref country code: AT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 19980701 |
|
REF | Corresponds to: |
Ref document number: 167967 Country of ref document: AT Date of ref document: 19980715 Kind code of ref document: T |
|
REG | Reference to a national code |
Ref country code: CH Ref legal event code: EP |
|
REF | Corresponds to: |
Ref document number: 69129699 Country of ref document: DE Date of ref document: 19980806 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: SE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 19981001 Ref country code: DK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 19981001 |
|
EN | Fr: translation not filed | ||
NLV1 | Nl: lapsed or annulled due to failure to fulfill the requirements of art. 29p and 29m of the patents act | ||
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: LU Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 19981213 |
|
REG | Reference to a national code |
Ref country code: CH Ref legal event code: PL |
|
PLBE | No opposition filed within time limit |
Free format text: ORIGINAL CODE: 0009261 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT |
|
26N | No opposition filed | ||
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: MC Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 19990630 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: GB Payment date: 19991220 Year of fee payment: 9 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: DE Payment date: 19991229 Year of fee payment: 9 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: GB Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20001213 |
|
GBPC | Gb: european patent ceased through non-payment of renewal fee |
Effective date: 20001213 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: DE Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20011002 |