# Patent

## Title
Audio gain control using specific-loudness-based auditory event detection

## Abstract

In one disclosed aspect, dynamic gain modifications are applied to an audio signal at least partly in response to auditory events and/or the degree of change in signal characteristics associated with said auditory event boundaries. In another aspect, an audio signal is divided into auditory events by comparing the difference in specific loudness between successive time blocks of the audio signal.

## Background
TECHNICAL FIELD
The present invention relates to audio dynamic range control methods and apparatus in which an audio processing device analyzes an audio signal and changes the level, gain or dynamic range of the audio, and all or some of the parameters of the audio gain and dynamics processing are generated as a function of auditory events. The invention also relates to computer programs for practicing such methods or controlling such apparatus.
The present invention also relates to methods and apparatus using a specific-loudness-based detection of auditory events. The invention also relates to computer programs for practicing such methods or controlling such apparatus.
BACKGROUND ART
Dynamics Processing of Audio
The techniques of automatic gain control (AGC) and dynamic range control (DRC) are well known and are a common element of many audio signal paths. In an abstract sense, both techniques measure the level of an audio signal in some manner and then gain-modify the signal by an amount that is a function of the measured level. In a linear, 1:1 dynamics processing system, the input audio is not processed and the output audio signal ideally matches the input audio signal. Additionally, if one has an audio dynamics processing system that automatically measures characteristics of the input signal and uses that measurement to control the output signal, if the input signal rises in level by 6 dB and the output signal is processed such that it only rises in level by 3 dB, then the output signal has been compressed by a ratio of 2:1 with respect to the input signal. International Publication Number WO 2006/047600 A1 (“Calculating and Adjusting the Perceived Loudness and/or the Perceived Spectral Balance of an Audio Signal” by Alan Jeffrey Seefeldt) provides a detailed overview of the five basic types of dynamics processing of audio: compression, limiting, automatic gain control (AGC), expansion and gating.
Auditory Events and Auditory Event Detection
The division of sounds into units or segments perceived as separate and distinct is sometimes referred to as “auditory event analysis” or “auditory scene analysis” (“ASA”) and the segments are sometimes referred to as “auditory events” or “audio events.” An extensive discussion of auditory scene analysis is set forth by Albert S. Bregman in his book Auditory Scene Analysts—The Perceptual Organization of Sound, Massachusetts Institute of Technology, 1991, Fourth printing, 2001, Second MIT Press paperback edition). In addition, U.S. Pat. No. 6,002,776 to Bhadkamkar, et al, Dec. 14, 1999 cites publications dating back to 1976 as “prior art work related to sound separation by auditory scene analysis.” However, the Bhadkamkar, et al patent discourages the practical use of auditory scene analysis, concluding that “[t]echniques involving auditory scene analysis, although interesting from a scientific point of view as models of human auditory processing, are currently far too computationally demanding and specialized to be considered practical techniques for sound separation until fundamental progress is made.”
A useful way to identify auditory events is set forth by Crockett and Crocket et al in various patent applications and papers listed below under the heading “Incorporation by Reference.” According to those documents, an audio signal is divided into auditory events, each of which tends to be perceived as separate and distinct, by detecting changes in spectral composition (amplitude as a function of frequency) with respect to time. This may be done, for example, by calculating the spectral content of successive time blocks of the audio signal, calculating the difference in spectral content between successive time blocks of the audio signal, and identifying an auditory event boundary as the boundary between successive time blocks when the difference in the spectral content between such successive time blocks exceeds a threshold. Alternatively, changes in amplitude with respect to time may be calculated instead of or in addition to changes in spectral composition with respect to time.
In its least computationally demanding implementation, the process divides audio into time segments by analyzing the entire frequency band (full bandwidth audio) or substantially the entire frequency band (in practical implementations, band limiting filtering at the ends of the spectrum is often employed) and giving the greatest weight to the loudest audio signal components. This approach takes advantage of a psychoacoustic phenomenon in which at smaller time scales (20 milliseconds (ms) and less) the ear may tend to focus on a single auditory event at a given time. This implies that while multiple events may be occurring at the same time, one component tends to be perceptually most prominent and may be processed individually as though it were the only event taking place. Taking advantage of this effect also allows the auditory event detection to scale with the complexity of the audio being processed. For example, if the input audio signal being processed is a solo instrument, the audio events that are identified will likely be the individual notes being played. Similarly for an input voice signal, the individual components of speech, the vowels and consonants for example, will likely be identified as individual audio elements. As the complexity of the audio increases, such as music with a drumbeat or multiple instruments and voice, the auditory event detection identifies the “most prominent” (i.e., the loudest) audio element at any given moment.
At the expense of greater computational complexity, the process may also take into consideration changes in spectral composition with respect to time in discrete frequency subbands (fixed or dynamically determined or both fixed and dynamically determined subbands) rather than the full bandwidth. This alternative approach takes into account more than one audio stream in different frequency subbands rather than assuming that only a single stream is perceptible at a particular time.
Auditory event detection may be implemented by dividing a time domain audio waveform into time intervals or blocks and then converting the data in each block to the frequency domain, using either a filter bank or a time-frequency transformation, such as the FFT. The amplitude of the spectral content of each block may be normalized in order to eliminate or reduce the effect of amplitude changes. Each resulting frequency domain representation provides an indication of the spectral content of the audio in the particular block. The spectral content of successive blocks is compared and changes greater than a threshold may be taken to indicate the temporal start or temporal end of an auditory event.
Preferably, the frequency domain data is normalized, as is described below. The degree to which the frequency domain data needs to be normalized gives an indication of amplitude. Hence, if a change in this degree exceeds a predetermined threshold that too may be taken to indicate an event boundary. Event start and end points resulting from spectral changes and from amplitude changes may be ORed together so that event boundaries resulting from either type of change are identified.
Although techniques described in said Crockett and Crockett at al applications and papers are particularly useful in connection with aspects of the present invention, other techniques for identifying auditory events and event boundaries may be employed in aspects of the present invention.
DISCLOSURE OF THE INVENTION
Conventional prior-art dynamics processing of audio involves multiplying the audio by a time-varying control signal that adjusts the gain of the audio producing a desired result. “Gain” is a scaling factor that scales the audio amplitude. This control signal may be generated on a continuous basis or from blocks of audio data, but it is generally derived by some form of measurement of the audio being processed, and its rate of change is determined by smoothing filters, sometimes with fixed characteristics and sometimes with characteristics that vary with the dynamics of the audio. For example, response times may be adjustable in accordance with changes in the magnitude or the power of the audio. Prior art methods such as automatic gain control (AGC) and dynamic range compression (DRC) do not assess in any psychoacoustically-based way the time intervals during which gain changes may be perceived as impairments and when they can be applied without imparting audible artifacts. Therefore, conventional audio dynamics processes can often introduce audible artifacts, i.e., the effects of the dynamics processing can introduce unwanted perceptible changes in the audio.
Auditory scene analysis identifies perceptually discrete auditory events, with each event occurring between two consecutive auditory event boundaries. The audible impairments caused by a gain change can be greatly reduced by ensuring that within an auditory event the gain is more nearly constant and by confining much of the change to the neighborhood of an event boundary. In the context of compressors or expanders, the response to an increase in audio level (often called the attack) may be rapid, comparable with or shorter than the minimum duration of auditory events, but the response to a decrease (the release or recovery) may be slower so that sounds that ought to appear constant or to decay gradually may be audibly disturbed. Under such circumstances, it is very beneficial to delay the gain recovery until the next boundary or to slow down the rate of change of gain during an event. For automatic gain control applications, where the medium- to long-term level or loudness of the audio is normalized and both attack and release times may therefore be long compared with the minimum duration of an auditory event, it is beneficial during events to delay changes or slow down rates of change in gain until the next event boundary for both increasing and decreasing gains.
According to one aspect of the present invention, an audio processing system receives an audio signal and analyzes and alters the gain and/or dynamic range characteristics of the audio. The dynamic range modification of the audio is often controlled by parameters of a dynamics processing system (attack and release time, compression ratio, etc.) that have significant effects on the perceptual artifacts introduced by the dynamics processing. Changes in signal characteristics with respect to time in the audio signal are detected and identified as auditory event boundaries, such that an audio segment between consecutive boundaries constitutes an auditory event in the audio signal. The characteristics of the auditory events of interest may include characteristics of the events such as perceptual strength or duration. Some of said one or more dynamics processing parameters are generated at least partly in response to auditory events and/or the degree of change in signal characteristics associated with said auditory event boundaries.
Typically, an auditory event is a segment of audio that tends to be perceived as separate and distinct. One usable measure of signal characteristics includes a measure of the spectral content of the audio, for example, as described in the cited Crockett and Crockett et al documents. All or some of the one or more audio dynamics processing parameters may be generated at least partly in response to the presence or absence and characteristics of one or more auditory events. An auditory event boundary may be identified as a change in signal characteristics with respect to time that exceeds a threshold. Alternatively, all or some of the one or more parameters may be generated at least partly in response to a continuing measure of the degree of change in signal characteristics associated with said auditory event boundaries. Although, in principle, aspects of the invention may be implemented in analog and/or digital domains, practical implementations are likely to be implemented in the digital domain in which each of the audio signals are represented by individual samples or samples within blocks of data. In this case, the signal characteristics may be the spectral content of audio within a block, the detection of changes in signal characteristics with respect to time may be the detection of changes in spectral content of audio from block to block, and auditory event temporal start and stop boundaries each coincide with a boundary of a block of data. It should be noted that for the more traditional case of performing dynamic gain changes on a sample-by-sample basis, that the auditory scene analysis described could be performed on a block basis and the resulting auditory event information being used to perform dynamic gain changes that are applied sample-by-sample.
By controlling key audio dynamics processing parameters using the results of auditory scene analysis, a dramatic reduction of audible artifacts introduced by dynamics processing may be achieved.
The present invention presents two ways of performing auditory scene analysis. The first performs spectral analysis and identifies the location of perceptible audio events that are used to control the dynamic gain parameters by identifying changes in spectral content. The second way transforms the audio into a perceptual loudness domain (that may provide more psychoacoustically relevant information than the first way) and identifies the location of auditory events that are subsequently used to control the dynamic gain parameters. It should be noted that the second way requires that the audio processing be aware of absolute acoustic reproduction levels, which may not be possible in some implementations. Presenting both methods of auditory scene analysis allows implementations of ASA-controlled dynamic gain modification using processes or devices that may or may not be calibrated to take into account absolute reproduction levels.
Aspects of the present invention are described herein in an audio dynamics processing environment that includes aspects of other inventions. Such other inventions are described in various pending United States and International patent applications of Dolby Laboratories Licensing Corporation, the owner of the present application, which applications are identified herein.



