WO2001041320A1 - Oversampled filter bank for subband processing - Google Patents

Oversampled filter bank for subband processing Download PDF

Info

Publication number
WO2001041320A1
WO2001041320A1 PCT/US1999/028343 US9928343W WO0141320A1 WO 2001041320 A1 WO2001041320 A1 WO 2001041320A1 US 9928343 W US9928343 W US 9928343W WO 0141320 A1 WO0141320 A1 WO 0141320A1
Authority
WO
WIPO (PCT)
Prior art keywords
subbands
frequency
filter bank
filter
path
Prior art date
Application number
PCT/US1999/028343
Other languages
French (fr)
Inventor
Douglas Wong
Original Assignee
Verance Corporation
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Verance Corporation filed Critical Verance Corporation
Priority to US10/148,161 priority Critical patent/US6947509B1/en
Priority to PCT/US1999/028343 priority patent/WO2001041320A1/en
Priority to AU18372/00A priority patent/AU1837200A/en
Publication of WO2001041320A1 publication Critical patent/WO2001041320A1/en

Links

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0248Filters characterised by a particular frequency response or filtering method
    • H03H17/0264Filter sets with mutual related characteristics
    • H03H17/0266Filter banks

Definitions

  • the present invention relates to a filter bank for use in digital signal processing, and in particular, to an efficient filter bank structure that uses oversampling and separate odd and even subband processing paths to reduce aliasing between subbands.
  • FIG. 1 illustrates a typical conventional filter bank structure.
  • a circuit 100 processes an input digital signal x(n) in N separate paths.
  • An analysis filter bank 120 includes respective analysis filters that decompose and transform the input signal x(n) into its frequency-domain subband components xb(0) , ...,xb(N-l) according to respective transfer functions H(0), H(l)...H(i),..., H(N-l) .
  • the subband components can be processed by distinct subband processors SB(0),
  • SB(1) ... SB(i) , ... , SB(N-l) of a subband processor 140 to be transformed into respective components yb(0), yb(l)..., yb(i ),..., yb(N-l).
  • a synthesis filter bank 160 reassembles and transforms the processed components into an output signal y(n) using synthesis filters F(0), F(l), ..., F(i) , ... F(N-l) .
  • Many text books provide a good introduction to the theory of digital filter banks.
  • DCT Discrete Cosine Transform
  • TDAC Time Domain Aliasing Cancellation
  • Cosine- modulated filter banks have been standardized into the MPEG-2 audio compression algorithms and Dolby Lab's(tm) AC3 compression algorithm.
  • a gain factor applied to one of the subband processors SB(i) will cause aliasing in the neighboring subbands, SB (i-l) and SB(i+l), during synthesis.
  • a higher order filter can be used. However, this decreases the temporal resolution of the filter banks as well as the computational efficiency since additional calculations are required.
  • Another approach is to try to avoid the -aliasing in the first place, by oversampling the subband components of interest. However, aliasing is still present between the subbands unless the oversampling ratio approaches the number of subbands M. This approach also is not computationally efficient.
  • a filter bank structure for subband processing that avoids constraining the type of subband processing, is computationally efficient, avoids aliasing between subbands, can be implemented using fast filter banks, and which has a relatively low filter order, and good temporal response and stop band rejection.
  • the filter bank structure should be suitable for use with any type of digital input signal, including one-dimensional signals such as audio signals, and two- dimensional signals such as video signals, and should allow any type of subband processing, including embedding of auxiliary data, noise shaping, and equalizing.
  • the filter bank structure should allow the use of real or complex filter banks.
  • the present invention provides a filter structure having the above and other advantages .
  • the present invention relates to an oversampled filter bank structure for subband processing that can be implemented using popular fast filter banks, such as FFT and cosine-modulated filter banks.
  • FFT fast filter banks
  • cosine-modulated filter banks are popular fast filter banks, such as FFT and cosine-modulated filter banks.
  • aliasing between subbands is substantially reduced, thus allowing independent processing of the subband signals.
  • the filter design criteria is significantly relaxed, thus resulting in a lower prototype filter order, better temporal response and higher stop band rejection.
  • a digital filter apparatus processes an digital input signal to provide a corresponding output signal in which subband processing has been performed.
  • a first frequency shifter frequency-shifts a first plurality of subbands of the digital input signal (such as the even subbands SBO, SB2, ...) to provide frequency- shifted subbands.
  • An oversampled analysis filter bank means filters the frequency-shifted subbands, along with a second plurality of subbands of the digital input signal (such as the odd subbands SB1, SB3 , ...), according to an oversampling ratio that is selected to reduce aliasing between adjacent ones of the frequency-shifted subbands, and between adjacent ones of the second plurality of subbands.
  • the oversampling ratio may be 2:1.
  • the oversampled analysis filter bank means may be implemented as a comb filter bank.
  • the first frequency shifter provides a ⁇ /M radian frequency shift, where M is the total number of subbands output from the analysis filter bank means.
  • M is the total number of subbands output from the analysis filter bank means.
  • the oversampled analysis filter means may be implemented as separate real filters for processing the even and odd subbands separately in parallel paths, or as a single complex filter bank which processes the even and odd subbands together.
  • the filter banks may be M/2-band filter banks.
  • a delay is provided for delaying the second plurality of subbands to compensate for a delay of the first frequency shifter.
  • the frequency shifter may be implemented using a critically-sampled cosine-modulated filter bank, or a Hubert transformer, for example.
  • a subband processor may be provided for processing the frequency- shifted subbands and the second plurality of subbands after filtering thereof.
  • the subband processing may involve any known technique, such as embedding of auxiliary data, noise shaping and so forth.
  • an oversampled synthesis filter means is provided for filtering the frequency- shifted subbands (e.g., even subbands) and the second plurality of subbands (e.g., odd subbands) after processing thereof at the subband processor.
  • a second frequency shifter shifts the already frequency- shifted subbands after filtering thereof at the synthesis filter means in a shift direction opposite to a shift direction of the first frequency shifter. Moreover, a second delay delays the second plurality of subbands to compensate for a delay of the second frequency shifter. Subsequently, a combiner combines all subbands to provide a digital output signal corresponding to the digital input signal.
  • the synthesis filter bank means may include separate real oversampled synthesis filter banks for the even and odd subbands, or a single complex synthesis filter bank.
  • the frequency shifters use cosine-modulated analysis and synthesis filter banks, proper frequency axis reversal is applied to each subband to maintain the orientation of the subband data.
  • a frequency rotation scheme can be used where this data is fed to the lowest subband of the synthesis filter of the frequency shifter.
  • a digital filter apparatus provides 4x oversampling to further reduce the likelihood of aliasing.
  • An input signal is processed in four parallel paths, three of which include frequency shifters for providing shifts of ⁇ /M, ⁇ /2M, and - ⁇ /2V radians, and one of which includes a delay for compensating for a processing time of the frequency-shifters .
  • the four sets of subbands are processed at an oversampled analysis filter bank means with an oversampling ratio (e.g., 4:1) that is selected to reduce aliasing between adjacent subbands in each set
  • filtering of two-dimensional data is provided.
  • FIG. 1 illustrates a typical conventional filter bank structure.
  • FIG. 2(a) illustrates a generalized filter structure with 2x-oversampled filter banks in odd and even subband paths in accordance with the present invention.
  • FIG. 2(b) illustrates a filter structure with 2x- oversampled real filter banks in odd and even subband paths in accordance with the present invention.
  • FIG. 2(c) illustrates a filter structure with a 2x-oversampled complex filter bank in accordance with the present invention.
  • FIG. 3(a) illustrates the spectrum for even subbands in accordance with the present invention.
  • FIG. 3(b) illustrates the spectrum for odd subbands in accordance with the present invention.
  • FIG. 4(a) illustrates a critically-sampled cosine- modulated frequency upshifter in accordance with the present invention.
  • FIG. 4(b) illustrates frequency axis reversal for a frequency shifter in accordance with the present invention.
  • FIG. 5 illustrates a Hubert Transform frequency shifter in accordance with the present invention.
  • FIG. 6 illustrates a critically-sampled cosine- modulated frequency downshifter in accordance with the present invention.
  • FIG. 7 illustrates a prototype filter design for analysis and synthesis filter banks in accordance with the present invention.
  • FIG. 8 illustrates a comb filter bank structure in accordance with the present invention.
  • FIG. 9 illustrates a graphical realization of one- dimensional filter bank processing in accordance with the present invention.
  • FIG. 10 illustrates a filter structure for one- dimensional filter bank processing with 4x-oversampled filter banks in four subband paths in accordance with the present invention.
  • FIG. 11(a) illustrates the subband structure for Path 0 of the filter of FIG. 10 in accordance with the present invention.
  • FIG. 11(b) illustrates the subband structure for Path 1 of the filter of FIG. 10 in accordance with the present invention.
  • FIG. 11(c) illustrates the subband structure for Path 2 of the filter of FIG. 10 in accordance with the present invention.
  • FIG. 11(d) illustrates the subband structure for Path 3 of the filter of FIG. 10 in accordance with the present invention.
  • FIG. 12(a) illustrates an original 2-D signal
  • FIG. 12(b) provides a graphical depiction of a conventional 2-D, critically-sampled filter bank.
  • FIG. 13(a) provides a graphical depiction of a path with no frequency shift (Path 0) in a 2-D, over- sampled filter bank in accordance with the present invention.
  • FIG. 13(b) provides a graphical depiction of a path with a horizontal frequency shift (Path 1) in a 2- D, over-sampled filter bank in accordance with the present invention.
  • FIG. 13(c) provides a graphical depiction of a path with a vertical frequency shift (Path 2) in a 2-D, over-sampled filter bank in accordance with the present invention.
  • FIG. 13(d) provides a graphical depiction of a path with horizontal and vertical frequency shifts (Path-3) in a 2-D, over-sampled filter bank in accordance with the present invention.
  • FIG. 13(e) provides a graphical depiction of the frequency regions that are processed in FIGs 13 (a) - 13(d) in a 2-D, over-sampled filter bank in accordance with the present invention.
  • FIG. 14 illustrates a real filter structure for two-dimensional filter bank processing with 2x- oversampled real filter banks in four subband paths in accordance with the present invention.
  • FIG. 2(a) illustrates a generalized filter structure with 2x-oversampled filter banks in odd and even subband paths in accordance with the present invention.
  • the filter apparatus 165 processes the input signal x(n) in an even subband path 170 and an odd subband path 180.
  • the even path 170 includes an M/2- band even subband oversampled analysis filter bank 172, and an M/2 -band even subband oversampled synthesis filter bank 174.
  • the odd path 180 includes an M/2 -band odd subband oversampled analysis filter bank 182, and an M/2 -band odd subband oversampled synthesis filter bank 184.
  • a subband processor can be provided between the filters 172 and 174, and between the filters 182 and 184.
  • a combiner 175 receives the output of the filter 174, x odd , and the output of the filter 184, x ⁇ v ⁇ n , to provide the output signal y(n) .
  • FIG. 2(b) illustrates a filter structure with 2x- oversampled real filter banks in odd and even subband paths in accordance with the present invention.
  • the filter apparatus 200 processes the input signal x(n) in an even subband path 210 and an odd subband path 250.
  • the even path 210 includes a ⁇ /M frequency shifter 212 and an M/2-band oversampled analysis filter bank 214.
  • the odd path 250 includes a delay 252, which compensates for the processing time of the frequency shifter 212, and an M/2 -band oversampled analysis filter bank 214', which is the complement of the filter bank 214.
  • a subband processor 240 receives the respective even and odd subband data, s ⁇ ven and s odd , respectively, of the input signal x(n) from the filters 214 and 214' , respectively.
  • the subband processing can include providing auxiliary data in each subband, harmonic analysis/synthesis as used in pitch-preserving time compression/expansion of audio, audio compression, and OFDM and transmultiplexers (transcoding multiplexers) as used in digital communications.
  • advantageous subband processing techniques are disclosed in commonly-assigned U.S. Patents 5,822,360 and 5,937,000.
  • the E-DNA ® data processing system of Solana Technology Development Corporation, San Diego, California, USA, is an example of such technology.
  • the even and odd subbands may be processed using separate and/or common circuitry at the subband processor 240.
  • the data output from the subband processor 240 in the even path 210, s ⁇ V ⁇ n ' is provided to an M/2-band synthesis filter bank 216 that is properly matched to the analysis filter bank 214, and a - ⁇ /M frequency shifter 218.
  • the data output from the subband processor 240 in the odd path 250, s odd ' is provided to an M/2 -band synthesis filter bank 216' that is properly matched to the analysis filter bank 214', and a delay 258, which compensates for the processing time of the frequency shifter 218.
  • a combiner 260 receives x ⁇ V ⁇ 1] and x odd to provide the final output signal y(n) .
  • frequency shifter 212 is shown as an upshifter and the frequency shifter 218 is shown as a downshifter, the opposite may be used.
  • M/2 is the downsampling ratio at the filters 214 and 214', and the upsampling ratio at the filters 216 and 216' .
  • M is the total number of subbands output from the analysis filter banks 214 and 214' together when 2x oversampling is used.
  • FIG. 2(c) illustrates a filter structure with a 2x-oversampled complex filter bank in accordance with the present invention. Like-numbered elements correspond to one another in the figures. Additionally, the prime symbol (ex: 200') denotes related elements.
  • a combiner 245 combines the even subband data x' with the delayed, odd subband data x d to form the complex value x'+jx d .
  • a corresponding synthesis filter 286 is used to recover the even subband data x ⁇ v ⁇ personally' and odd subband data x odd ' .
  • x ⁇ V ⁇ n ' and x odd ' are full-band data signals
  • s ⁇ ven ' and s odd ' are subband data signals.
  • the subband processing is the same in both processors 240' and 240.
  • FIG. 3 (a) illustrates the subband spectrum for even subbands in accordance with the present invention.
  • a horizontal axis 310' denotes frequency, ⁇
  • a vertical axis 320' denotes the magnitude of the spectrum X ⁇ V ⁇ n ( ⁇ ) of the even subband values just prior to recombining with the odd subband values to form the output signal y(n) .
  • the end-to-end filter response or transfer function is X ⁇ ven ( ⁇ ) /X ( ⁇ ) using the notation of FIGs 2(a)-2(c) .
  • SBO, SB2 and SB4 denote the zeroeth, second and fourth subbands respectively.
  • SBO extends between 0 and ⁇ / M ' while its complex conjugate, SBO*, extends from - ⁇ r,/M to 0.
  • SB2 extends between ⁇ /M and 3 7I: /M, and SB3 extends between 3 ⁇ /M and 5 ⁇ /M.
  • SB2* extends between - /M and -3- ⁇ /M, although part of it is not shown .
  • the actual passbands of the even subbands are shown at 337, 338, 340 and 342, while the corresponding idealized passbands, having a width of ⁇ /M, are shown at 337', 338', 340' and 342', respectively.
  • FIG. 3(b) illustrates the subband spectrum for odd subbands in accordance with the present invention.
  • a horizontal axis 310 denotes frequency, ⁇
  • a vertical axis 320 denotes the magnitude of the spectrum X odd ( ⁇ ) of the odd subband values just prior to recombining with the even subband values to form the output signal y(n) .
  • the end-to-end response or transfer function is X odd ( ⁇ ) / ⁇ ( ⁇ ) using the notation of FIGs 2 (a) - (c) .
  • SBl, SB3 and SB5 denote the first, third and fifth subbands, respectively.
  • SBl extends between 0 and 2 ⁇ /M
  • SB3 extends between 2 ⁇ /M and 4 ⁇ /M
  • SB5 extends between 4 ⁇ /M and 6 ⁇ /M.
  • the actual passbands of the odd subbands are shown at 330, 332, 334 and 336, while the corresponding idealized passbands, having a width of ⁇ /M and a sharp cut-off response, are shown at 330', 332', 334', and 336', respectively.
  • the pass bands 330, 332, 334, 336, 338, 340, 342,. . . of the oversampled subbands SBl*, SBl, SB3 , SB5 , SBO, SB2 and SB4 , . . ., respectively occupy less than the available bandwidth, 2 ⁇ /M.
  • a critically-sampled signal utilizes the entire subband bandwidth.
  • each subband has very little signal energy around its band edges 331, 333,
  • the oversampling ratio is determined by the ratio of the available bandwidth to the utilized passband, and is usually chosen to be two to minimize computation and implementation complexity, although other ratios can be used.
  • each of the subband paths 210, 250 has only M/2 subbands, but when we combine the two paths, there are a total of M subbands (output from the filter bank 284) , which cover the full bandwidth of the signal space, and thus the filters 200 or 200' overall each comprise an M-band oversampled filter bank.
  • the analysis filter banks 214, 214' are complements of one other, as are the synthesis filter banks 216, 216'.
  • the following design considerations should be taken into account. Namely, the difference in band alignment in the even and odd subband paths implies that two different filter banks need to be realized. Moreover, the even subband path is not easily realizable conventionally, since the subband centers are at multiples of ⁇ /M, which is where the downsampling/aliasing occurs with most fast filter banks algorithm.
  • the 2x oversampled filter banks of FIGs 2(b) and 2(c) solve the aforementioned problems.
  • the M/2 samples of x(n) are replicated in two paths, with one path (even path 210) going to the frequency shifter 212 of ⁇ /M radians, and the other path (odd path 250) going through the delay line 252, which compensates for the delay introduced by the frequency shifter 212.
  • the subbands samples are reassembled/synthesized into the odd and even time- domain sequences with M/2 samples using two sets of identical oversampled M/2 -band synthesis filter banks 216, 216' .
  • the synthesized even sequence is frequency shifted back by - ⁇ /M radians at the shifter 218.
  • the synthesized odd sequence is delayed at the delay 258 to align with the even path 210.
  • the odd and even sequences are added together at the combiner 260 to form M/2 time domain samples, y(n) .
  • x(n) consists of a sequence of M/2 samples, for which there will be a corresponding M/2 samples of y(n), which is the reconstructed version of x(n) .
  • Oversampling occurs since M/2 samples of x(n) are transformed into M samples of s ⁇ v ⁇ n and s odd combined in the subband domain. Thus, the transformation is considered to be over-complete since it uses more samples than is fundamentally necessary to represent the original sequence .
  • the frequency shifters 212, 218 can be realized in a variety of ways, as discussed below.
  • FIG. 4(a) illustrates a critically-sampled cosine- modulated frequency upshifter in accordance with the present invention.
  • This approach to frequency shifting uses an M-band (or multiple thereof) cosine-modulated filter bank to perform the frequency shifting.
  • a time- domain transfer function h(n) 410 decomposes the time sequence x(n) into M subband samples, e.g., in paths 420, 430, 440,... , 450.
  • SBO data is provided in path 420
  • SB2 data is provided in path 430, and so forth.
  • the subband samples are downsampled by a factor of M at a corresponding downsampler 422, 432, 442,..., 452, and cosine-modulated at a corresponding modulator 424, 434, 444, ...454.
  • a frequency rotation scheme can be employed where the highest subband output from the analysis filter 405, in path 450 is fed to the modulator 426 of the lowest subband of the synthesis filter 406, in path 420. All information in x(n) can therefore be retained.
  • Frequency axis reversal functions 461-463 are used due to the structure of the cosine-modulated filter banks, as discussed further in connection with FIG. 4(b) .
  • the samples of SBO, in path 420, are provided to modulator 436, in path 430.
  • the samples of SB2 , in path 430, are provided to modulator 446, in path 440, and so forth.
  • the samples in each subband in the analysis filter bank 405 are fed to the next higher subband in the synthesis filter bank 406.
  • Cosine modulators 426, 436, 446, ... , 448 demodulate the subband samples from the next lower subband of the analysis filter 405 and provide them to respective upsamplers 428, 438, 448, ... 458, for upsampling by a factor of M.
  • the resulting data is provided to a function f (n) 460, which re-assembles the samples in proper sequence to obtain the time sequence x' (n) , which is shifted by ** ⁇ /M radians relative to x(n) .
  • This method is fairly efficient, but the delay is limited by the delay of the M-band filter banks.
  • a second approach for frequency shifting uses an M-band (or multiple thereof, e.g., 2M, 3M, ...) Discrete Fourier Transform (DFT) type of filter bank such as used in Time Domain Aliasing Cancellation
  • TDAC time domain aliasing cancellation
  • FIG. 4(b) illustrates frequency axis reversal for a frequency shifter in accordance with the present invention.
  • Frequency axis reversal is an inherent characteristic of cosine-modulated filter banks.
  • the spectrum of X( ⁇ ) is shown at 470.
  • each (non-discarded) subband undergoes frequency axis reversal, as shown at 475.
  • a corresponding reversal must be implemented, as shown at 480.
  • the corrected subbands can then be input to the synthesis filter 406 of the frequency shifter
  • a first set of every other subband output from the analysis filter bank 405 contains signals that are frequency axis flipped, i.e., 0 ⁇ and ⁇ _- > 0.
  • a second set of the remaining alternate subbands will be frequency axis flipped after processing at the synthesis filter 406.
  • Proper frequency axis reversal for each subband avoids this problem.
  • FIG. 5 illustrates a Hilbert Transform frequency shifter in accordance with the present invention.
  • a third approach for implementing the frequency shifters 212 or 218 of FIGs 2(b) and 2(c) uses a Hilbert transformer with a complex modulator.
  • the frequency shifter 212'' includes a Hilbert Transformer 510, a complex modulator 520, and a function 530 for obtain the real portion of the output of the modulator 520.
  • An ideal Hilbert transformer is an all-pass filter than imparts a 90° phase shift on the input signal.
  • the frequency response is specified as:
  • the Hilbert transformer can be implemented using a Fast Fourier Transform (FFT) or a Finite Impulse Response (FIR) filter, for example.
  • FFT Fast Fourier Transform
  • FIR Finite Impulse Response
  • the Hilbert transform of the sequence is then modulated with a complex modulator 520 with center frequency at 71 / , to achieve a frequency upshift.
  • the real part of the modulated sequence represents the frequency- shifted version of the original time-domain sequence.
  • the frequency shifter 212 ' ' can be used to achieve a frequency downshift by using exp(-j ⁇ /M) instead of exp(j ⁇ /M) at the modulator 520.
  • the input signal to the frequency shifter is then x ⁇ v ⁇ Budapest ' (n) , and the output signal is x ⁇ V ⁇ n (n) .
  • FIG. 6 illustrates a critically-sampled cosine- modulated frequency downshifter in accordance with the present invention.
  • the frequency downshifter 218' which includes subband paths 620, 630, 640,..., 650, is analogous to the frequency upshifter 212'.
  • An analysis filter bank 605 and synthesis filter bank 606 correspond to the analysis filter bank 405 and synthesis filter bank 406, respectively, of FIG. 4(a), except each modulator 426, 436, 446, ..., 448 receives subband samples from the next higher (instead of lower) subband from the analysis filter bank 605.
  • Frequency axis reversal functions 661-663 are also provided. They are analogous to the reversal functions 461-463, as discussed in connection with FIGs 4(a) and 4(b) .
  • FIG. 7 illustrates a prototype filter design for analysis and synthesis filter banks in accordance with the present invention.
  • the oversampled M/2-band filter banks 214, 214', 216 and 216', along with the complex filters 284 and 286, may be implemented as comb filters.
  • a comb filter may be viewed as a filter in which nulls occur periodically across the frequency band.
  • the comb filter structure of the present invention is realized using a modified form of cosine- modulated filter bank, which can be implemented with fast transforms such as the DCT or FFT.
  • Other uniform filter banks may be used, such as a pseudo-Quadrature Mirror Filter (QMF), as described in J . H. Rothweiler, Polyphase Quadrature Fil ters - A New Subband Coding
  • TDAC Transimeter-to-Digital Converter
  • the prototype filter is designed to have a bandwidth of ⁇ .
  • the prototype filter is designed to have a bandwidth of ⁇ /2.
  • the response 730 of the corresponding comb filter is shown in FIG 7.
  • a horizontal axis 710 indicates frequency, while a vertical axis 720 indicates the magnitude of the response, in decibels (dB) .
  • dB decibels
  • a series of non-overlapping band-pass filters result which have band edges at multiples of ⁇ /4 and 3 ⁇ /4, respectively.
  • the effective bandwidth of the pass band 750 is therefore ⁇ /2. Note that an attenuation of -40 dB can be considered to represent an amplitude of essentially zero.
  • the stop band rejection of the filter is designed to be low enough so that aliasing is small or negligible at 0 and ⁇ .
  • Stop bands 740 and 760 are provided between 0 and ⁇ /4 , and between 3 ⁇ /4 and ⁇ » respectively.
  • An additional benefit when designing this prototype filter is that we are no longer constrained by the severe aliasing cancellation matrix, discussed in P.P. Vaidyanathan, Mul tirate System And Fil ter Banks . Prentice Hall 1993, as with the critically- sampled case, thus yielding additional freedom in the design. Effectively, the filter order can be reduced substantially to provide better temporal resolution, and to achieve improved stop band rejection characteristics .
  • the response 730 applies to both the analysis and synthesis filter banks since, in the design of a cosine-modulated filter bank, one only needs to specify the prototype filter characteristics, and the response of the analysis and synthesis filters will be dictated by the prototype filter.
  • FIG. 8 illustrates a comb filter bank structure in accordance with the present invention.
  • the circuit 800 can perform the functions of the filter 214, subband processor 240, and filter 216 of FIG. 2(b), or the filter 214', subband processor 240, and filter 216' of FIG. 2(b).
  • An analysis comb filter bank 214' ' includes a transfer function h(n) 810 for providing the input signal in M/2 subbands, downsamplers 820, 822, ...,
  • FIG. 9 illustrates a graphical realization of one- dimensional filter bank processing in accordance with the present invention.
  • the graphical illustration shows how the oversampled filter banks are realized with the frequency shifter and comb filter banks already discussed.
  • X d ( ⁇ ) (920) is formed by delaying X( ⁇ ) .
  • the passbands of the oversampled filter banks are shown at 922.
  • X odd (co) (926) is formed, whose spectrum is essentially the same as S odd ( ⁇ ) (924) .
  • X' ( ⁇ ) (930) is formed by frequency- shifting X( ⁇ ) by ⁇ /M.
  • the passbands of the oversampled filter banks are shown at 932.
  • X ⁇ v ⁇ n ( ⁇ ) (936) is formed.
  • FIG. 10 illustrates a 2M-band filter (e.g., four, M/2-band filters) with 4x-oversampled filter banks with four subband paths in accordance with the present invention.
  • 2M-band filter e.g., four, M/2-band filters
  • the filter structure 1000 includes four subband processing paths, i.e., Path 0 (1010), Path 1 (1020), Path 2 (1030) , and Path 3 (1040) .
  • Path 0 (1010) includes a ⁇ /M frequency shifter 1012, an M/2 -band analysis filter bank 1014, an M/2 -band synthesis filter
  • Path 1 includes a ⁇ /2M frequency shifter 1012', an M/2-band analysis filter bank 1014', an M/2-band synthesis filter 1016', and a - ⁇ /2M frequency shifter 1018'. Note that the frequency shift functions described herein can be implemented as discussed previously (e.g., in connection with FIGs 4(a), 4(b), 5 and 6) with modification as required for the magnitude and direction of frequency shift.
  • Path 2 (1030) includes a delay 252, an M/2-band analysis filter bank 1014'', an M/2 -band synthesis filter 1016' ' , and a delay 252' .
  • the delay 252 compensates for the processing time of the frequency shifters 1012, 1012', and 1012'', while the delay 252' compensates for the processing time of the frequency shifters 1018, 1018', and 1018''.
  • Path 3 (1040) includes a - ⁇ /2M frequency shifter 1012' ' , an M/2-band analysis filter bank 1014' ' , an M/2 -band synthesis filter 1016'', and a ⁇ / 2M frequency shifter 1018 ' ' .
  • each successive subband is processed in a successive one of the Paths 1010, 1020, 1030, or 1040.
  • Paths 1010, 1020, 1030, or 1040 For example, assuming 32 subbands are used, i.e., numbered SB0-SB31, SBO, SB4 , SB8, ... are processed in Path 0 (1010), SBl, SB5 , SB9 , ... are processed in Path 1 (1020), SB2 , SB6, SB10, ... are processed in Path 2 (1030), and SB3 , SB7, SB11, ... are processed in Path 3 (1040) .
  • the filter 1000 could optionally be modified to use one or two complex filters, instead of four real filters, based on the discussion of FIG. 2(c) .
  • FIG. 11(a) illustrates the subband structure for Path 0 of the filter bank of FIG. 10 in accordance with the present invention.
  • the pass bands 1112, 1116, 1120, ... of the 4x oversampled subbands SBO*, SBO, SB4, SB8, ... occupy only about one-fourth of the available bandwidth, 2 ⁇ /M.
  • the corresponding idealized pass bands are shown at 1112', 1116', 1120', ... , respectively.
  • X Path - o ( ⁇ ) designates the signal in Path 0 (1010), as shown in FIG. 10.
  • FIG. 11(b) illustrates the subband structure for Path 1 of the filter bank of FIG. 10 in accordance with the present invention.
  • pass bands 1132, 1136, 1140, 1144, ... of the 4x oversampled subbands SBl*, SBl, SB5 , SB9 , ... are provided, with band edges 1130, 1134, 1138, ...
  • the corresponding idealized pass bands 1132', 1136', 1140', 1144', ... are also shown.
  • ⁇ p th - ⁇ ( ⁇ ) designates the signal in Path 1 (1020), as shown in FIG. 10.
  • FIG. 11(c) illustrates the subband structure for Path 2 of the filter bank of FIG. 10 in accordance with the present invention.
  • pass bands 1152, 1156, 1160, ... of the 4x oversampled subbands SB2 , SB6, SB10, ... are provided, with band edges 1150, 1154, 1158, ...
  • the corresponding idealized pass bands 1152', 1156', 1160', ... are also shown.
  • X Pth - 2 (oc>) designates the signal in Path 2 (1030), as shown in FIG. 10.
  • FIG. 11(d) illustrates the subband structure for Path 3 of the filter bank of FIG. 10 in accordance with the present invention.
  • pass bands 1182, 1186, ... of the 4x oversampled subbands SB2 , SB6, SB10, ... are provided, with band edges 1180, 1184, 1188, ...
  • the corresponding idealized pass bands 1182' , 1186' , ... are also shown.
  • ⁇ p ah - 3 ( ⁇ ) designates the signal in Path 3 (1040), as shown in FIG. 10.
  • FIG. 12(a) illustrates an original 2-D signal (image) region with several subregions for processing in accordance with the present invention.
  • the image 1200 is shown as including a number of subregions A-P.
  • image data is processed by subdividing it into several blocks, each of which includes many pixels.
  • a block-based spatial transform such as the Discrete Cosine Transform (DCT) is often used, to provide frequency domain coefficients.
  • DCT Discrete Cosine Transform
  • Other possible transforms include the Discrete Fourier
  • FIG. 12(b) provides a graphical depiction of a conventional 2-D, critically-sampled filter bank.
  • FIGs 12(b) -13(e) are diagrammatic illustrations that show a composite of image data in the frequency domain (such as DCT frequency coefficients or the like) and associated filters.
  • Each filter in the image/filter composite 1200' is labeled with (x,y) coordinates, and is overlaid on the corresponding image data subregion which it processes.
  • x and y are used to denote horizontal and vertical direction, respectively, and should not be confused with the input or output signals x(n) and y(n) . The meaning should be clear from the context.
  • An 8x8 2-D filter bank is assumed for illustration purposes, although different sizes can be used.
  • an image subregion 1210 has data that is processed by a filter labeled (7,7) .
  • the image data within the subregion 1210, which is part of the larger image subregion "P" is passed through a pass band 1212 of the corresponding filter, while image data beyond the band edges 1214 will not be passed by this filter.
  • the 2-D, critically-sampled filter bank suffers from the same problem of aliasing between subbands as discussed for the 1-D critically sampled filter bank.
  • FIG. 13(a) provides a graphical depiction of a path with no frequency shift (Path 0) in a 2-D, over- sampled filter bank in accordance with the present invention .
  • a 2x oversampled 2-D filter bank structure in accordance with the present invention processes image or other 2-D data in four processing paths, namely Paths 0-3.
  • the 2-D data toward the center of each image region is processed in the first path, Path 0, while data at vertical boundaries between regions (x-axis shifted by ⁇ /M) is processed in a second path, Path-1, data at horizontal boundaries between regions (y-axis shifted by ⁇ / M ) 1S processed in a third path, Path-2, and data at corners of the boundaries between regions (x-axis and y-axis shifted by ⁇ /M) 1S processed in a fourth path, Path-3.
  • each filter is labeled with (P,x,y) coordinates in the image/filter composite 1300-0, where the "P" coordinate denotes the path, and is shown overlaid with the corresponding image data subregion.
  • an image subregion 1310, labeled "H” has data that is processed by an oversampled filter (0,1,3) .
  • the image data toward the center of the subregion 1310 will pass through a pass band 1312 of the corresponding filter essentially unchanged, while data between the pass band 1312 and band edges 1316 is stopped by a stop band 1314.
  • FIG. 13(b) provides a graphical depiction of a path with a horizontal frequency shift (Path 1) in a 2- D, over-sampled filter bank in accordance with the present invention.
  • the image/filter composite 1300-1 shows that data at vertical boundaries between regions (x-axis shifted by ⁇ /M) 1S processed using the corresponding filters as shown. For example, data between the regions "A” and "B” is processed using a filter with index (1,0,1).
  • FIG. 13(c) provides a graphical depiction of a path with a vertical frequency shift (Path 2) in a 2-D, over-sampled filter bank in accordance with the present invention.
  • the image/filter composite 1300-2 shows that data at horizontal boundaries between regions (y-axis shifted by ⁇ /M) is processed using the corresponding filters as shown. For example, data between the regions "A” and “E” is processed using a filter with index (2,1,0).
  • FIG. 13(d) provides a graphical depiction of a path with horizontal and vertical frequency shifts (Path-3) in a 2-D, over-sampled filter bank in accordance with the present invention.
  • the image/filter composite 1300-3 shows that data at corners of the boundaries between regions (x-axis and y-axis shifted by ⁇ /M) is processed in a fourth path.
  • data at the corner of regions "A”, “B”, “E” and “F” is processed using a filter with index (3,1,1).
  • FIG. 13(e) provides a graphical depiction of the freq ⁇ ency regions that are processed in FIGs 13 (a) - 13(d) in a 2-D, over-sampled filter bank in-' accordance with the present invention.
  • the image/filter composite 1300 shows the aggregation of the 2-D data regions that are processed by each of the four paths. In the example shown, sixteen regions are processed in each path, for a total of sixty-four regions.
  • a region 1385 illustrates an area of imperfect reconstruction by the filter banks.
  • the area 1385 is not covered by any of the four paths. To obtain complete coverage, additional processing paths can be introduced. This is similar to the 1-D case, where we decided to drop the last subband. Notice that the area 1385 corresponds to high- requency regions of the 2-D signal, so, for most applications, this data is not used at all and can be discarded without noticeably impairing the resulting output signal quality.
  • the present invention can be extended to an oversampled filter bank structure for use with N- dimensions .
  • an example of a 2-D oversampled filter bank has been shown.
  • Such a filter bank can be used, e.g., in applications relating to video or image processing.
  • FIG. 14 illustrates a real filter structure for two-dimensional filter bank processing with 2x- oversampled real filter banks in four subband paths in accordance with the present invention.
  • frequency shift functions described herein can be implemented as discussed previously (e.g., in connection with FIGs 4(a), 4(b), -5 and 6) with modification as required for the magnitude and direction of frequency shift.
  • the filter structure can be modified to used one or more complex filters as discussed in connection with FIG. 2(c) .
  • the filter structure 1400 includes processing paths 1410, 1420, 1430 and 1440. Generally, a number of delay elements (in path 0: 252-0, 252-0", 252-0'',
  • 252-0''*; in path 1: 252-1, 252-1'; and, in path 2: 252-2, 252-2') are provided to compensate for the processing time of the functions that are shown as vertically aligned in the figure.
  • Each path includes an analysis filter bank 1414-0,
  • Path 1 (1420) includes counterpart x-axis frequency shifters 1412-1 and 1418-1, while Path 2 (1430) includes counterpart y-axis frequency shifters
  • Path 3 (1440) includes both counterpart x-axis frequency shifters 1412-3 and 1418-
  • a combiner 260' receives the output in each path to provide the final output signal y(n) .
  • the x- axis and y-axis terminology should not be confused with the input signal x(n) and output signal y(n) .
  • x-axis (1412-1, 1412-3) and y-axis (1413-2, 1413-3) shift functions can use common hardware and/or software.
  • the present invention provides an oversampled filter bank structure that can be implemented using popular and efficient fast filter banks.
  • aliasing between subbands is substantially reduced, thus allowing the independent processing of the subband signals using any known subband processing technique.
  • an input audio or video signal may be processed to embed auxiliary data in its subbands.
  • the auxiliary data may be used for any number of purposes, such as copyright protection, broadcast verification, identification of soundtrack, carrying electronic coupon data, and so forth.
  • equalizing and noise shaping of the subbands as well as other processing may be provided.
  • the oversampled filter bank can be used in OFDM, which has found application in Digital Subscriber Loop (DSL) and digital broadcast systems.
  • DSL Digital Subscriber Loop
  • even subbands of an input signal are frequency- shifted prior to filtering and subband processing, while the odd subbands are delayed to compensate for the processing time of the frequency shifting.
  • Separate analysis and synthesis filter banks may be provided for the even and odd subbands, or the subbands may be processed together in a single complex filter bank. 2x oversampling may be advantageously used in the filter banks.
  • Frequency shifting may be achieved using critically-sampled cosine filter banks, or a Hilbert transform, for example.
  • the subbands are processed in four subband paths, and 4x oversampling is used at the filter banks.
  • the invention may also be extended to two-dimensional or higher filtering applications.
  • a 2-D oversampled filter bank structure that processes data in four paths, three of which undergo frequency shifts.
  • the structure is suited for processing image or other 2-D data, such as seismic or other vibrational data, temperature data and the like.
  • the digital input signal discussed herein may only be a portion of an overall signal that is provided at some upstream location. Only a portion of the overall signal need be processed in accordance with the invention. Additionally, the system of the invention may be aggregated to process additional input signals, or additional portions of an input signal.

Abstract

An oversampled filter bank structure that can be implemented using popular and efficient fast filter banks to allow subband processing of an input signal with substantially reduced aliasing between subbands. Even subbands (SB0, SB2, SB4,...) of an input signal (x(n)) are frequency-shifted (212, 1012, 1012', 1012'') prior to analysis filtering (214, 214', 214'') at a 2x oversampled filter bank, subband processing (240, 240', 240''), and synthesis filtering (216, 216', 216''). A subsequent frequency-shift (218, 218') returns the even subbands to their original band positions. The odd subbands (SB1, SB3, SB5,...) are delayed (252) to compensate for the processing time of the frequency shifting. Separate analysis (214, 214') and synthesis (216, 216') filter banks may be provided for the even and odd subbands, or common complex analysis (284) and synthesis (286) filter banks may be used. In another embodiment, the subbands are processed in four subband paths (Paths 0, 1, 2, 3), and 4x oversampling is used. A filter bank structure (1400) for 2-D data is also provided.

Description

OVERSA PLED FILTER BANK FOR SUBBAND PROCESSING
BACKGROUND OF THE INVENTION
The present invention relates to a filter bank for use in digital signal processing, and in particular, to an efficient filter bank structure that uses oversampling and separate odd and even subband processing paths to reduce aliasing between subbands.
The digital filter bank is an enabling technology for modern audio and video processing systems, and, more recently, digital data communication systems such as Orthogonal Frequency-Division Multiplexing (OFDM) . FIG. 1 illustrates a typical conventional filter bank structure. A circuit 100 processes an input digital signal x(n) in N separate paths. An analysis filter bank 120 includes respective analysis filters that decompose and transform the input signal x(n) into its frequency-domain subband components xb(0) , ...,xb(N-l) according to respective transfer functions H(0), H(l)...H(i),..., H(N-l) .* The subband components can be processed by distinct subband processors SB(0),
SB(1) ... SB(i) , ... , SB(N-l) of a subband processor 140 to be transformed into respective components yb(0), yb(l)..., yb(i ),..., yb(N-l). Various types of subband processing can be performed. A synthesis filter bank 160 reassembles and transforms the processed components into an output signal y(n) using synthesis filters F(0), F(l), ..., F(i) , ... F(N-l) . Many text books provide a good introduction to the theory of digital filter banks.
Of particular interest is the class of critically- sampled uniform filter banks, which have found wide application in the areas of audio and video 'processing. For example, the DCT (Discrete Cosine Transform) is used in the MPEG-2 video compression engine, whereas TDAC (Time Domain Aliasing Cancellation) and Cosine- modulated filter banks have been standardized into the MPEG-2 audio compression algorithms and Dolby Lab's(tm) AC3 compression algorithm.
While the efficiency of these filter banks makes them suitable for many signal processing applications, they suffer from aliasing between the subbands. Most practical filters have finite rejection at the Nyquist frequency so the signals beyond the Nyquist frequency are not sufficiently attenuated prior to downsampling, and appear as aliased components in the downsampled signal . Most of these filter banks are designed to be "aliasing canceling", which means that the synthesis filter bank is specifically designed to account for, and cancel, the aliasing caused by the analysis filter banks. However, this cancellation property severely constrains the kind of processing that can be introduced by the subband processors SB(0) , ..., SB(N- 1) . For example, a gain factor applied to one of the subband processors SB(i) will cause aliasing in the neighboring subbands, SB (i-l) and SB(i+l), during synthesis. This makes the use of such filter banks unsuitable for applications such as subband equalizers or noise shaper. To reduce aliasing, a higher order filter can be used. However, this decreases the temporal resolution of the filter banks as well as the computational efficiency since additional calculations are required. Another approach is to try to avoid the -aliasing in the first place, by oversampling the subband components of interest. However, aliasing is still present between the subbands unless the oversampling ratio approaches the number of subbands M. This approach also is not computationally efficient.
Accordingly, it would be desirable to provide a filter bank structure for subband processing that avoids constraining the type of subband processing, is computationally efficient, avoids aliasing between subbands, can be implemented using fast filter banks, and which has a relatively low filter order, and good temporal response and stop band rejection.
The filter bank structure should be suitable for use with any type of digital input signal, including one-dimensional signals such as audio signals, and two- dimensional signals such as video signals, and should allow any type of subband processing, including embedding of auxiliary data, noise shaping, and equalizing. The filter bank structure should allow the use of real or complex filter banks.
The present invention provides a filter structure having the above and other advantages . SUMMARY OF THE INVENTION
The present invention relates to an oversampled filter bank structure for subband processing that can be implemented using popular fast filter banks, such as FFT and cosine-modulated filter banks. Advantageously, aliasing between subbands is substantially reduced, thus allowing independent processing of the subband signals. Also, the filter design criteria is significantly relaxed, thus resulting in a lower prototype filter order, better temporal response and higher stop band rejection.
In a particular embodiment, a digital filter apparatus processes an digital input signal to provide a corresponding output signal in which subband processing has been performed. A first frequency shifter frequency-shifts a first plurality of subbands of the digital input signal (such as the even subbands SBO, SB2, ...) to provide frequency- shifted subbands. An oversampled analysis filter bank means filters the frequency-shifted subbands, along with a second plurality of subbands of the digital input signal (such as the odd subbands SB1, SB3 , ...), according to an oversampling ratio that is selected to reduce aliasing between adjacent ones of the frequency-shifted subbands, and between adjacent ones of the second plurality of subbands.
The oversampling ratio may be 2:1. The oversampled analysis filter bank means may be implemented as a comb filter bank. The first frequency shifter provides a π/M radian frequency shift, where M is the total number of subbands output from the analysis filter bank means. Preferably, M>8. The oversampled analysis filter means may be implemented as separate real filters for processing the even and odd subbands separately in parallel paths, or as a single complex filter bank which processes the even and odd subbands together. The filter banks may be M/2-band filter banks.
In any case, a delay is provided for delaying the second plurality of subbands to compensate for a delay of the first frequency shifter.
The frequency shifter may be implemented using a critically-sampled cosine-modulated filter bank, or a Hubert transformer, for example.
A subband processor may be provided for processing the frequency- shifted subbands and the second plurality of subbands after filtering thereof. The subband processing may involve any known technique, such as embedding of auxiliary data, noise shaping and so forth.
After subband processing, an oversampled synthesis filter means is provided for filtering the frequency- shifted subbands (e.g., even subbands) and the second plurality of subbands (e.g., odd subbands) after processing thereof at the subband processor.
A second frequency shifter shifts the already frequency- shifted subbands after filtering thereof at the synthesis filter means in a shift direction opposite to a shift direction of the first frequency shifter. Moreover, a second delay delays the second plurality of subbands to compensate for a delay of the second frequency shifter. Subsequently, a combiner combines all subbands to provide a digital output signal corresponding to the digital input signal.
Analogous to the analysis filter bank means, the synthesis filter bank means may include separate real oversampled synthesis filter banks for the even and odd subbands, or a single complex synthesis filter bank. Importantly, when the frequency shifters use cosine-modulated analysis and synthesis filter banks, proper frequency axis reversal is applied to each subband to maintain the orientation of the subband data.
Also, in the frequency shifter, to avoid discarding data from the highest even subband at the analysis filter of the frequency shifter, a frequency rotation scheme can be used where this data is fed to the lowest subband of the synthesis filter of the frequency shifter.
In an optional embodiment, a digital filter apparatus provides 4x oversampling to further reduce the likelihood of aliasing. An input signal is processed in four parallel paths, three of which include frequency shifters for providing shifts of π/M, π/2M, and -π/2V radians, and one of which includes a delay for compensating for a processing time of the frequency-shifters . The four sets of subbands are processed at an oversampled analysis filter bank means with an oversampling ratio (e.g., 4:1) that is selected to reduce aliasing between adjacent subbands in each set
In a further optional embodiment, filtering of two-dimensional data is provided.
Corresponding methods are also presented.
BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 illustrates a typical conventional filter bank structure.
FIG. 2(a) illustrates a generalized filter structure with 2x-oversampled filter banks in odd and even subband paths in accordance with the present invention.
FIG. 2(b) illustrates a filter structure with 2x- oversampled real filter banks in odd and even subband paths in accordance with the present invention.
FIG. 2(c) illustrates a filter structure with a 2x-oversampled complex filter bank in accordance with the present invention.
FIG. 3(a) illustrates the spectrum for even subbands in accordance with the present invention.
FIG. 3(b) illustrates the spectrum for odd subbands in accordance with the present invention.
FIG. 4(a) illustrates a critically-sampled cosine- modulated frequency upshifter in accordance with the present invention.
FIG. 4(b) illustrates frequency axis reversal for a frequency shifter in accordance with the present invention.
FIG. 5 illustrates a Hubert Transform frequency shifter in accordance with the present invention.
FIG. 6 illustrates a critically-sampled cosine- modulated frequency downshifter in accordance with the present invention. FIG. 7 illustrates a prototype filter design for analysis and synthesis filter banks in accordance with the present invention.
FIG. 8 illustrates a comb filter bank structure in accordance with the present invention.
FIG. 9 illustrates a graphical realization of one- dimensional filter bank processing in accordance with the present invention.
FIG. 10 illustrates a filter structure for one- dimensional filter bank processing with 4x-oversampled filter banks in four subband paths in accordance with the present invention.
FIG. 11(a) illustrates the subband structure for Path 0 of the filter of FIG. 10 in accordance with the present invention.
FIG. 11(b) illustrates the subband structure for Path 1 of the filter of FIG. 10 in accordance with the present invention.
FIG. 11(c) illustrates the subband structure for Path 2 of the filter of FIG. 10 in accordance with the present invention.
FIG. 11(d) illustrates the subband structure for Path 3 of the filter of FIG. 10 in accordance with the present invention. FIG. 12(a) illustrates an original 2-D signal
(image) region with several subregions for processing in accordance with the present invention.
FIG. 12(b) provides a graphical depiction of a conventional 2-D, critically-sampled filter bank. FIG. 13(a) provides a graphical depiction of a path with no frequency shift (Path 0) in a 2-D, over- sampled filter bank in accordance with the present invention.
FIG. 13(b) provides a graphical depiction of a path with a horizontal frequency shift (Path 1) in a 2- D, over-sampled filter bank in accordance with the present invention.
FIG. 13(c) provides a graphical depiction of a path with a vertical frequency shift (Path 2) in a 2-D, over-sampled filter bank in accordance with the present invention.
FIG. 13(d) provides a graphical depiction of a path with horizontal and vertical frequency shifts (Path-3) in a 2-D, over-sampled filter bank in accordance with the present invention. FIG. 13(e) provides a graphical depiction of the frequency regions that are processed in FIGs 13 (a) - 13(d) in a 2-D, over-sampled filter bank in accordance with the present invention.
FIG. 14 illustrates a real filter structure for two-dimensional filter bank processing with 2x- oversampled real filter banks in four subband paths in accordance with the present invention.
DETAILED DESCRIPTION OF THE INVENTION
The present invention relates to a filter bank for use in digital signal processing to prepare "a signal for subband processing. FIG. 2(a) illustrates a generalized filter structure with 2x-oversampled filter banks in odd and even subband paths in accordance with the present invention.
The filter apparatus 165 processes the input signal x(n) in an even subband path 170 and an odd subband path 180. The even path 170 includes an M/2- band even subband oversampled analysis filter bank 172, and an M/2 -band even subband oversampled synthesis filter bank 174. Similarly, the odd path 180 includes an M/2 -band odd subband oversampled analysis filter bank 182, and an M/2 -band odd subband oversampled synthesis filter bank 184.
A subband processor can be provided between the filters 172 and 174, and between the filters 182 and 184.
A combiner 175 receives the output of the filter 174, xodd, and the output of the filter 184, xβvβn, to provide the output signal y(n) .
FIG. 2(b) illustrates a filter structure with 2x- oversampled real filter banks in odd and even subband paths in accordance with the present invention.
Here, a more specific implementation of the generalized structure of FIG. 2(a) is provided. The filter apparatus 200 processes the input signal x(n) in an even subband path 210 and an odd subband path 250. The even path 210 includes a π/M frequency shifter 212 and an M/2-band oversampled analysis filter bank 214. The odd path 250 includes a delay 252, which compensates for the processing time of the frequency shifter 212, and an M/2 -band oversampled analysis filter bank 214', which is the complement of the filter bank 214. A subband processor 240 receives the respective even and odd subband data, sθven and sodd, respectively, of the input signal x(n) from the filters 214 and 214' , respectively. The subband processing can include providing auxiliary data in each subband, harmonic analysis/synthesis as used in pitch-preserving time compression/expansion of audio, audio compression, and OFDM and transmultiplexers (transcoding multiplexers) as used in digital communications. For example, advantageous subband processing techniques are disclosed in commonly-assigned U.S. Patents 5,822,360 and 5,937,000. The E-DNA® data processing system of Solana Technology Development Corporation, San Diego, California, USA, is an example of such technology.
The even and odd subbands may be processed using separate and/or common circuitry at the subband processor 240.
The data output from the subband processor 240 in the even path 210, sΘVθn' , is provided to an M/2-band synthesis filter bank 216 that is properly matched to the analysis filter bank 214, and a -π/M frequency shifter 218. The data output from the subband processor 240 in the odd path 250, sodd' , is provided to an M/2 -band synthesis filter bank 216' that is properly matched to the analysis filter bank 214', and a delay 258, which compensates for the processing time of the frequency shifter 218.
A combiner 260 receives xΘVΘ1] and xodd to provide the final output signal y(n) .
Note that the frequency shifting is provided in the even subband path to allow the use of fast filter banks. Without this frequency shifting, the downsampling edges would land on the center of the passband of interest.
Additionally, while the frequency shifter 212 is shown as an upshifter and the frequency shifter 218 is shown as a downshifter, the opposite may be used.
M/2 is the downsampling ratio at the filters 214 and 214', and the upsampling ratio at the filters 216 and 216' . M is the total number of subbands output from the analysis filter banks 214 and 214' together when 2x oversampling is used.
M=32 (M/2=16) has been used satisfactorily by the present inventors. M should preferably be selected as a multiple of two to allow the use of efficient algorithms, such as the DFT and FFT, in the filter implementation. For practical purposes, M=8 as a minimum. The efficiency of the fast filter banks does not become substantial unless M>>log2(M).
FIG. 2(c) illustrates a filter structure with a 2x-oversampled complex filter bank in accordance with the present invention. Like-numbered elements correspond to one another in the figures. Additionally, the prime symbol (ex: 200') denotes related elements.
Here, a single complex M/2 -band analysis filter bank 284 is used instead of having separate filter banks for the odd and even subbands . A combiner 245 combines the even subband data x' with the delayed, odd subband data xd to form the complex value x'+jxd.
After subband processing at the processor 240', a corresponding synthesis filter 286 is used to recover the even subband data xβvβ„' and odd subband data xodd' .
In Figs. 2(b) and 2(c), note that xΘVΘn' and xodd' are full-band data signals, whereas sβven' and sodd' are subband data signals. The subband processing is the same in both processors 240' and 240.
FIG. 3 (a) illustrates the subband spectrum for even subbands in accordance with the present invention.
A horizontal axis 310' denotes frequency, ω, while a vertical axis 320' denotes the magnitude of the spectrum XβVθn(ω) of the even subband values just prior to recombining with the odd subband values to form the output signal y(n) .
The end-to-end filter response or transfer function is Xθven (ω) /X (ω) using the notation of FIGs 2(a)-2(c) .
SBO, SB2 and SB4 denote the zeroeth, second and fourth subbands respectively. SBO extends between 0 and π/M' while its complex conjugate, SBO*, extends from -τr,/M to 0. SB2 extends between π/M and 37I:/M, and SB3 extends between 3π/M and 5π/M. SB2* extends between - /M and -3-π/M, although part of it is not shown .
The actual passbands of the even subbands are shown at 337, 338, 340 and 342, while the corresponding idealized passbands, having a width of π /M, are shown at 337', 338', 340' and 342', respectively.
FIG. 3(b) illustrates the subband spectrum for odd subbands in accordance with the present invention. A horizontal axis 310 denotes frequency, ω, while a vertical axis 320 denotes the magnitude of the spectrum Xodd(ω) of the odd subband values just prior to recombining with the even subband values to form the output signal y(n) . The end-to-end response or transfer function is Xodd (ω) /χ (ω) using the notation of FIGs 2 (a) - (c) .
SBl, SB3 and SB5 denote the first, third and fifth subbands, respectively. For example, SBl extends between 0 and 2π/M, SB3 extends between 2π/M and 4π/M, and SB5 extends between 4π/M and 6π/M. SBl* denotes the complex conjugate of SBl, and extends between 0 and -2π/M. Not all subbands are shown. For example, with M/2=16, subbands SBl, SB3 , SB5 , SB7 , SB9, SB11, SB13 and SB15 are present in the odd subband path. The actual passbands of the odd subbands are shown at 330, 332, 334 and 336, while the corresponding idealized passbands, having a width of π /M and a sharp cut-off response, are shown at 330', 332', 334', and 336', respectively. Referring to FIGs 3 (a) and 3 (b) , there are two important points to note. First, unlike the critically-sampled counterpart, the pass bands 330, 332, 334, 336, 338, 340, 342,. . . of the oversampled subbands SBl*, SBl, SB3 , SB5 , SBO, SB2 and SB4 , . . ., respectively, occupy less than the available bandwidth, 2π/M. In contrast, a critically-sampled signal utilizes the entire subband bandwidth.
With the present invention, each subband has very little signal energy around its band edges 331, 333,
335, 350, 352, 354 ... Therefore, little or negligible aliasing between subbands is introduced by the filter banks 214, 214' and 284. The oversampling ratio is determined by the ratio of the available bandwidth to the utilized passband, and is usually chosen to be two to minimize computation and implementation complexity, although other ratios can be used.
Second, each of the subband paths 210, 250 has only M/2 subbands, but when we combine the two paths, there are a total of M subbands (output from the filter bank 284) , which cover the full bandwidth of the signal space, and thus the filters 200 or 200' overall each comprise an M-band oversampled filter bank.
Referring again to FIG. 2(b), in the one- dimensional (1-D) filter 200, the analysis filter banks 214, 214' are complements of one other, as are the synthesis filter banks 216, 216'. For the 1-D oversampled filter banks, the following design considerations should be taken into account. Namely, the difference in band alignment in the even and odd subband paths implies that two different filter banks need to be realized. Moreover, the even subband path is not easily realizable conventionally, since the subband centers are at multiples of π/M, which is where the downsampling/aliasing occurs with most fast filter banks algorithm.
However, the 2x oversampled filter banks of FIGs 2(b) and 2(c) solve the aforementioned problems. In particular, the M/2 samples of x(n) are replicated in two paths, with one path (even path 210) going to the frequency shifter 212 of π/M radians, and the other path (odd path 250) going through the delay line 252, which compensates for the delay introduced by the frequency shifter 212. Two sets of identical oversampled M/2 -band analysis filter banks 214, 214' are used for decomposing the even and odd sequences . For M/2 time domain samples, this results in 2*M/2 = M subband samples total . After the subband processing at the subband processor 240, the subbands samples are reassembled/synthesized into the odd and even time- domain sequences with M/2 samples using two sets of identical oversampled M/2 -band synthesis filter banks 216, 216' .
The synthesized even sequence is frequency shifted back by -π/M radians at the shifter 218. The synthesized odd sequence is delayed at the delay 258 to align with the even path 210. The odd and even sequences are added together at the combiner 260 to form M/2 time domain samples, y(n) . Note that x(n) consists of a sequence of M/2 samples, for which there will be a corresponding M/2 samples of y(n), which is the reconstructed version of x(n) . Oversampling occurs since M/2 samples of x(n) are transformed into M samples of sβvβn and sodd combined in the subband domain. Thus, the transformation is considered to be over-complete since it uses more samples than is fundamentally necessary to represent the original sequence .
The frequency shifters 212, 218 can be realized in a variety of ways, as discussed below.
FIG. 4(a) illustrates a critically-sampled cosine- modulated frequency upshifter in accordance with the present invention.
The frequency upshifter 212 ' receives the input signal x(n), having the spectral representation X(ω), and outputs the signal x' (n) , having the spectral representation X' (ω)=χ_π/M) •
This approach to frequency shifting uses an M-band (or multiple thereof) cosine-modulated filter bank to perform the frequency shifting. At an analysis filter portion 405 of the frequency shifter 212', a time- domain transfer function h(n) 410 decomposes the time sequence x(n) into M subband samples, e.g., in paths 420, 430, 440,... , 450. For example, SBO data is provided in path 420, SB2 data is provided in path 430, and so forth. In each path, the subband samples are downsampled by a factor of M at a corresponding downsampler 422, 432, 442,..., 452, and cosine-modulated at a corresponding modulator 424, 434, 444, ...454. The samples from the last (Mth) subband*' (cosine modulation index M-l) in path 450 are discarded. This is acceptable for many applications, such as wide-band audio applications, because the last subband usually contains noise with little signal energy. For example, if an audio signal having a spectrum from 0 to 24 KHz is represented by M=32 subbands, the last band that is discarded has components near 24 KHz, which do not carry significant information for the human listener. Moreover, the discarded subband provides a guard band that can avoid channel interference.
However, if desired, a frequency rotation scheme can be employed where the highest subband output from the analysis filter 405, in path 450 is fed to the modulator 426 of the lowest subband of the synthesis filter 406, in path 420. All information in x(n) can therefore be retained.
The subband samples in each path are fed to the next higher subband at the synthesis filter 406 of the frequency shifter 212', and the zeroeth subband is set to zero at modulator 426. Frequency axis reversal functions 461-463 are used due to the structure of the cosine-modulated filter banks, as discussed further in connection with FIG. 4(b) . Frequency axis reversal for a band-limited sequence can be expressed as: x' (n) = x(n) *exp ( j*π*n) . In an actual implementation, this is equivalent to: x' (n) = x(n)*(-l)n, which flips the sign (polarity) of every other sample. This function can be incorporated into the synthesis filter bank 406, but is shown as being separate in FIG. 4(a) for clarity.
The samples of SBO, in path 420, are provided to modulator 436, in path 430. The samples of SB2 , in path 430, are provided to modulator 446, in path 440, and so forth. Generally, the samples in each subband in the analysis filter bank 405 are fed to the next higher subband in the synthesis filter bank 406.
Cosine modulators 426, 436, 446, ... , 448 demodulate the subband samples from the next lower subband of the analysis filter 405 and provide them to respective upsamplers 428, 438, 448, ... 458, for upsampling by a factor of M. The resulting data is provided to a function f (n) 460, which re-assembles the samples in proper sequence to obtain the time sequence x' (n) , which is shifted by **χ/M radians relative to x(n) . This method is fairly efficient, but the delay is limited by the delay of the M-band filter banks.
A second approach for frequency shifting uses an M-band (or multiple thereof, e.g., 2M, 3M, ...) Discrete Fourier Transform (DFT) type of filter bank such as used in Time Domain Aliasing Cancellation
(TDAC) . See J.P. Princen et al . , Analysis /Synthesis filter bank design based on time domain aliasing cancellation , ASSP-34, Oct 1986. The structure is similar to the first approach. Specifically, this approach is slightly slower than the cosine-modulated filter bank, but has the advantage that we can choose the number of subbands without worrying about frequency reversal, discussed above.
With the second approach, we assume a N*M band TDAC filter bank, where N is an integer. To*' realize a frequency shift of π/M, analysis/synthesis filters are provided back-to-back with the following relationships: sa(0) --> ss(N) sa(l) --> ss(N+l) sa(2) --> ss(N +2)
sa (i) --> ss (N + i) , where sa(i) is the analysis subband output at index i, and ss(i) is the synthesis subband input at index i. FIG. 4(b) illustrates frequency axis reversal for a frequency shifter in accordance with the present invention. Frequency axis reversal is an inherent characteristic of cosine-modulated filter banks. The spectrum of X(ω) is shown at 470. After processing at the analysis filter 405 of the frequency shifter 212', each (non-discarded) subband undergoes frequency axis reversal, as shown at 475. To correct this, a corresponding reversal must be implemented, as shown at 480. The corrected subbands can then be input to the synthesis filter 406 of the frequency shifter
212' to provide the frequency-shifted subbands X {ω-π/M) as shown at 480.
Essentially, a first set of every other subband output from the analysis filter bank 405 contains signals that are frequency axis flipped, i.e., 0→π and π_->0. A second set of the remaining alternate subbands will be frequency axis flipped after processing at the synthesis filter 406. However, by routing the subband data into the next higher subband for frequency upshifting, we are feeding non- flipped signa'ls into the synthesis filter bank, which would attempt to flip these subband signals. Proper frequency axis reversal for each subband avoids this problem.
An analogous result is achieved for frequency downshifting, as discussed in connection with FIG. 6, with frequency axis reversal functions 661-663.
FIG. 5 illustrates a Hilbert Transform frequency shifter in accordance with the present invention.
A third approach for implementing the frequency shifters 212 or 218 of FIGs 2(b) and 2(c) uses a Hilbert transformer with a complex modulator. Specifically, the frequency shifter 212'' includes a Hilbert Transformer 510, a complex modulator 520, and a function 530 for obtain the real portion of the output of the modulator 520.
An ideal Hilbert transformer is an all-pass filter than imparts a 90° phase shift on the input signal. The frequency response is specified as:
H(ω)= -j for 0 < ω < π = j for -π < ω < 0.
The Hilbert transformer can be implemented using a Fast Fourier Transform (FFT) or a Finite Impulse Response (FIR) filter, for example. The Hilbert transform of the sequence is then modulated with a complex modulator 520 with center frequency at 71/ , to achieve a frequency upshift. The real part of the modulated sequence represents the frequency- shifted version of the original time-domain sequence.
The frequency shifter 212 ' ' can be used to achieve a frequency downshift by using exp(-jπ/M) instead of exp(jπ/M) at the modulator 520. The input signal to the frequency shifter is then xβvβ„ ' (n) , and the output signal is xΘVθn(n) .
FIG. 6 illustrates a critically-sampled cosine- modulated frequency downshifter in accordance with the present invention.
The frequency downshifter 218' receives the input signal xβVθn' (n) , having the spectral representation χ even' (ω) > and outputs the signal xβvβn(n), having the spectral representation Xθven (ω) =XΘVβn ' (ω+π/M) .
The frequency downshifter 218', which includes subband paths 620, 630, 640,..., 650, is analogous to the frequency upshifter 212'. An analysis filter bank 605 and synthesis filter bank 606 correspond to the analysis filter bank 405 and synthesis filter bank 406, respectively, of FIG. 4(a), except each modulator 426, 436, 446, ..., 448 receives subband samples from the next higher (instead of lower) subband from the analysis filter bank 605. Frequency axis reversal functions 661-663 are also provided. They are analogous to the reversal functions 461-463, as discussed in connection with FIGs 4(a) and 4(b) . FIG. 7 illustrates a prototype filter design for analysis and synthesis filter banks in accordance with the present invention.
In accordance with a critical aspect of the present invention, the oversampled M/2-band filter banks 214, 214', 216 and 216', along with the complex filters 284 and 286, may be implemented as comb filters. A comb filter may be viewed as a filter in which nulls occur periodically across the frequency band. The comb filter structure of the present invention is realized using a modified form of cosine- modulated filter bank, which can be implemented with fast transforms such as the DCT or FFT. Other uniform filter banks may be used, such as a pseudo-Quadrature Mirror Filter (QMF), as described in J . H. Rothweiler, Polyphase Quadrature Fil ters - A New Subband Coding
Technique, International Conference on Acoustics, Speech and Signal Processing (ICASSP) , 1983. TDAC can also be used with appropriate modification. Consider an M/2 -band cosine-modulated filter bank.
In a conventional, critically-sampled implementation, the prototype filter is designed to have a bandwidth of π. In the present invention, the prototype filter is designed to have a bandwidth of π/2. The response 730 of the corresponding comb filter is shown in FIG 7. A horizontal axis 710 indicates frequency, while a vertical axis 720 indicates the magnitude of the response, in decibels (dB) . After the cosine modulation, a series of non-overlapping band-pass filters result which have band edges at multiples of π/4 and 3π/4, respectively. The effective bandwidth of the pass band 750 is therefore π/2. Note that an attenuation of -40 dB can be considered to represent an amplitude of essentially zero. The stop band rejection of the filter is designed to be low enough so that aliasing is small or negligible at 0 and π. Stop bands 740 and 760 are provided between 0 and π/4 , and between 3π/4 and π» respectively. An additional benefit when designing this prototype filter is that we are no longer constrained by the severe aliasing cancellation matrix, discussed in P.P. Vaidyanathan, Mul tirate System And Fil ter Banks . Prentice Hall 1993, as with the critically- sampled case, thus yielding additional freedom in the design. Effectively, the filter order can be reduced substantially to provide better temporal resolution, and to achieve improved stop band rejection characteristics . Note that the response 730 applies to both the analysis and synthesis filter banks since, in the design of a cosine-modulated filter bank, one only needs to specify the prototype filter characteristics, and the response of the analysis and synthesis filters will be dictated by the prototype filter.
FIG. 8 illustrates a comb filter bank structure in accordance with the present invention.
The circuit 800 can perform the functions of the filter 214, subband processor 240, and filter 216 of FIG. 2(b), or the filter 214', subband processor 240, and filter 216' of FIG. 2(b).
An analysis comb filter bank 214' ' includes a transfer function h(n) 810 for providing the input signal in M/2 subbands, downsamplers 820, 822, ...,
824, 826 for each of the M/2 even or odd subbands, and corresponding cosine modulators 830, 832, ... , 834, 836. After processing at the subband processor 240', the even or odd subband data is processed at a synthesis comb filter bank 216'', which includes cosine modulators 840, 842, ..., 844, 846, corresponding M/2 upsamplers 850, 852, ..., 854, 856, and a transfer function f (n) 860 that assembles the subband data to provide the output xβvβn' (n) or xodd' (n) . FIG. 9 illustrates a graphical realization of one- dimensional filter bank processing in accordance with the present invention.
The graphical illustration shows how the oversampled filter banks are realized with the frequency shifter and comb filter banks already discussed.
For simplicity, only a portion of the actual spectrum X(©) (915) of the input signal x(n), i.e., from 0 to 5π/M is shown. A number of subband portions A-J are shown.
For the processing of the odd subbands, Xd(ω) (920) is formed by delaying X(ω) . The passbands of the oversampled filter banks are shown at 922. The resulting odd subband samples Sodd(ω) (924), which includes subband portions B, C, F, G and J, occupy only one-half of the available subband width, 2π/M. After a delay, Xodd(co) (926) is formed, whose spectrum is essentially the same as Sodd(ω) (924) .
For the processing of the even subbands,. X' (ω) (930) is formed by frequency- shifting X(ω) by π/M. The passbands of the oversampled filter banks are shown at 932. The resulting even subband samples SΘVen(ω) (934), including subband portions A, D, E and H, occupy only one-half of the available subband width, 2π/M. After a -JC/M frequency downshift, Xβvβn(ω) (936) is formed.
The subband values Xodd(ω) and Xβvβn(ω) are then combined to form the full -band digital output signal
Y(ω) (950).
FIG. 10 illustrates a 2M-band filter (e.g., four, M/2-band filters) with 4x-oversampled filter banks with four subband paths in accordance with the present invention.
The filter structure 1000 includes four subband processing paths, i.e., Path 0 (1010), Path 1 (1020), Path 2 (1030) , and Path 3 (1040) . Path 0 (1010) includes a π/M frequency shifter 1012, an M/2 -band analysis filter bank 1014, an M/2 -band synthesis filter
1016, and a ~π/M frequency shifter 1018. Each synthesis filter bank receives corresponding data from the subband processor 240''. Path 1 (1020) includes a π/2M frequency shifter 1012', an M/2-band analysis filter bank 1014', an M/2-band synthesis filter 1016', and a -π/2M frequency shifter 1018'. Note that the frequency shift functions described herein can be implemented as discussed previously (e.g., in connection with FIGs 4(a), 4(b), 5 and 6) with modification as required for the magnitude and direction of frequency shift.
Path 2 (1030) includes a delay 252, an M/2-band analysis filter bank 1014'', an M/2 -band synthesis filter 1016' ' , and a delay 252' . The delay 252 compensates for the processing time of the frequency shifters 1012, 1012', and 1012'', while the delay 252' compensates for the processing time of the frequency shifters 1018, 1018', and 1018''.
Path 3 (1040) includes a -π/2M frequency shifter 1012' ' , an M/2-band analysis filter bank 1014' ' , an M/2 -band synthesis filter 1016'', and a π/2M frequency shifter 1018 ' ' .
With this scheme, each successive subband is processed in a successive one of the Paths 1010, 1020, 1030, or 1040. For example, assuming 32 subbands are used, i.e., numbered SB0-SB31, SBO, SB4 , SB8, ... are processed in Path 0 (1010), SBl, SB5 , SB9 , ... are processed in Path 1 (1020), SB2 , SB6, SB10, ... are processed in Path 2 (1030), and SB3 , SB7, SB11, ... are processed in Path 3 (1040) . The filter 1000 could optionally be modified to use one or two complex filters, instead of four real filters, based on the discussion of FIG. 2(c) .
FIG. 11(a) illustrates the subband structure for Path 0 of the filter bank of FIG. 10 in accordance with the present invention. The pass bands 1112, 1116, 1120, ... of the 4x oversampled subbands SBO*, SBO, SB4, SB8, ... occupy only about one-fourth of the available bandwidth, 2π/M. The corresponding idealized pass bands are shown at 1112', 1116', 1120', ... , respectively. Specifically, only about π/2M radians is occupied by each subband, which leaves a substantial region with little signal energy around the band edges 1110, 1114, 1118, ... Therefore, little or negligible aliasing between subbands is introduced by the filter structure 1000. XPath-o(ω) designates the signal in Path 0 (1010), as shown in FIG. 10.
A similar result is achieved for the other subbands, as discussed below. Generally, a primary benefit of 4x oversampling that there is even less of a chance of aliasing than with 2x oversampling. The present inventor has found that -40dB attenuation at the band edges is more than enough for most practical purposes. However, there may be cases where the application demands even more attenuation. One way to achieve this is to increase the prototype filter order and trade off temporal resolution. The second way, discussed here, is to use an even higher oversampling ratio, which results in more significant filter rolloff at the band edges.
FIG. 11(b) illustrates the subband structure for Path 1 of the filter bank of FIG. 10 in accordance with the present invention.
Here, pass bands 1132, 1136, 1140, 1144, ... of the 4x oversampled subbands SBl*, SBl, SB5 , SB9 , ... are provided, with band edges 1130, 1134, 1138, ... The corresponding idealized pass bands 1132', 1136', 1140', 1144', ... are also shown. χpth-ι(ω) designates the signal in Path 1 (1020), as shown in FIG. 10.
FIG. 11(c) illustrates the subband structure for Path 2 of the filter bank of FIG. 10 in accordance with the present invention.
Here, pass bands 1152, 1156, 1160, ... of the 4x oversampled subbands SB2 , SB6, SB10, ... are provided, with band edges 1150, 1154, 1158, ... The corresponding idealized pass bands 1152', 1156', 1160', ... are also shown.
X Pth-2(oc>) designates the signal in Path 2 (1030), as shown in FIG. 10.
FIG. 11(d) illustrates the subband structure for Path 3 of the filter bank of FIG. 10 in accordance with the present invention.
Here, pass bands 1182, 1186, ... of the 4x oversampled subbands SB2 , SB6, SB10, ... are provided, with band edges 1180, 1184, 1188, ... The corresponding idealized pass bands 1182' , 1186' , ... are also shown. χpah-3(ω) designates the signal in Path 3 (1040), as shown in FIG. 10.
FIG. 12(a) illustrates an original 2-D signal (image) region with several subregions for processing in accordance with the present invention.
The image 1200 is shown as including a number of subregions A-P. Typically, image data is processed by subdividing it into several blocks, each of which includes many pixels. A block-based spatial transform, such as the Discrete Cosine Transform (DCT) is often used, to provide frequency domain coefficients. Other possible transforms include the Discrete Fourier
Transform, Karhunen-Loeve Transform, Walsh Hadamard Transform, and wavelet transform, as well as other known transforms.
FIG. 12(b) provides a graphical depiction of a conventional 2-D, critically-sampled filter bank. Note that FIGs 12(b) -13(e) are diagrammatic illustrations that show a composite of image data in the frequency domain (such as DCT frequency coefficients or the like) and associated filters. Each filter in the image/filter composite 1200' is labeled with (x,y) coordinates, and is overlaid on the corresponding image data subregion which it processes. Note that x and y are used to denote horizontal and vertical direction, respectively, and should not be confused with the input or output signals x(n) and y(n) . The meaning should be clear from the context.
An 8x8 2-D filter bank is assumed for illustration purposes, although different sizes can be used.
For example, an image subregion 1210 has data that is processed by a filter labeled (7,7) . The image data within the subregion 1210, which is part of the larger image subregion "P", is passed through a pass band 1212 of the corresponding filter, while image data beyond the band edges 1214 will not be passed by this filter. However, the 2-D, critically-sampled filter bank suffers from the same problem of aliasing between subbands as discussed for the 1-D critically sampled filter bank.
FIG. 13(a) provides a graphical depiction of a path with no frequency shift (Path 0) in a 2-D, over- sampled filter bank in accordance with the present invention .
As explained further in connection with FIG. 14, a 2x oversampled 2-D filter bank structure in accordance with the present invention processes image or other 2-D data in four processing paths, namely Paths 0-3.
In a particular embodiment, the 2-D data toward the center of each image region is processed in the first path, Path 0, while data at vertical boundaries between regions (x-axis shifted by π/M) is processed in a second path, Path-1, data at horizontal boundaries between regions (y-axis shifted by π/M) 1S processed in a third path, Path-2, and data at corners of the boundaries between regions (x-axis and y-axis shifted by π/M) 1S processed in a fourth path, Path-3. Here, each filter is labeled with (P,x,y) coordinates in the image/filter composite 1300-0, where the "P" coordinate denotes the path, and is shown overlaid with the corresponding image data subregion. For example, an image subregion 1310, labeled "H", has data that is processed by an oversampled filter (0,1,3) .
The image data toward the center of the subregion 1310 will pass through a pass band 1312 of the corresponding filter essentially unchanged, while data between the pass band 1312 and band edges 1316 is stopped by a stop band 1314.
FIG. 13(b) provides a graphical depiction of a path with a horizontal frequency shift (Path 1) in a 2- D, over-sampled filter bank in accordance with the present invention.
In this path, the image/filter composite 1300-1 shows that data at vertical boundaries between regions (x-axis shifted by π/M) 1S processed using the corresponding filters as shown. For example, data between the regions "A" and "B" is processed using a filter with index (1,0,1).
FIG. 13(c) provides a graphical depiction of a path with a vertical frequency shift (Path 2) in a 2-D, over-sampled filter bank in accordance with the present invention.
In this path, the image/filter composite 1300-2 shows that data at horizontal boundaries between regions (y-axis shifted by π/M) is processed using the corresponding filters as shown. For example, data between the regions "A" and "E" is processed using a filter with index (2,1,0).
FIG. 13(d) provides a graphical depiction of a path with horizontal and vertical frequency shifts (Path-3) in a 2-D, over-sampled filter bank in accordance with the present invention.
In this path, the image/filter composite 1300-3 shows that data at corners of the boundaries between regions (x-axis and y-axis shifted by π/M) is processed in a fourth path. For example, data at the corner of regions "A", "B", "E" and "F" is processed using a filter with index (3,1,1).
FIG. 13(e) provides a graphical depiction of the freqμency regions that are processed in FIGs 13 (a) - 13(d) in a 2-D, over-sampled filter bank in-' accordance with the present invention.
The image/filter composite 1300 shows the aggregation of the 2-D data regions that are processed by each of the four paths. In the example shown, sixteen regions are processed in each path, for a total of sixty-four regions.
A region 1385 illustrates an area of imperfect reconstruction by the filter banks. Thus, by summing the regions from FIGs 13(a) , (b) , (c) and (d) , we find that the area 1385 is not covered by any of the four paths. To obtain complete coverage, additional processing paths can be introduced. This is similar to the 1-D case, where we decided to drop the last subband. Notice that the area 1385 corresponds to high- requency regions of the 2-D signal, so, for most applications, this data is not used at all and can be discarded without noticeably impairing the resulting output signal quality.
Generally, the present invention can be extended to an oversampled filter bank structure for use with N- dimensions . For illustration purposes, an example of a 2-D oversampled filter bank has been shown. Such a filter bank can be used, e.g., in applications relating to video or image processing. FIG. 14 illustrates a real filter structure for two-dimensional filter bank processing with 2x- oversampled real filter banks in four subband paths in accordance with the present invention.
Note that the frequency shift functions described herein can be implemented as discussed previously (e.g., in connection with FIGs 4(a), 4(b), -5 and 6) with modification as required for the magnitude and direction of frequency shift.
Additionally, the filter structure can be modified to used one or more complex filters as discussed in connection with FIG. 2(c) .
The filter structure 1400 includes processing paths 1410, 1420, 1430 and 1440. Generally, a number of delay elements (in path 0: 252-0, 252-0", 252-0'',
252-0''*; in path 1: 252-1, 252-1'; and, in path 2: 252-2, 252-2') are provided to compensate for the processing time of the functions that are shown as vertically aligned in the figure.
Each path includes an analysis filter bank 1414-0,
1414-1, 1414-2 and 1414-3, and a corresponding synthesis filter bank 1416-0, 1416-1, 1416-2 and 1416-
3 , which receives the corresponding data from the subband processor 240' ' ' .
Path 1 (1420) includes counterpart x-axis frequency shifters 1412-1 and 1418-1, while Path 2 (1430) includes counterpart y-axis frequency shifters
1413-2 and 1417-2, and Path 3 (1440) includes both counterpart x-axis frequency shifters 1412-3 and 1418-
3, and counterpart y-axis frequency shifters 1413-3 and
1417-3. A combiner 260' receives the output in each path to provide the final output signal y(n) . Again, the x- axis and y-axis terminology should not be confused with the input signal x(n) and output signal y(n) .
Note that the x-axis (1412-1, 1412-3) and y-axis (1413-2, 1413-3) shift functions can use common hardware and/or software.
Accordingly, it can be seen that the present invention provides an oversampled filter bank structure that can be implemented using popular and efficient fast filter banks. Advantageously, aliasing between subbands is substantially reduced, thus allowing the independent processing of the subband signals using any known subband processing technique. For example, an input audio or video signal may be processed to embed auxiliary data in its subbands. The auxiliary data may be used for any number of purposes, such as copyright protection, broadcast verification, identification of soundtrack, carrying electronic coupon data, and so forth. Moreover, equalizing and noise shaping of the subbands as well as other processing may be provided. Additionally, the oversampled filter bank can be used in OFDM, which has found application in Digital Subscriber Loop (DSL) and digital broadcast systems. In a particular embodiment, even subbands of an input signal are frequency- shifted prior to filtering and subband processing, while the odd subbands are delayed to compensate for the processing time of the frequency shifting. Separate analysis and synthesis filter banks may be provided for the even and odd subbands, or the subbands may be processed together in a single complex filter bank. 2x oversampling may be advantageously used in the filter banks.
Frequency shifting may be achieved using critically-sampled cosine filter banks, or a Hilbert transform, for example.
In an optional embodiment, the subbands are processed in four subband paths, and 4x oversampling is used at the filter banks. The invention may also be extended to two-dimensional or higher filtering applications.
Moreover, the invention can be extended to multidimensional applications. For example, a 2-D oversampled filter bank structure is shown that processes data in four paths, three of which undergo frequency shifts. The structure is suited for processing image or other 2-D data, such as seismic or other vibrational data, temperature data and the like.
The invention can be extended using the principles shown to additional dimensions. Although the invention has been described in connection with various specific embodiments, those skilled in the art will appreciate that numerous adaptations and modifications may be made thereto without departing from the spirit and scope of the invention as set forth in the claims.
For example, known computer hardware, firmware and/or software techniques may be used to implement the invention.
Also, the digital input signal discussed herein may only be a portion of an overall signal that is provided at some upstream location. Only a portion of the overall signal need be processed in accordance with the invention. Additionally, the system of the invention may be aggregated to process additional input signals, or additional portions of an input signal.

Claims

What is claimed is :
1. A digital filter apparatus for processing a digital input signal, comprising: oversampled analysis filter bank means"' for filtering at least a first and second plurality of subbands of said digital input signal according to an oversampling ratio; wherein: said second plurality of subbands are distinct from said first plurality of subbands; and the oversampling ratio is selected to reduce aliasing between adjacent ones of the subbands.
2. The apparatus of claim 1, wherein: the oversampling ratio is 2:1.
3. The apparatus of claim 1, wherein: said first plurality of subbands comprise even- numbered subbands of said digital input signal; and said second plurality of subbands comprise odd- numbered subbands of said digital input signal .
4. The apparatus of claim 1, wherein: said filter bank means comprises at least one comb filter.
5. The apparatus of claim 1, wherein: said filter bank means outputs M subbands, where M
> 8.
6. The apparatus of claim 1, further comprising: a first frequency shifter for frequency-shifting the first plurality of subbands of said digital input signal to provide corresponding frequency-shifted subbands for filtering by said filter bank "means .
7. The apparatus of claim 6, wherein: said first frequency shifter comprises critically- sampled cosine-modulated analysis and synthesis filter banks .
8. The apparatus of claim 6, wherein: said first frequency shifter comprises a Hilbert trans ormer.
9. The apparatus of claim 6, further comprising: first and second parallel paths for carrying said digital input signal to said filter bank means; wherein said first frequency shifter is provided in said first path; and a first delay is provided in said second path for delaying the second plurality of subbands to compensate for a delay caused by a processing time of said first frequency shifter.
10. The apparatus of claim 9, wherein: said filter bank means comprises a first oversampled analysis filter bank arranged in said first path for said filtering of said frequency-shifted subbands, and a second oversampled analysis filter bank arranged in said second path for said filtering of the delayed subbands .
11. The apparatus of claim 10, wherein: said first and second analysis filter ioanks comprise respective M/2 -band filter banks, and M is the number of subbands output from the filter bank means.
12. The apparatus of claim 6, further comprising: a subband processor for processing the frequency- shifted subbands and the second plurality of subbands after filtering thereof.
13. The apparatus of claim 12, further comprising: oversampled synthesis filter means for filtering said frequency-shifted subbands after processing thereof at said subband processor, and for filtering said second plurality of subbands after processing thereof at said subband processor.
14. The apparatus of claim 13, further comprising: a second frequency shifter for frequency-shifting said frequency-shifted subbands after filtering thereof at said synthesis filter means in a shift direction opposite to a shift direction of said first frequency shifter.
15. The apparatus of claim 14, wherein: said first frequency shifter provides a frequency upshift, and said second frequency shifter provides a corresponding frequency down shift .
16. The apparatus of claim 14, further comprising: a second delay for delaying the second plurality of subbands to compensate for a delay of said second frequency shifter.
17. The apparatus of claim 13, further comprising: a combiner for combining said frequency-shifted subbands after filtering thereof at said synthesis filter means, with said second plurality of subbands after filtering thereof at said synthesis filter means, to provide a digital output signal corresponding to said digital input signal.
18. The apparatus of claim 13, wherein: said synthesis filter means comprises a first oversampled synthesis filter bank for said synthesis filtering of said frequency-shifted subbands, and a second oversampled synthesis filter bank for said synthesis filtering of said second plurality of subbands .
19. The apparatus of claim 18, wherein: said first and second oversampled synthesis filter banks comprise respective comb filters.
20. The apparatus of claim 6, wherein: said first frequency shifter provides a π/M radian frequency shift; and
M is the number of subbands output from said filter bank means.
21. The apparatus of claim 6, wherein said first frequency shifter comprises: a cosine-modulated analysis filter bank for processing the first plurality of subbands; means for providing frequency axis reversal for the subbands output from the cosine-modulated analysis filter bank; and a cosine-modulated synthesis filter bank for processing the subbands output from the frequency axis reversal means to provide the frequency- shifted subbands .
22. The apparatus of claim 6, wherein said first frequency shifter comprises: a cosine-modulated analysis filter bank for processing the first plurality of subbands; a cosine-modulated synthesis filter bank for processing the subbands output from the cosine- modulated analysis filter bank to provide the frequency- shifted subbands; and means for feeding a highest subband that is output from the cosine-modulated analysis filter bank to a lowest subband input of the cosine-modulated synthesis filter bank.
23. The apparatus of claim 6, wherein said filter bank means comprises a complex filter bank, further comprising : a delay for delaying the second plurality of subbands to compensate for a delay of said first frequency shifter; and means for forming a complex signal comprising a real portion obtained from the frequency-shifted subbands, and an imaginary portion obtained from the delayed subbands ; and means for providing the complex signal to the complex filter bank for filtering thereat.
24. The apparatus of claim 1, wherein: said oversampled analysis filter bank means is adapted to filter additional subbands of said digital input signal, including at least a third and fourth plurality of subbands of said digital input signal according to the oversampling ratio,- wherein: said third plurality of subbands are distinct from said fourth plurality of subbands .
25. The apparatus of claim 24, wherein: said oversampling ratio is 4:1.
26. The apparatus of claim 24, further comprising: first, second, third and forth parallel paths for providing said digital input signal to said analysis filter bank means,- a first frequency shifter in said first path for frequency-shifting said first plurality of subbands,- a second frequency shifter in said second path for frequency-shifting said second plurality of subbands,- a third frequency shifter in said fourth path for frequency-shifting said third plurality of subbands,- and a delay in said third path for delaying said fourth plurality of subbands to compensate for a processing time of said first, second and third frequency-shif ers .
27. The apparatus of claim 26, wherein: said first frequency shifter provides a π/M radian frequency shift; said second frequency shifter provides a π/2M radian frequency shift; said third frequency shifter provides a -π/2M radian frequency shift; and
2M is the number of subbands output from the oversampled analysis filter bank means.
28. The apparatus of claim 26, wherein.- said analysis filter bank means comprises M/2 -band oversampled filter banks for each of said first, second, third and fourth paths, and 2M is the number of subbands output from the oversampled analysis filter bank means .
29. The apparatus of claim 28, wherein: said M/2 -band oversampled filter banks comprise comb filters.
30. The apparatus of claim 24, wherein the digital input signal comprises two-dimensional data extending along an x-axis and a perpendicular y-axis, further comprising: first, second, third and forth parallel paths for providing said digital input signal to said analysis filter bank means; x-axis frequency shifting means for frequency- shifting said second plurality of subbands in said second path, and said fourth plurality of subbands in said fourth path; y-axis frequency shifting means for frequency- shifting said third plurality of subbands in said third path, and said fourth plurality of subbands in said fourth path,- and delay means for delaying said first plurality of subbands in said first path to compensate for a processing time of said x-axis and y-axis frequency shifting means.
31. The apparatus of claim 30, wherein: said x-axis and y-axis frequency shifting means provide a magnitude π/M radian frequency shift; and
4xMxM is the number of subbands output from the oversampled analysis filter bank means .
32. The apparatus of claim 30, wherein: said analysis filter bank means comprises M/2 -band oversampled filter banks for each of said first, second, third and forth paths, and 4xMxM is the number of subbands output from the oversampled analysis filter bank means .
33. The apparatus of claim 32, wherein: said M/2 -band oversampled filter banks comprise comb filters.
34. A digital filtering method for a digital input signal, comprising the steps of: filtering at least a first and second plurality of subbands of said digital input signal at an oversampled analysis filter bank means according to an oversampling ratio,- wherein: said second plurality of subbands are distinct from said first plurality of subbands; and the oversampling ratio is selected to reduce aliasing between adjacent ones of the subbands.
PCT/US1999/028343 1999-11-30 1999-11-30 Oversampled filter bank for subband processing WO2001041320A1 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US10/148,161 US6947509B1 (en) 1999-11-30 1999-11-30 Oversampled filter bank for subband processing
PCT/US1999/028343 WO2001041320A1 (en) 1999-11-30 1999-11-30 Oversampled filter bank for subband processing
AU18372/00A AU1837200A (en) 1999-11-30 1999-11-30 Oversampled filter bank for subband processing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/US1999/028343 WO2001041320A1 (en) 1999-11-30 1999-11-30 Oversampled filter bank for subband processing

Publications (1)

Publication Number Publication Date
WO2001041320A1 true WO2001041320A1 (en) 2001-06-07

Family

ID=22274187

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US1999/028343 WO2001041320A1 (en) 1999-11-30 1999-11-30 Oversampled filter bank for subband processing

Country Status (2)

Country Link
AU (1) AU1837200A (en)
WO (1) WO2001041320A1 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002080362A1 (en) * 2001-04-02 2002-10-10 Coding Technologies Sweden Ab Aliasing reduction using complex-exponential modulated filterbanks
EP1494353A2 (en) * 2003-06-09 2005-01-05 Northrop Grumman Corporation Oversampled filterbank with near perfect reconstruction constraint
US8484272B2 (en) 2004-08-20 2013-07-09 Qualcomm Incorporated Unified pulse shaping for multi-carrier and single-carrier waveforms
CN108140396A (en) * 2015-09-22 2018-06-08 皇家飞利浦有限公司 Audio signal processing

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4691292A (en) * 1983-04-13 1987-09-01 Rca Corporation System for digital multiband filtering
US5568142A (en) * 1994-10-20 1996-10-22 Massachusetts Institute Of Technology Hybrid filter bank analog/digital converter
US5732189A (en) * 1995-12-22 1998-03-24 Lucent Technologies Inc. Audio signal coding with a signal adaptive filterbank
US5745392A (en) * 1995-10-05 1998-04-28 Chevron U.S.A. Inc. Method for reducing data storage and transmission requirements for seismic data

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4691292A (en) * 1983-04-13 1987-09-01 Rca Corporation System for digital multiband filtering
US5568142A (en) * 1994-10-20 1996-10-22 Massachusetts Institute Of Technology Hybrid filter bank analog/digital converter
US5745392A (en) * 1995-10-05 1998-04-28 Chevron U.S.A. Inc. Method for reducing data storage and transmission requirements for seismic data
US5732189A (en) * 1995-12-22 1998-03-24 Lucent Technologies Inc. Audio signal coding with a signal adaptive filterbank

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
RAMSTAD T.A. ET AL.: "Cosine-Modulated Analysis-Synthesis Filtering with Critical Sampling and Perfect Reconstruction", IEEE INTERNATIONAL CONFERENCE ON ACOUSTICS, SPEECH AND SIGNAL PROCESSING, ICASSP-91, vol. 3, 1991, pages 1789 - 1792, XP002927264 *
SMITH M.J.T. ET AL.: "Analysis/Synthesis Techniques for Subband Image Coding", IEEE TRANSACTIONS ON ACOUSTICS, SPEECH AND SIGNAL PROCESSING, vol. 38, no. 8, August 1990 (1990-08-01), pages 1446 - 1456, XP002927263 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002080362A1 (en) * 2001-04-02 2002-10-10 Coding Technologies Sweden Ab Aliasing reduction using complex-exponential modulated filterbanks
US7242710B2 (en) 2001-04-02 2007-07-10 Coding Technologies Ab Aliasing reduction using complex-exponential modulated filterbanks
EP1494353A2 (en) * 2003-06-09 2005-01-05 Northrop Grumman Corporation Oversampled filterbank with near perfect reconstruction constraint
EP1494353A3 (en) * 2003-06-09 2006-01-11 Northrop Grumman Corporation Oversampled filterbank with near perfect reconstruction constraint
US7372907B2 (en) 2003-06-09 2008-05-13 Northrop Grumman Corporation Efficient and flexible oversampled filterbank with near perfect reconstruction constraint
US8484272B2 (en) 2004-08-20 2013-07-09 Qualcomm Incorporated Unified pulse shaping for multi-carrier and single-carrier waveforms
TWI425797B (en) * 2004-08-20 2014-02-01 Qualcomm Inc Unified pulse shaping for multi-carrier and single-carrier waveforms
CN108140396A (en) * 2015-09-22 2018-06-08 皇家飞利浦有限公司 Audio signal processing

Also Published As

Publication number Publication date
AU1837200A (en) 2001-06-12

Similar Documents

Publication Publication Date Title
US6947509B1 (en) Oversampled filter bank for subband processing
Karp et al. Modified DFT filter banks with perfect reconstruction
AU733587B2 (en) Method and apparatus for sub-band modulation of signals for transmission and/or storage
EP0126148B1 (en) Image processing method using a block overlap transformation procedure
Fernandes et al. Complex wavelet transforms with allpass filters
US7620673B2 (en) Complimentary discrete fourier transform processor
Diniz et al. Design of high-resolution cosine-modulated transmultiplexers with sharp transition band
Viholainen et al. Complex modulated critically sampled filter banks based on cosine and sine modulation
KR20230124676A (en) High-Performance Filter Bank Channelizer
Viholainen et al. Implementation of parallel cosine and sine modulated filter banks for equalized transmultiplexer systems
Furtado et al. Numerically efficient optimal design of cosine-modulated filter banks with peak-constrained least-squares behavior
WO2001041320A1 (en) Oversampled filter bank for subband processing
KR20060003070A (en) A multi-channel tuner using a discrete cosine transform
Harris et al. Interleaving different bandwidth narrowband channels in perfect reconstruction cascade polyphase filter banks for efficient flexible variable bandwidth filters in wideband digital transceivers
Fliege Closed form design of prototype filters for linear phase DFT polyphase filter banks
Harris et al. Cascade non-maximally decimated filter banks form efficient variable bandwidth filters for wideband digital transceivers
Ramachandran et al. Bandwidth efficient transmultiplexers, part 2: Subband complements and performance aspects
Zhong et al. Design of oversampled nonuniform filter banks with arbitrary rational frequency partitioning
Pun et al. A novel and efficient design of multidimensional PR two-channel filter banks with hourglass-shaped passband support
Coulombe et al. Nonuniform perfect reconstruction filter banks over lattices with application to transmultiplexers
Johansson et al. Two-channel FIR filter banks based on the frequency-response masking approach
Venosa et al. Implement Even and Odd Stacked Frequency Bins in Cascade Non-Maximally Decimated Analysis and Synthesis Filter Banks
Sablatash et al. Polyphase Realization of VSB Filter Banks for Spectrally Efficient Transmitter Multiplexers and Receiver Demultiplexers Using Spectral Factorization
Harris et al. How to Perform Very Wideband Digital Filtering in Modern Software Defined Radios
Lin et al. On the study of four-parallelogram filter banks

Legal Events

Date Code Title Description
ENP Entry into the national phase

Ref country code: AU

Ref document number: 2000 18372

Kind code of ref document: A

Format of ref document f/p: F

AK Designated states

Kind code of ref document: A1

Designated state(s): AE AL AM AT AU AZ BA BB BG BR BY CA CH CN CU CZ DE DK EE ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MD MG MK MN MW MX NO NZ PL PT RO RU SD SE SG SI SK SL TJ TM TR TT UA UG US UZ VN YU ZA ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): GH GM KE LS MW SD SL SZ TZ UG ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE BF BJ CF CG CI CM GA GN GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
WWE Wipo information: entry into national phase

Ref document number: 10148161

Country of ref document: US

REG Reference to national code

Ref country code: DE

Ref legal event code: 8642

122 Ep: pct application non-entry in european phase