xref: /aosp_15_r20/external/liblc3/test/encoder.py (revision 49fe348c0058011ee60b6957cdd9d52742df84bc)
1*49fe348cSAndroid Build Coastguard Worker#
2*49fe348cSAndroid Build Coastguard Worker# Copyright 2022 Google LLC
3*49fe348cSAndroid Build Coastguard Worker#
4*49fe348cSAndroid Build Coastguard Worker# Licensed under the Apache License, Version 2.0 (the "License");
5*49fe348cSAndroid Build Coastguard Worker# you may not use this file except in compliance with the License.
6*49fe348cSAndroid Build Coastguard Worker# You may obtain a copy of the License at
7*49fe348cSAndroid Build Coastguard Worker#
8*49fe348cSAndroid Build Coastguard Worker#     http://www.apache.org/licenses/LICENSE-2.0
9*49fe348cSAndroid Build Coastguard Worker#
10*49fe348cSAndroid Build Coastguard Worker# Unless required by applicable law or agreed to in writing, software
11*49fe348cSAndroid Build Coastguard Worker# distributed under the License is distributed on an "AS IS" BASIS,
12*49fe348cSAndroid Build Coastguard Worker# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13*49fe348cSAndroid Build Coastguard Worker# See the License for the specific language governing permissions and
14*49fe348cSAndroid Build Coastguard Worker# limitations under the License.
15*49fe348cSAndroid Build Coastguard Worker#
16*49fe348cSAndroid Build Coastguard Worker
17*49fe348cSAndroid Build Coastguard Workerimport numpy as np
18*49fe348cSAndroid Build Coastguard Workerimport scipy.signal as signal
19*49fe348cSAndroid Build Coastguard Workerimport scipy.io.wavfile as wavfile
20*49fe348cSAndroid Build Coastguard Workerimport struct
21*49fe348cSAndroid Build Coastguard Workerimport argparse
22*49fe348cSAndroid Build Coastguard Worker
23*49fe348cSAndroid Build Coastguard Workerimport lc3
24*49fe348cSAndroid Build Coastguard Workerimport tables as T, appendix_c as C
25*49fe348cSAndroid Build Coastguard Worker
26*49fe348cSAndroid Build Coastguard Workerimport attdet, ltpf
27*49fe348cSAndroid Build Coastguard Workerimport mdct, energy, bwdet, sns, tns, spec
28*49fe348cSAndroid Build Coastguard Workerimport bitstream
29*49fe348cSAndroid Build Coastguard Worker
30*49fe348cSAndroid Build Coastguard Workerclass Encoder:
31*49fe348cSAndroid Build Coastguard Worker
32*49fe348cSAndroid Build Coastguard Worker    def __init__(self, dt_ms, sr_hz):
33*49fe348cSAndroid Build Coastguard Worker
34*49fe348cSAndroid Build Coastguard Worker        dt = { 7.5: T.DT_7M5, 10: T.DT_10M }[dt_ms]
35*49fe348cSAndroid Build Coastguard Worker
36*49fe348cSAndroid Build Coastguard Worker        sr = {  8000: T.SRATE_8K , 16000: T.SRATE_16K, 24000: T.SRATE_24K,
37*49fe348cSAndroid Build Coastguard Worker               32000: T.SRATE_32K, 48000: T.SRATE_48K }[sr_hz]
38*49fe348cSAndroid Build Coastguard Worker
39*49fe348cSAndroid Build Coastguard Worker        self.ne = T.NE[dt][sr]
40*49fe348cSAndroid Build Coastguard Worker
41*49fe348cSAndroid Build Coastguard Worker        self.attdet = attdet.AttackDetector(dt, sr)
42*49fe348cSAndroid Build Coastguard Worker        self.ltpf = ltpf.LtpfAnalysis(dt, sr)
43*49fe348cSAndroid Build Coastguard Worker
44*49fe348cSAndroid Build Coastguard Worker        self.mdct = mdct.MdctForward(dt, sr)
45*49fe348cSAndroid Build Coastguard Worker        self.energy = energy.EnergyBand(dt, sr)
46*49fe348cSAndroid Build Coastguard Worker        self.bwdet = bwdet.BandwidthDetector(dt, sr)
47*49fe348cSAndroid Build Coastguard Worker        self.sns = sns.SnsAnalysis(dt, sr)
48*49fe348cSAndroid Build Coastguard Worker        self.tns = tns.TnsAnalysis(dt)
49*49fe348cSAndroid Build Coastguard Worker        self.spec = spec.SpectrumAnalysis(dt, sr)
50*49fe348cSAndroid Build Coastguard Worker
51*49fe348cSAndroid Build Coastguard Worker    def analyse(self, x, nbytes):
52*49fe348cSAndroid Build Coastguard Worker
53*49fe348cSAndroid Build Coastguard Worker        att = self.attdet.run(nbytes, x)
54*49fe348cSAndroid Build Coastguard Worker
55*49fe348cSAndroid Build Coastguard Worker        pitch_present = self.ltpf.run(x)
56*49fe348cSAndroid Build Coastguard Worker
57*49fe348cSAndroid Build Coastguard Worker        x = self.mdct.run(x)[:self.ne]
58*49fe348cSAndroid Build Coastguard Worker
59*49fe348cSAndroid Build Coastguard Worker        (e, nn_flag) = self.energy.compute(x)
60*49fe348cSAndroid Build Coastguard Worker        if nn_flag:
61*49fe348cSAndroid Build Coastguard Worker            self.ltpf.disable()
62*49fe348cSAndroid Build Coastguard Worker
63*49fe348cSAndroid Build Coastguard Worker        bw = self.bwdet.run(e)
64*49fe348cSAndroid Build Coastguard Worker
65*49fe348cSAndroid Build Coastguard Worker        x = self.sns.run(e, att, x)
66*49fe348cSAndroid Build Coastguard Worker
67*49fe348cSAndroid Build Coastguard Worker        x = self.tns.run(x, bw, nn_flag, nbytes)
68*49fe348cSAndroid Build Coastguard Worker
69*49fe348cSAndroid Build Coastguard Worker        (xq, lastnz, x) = self.spec.run(bw, nbytes,
70*49fe348cSAndroid Build Coastguard Worker            self.bwdet.get_nbits(), self.ltpf.get_nbits(),
71*49fe348cSAndroid Build Coastguard Worker            self.sns.get_nbits(), self.tns.get_nbits(), x)
72*49fe348cSAndroid Build Coastguard Worker
73*49fe348cSAndroid Build Coastguard Worker        return pitch_present
74*49fe348cSAndroid Build Coastguard Worker
75*49fe348cSAndroid Build Coastguard Worker    def encode(self, pitch_present, nbytes):
76*49fe348cSAndroid Build Coastguard Worker
77*49fe348cSAndroid Build Coastguard Worker        b = bitstream.BitstreamWriter(nbytes)
78*49fe348cSAndroid Build Coastguard Worker
79*49fe348cSAndroid Build Coastguard Worker        self.bwdet.store(b)
80*49fe348cSAndroid Build Coastguard Worker
81*49fe348cSAndroid Build Coastguard Worker        self.spec.store(b)
82*49fe348cSAndroid Build Coastguard Worker
83*49fe348cSAndroid Build Coastguard Worker        self.tns.store(b)
84*49fe348cSAndroid Build Coastguard Worker
85*49fe348cSAndroid Build Coastguard Worker        b.write_bit(pitch_present)
86*49fe348cSAndroid Build Coastguard Worker
87*49fe348cSAndroid Build Coastguard Worker        self.sns.store(b)
88*49fe348cSAndroid Build Coastguard Worker
89*49fe348cSAndroid Build Coastguard Worker        if pitch_present:
90*49fe348cSAndroid Build Coastguard Worker            self.ltpf.store(b)
91*49fe348cSAndroid Build Coastguard Worker
92*49fe348cSAndroid Build Coastguard Worker        self.spec.encode(b)
93*49fe348cSAndroid Build Coastguard Worker
94*49fe348cSAndroid Build Coastguard Worker        return b.terminate()
95*49fe348cSAndroid Build Coastguard Worker
96*49fe348cSAndroid Build Coastguard Worker    def run(self, x, nbytes):
97*49fe348cSAndroid Build Coastguard Worker
98*49fe348cSAndroid Build Coastguard Worker        pitch_present = self.analyse(x, nbytes)
99*49fe348cSAndroid Build Coastguard Worker
100*49fe348cSAndroid Build Coastguard Worker        data = self.encode(pitch_present, nbytes)
101*49fe348cSAndroid Build Coastguard Worker
102*49fe348cSAndroid Build Coastguard Worker        return data
103*49fe348cSAndroid Build Coastguard Worker
104*49fe348cSAndroid Build Coastguard Workerdef check_appendix_c(dt):
105*49fe348cSAndroid Build Coastguard Worker
106*49fe348cSAndroid Build Coastguard Worker    i0 = dt - T.DT_7M5
107*49fe348cSAndroid Build Coastguard Worker
108*49fe348cSAndroid Build Coastguard Worker    enc_c = lc3.setup_encoder(int(T.DT_MS[dt] * 1000), 16000)
109*49fe348cSAndroid Build Coastguard Worker    ok = True
110*49fe348cSAndroid Build Coastguard Worker
111*49fe348cSAndroid Build Coastguard Worker    for i in range(len(C.X_PCM[i0])):
112*49fe348cSAndroid Build Coastguard Worker
113*49fe348cSAndroid Build Coastguard Worker        data = lc3.encode(enc_c, C.X_PCM[i0][i], C.NBYTES[i0])
114*49fe348cSAndroid Build Coastguard Worker        ok = ok and data == C.BYTES_AC[i0][i]
115*49fe348cSAndroid Build Coastguard Worker
116*49fe348cSAndroid Build Coastguard Worker    return ok
117*49fe348cSAndroid Build Coastguard Worker
118*49fe348cSAndroid Build Coastguard Workerdef check():
119*49fe348cSAndroid Build Coastguard Worker
120*49fe348cSAndroid Build Coastguard Worker    ok = True
121*49fe348cSAndroid Build Coastguard Worker
122*49fe348cSAndroid Build Coastguard Worker    for dt in ( T.DT_7M5, T.DT_10M ):
123*49fe348cSAndroid Build Coastguard Worker        ok = ok and check_appendix_c(dt)
124*49fe348cSAndroid Build Coastguard Worker
125*49fe348cSAndroid Build Coastguard Worker    return ok
126