Abstract:
The present invention relates to a motion classification device  14  comprising sensor means  141  for monitoring movements of a targeted object, said sensor means  141  being operable to output motion signals based on said monitored movements, a processing means  147  for processing said motion signals, said processing means being operable to compare said processed motion signals with a set of predetermined values such that an event of said targeted object  12  being moved from a first mode of operation to a second mode of operation can be detected by determining the difference between said processed motion signals and said set of predetermined values.

Description:
FIELD OF THE INVENTION 
       [0001]    The present invention relates to a motion classification device. 
       BACKGROUND OF THE INVENTION 
       [0002]    Motion sensor systems exist in numerous forms and are implemented in many applications. For instance, Apple&#39;s computer has implemented a Sudden Motion Sensor (SMS) in their notebook computer system which is a motion-based hardware and data protection system. Essentially, the SMS uses a triaxial accelerometer to detect sudden acceleration (for example when the computer is dropped) and prepares the relatively fragile hard disk drive mechanism for impact. The system then disengages the hard disk drive heads from the hard disk platters, preventing data loss and drive damage from a disk head crash. The drive resumes its normal operation when the SMS detects that the computer is stable again. 
         [0003]    Another form of motion sensor system, more commonly referred to as a motion detection system, is often applied to sense the presence of an intruder within a defined area. Such a system usually includes motion detection devices that are sensitive to infra-red radiation as it is propagated by objects, such as the human or other animal bodies. These bodies cause the motion detection devices to trigger an alarm when the movement of such a body is passed from one area of detection to another area of detection. 
         [0004]    The present invention discloses a motion classification device that is suitable for implementation in many applications. For example, the device could be implemented as part of monitoring equipment for monitoring mechanical motion and diagnosing any faults in engines of vehicles. The device could also be used for monitoring stationary engines including pumps and generators used in power stations. Another possible use of a motion classification device is for monitoring behaviour of animals in the wild. 
         [0005]    One further example is to implement the motion classification device in a cargo container security system. Essentially, the motion classification device can be implemented to classify different modes of transport of a cargo container so as to identify any potential threats to or theft of such a container. 
         [0006]    Cargo containers are commonly used for shipping large quantities of goods from one country to another. Cargo container based shipping is viewed as vulnerable to theft. Theft of and tampering with such shipping containers is a significant source of lost revenue to shippers and merchants. The International Cargo Security Council (ICSC) has estimated about $60 billion worth of cargo is stolen each year. In addition, there are also indirect costs such as lost sales, expedited shipments of replacement goods, disrupted customer service and so on. 
         [0007]    The transportation and shipping industries become viable targets not only for theft of cargo, but also as potential conveyors of illicit goods, for instance explosives and dangerous chemicals, and for poisoning of food carried therein. There are also reported incidents where illegal immigrants were found inside cargo containers. Therefore, measures to improve security are needed. 
         [0008]    US2005231365 relates to a battery operated cable security seal for securing the door of a cargo container using a stranded metal cable which includes monitoring electronics for monitoring the locked and tampered states of the cable. The monitoring electronics in this case is an RFID tag transmission device which transmits a fault condition (that is when an intruder tries to break the seal) to a host administrator. 
         [0009]    WO2004083078 is directed to a cargo transport system which prevents access to or operation of the container unless the container has been transported to a prescribed location. The cargo transport system contains a lock and a sensor that respectively lock and sense the environment surrounding the contents and/or of the cargo transport container. The lock and sensor are in communication with a remote monitoring location where appropriate responses can be directed to breaches in cargo transport security or in the cargo environment. 
         [0010]    WO2004022434 discloses a detector disposed within a container, the detector being capable of detecting deviations that could be indicative of possible security threats. A communication device is also disposed within the container housing. The communication device is capable of transmitting possible threat information to a central cargo data collection location. 
         [0011]    These prior art examples provide methods of detecting an intrusion or unauthorised access on a cargo container. However, none of these prior art examples provide a method of detecting a cargo container being stolen and transported from its original location to another location. 
         [0012]    One common method of providing external security measures is the surveillance of the area in which a shipping container is stored. However, this approach has met with limited success. 
       SUMMARY OF THE INVENTION 
       [0013]    In a first aspect of the present invention, there is provided a motion classification device comprising sensor means for monitoring movements of a targeted object, said sensor means being operable to output motion signals based on said monitored movements, processing means for processing said motion signals, said processing means being operable to compare said processed motion signals with a set of predetermined values such that an event of said targeted object being moved from a first mode of operation to a second mode of operation can be detected by determining the difference between said processed motion signals and said set of predetermined values. 
         [0014]    Preferably, said first mode of operation is any one of:
       (i) when said motion classification device is first fitted onto said targeted object; and   (ii) current mode of operation of said motion classification device.       
 
         [0017]    The second mode of operation may be any other mode which excludes said first mode of operation. 
         [0018]    Preferably, said sensor means includes at least one multiple-axis motion sensor. 
         [0019]    Said multiple-axis motion sensor may comprise a plurality of axis-specific motion sensors, the axis-specific motion sensors being orthogonally positioned relative to each other. 
         [0020]    The motion sensor may be any one of:
       i. an accelerometer in an inertial measurement unit;   ii. a gyroscope in an inertial measurement unit; and   iii. an accelerometer in a vibration unit.       
 
         [0024]    Preferably, said processing means includes input terminals connected to said multiple-axis motion sensor, said processing means being operable to receive said motion signals from each motion sensor independently. 
         [0025]    Preferably, said processing means is operable to sample said motion signals received from said multiple-axis motion sensor into a block of signal samples for each axis. 
         [0026]    Preferably, said processing means is further operable to generate a spectral magnitude of vector signals for at least one frequency component based on said block of signal samples from each axis. 
         [0027]    Preferably, said processing means is further operable to generate a set of classification coefficients based on said generated spectral magnitude of vector signals. 
         [0028]    Preferably, said processing means is further operable to determine statistics of said output motion signals over a period of time. 
         [0029]    The classification coefficients may be any one of:
       i. a set of general spectral classification coefficients;   ii. a set of general spectral modulation classification coefficients; and   iii. a set of determined statistics of said output motion signals.       
 
         [0033]    Said set of predetermined values may be any one of:
       i. a set of fixed predetermined values; and   ii. a range of predetermined values.       
 
         [0036]    Said processing means may be operable to further process said classification coefficients. 
         [0037]    Said processing means may be operable to perform a pre-alert routine when said set of classification coefficients is not equivalent to said set of fixed predetermined values. 
         [0038]    Said pre-alert routine may include recording temporal and/or spatial information of said motion classification device. 
         [0039]    Said processing means may be operable to trigger an alert signal when said set of classification coefficients is not within said range of predetermined values. 
         [0040]    In a second aspect of the present invention, there is provided a motion classification device comprising sensor means for monitoring movements of a targeted object, sensor means being operable to output motion signals based on said monitored movements, a processing means for processing motion signals, and to determine statistics of said motion signals over a period of time when said motion classification device is in a first mode of operation, wherein said processing means is further operable to compare said processed motion signals with said determined statistics of said motion signals such that an event of said targeted object being moved from a first mode of operation to a second mode of operation can be detected by determining the difference between said processed motion signals and said determined statistics of said motion signals. 
         [0041]    In a third aspect of the present invention, there is provided a cargo container including a motion classification device as defined in any of the statements set out above. 
         [0042]    In a fourth aspect of the present invention, there is provided a method of retrofitting a targeted object with a motion classification device comprising the steps of installing movement monitoring means for monitoring movement of said targeted object, and configuring a processing means operable to determine change in monitored movement such that an event of said targeted object being moved from a first mode of operation to a second mode of operation can be detected. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0043]    Embodiments of the present invention will now be described with reference to the accompanying drawings, wherein: 
           [0044]      FIG. 1  illustrates an overview of a motion classification device implementation in accordance with an embodiment of the present invention; 
           [0045]      FIG. 2  illustrates a schematic diagram of a motion classification device in accordance with an embodiment of the present invention; 
           [0046]      FIG. 3  illustrates a three-axis sensor in accordance with an embodiment of the present invention; 
           [0047]      FIG. 4  illustrates the steps of processing motion signals and generating classification coefficients in accordance with an embodiment of the present invention; and 
           [0048]      FIG. 5  illustrates the steps of processing motion signals and generating classification coefficients in accordance with an alternative embodiment of the present invention. 
       
    
    
     DETAILED DESCRIPTION 
       [0049]    A motion classification device for detecting an event in which a cargo container is being stolen and transferred from one mode of transport to another mode of transport is described in the following paragraphs. 
         [0050]    An overview of the implementation  10  of a motion classification device  14  in accordance with an embodiment of the present invention is illustrated in  FIG. 1 . Referring to  FIG. 1 , a cargo container  12  having a motion classification device  14  fitted therein is shown.  FIG. 1  illustrates the situation when the cargo container  12  is being stolen and transferred from one mode of transport (a truck  16 ) to another mode of transport (a ship  18 ). 
         [0051]    It is appreciated by the man skilled in the art that the motion of the cargo container  12  positioned on the back of a truck  16  will be different from the motion of the cargo container  12  positioned on a ship  14 . The motion classification device  14  continuously monitors any changes in motion of the cargo container  12 , and triggers an alert signal if it detects such a change. The alert signal may be an audible burglar alarm on board the cargo container  12 . Alternatively, instead of triggering the burglar alarm on board a cargo container  12 , a communication unit may also be used to communicate with a remote monitoring base  20  to alert the security guard that the cargo container  12  has been moved or transferred from its present location. Furthermore, the fact that an alert signal has been triggered, together with other relevant information (such as place and time when the alert signal is triggered) can be stored into a mass storage device. This information may be useful to the owner or the criminal investigation authority when the container is recovered. 
         [0052]    The motion classification device  14  will now be described in more detail with respect to  FIGS. 2 to 5 . 
         [0053]      FIG. 2  shows schematically the components of the motion classification device  14  described above. The motion classification device  14  includes a sensor system  141 , a microcontroller  143 , including an Input/Output (I/O) interface  144 , a communication unit  145 , a working memory  148 , a processor  147 , a mass storage unit  146 , and an alert system  150 . 
         [0054]    A sensor system  141  is provided at the input of the motion classification device  14  to monitor movements of the cargo container  12 . As shown in  FIG. 3 , the sensor system  141  takes the form of a three-axis motion sensor  160 . Essentially, the three-axis motion sensor  160  comprises three independent motion sensors  162 ,  164 ,  166  positioned orthogonally relative to each other. In the described embodiment, each motion sensor includes an Inertial Measurement Unit (IMU) (not shown) for providing direct measures of acceleration of an axis and rate of rotation about the same axis. It is known to the person skilled in the art that a MEMS (Micro-Electro-Mechanical System) accelerometer will detect acceleration and a MEMS gyroscope will detect the rate of rotation. Furthermore, a vibration unit (not shown) is also implemented to provide a direct measure of acceleration along three orthogonal axes. The main difference between the accelerometers in the vibration unit and those in the IMU is that the vibration unit accelerometers respond to vibrations at a higher frequency than those in the IMU. 
         [0055]    The outputs of the sensor system  141  are connected to the signal processor  147  via the I/O interface  144  of the microcontroller  143 . By this connection, the measured motion signals can be input to the signal processor  147 . The I/O interface  144  also includes an analogue-to-digital converter (ADC) (not shown) which converts the analogue output signals from the sensor system  141  into digital input signals. By means of a general purpose bus  142 , external devices (such as the sensor system  141  and the alert system  150 ) through the I/O interface  144  are in communication with the signal processor  147 . 
         [0056]    The signal processor  147  is operable to execute machine code instructions stored in a working memory  148  and/or retrievable from a mass storage unit  146 . The signal processor  147  processes the incoming signals in accordance with the method described in the forthcoming paragraphs. For clarity, a flow diagram is also included in  FIG. 4 . 
         [0057]    Conventionally, the spectral magnitude of a signal is determined by considering the samples of the three orthogonal components of a three-axis sensor measurement and is denoted by the vector sample sequence:
       s(n)={{x(n), y(n), z(n)}: n=1, 2, . . . }       
 
         [0059]    The magnitude of each vector signal sample s(n) is then given by 
         [0000]      | s ( n )|=√{square root over ( x   2 ( n )+ y   2 ( n )+ z   2 ( n ))}{square root over ( x   2 ( n )+ y   2 ( n )+ z   2 ( n ))}{square root over ( x   2 ( n )+ y   2 ( n )+ z   2 ( n ))} 
         [0060]    If the scalar signal in any or all of the three directions is a pure sinusoid, the magnitude signal sequence |s(n)|n=1, 2, . . . and so on will be a rectified sinusoid, and its spectral content will contain harmonics of the signal. These harmonics constitute unwanted artefacts of the detected motion. If there are multiple sinusoids in the original signal, then the resultant magnitude signal will not only contain harmonics, but also any inter-modulation products and harmonics thereof. The problem of the unwanted artefacts can be partially reduced by adding sufficient bias to each of the individual components x(n), y(n) and z(n) such that the polarity of each signal component does not change. This will remove the harmonic problem, and the resulting inter-modulation products. However, if the actual motion contains multiple sinusoids, adding bias will not remove any inter-modulation products between them. Moreover, even for a single pure sinusoid with a dc bias, the magnitude signal will contain components at both the fundamental frequency and twice the fundamental frequency. 
         [0061]    In the present embodiment of the invention, as shown in step S 14  of  FIG. 4 , the input signals for each axis of the sensor are read separately. This is achieved by loop function found in step S 18 , as illustrated. As noted in step S 16 , these signals are also processed separately. In detail, the Discrete Fourier Transform (DFT) of a block of N signal samples starting at sample n=n 0 +μk is denoted by 
         [0000]    
       
         
           
             
               
                 
                   
                     X 
                     m 
                   
                    
                   
                     ( 
                     
                       
                         n 
                         0 
                       
                       + 
                       
                         k 
                          
                         
                             
                         
                          
                         μ 
                       
                     
                     ) 
                   
                 
                 = 
                 
                   
                     
                       1 
                       N 
                     
                      
                     
                       
                         ∑ 
                         
                           n 
                           = 
                           0 
                         
                         
                           N 
                           - 
                           1 
                         
                       
                        
                       
                         
                           x 
                            
                           
                             ( 
                             
                               n 
                               + 
                               
                                 n 
                                 0 
                               
                               + 
                               
                                 k 
                                  
                                 
                                     
                                 
                                  
                                 μ 
                               
                             
                             ) 
                           
                         
                          
                         
                           w 
                           
                             - 
                             mn 
                           
                         
                          
                         m 
                       
                     
                   
                   = 
                   1 
                 
               
               , 
               2 
               , 
               … 
                
               
                   
               
               , 
               
                 N 
                 / 
                 2 
               
               , 
               … 
                
               
                   
               
               , 
               
                 N 
                 / 
                 2 
               
             
              
             
                 
             
           
         
       
       
         
           
             
                 
             
              
             where 
              
             
                 
             
           
         
       
       
         
           
             
                 
             
              
             
               w 
               = 
               
                 exp 
                  
                 
                   ( 
                   
                     j 
                      
                     
                         
                     
                      
                     2 
                      
                     
                       π 
                       / 
                       N 
                     
                   
                   ) 
                 
               
             
           
         
       
     
         [0062]    Similar equations are also applied for y(n) and z(n) as noted in step S 18 . 
         [0063]    Alternatively, the signals for the three-axis sensor could also be processed in parallel. As shown in  FIG. 5 , the inner loop function found in step S 18  of  FIG. 4  can be eliminated. 
         [0064]    In step S 20 , the spectral magnitude of the vector signal s(n) at the frequency corresponding to DFT cell m is calculated 
         [0000]        S   m ( n   0   +k   μ )=√{square root over (| x   m ( n   0   +k μ)| 2   +|Y   m ( n   0   +kμ )| 2   +|Z   m ( n   0   +k μ)| 2 )}{square root over (| x   m ( n   0   +k μ)| 2   +|Y   m ( n   0   +kμ )| 2   +|Z   m ( n   0   +k μ)| 2 )}{square root over (| x   m ( n   0   +k μ)| 2   +|Y   m ( n   0   +kμ )| 2   +|Z   m ( n   0   +k μ)| 2 )} m=1, 2, . . . , N/2 
         [0065]    Alternatively, the signal x(n) can be windowed prior to forming the DFT so as to remove spectral leakage effects. If a Hamming window is used, for which h(n)=1−cos(2πn/N) n−0, 2, . . . , N−1, then: 
         [0000]    
       
         
           
             
               
                 
                   X 
                   m 
                 
                  
                 
                   ( 
                   
                     
                       n 
                       0 
                     
                     + 
                     
                       k 
                        
                       
                           
                       
                        
                       μ 
                     
                   
                   ) 
                 
               
               = 
               
                 
                   
                     1 
                     N 
                   
                    
                   
                     
                       ∑ 
                       
                         n 
                         = 
                         0 
                       
                       
                         N 
                         - 
                         1 
                       
                     
                      
                     
                       
                         h 
                          
                         
                           ( 
                           n 
                           ) 
                         
                       
                        
                       
                         x 
                          
                         
                           ( 
                           
                             n 
                             + 
                             
                               n 
                               0 
                             
                             + 
                             
                               k 
                                
                               
                                   
                               
                                
                               μ 
                             
                           
                           ) 
                         
                       
                        
                       
                         w 
                         
                           - 
                           mn 
                         
                       
                        
                       m 
                     
                   
                 
                 = 
                 1 
               
             
             , 
             2 
             , 
             … 
              
             
                 
             
             , 
             
               N 
               / 
               2 
             
           
         
       
     
         [0066]    In step S 22 , a set of such magnitude spectra for k=0, 1, . . . , K−1 is derived by plotting a “waterfall display” where each spectrum is taken at intervals of μ samples. Alternatively, it can be plotted as an intensity plot (spectrogram) with frequency m along one axis and time k along the other axis. Essentially, the “display” will simply be a set of classification coefficients to be compared with a set of predetermined values so as to identify the current mode of transport and to detect whether the current mode of transport has changed from its initial mode of transport. 
         [0067]    In the present embodiment, the spectrogram information is used to identify the mode of transport of a vehicle on which the motion sensors are mounted. 
         [0068]    For steady vehicle motion, the spectrogram will show features that tend to be either constant or periodic over time. A further stage of processing is thus provided, namely to take the spectrum of the intensity in each frequency cell. With K time samples for each spectral intensity, the “modulation spectrum” at time n 0  is given by 
         [0000]    
       
         
           
             
               
                 
                   P 
                   m 
                 
                  
                 
                   ( 
                   
                     
                       n 
                       0 
                     
                     , 
                     τ 
                   
                   ) 
                 
               
               = 
               
                 
                   
                     1 
                     K 
                   
                    
                   
                     
                       ∑ 
                       
                         k 
                         = 
                         0 
                       
                       
                         K 
                         - 
                         1 
                       
                     
                      
                     
                       
                         
                           S 
                           m 
                           2 
                         
                          
                         
                           ( 
                           
                             
                               n 
                               0 
                             
                             + 
                             
                               k 
                                
                               
                                   
                               
                                
                               μ 
                             
                           
                           ) 
                         
                       
                        
                       
                         w 
                         
                           
                             - 
                             k 
                           
                            
                           
                               
                           
                            
                           τ 
                         
                       
                        
                       
                           
                       
                        
                       τ 
                     
                   
                 
                 = 
                 0 
               
             
             , 
             1 
             , 
             … 
              
             
                 
             
             , 
             
               K 
               / 
               2 
             
           
         
       
       
         
           where 
         
       
       
         
           
             w 
             = 
             
               exp 
                
               
                 ( 
                 
                   j 
                    
                   
                       
                   
                    
                   2 
                    
                   
                     π 
                     / 
                     K 
                   
                 
                 ) 
               
             
           
         
       
     
         [0069]    The set of coefficients {P m (n 0 , τ): m=1, 2, . . . , N/2; τ=0, 1 . . . K/2} (spectral modulation coefficients) is then used in a pattern recognition process to identify the type of vehicle on which the motions sensors are mounted. 
         [0070]    In step S 26 , the signal processor compares the coefficients derived from the above equations with a set of predetermined values. The set of predetermined values may be (1) a fixed set of predetermined values, and/or (2) a range of predetermined values. If the set of classification coefficients is not equal or close to the fixed set of predetermined values, or in other words the mode of transport of the container has changed from its initial mode of transport, the signal processor  147  will then record the temporal and spatial information of the motion classification device at which the mode of transport has changed. In step S 28 , the signal processor checks whether the coefficients fall within the predetermined range of values. If the set of coefficients is not within the range of predetermined values, the signal processor  147  will then send a signal through the I/O interface  144  to the alert system  150 . Alternatively, the signal processor  147  may send a signal to the communication unit  145  which in turn sends a signal to the remote monitoring base  20 . 
         [0071]    The set of predetermined values can also be updated during normal operation of the motion classification device  14 . This can be achieved by transmitting an updated set of predetermined values from a remote source to the motion classification device  14 . The motion classification device  14  receives the updated predetermined values via the communication unit  145  which in turns passes the updated predetermined values to the mass storage unit  146  so as to update the currently stored predetermined values. 
         [0072]    In an alternative embodiment of the present invention, in addition to spectral modulation coefficients, the statistics of the three-axis sensor output waveforms is also included to aid the classification process. 
         [0073]    For a sampled signal x n  n=1, 2, . . . , N, the mean μ 1  is given by 
         [0000]    
       
         
           
             
               
                 μ 
                 1 
               
                
               
                 ( 
                 x 
                 ) 
               
             
             = 
             
               
                 x 
                 _ 
               
               = 
               
                 
                   1 
                   N 
                 
                  
                 
                   
                     ∑ 
                     
                       n 
                       = 
                       1 
                     
                     N 
                   
                    
                   
                     x 
                     n 
                   
                 
               
             
           
         
       
     
         [0074]    The variance μ 2  is 
         [0000]    
       
         
           
             
               
                 μ 
                 2 
               
                
               
                 ( 
                 x 
                 ) 
               
             
             = 
             
               
                 σ 
                 x 
                 2 
               
               = 
               
                 
                   1 
                   
                     N 
                     - 
                     1 
                   
                 
                  
                 
                   
                     ∑ 
                     
                       n 
                       = 
                       1 
                     
                     N 
                   
                    
                   
                     
                       ( 
                       
                         
                           x 
                           n 
                         
                         - 
                         
                           x 
                           _ 
                         
                       
                       ) 
                     
                     2 
                   
                 
               
             
           
         
       
     
         [0075]    The skewness μ 3 * is defined as 
         [0000]    
       
         
           
             
               
                 μ 
                 3 
                 * 
               
                
               
                 ( 
                 x 
                 ) 
               
             
             = 
             
               
                 1 
                 
                   N 
                   - 
                   2 
                 
               
                
               
                 
                   ∑ 
                   
                     n 
                     = 
                     1 
                   
                   N 
                 
                  
                 
                   
                     ( 
                     
                       
                         
                           x 
                           n 
                         
                         - 
                         
                           x 
                           _ 
                         
                       
                       
                         σ 
                         x 
                       
                     
                     ) 
                   
                   3 
                 
               
             
           
         
       
     
         [0000]    and the kurtosis μ 4 * by 
         [0000]    
       
         
           
             
               
                 μ 
                 4 
                 * 
               
                
               
                 ( 
                 x 
                 ) 
               
             
             = 
             
               
                 1 
                 
                   N 
                   - 
                   2 
                 
               
                
               
                 
                   ∑ 
                   
                     n 
                     = 
                     1 
                   
                   N 
                 
                  
                 
                   
                     ( 
                     
                       
                         
                           x 
                           n 
                         
                         - 
                         
                           x 
                           _ 
                         
                       
                       
                         σ 
                         x 
                       
                     
                     ) 
                   
                   4 
                 
               
             
           
         
       
     
         [0000]    a set of 3 statistical coefficients (σ, μ 3 * and μ 4 *) will be associated with each “spectral modulation gram” for each of the three axes, along with the overall standard deviation σ xyz  defined by 
         [0000]      σ xyz =√{square root over (σ x   2 +σ y   2 +σ z   2 )} 
         [0076]    These statistical coefficients can be collected over time and stored as a set of predetermined values in the mass storage unit, which are then used for comparison with spectral coefficients for identifying any changes in the mode of transport. 
         [0077]    It will be appreciated by the person skilled in the art that further processing of the derived coefficients (such as general spectral classification coefficients, general spectral modulation classification coefficients, and determined statistics of said output motion signals) may be performed prior to the step of comparing the derived coefficients with the set of predetermined values. Examples of such further processing may include, but is not limited to, a filtering process, determining a confidence interval between the derived coefficients and the predetermined values, and scaling the derived coefficients in the frequency domain and/or time domain. 
         [0078]    In an alternative embodiment of the present invention, there is an option to average the spectra over a number of spectral “slices” prior to computing the spectrogram and the subsequent spectral modulation coefficients so as to minimise computational load. 
         [0079]    It will be appreciated that the foregoing provides description of specific embodiments of the invention and that no limitation on the scope of protection sought herein is to be implied therefrom. The scope of protection sought is to be determined from the claims, read with reference to, but not bound by, the description and drawings.