Binaural synthesis (confbridge): Adds utils/conf_bridge_binaural_hrir_importer

Adds the import tool for converting a HRIR database to hrirs.h

ASTERISK-26292

Change-Id: I51eb31b54c23ffd9b544bdc6a09d20c112c8a547

1 files changed, 148 insertions, 0 deletions

diff --git a/utils/conf_bridge_binaural_hrir_importer.c b/utils/conf_bridge_binaural_hrir_importer.c
new file mode 100644
index 000000000..5690d86ea
--- /dev/null
+++ b/utils/conf_bridge_binaural_hrir_importer.c
@@ -0,0 +1,148 @@
+/*
+ * Asterisk -- An open source telephony toolkit.
+ *
+ * Copyright (C) 2016, Digium, Inc.
+ *
+ * Frank Haase <fra.haase@gmail.com>
+ * Dennis Guse <dennis.guse@alumni.tu-berlin.de>
+ *
+ * See http://www.asterisk.org for more information about
+ * the Asterisk project. Please do not directly contact
+ * any of the maintainers of this project for assistance;
+ * the project provides a web site, mailing lists and IRC
+ * channels for your use.
+ *
+ * This program is free software, distributed under the terms of
+ * the GNU General Public License Version 2. See the LICENSE file
+ * at the top of the source tree.
+ */
+
+/*!
+ * \file
+ * Converts a Head Related Impulse Response (HRIR) database (a multi-channel wave) into a C header file.
+ * HRIR for the left ear and HRIR for right ear have to be interleaved.
+ * No further signal processing is applied (e.g., resampling).
+ *
+ * Info messages are printed to stderror and the generated header file to output.
+ */
+
+#include <stdio.h>
+#include <string.h>
+#include <stdlib.h>
+#include <sndfile.h>
+
+int main (int argc, char **argv)
+{
+	char *hrir_filename;
+	unsigned int binaural_index_start;
+	unsigned int binaural_index_end;
+
+	SNDFILE *hrir_file;
+	SF_INFO hrir_info;
+	float *hrir_data;
+
+	unsigned int impulse_response_index_start;
+	unsigned int impulse_response_index_end;
+
+	int j;
+	int ir_current;
+
+	if(argc != 4) {
+		puts("HRIR database to C header file converter.");
+		puts("Usage: conf_bridge_binaural_hrir_importer HRIR.wav INDEX_START INDEX_END > OUTPUT.h");
+		puts("Example: conf_bridge_binaural_hrir_importer hrirs.wav 0 180 > ../bridges/bridge_softmix/include/hrirs.h");
+
+		return -1;
+	}
+
+	/* Parse arguments */
+	hrir_filename = argv[1];
+	binaural_index_start = atoi(argv[2]);
+	binaural_index_end = atoi(argv[3]);
+
+	/* Read HRIR database */
+	hrir_file = sf_open(hrir_filename, SFM_READ, &hrir_info);
+	if(hrir_file == NULL) {
+		fprintf(stderr, "ERROR: Could not open HRIR database (%s).\n", hrir_filename);
+
+		return -1;
+	}
+	fprintf(stderr, "INFO: Opened HRIR database (%s) with: number channels: %d; samplerate: %d; samples per channel: %ld\n", hrir_filename, hrir_info.channels, hrir_info.samplerate, hrir_info.frames);
+
+	hrir_data = (float *)malloc(hrir_info.channels * hrir_info.frames * sizeof(float));
+	if(hrir_data == NULL) {
+		fprintf(stderr, "ERROR: Out of memory!");
+
+		return -1;
+	}
+
+	/* Channels are interleaved */
+	sf_read_float(hrir_file, hrir_data, hrir_info.channels * hrir_info.frames);
+	sf_close(hrir_file);
+
+	if(binaural_index_start >= binaural_index_end) {
+		fprintf(stderr, "ERROR: INDEX_START (%d) must be smaller than INDEX_END (%d).", binaural_index_start, binaural_index_end);
+		free(hrir_data);
+
+		return -1;
+	}
+
+	if (binaural_index_end * 2 >= hrir_info.channels) {
+		fprintf(stderr, "ERROR: END_INDEX (%d) is out of range for HRIR database (%s).\n", binaural_index_end, hrir_filename);
+		free(hrir_data);
+
+		return -1;
+	}
+
+	/* Convert indices */
+	impulse_response_index_start = 2 * binaural_index_start;
+	impulse_response_index_end = (binaural_index_end + 1) * 2;
+
+	/* Write header */
+	printf("//Used hrirs database: %s\n", hrir_filename);
+	printf("//Start index in database: %d\n", impulse_response_index_start);
+	printf("//End index in database: %d\n", impulse_response_index_end);
+
+	printf("#define HRIRS_IMPULSE_LEN %ld\n", hrir_info.frames);
+	printf("#define HRIRS_IMPULSE_SIZE %d\n", binaural_index_end - binaural_index_start + 1);
+	printf("#define HRIRS_SAMPLE_RATE %d\n", hrir_info.samplerate);
+
+	printf("float hrirs_left[HRIRS_IMPULSE_SIZE][HRIRS_IMPULSE_LEN] = {\n");
+	for (ir_current = impulse_response_index_start; ir_current < impulse_response_index_end; ir_current += 2) {
+			printf("{");
+
+			for (j = 0; j < hrir_info.frames - 1; j++) {
+				printf("%.16f,", hrir_data[ir_current * hrir_info.frames + j]);
+			}
+			/* Write last without trailing "," */
+			printf("%.16f", hrir_data[ir_current * hrir_info.frames + hrir_info.frames - 1]);
+
+			if (ir_current + 2 < impulse_response_index_end) {
+				printf("},\n");
+			}	else {
+				printf("}};");
+			}
+	}
+
+	printf("\nfloat hrirs_right[HRIRS_IMPULSE_SIZE][HRIRS_IMPULSE_LEN] = {\n");
+	for (ir_current = impulse_response_index_start + 1; ir_current < impulse_response_index_end + 1; ir_current += 2) {
+			printf("{");
+
+			for (j = 0; j < hrir_info.frames - 1; j++) {
+				printf("%.16f,", hrir_data[ir_current * hrir_info.frames + j]);
+			}
+			 /* Write last without trailing "," */
+			printf("%.16f", hrir_data[ir_current * hrir_info.frames + hrir_info.frames - 1]);
+
+			if (ir_current + 2 < impulse_response_index_end) {
+				printf("},\n");
+			}	else {
+				printf("}};");
+			}
+	}
+
+	fprintf(stderr, "INFO: Successfully converted: imported %d impulse responses.\n", impulse_response_index_end - impulse_response_index_start);
+	free(hrir_data);
+
+	return 0;
+}