DESCRIPTION OF THE DRAWINGS
FIG. 1 is a flow chart showing an example of processing steps for performing auditory scene analysis.
FIG. 2 shows an example of block processing, windowing and performing the DFT on audio while performing the auditory scene analysis.
FIG. 3 is in the nature of a flow chart or functional block diagram, showing parallel processing in which audio is used to identify auditory events and to identify the characteristics of the auditory events such that the events and their characteristics are used to modify dynamics processing parameters.
FIG. 4 is in the nature of a flow chart or functional block diagram, showing processing in which audio is used only to identify auditory events and the event characteristics are determined from the audio event detection such that the events and their characteristics are used to modify the dynamics processing parameters.
FIG. 5 is in the nature of a flow chart or functional block diagram, showing processing in which audio is used only to identify auditory events and the event characteristics are determined from the audio event detection and such that only the characteristics of the auditory events are used to modify the dynamics processing parameters.
FIG. 6 shows a set idealized auditory filter characteristic responses that approximate critical banding on the ERB scale. The horizontal scale is frequency in Hertz and the vertical scale is level in decibels.
FIG. 7 shows the equal loudness contours of ISO 226. The horizontal scale is frequency in Hertz (logarithmic base 10 scale) and the vertical scale is sound pressure level in decibels.
FIGS. 8a-c shows idealized input/output characteristics and input gain characteristics of an audio dynamic range compressor.
FIGS. 9a-f show an example of the use of auditory events to control the release time in a digital implementation of a traditional Dynamic Range Controller (DRC) in which the gain control is derived from the Root Mean Square (RMS) power of the signal.
FIGS. 10a-f show an example of the use of auditory events to control the release time in a digital implementation of a traditional Dynamic Range Controller (DRC) in which the gain control is derived from the Root Mean Square (RMS) power of the signal for an alternate signal to that used in FIG. 9.
FIG. 11 depicts a suitable set of idealized AGC and DRC curves for the application of AGC followed by DRC in a loudness domain dynamics processing system. The goal of the combination is to make all processed audio have approximately the same perceived loudness while still maintaining at least some of the original audio's dynamics.



BEST MODE FOR CARRYING OUT THE INVENTION
Auditory Scene Analysis
Original, Non-Loudness Domain Method
In accordance with an embodiment of one aspect of the present invention, auditory scene analysis may be composed of four general processing steps as shown in a portion of FIG. 1. The first step 1-1 (“Perform Spectral Analysis”) takes a time-domain audio signal, divides it into blocks and calculates a spectral profile or spectral content for each of the blocks. Spectral analysis transforms the audio signal into the short-term frequency domain. This may be performed using any filterbank, either based on transforms or banks of bandpass filters, and in either linear or warped frequency space (such as the Bark scale or critical band, which better approximate the characteristics of the human ear). With any filterbank there exists a tradeoff between time and frequency. Greater time resolution, and hence shorter time intervals, leads to lower frequency resolution. Greater frequency resolution, and hence narrower subbands, leads to longer time intervals.
The first step, illustrated conceptually in FIG. 1 calculates the spectral content of successive time segments of the audio signal. In a practical embodiment, the ASA block size may be from any number of samples of the input audio signal, although 512 samples provide a good tradeoff of time and frequency resolution. In the second step 1-2, the differences in spectral content from block to block are determined (“Perform spectral profile difference measurements”). Thus, the second step calculates the difference in spectral content between successive time segments of the audio signal. As discussed above, a powerful indicator of the beginning or end of a perceived auditory event is believed to be a change in spectral content. In the third step 1-3 (“Identify location of auditory event boundaries”), when the spectral difference between one spectral-profile block and the next is greater than a threshold, the block boundary is taken to be an auditory event boundary. The audio segment between consecutive boundaries constitutes an auditory event. Thus, the third step sets an auditory event boundary between successive time segments when the difference in the spectral profile content between such successive time segments exceeds a threshold, thus defining auditory events. In this embodiment, auditory event boundaries define auditory events having a length that is an integral multiple of spectral profile blocks with a minimum length of one spectral profile block (512 samples in this example). In principle, event boundaries need not be so limited. As an alternative to the practical embodiments discussed herein, the input block size may vary, for example, so as to be essentially the size of an auditory event.
Following the identification of the event boundaries, key characteristics of the auditory event are identified, as shown in step 1-4.
Either overlapping or non-overlapping segments of the audio may be windowed and used to compute spectral profiles of the input audio. Overlap results in finer resolution as to the location of auditory events and, also, makes it less likely to miss an event, such as a short transient. However, overlap also increases computational complexity. Thus, overlap may be omitted. FIG. 2 shows a conceptual representation of non-overlapping N sample blocks being windowed and transformed into the frequency domain by the Discrete Fourier Transform (DFT). Each block may be windowed and transformed into the frequency domain, such as by using the DFT, preferably implemented as a Fast Fourier Transform (FFT) for speed.
The following variables may be used to compute the spectral profile of the input block:

    
    
        M=number of windowed samples in a block used to compute spectral profile
        P=number of samples of spectral computation overlap
    
    


In general, any integer numbers may be used for the variables above. However, the implementation will be more efficient if M is set equal to a power of 2 so that standard FFTs may be used for the spectral profile calculations. In a practical embodiment of the auditory scene analysis process, the parameters listed may be set to:

    
    
        M=512 samples (or 11.6 ms at 44.1 kHz)
        p=0 samples (no overlap)
    
    


