1. Technical Field
The present invention relates to low bitrate audio source coding systems. Different parametric representations of stereo properties of an input signal are introduced, and the application thereof at the decoder side is explained, ranging from pseudo-stereo to full stereo coding of spectral envelopes, the latter of which is especially suited for HFR based codecs.
2. Description of the Related Art
Audio source coding techniques can be divided into two classes: natural audio coding and speech coding. At medium to high bitrates, natural audio coding is commonly used for speech and music signals, and stereo transmission and reproduction is possible. In applications where only low bitrates are available, e.g. Internet streaming audio targeted at users with slow telephone modem connections, or in the emerging digital AM broadcasting systems, mono coding of the audio program material is unavoidable. However, a stereo impression is still desirable, in particular when listening with headphones, in which case a pure mono signal is perceived as originating from “within the head”, which can be an unpleasant experience.
One approach to address this problem is to synthesize a stereo signal at the decoder side from a received pure mono signal. Throughout the years, several different “pseudo-stereo” generators have been proposed. For example in [U.S. Pat. No. 5,883,962], enhancement of mono signals by means of adding delayed/phase shifted versions of a signal to the unprocessed signal, thereby creating a stereo illusion, is described. Hereby the processed signal is added to the original signal for each of the two outputs at equal levels but with opposite signs, ensuring that the enhancement signals cancel if the two channels are added later on in the signal path. In [PCT WO 98/57436] a similar system is shown, albeit without the above mono-compatibility of the enhanced signal. Prior art methods have in common that they are applied as pure post-processes. In other words, no information on the degree of stereo-width, let alone position in the stereo sound stage, is available to the decoder. Thus, the pseudo-stereo signal mayor may not have a resemblance of the stereo character of the original signal. A particular situation where prior art systems fall short, is when the original signal is a pure mono signal, which often is the case for speech recordings. This mono signal is blindly converted to a synthetic stereo signal at the decoder, which in the speech case often causes annoying artifacts, and may reduce the clarity and speech intelligibility.
Other prior art systems, aiming at true stereo transmission at low bitrates, typically employ a sum and difference coding scheme. Thus, the original left (L) and right (R) signals are converted to a sum signal, S=(L+R)I2, and a difference signal, D=(L−R)/2, and subsequently encoded and transmitted. The receiver decodes the Sand D signals, whereupon the original UR-signal is recreated through the operations L=S+D, and R=S−D. The advantage of this, is that very often a redundancy between L and R is at hand, whereby the information in D to be encoded is less, requiring fewer bits, than in S. Clearly, the extreme case is a pure mono signal, i.e. L and R are identical. A traditional L/R-codec encodes this mono signal twice, whereas a S/D codec detects this redundancy, and the D signal does (ideally) not require any bits at all. Another extreme is represented by the situation where R=−L, corresponding to “out of phase” signals. Now, the S signal is zero, whereas the D signal computes to L. Again, the S/D scheme has a clear advantage to standard L/R-coding. However, consider the situation where e.g. R=0 during a passage, which was not uncommon in the early days of stereo recordings. Both S and D equal L/2, and the S/D-scheme does not offer any advantage. On the contrary, L/R-coding handles this very well: The R signal does not require any bits. For this reason, prior art codecs employ adaptive switching between those two coding schemes, depending on what method that is most beneficial to use at a given moment. The above examples are merely theoretical (except for the dual mono case, which is common in speech only programs). Thus, real world stereo program material contains significant amounts of stereo information, and even if the above switching is implemented, the resulting bitrate is often still too high for many applications. Furthermore, as can be seen from the resynthesis relations above, very coarse quantization of the D signal in an attempt to further reduce the bitrate is not feasible, since the quantization errors translate to non-neglectable level errors in the L and R signals.
The present invention employs detection of signal stereo properties prior to coding and transmission. In the simplest form, a detector measures the amount of stereo perspective that is present in the input stereo signal. This amount is then transmitted as a stereo width parameter, together with an encoded mono sum of the original signal. The receiver decodes the mono signal, and applies the proper amount of stereo-width, using a pseudo-stereo generator, which is controlled by said parameter. As a special case, a mono input signal is signaled as zero stereo width, and correspondingly no stereo synthesis is applied in the decoder. According to the invention, useful measures of the stereo-width can be derived e.g. from the difference signal or from the cross-correlation of the original left and right channel. The value of such computations can be mapped to a small number of states, which are transmitted at an appropriate fixed rate in time, or on an as-needed basis. The invention also teaches how to filter the synthesized stereo components, in order to reduce the risk of unmasking coding artifacts which typically are associated with low bitrate coded signals.
Alternatively, the overall stereo-balance or localization in the stereo field is detected in the encoder. This information, optionally together with the above width-parameter, is efficiently transmitted as a balance-parameter, along with the encoded mono signal. Thus, displacements to either side of the sound stage can be recreated at the decoder, by correspondingly altering the gains of the two output channels. According to the invention, this stereo-balance parameter can be derived from the quotient of the left and right signal powers. The transmission of both types of parameters requires very few bits compared to full stereo coding, whereby the total bitrate demand is kept low. In a more elaborate version of the invention, which offers a more accurate parametric stereo depiction, several balance and stereo-width parameters are used, each one representing separate frequency bands.
The balance-parameter generalized to a per frequency-band operation, together with a corresponding per band operation of a level-parameter, calculated as the sum of the left and right signal powers, enables a new, arbitrary detailed, representation of the power spectral density of a stereo signal. A particular benefit of this representation, in addition to the benefits from stereo redundancy that also S/D-systems take advantage of, is that the balance-signal can be quantized with less precision than the level ditto, since the quantization error, when converting back to a stereo spectral envelope, causes an “error in space”, i.e. perceived localization in the stereo panorama, rather than an error in level. Analogous to a traditional switched L/R- and S/D-system, the level/balance-scheme can be adaptively switched off, in favor of a levelL/levelR-signal, which is more efficient when the overall signal is heavily offset towards either channel. The above spectral envelope coding scheme can be used whenever an efficient coding of power spectral envelopes is required, and can be incorporated as a tool in new stereo source codecs. A particularly interesting application is in HFR systems that are guided by information about the original signal highband envelope. In such a system, the lowband is coded and decoded by means of an arbitrary codec, and the highband is regenerated at the decoder using the decoded lowband signal and the transmitted highband envelope information [PCT WO 98/57436]. Furthermore, the possibility to build a scalable HFR-based stereo codec is offered, by locking the envelope coding to level/balance operation. Hereby the level values are fed into the primary bitstream, which, depending on the implementation, typically decodes to a mono signal. The balance values are fed into the secondary bitstream, which in addition to the primary bitstream is available-to receivers close to the transmitter, taking an IBOC (In-Band On-Channel) digital AM-broadcasting system as an example. When the two bitstreams are combined, the decoder produces a stereo output signal. In addition to the level values, the primary bitstream can contain stereo parameters, e.g. a width parameter. Thus, decoding of this bitstream alone already yields a stereo output, which is improved when both bit streams are available.
The present invention will now be described by way of illustrative examples, not limiting the scope or spirit of the invention, with reference to the accompanying drawings, in which:
a is a block schematic of a parametric stereo decoder module,
b is a block schematic of a pseudo-stereo generator with control parameter inputs,
c is a block schematic of a balance adjuster with control parameter inputs,
a is a block schematic of the encoder side of a scalable HFR-based stereo codec, employing level/balance-coding of the spectral envelope,
b is a block schematic of the corresponding decoder side.
The below-described embodiments are merely illustrative for the principles of the present invention. It is understood that modifications and variations of the arrangements and the details described herein will be apparent to others skilled in the art. It is the intent therefore, to be limited only by the scope of the impending patent claims, and not by the specific details presented by way of description and explanation of the embodiments herein. For the sake of clarity, all below examples assume two channel systems, but apparent to others skilled in the art, the methods can be applied to multichannel systems, such as a 5.1 system.
One method of parameterization of stereo properties according to the present invention, is to determine the original signal stereo-width at the encoder side. A first approximation of the stereo-width is the difference signal, D=L−R, since, roughly put, a high degree of similarity between L and R computes to a small value of D, and vice versa. A special case is dual mono, where L=R and thus D=O. Thus, even this simple algorithm is capable of detecting the type of mono input signal commonly associated with news broadcasts, in which case pseudo-stereo is not desired. However, a mono signal that is fed to Land R at different levels does not yield a zero D signal, even though the perceived width is zero. Thus, in practice more elaborate detectors might be required, employing for example cross-correlation methods. One should make sure that the value describing the left-right difference or correlation in some way is normalized with the total signal level, in order to achieve a level independent detector. A problem with the aforementioned detector is the case when mono speech is mixed with a much weaker stereo signal e.g. stereo noise or background music during speech-to-music/music-to-speech transitions. At the speech pauses the detector will then indicate a wide stereo signal. This is solved by normalizing the stereo-width value with a signal containing information of previous total energy level e.g., a peak decay signal of the total energy. Furthermore, to prevent the stereo-width detector from being trigged by high frequency noise or channel different high frequency distortion, the detector signals should be pre-filtered by a lowpass filter, typically with a cutoff frequency somewhere above a voice's second format, and optionally also by a high-pass filter to avoid unbalanced signal-offsets or hum. Regardless of detector type, the calculated stereo-width is mapped to a finite set of values, covering the entire range, from mono to wide stereo.
a gives an example of the contents of the parametric stereo decoder introduced in
Any prior art pseudo-stereo generator can be used for the width block, such as those mentioned in the background section, or a Schroeder-type early reflection simulating unit (multitap delay) or reverberator.
An alternative method of detecting stereo-properties according to the invention, is described as follows. Again, let L and R denote the left and right input signals. The corresponding signal powers are then given by PL˜L2 and PR˜R2. Now, a measure of the stereo-balance can be calculated as the quotient of the two signal powers, or more specifically as B=(PL+e)/(PR+e), where e is an arbitrary, very small number, which eliminates division by zero. The balance parameter, B, can be expressed in dB given by the relation BdB=10 log10(B). As an example, the three cases PL=10 PR, PL=PR, and PL=O.1PR correspond to balance values of +10 dB, 0 dB, and −10 dB respectively. Clearly, those values map to the locations “left”, “center”, and “right”. Experiments have shown that the span of the balance parameter can be limited to for example +/−40 dB, since those extreme values are already perceived as if the sound originates entirely from one of the two loudspeakers or headphone drivers. This limitation reduces the signal space to cover in the transmission, thus offering bitrate reduction. Furthermore, a progressive quantization scheme can be used, whereby smaller quantization steps are used around zero, and larger steps towards the outer limits, which further reduces the bitrate. Often the balance is constant over time for extended passages. Thus, a last step to significantly reduce the number of average bits needed can be taken: After transmission of an initial balance value, only the differences between consecutive balance values are transmitted, whereby entropy coding is employed. Very commonly, this difference is zero, which thus is signaled by the shortest possible codeword. Clearly, in applications where bit errors are possible, this delta coding must be reset at an appropriate time interval, in order to eliminate uncontrolled error propagation.
The most rudimental decoder usage of the balance parameter, is simply to offset the mono signal towards either of the two reproduction channels, by feeding the mono signal to both outputs and adjusting the gains correspondingly, as illustrated in
The methods described so far, are intended for very low bitrate applications. In applications where higher bitrates are available, it is possible to use more elaborate versions of the above width and balance methods. Stereo-width detection can be made in several frequency bands, resulting in individual stereowidth values for each frequency band. Similarly, balance calculation can operate in a multiband fashion, which is equivalent to applying different filter-curves to two channels that are fed by a mono signal.
The parametric balance coding method can, especially for lower frequency bands, give a somewhat unstable behavior, due to lack of frequency resolution, or due to too many sound events occurring in one frequency band at the same time but at different balance positions. Those balance-glitches are usually characterized by a deviant balance value during just a short period of time, typically one or a few consecutive values calculated, dependent on the update rate. In order to avoid disturbing balance-glitches, a stabilization process can be applied on the balance data. This process may use a number of balance values before and after current time position, to calculate the median value of those. The median value can subsequently be used as a limiter value for the current balance value i.e., the current balance value should not be allowed to go beyond the median value. The current value is then limited by the range between the last value and the median value. Optionally, the current balance value can be allowed to pass the limited values by a certain overshoot factor. Furthermore, the overshoot factor, as well as the number of balance values used for calculating the median, should be seen as frequency dependent properties and hence be individual for each frequency band.
At low update ratios of the balance information, the lack of time resolution can cause failure in synchronization between motions of the stereo image and the actual sound events. To improve this behavior in terms of synchronization, an interpolation scheme based on identifying sound events can be used. Interpolation here refers to interpolations between two, in time consecutive balance values. By studying the mono signal at the receiver side, information about beginnings and ends of different sound events can be obtained. One way is to detect a sudden increase or decrease of signal energy in a particular frequency band. The interpolation should after guidance from that energy envelope in time make sure that the changes in balance position should be performed preferably during time segments containing little signal energy. Since human ear is more sensitive to entries than trailing parts of a sound, the interpolation scheme benefits from finding the beginning of a sound by e.g., applying peak-hold to the energy and then let the balance value increments be a function of the peak-holded energy, where a small energy value gives a large increment and vice versa. For time segments containing uniformly distributed energy in time i.e., as for some stationary signals, this interpolation method equals linear interpolation between the two balance values. If the balance values are quotients of left and right energies, logarithmic balance values are preferred, for left-right symmetry reasons. Another advantage of applying the whole interpolation algorithm in the logarithmic domain is the human ear's tendency of relating levels to a logarithmic scale.
Also, for low update ratios of the stereo-width gain values, interpolation can be applied to the same. A simple way is to interpolate linearly between two in time consecutive stereo-width values. More stable behavior of the stereo-width can be achieved by smoothing the stereo-width gain values over a longer time segment containing several stereo-width parameters. By utilizing smoothing with different attack and release time constants, a system well suited for program material containing mixed or interleaved speech and music is achieved. An appropriate design of such smoothing filter is made using a short attack time constant, to get a short rise-time and hence an immediate response to music entries in stereo, and a long release time, to get a long fall-time. To be able to fast switch from a wide stereo mode to mono, which can be desirable for sudden speech entries, there is a possibility to bypass or reset the smoothing filter by signaling this event. Furthermore, attack time constants, release time constants and other smoothing filter characteristics can also be signaled by an encoder.
For signals containing masked distortion from a psycho-acoustical codec, one common problem with introducing stereo information based on the coded mono signal is an unmasking effect of the distortion. This phenomenon usually referred as “stereo-unmasking” is the result of non-centered sounds that do not fulfill the masking criterion. The problem with stereo-unmasking might be solved or partly solved by, at the decoder side, introducing a detector aimed for such situations. Known technologies for measuring signal to mask ratios can be used to detect potential stereo-unmasking. Once detected, it can be explicitly signaled or the stereo parameters can just simply be decreased.
At the encoder side, one option, as taught by the invention, is to employ a Hilbert transformer to the input signal, i.e. a 90 degree phase shift between the two channels is introduced. When subsequently forming the mono signal by addition of the two signals, a better balance between a center-panned mono signal and “true” stereo signals is achieved, since the Hilbert transformation introduces a 3 dB attenuation for center information. In practice, this improves mono coding of e.g. contemporary pop music, where for instance the lead vocals and the bass guitar commonly is recorded using a single mono source.
The multiband balance-parameter method is not limited to the type of application described in
One particularly interesting application of the above envelope coding method is coding of highband spectral envelopes for HFR-based codecs. In this case no highband residual signal is transmitted. Instead this residual is derived from the lowband. Thus, there is no strict relation between residual and envelope representation, and envelope quantization is more crucial. In order to study the effects of quantization, let Pq and Bq denote the quantized values of P and B respectively. Pq and Bq are then inserted into the above relations, and the sum is formed: PL q+PR q=BqPq/(Bq+1)+Pq/(Bq+1)=Pq(Bq+1)/(Bq+1)=Pq. The interesting feature here is that Bq is eliminated, and the error in total power is solely determined by the quantization error in P. This implies that even though B is heavily quantized, the perceived level is correct, assuming that sufficient precision in the quantization of P is used. In other words, distortion in B maps to distortion in space, rather than in level. As long as the sound sources are stationary in the space over time, this distortion in the stereo perspective is also stationary, and hard to notice. As already stated, the quantization of the stereo-balance can also be coarser towards the outer extremes, since a given error in dB corresponds to a smaller error in perceived angle when the angle to the centerline is large, due to properties of human hearing.
When quantizing frequency dependent data e.g., multi band stereo-width gain values or multi band balance values, resolution and range of the quantization method can advantageously be selected to match the properties of a perceptual scale. If such scale is made frequency dependent, different quantization methods, or so called quantization classes, can be chosen for the different frequency bands. The encoded parameter values representing the different frequency bands, should then in some cases, even if having identical values, be interpreted in different ways i.e., be decoded into different values.
Analogous to a switched L/R- to S/D-coding scheme, the P and B signals may be adaptively substituted by the PL and PR signals, in order to better cope with extreme signals. As taught by [PCf/SEOOI00158], delta coding of envelope samples can be switched from delta-in-time to delta-in-frequency, depending on what direction is most efficient in terms of number of bits at a particular moment. The balance parameter can also take advantage of this scheme: Consider for example a source that moves in stereo field over time. Clearly, this corresponds to a successive change of balance values over time, which depending on the speed of the source versus the update rate of the parameters, may correspond to large delta-in-time values, corresponding to large codewords when employing entropy coding. However, assuming that the source has uniform sound radiation versus frequency, the delta-in-frequency values of the balance parameter are zero at every point in time, again corresponding to small codewords. Thus, a lower bitrate is achieved in this case, when using the frequency delta coding direction. Another example is a source that is stationary in the room, but has a non-uniform radiation. Now the delta-in-frequency values are large, and delta-in-time is the preferred choice.
The P/B-coding scheme offers the possibility to build a scalable HFR-codec, see
For the lowband codec different possibilities exist: It may constantly operate in S/D-mode, and the S and D signals be sent to primary and secondary bitstreams respectively. In this case, a decoding of the primary bitstream results in a full band mono signal. Of course, this mono signal can be enhanced by parametric stereo methods according to the invention, in which case the stereo-parameter(s) also must be located in the primary bitstream. Another possibility is to feed a stereo coded lowband signal to the primary bitstream, optionally together with highband width- and balance-parameters. Now decoding of the primary bitstream results in true stereo for the lowband, and very realistic pseudo-stereo for the highband, since the stereo properties of the lowband are reflected in the high frequency reconstruction. Stated in another way: Even though the available highband envelope representation or spectral coarse structure is in mono, the synthesized highband residual or spectral fine structure is not. In this type of implementation, the secondary bitstream may contain more lowband information, which when combined with that of the primary bitstream, yields a higher quality lowband reproduction. The topology of
The bitstreams are transmitted or stored, and either only 419 or both 419 and 417 are fed to the decoder,
Stated in other words, a method for coding of stereo properties of an input signal, includes at an encoder, the step of calculating a width-parameter that signals a stereo-width of said input signal, and at a decoder, a step of generating a stereo output signal, using said width-parameter to control a stereo-width of said output signal. The method further comprises at said encoder, forming a mono signal from said input signal, wherein, at said decoder, said generation implies a pseudo-stereo method operating on said mono signal. The method further implies splitting of said mono signal into two signals as well as addition of delayed version(s) of said mono signal to said two signals, at level(s) controlled by said width-parameter. The method further includes that said delayed version(s) are high-pass filtered and progressively attenuated at higher frequencies prior to being added to said two signals. The method further includes that said width-parameter is a vector, and the elements of said vector correspond to separate frequency bands. The method further includes that if said input signal is of type dual mono, said output signal is also of type dual mono.
A method for coding of stereo properties of an input signal, includes at an encoder, calculating a balance parameter that signals a stereo-balance of said input signal, and at a decoder, generate a stereo output signal, using said balance-parameter to control a stereo-balance of said output signal.
In this method, at said encoder, a mono signal from said input signal is formed, and at said decoder, said generation implies splitting of said mono signal into two signals, and said control implies adjustment of levels of said two signals. The method further includes that a power for each channel of said input signal is calculated, and said balance-parameter is calculated from a quotient between said powers. The method further includes that said powers and said balance-parameter are vectors where every element corresponds to a specific frequency band. The method further includes that at said decoder it is interpolated between two in time consecutive values of said balance-parameters in a way that the momentary value of the corresponding power of said mono signal controls how steep the momentary interpolation should be. The method further includes that said interpolation method is performed on balance values represented as logarithmic values. The method further includes that said values of balance parameters are limited to a range between a previous balance value, and a balance value extracted from other balance values by a median filter or other filter process, where said range can be further extended by moving the borders of said range by a certain factor. The method further includes that said method of extracting limiting borders for balance values, is, for a multi band system, frequency dependent. The method further includes that an additional level-parameter is calculated as a vector sum of said powers and sent to said decoder, thereby providing said decoder a representation of a spectral envelope of said input signal. The method further includes that said level-parameter and said balance-parameter adaptively are replaced by said powers. The method further includes that said spectral envelope is used to control a HFR-process in a decoder. The method further includes that said level-parameter is fed into a primary bitstream of a scalable HFR-based stereo codec, and said balance-parameter is fed into a secondary bitstream of said codec. Said mono signal and said width-parameter are fed into said primary bitstream. Furthermore, said width-parameters are processed by a function that gives smaller values for a balance value that corresponds to a balance position further from the center position. The method further includes that a quantization of said balance-parameter employs smaller quantization steps around a center position and larger steps towards outer positions. The method further includes that said width-parameters and said balance-parameters are quantized using a quantization method in terms of resolution and range which, for a multiband system, is frequency dependent. The method further includes that said balance parameter adaptively is delta-coded either in time or in frequency. The method further includes that said input signal is passed though a Hilbert transformer prior to forming said mono signal.
An apparatus for parametric stereo coding, includes, at an encoder, means for calculation of a width-parameter that signals a stereo-width of an input signal, and means for forming a mono signal from said input signal, and, at a decoder, means for generating a stereo output signal from said mono signal, using said width-parameter to control a stereo-width of said output signal.
Number | Date | Country | Kind |
---|---|---|---|
0102481 | Jul 2001 | SE | national |
0200796 | Mar 2002 | SE | national |
0202159 | Jul 2002 | SE | national |
PCT/SE02/01372 | Jul 2002 | WO | international |
This application is a divisional of U.S. patent application Ser. No. 11/238,982 filed on Sep. 28, 2005, which is a divisional of U.S. patent application Ser. No. 10/483,453 filed on Jan. 8, 2004, now U.S. Pat. No. 7,382,886, each of which are hereby incorporated in their entireties by reference.
Number | Name | Date | Kind |
---|---|---|---|
4053711 | DeFreitas et al. | Oct 1977 | A |
4166924 | Berkley et al. | Sep 1979 | A |
4706287 | Blackmer et al. | Nov 1987 | A |
5463424 | Dressler | Oct 1995 | A |
5559891 | Kuusama et al. | Sep 1996 | A |
5613035 | Kim | Mar 1997 | A |
5671287 | Gerzon | Sep 1997 | A |
5701346 | Herre et al. | Dec 1997 | A |
5862228 | Davis | Jan 1999 | A |
5873065 | Akagiri et al. | Feb 1999 | A |
5883962 | Hawks | Mar 1999 | A |
5890108 | Yeldener | Mar 1999 | A |
5890125 | Davis et al. | Mar 1999 | A |
6507658 | Abel et al. | Jan 2003 | B1 |
6771777 | Gbur et al. | Aug 2004 | B1 |
20020037086 | Irwan et al. | Mar 2002 | A1 |
Number | Date | Country |
---|---|---|
19947098 | Nov 2000 | DE |
0478096 | Jan 1987 | EP |
0273567 | Jul 1988 | EP |
0858067 | Aug 1998 | EP |
0858067 | Aug 1998 | EP |
0989543 | Mar 2000 | EP |
1107232 | Jun 2001 | EP |
2100430 | Dec 1982 | GB |
H02-012299 | Jan 1990 | JP |
02-177782 | Jul 1990 | JP |
03-214956 | Sep 1991 | JP |
H04-301688 | Oct 1992 | JP |
H05-165500 | Jul 1993 | JP |
H06-090209 | Mar 1994 | JP |
H06-202629 | Jul 1994 | JP |
H06-215482 | Aug 1994 | JP |
H08-254994 | Oct 1996 | JP |
H08-305398 | Nov 1996 | JP |
H09-505193 | May 1997 | JP |
H09-261064 | Oct 1997 | JP |
H11-262100 | Sep 1999 | JP |
11-317672 | Nov 1999 | JP |
2000-083014 | Mar 2000 | JP |
2000-505266 | Apr 2000 | JP |
2001-184090 | Jul 2001 | JP |
2004-535145 | Nov 2004 | JP |
96-003455 | Mar 1996 | KR |
96-0012475 | Sep 1996 | KR |
WO 9504442 | Feb 1995 | WO |
WO 9516333 | Jun 1995 | WO |
WO9857436 | Dec 1998 | WO |
WO 0045378 | Aug 2000 | WO |
WO 03007656 | Jan 2003 | WO |
Number | Date | Country | |
---|---|---|---|
20090316914 A1 | Dec 2009 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 11238982 | Sep 2005 | US |
Child | 12496926 | US | |
Parent | 10483453 | Jan 2004 | US |
Child | 11238982 | US |