Abstract:
A system is disclosed that provides a sonar system featuring combined spatial and temporal correlation array processing. The disclosed system merges spatial-based velocity errors with temporal-based velocity errors by using a combining algorithm that is based on the optimization criteria that the output after combining has no bias error and has a minimized mean square error. With these criteria, the combining algorithm takes the form of a weighted summation of the spatial data stream and temporal data stream inputs with the weights being a function of the relative velocity errors. In combining spatial and temporal processing in this manner, the disclosed system provides improved velocity measurements for a wider range of ship&#39;s speeds. In particular, the system achieves performance levels of a spatial correlation sonar for low ship&#39;s speeds for which a temporal correlation sonar solution might not be available. The system also achieves performance levels of a temporal correlation sonar for high ship&#39;s speeds for which a spatial correlation sonar solution might not be available.

Description:
CROSS REFERENCE TO RELATED APPLICATIONS 
       [0001]    The underlying concepts, but not necessarily the language, of the following case are incorporated by reference: 
         [0002]    (1) U.S. Patent 60/894888, filed Mar. 14, 2007. 
         [0000]    If there are any contradictions or inconsistencies in language between this application and the case incorporated by reference that might affect the interpretation of the claims in this application, the claims in this application should be interpreted to be consistent with the language in this application. 
     
    
     FIELD OF THE INVENTION 
       [0003]    The present invention relates generally to SONAR systems, and, more particularly, to correlation SONAR systems. 
       BACKGROUND OF THE INVENTION 
       [0004]    SONAR (SOund Navigation And Ranging)—or sonar—is the generic name of the technology that is used to locate objects underwater. Sonar is used in marine, geological, and biological research, undersea mapping and navigation, and various commercial and military applications. 
         [0005]    An “active” sonar system is a type of sonar system in which a “projector” emits a pulse of sound and underwater microphones called “hydrophones” receive underwater sounds to be signal processed. If the transmitted pulses encounter an underwater object (a “target”), a portion of the sound is reflected as an “echo.” Knowing the speed of sound in water and the time for the sound wave to travel to the target and back, the distance between the listening-post vessel (e.g., ship, etc.) and the target can be calculated. Active sonar systems generally use highly directional beams of sound when searching for targets, which enable them to determine direction to the target, as well as the distance. 
         [0006]    Another application of active sonar processing is for measuring the velocity of the sound-projecting vessel itself. The sonar source of the vessel directs sonic pulses towards the ocean floor, and the receivers detect echoes of those pulses. The velocity of the vessel is then calculated based upon the distance traveled by the vessel between the transmission and reception of a first pulse and a second pulse. Examples of velocity-measurement sonars are spatial correlation sonar and temporal correlation sonar, which rely on selecting a maximum “correlation” between hydrophones in the case of spatial correlation or pulses in the case of temporal correlation. 
         [0007]    A spatial correlation sonar requires an array of many hydrophones with associated electronics components; consequently, spatial correlation sonar might be burdened with major ship&#39;s array installation costs. In contrast, a basic temporal correlation sonar requires only three hydrophones with reduced support electronics; as a result, temporal correlation sonar is likely to have lower ship&#39;s installation costs. A spatial system requires substantial processing to form arrays and extract ship&#39;s velocity via time-consuming array processing, whereas processor throughput requirements are reduced for a temporal system which does not need to deal with multiple velocity vectors and time-consuming array processing. 
         [0008]    With respect to operating conditions, a spatial correlation system performs well for low ship&#39;s speeds, but its performance steadily degrades as speed increases. This is because it must reduce the time between correlated pulses (i.e., the correlation time) in order to keep the correlation solution limited to involving hydrophone pairs inside the limited-size hydrophone array while allowing for ship&#39;s dynamics. A spatial system&#39;s correlation time is a compromise value to allow for (i) a maximum ship&#39;s speed in one chosen direction (e.g., forward), (ii) the use of as large a value as possible consistent with the hydrophone array dimensions, and (iii) the use of a value that allows for variations in ship&#39;s velocity estimates in a sonar cycle. These conditions may result in a spatial correlation sonar being unable to make use of the full hydrophone array size. 
         [0009]    In contrast, a temporal correlation system cannot function for low ship&#39;s speeds because this requires the use of an excessively long and often unachievable time between correlated pulses (i.e., the time separation between correlated pulses which exceeds the transmit-to-echo roundtrip time). A temporal correlation sonar, however, provides a performance enhancement for high ship&#39;s speeds by making full use of the size of a hydrophone array in all directions and permits further high speed performance improvements by use of available hydrophones that are external to a correlation sonar hydrophone array. 
         [0010]    What is needed is an integrated sonar system that provides the performance benefits of both spatial and temporal correlation sonars, as well as a means of detecting whether one of these correlation sonar systems is performing poorly, without some of the disadvantages in the prior art. 
       SUMMARY OF THE INVENTION 
       [0011]    The present invention provides a sonar system that features combined spatial and temporal correlation array processing. In accordance with the illustrative embodiment of the present invention, the disclosed sonar system merges spatial-based velocity errors with temporal-based velocity errors. The merging is accomplished by using a combining algorithm that is based on the optimization criteria that the output after combining has no bias error and has a minimized mean square error. With these criteria, the combining algorithm takes the form of a weighted summation of the spatial data stream and temporal data stream inputs with the weights being a function of the relative velocity errors. 
         [0012]    In combining spatial and temporal processing in this manner, the system of the illustrative embodiment provides improved velocity measurements for a wider range of ship&#39;s speeds. In particular, the illustrative-embodiment system achieves performance levels of a spatial correlation sonar for low ship&#39;s speeds for which a temporal correlation sonar solution might not be available. The system also achieves performance levels of a temporal correlation sonar for high ship&#39;s speeds for which a spatial correlation sonar solution might not be available. 
         [0013]    The system of the illustrative embodiment also provides an improvement in sonar system reliability by providing ship&#39;s velocity measurements in cases where one or more hydrophone-channel failures would render either the spatial or temporal systems individually inoperative. Additionally, because of the different characteristics of the spatial and temporal correlation sonar systems, each correlation system is used to monitor the performance of the other, as part of an integrity improvement that was previously unavailable with systems in the prior art. As part of the monitoring, in some embodiments an alarm is provided to indicate that the spatial or temporal correlation velocity solution is inaccurate at a particular moment. 
         [0014]    There are a number of performance benefits to being able to improve correlation system accuracy, availability, reliability and integrity. One benefit is that a velocity estimation with greater performance results in improved estimates of a ship&#39;s position (when those estimates are velocity-based) for submerged vehicles (e.g., submarines, AUVs, UUVs, etc.). This is particularly important for submerged vehicles for which Global Positioning System (GPS) fixes are not available or otherwise kept to a minimum to maintain the ship&#39;s covertness. Another benefit of having high-performance velocity estimates is that they improve the accuracy of certain on-board missile-delivery systems that employ a technique to obtain a velocity fix, which is then provided as initialization data to the missile before launch. In other words, high-availability estimates reduce the missile Circular Error Probable (CEP). 
         [0015]    The illustrative embodiment features a hydrophone-based system with correlation sonar processing. After reading this specification, those who are skilled in the art will further appreciate that the techniques described herein can be applied to other types of systems such as RADAR (RAdio Detection And Ranging). 
         [0016]    The illustrative embodiment of the present invention comprises: receiving i) a first spatial velocity estimate obtained from a spatial correlation sonar process and ii) a first temporal velocity estimate obtained from a temporal correlation sonar process; and generating a composite velocity error that is based on the first spatial velocity estimate and the first temporal velocity estimate. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS  
         [0017]      FIG. 1  depicts a block diagram of sonar processing system  100 , in accordance with the illustrative embodiment of the present invention. 
           [0018]      FIG. 2  depicts a schematic diagram of prior-art, spatial correlation sonar transceiver system  200 . 
           [0019]      FIG. 3  depicts a prior-art method for performing spatial correlation sonar processing. 
           [0020]      FIG. 4  depicts a schematic diagram of prior-art, temporal correlation sonar transceiver system  400 . 
           [0021]      FIG. 5  depicts pulse train  500  that comprises pulse pairs used for temporal correlation processing. 
           [0022]      FIG. 6  depicts the processing associated with a single processing channel for a single directional component. 
       
    
    
     DETAILED DESCRIPTION  
       [0023]      FIG. 1  depicts a block diagram of sonar processing system  100 , in accordance with the illustrative embodiment of the present invention. System  100  comprises inertial navigation system  101 , spatial correlation sonar system  102 , temporal correlation sonar system  103 , and spatial-temporal sonar combining system  104 , interconnected as shown. 
         [0024]    Inertial navigation system  101  is a navigation aid that continuously measures velocity without the need for external references, in well-known fashion. Because system  101  computes velocity independently of any external references, it is subject to integration drift, which results in progressively larger errors over time. In order to correct for the integration drift, system  101  accepts velocity corrections from an outside source, in this case from system  104 . In turn, system  101  provides its velocity measurements to system  104 . 
         [0025]    Spatial correlation sonar system  102  provides to system  104  velocity measurements that are the result of spatial correlation processing. System  102  comprises an interface to sonar transceiver system  200 , as described below and with respect to  FIG. 2 . System  102  further comprises a processor for performing the tasks described below and with respect to  FIG. 3  for the purposes of exchanging signals with transceiver system  200  and for processing signals received from transceiver system  200  in order to determine velocity measurements. 
         [0026]    Temporal correlation sonar system  103  provides to system  104  velocity measurements that are the result of temporal correlation processing. System  103  comprises an interface to sonar transceiver system  400 , as described below and with respect to  FIG. 4 . System  103  further comprises a processor for performing the tasks associated with temporal correlation array processing, as described below and with respect to  FIG. 5 . 
         [0027]    Spatial-temporal sonar combining system  104  comprises a receive interface for receiving velocity measurements from inertial navigation system  101 , spatial correlation sonar system  102 , and temporal correlation sonar system  103 , as well as a transmit interface to provide inertial navigation system  101  with velocity measurement corrections. System  104  further comprises a processor for performing the tasks associated with processing the velocity measurements from systems  101 ,  102 , and  103  for the purpose of generating the velocity measurement corrections. The salient tasks for generating the corrections are described below and with respect to  FIG. 6 . 
         [0028]      FIG. 2  depicts a schematic diagram of prior-art, spatial correlation sonar transceiver system  200  for transmitting and receiving signals for the purpose of measuring velocity. The depicted transceiver system is considered here for pedagogical purposes; after reading this specification, those who are skilled in the art will appreciate that the embodiments described herein are applicable to other transceiver systems as well. System  200 , which may be mounted on the underside of a ship, comprises transmitter  202  and receiver array  204 . Receiver array  204  comprises receivers  1  through  16 , wherein each receiver is commonly referred to as a “hydrophone.” (The terms “receiver” and “hydrophone” are used interchangeably herein.) 
         [0029]    The hydrophones in array  204  are used to perform sonar signal processing. One type of sonar signal processing is velocity-measuring, correlation sonar. A prior-art method for performing spatial correlation sonar processing is discussed with respect to  FIG. 3 . The processing that takes place as part of method  300  in  FIG. 3  is discussed, in order to provide a basis of understanding for some of the concepts in the later-discussed technique of the illustrative embodiment. 
         [0030]    At task  302  of method  300 , a series of pulses are transmitted vertically towards the ocean bottom via transmitter  202 . 
         [0031]    At task  304 , echoes are detected at each hydrophone in receiver array  204 . For a ground-referenced correlation sonar, the echoes are returned from the ocean floor. For a water-referenced correlation sonar, the echoes are returned from the water volume beneath the ship. 
         [0032]    At task  306 , pulse echo data is amplified, converted from analog to digital, and then digitally filtered to yield in-phase (“I”) and quadrature (“Q”) data for each hydrophone channel. This I and Q data contains all of the amplitude and phase information contained in the echo pulses, but is base banded and thus vastly reduced in data rate from the A-to-D-converted echo signals. 
         [0033]    In accordance with task  308 , a pulse location algorithm is employed to define the leading edge of each pulse. 
         [0034]    At task  310 , a pair of prime hydrophones is identified for each of the ship&#39;s various possible velocity vectors, given the arrangement of receiver array  204 . A channel-selection pair matrix, which includes all of the possible non-redundant ship&#39;s velocity vectors, is created as the result of performing task  310 . 
         [0035]    At task  312 , I and Q data from a first echo pulse in a reference channel is correlated with I and Q data from a later echo pulse for each of the other channels, thereby forming a correlation product for each channel-selection pair (i.e., hydrophone pair). These individual correlation products together describe a measured correlation function. 
         [0036]    At task  314 , the hydrophone pair having the best correlation product (the “best-correlated” hydrophone pair) is identified. 
         [0037]    At task  316 , an M-by-M array of hydrophone-pair correlation products is formed (e.g., M equals three, etc.) in the channel-selection pair matrix, wherein the array is centered about the best-correlated hydrophone pair from task  314 . A search for the peak of the correlation function is performed, which can possibly lie somewhere between the best-correlated hydrophone pair and another hydrophone pair in the M-by-M array. For example, an interpolation algorithm can be used on the correlation products, in order to find the location of the peak in relation to the hydrophone pairs. 
         [0038]    The correlation function that is described by the correlation products and the peak is a relationship between i) the correlation between hydrophone pairs and ii) their displacement in the x and y directions, where “x” and “y” correspond to the fore/aft and athwart-ship directions, respectively. The location of the correlation peak provides “maximum-correlation distance” components in the fore/aft and athwart-ship directions. The velocity is determined at task  318  by dividing the maximum-correlation distance for each directional component by twice the time differential between the pulses. 
         [0039]    Method  300  is repeated, periodically or sporadically, in order to provide successive velocity measurements. 
         [0040]      FIG. 4  depicts a schematic diagram of prior-art, temporal correlation sonar transceiver system  400  for measuring velocity. The depicted transceiver system is considered here for pedagogical purposes; after reading this specification, those who are skilled in the art will appreciate that the embodiments described herein are applicable to other transceiver systems as well.  FIG. 4  also depicts measured velocity vector  401 , consisting of forward-aft component velocity component  402  and athwart-ships component  403 . 
         [0041]    In accordance with the illustrative embodiment, transceiver system  400  is identical to spatial correlation sonar transceiver system  200 , shown in  FIG. 2 , except for the hydrophones that are used for velocity estimation and the velocity processing algorithms. No hardware changes are required to adapt transceiver system  200  to provide the hydrophone configuration of transceiver system  400 , in that system  400  is able to use transmitter  202  and a subset of the hydrophones from receiver array  204 , in accordance with the illustrative embodiment. The front-end processing, which includes preamplifiers, band pass filtering, and I-Q sampling, is the same for both system  200  and  400 , as well as the transmit power levels and pulse patterns. In some alternative embodiments, as those who are skilled in the art will appreciate, physically distinct transmitters and hydrophones can be provided, as part of separate spatial and temporal array systems, or additional hydrophones can be used for one type of correlation processing (e.g., temporal, etc.) in relation to the other type (e.g., spatial, etc.). 
         [0042]      FIG. 5  depicts pulse train  500  that comprises pulse pairs used for temporal correlation processing. The accompanying temporal processing uses an echo received from a first pulse at hydrophone  1  of  FIG. 4 , namely echo  501 , with a number of echoes received at hydrophone  13  of  FIG. 4 , namely echoes  502 - 1  through  502 -N for forward velocity estimation, wherein N is a positive integer greater than one. Similarly, the temporal processing uses an echo received from a first pulse at hydrophone  1  (such as echo  501 ) with a number of echoes received in hydrophone  5  (such as echoes  502 - 1  through  502 -N) for athwart-ships velocity estimation. For each directional component, the time separation (such as interval  503 ) that corresponds to maximum pulse-pair correlation is related to the ship&#39;s velocity along that component. A detailed description of a temporal correlation sonar system can be found in “Bi-Directional Temporal Correlation Sonar” U.S. Pat. No. 6,804,167, which is incorporated herein by reference. 
         [0043]      FIG. 6  depicts spatial-temporal correlation sonar combining system  104 , in accordance with the illustrative embodiment of the present invention. As discussed above and with respect to  FIG. 1 , combining system  104  receives velocity measurements from inertial navigation system  101 , spatial correlation array processing system  102 , and temporal correlation array processing system  103 , and uses those measurements to continuously generate corrections for the velocity measurements determined by inertial navigation system  101 . As those who are skilled in the art will appreciate, after reading this specification, some of the tasks depicted in  FIG. 6  can be performed in parallel or in a different order from that depicted. 
         [0044]      FIG. 6  depicts the processing associated with a single processing channel for a single directional component such as north-south. The processing depicted in  FIG. 6  can be replicated for a second processing channel, in order to cover a second directional component that is orthogonal to the first, such as east-west, as those who are skilled in the art will appreciate. 
         [0045]    System  104  relies on spatial and temporal correlation processing. As part of the correlation processing, the sonar system of the illustrative embodiment transmits multiple pairs of pulses during each period referred to as a “sonar cycle,” which refers to the interval between i) the first transmitted pulse of a burst of pulses and ii) the first transmitted pulse of the next burst of pulses. Pairs of echoes are then processed by spatial sonar processing system  102 , and the same or other pairs of echoes are also processed by temporal sonar processing system  103 . The outputs of processing systems  102  and  103  are sonar pulse-pair velocity estimates for a given directional component, wherein the estimates are used by tasks  601  and  604 , respectively. 
         [0046]    The spatial correlation processing path (i.e., path  620 ) is now described. At task  601 , a velocity estimate (i.e., the “navigator velocity”) received from inertial navigation system  101  is subtracted from a pulse-pair velocity estimate received from spatial correlation sonar processing system  102 . The difference is referred to as the raw spatial velocity error for the particular directional component being evaluated (e.g., north-south). Note that the navigator velocities for these pulse pair differences are evaluated at the correlation sonar times; in other words, the navigator velocity represents the mean velocity between the two transmit pulses and two receive pulses for the pulse pair. The resultant raw spatial velocity error is provided to task  602 , as part of a data stream of raw spatial velocity errors. 
         [0047]    At task  602 , the received data stream of raw spatial velocity errors is edited to eliminate values that are unacceptable due to excess noise, ocean bottom conditions, and other effects. In accordance with the illustrative embodiment, the raw spatial velocity error is subjected to an editing function, in which if the absolute value of an error value is in excess of an editing threshold (e.g., in velocity units of feet-per-second or knots, etc.) then the error value is discarded. The edited data stream is then provided to task  603 , and data-editing results (e.g., statistics, etc.) are provided to task  608 . 
         [0048]    At task  603 , the data within the received edited data stream is then averaged over one sonar cycle, and the result is provided to task  607 , as part of a data stream of averaged spatial velocity errors. 
         [0049]    The temporal correlation processing path (i.e., path  630 ) is described here and is similar to the spatial-related processing of tasks  601  through  603 . At task  604 , a velocity estimate (i.e., the “navigator velocity”) received from inertial navigation system  101  is subtracted from a pulse-pair velocity estimate received from temporal correlation sonar processing system  103 . The difference is referred to as the raw temporal velocity error for the particular directional component being evaluated (e.g., north-south). Note that the navigator velocities for these pulse pair differences are evaluated at the correlation sonar times; in other words, the navigator velocity represents the mean velocity between the two transmit pulses and two receive pulses for the pulse pair. The resultant raw temporal velocity error is provided to task  605 , as part of a data stream of raw temporal velocity errors. 
         [0050]    At task  605 , the received data stream of raw temporal velocity errors is edited to eliminate values that are unacceptable due to excess noise, ocean bottom conditions, and other effects. In accordance with the illustrative embodiment, the raw temporal velocity error is subjected to an editing function, in which if the absolute value of an error value is in excess of an editing threshold (e.g., in velocity units of feet-per-second or knots, etc.) then the error value is discarded. The edited data stream is then provided to task  606 , and data-editing results (e.g., statistics, etc.) are provided to task  608 . 
         [0051]    At task  606 , the data within the received edited data stream is then averaged over one sonar cycle, and the result is provided to task  607 , as part of a data stream of averaged temporal velocity errors. 
         [0052]    At task  607 , the averaged spatial and temporal velocity errors are merged. The key to the merging of the temporal and spatial data is a combining filter, which is now described. The combining filter takes as inputs the two data streams of the edited and averaged temporal and spatial velocity errors, and outputs an optimized combination of the two data streams. The selected optimization criteria are that the output has no bias error and a minimized mean square error. With these requirements, the filter takes the form of a weighted sum of the two data streams (spatial and temporal) with the weights being a function of the relative velocity errors. Thus, the filter output velocity correction estimate, V e , can be represented as: 
         [0000]        V   e =[σ 2   T /(σ 2   T +σ 2   S )]  v   S +[σ 2   S/(σ   2   T +σ 2   S )]  v   T , wherein   (1) 
         [0053]    v S =the spatial velocity error from task  603 , 
         [0054]    v T =the temporal velocity error from task  606 , 
         [0055]    σ 2   S =variance of the spatial velocity error, and 
         [0056]    σ 2   T =variance of the temporal velocity error. 
         [0000]    As already mentioned, this equation is for a single velocity component, such as north-south. 
         [0057]    Rearranging equation (1) yields: 
         [0000]        V   e =[1/(1+σ 2   S/σ   2   T )]  v   S+[ 1/(1+σ 2   T /σ 2   S )]  v   T    (2) 
         [0000]    Equation (2) shows that the filter coefficients are only a function of the relative spatial and temporal errors and not of their absolute values. 
         [0058]    It has been shown (such as in Bradley, S., Deines, K., and Rowe, F. “Acoustic Correlation Current Profiler”  IEEE Journal of Oceanic Engineering  Vol. 16 Issue 4 October 1991 pp.408-414, incorporated herein by reference) that the standard deviation of the velocity estimate for a spatial correlation sonar, σ S , is given by: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       σ 
                       S 
                     
                     = 
                     
                       
                         1 
                         
                           4 
                            
                           
                             
                               π 
                                
                               
                                   
                               
                                
                               
                                 N 
                                 S 
                               
                             
                           
                         
                       
                        
                       
                         ( 
                         
                           
                             W 
                             S 
                           
                           
                             t 
                             c 
                           
                         
                         ) 
                       
                        
                       
                         ( 
                         
                           1 
                           + 
                           
                             1 
                             SNR 
                           
                         
                         ) 
                       
                     
                   
                   , 
                   
                     wherein 
                      
                     
                       : 
                     
                   
                 
               
               
                 
                   ( 
                   3 
                   ) 
                 
               
             
           
         
       
     
         [0059]    N S =the number of independent spatial samples, 
         [0060]    W S =the width of the spatial correlation function, 
         [0061]    t c =the correlation time delay, and 
         [0062]    SNR=the signal-to-noise ratio. 
         [0063]    Additionally, since t c =RS/2V, wherein V=velocity and RS=the hydrophone separation for maximum spatial correlation, equation (3) becomes 
         [0000]    
       
         
           
             
               
                 
                   
                     σ 
                     S 
                   
                   = 
                   
                     
                       V 
                       
                         2 
                          
                         
                           
                             π 
                              
                             
                                 
                             
                              
                             
                               N 
                               S 
                             
                           
                         
                       
                     
                      
                     
                       ( 
                       
                         
                           W 
                           S 
                         
                         RS 
                       
                       ) 
                     
                      
                     
                       ( 
                       
                         1 
                         + 
                         
                           1 
                           SNR 
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   ( 
                   4 
                   ) 
                 
               
             
           
         
       
     
         [0064]    It can be assumed that the equation for the standard deviation of the velocity estimate for a temporal correlation sonar, σ T , has a similar form: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       σ 
                       T 
                     
                     = 
                     
                       
                         V 
                         
                           2 
                            
                           
                             
                               π 
                                
                               
                                   
                               
                                
                               
                                 N 
                                 T 
                               
                             
                           
                         
                       
                        
                       
                         ( 
                         
                           
                             W 
                             T 
                           
                           T 
                         
                         ) 
                       
                        
                       
                         ( 
                         
                           1 
                           + 
                           
                             1 
                             SNR 
                           
                         
                         ) 
                       
                     
                   
                   , 
                   
                     wherein 
                      
                     
                       : 
                     
                   
                 
               
               
                 
                   ( 
                   5 
                   ) 
                 
               
             
           
         
       
     
         [0065]    W T =the width of temporal correlation function, 
         [0066]    T=the correlation time for maximum temporal correlation, 
         [0067]    N T =the number of independent temporal samples, and 
         [0068]    SNR=the signal-to-noise ratio. 
         [0000]    Additionally, V is assumed to be the same for spatial- and temporal-related standard deviation, since both processing paths are measuring the same velocity. SNR is the same for both processing paths since they share the same transmitting transducer and receiving hydrophones. 
         [0069]    It is known that W T =(T/RT)W S , wherein RT is equal to the hydrophone separation for maximum temporal correlation. Thus, dividing equation (5) by equation (4) yields: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       σ 
                       T 
                     
                     / 
                     
                       σ 
                       S 
                     
                   
                   = 
                   
                     
                       
                         
                           N 
                           S 
                         
                         
                           N 
                           T 
                         
                       
                     
                      
                     
                       ( 
                       
                         RS 
                         / 
                         RT 
                       
                       ) 
                     
                   
                 
               
               
                 
                   ( 
                   6 
                   ) 
                 
               
             
           
         
       
     
         [0070]    This result shows that the filter coefficients depend on two relative parameters: 1) the number of independent samples used in the velocity calculations, and 2) the correlation distances used. In practice, the number of independent samples can be taken as the number of pulse pairs that are used in the each sonar cycle. The spatial correlation distance, RS, is the equivalent hydrophone separation that yields the maximum correlation. The temporal correlation distance, RT, is the actual hydrophone separation of the two hydrophones used in the temporal processing. This ratio (temporal/spatial) varies as a function of the relative speeds in the fore-aft and athwart-ships directions. For the spatial correlation sonar hydrophone box array illustrated in  FIG. 2  for mostly fore-aft speed, the ratio is about 4/3 (1.33) while for equal fore-aft and athwart-ships speeds, the ratio increases to about 5.6/2.8 (2). The ratio RS/RT also corresponds to a direct proportion relationship in velocity measuring bias error between the spatial and temporal correlation sonars. 
         [0071]    The filter output velocity correction estimate, V e , is then provided to the integrity filter of task  608 . 
         [0072]    At task  608 , integrity tests are performed on the velocity correction data, in accordance with the illustrative embodiment. The integrity tests are performed every sonar cycle to determine if there should or should not be an output provided from the Kalman filter of task  609 . This processing requires concurrent execution at a rate of once per sonar cycle of two integrity filters, one each for the spatial and temporal data streams. Each integrity filter receives the following data from its corresponding processing path (i.e., spatial processing path  620  and temporal processing path  630 ) every sonar cycle:
       i) the number of pulse-pair velocity estimation attempts, as made by either system  102  (for spatial processing) or system  103  (for temporal processing);   ii) the number of pulse-pair velocity estimates computed (there may have been some failed velocity estimation attempts made by system  102  or  103 ); and   iii) the raw spatial or temporal velocity errors for each pulse pair.       
 
         [0076]    Each integrity filter generates a result based in part on the “filter fullness” value for the particular spatial or temporal processing path, which value is defined as follows. As described above and with respect to tasks  602  or  605 , the raw spatial velocity data or raw temporal velocity data are subjected to an editing function that results in possibly some data being discarded. The ratio of i) pulse-pair errors surviving the editing to ii) the total number of pulse-pair velocity estimation attempts over a certain number of sonar cycles (or over a given time interval) is referred to as the “filter fullness” for the spatial data (i.e., ratio R S ) or temporal data (i.e., ratio R T ) and can be a value between 0 and 1, inclusive. 
         [0077]    Task  608  receives editing data from tasks  602  and  605 , the composite velocity error data from task  607 , and the cycle-averaged data from tasks  603  and  606 , in which the mean error values of all surviving pulse-pair velocity errors are computed over the same number of cycles or time (e.g., one sonar cycle, etc.) for each processing path. If the filter fullness of both the spatial and temporal paths is below a predetermined ratio value (for example 0.50) and are therefore unacceptable, then the Kalman filter will be instructed to provide no output, an alarm will be issued, and processing will continue with the next sonar cycle. If the filter fullness of both filters is greater than the predetermined ratio value, then the integrity testing will compare the mean error value of one path to the other. If the absolute values of these two long-term averages (for spatial and temporal) differ by more than a predetermined threshold value (in velocity units), then the Kalman filter will be instructed to provide no output and an alarm will be issued. Otherwise, the Kalman filter will be provided with the composite velocity error data from task  607 , the Kalman filter&#39;s output will be enabled, and no alarm will be issued. 
         [0078]    In some embodiments, the results from the integrity tests are presented to a status display, for the purpose of notifying an equipment operator of the velocity measurement statistics. Moreover, in some embodiments, integrity tests can be done as part of the Kalman filter processing (for example, by evaluating the Kalman input data variability or Kalman output levels) instead of, or in addition to, the pre-Kalman filter integrity tests described herein. 
         [0079]    As those who are skilled in the art will appreciate different numerical methods can be used to compute the long-term averages. For example, the edited velocity corrections provided by tasks  602  and  605  can be used in the averaging, or the averaged velocity corrections provided by tasks  603  and  606  can be used in order to reduce the number of computations needed. 
         [0080]    At task  609 , Kalman filtering, as is known in the art, is performed on the data stream of velocity corrections received from task  608 . As is well-known in the art, a spatial correlation sonar or a temporal correlation sonar is most effective when its noisy but unbiased velocity estimates are used to correct low-noise, slowly-varying inertial navigator velocity errors (i.e., varying at the Schuler rate or less). To account for the noisy component attributable to the correlation sonar velocity measurement data, the velocity correction estimates are filtered at task  609  by a Kalman filter. This results in an improved estimate of the error in the navigator velocity measurements. As those who are skilled in the art will appreciate, a type of filter different than a Kalman filter can be used, in some alternative embodiments. 
         [0081]    The Kalman-filtered velocity corrections are then provided to inertial navigation system  101  for the purpose of correcting the inertial-based, velocity estimate. 
         [0082]    It is to be understood that the disclosure teaches just one example of the illustrative embodiment and that many variations of the invention can easily be devised by those skilled in the art after reading this disclosure and that the scope of the present invention is to be determined by the following claims.