The above-listed values were determined experimentally and were found generally to identify with sufficient accuracy the location and duration of auditory events. However, setting the value of P to 256 samples (50% overlap) rather than zero samples (no overlap) has been found to be useful in identifying some hard-to-find events. While many different types of windows may be used to minimize spectral artifacts due to windowing, the window used in the spectral profile calculations is an M-point Hanning, Kaiser-Bessel or other suitable, preferably non-rectangular, window. The above-indicated values and a Hanning window type were selected after extensive experimental analysis as they have shown to provide excellent results across a wide range of audio material. Non-rectangular windowing is preferred for the processing of audio signals with predominantly low frequency content Rectangular windowing produces spectral artifacts that may cause incorrect detection of events. Unlike certain encoder/decoder (codec) applications where an overall overlap/add process must provide a constant level, such a constraint does not apply here and the window may be chosen for characteristics such as its time/frequency resolution and stop-band rejection.
In step 1-1 (FIG. 1), the spectrum of each M-sample block may be computed by windowing the data with an M-point Hanning, Kaiser-Bessel or other suitable window, converting to the frequency domain using an M-point Fast Fourier Transform, and calculating the magnitude of the complex FFT coefficients. The resultant data is normalized so that the largest magnitude is set to unity, and the normalized array of M numbers is converted to the log domain. The data may also be normalized by some other metric such as the mean magnitude value or mean power value of the data. The array need not be converted to the log domain, but the conversion simplifies the calculation of the difference measure in step 1-2. Furthermore, the log domain more closely matches the nature of the human auditory system. The resulting log domain values have a range of minus infinity to zero. In a practical embodiment, a lower limit may be imposed on the range of values; the limit may be fixed, for example −60 dB, or be frequency-dependent to reflect the lower audibility of quiet sounds at low and very high frequencies. (Note that it would be possible to reduce the size of the array to M/2 in that the FFT represents negative as well as positive frequencies).
Step 1-2 calculates a measure of the difference between the spectra of adjacent blocks. For each block, each of the M (log) spectral coefficients from step 1-1 is subtracted from the corresponding coefficient for the preceding block, and the magnitude of the difference calculated (the sign is ignored). These M differences are then summed to one number. This difference measure may also be expressed as an average difference per spectral coefficient by dividing the difference measure by the number of spectral coefficients used in the sum (in this case M coefficients).
Step 1-3 identifies the locations of auditory event boundaries by applying a threshold to the array of difference measures from step 1-2 with a threshold value. When a difference measure exceeds a threshold, the change in spectrum is deemed sufficient to signal a new event and the block number of the change is recorded as an event boundary. For the values of M and P given above and for log domain values (in step 1-1) expressed in units of dB, the threshold may be set equal to 2500 if the whole magnitude FFT (including the mirrored part) is compared or 1250 if half the FFT is compared (as noted above, the FFT represents negative as well as positive frequencies—for the magnitude of the FFT, one is the mirror image of the other). This value was chosen experimentally and it provides good auditory event boundary detection. This parameter value may be changed to reduce (increase the threshold) or increase (decrease the threshold) the detection of events.
The process of FIG. 1 may be represented more generally by the equivalent arrangements of FIGS. 3, 4 and 5. In FIG. 3, an audio signal is applied in parallel to an “Identify Auditory Events” function or step 3-1 that divides the audio signal into auditory events, each of which tends to be perceived as separate and distinct and to an optional “Identify Characteristics of Auditory Events” function or step 3-2. The process of FIG. 1 may be employed to divide the audio signal into auditory events and their characteristics identified or some other suitable process may be employed. The auditory event information, which may be an identification of auditory event boundaries, determined by function or step 3-1 is then used to modify the audio dynamics processing parameters (such as attack, release, ratio, etc.), as desired, by a “Modify Dynamics Parameters” function or step 3-3. The optional “Identify Characteristics” function or step 3-3 also receives the auditory event information. The “Identify Characteristics” function or step 3-3 may characterize some or all of the auditory events by one or more characteristics. Such characteristics may include an identification of the dominant subband of the auditory event, as described in connection with the process of FIG. 1. The characteristics may also include one or more audio characteristics, including, for example, a measure of power of the auditory event, a measure of amplitude of the auditory event, a measure of the spectral flatness of the auditory event, and whether the auditory event is substantially silent, or other characteristics that help modify dynamics parameters such that negative audible artifacts of the processing are reduced or removed. The characteristics may also include other characteristics such as whether the auditory event includes a transient.
Alternatives to the arrangement of FIG. 3 are shown in FIGS. 4 and 5. In FIG. 4, the audio input signal is not applied directly to the “Identify Characteristics” function or step 4-3, but it does receive information from the “Identify Auditory Events” function or step 4-1. The arrangement of FIG. 1 is a specific example of such an arrangement. In FIG. 5, the functions or steps 5-1, 5-2 and 5-3 are arranged in series.
The details of this practical embodiment are not critical. Other ways to calculate the spectral content of successive time segments of the audio signal, calculate the differences between successive time segments, and set auditory event boundaries at the respective boundaries between successive time segments when the difference in the spectral profile content between such successive time segments exceeds a threshold may be employed.
Auditory Scene Analysis
New, Loudness Domain Method
International application under the Patent Cooperation Treaty S.N. PCT/US2005/038579, filed Oct. 25, 2005, published as International Publication Number WO 2006/047600 A1, entitled “Calculating and Adjusting the Perceived Loudness and/or the Perceived Spectral Balance of an Audio Signal” by Alan Jeffrey Seefeldt discloses, among other things, an objective measure of perceived loudness based on a psychoacoustic model. Said application is hereby incorporated by reference in its entirety. As described in said application, from an audio signal, x[n], an excitation signal E[b,t] is computed that approximates the distribution of energy along the basilar membrane of the inner ear at critical band b during time block t. This excitation may be computed from the Short-time Discrete Fourier Transform (STDFT) of the audio signal as follows:




  
    
      
        
          E
          ⁡
          
            [
            
              b
              ,
              t
            
            ]
          
        
        =
        
          
            
              λ
              b
            
            ⁢
            
              E
              ⁡
              
                [
                
                  b
                  ,
                  
                    t
                    -
                    1
                  
                
                ]
              
            
          
          +
          
            
              (
              
                1
                -
                
                  λ
                  b
                
              
              )
            
            ⁢
            
              
                ∑
                k
              
              ⁢
              
                
                  
                    
                    
                      T
                      ⁡
                      
                        [
                        k
                        ]
                      
                    
                    
                  
                  2
                
                ⁢
                
                  
                    
                    
                      
                        C
                        b
                      
                      ⁡
                      
                        [
                        k
                        ]
                      
                    
                    
                  
                  2
                
                ⁢
                
                  
                    
                    
                      X
                      ⁡
                      
                        [
                        
                          k
                          ,
                          t
                        
                        ]
                      
                    
                    
                  
                  2
                
              
            
          
        
      
    
    
      
        (
        1
        )
      
    
  




where X[k,t] represents the STDFT of x[n] at time block t and bin k. Note that in equation 1 t represents time in discrete units of transform blocks as opposed to a continuous measure, such as seconds. T[k] represents the frequency response of a filter simulating the transmission of audio through the outer and middle ear, and Cb[k] represents the frequency response of the basilar membrane at a location corresponding to critical band b. FIG. 6 depicts a suitable set of critical band filter responses in which 40 bands are spaced uniformly along the Equivalent Rectangular Bandwidth (ERB) scale, as defined by Moore and Glasberg. Each filter shape is described by a rounded exponential function and the bands are distributed using a spacing of 1 ERB. Lastly, the smoothing time constant λb in equation 1 may be advantageously chosen proportionate to the integration time of human loudness perception within band b.

Using equal loudness contours, such as those depicted in FIG. 7, the excitation at each band is transformed into an excitation level that would generate the same perceived loudness at 1 kHz. Specific loudness, a measure of perceptual loudness distributed across frequency and time, is then computed from the transformed excitation, E1kHz[b,t], through a compressive non-linearity. One such suitable function to compute the specific loudness N[b,t] is given by:




  
    
      
        
          N
          ⁡
          
            [
            
              b
              ,
              t
            
            ]
          
        
        =
        
          β
          ⁡
          
            (
            
              
                
                  (
                  
                    
                      
                        E
                        
                          1
                          ⁢
                          
                            
                          
                          ⁢
                          kHz
                        
                      
                      ⁡
                      
                        [
                        
                          b
                          ,
                          t
                        
                        ]
                      
                    
                    
                      TQ
                      
                        1
                        ⁢
                        
                          
                        
                        ⁢
                        kHz
                      
                    
                  
                  )
                
                a
              
              -
              1
            
            )
          
        
      
    
    
      
        (
        2
        )
      
    
  




where TQ1kHz is the threshold in quiet at 1 kHz and the constants β and α are chosen to match growth of loudness data as collected from listening experiments. Abstractly, this transformation from excitation to specific loudness may be presented by the function Ψ{ } such that:

N[b,t]=Ψ{E[b,t]}

Finally, the total loudness, L[t], represented in units of sone, is computed by summing the specific loudness across bands:





  
    
      
        
          L
          ⁡
          
            [
            t
            ]
          
        
        =
        
          
            ∑
            b
          
          ⁢
          
            N
            ⁡
            
              [
              
                b
                ,
                t
              
              ]
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




