Digital audio transmission.
Abstract
@ A digital sub-band coder for audio signals provides for adaptive allocation of bits to (eg) two sub-bands. Estimation of the relative energy content of the bands is determined by reference only to those bits which are always transmitted, so that decoding can take place without the need for transmission of side information. Accompanying data may be transmitted in a fixed bit position, or may be by insertion of the data into the lower energy band. In the latter case, the data may optionally be duplicated in both bands when their energy contents are random. The invention can be applied to embedded ADPCM coders.

Term
Term ended
Projected expiry passed 15 November 2004, 21.9 years ago.
- Priority
- Filed
- Published
- Projected expiry
- Today
19 claims: 10 independent, 9 dependent
- 1A method of coding a digital sub-band coded signal, in which the number of bits allocated to the respective sub-bands is adaptively determined as a function of those bits which are always transmitted, whereby the coded signal may be decoded by reference only to the transmitted bits.
- 5A method as claimed in any preceding claim in which the energy content of two adjacent sub-bands is compared.
- 8A method as claimed in anyone of the preceding claims in which each sub-band is coded by an embedded adaptive differential pulse code modulation (ADPCM) encoder and the comparison of the energy content of the sub-bands is effected by comparing the stepsizes of the ADPCM encoders.
- 9A method as claimed in any preceding claim in which the sub-bands are two adjoining half-bands of 64kbit/s voice transmission system of a 7kHz speech bandwidth, the total transmission capacity is 8 bits, the lower sub-band is encoded by 4 or 5 bits and the upper sub-band by 2 or 3 bits.
- 10A method of transmitting a digital data signal together with a digital sub-band coded signal which comprises comparing the energy contents different sub-bands and adaptively inserting data bits into a sub-band with a low energy content.
- 11A digital sub-band coder comprising filter means for dividing the input signal into two or more frequency sub-bands, and adaptive bit allocation means for adaptively determining the number of bits to be used for encoding each sub-band in dependence on the relative energy content of the sub-bands, characterised in that the bit allocation means is responsive only to those bits of each sub-band which are always transmitted.
- 19Apparatus to transmit digital sub-band encoded signals and digital data signals, comprising means to divide an incoming frequency band into a plurality of sub-bands and means for encoding each sub-band separately, characterised in that the apparatus further comprises means to measure and compare the energy contents of the sub-bands and selectively insert data bits into a sub-band having a low energy envelope.
Independent claims10
43 paragraphs, as filed
This invention relates to transmission using sub-band coding, and particularly, though not exclusively, to multiplexing data and digitally coded audio for transmission.
Sub-band coding (SBC) is a coding technique in which a parent frequency band, usually speech, is divided into two or more sub-bands, the frequency range within each of which is separately encoded and transmitted. The main advantages of sub-band coding, for example, over pulse code modulation (PCM) or adaptive differential pulse code modulation (ADPCM) are generally considered to be that: <ul id="ul0001" list-style="none"><li>a) the quantisation noise generated within each sub-band remains confined to that frequency range and on decoding does not mask low level speech sounds in other bands; and</li><li>b) each sub-band signal may be coded according to the perceptual contribution that sub-band makes to the overall subjective quality.</li></ul>
According to one aspect of the present invention, there is provided a method of coding a digital sub-band coded signal, in which the number of bits allocated to the respective sub-bands is adaptively determined as a function of those bits which are always transmitted, whereby the coded signal may be decoded by reference only to the transmitted bits.
In another aspect \'!e provide a digital sub-band coder comprising filter means for dividing the input signal into two or more frequency sub-bands, and adaptive bit allocation means for adaptively determining the number of bits to be used for encoding each sub-band in dependence on the relative energy content of the sub-bands, characterised in that the bit allocation means is responsive only to those bits of each sub-band which are always transmitted.
A further aspect of the invention provides a method of transmitting a digital data signal together with a digital sub-band coded signal which comprises comparing the energy contents different sub-bands and adaptively inserting data bits into a sub-band with a low energy content.
In yet another aspect, we provide apparatus to transmit digital sub-band encoded signals and digital data signals, comprising means to divide an incoming frequency band into a plurality of sub-bands and means for encoding each sub-band separately, characterised in that the apparatus further comprises means to measure and compare the energy contents of the sub-bands and selectively insert data bits into a sub-band having a low energy envelope.
Conveniently the comparison is made between two adjoining sub-bands. The sub-bands may themselves be sub-bands of lower order sub-bands, or may be the only sub-bands into which the parent freauency band has been divided.
The comparison may be made, for example, by subtracting the respective sub-band stepsizes, but is preferably made by determining their ratio.
The present invention will now be explained further by way of example and with reference to the accompanying drawings of which: <ul id="ul0002" list-style="none"><li>Figures 1, 3 and 5 are schematic diagrams of SBC encoders incorporating the present invention;</li><li>Figure 2, 4 and 6 are schematic diagrams of SBC decoders incorporating the present invention; and</li><li>Figure 7 is a schematic diagram of an embedded ADPCM coder and decoder;</li></ul>
Referring to Figure 1, a sub-band coder for speech signals receives at an input 1, a speech signal which is assumed to be PCM coded at a sampling rate of 16kHz. The input frequency spectrum is divided into upper and lower sub-bands by high-pass and low-pass filters 2, 3, the two bands then being down-sampled to 8kHz in "decimators" 4,5. For the purposes of illustration, the signals at this point are assumed to be quantised to four-bit accuracy, in both bands, although in practice more bits would be used, not necessarily the same number in each band.
An output multiplexer 6 is provided to assemble the output bits into suitable form for onward transmission from its output; it could be part of a t.d.m. system, or simply a parallel-to-serial converter. It has, by way of example, an 8-bit input. One bit of these eight is permanently assigned to transmission of an auxiliary data signal D; note however that this is incidental to the principle being illustrated in Figure 1.
An adaptive allocation of the remaining seven available bits to the two channels is effected on the basis that one channel will be transmitted with 3-bit accuracy, the other 4-bit.· Thus the 3 most significant bits h<sub>3</sub>-h<sub>1</sub>, 1<sub>3</sub>-1<sub>1</sub> of the two channels are permanently assigned to respective bit inputs of the multiplexer 6, whilst the eighth input of the multiplexer receives, via a selector 8, the least significant bit h<sub>0</sub> or 1<sub>0</sub> of a selected one of the two channels.
The corresponding decoder (Figure 2) is essentially a mirror image of the transmitter, a demultiplexer 16 feeding the three most significant inputs of up-samplers 14, 15 which in turn supply interpolation filters 12,13 the outputs of which are added at 10 to form the decoded output 11.
Returning to the coder of Figure 1, the selector 8 is controlled by a decision circuit 9 which determines which of the two least significant bits h<sub>O</sub> or 1<sub>0</sub> is to be transmitted. This decision is made as a function of the relative signal or energy content of the two channels, and the corresponding deselection (by a one-bit demultiplexer 18) at the decoder is controlled by a functionally identical decision circuit 19.
So that the decoder decision circuit can accurately "track" operation of that in the coder, the decision circuits are in each case responsive only to the three most significant bits h<sub>3</sub>-h<sub>1</sub>, 1<sub>3</sub>-1<sub>1</sub> of the two channels. The decoder decision circuit does not, of course, have access to both h<sub>O</sub> and 1<sub>0</sub>; nor does it "know" prior to making its decision, which of them it does have access to.
The decision circuits 9, 19 serve to assess which channel has the higher energy content, so that the dropping of the least significant bit will occur in a channel where it makes the lesser noise contribution. Although it would be possible to do this on a sample-by-sample basis (where a simple comparator would suffice), for a practical system it would compare and average energy estimates formed over a period; this will be discussed further below.
Figure 3 shows a modified version of the coder of Figure 1, where, of the eight available bits, four are allocated to each of the two channels, the auxiliary data D being adaptively inserted in the least significant bit position h<sub>O</sub> or l<sub>0</sub> of one of the two channels. As shown, however, the decision circuit 9' has two outputs x, y which are active when the ratio of the energy in the upper channel to that in the lower respectively exceeds a lower threshold or falls below a higher one, so that when it falls between the two thresholds, the data is inserted into both channels. In the receiver, a single decision threshold is used, as before. In this way, corruption of the data by errors in the decision, due to finite computational accuracy, or transmission errors, are reduced.
A more practical version of the coder and decoder of Figures 2 and 3 will now be described; note however that much of what follows could also be applied to the coder and decoder of Figures 1 and 2.
The coder, suitable for receiving a wide band speech input signal of 7kHz bandwidth sampled at 16kHz and transmitting it at 64kbit/s uses ADPCM (adaptive differential pulse code modulation). Differential coders normally suffer from an excessive signal-to-noise penalty if bits are deleted and to rectify this, an embedded ADPCM coding is used. Embedded coders are described in "Embedded DPCM for variable bit rate transmission", David J Goodman, IEEE Trans. Comm. Vol COM-28, No.7, July 1980, pp1040 to 1046. The encoder described there is not an adaptive coder (nor is it used for sub-band coding). Essentially, an embedded DPCM coder solves the problem (with a small SNR penalty due to less accurate prediction and hence a higher quantiser step size and thus quantisation noise) in that any bits liable to be omitted on transmission are excluded from the predictor loop at both coder and decoder.
The sub-band coder 1 comprises a digital filter 22 consisting of a high pass filter 22a and a low pass filter 22b which divides the digitally (eg PCM) encoded input signal at terminal 3 into two sub-bands, which are down sampled in respective down sampling circuits 24a and 24b of a "decimator" 24 from where they are supplied to ADPCM encoder 5 connected to the output of down sampler 4a, and ADPCM encoder 6 connected to the output of down sampler 4b. The ADPCM encoders 5 and 6 have code word outputs 7 and 8 respectively. The ADPCM encoders 5 and 6 also have outputs 9 and 10 providing a signal indicative of the stepsize to a decision circuit 31. The output of the decision circuit 31 is connected to data insertion circuits 32 and 33, and controls their operation in the manner explained below. The outputs of data insertion circuits 32 and 33 form the inputs of an 8bit multiplexer 34 whose output is the signal to be transmitted.
The digital filter 2 comprises quadrature mirror filters (QMFs) which are used to split the input spectrum into two overlapping bands (0-4 kHz, and 4-8 kHz). Suitable filter constants are listed in Table 1. The signals contained in each of the two sub-bands are then applied to the respective down sampler 24a (4-8 kHz) and 24b (0-4 kHz) where they are down sampled to 8 kilo samples per second and applied to the two independent ADPCM coders 25 and 26 respectively. ADPCM coders per se are known, and for the present purposes it will be sufficient to note that the encoding of the two sub-bands in ADPCM encoders 5 and 6 is not symmetric, in that the lower sub-band is encoded using ADPCM with a 5 bit robust Jayant quantiser (AQJ) and a fixed fourth order predictor optimised for 4 kHz band width speech. The higher sub-band (4-8 kHz) is encoded by a 3 bit robust AQJ with a fixed first order predictor.
In the receiver the lower and upper sub-bands are demultiplexed in demultiplexer 55 and are then processed in essentially the inverse order of that performed on encoding, with the components 42 to 53 being essentially the mirror images of the corresponding components 22 to 33 of the encoder.
The 64 Kbit/s signal is demultiplexed, decoded, interpolated and processed by the receiver QMF bank to recover the 7 kHz bandwidth signal, the two channels being combined in an adder 43.
The SBC encoder and decoder additionally include data insertion and data recovery circuits 32, 33 and 52, 53 respectively, as well as the data insertion and data recovery decision circuits 31 and 51.
The 3 bit and 5 bit outputs of the ADPCM encoders 5 and 6 respectively are applied to the data insertion circuits 12 and 13 where a single data bit is inserted - into, and thus replaces the existing bit in, the fifth encoded (least significant) bit position for the lower channel output, or into the third encoded bit position for the upper channel, or into both, in accordance with the data insertion strategy described below.
Thus the output from the SBC coder 1 is configured either as 5 <sup>+</sup> (2 <sup>+</sup> D) or (4 <sup>+</sup> D) <sup>+</sup> 3, or (4 <sup>+</sup> D) <sup>+</sup> (2 + D), where D is the data bit.
In this embodiment, the decision circuits 31,51 are comparators which compare the adaptor stepsize (scaling factor) generated in the stepsize adapting circuits 46 and 47 of Figure 7, which is a block diagram of the ADPCM coders 25,26 and decoders 45,46. The data insertion decision strategy and the corresponding data recovery strategy operate on the following principles.
It is clearly desirable that the decision as to which channel receives the data bit should not require any transmission of side information to the receiver to provide correct demultiplexing. For this reason a backward adaptive data multiplexing strategy is employed where the multiplexing and demultiplexing data channel selection decision is derived from the transmit quantiser output code words and receiver input code words respectively.
The data channel selection strategy is determined as follows. A spectral envelope estimate, that is essentially a measure for the energy content, is obtained for each sub-band signal and the estimates for the two channels are then compared to give a three way decision voiced (such as vowels), unvoiced (such as the syllable s) and intermediate.
When it is clear that voiced sounds are present the upper channel (4-8 kHz) is selected for the data bit, and 5 bits are then available for accurate representation of the predominantly lower channel signal.
When unvoiced sounds are present, or alternatively, when it is not clear whether the speech is voiced or unvoiced, the lower channel (0-4 kHz) is arranged to carry the data bit thus leaving 4 bits available for coding the speech signal in the lower band. The full 3 bits are then available in the upper channel for accurately coding the higher frequencies of the speech signal.
To avoid the need to transmit any side information the short term spectral envelope information for the data insertion decision is derived from the quantiser stepsize parameters for the two channels. This information is derived in a backward adaptive mode, using the robust Jayant algorithm, ie the stepsize parameter for the lower and upper channel quantisers respectively is given by<maths id="math0001" num=""><img file="EP0145332A2_D0001.tif" /></maths>and<maths id="math0002" num=""><img file="EP0145332A2_D0002.tif" /></maths>where Il(n) is the nth output codeword for the lower channel quantiser and Iu(n) is the output codeword from the upper channel quantiser. Ml(.) and Mu(.) are multiplier functions shown in table 1b for quantisers with various numbers of bits. y is a 'leakage factor' which is required in order to dissipate the effect of transmission errors. Typically, γ is chosen to be 0.984.
In practice it is more convenient in the coder implementation to re-express equations 1 and 2 in logarithmic form to eliminate the complex exponentiation operations:<maths id="math0003" num=""><img file="EP0145332A2_D0003.tif" /></maths>and<maths id="math0004" num=""><img file="EP0145332A2_D0004.tif" /></maths>where, dl(.), du(.), ml(.) and mu(.) are respectively the logarithmic values of Δl(.), Δu(.),. Ml(.) and Mu(.).
It is easy to show 'that the above recursive algorithm produces a measure of the short-term envelope of the bandlimited speech within each sub-band.
The decision as to which sub-band is to receive the data bit is determined as follows. First the scaled ratio of the lower and upper quantiser step-sizes is obtained as:<maths id="math0005" num=""><img file="EP0145332A2_D0005.tif" /></maths>where is a constant and V(n) is the time varying signal applied to the subsequent data bit assignment strategy.
To simplify the digital processing, equation 5 is expressed in a logarithmic form:<maths id="math0006" num=""><img file="EP0145332A2_D0006.tif" /></maths>
By experiment the value of in equation 7 was set to 2.
It is apparent that a voiced/unvoiced decision could be made if the decision threshold is set to v(n) = 0. Unfortunately the use of such a hard decision process creates the possibility of systematic data channel selection errors occurring when the value of v(n) lies close to the decision level. This may happen due to the effects of finite precision limitations or as a result of the residual effect of channel errors and would cause the v(n) calculated from the decoded signals in the receiver not being perfectly in synchronism with the function v(n) determined at the transmitter.
For this reason the transmitter implies (soft) decision thresholds, which cause the insertion of the data bit simultaneously into both sub-bands when v(n) falls between an upper and a lower threshold, u and 1. The decoder forms the comparison with a hard middle threshold, . Although the use of both channels momentarily reduces the performance on the speech channel, it makes the data channel more robust to transmission errors and eliminates systematic errors due to misalignment of the stepsize parameters derived in the encoder and the decoder respectively.
In a pilot experiment, the adaptive data insertion technique of the present invention was simulated in floating point Fortran on a mini computer and compared with a fixed data bit assignment technique, that is the technique where the data are multiplexed at all times into either the upper or the lower sub-band code words. The finite precision limitations of an NEC µP0 7720 microprocessor such as is commonly used in speech coding transmitters and receivers, were simulated and facilities were also provided to inject random errors into the serialised output (the channel) of the sub-band coder.
The afore-described is a simple yet potentially powerful technique for insertion, for example, 8 Kbit/s of data into a 56 Kbit/s SBC coded wideband speech channel. The technique described may be combined with for example adaptive prediction in the lower band, to increase the data carrying capacity beyond this.
The main advantages of the technique are: <ul id="ul0003" list-style="none"><li>1. The scheme is relatively simple to implement and uses processing capacity available in the coder DSP (digital signal processing) micro computer.</li><li>2. There is no need for transmitting side information to control the data demultiplexing function.</li><li>3. In the example described above the scheme offers a constant data carrying capacity of 8 Kbit/s.</li><li>4. The technique may be combined with other techniques to increase the data carrying capacity, in the example 2 16 Kbit/s or more.</li></ul>
Compared to a fixed data multiplexing strategy, where the data bit is assigned constantly to eg the lower channel, the present technique results in a noticeable reduction in the perceivable distortion. <tables id="tabl0001" num="0001"><img file="EP0145332A2_D0007.tif" /></tables>
12 sheets
Sheet 1 Sheet 2 Sheet 3 Sheet 4 Sheet 5 Sheet 6 Sheet 7 Sheet 8 Sheet 9 Sheet 10 Sheet 11 Sheet 12
Every citation, both ways
| Document | Relation | Office | Cited during |
|---|---|---|---|
| US4803727A | Cited by | United States of America | Search report |
| DE3703143A1 | Cited by | Germany | Search report |
| US8233924B2 | Cited by | United States of America | Applicant |
| EP0176243A3 | Cited by | European Patent Office (EPO) | Search report |
| US4949383A | Cited by | United States of America | Search report |
| EP0273567A1 | Cited by | European Patent Office (EPO) | Search report |
| EP0273567A1 | Cited by | European Patent Office (EPO) | Search report |
| EP0176243A2 | Cited by | European Patent Office (EPO) | Search report |
| US4942607A | Cited by | United States of America | Search report |
| DE3721478A1 | Cited by | Germany | Search report |
| DE3642982A1 | Cited by | Germany | Search report |
| EP0372601A1 | Cited by | European Patent Office (EPO) | Search report |
| EP0070949A1 | Cites | European Patent Office (EPO) | Search report |
| GB2063018A | Cites | United Kingdom | Search report |
9 members in 6 offices
Priority claims4
| Document | Office | Kind | Date |
|---|---|---|---|
| 8330885 | United Kingdom | A | |
| 8330885 | United Kingdom | – | |
| 8330885 | – | – | – |
| GB19830030885 | – | – | – |
Members9
| Document | Office | Kind | |
|---|---|---|---|
| GB8330885D0 | United Kingdom | D0 | |
| EP0145332A2This record | European Patent Office (EPO) | A2 | |
| EP0145332A3 | European Patent Office (EPO) | A3 | |
| JPS60169249A | Japan | A | |
| US4703480A | United States of America | A | |
| EP0145332B1 | European Patent Office (EPO) | B1 | |
| AT35883T | Austria | T | |
| DE3472873D1 | Germany | D1 | |
| JPH0243382B2 | Japan | B2 |
51 legal events, as 4 offices reported them to INPADOC
Over the term
Point at a mark for the eventEvents
| Event | Code | Office | |
|---|---|---|---|
| Appeal reference modifiedAppealORIGINAL CODE: EPIDOSCREFNOAPAH | APAH | EP | |
| Notification of lapseLapsedST | ST | FR | |
| Nl: lapsed or anulled due to non-payment of the annual feeLapsedNLV4 | NLV4 | EP | |
| Lapsed in a contracting state [announced via postgrant information from national office to epo]LapsedPG25 | PG25 | EP | |
| Gb: european patent ceased through non-payment of renewal feeCeasedGBPC | GBPC | EP | |
| Lapsed in a contracting state [announced via postgrant information from national office to epo]LapsedPG25 | PG25 | EP | |
| Lapsed in a contracting state [announced via postgrant information from national office to epo]LapsedPG25 | PG25 | EP | |
| Lapsed in a contracting state [announced via postgrant information from national office to epo]LapsedPG25 | PG25 | EP | |
| European patent in force as of 2002-01-01IF02 | IF02 | GB | |
| Annual fee paid to national office [announced via postgrant information from national office to epo]GrantedPGFP | PGFP | EP | |
| Annual fee paid to national office [announced via postgrant information from national office to epo]GrantedPGFP | PGFP | EP | |
| Annual fee paid to national office [announced via postgrant information from national office to epo]GrantedPGFP | PGFP | EP | |
| Annual fee paid to national office [announced via postgrant information from national office to epo]GrantedPGFP | PGFP | EP | |
| Patent ceasedCeasedPL | PL | CH | |
| Be: lapsedLapsedBERE | BERE | EP | |
| Lapsed in a contracting state [announced via postgrant information from national office to epo]LapsedPG25 | PG25 | EP | |
| Lapsed in a contracting state [announced via postgrant information from national office to epo]LapsedPG25 | PG25 | EP | |
| Lapsed in a contracting state [announced via postgrant information from national office to epo]LapsedPG25 | PG25 | EP | |
| Annual fee paid to national office [announced via postgrant information from national office to epo]GrantedPGFP | PGFP | EP | |
| Annual fee paid to national office [announced via postgrant information from national office to epo]GrantedPGFP | PGFP | EP | |
| Se: european patent has lapsedLapsedEUG | EUG | EP | |
| Lapsed in a contracting state [announced via postgrant information from national office to epo]LapsedPG25 | PG25 | EP | |
| Lapsed in a contracting state [announced via postgrant information from national office to epo]LapsedPG25 | PG25 | EP | |
| Lapsed in a contracting state [announced via postgrant information from national office to epo]LapsedPG25 | PG25 | EP | |
| Annual fee paid to national office [announced via postgrant information from national office to epo]GrantedPGFP | PGFP | EP | |
| Se: european patent in force in swedenEAL | EAL | EP | |
| Annual fee paid to national office [announced via postgrant information from national office to epo]GrantedPGFP | PGFP | EP | |
| Annual fee paid to national office [announced via postgrant information from national office to epo]GrantedPGFP | PGFP | EP | |
| Lu: last paid annual feeEPTA | EPTA | EP | |
| It: last paid annual feeITTA | ITTA | EP | |
| Nl: decision of oppositionOppositionNLR2 | NLR2 | EP | |
| Opposition rejectedOpposition27O | 27O | EP | |
| Opposition rejectedOppositionORIGINAL CODE: 0009273PLBN | PLBN | EP | |
| Information on the status of an ep patent application or granted ep patentGrantedSTATUS: OPPOSITION REJECTEDSTAA | STAA | EP | |
| Nl: opposition has been filed with the epoOppositionNLR1 | NLR1 | EP | |
| Opposition filedOpposition26 | 26 | EP | |
| Opposition filedOppositionORIGINAL CODE: 0009260PLBI | PLBI | EP | |
| Fr: translation filedET | ET | EP | |
| Corresponds to:REF | REF | EP | |
| It: translation for a ep patent filedITF | ITF | EP | |
| It: translation for a ep patent filedITF | ITF | EP | |
| Designated contracting statesAK | AK | EP | |
| Corresponds to:REF | REF | EP | |
| (expected) grantORIGINAL CODE: 0009210GRAA | GRAA | EP | |
| Party data changed (applicant data changed or rights of an application transferred)RAP1 | RAP1 | EP | |
| First examination report despatched17Q | 17Q | EP | |
| Request for examination filed17P | 17P | EP | |
| Designated contracting statesAK | AK | EP | |
| Designated contracting statesAK | AK | EP | |
| Search report despatchedORIGINAL CODE: 0009013PUAL | PUAL | EP | |
| Public reference made under article 153(3) epc to a published international application that has entered the european phaseORIGINAL CODE: 0009012PUAI | PUAI | EP |
Numbers
- Publication
- 0145332
- Publication, DOCDB
- 0145332
- Publication, EPODOC
- EP0145332
- Application
- 84307933
- Application, DOCDB
- 84307933
- Application, EPODOC
- EP19840307933
Titles3
- German
- Digitale Übertragung von Audiosignalen
- English
- Digital audio transmission
- French
- Transmission numérique de signaux audio
Classification
- CPC, 2
- H04M11/064
- H04B1/667
- IPC, 4
- H04B14 04
- H03M7 30
- H04B1 66
- H04M11 06
Designated states11
- Contracting states, 11
- Austria
- Belgium
- Switzerland
- Germany
- France
- United Kingdom
- Italy
- Liechtenstein
- Luxembourg
- Netherlands (Kingdom of the)
- Sweden