| /* |
| * NoiseLevelEstimation.cpp |
| * |
| * Copyright 2021 HIMSA II K/S - www.himsa.com. Represented by EHIMA - |
| * www.ehima.com |
| * |
| * Licensed under the Apache License, Version 2.0 (the "License"); |
| * you may not use this file except in compliance with the License. |
| * You may obtain a copy of the License at |
| * |
| * http://www.apache.org/licenses/LICENSE-2.0 |
| * |
| * Unless required by applicable law or agreed to in writing, software |
| * distributed under the License is distributed on an "AS IS" BASIS, |
| * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. |
| * See the License for the specific language governing permissions and |
| * limitations under the License. |
| */ |
| |
| #include "NoiseLevelEstimation.hpp" |
| |
| #include <algorithm> |
| #include <cmath> |
| |
| namespace Lc3Enc { |
| |
| NoiseLevelEstimation::NoiseLevelEstimation(const Lc3Config& lc3Config_) |
| : lc3Config(lc3Config_), F_NF(0), datapoints(nullptr) {} |
| |
| NoiseLevelEstimation::~NoiseLevelEstimation() {} |
| |
| void NoiseLevelEstimation::run(const double* const X_f, |
| const int16_t* const X_q, uint8_t P_bw, |
| double gg) { |
| const uint16_t bw_stop_table_10ms[5] = {80, 160, 240, 320, 400}; |
| const uint16_t bw_stop_table_7p5ms[5] = {60, 120, 180, 240, 300}; |
| const uint16_t bw_stop = (Lc3Config::FrameDuration::d10ms == lc3Config.N_ms) |
| ? bw_stop_table_10ms[P_bw] |
| : bw_stop_table_7p5ms[P_bw]; |
| uint16_t NF_start = |
| (Lc3Config::FrameDuration::d10ms == lc3Config.N_ms) ? 24 : 18; |
| uint16_t NF_width = |
| (Lc3Config::FrameDuration::d10ms == lc3Config.N_ms) ? 3 : 2; |
| |
| double L_NF_numerator = 0; |
| uint16_t L_NF_denominator = 0; |
| for (uint16_t k = NF_start; (k < lc3Config.NE) && (k < bw_stop); k++) { |
| bool I_NF_k = true; |
| for (uint16_t i = k - NF_width; (i < bw_stop) && (i <= k + NF_width); i++) { |
| if (X_q[i] != 0) { |
| // -> I_NF[k] = 0; |
| I_NF_k = false; |
| break; |
| } |
| } |
| if (I_NF_k) { |
| // L_NF_numerator += fabs(X_f[k]); |
| L_NF_numerator += (X_f[k] > 0) ? X_f[k] : -X_f[k]; |
| L_NF_denominator++; |
| } |
| } |
| double L_NF = 0; |
| if ((L_NF_denominator > 0) && (gg > 0)) { |
| // L_NF_numerator /= gg; |
| L_NF = L_NF_numerator / gg / L_NF_denominator; |
| } |
| |
| double diff = 8 - 16 * L_NF; |
| if (diff >= 0) { |
| F_NF = std::min(static_cast<int>(diff + 0.5), |
| 7); // specification says rounding, |
| // F_NF = std::min( static_cast<int>(ceil(diff)), 7 ); // but sometimes it |
| // seems like the reference encoder computes ceil(.) |
| } else { |
| F_NF = 0; |
| } |
| |
| if (nullptr != datapoints) { |
| datapoints->log("L_NF", &L_NF, sizeof(L_NF)); |
| datapoints->log("L_NF_numerator", &L_NF_numerator, sizeof(L_NF_numerator)); |
| datapoints->log("L_NF_denominator", &L_NF_denominator, |
| sizeof(L_NF_denominator)); |
| } |
| } |
| |
| void NoiseLevelEstimation::registerDatapoints(DatapointContainer* datapoints_) { |
| datapoints = datapoints_; |
| if (nullptr != datapoints) { |
| datapoints->addDatapoint("F_NF", &F_NF, sizeof(F_NF)); |
| } |
| } |
| |
| } // namespace Lc3Enc |