The specific loudness N[b,t] is a spectral representation meant to simulate the manner in which a human perceives audio as a function of frequency and time. It captures variations in sensitivity to different frequencies, variations in sensitivity to level, and variations in frequency resolution. As such, it is a spectral representation well matched to the detection of auditory events. Though more computationally complex, comparing the difference of N[b,t] across bands between successive time blocks may in many cases result in more perceptually accurate detection of auditory events in comparison to the direct use of successive FFT spectra described above.
In said patent application, several applications for modifying the audio based on this psychoacoustic loudness model are disclosed. Among these are several dynamics processing algorithms, such as AGC and DRC. These disclosed algorithms may benefit from the use of auditory events to control various associated parameters. Because specific loudness is already computed, it is readily available for the purpose of detecting said events. Details of a preferred embodiment are discussed below.
Audio Dynamics Processing Parameter Control with Auditory Events
Two examples of embodiments of the invention are now presented. The first describes the use of auditory events to control the release time in a digital implementation of a Dynamic Range Controller (DRC) in which the gain control is derived from the Root Mean Square (RMS) power of the signal. The second embodiment describes the use of auditory events to control certain aspects of a more sophisticated combination of AGC and DRC implemented within the context of the psychoacoustic loudness model described above. These two embodiments are meant to serve as examples of the invention only, and it should be understood that the use of auditory events to control parameters of a dynamics processing algorithm is not restricted to the specifics described below.
Dynamic Range Control
The described digital implementation of a DRC segments an audio signal x[n] into windowed, half-overlapping blocks, and for each block a modification gain based on a measure of the signal's local power and a selected compression curve is computed. The gain is smoothed across blocks and then multiplied with each block. The modified blocks are finally overlap-added to generate the modified audio signal y[n].
It should be noted, that while the auditory scene analysis and digital implementation of DRC as described here divides the time-domain audio signal into blocks to perform analysis and processing, the DRC processing need not be performed using block segmentation. For example the auditory scene analysis could be performed using block segmentation and spectral analysis as described above and the resulting auditory event locations and characteristics could be used to provide control information to a digital implementation of a traditional DRC implementation that typically operates on a sample-by-sample basis. Here, however, the same blocking structure used for auditory scene analysis is employed for the DRC to simplify the description of their combination.
Proceeding with the description of a block based DRC implementation, the overlapping blocks of the audio signal may be represented as:

x[n,t]=w[n]x[n+tM/2] for 0<n<M−1  (4)

where M is the block length and the hopsize is M/2, w[n] is the window, n is the sample index within the block, and t is the block index (note that here t is used in the same way as with the STDFT in equation 1; it represents time in discrete units of blocks rather than seconds, for example). Ideally, the window w[n] tapers to zero at both ends and sums to unity when half-overlapped with itself; the commonly used sine window meets these criteria, for example.

For each block, one may then compute the RMS power to generate a power measure P[t] in dB per block:




  
    
      
        
          P
          ⁡
          
            [
            t
            ]
          
        
        =
        
          10
          *
          log
          ⁢
          
            
          
          ⁢
          10
          ⁢
          
            
              (
              
                
                  1
                  M
                
                ⁢
                
                  
                    ∑
                    
                      n
                      =
                      1
                    
                    M
                  
                  ⁢
                  
                    
                      x
                      2
                    
                    ⁡
                    
                      [
                      
                        n
                        ,
                        t
                      
                      ]
                    
                  
                
              
              )
            
            .
          
        
      
    
    
      
        (
        5
        )
      
    
  




As mentioned earlier, one could smooth this power measure with a fast attack and slow release prior to processing with a compression curve, but as an alternative the instantaneous power P[t] is processed and the resulting gain is smoothed. This alternate approach has the advantage that a simple compression curve with sharp knee points may be used, but the resulting gains are still smooth as the power travels through the knee-point. Representing a compression curve as shown in FIG. 8c as a function F of signal level that generates a gain, the block gain G[t] is given by:

G[t]=F{P[t]}  (6)

