|
@@ -1,5 +1,5 @@
|
|
|
/**
|
|
|
- * Copyright 2021 Xiaomi Corporation (authors: Fangjun Kuang)
|
|
|
+ * Copyright (c) 2022-2024 Xiaomi Corporation (authors: Fangjun Kuang)
|
|
|
*
|
|
|
* See LICENSE for clarification regarding multiple authors
|
|
|
*
|
|
@@ -20,10 +20,11 @@
|
|
|
|
|
|
#include <cassert>
|
|
|
#include <fstream>
|
|
|
-#include <iostream>
|
|
|
#include <utility>
|
|
|
#include <vector>
|
|
|
|
|
|
+#include "platform.h" // NOLINT
|
|
|
+
|
|
|
namespace sherpa_ncnn {
|
|
|
namespace {
|
|
|
// see http://soundfile.sapp.org/doc/WaveFormat/
|
|
@@ -31,53 +32,12 @@ namespace {
|
|
|
// Note: We assume little endian here
|
|
|
// TODO(fangjun): Support big endian
|
|
|
struct WaveHeader {
|
|
|
- bool Validate() const {
|
|
|
- // F F I R
|
|
|
- if (chunk_id != 0x46464952) {
|
|
|
- return false;
|
|
|
- }
|
|
|
- // E V A W
|
|
|
- if (format != 0x45564157) {
|
|
|
- return false;
|
|
|
- }
|
|
|
-
|
|
|
- // t m f
|
|
|
- if (subchunk1_id != 0x20746d66) {
|
|
|
- return false;
|
|
|
- }
|
|
|
-
|
|
|
- if (subchunk1_size != 16) { // 16 for PCM
|
|
|
- return false;
|
|
|
- }
|
|
|
-
|
|
|
- if (audio_format != 1) { // 1 for PCM
|
|
|
- return false;
|
|
|
- }
|
|
|
-
|
|
|
- if (num_channels != 1) { // we support only single channel for now
|
|
|
- return false;
|
|
|
- }
|
|
|
- if (byte_rate != (sample_rate * num_channels * bits_per_sample / 8)) {
|
|
|
- return false;
|
|
|
- }
|
|
|
-
|
|
|
- if (block_align != (num_channels * bits_per_sample / 8)) {
|
|
|
- return false;
|
|
|
- }
|
|
|
-
|
|
|
- if (bits_per_sample != 16) { // we support only 16 bits per sample
|
|
|
- return false;
|
|
|
- }
|
|
|
-
|
|
|
- return true;
|
|
|
- }
|
|
|
-
|
|
|
// See
|
|
|
// https://en.wikipedia.org/wiki/WAV#Metadata
|
|
|
// and
|
|
|
// https://www.robotplanet.dk/audio/wav_meta_data/riff_mci.pdf
|
|
|
void SeekToDataChunk(std::istream &is) {
|
|
|
- // a t a d
|
|
|
+ // a t a d
|
|
|
while (is && subchunk2_id != 0x61746164) {
|
|
|
// const char *p = reinterpret_cast<const char *>(&subchunk2_id);
|
|
|
// printf("Skip chunk (%x): %c%c%c%c of size: %d\n", subchunk2_id, p[0],
|
|
@@ -104,46 +64,255 @@ struct WaveHeader {
|
|
|
};
|
|
|
static_assert(sizeof(WaveHeader) == 44, "");
|
|
|
|
|
|
+/*
|
|
|
+sox int16-1-channel-zh.wav -b 8 int8-1-channel-zh.wav
|
|
|
+
|
|
|
+sox int16-1-channel-zh.wav -c 2 int16-2-channel-zh.wav
|
|
|
+
|
|
|
+we use audacity to generate int32-1-channel-zh.wav and float32-1-channel-zh.wav
|
|
|
+because sox uses WAVE_FORMAT_EXTENSIBLE, which is not easy to support
|
|
|
+in sherpa-ncnn.
|
|
|
+ */
|
|
|
+
|
|
|
// Read a wave file of mono-channel.
|
|
|
// Return its samples normalized to the range [-1, 1).
|
|
|
-std::vector<float> ReadWaveImpl(std::istream &is, float expected_sample_rate,
|
|
|
+std::vector<float> ReadWaveImpl(std::istream &is, int32_t *sampling_rate,
|
|
|
bool *is_ok) {
|
|
|
- WaveHeader header;
|
|
|
- is.read(reinterpret_cast<char *>(&header), sizeof(header));
|
|
|
- if (!is) {
|
|
|
+ WaveHeader header{};
|
|
|
+ is.read(reinterpret_cast<char *>(&header.chunk_id), sizeof(header.chunk_id));
|
|
|
+
|
|
|
+ // F F I R
|
|
|
+ if (header.chunk_id != 0x46464952) {
|
|
|
+ NCNN_LOGE("Expected chunk_id RIFF. Given: 0x%08x\n", header.chunk_id);
|
|
|
+ *is_ok = false;
|
|
|
+ return {};
|
|
|
+ }
|
|
|
+
|
|
|
+ is.read(reinterpret_cast<char *>(&header.chunk_size),
|
|
|
+ sizeof(header.chunk_size));
|
|
|
+
|
|
|
+ is.read(reinterpret_cast<char *>(&header.format), sizeof(header.format));
|
|
|
+
|
|
|
+ // E V A W
|
|
|
+ if (header.format != 0x45564157) {
|
|
|
+ NCNN_LOGE("Expected format WAVE. Given: 0x%08x\n", header.format);
|
|
|
*is_ok = false;
|
|
|
return {};
|
|
|
}
|
|
|
|
|
|
- if (!header.Validate()) {
|
|
|
+ is.read(reinterpret_cast<char *>(&header.subchunk1_id),
|
|
|
+ sizeof(header.subchunk1_id));
|
|
|
+
|
|
|
+ is.read(reinterpret_cast<char *>(&header.subchunk1_size),
|
|
|
+ sizeof(header.subchunk1_size));
|
|
|
+
|
|
|
+ if (header.subchunk1_id == 0x4b4e554a) {
|
|
|
+ // skip junk padding
|
|
|
+ is.seekg(header.subchunk1_size, std::istream::cur);
|
|
|
+
|
|
|
+ is.read(reinterpret_cast<char *>(&header.subchunk1_id),
|
|
|
+ sizeof(header.subchunk1_id));
|
|
|
+
|
|
|
+ is.read(reinterpret_cast<char *>(&header.subchunk1_size),
|
|
|
+ sizeof(header.subchunk1_size));
|
|
|
+ }
|
|
|
+
|
|
|
+ if (header.subchunk1_id != 0x20746d66) {
|
|
|
+ NCNN_LOGE("Expected subchunk1_id 0x20746d66. Given: 0x%08x\n",
|
|
|
+ header.subchunk1_id);
|
|
|
*is_ok = false;
|
|
|
return {};
|
|
|
}
|
|
|
|
|
|
- header.SeekToDataChunk(is);
|
|
|
- if (!is) {
|
|
|
+ // NAudio uses 18
|
|
|
+ // See https://github.com/naudio/NAudio/issues/1132
|
|
|
+ if (header.subchunk1_size != 16 &&
|
|
|
+ header.subchunk1_size != 18) { // 16 for PCM
|
|
|
+ NCNN_LOGE("Expected subchunk1_size 16. Given: %d\n", header.subchunk1_size);
|
|
|
*is_ok = false;
|
|
|
return {};
|
|
|
}
|
|
|
|
|
|
- if (expected_sample_rate != header.sample_rate) {
|
|
|
+ is.read(reinterpret_cast<char *>(&header.audio_format),
|
|
|
+ sizeof(header.audio_format));
|
|
|
+
|
|
|
+ if (header.audio_format != 1 && header.audio_format != 3) {
|
|
|
+ // 1 for integer PCM
|
|
|
+ // 3 for floating point PCM
|
|
|
+ // see https://www.mmsp.ece.mcgill.ca/Documents/AudioFormats/WAVE/WAVE.html
|
|
|
+ // and https://github.com/microsoft/DirectXTK/wiki/Wave-Formats
|
|
|
+ NCNN_LOGE("Expected audio_format 1. Given: %d\n", header.audio_format);
|
|
|
+
|
|
|
+ if (header.audio_format == static_cast<int16_t>(0xfffe)) {
|
|
|
+ NCNN_LOGE("We don't support WAVE_FORMAT_EXTENSIBLE files.");
|
|
|
+ }
|
|
|
+
|
|
|
*is_ok = false;
|
|
|
return {};
|
|
|
}
|
|
|
|
|
|
- // header.subchunk2_size contains the number of bytes in the data.
|
|
|
- // As we assume each sample contains two bytes, so it is divided by 2 here
|
|
|
- std::vector<int16_t> samples(header.subchunk2_size / 2);
|
|
|
+ is.read(reinterpret_cast<char *>(&header.num_channels),
|
|
|
+ sizeof(header.num_channels));
|
|
|
+
|
|
|
+ if (header.num_channels != 1) { // we support only single channel for now
|
|
|
+ NCNN_LOGE(
|
|
|
+ "Warning: %d channels are found. We only use the first channel.\n",
|
|
|
+ header.num_channels);
|
|
|
+ }
|
|
|
+
|
|
|
+ is.read(reinterpret_cast<char *>(&header.sample_rate),
|
|
|
+ sizeof(header.sample_rate));
|
|
|
|
|
|
- is.read(reinterpret_cast<char *>(samples.data()), header.subchunk2_size);
|
|
|
+ is.read(reinterpret_cast<char *>(&header.byte_rate),
|
|
|
+ sizeof(header.byte_rate));
|
|
|
+
|
|
|
+ is.read(reinterpret_cast<char *>(&header.block_align),
|
|
|
+ sizeof(header.block_align));
|
|
|
+
|
|
|
+ is.read(reinterpret_cast<char *>(&header.bits_per_sample),
|
|
|
+ sizeof(header.bits_per_sample));
|
|
|
+
|
|
|
+ if (header.byte_rate !=
|
|
|
+ (header.sample_rate * header.num_channels * header.bits_per_sample / 8)) {
|
|
|
+ NCNN_LOGE("Incorrect byte rate: %d. Expected: %d", header.byte_rate,
|
|
|
+ (header.sample_rate * header.num_channels *
|
|
|
+ header.bits_per_sample / 8));
|
|
|
+ *is_ok = false;
|
|
|
+ return {};
|
|
|
+ }
|
|
|
+
|
|
|
+ if (header.block_align !=
|
|
|
+ (header.num_channels * header.bits_per_sample / 8)) {
|
|
|
+ NCNN_LOGE("Incorrect block align: %d. Expected: %d\n", header.block_align,
|
|
|
+ (header.num_channels * header.bits_per_sample / 8));
|
|
|
+ *is_ok = false;
|
|
|
+ return {};
|
|
|
+ }
|
|
|
+
|
|
|
+ if (header.bits_per_sample != 8 && header.bits_per_sample != 16 &&
|
|
|
+ header.bits_per_sample != 32) {
|
|
|
+ NCNN_LOGE("Expected bits_per_sample 8, 16 or 32. Given: %d\n",
|
|
|
+ header.bits_per_sample);
|
|
|
+ *is_ok = false;
|
|
|
+ return {};
|
|
|
+ }
|
|
|
+
|
|
|
+ if (header.subchunk1_size == 18) {
|
|
|
+ // this is for NAudio. It puts extra bytes after bits_per_sample
|
|
|
+ // See
|
|
|
+ // https://github.com/naudio/NAudio/blob/master/NAudio.Core/Wave/WaveFormats/WaveFormat.cs#L223
|
|
|
+
|
|
|
+ int16_t extra_size = -1;
|
|
|
+ is.read(reinterpret_cast<char *>(&extra_size), sizeof(int16_t));
|
|
|
+ if (extra_size != 0) {
|
|
|
+ NCNN_LOGE(
|
|
|
+ "Extra size should be 0 for wave from NAudio. Current extra size "
|
|
|
+ "%d\n",
|
|
|
+ extra_size);
|
|
|
+ *is_ok = false;
|
|
|
+ return {};
|
|
|
+ }
|
|
|
+ }
|
|
|
+
|
|
|
+ is.read(reinterpret_cast<char *>(&header.subchunk2_id),
|
|
|
+ sizeof(header.subchunk2_id));
|
|
|
+
|
|
|
+ is.read(reinterpret_cast<char *>(&header.subchunk2_size),
|
|
|
+ sizeof(header.subchunk2_size));
|
|
|
+
|
|
|
+ header.SeekToDataChunk(is);
|
|
|
if (!is) {
|
|
|
*is_ok = false;
|
|
|
return {};
|
|
|
}
|
|
|
|
|
|
- std::vector<float> ans(samples.size());
|
|
|
- for (int32_t i = 0; i != ans.size(); ++i) {
|
|
|
- ans[i] = samples[i] / 32768.;
|
|
|
+ *sampling_rate = header.sample_rate;
|
|
|
+
|
|
|
+ std::vector<float> ans;
|
|
|
+
|
|
|
+ if (header.bits_per_sample == 16 && header.audio_format == 1) {
|
|
|
+ // header.subchunk2_size contains the number of bytes in the data.
|
|
|
+ // As we assume each sample contains two bytes, so it is divided by 2 here
|
|
|
+ std::vector<int16_t> samples(header.subchunk2_size / 2);
|
|
|
+
|
|
|
+ is.read(reinterpret_cast<char *>(samples.data()), header.subchunk2_size);
|
|
|
+ if (!is) {
|
|
|
+ NCNN_LOGE("Failed to read %d bytes", header.subchunk2_size);
|
|
|
+ *is_ok = false;
|
|
|
+ return {};
|
|
|
+ }
|
|
|
+
|
|
|
+ ans.resize(samples.size() / header.num_channels);
|
|
|
+
|
|
|
+ // samples are interleaved
|
|
|
+ for (int32_t i = 0; i != static_cast<int32_t>(ans.size()); ++i) {
|
|
|
+ ans[i] = samples[i * header.num_channels] / 32768.;
|
|
|
+ }
|
|
|
+ } else if (header.bits_per_sample == 8 && header.audio_format == 1) {
|
|
|
+ // number of samples == number of bytes for 8-bit encoded samples
|
|
|
+ //
|
|
|
+ // For 8-bit encoded samples, they are unsigned!
|
|
|
+ std::vector<uint8_t> samples(header.subchunk2_size);
|
|
|
+
|
|
|
+ is.read(reinterpret_cast<char *>(samples.data()), header.subchunk2_size);
|
|
|
+ if (!is) {
|
|
|
+ NCNN_LOGE("Failed to read %d bytes", header.subchunk2_size);
|
|
|
+ *is_ok = false;
|
|
|
+ return {};
|
|
|
+ }
|
|
|
+
|
|
|
+ ans.resize(samples.size() / header.num_channels);
|
|
|
+ for (int32_t i = 0; i != static_cast<int32_t>(ans.size()); ++i) {
|
|
|
+ // Note(fangjun): We want to normalize each sample into the range [-1, 1]
|
|
|
+ // Since each original sample is in the range [0, 256], dividing
|
|
|
+ // them by 128 converts them to the range [0, 2];
|
|
|
+ // so after subtracting 1, we get the range [-1, 1]
|
|
|
+ //
|
|
|
+ ans[i] = samples[i * header.num_channels] / 128. - 1;
|
|
|
+ }
|
|
|
+ } else if (header.bits_per_sample == 32 && header.audio_format == 1) {
|
|
|
+ // 32 here is for int32
|
|
|
+ //
|
|
|
+ // header.subchunk2_size contains the number of bytes in the data.
|
|
|
+ // As we assume each sample contains 4 bytes, so it is divided by 4 here
|
|
|
+ std::vector<int32_t> samples(header.subchunk2_size / 4);
|
|
|
+
|
|
|
+ is.read(reinterpret_cast<char *>(samples.data()), header.subchunk2_size);
|
|
|
+ if (!is) {
|
|
|
+ NCNN_LOGE("Failed to read %d bytes", header.subchunk2_size);
|
|
|
+ *is_ok = false;
|
|
|
+ return {};
|
|
|
+ }
|
|
|
+
|
|
|
+ ans.resize(samples.size() / header.num_channels);
|
|
|
+ for (int32_t i = 0; i != static_cast<int32_t>(ans.size()); ++i) {
|
|
|
+ ans[i] = static_cast<float>(samples[i * header.num_channels]) / (1 << 31);
|
|
|
+ }
|
|
|
+ } else if (header.bits_per_sample == 32 && header.audio_format == 3) {
|
|
|
+ // 32 here is for float32
|
|
|
+ //
|
|
|
+ // header.subchunk2_size contains the number of bytes in the data.
|
|
|
+ // As we assume each sample contains 4 bytes, so it is divided by 4 here
|
|
|
+ std::vector<float> samples(header.subchunk2_size / 4);
|
|
|
+
|
|
|
+ is.read(reinterpret_cast<char *>(samples.data()), header.subchunk2_size);
|
|
|
+ if (!is) {
|
|
|
+ NCNN_LOGE("Failed to read %d bytes", header.subchunk2_size);
|
|
|
+ *is_ok = false;
|
|
|
+ return {};
|
|
|
+ }
|
|
|
+
|
|
|
+ ans.resize(samples.size() / header.num_channels);
|
|
|
+ for (int32_t i = 0; i != static_cast<int32_t>(ans.size()); ++i) {
|
|
|
+ ans[i] = samples[i * header.num_channels];
|
|
|
+ }
|
|
|
+ } else {
|
|
|
+ NCNN_LOGE(
|
|
|
+ "Unsupported %d bits per sample and audio format: %d. Supported values "
|
|
|
+ "are: 8, 16, 32.",
|
|
|
+ header.bits_per_sample, header.audio_format);
|
|
|
+ *is_ok = false;
|
|
|
+ return {};
|
|
|
}
|
|
|
|
|
|
*is_ok = true;
|
|
@@ -152,15 +321,41 @@ std::vector<float> ReadWaveImpl(std::istream &is, float expected_sample_rate,
|
|
|
|
|
|
} // namespace
|
|
|
|
|
|
-std::vector<float> ReadWave(const std::string &filename,
|
|
|
- float expected_sample_rate, bool *is_ok) {
|
|
|
+std::vector<float> ReadWave(const std::string &filename, int32_t *sampling_rate,
|
|
|
+ bool *is_ok) {
|
|
|
std::ifstream is(filename, std::ifstream::binary);
|
|
|
- return ReadWave(is, expected_sample_rate, is_ok);
|
|
|
+ return ReadWave(is, sampling_rate, is_ok);
|
|
|
}
|
|
|
|
|
|
-std::vector<float> ReadWave(std::istream &is, float expected_sample_rate,
|
|
|
+std::vector<float> ReadWave(const std::string &filename,
|
|
|
+ int32_t expected_sampling_rate, bool *is_ok) {
|
|
|
+ int32_t sampling_rate = -1;
|
|
|
+ auto samples = ReadWave(filename, &sampling_rate, is_ok);
|
|
|
+ if (*is_ok && expected_sampling_rate != sampling_rate) {
|
|
|
+ *is_ok = false;
|
|
|
+ NCNN_LOGE("Expected sample rate: %d, actual sample rate: %d",
|
|
|
+ expected_sampling_rate, sampling_rate);
|
|
|
+ samples.clear();
|
|
|
+ }
|
|
|
+ return samples;
|
|
|
+}
|
|
|
+
|
|
|
+std::vector<float> ReadWave(std::istream &is, int32_t *sampling_rate,
|
|
|
bool *is_ok) {
|
|
|
- auto samples = ReadWaveImpl(is, expected_sample_rate, is_ok);
|
|
|
+ auto samples = ReadWaveImpl(is, sampling_rate, is_ok);
|
|
|
+ return samples;
|
|
|
+}
|
|
|
+
|
|
|
+std::vector<float> ReadWave(std::istream &is, int32_t expected_sampling_rate,
|
|
|
+ bool *is_ok) {
|
|
|
+ int32_t sampling_rate = -1;
|
|
|
+ auto samples = ReadWave(is, &sampling_rate, is_ok);
|
|
|
+ if (*is_ok && expected_sampling_rate != sampling_rate) {
|
|
|
+ *is_ok = false;
|
|
|
+ NCNN_LOGE("Expected sample rate: %d, actual sample rate: %d",
|
|
|
+ expected_sampling_rate, sampling_rate);
|
|
|
+ samples.clear();
|
|
|
+ }
|
|
|
return samples;
|
|
|
}
|
|
|
|