Assuming that the compression curve applies greater attenuation as signal level increases, the gain will be decreasing when the signal is in “attack mode” and increasing when in “release mode”. Therefore, a smoothed gain G[t] may be computed according to:





  
    
      
        
          
            
              G
              _
            
            ⁡
            
              [
              t
              ]
            
          
          =
          
            
              
                α
                ⁡
                
                  [
                  t
                  ]
                
              
              ·
              
                
                  G
                  _
                
                ⁡
                
                  [
                  
                    t
                    -
                    1
                  
                  ]
                
              
            
            +
            
              
                (
                
                  1
                  -
                  
                    α
                    ⁡
                    
                      [
                      t
                      ]
                    
                  
                
                )
              
              ⁢
              
                G
                ⁡
                
                  [
                  t
                  ]
                
              
            
          
        
        ⁢
        
          

        
        ⁢
        where
      
    
    
      
        (
        
          7
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          α
          ⁡
          
            [
            t
            ]
          
        
        =
        
          {
          
            
              
                
                  
                    α
                    attach
                  
                
                
                  
                    
                      G
                      ⁡
                      
                        [
                        t
                        ]
                      
                    
                    <
                    
                      
                        G
                        _
                      
                      ⁡
                      
                        [
                        
                          t
                          -
                          1
                        
                        ]
                      
                    
                  
                
              
              
                
                  
                    α
                    release
                  
                
                
                  
                    
                      G
                      ⁡
                      
                        [
                        t
                        ]
                      
                    
                    ≥
                    
                      
                        
                          G
                          _
                        
                        ⁡
                        
                          [
                          
                            t
                            -
                            1
                          
                          ]
                        
                      
                      .
                    
                  
                
              
            
            ⁢
            
              

            
            ⁢
            and
          
        
      
    
    
      
        (
        
          7
          ⁢
          b
        
        )
      
    
  
  
    
      
        
          α
          release
        
        >>
        
          α
          attach
        
      
    
    
      
        (
        
          7
          ⁢
          c
        
        )
      
    
  




Finally, the smoothed gain G[t], which is in dB, is applied to each block of the signal, and the modified blocks are overlap-added to produce the modified audio:

y[n+tM/2]=(10 G[t]/20)x[n,t]+(10 G[t−1]/20)x[n+M/2,t−1] for 0<n<M/2  (8)

Note that because the blocks have been multiplied with a tapered window, as shown in equation 4, the overlap-add synthesis shown above effectively smooths the gains across samples of the processed signal y[n]. Thus, the gain control signal receives smoothing in addition to that in shown in equation 7a. In a more traditional implementation of DRC operating sample-by-sample rather than block-by-block, gain smoothing more sophisticated than the simple one-pole filter shown in equation 7a might be necessary in order to prevent audible distortion in the processed signal. Also, the use of block based processing introduces an inherent delay of M/2 samples into the system, and as long as the decay time associated with αattack is close to this delay, the signal x[n] does not need to be delayed further before the application of the gains for the purposes of preventing overshoot.

FIGS. 9a through 9c depict the result of applying the described DRC processing to an audio signal. For this particular implementation, a block length of M=512 is used at a sampling rate of 44.1 kHz. A compression curve similar to the one shown in FIG. 8b is used: above −20 dB relative to full scale digital the signal is attenuated with a ratio of 5:1, and below −30 dB the signal is boosted with a ratio of 5:1. The gain is smoothed with an attack coefficient αattack corresponding to a half-decay time of 10 ms and a release coefficient αrelease corresponding to a half-decay time of 500 ms. The original audio signal depicted in FIG. 9a consists of six consecutive piano chords, with the final chord, located around sample 1.75×105, decaying into silence. Examining a plot of the gain G[t] in FIG. 9b, it should be noted that the gain remains close to 0 dB while the six chords are played. This is because the signal energy remains, for the most part, between −30 dB and −20 dB, the region within which the DRC curve calls for no modification. However, after the hit of the last chord, the signal energy falls below −30 dB, and the gain begins to rise, eventually beyond 15 dB, as the chord decays. FIG. 9c depicts the resulting modified audio signal, and one can see that the tail of the final chord is boosted significantly. Audibly, this boosting of the chord's natural, low-level decay sound creates an extremely unnatural result. It is the aim of the present invention to prevent problems of this type that are associated with a traditional dynamics processor.
FIGS. 10a through 10c depict the results of applying the exact same DRC system to a different audio signal. In this case the first half of the signal consists of an up-tempo music piece at a high level, and then at approximately sample 10×104 the signal switches to a second up-tempo music piece, but at a significantly lower level. Examining the gain in FIG. 6b, one sees that the signal is attenuated by approximately 10 dB during the first half, and then the gain rises back up to 0 dB during the second half when the softer piece is playing. In this case, the gain behaves as desired. One would like the second piece to be boosted relative to the first, and the gain should increase quickly after the transition to the second piece to be audibly unobtrusive. One sees a gain behavior that is similar to that for the first signal discussed, but here the behavior is desirable. Therefore, one would like to fix the first case without affecting the second. The use of auditory events to control the release time of this DRC system provides such a solution.
In the first signal that was examined in FIG. 9, the boosting of the last chord's decay seems unnatural because the chord and its decay are perceived as a single auditory event whose integrity is expected to be maintained. In the second case, however, many auditory events occur while the gain increases, meaning that for any individual event, little change is imparted. Therefore the overall gain change is not as objectionable. One may therefore argue that a gain change should be allowed only in the near temporal vicinity of an auditory event boundary. One could apply this principal to the gain while it is in either attack or release mode, but for most practical implementations of a DRC, the gain moves so quickly in attack mode in comparison to the human temporal resolution of event perception that no control is necessary. One may therefore use events to control smoothing of the DRC gain only when it is in release mode.
A suitable behavior of the release control is now described. In qualitative terms, if an event is detected, the gain is smoothed with the release time constant as specified above in Equation 7a. As time evolves past the detected event, and if no subsequent events are detected, the release time constant continually increases so that eventually the smoothed gain is “frozen” in place. If another event is detected, then the smoothing time constant is reset to the original value and the process repeats. In order to modulate the release time, one may first generate a control signal based on the detected event boundaries.
As discussed earlier, event boundaries may be detected by looking for changes in successive spectra of the audio signal. In this particular implementation, the DFT of each overlapping block x[n,t] may be computed to generate the STDFT of the audio signal x[n]:




  
    
      
        
          X
          ⁡
          
            [
            
              k
              ,
              t
            
            ]
          
        
        =
        
          
            ∑
            
              n
              =
              0
            
            
              M
              -
              1
            
          
          ⁢
          
            
              x
              ⁡
              
                [
                
                  n
                  ,
                  t
                
                ]
              
            
            ⁢
            
              e
              
                
                  -
                  j
                
                ⁢
                
                  
                    2
                    ⁢
                    π
                    ⁢
                    
                      
                    
                    ⁢
                    kn
                  
                  M
                
              
            
          
        
      
    
    
      
        (
        9
        )
      
    
  




Next, the difference between the normalized log magnitude spectra of successive blocks may be computed according to:





  
    
      
        
          
            D
            ⁡
            
              [
              t
              ]
            
          
          =
          
            
              ∑
              k
            
            ⁢
            
              
              
                
                  
                    X
                    NORM
                  
                  ⁡
                  
                    [
                    
                      k
                      ,
                      t
                    
                    ]
                  
                
                -
                
                  
                    X
                    NORM
                  
                  ⁡
                  
                    [
                    
                      k
                      ,
                      
                        t
                        -
                        1
                      
                    
                    ]
                  
                
              
              
            
          
        
        ⁢
        
          

        
        ⁢
        where
      
    
    
      
        (
        
          10
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            X
            NORM
          
          ⁡
          
            [
            
              k
              ,
              t
            
            ]
          
        
        =
        
          log
          ⁡
          
            (
            
              
                
                
                  X
                  ⁡
                  
                    [
                    
                      k
                      ,
                      t
                    
                    ]
                  
                
                
              
              
                
                  max
                  k
                
                ⁢
                
                  {
                  
                    
                    
                      X
                      ⁡
                      
                        [
                        
                          k
                          ,
                          t
                        
                        ]
                      
                    
                    
                  
                  }
                
              
            
            )
          
        
      
    
    
      
        (
        
          10
          ⁢
          b
        
        )
      
    
  




Here the maximum of |X[k,t]| across bins k is used for normalization, although one might employ other normalization factors; for example, the average of |X[k,t]| across bins. If the difference D[t] exceeds a threshold Dmin, then an event is considered to have occurred. Additionally, one may assign a strength to this event, lying between zero and one, based on the size of D[t] in comparison to a maximum threshold Dmax. The resulting auditory event strength signal A[t] may be computed as:





  
    
      
        
          A
          ⁡
          
            [
            t
            ]
          
        
        =
        
          {
          
            
              
                
                  
                  ⁢
                  0
                
              
              
                
                  
                  ⁢
                  
                    
                      D
                      ⁡
                      
                        [
                        t
                        ]
                      
                    
                    ≤
                    
                      D
                      min
                    
                  
                
              
            
            
              
                
                  
                    
                      D
                      ⁡
                      
                        [
                        t
                        ]
                      
                    
                    -
                    
                      D
                      
                        mi
                        ⁢
                        n
                      
                    
                  
                  
                    
                      D
                      max
                    
                    -
                    
                      D
                      min
                    
                  
                
              
              
                
                  
                    D
                    min
                  
                  <
                  
                    D
                    ⁡
                    
                      [
                      t
                      ]
                    
                  
                  <
                  
                    D
                    max
                  
                
              
            
            
              
                
                  
                  ⁢
                  1
                
              
              
                
                  
                  ⁢
                  
                    
                      D
                      ⁡
                      
                        [
                        t
                        ]
                      
                    
                    ≥
                    
                      D
                      max
                    
                  
                
              
            
          
        
      
    
    
      
        (
        11
        )
      
    
  




By assigning a strength to the auditory event proportional to the amount of spectral change associated with that event, greater control over the dynamics processing is achieved in comparison to a binary event decision. The inventors have found that larger gain changes are acceptable during stronger events, and the signal in equation 11 allows such variable control.

The signal A[t] is an impulsive signal with an impulse occurring at the location of an event boundary. For the purposes of controlling the release time, one may further smooth the signal A[t] so that it decays smoothly to zero after the detection of an event boundary. The smoothed event control signal Ā[t] may be computed from A[t] according to:




  
    
      
        
          
            A
            _
          
          ⁡
          
            [
            t
            ]
          
        
        =
        
          {
          
            
              
                
                  
                  ⁢
                  
                    A
                    ⁡
                    
                      [
                      t
                      ]
                    
                  
                
              
              
                
                  
                  ⁢
                  
                    
                      A
                      ⁡
                      
                        [
                        t
                        ]
                      
                    
                    >
                    
                      
                        α
                        event
                      
                      ⁢
                      
                        
                          A
                          _
                        
                        ⁡
                        
                          [
                          
                            t
                            -
                            1
                          
                          ]
                        
                      
                    
                  
                
              
            
            
              
                
                  
                    α
                    event
                  
                  ⁢
                  
                    
                      A
                      _
                    
                    ⁡
                    
                      [
                      
                        t
                        -
                        1
                      
                      ]
                    
                  
                
              
              
                
                  
                  ⁢
                  otherwise
                
              
            
          
        
      
    
    
      
        (
        12
        )
      
    
  




Here αevent controls the decay time of the event control signal. FIGS. 9d and 10d depict the event control signal Ā[t] for the two corresponding audio signals, with the half-decay time of the smoother set to 250 ms. In the first case, one sees that an event boundary is detected for each of the six piano chords, and that the event control signal decays smoothly towards zero after each event. For the second signal, many events are detected very close to each other in time, and therefore the event control signal never decays fully to zero.

One may now use the event control signal Ā[t] to vary the release time constant used for smoothing the gain. When the control signal is equal to one, the smoothing coefficient α[t] from Equation 7a equals αrelease, as before, and when the control signal is equal to zero, the coefficient equals one so that the smoothed gain is prevented from changing. The smoothing coefficient is interpolated between these two extremes using the control signal according to:




  
    
      
        
          α
          ⁡
          
            [
            t
            ]
          
        
        =
        
          {
          
            
              
                
                  
                  ⁢
                  
                    α
                    attack
                  
                
              
              
                
                  
                    G
                    ⁡
                    
                      [
                      t
                      ]
                    
                  
                  <
                  
                    
                      G
                      _
                    
                    ⁡
                    
                      [
                      
                        t
                        -
                        1
                      
                      ]
                    
                  
                
              
            
            
              
                
                  
                    
                      
                        
                          
                          ⁢
                          A
                        
                        _
                      
                      ⁡
                      
                        [
                        t
                        ]
                      
                    
                    ⁢
                    
                      α
                      release
                    
                  
                  +
                  
                    (
                    
                      1
                      -
                      
                        
                          A
                          _
                        
                        ⁡
                        
                          [
                          t
                          ]
                        
                      
                    
                    )
                  
                
              
              
                
                  
                    G
                    ⁡
                    
                      [
                      t
                      ]
                    
                  
                  ≥
                  
                    
                      G
                      _
                    
                    ⁡
                    
                      [
                      
                        t
                        -
                        1
                      
                      ]
                    
                  
                
              
            
          
        
      
    
    
      
        (
        13
        )
      
    
  




By interpolating the smoothing coefficient continuously as a function of the event control signal, the release time is reset to a value proportionate to the event strength at the onset of an event and then increases smoothly to infinity after the occurrence of an event. The rate of this increase is dictated by the coefficient αevent used to generate the smoothed event control signal.

FIGS. 9e and 10e show the effect of smoothing the gain with the event-controlled coefficient from Equation 13 as opposed to non-event-controlled coefficient from Equation 7b. In the first case, the event control signal falls to zero after the last piano chord, thereby preventing the gain from moving upwards. As a result, the corresponding modified audio in FIG. 9f does not suffer from an unnatural boost of the chord's decay. In the second case, the event control signal never approaches zero, and therefore the smoothed gain signal is inhibited very little through the application of the event control. The trajectory of the smoothed gain is nearly identical to the non-event-controlled gain in FIG. 10b. This is exactly the desired effect.
Loudness Based AGC and DRC
As an alternative to traditional dynamics processing techniques where signal modifications are a direct function of simple signal measurements such as Peak or RMS power, International Patent Application S.N. PCT/US2005/038579 discloses use of the psychoacoustic based loudness model described earlier as a framework within which to perform dynamics processing. Several advantages are cited. First, measurements and modifications are specified in units of sone, which is a more accurate measure of loudness perception than more basic measures such as Peak or RMS power. Secondly, the audio may be modified such that the perceived spectral balance of the original audio is maintained as the overall loudness is changed. This way, changes to the overall loudness become less perceptually apparent in comparison to a dynamics processor that utilizes a wideband gain, for example, to modify the audio. Lastly, the psychoacoustic model is inherently multi-band, and therefore the system is easily configured to perform multi-band dynamics processing in order to alleviate the well-known cross-spectral pumping problems associated with a wideband dynamics processor.
Although performing dynamics processing in this loudness domain already holds several advantages over more traditional dynamics processing, the technique may be further improved through the use of auditory events to control various parameters. Consider the audio segment containing piano chords as depicted in 27a and the associated DRC shown in FIGS. 10b and c. One could perform a similar DRC in the loudness domain, and in this case, when the loudness of the final piano chord's decay is boosted, the boost would be less apparent because the spectral balance of the decaying note would be maintained as the boost is applied. However, a better solution is to not boost the decay at all, and therefore one may advantageously apply the same principle of controlling attack and release times with auditory events in the loudness domain as was previously described for the traditional DRC.
The loudness domain dynamics processing system that is now described consists of AGC followed by DRC. The goal of this combination is to make all processed audio have approximately the same perceived loudness while still maintaining at least some of the original audio's dynamics. FIG. 11 depicts a suitable set of AGC and DRC curves for this application. Note that the input and output of both curves is represented in units of sone since processing is performed in the loudness domain. The AGC curve strives to bring the output audio closer to some target level, and, as mentioned earlier, does so with relatively slow time constants. One may think of the AGC as making the long-term loudness of the audio equal to the target, but on a short-term basis, the loudness may fluctuate significantly around this target. Therefore, one may employ faster acting DRC to limit these fluctuations to some range deemed acceptable for the particular application. FIG. 11 shows such a DRC curve where the AGC target falls within the “null band” of the DRC, the portion of the curve that calls for no modification. With this combination of curves, the AGC places the long-term loudness of the audio within the null-band of the DRC curve so that minimal fast-acting DRC modifications need be applied. If the short-term loudness still fluctuates outside of the null-band, the DRC then acts to move the loudness of the audio towards this null-band. As a final general note, one may apply the slow acting AGC such that all bands of the loudness model receive the same amount of loudness modification, thereby maintaining the perceived spectral balance, and one may apply the fast acting DRC in a manner that allows the loudness modification to vary across bands in order alleviate cross-spectral pumping that might otherwise result from fast acting band-independent loudness modification.
Auditory events may be utilized to control the attack and release of both the AGC and DRC. In the case of AGC, both the attack and release times are large in comparison to the temporal resolution of event perception, and therefore event control may be advantageously employed in both cases. With the DRC, the attack is relatively short, and therefore event control may be needed only for the release as with the traditional DRC described above.
As discussed earlier, one may use the specific loudness spectrum associated with the employed loudness model for the purposes of event detection. A difference signal D[t], similar to the one in Equations 10a and b may be computed from the specific loudness N[b,t], defined in Equation 2, as follows:




  
    
      
        
          
            D
            ⁡
            
              [
              t
              ]
            
          
          =
          
            
              ∑
              b
            
            ⁢
            
              
              
                
                  
                    N
                    NORM
                  
                  ⁡
                  
                    [
                    
                      b
                      ,
                      t
                    
                    ]
                  
                
                -
                
                  
                    N
                    NORM
                  
                  ⁡
                  
                    [
                    
                      b
                      ,
                      
                        t
                        -
                        1
                      
                    
                    ]
                  
                
              
              
            
          
        
        ⁢
        
          

        
        ⁢
        where
      
    
    
      
        (
        
          14
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            N
            NORM
          
          ⁡
          
            [
            
              b
              ,
              t
            
            ]
          
        
        =
        
          
            N
            ⁡
            
              [
              
                b
                ,
                t
              
              ]
            
          
          
            
              max
              b
            
            ⁢
            
              {
              
                N
                ⁡
                
                  [
                  
                    b
                    ,
                    t
                  
                  ]
                
              
              }
            
          
        
      
    
    
      
        (
        
          14
          ⁢
          b
        
        )
      
    
  




Here the maximum of |N[b,t]| across frequency bands b is used for normalization, although one might employ other normalization factors; for example, the average of |N[b,t]| across frequency bands. If the difference D[t] exceeds a threshold Dmin, then an event is considered to have occurred. The difference signal may then be processed in the same way shown in Equations 11 and 12 to generate a smooth event control signal Ā[t] used to control the attack and release times.

The AGC curve depicted in FIG. 11 may be represented as a function that takes as its input a measure of loudness and generates a desired output loudness:

Lo=FAGC{L1}  (15a)

The DRC curve may be similarly represented:

Lo=FDRC{L1}  (15b)

For the AGC, the input loudness is a measure of the audio's long-term loudness. One may compute such a measure by smoothing the instantaneous loudness L[t], defined in Equation 3, using relatively long time constants (on the order of several seconds). It has been shown that in judging an audio segment's long term loudness, humans weight the louder portions more heavily than the softer, and one may use a faster attack than release in the smoothing to simulate this effect. With the incorporation of event control for both the attack and release, the long-term loudness used for determining the AGC modification may therefore be computed according to:





  
    
      
        
          
            
              L
              AGC
            
            ⁡
            
              [
              t
              ]
            
          
          =
          
            
              
                
                  α
                  AGC
                
                ⁡
                
                  [
                  t
                  ]
                
              
              ⁢
              
                
                  L
                  AGC
                
                ⁡
                
                  [
                  
                    t
                    -
                    1
                  
                  ]
                
              
            
            +
            
              
                (
                
                  1
                  -
                  
                    
                      α
                      AGC
                    
                    ⁡
                    
                      [
                      t
                      ]
                    
                  
                
                )
              
              ⁢
              
                L
                ⁡
                
                  [
                  t
                  ]
                
              
            
          
        
        ⁢
        
          

        
        ⁢
        where
      
    
    
      
        (
        
          16
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            α
            AGC
          
          ⁡
          
            [
            t
            ]
          
        
        =
        
          {
          
            
              
                
                  
                    
                      
                        A
                        _
                      
                      ⁡
                      
                        [
                        t
                        ]
                      
                    
                    ⁢
                    
                      α
                      AGCattach
                    
                  
                  +
                  
                    (
                    
                      1
                      -
                      
                        
                          A
                          _
                        
                        ⁡
                        
                          [
                          t
                          ]
                        
                      
                    
                    )
                  
                
              
              
                
                  
                    L
                    ⁡
                    
                      [
                      t
                      ]
                    
                  
                  >
                  
                    
                      L
                      AGC
                    
                    ⁡
                    
                      [
                      
                        t
                        -
                        1
                      
                      ]
                    
                  
                
              
            
            
              
                
                  
                    
                      
                        A
                        _
                      
                      ⁡
                      
                        [
                        t
                        ]
                      
                    
                    ⁢
                    
                      α
                      AGCrelease
                    
                  
                  +
                  
                    (
                    
                      1
                      -
                      
                        
                          A
                          _
                        
                        ⁡
                        
                          [
                          t
                          ]
                        
                      
                    
                    )
                  
                
              
              
                
                  
                    L
                    ⁡
                    
                      [
                      t
                      ]
                    
                  
                  ≤
                  
                    
                      
                        L
                        AGC
                      
                      ⁡
                      
                        [
                        
                          t
                          -
                          1
                        
                        ]
                      
                    
                    .
                  
                
              
            
          
        
      
    
    
      
        (
        
          16
          ⁢
          b
        
        )
      
    
  




In addition, one may compute an associated long-term specific loudness spectrum that will later be used for the multi-band DRC:

NAGC[b,t]=αAGC[t]NAGC[b,t−1]+(1−αAGC[t])N[b,t]  (16c)

In practice one may choose the smoothing coefficients such that the attack time is approximately half that of the release. Given the long-term loudness measure, one may then compute the loudness modification scaling associated with the AGC as the ratio of the output loudness to input loudness:





  
    
      
        
          
            S
            AGC
          
          ⁡
          
            [
            t
            ]
          
        
        =
        
          
            
              F
              AGC
            
            ⁢
            
              {
              
                
                  L
                  AGC
                
                ⁡
                
                  [
                  t
                  ]
                
              
              }
            
          
          
            
              L
              AGC
            
            ⁡
            
              [
              t
              ]
            
          
        
      
    
    
      
        (
        17
        )
      
    
  




The DRC modification may now be computed from the loudness after the application of the AGC scaling. Rather than smooth a measure of the loudness prior to the application of the DRC curve, one may alternatively apply the DRC curve to the instantaneous loudness and then subsequently smooth the resulting modification. This is similar to the technique described earlier for smoothing the gain of the traditional DRC. In addition, the DRC may be applied in a multi-band fashion, meaning that the DRC modification is a function of the specific loudness N[b,t] in each band b, rather than the overall loudness L[t]. However, in order to maintain the average spectral balance of the original audio, one may apply DRC to each band such that the resulting modifications have the same average effect as would result from applying DRC to the overall loudness. This may be achieved by scaling each band by the ratio of the long-term overall loudness (after the application of the AGC scaling) to the long-term specific loudness, and using this value as the argument to the DRC function. The result is then rescaled by the inverse of said ratio to produce the output specific loudness. Thus, the DRC scaling in each band may be computed according to:




  
    
      
        
          
            S
            DRC
          
          ⁡
          
            [
            
              b
              ,
              t
            
            ]
          
        
        =
        
          
            
              
                
                  N
                  AGC
                
                ⁡
                
                  [
                  
                    b
                    ,
                    t
                  
                  ]
                
              
              
                
                  
                    S
                    AGC
                  
                  ⁡
                  
                    [
                    t
                    ]
                  
                
                ⁢
                
                  
                    L
                    AGC
                  
                  ⁡
                  
                    [
                    t
                    ]
                  
                
              
            
            ⁢
            
              F
              DRC
            
            ⁢
            
              {
              
                
                  
                    
                      
                        S
                        AGC
                      
                      ⁡
                      
                        [
                        t
                        ]
                      
                    
                    ⁢
                    
                      
                        L
                        AGC
                      
                      ⁡
                      
                        [
                        t
                        ]
                      
                    
                  
                  
                    
                      N
                      AGC
                    
                    ⁡
                    
                      [
                      t
                      ]
                    
                  
                
                ⁢
                
                  N
                  ⁡
                  
                    [
                    
                      b
                      ,
                      t
                    
                    ]
                  
                
              
              }
            
          
          
            N
            ⁡
            
              [
              
                b
                ,
                t
              
              ]
            
          
        
      
    
    
      
        (
        18
        )
      
    
  




The AGC and DRC modifications may then be combined to form a total loudness scaling per band:

STOT[b,t]=SAGC[t]SDRC[b,t]  (19)

This total scaling may then be smoothed across time independently for each band with a fast attack and slow release and event control applied to the release only. Ideally smoothing is performed on the logarithm of the scaling analogous to the gains of the traditional DRC being smoothed in their decibel representation, though this is not essential. To ensure that the smoothed total scaling moves in sync with the specific loudness in each band, attack and release modes may by determined through the simultaneous smoothing of specific loudness itself:





  
    
      
        
          
            
              S
              _
            
            TOT
          
          ⁡
          
            [
            
              b
              ,
              t
            
            ]
          
        
        =
        
          exp
          ⁡
          
            (
            
              
                
                  
                    
                      
                        
                          α
                          TOT
                        
                        ⁡
                        
                          [
                          
                            b
                            ,
                            t
                          
                          ]
                        
                      
                      ⁢
                      
                        log
                        ⁡
                        
                          (
                          
                            
                              
                                S
                                _
                              
                              TOT
                            
                            ⁡
                            
                              [
                              
                                b
                                ,
                                
                                  t
                                  -
                                  1
                                
                              
                              ]
                            
                          
                          )
                        
                      
                    
                    +
                  
                
              
              
                
                  
                    
                      (
                      
                        1
                        -
                        
                          
                            α
                            TOT
                          
                          ⁡
                          
                            [
                            
                              b
                              ,
                              t
                            
                            ]
                          
                        
                      
                      )
                    
                    ⁢
                    
                      log
                      ⁡
                      
                        (
                        
                          
                            S
                            TOT
                          
                          ⁡
                          
                            [
                            
                              b
                              ,
                              t
                            
                            ]
                          
                        
                        )
                      
                    
                  
                
              
            
            )
          
        
      
    
    
      
        (
        
          20
          ⁢
          a
        
        )
      
    
  
  
    
      
        
          
            
              N
              _
            
            ⁡
            
              [
              
                b
                ,
                t
              
              ]
            
          
          =
          
            
              
                
                  α
                  TOT
                
                ⁡
                
                  [
                  
                    b
                    ,
                    t
                  
                  ]
                
              
              ⁢
              
                
                  N
                  _
                
                ⁡
                
                  [
                  
                    b
                    ,
                    
                      t
                      -
                      1
                    
                  
                  ]
                
              
            
            +
            
              
                (
                
                  1
                  -
                  
                    
                      α
                      TOT
                    
                    ⁡
                    
                      [
                      
                        b
                        ,
                        t
                      
                      ]
                    
                  
                
                )
              
              ⁢
              
                N
                ⁡
                
                  [
                  
                    b
                    ,
                    t
                  
                  ]
                
              
            
          
        
        ⁢
        
          

        
        ⁢
        where
      
    
    
      
        (
        
          20
          ⁢
          b
        
        )
      
    
  
  
    
      
        
          
            α
            TOT
          
          ⁡
          
            [
            
              b
              ,
              t
            
            ]
          
        
        =
        
          {
          
            
              
                
                  
                  ⁢
                  
                    α
                    TOTattack
                  
                
              
              
                
                  
                    N
                    ⁡
                    
                      [
                      
                        b
                        ,
                        t
                      
                      ]
                    
                  
                  >
                  
                    
                      N
                      _
                    
                    ⁡
                    
                      [
                      
                        b
                        ,
                        
                          t
                          -
                          1
                        
                      
                      ]
                    
                  
                
              
            
            
              
                
                  
                    
                      
                        
                          
                          ⁢
                          A
                        
                        _
                      
                      ⁡
                      
                        [
                        t
                        ]
                      
                    
                    ⁢
                    
                      α
                      TOTrelease
                    
                  
                  +
                  
                    (
                    
                      1
                      -
                      
                        
                          A
                          _
                        
                        ⁡
                        
                          [
                          t
                          ]
                        
                      
                    
                    )
                  
                
              
              
                
                  
                    N
                    ⁡
                    
                      [
                      
                        b
                        ,
                        t
                      
                      ]
                    
                  
                  ≤
                  
                    
                      N
                      _
                    
                    ⁡
                    
                      [
                      
                        b
                        ,
                        
                          t
                          -
                          1
                        
                      
                      ]
                    
                  
                
              
            
          
        
      
    
    
      
        (
        
          20
          ⁢
          c
        
        )
      
    
  




Finally one may compute a target specific loudness based on the smoothed scaling applied to the original specific loudness

{circumflex over (N)}[b,t]= STOT[b,t]N[b,t]  (21)

and then solve for gains G[b,t] that when applied to the original excitation result in a specific loudness equal to the target:

{circumflex over (N)}[b,t]=Ψ{G2[b,t]E[b,t]}  (22)

The gains may be applied to each band of the filterbank used to compute the excitation, and the modified audio may then be generated by inverting the filterbank to produce a modified time domain audio signal.

Additional Parameter Control
While the discussion above has focused on the control of AGC and DRC attack and release parameters via auditory scene analysis of the audio being processed, other important parameters may also benefit from being controlled via the ASA results. For example, the event control signal Ā[t] from Equation 12 may be used to vary the value of the DRC ratio parameter that is used to dynamically adjust the gain of the audio. The Ratio parameter, similarly to the attack and release time parameters, may contribute significantly to the perceptual artifacts introduced by dynamic gain adjustments.
Implementation
The invention may be implemented in hardware or software, or a combination of both (e.g., programmable logic arrays). Unless otherwise specified, the algorithms included as part of the invention are not inherently related to any particular computer or other apparatus. In particular, various general-purpose machines may be used with programs written in accordance with the teachings herein, or it may be more convenient to construct more specialized apparatus (e.g., integrated circuits) to perform the required method steps. Thus, the invention may be implemented in one or more computer programs executing on one or more programmable computer systems each comprising at least one processor, at least one data storage system (including volatile and non-volatile memory and/or storage elements), at least one input device or port, and at least one output device or port. Program code is applied to input data to perform the functions described herein and generate output information. The output information is applied to one or more output devices, in known fashion.
Each such program may be implemented in any desired computer language (including machine, assembly, or high level procedural, logical, or object oriented programming languages) to communicate with a computer system. In any case, the language may be a compiled or interpreted language.
Each such computer program is preferably stored on or downloaded to a storage media or device (e.g., solid state memory or media, or magnetic or optical media) readable by a general or special purpose programmable computer, for configuring and operating the computer when the storage media or device is read by the computer system to perform the procedures described herein. The inventive system may also be considered to be implemented as a computer-readable storage medium, configured with a computer program, where the storage medium so configured causes a computer system to operate in a specific and predefined manner to perform the functions described herein.
A number of embodiments of the invention have been described. Nevertheless, it will be understood that various modifications may be made without departing from the spirit and scope of the invention. For example, some of the steps described herein may be order independent, and thus may be performed in an order different from that described.
It should be understood that implementation of other variations and modifications of the invention and its various aspects will be apparent to those skilled in the art, and that the invention is not limited by these specific embodiments described. It is therefore contemplated to cover by the present invention any and all modifications, variations, or equivalents that fall within the true spirit and scope of the basic underlying principles disclosed and claimed herein.
INCORPORATION BY REFERENCE
The following patents, patent applications and publications are hereby incorporated by reference, each in their entirety.
Audio Dynamics Processing


    Audio Engineer's Reference Book, edited by Michael Talbot-Smith, 2nd edition. Limiters and Compressors, Alan Tutton, 2-1492-165. Focal Press, Reed Educational and Professional Publishing, Ltd., 1999.


Detecting and Using Auditory Events


    U.S. patent application Ser. No. 10/474,387, “High Quality Time-Scaling and Pitch-Scaling of Audio Signals” of Brett Graham Crockett, published Jun. 24, 2004 as US 2004/0122662 A1.
    U.S. patent application Ser. No. 10/478,398, “Method for Time Aligning Audio Signals Using Characterizations Based on Auditory Events” of Brett G. Crockett et al, published Jul. 29, 2004 as US 2004/0148159 A1.
    U.S. patent application Ser. No. 10/478,538, “Segmenting Audio Signals Into Auditory Events” of Brett G. Crockett, published Aug. 26, 2004 as US 2004/0165730 A1. Aspects of the present invention provide a way to detect auditory events in addition to those disclosed in said application of Crockett.
    U.S. patent application Ser. No. 10/478,397, “Comparing Audio Using Characterizations Based on Auditory Events” of Brett G. Crockett et al, published Sep. 2, 2004 as US 2004/0172240 A1.
    International Application under the Patent Cooperation Treaty S.N. PCT/US 05/24630 filed Jul. 13, 2005, entitled “Method for Combining Audio Signals Using Auditory Scene Analysis,” of Michael John Smithers, published Mar. 9, 2006 as WO 2006/026161.
    International Application under the Patent Cooperation Treaty S.N. PCT/US 2004/016964, filed May 27, 2004, entitled “Method, Apparatus and Computer Program for Calculating and Adjusting the Perceived Loudness of an Audio Signal” of Alan Jeffrey Seefeldt et al, published Dec. 23, 2004 as WO 2004/111994 A2.
    International application under the Patent Cooperation Treaty S.N. PCT/US2005/038579, filed Oct. 25, 2005, entitled “Calculating and Adjusting the Perceived Loudness and/or the Perceived Spectral Balance of an Audio Signal” by Alan Jeffrey Seefeldt and published as International Publication Number WO 2006/047600.
    “A Method for Characterizing and Identifying Audio Based on Auditory Scene Analysis,” by Brett Crockett and Michael Smithers, Audio Engineering Society Convention Paper 6416, 118th Convention, Barcelona, May 28-31, 2005.
    “High Quality Multichannel Time Scaling and Pitch-Shifting using Auditory Scene Analysis,” by Brett Crockett, Audio Engineering Society Convention Paper 5948, New York, October 2003.
    “A New Objective Measure of Perceived Loudness” by Alan Seefeldt et al, Audio Engineering Society Convention Paper 6236, San Francisco, Oct. 28, 2004.
    Handbook for Sound Engineers, The New Audio Cyclopedia, edited by Glen M. Ballou, 2nd edition. Dynamics, 850-851. Focal Press an imprint of Butterworth-Heinemann, 1998.
    Audio Engineer's Reference Book, edited by Michael Talbot-Smith, 2nd edition, Section 2.9 (“Limiters and Compressors” by Alan Tutton), pp. 2.149-2.165, Focal Press, Reed Educational and Professional Publishing, Ltd., 1999.

## Claims

1. A method for modifying a parameter of an audio dynamics processor, comprising
detecting changes in spectral characteristics with respect to time in an audio signal,
identifying as auditory event boundaries changes greater than a threshold in spectral characteristics with respect to time in said audio signal, wherein an audio segment between consecutive boundaries constitutes an auditory event,
generating a parameter-modifying dynamically-varying control signal based on said identified event boundaries, and
modifying the parameter of the audio dynamics processor as a function of the control signal.

2. A method according to claim 1 wherein the parameter is one of attack time, release time, and ratio.

3. A method according to claim 1 wherein the parameter modified is a gain-smoothing time constant.

4. A method according to claim 3 wherein the gain-smoothing time constant is a gain-smoothing attack time constant.

5. A method according to claim 3 wherein the gain-smoothing time constant is a gain-smoothing release time constant.

6. A method according to claim 1 wherein said parameter-modifying control signal is based on the location of said identified auditory event boundaries and the degree of change in spectral characteristics associated with each of said auditory event boundaries.

7. A method according to claim 6 wherein generating a parameter modifying control signal comprises
providing an impulse at each of auditory event boundaries, each such impulse having an amplitude proportional to the degree of said changes in spectral characteristics, and
time smoothing each such impulse so that its amplitude decays smoothly toward zero.

8. A method according to claim 1 wherein changes in spectral characteristics with respect to time are detected by comparing differences in specific loudness.

9. A method according to claim 8 wherein said audio signal is represented by a discrete time sequence x[n] that has been sampled from an audio source at a sampling frequency fs and the changes in spectral characteristics with respect to time are calculated by comparing the difference in specific loudness N[b,t] across frequency bands b between successive time blocks t.

10. A method according to claim 9 wherein the difference in spectral content between successive time blocks of the audio signal is calculated according to



  
    
      D
      ⁡
      
        [
        t
        ]
      
    
    =
    
      
        ∑
        b
      
      ⁢
      
        
        
          
            
              N
              NORM
            
            ⁡
            
              [
              
                b
                ,
                t
              
              ]
            
          
          -
          
            
              N
              NORM
            
            ⁡
            
              [
              
                b
                ,
                
                  t
                  -
                  1
                
              
              ]
            
          
        
        
      
    
  




  where




  
    
      
        N
        NORM
      
      ⁡
      
        [
        
          b
          ,
          t
        
        ]
      
    
    =
    
      
        
          N
          ⁡
          
            [
            
              b
              ,
              t
            
            ]
          
        
        
          
            max
            b
          
          ⁢
          
            {
            
              N
              ⁡
              
                [
                
                  b
                  ,
                  t
                
                ]
              
            
            }
          
        
      
      .

11. A method according to claim 9 wherein the difference in spectral content between successive time blocks of the audio signal is calculated according to



  
    
      D
      ⁡
      
        [
        t
        ]
      
    
    =
    
      
        ∑
        b
      
      ⁢
      
        
        
          
            
              N
              NORM
            
            ⁡
            
              [
              
                b
                ,
                t
              
              ]
            
          
          -
          
            
              N
              NORM
            
            ⁡
            
              [
              
                b
                ,
                
                  t
                  -
                  1
                
              
              ]
            
          
        
        
      
    
  




  where




  
    
      
        N
        NORM
      
      ⁡
      
        [
        
          b
          ,
          t
        
        ]
      
    
    =
    
      
        
          N
          ⁡
          
            [
            
              b
              ,
              t
            
            ]
          
        
        
          
            avg
            b
          
          ⁢
          
            {
            
              N
              ⁡
              
                [
                
                  b
                  ,
                  t
                
                ]
              
            
            }
          
        
      
      .

12. A non-transitory computer-readable storage medium encoded with a computer program for causing a computer to perform the method of claim 1.

13. A method for modifying a parameter of an audio dynamics processor, comprising
detecting changes in a signal characteristic with respect to time in an audio signal,
identifying as auditory event boundaries changes greater than a threshold in said signal characteristic with respect to time in said audio signal, wherein an audio segment between consecutive boundaries constitutes an auditory event,
generating a parameter-modifying dynamically-varying control signal based on said auditory event, and
modifying the parameter of the audio dynamics processor as a function of the control signal.

14. An audio processing method in which a processor receives an input channel and generates an output channel that is generated by applying dynamic gain modifications to the input channel, comprising
detecting changes in signal characteristics with respect to time in the audio input channel,
identifying as auditory event boundaries changes in signal characteristics with respect to time in said input channel, wherein an audio segment between consecutive boundaries constitutes an auditory event in the channel, and
generating all or some of one or more dynamically-varying parameters of the audio dynamic gain modification method in response to auditory events and/or the degree of change in signal characteristics associated with said auditory event boundaries.

15. Apparatus for modifying a parameter of an audio dynamics processor, comprising
means for detecting changes in spectral characteristics with respect to time in an audio signal,
means for identifying as auditory event boundaries changes greater than a threshold in spectral characteristics with respect to time in said audio signal, wherein an audio segment between consecutive boundaries constitutes an auditory event,
means for generating a parameter-modifying dynamically-varying control signal based on said identified event boundaries, and
means for modifying the parameter of the audio dynamics processor as a function of the control signal.