[libsigrok.git] / src / input / wav.c

/*
 * This file is part of the libsigrok project.
 *
 * Copyright (C) 2013 Bert Vermeulen <bert@biot.com>
 *
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program.  If not, see <http://www.gnu.org/licenses/>.
 */

#include <config.h>
#include <sys/types.h>
#include <sys/stat.h>
#include <unistd.h>
#include <fcntl.h>
#include <ctype.h>
#include <string.h>
#include <stdint.h>
#include <libsigrok/libsigrok.h>
#include "libsigrok-internal.h"

#define LOG_PREFIX "input/wav"

/* How many bytes at a time to process and send to the session bus. */
#define CHUNK_SIZE 4096

/* Minimum size of header + 1 8-bit mono PCM sample. */
#define MIN_DATA_CHUNK_OFFSET    45

/* Expect to find the "data" chunk within this offset from the start. */
#define MAX_DATA_CHUNK_OFFSET    1024

#define WAVE_FORMAT_PCM_         0x0001
#define WAVE_FORMAT_IEEE_FLOAT_  0x0003
#define WAVE_FORMAT_EXTENSIBLE_  0xfffe

struct context {
	gboolean started;
	int fmt_code;
	uint64_t samplerate;
	int samplesize;
	int num_channels;
	int unitsize;
	gboolean found_data;
};

static int parse_wav_header(GString *buf, struct context *inc)
{
	uint64_t samplerate;
	unsigned int fmt_code, samplesize, num_channels, unitsize;

	if (buf->len < MIN_DATA_CHUNK_OFFSET)
		return SR_ERR_NA;

	fmt_code = RL16(buf->str + 20);
	samplerate = RL32(buf->str + 24);

	samplesize = RL16(buf->str + 32);
	num_channels = RL16(buf->str + 22);
	if (num_channels == 0)
		return SR_ERR;
	unitsize = samplesize / num_channels;
	if (unitsize != 1 && unitsize != 2 && unitsize != 4) {
		sr_err("Only 8, 16 or 32 bits per sample supported.");
		return SR_ERR_DATA;
	}

	if (fmt_code == WAVE_FORMAT_PCM_) {
	} else if (fmt_code == WAVE_FORMAT_IEEE_FLOAT_) {
		if (unitsize != 4) {
			sr_err("only 32-bit floats supported.");
			return SR_ERR_DATA;
		}
	} else if (fmt_code == WAVE_FORMAT_EXTENSIBLE_) {
		if (buf->len < 70)
			/* Not enough for extensible header and next chunk. */
			return SR_ERR_NA;

		if (RL16(buf->str + 16) != 40) {
			sr_err("WAV extensible format chunk must be 40 bytes.");
			return SR_ERR;
		}
		if (RL16(buf->str + 36) != 22) {
			sr_err("WAV extension must be 22 bytes.");
			return SR_ERR;
		}
		if (RL16(buf->str + 34) != RL16(buf->str + 38)) {
			sr_err("Reduced valid bits per sample not supported.");
			return SR_ERR_DATA;
		}
		/* Real format code is the first two bytes of the GUID. */
		fmt_code = RL16(buf->str + 44);
		if (fmt_code != WAVE_FORMAT_PCM_ && fmt_code != WAVE_FORMAT_IEEE_FLOAT_) {
			sr_err("Only PCM and floating point samples are supported.");
			return SR_ERR_DATA;
		}
		if (fmt_code == WAVE_FORMAT_IEEE_FLOAT_ && unitsize != 4) {
			sr_err("only 32-bit floats supported.");
			return SR_ERR_DATA;
		}
	} else {
		sr_err("Only PCM and floating point samples are supported.");
		return SR_ERR_DATA;
	}

	if (inc) {
		inc->fmt_code = fmt_code;
		inc->samplerate = samplerate;
		inc->samplesize = samplesize;
		inc->num_channels = num_channels;
		inc->unitsize = unitsize;
		inc->found_data = FALSE;
	}

	return SR_OK;
}

static int format_match(GHashTable *metadata)
{
	GString *buf;
	int ret;

	buf = g_hash_table_lookup(metadata, GINT_TO_POINTER(SR_INPUT_META_HEADER));
	if (strncmp(buf->str, "RIFF", 4))
		return SR_ERR;
	if (strncmp(buf->str + 8, "WAVE", 4))
		return SR_ERR;
	if (strncmp(buf->str + 12, "fmt ", 4))
		return SR_ERR;
	/*
	 * Only gets called when we already know this is a WAV file, so
	 * this parser can log error messages.
	 */
	if ((ret = parse_wav_header(buf, NULL)) != SR_OK)
		return ret;

	return SR_OK;
}

static int init(struct sr_input *in, GHashTable *options)
{
	(void)options;

	in->sdi = g_malloc0(sizeof(struct sr_dev_inst));
	in->priv = g_malloc0(sizeof(struct context));

	return SR_OK;
}

static int find_data_chunk(GString *buf, int initial_offset)
{
	unsigned int offset, i;

	offset = initial_offset;
	while (offset < MIN(MAX_DATA_CHUNK_OFFSET, buf->len)) {
		if (!memcmp(buf->str + offset, "data", 4))
			/* Skip into the samples. */
			return offset + 8;
		for (i = 0; i < 4; i++) {
			if (!isalnum(buf->str[offset + i])
					&& !isblank(buf->str[offset + i]))
				/* Doesn't look like a chunk ID. */
				return -1;
		}
		/* Skip past this chunk. */
		offset += 8 + RL32(buf->str + offset + 4);
	}

	if (offset > MAX_DATA_CHUNK_OFFSET)
		return -1;

	return offset;
}

static void send_chunk(const struct sr_input *in, int offset, int num_samples)
{
	struct sr_datafeed_packet packet;
	struct sr_datafeed_analog_old analog;
	struct context *inc;
	float fdata[CHUNK_SIZE];
	int total_samples, samplenum;
	char *s, *d;

	inc = in->priv;

	s = in->buf->str + offset;
	d = (char *)fdata;
	memset(fdata, 0, CHUNK_SIZE);
	total_samples = num_samples * inc->num_channels;
	for (samplenum = 0; samplenum < total_samples; samplenum++) {
		if (inc->fmt_code == WAVE_FORMAT_PCM_) {
			switch (inc->unitsize) {
			case 1:
				/* 8-bit PCM samples are unsigned. */
				fdata[samplenum] = *(uint8_t*)(s) / (float)255;
				break;
			case 2:
				fdata[samplenum] = RL16S(s) / (float)INT16_MAX;
				break;
			case 4:
				fdata[samplenum] = RL32S(s) / (float)INT32_MAX;
				break;
			}
		} else {
			/* BINARY32 float */
#ifdef WORDS_BIGENDIAN
			int i;
			for (i = 0; i < inc->unitsize; i++)
				d[i] = s[inc->unitsize - 1 - i];
#else
			memcpy(d, s, inc->unitsize);
#endif
		}
		s += inc->unitsize;
		d += inc->unitsize;
	}
	packet.type = SR_DF_ANALOG_OLD;
	packet.payload = &analog;
	analog.channels = in->sdi->channels;
	analog.num_samples = num_samples;
	analog.mq = 0;
	analog.mqflags = 0;
	analog.unit = 0;
	analog.data = fdata;
	sr_session_send(in->sdi, &packet);
}

static int process_buffer(struct sr_input *in)
{
	struct context *inc;
	struct sr_datafeed_packet packet;
	struct sr_datafeed_meta meta;
	struct sr_config *src;
	int offset, chunk_samples, total_samples, processed, max_chunk_samples;
	int num_samples, i;

	inc = in->priv;
	if (!inc->started) {
		std_session_send_df_header(in->sdi, LOG_PREFIX);

		packet.type = SR_DF_META;
		packet.payload = &meta;
		src = sr_config_new(SR_CONF_SAMPLERATE, g_variant_new_uint64(inc->samplerate));
		meta.config = g_slist_append(NULL, src);
		sr_session_send(in->sdi, &packet);
		g_slist_free(meta.config);
		sr_config_free(src);

		inc->started = TRUE;
	}

	if (!inc->found_data) {
		/* Skip past size of 'fmt ' chunk. */
		i = 20 + RL32(in->buf->str + 16);
		offset = find_data_chunk(in->buf, i);
		if (offset < 0) {
			if (in->buf->len > MAX_DATA_CHUNK_OFFSET) {
				sr_err("Couldn't find data chunk.");
				return SR_ERR;
			}
		}
		inc->found_data = TRUE;
	} else
		offset = 0;

	/* Round off up to the last channels * unitsize boundary. */
	chunk_samples = (in->buf->len - offset) / inc->samplesize;
	max_chunk_samples = CHUNK_SIZE / inc->samplesize;
	processed = 0;
	total_samples = chunk_samples;
	while (processed < total_samples) {
		if (chunk_samples > max_chunk_samples)
			num_samples = max_chunk_samples;
		else
			num_samples = chunk_samples;
		send_chunk(in, offset, num_samples);
		offset += num_samples * inc->samplesize;
		chunk_samples -= num_samples;
		processed += num_samples;
	}

	if ((unsigned int)offset < in->buf->len) {
		/*
		 * The incoming buffer wasn't processed completely. Stash
		 * the leftover data for next time.
		 */
		g_string_erase(in->buf, 0, offset);
	} else
		g_string_truncate(in->buf, 0);

	return SR_OK;
}

static int receive(struct sr_input *in, GString *buf)
{
	struct context *inc;
	int ret;
	char channelname[8];

	g_string_append_len(in->buf, buf->str, buf->len);

	if (in->buf->len < MIN_DATA_CHUNK_OFFSET) {
		/*
		 * Don't even try until there's enough room
		 * for the data segment to start.
		 */
		return SR_OK;
	}

	inc = in->priv;
	if (!in->sdi_ready) {
		if ((ret = parse_wav_header(in->buf, inc)) == SR_ERR_NA)
			/* Not enough data yet. */
			return SR_OK;
		else if (ret != SR_OK)
			return ret;

		for (int i = 0; i < inc->num_channels; i++) {
			snprintf(channelname, 8, "CH%d", i + 1);
			sr_channel_new(in->sdi, i, SR_CHANNEL_ANALOG, TRUE, channelname);
		}

		/* sdi is ready, notify frontend. */
		in->sdi_ready = TRUE;
		return SR_OK;
	}

	ret = process_buffer(in);

	return ret;
}

static int end(struct sr_input *in)
{
	struct context *inc;
	int ret;

	if (in->sdi_ready)
		ret = process_buffer(in);
	else
		ret = SR_OK;

	inc = in->priv;
	if (inc->started)
		std_session_send_df_end(in->sdi, LOG_PREFIX);

	return ret;
}

SR_PRIV struct sr_input_module input_wav = {
	.id = "wav",
	.name = "WAV",
	.desc = "WAV file",
	.exts = (const char*[]){"wav", NULL},
	.metadata = { SR_INPUT_META_HEADER | SR_INPUT_META_REQUIRED },
	.format_match = format_match,
	.init = init,
	.receive = receive,
	.end = end,
};
Commit	Line	Data
	1	/*
	2	* This file is part of the libsigrok project.
	3	*
	4	* Copyright (C) 2013 Bert Vermeulen <bert@biot.com>
	5	*
	6	* This program is free software: you can redistribute it and/or modify
	7	* it under the terms of the GNU General Public License as published by
	8	* the Free Software Foundation, either version 3 of the License, or
	9	* (at your option) any later version.
	10	*
	11	* This program is distributed in the hope that it will be useful,
	12	* but WITHOUT ANY WARRANTY; without even the implied warranty of
	13	* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
	14	* GNU General Public License for more details.
	15	*
	16	* You should have received a copy of the GNU General Public License
	17	* along with this program. If not, see <http://www.gnu.org/licenses/>.
	18	*/
	19
	20	#include <config.h>
	21	#include <sys/types.h>
	22	#include <sys/stat.h>
	23	#include <unistd.h>
	24	#include <fcntl.h>
	25	#include <ctype.h>
	26	#include <string.h>
	27	#include <stdint.h>
	28	#include <libsigrok/libsigrok.h>
	29	#include "libsigrok-internal.h"
	30
	31	#define LOG_PREFIX "input/wav"
	32
	33	/* How many bytes at a time to process and send to the session bus. */
	34	#define CHUNK_SIZE 4096
	35
	36	/* Minimum size of header + 1 8-bit mono PCM sample. */
	37	#define MIN_DATA_CHUNK_OFFSET 45
	38
	39	/* Expect to find the "data" chunk within this offset from the start. */
	40	#define MAX_DATA_CHUNK_OFFSET 1024
	41
	42	#define WAVE_FORMAT_PCM_ 0x0001
	43	#define WAVE_FORMAT_IEEE_FLOAT_ 0x0003
	44	#define WAVE_FORMAT_EXTENSIBLE_ 0xfffe
	45
	46	struct context {
	47	gboolean started;
	48	int fmt_code;
	49	uint64_t samplerate;
	50	int samplesize;
	51	int num_channels;
	52	int unitsize;
	53	gboolean found_data;
	54	};
	55
	56	static int parse_wav_header(GString buf, struct context inc)
	57	{
	58	uint64_t samplerate;
	59	unsigned int fmt_code, samplesize, num_channels, unitsize;
	60
	61	if (buf->len < MIN_DATA_CHUNK_OFFSET)
	62	return SR_ERR_NA;
	63
	64	fmt_code = RL16(buf->str + 20);
	65	samplerate = RL32(buf->str + 24);
	66
	67	samplesize = RL16(buf->str + 32);
	68	num_channels = RL16(buf->str + 22);
	69	if (num_channels == 0)
	70	return SR_ERR;
	71	unitsize = samplesize / num_channels;
	72	if (unitsize != 1 && unitsize != 2 && unitsize != 4) {
	73	sr_err("Only 8, 16 or 32 bits per sample supported.");
	74	return SR_ERR_DATA;
	75	}
	76
	77	if (fmt_code == WAVE_FORMAT_PCM_) {
	78	} else if (fmt_code == WAVE_FORMAT_IEEE_FLOAT_) {
	79	if (unitsize != 4) {
	80	sr_err("only 32-bit floats supported.");
	81	return SR_ERR_DATA;
	82	}
	83	} else if (fmt_code == WAVE_FORMAT_EXTENSIBLE_) {
	84	if (buf->len < 70)
	85	/* Not enough for extensible header and next chunk. */
	86	return SR_ERR_NA;
	87
	88	if (RL16(buf->str + 16) != 40) {
	89	sr_err("WAV extensible format chunk must be 40 bytes.");
	90	return SR_ERR;
	91	}
	92	if (RL16(buf->str + 36) != 22) {
	93	sr_err("WAV extension must be 22 bytes.");
	94	return SR_ERR;
	95	}
	96	if (RL16(buf->str + 34) != RL16(buf->str + 38)) {
	97	sr_err("Reduced valid bits per sample not supported.");
	98	return SR_ERR_DATA;
	99	}
	100	/* Real format code is the first two bytes of the GUID. */
	101	fmt_code = RL16(buf->str + 44);
	102	if (fmt_code != WAVE_FORMAT_PCM_ && fmt_code != WAVE_FORMAT_IEEE_FLOAT_) {
	103	sr_err("Only PCM and floating point samples are supported.");
	104	return SR_ERR_DATA;
	105	}
	106	if (fmt_code == WAVE_FORMAT_IEEE_FLOAT_ && unitsize != 4) {
	107	sr_err("only 32-bit floats supported.");
	108	return SR_ERR_DATA;
	109	}
	110	} else {
	111	sr_err("Only PCM and floating point samples are supported.");
	112	return SR_ERR_DATA;
	113	}
	114
	115	if (inc) {
	116	inc->fmt_code = fmt_code;
	117	inc->samplerate = samplerate;
	118	inc->samplesize = samplesize;
	119	inc->num_channels = num_channels;
	120	inc->unitsize = unitsize;
	121	inc->found_data = FALSE;
	122	}
	123
	124	return SR_OK;
	125	}
	126
	127	static int format_match(GHashTable *metadata)
	128	{
	129	GString *buf;
	130	int ret;
	131
	132	buf = g_hash_table_lookup(metadata, GINT_TO_POINTER(SR_INPUT_META_HEADER));
	133	if (strncmp(buf->str, "RIFF", 4))
	134	return SR_ERR;
	135	if (strncmp(buf->str + 8, "WAVE", 4))
	136	return SR_ERR;
	137	if (strncmp(buf->str + 12, "fmt ", 4))
	138	return SR_ERR;
	139	/*
	140	* Only gets called when we already know this is a WAV file, so
	141	* this parser can log error messages.
	142	*/
	143	if ((ret = parse_wav_header(buf, NULL)) != SR_OK)
	144	return ret;
	145
	146	return SR_OK;
	147	}
	148
	149	static int init(struct sr_input in, GHashTable options)
	150	{
	151	(void)options;
	152
	153	in->sdi = g_malloc0(sizeof(struct sr_dev_inst));
	154	in->priv = g_malloc0(sizeof(struct context));
	155
	156	return SR_OK;
	157	}
	158
	159	static int find_data_chunk(GString *buf, int initial_offset)
	160	{
	161	unsigned int offset, i;
	162
	163	offset = initial_offset;
	164	while (offset < MIN(MAX_DATA_CHUNK_OFFSET, buf->len)) {
	165	if (!memcmp(buf->str + offset, "data", 4))
	166	/* Skip into the samples. */
	167	return offset + 8;
	168	for (i = 0; i < 4; i++) {
	169	if (!isalnum(buf->str[offset + i])
	170	&& !isblank(buf->str[offset + i]))
	171	/* Doesn't look like a chunk ID. */
	172	return -1;
	173	}
	174	/* Skip past this chunk. */
	175	offset += 8 + RL32(buf->str + offset + 4);
	176	}
	177
	178	if (offset > MAX_DATA_CHUNK_OFFSET)
	179	return -1;
	180
	181	return offset;
	182	}
	183
	184	static void send_chunk(const struct sr_input *in, int offset, int num_samples)
	185	{
	186	struct sr_datafeed_packet packet;
	187	struct sr_datafeed_analog_old analog;
	188	struct context *inc;
	189	float fdata[CHUNK_SIZE];
	190	int total_samples, samplenum;
	191	char s, d;
	192
	193	inc = in->priv;
	194
	195	s = in->buf->str + offset;
	196	d = (char *)fdata;
	197	memset(fdata, 0, CHUNK_SIZE);
	198	total_samples = num_samples * inc->num_channels;
	199	for (samplenum = 0; samplenum < total_samples; samplenum++) {
	200	if (inc->fmt_code == WAVE_FORMAT_PCM_) {
	201	switch (inc->unitsize) {
	202	case 1:
	203	/* 8-bit PCM samples are unsigned. */
	204	fdata[samplenum] = (uint8_t)(s) / (float)255;
	205	break;
	206	case 2:
	207	fdata[samplenum] = RL16S(s) / (float)INT16_MAX;
	208	break;
	209	case 4:
	210	fdata[samplenum] = RL32S(s) / (float)INT32_MAX;
	211	break;
	212	}
	213	} else {
	214	/* BINARY32 float */
	215	#ifdef WORDS_BIGENDIAN
	216	int i;
	217	for (i = 0; i < inc->unitsize; i++)
	218	d[i] = s[inc->unitsize - 1 - i];
	219	#else
	220	memcpy(d, s, inc->unitsize);
	221	#endif
	222	}
	223	s += inc->unitsize;
	224	d += inc->unitsize;
	225	}
	226	packet.type = SR_DF_ANALOG_OLD;
	227	packet.payload = &analog;
	228	analog.channels = in->sdi->channels;
	229	analog.num_samples = num_samples;
	230	analog.mq = 0;
	231	analog.mqflags = 0;
	232	analog.unit = 0;
	233	analog.data = fdata;
	234	sr_session_send(in->sdi, &packet);
	235	}
	236
	237	static int process_buffer(struct sr_input *in)
	238	{
	239	struct context *inc;
	240	struct sr_datafeed_packet packet;
	241	struct sr_datafeed_meta meta;
	242	struct sr_config *src;
	243	int offset, chunk_samples, total_samples, processed, max_chunk_samples;
	244	int num_samples, i;
	245
	246	inc = in->priv;
	247	if (!inc->started) {
	248	std_session_send_df_header(in->sdi, LOG_PREFIX);
	249
	250	packet.type = SR_DF_META;
	251	packet.payload = &meta;
	252	src = sr_config_new(SR_CONF_SAMPLERATE, g_variant_new_uint64(inc->samplerate));
	253	meta.config = g_slist_append(NULL, src);
	254	sr_session_send(in->sdi, &packet);
	255	g_slist_free(meta.config);
	256	sr_config_free(src);
	257
	258	inc->started = TRUE;
	259	}
	260
	261	if (!inc->found_data) {
	262	/* Skip past size of 'fmt ' chunk. */
	263	i = 20 + RL32(in->buf->str + 16);
	264	offset = find_data_chunk(in->buf, i);
	265	if (offset < 0) {
	266	if (in->buf->len > MAX_DATA_CHUNK_OFFSET) {
	267	sr_err("Couldn't find data chunk.");
	268	return SR_ERR;
	269	}
	270	}
	271	inc->found_data = TRUE;
	272	} else
	273	offset = 0;
	274
	275	/* Round off up to the last channels * unitsize boundary. */
	276	chunk_samples = (in->buf->len - offset) / inc->samplesize;
	277	max_chunk_samples = CHUNK_SIZE / inc->samplesize;
	278	processed = 0;
	279	total_samples = chunk_samples;
	280	while (processed < total_samples) {
	281	if (chunk_samples > max_chunk_samples)
	282	num_samples = max_chunk_samples;
	283	else
	284	num_samples = chunk_samples;
	285	send_chunk(in, offset, num_samples);
	286	offset += num_samples * inc->samplesize;
	287	chunk_samples -= num_samples;
	288	processed += num_samples;
	289	}
	290
	291	if ((unsigned int)offset < in->buf->len) {
	292	/*
	293	* The incoming buffer wasn't processed completely. Stash
	294	* the leftover data for next time.
	295	*/
	296	g_string_erase(in->buf, 0, offset);
	297	} else
	298	g_string_truncate(in->buf, 0);
	299
	300	return SR_OK;
	301	}
	302
	303	static int receive(struct sr_input in, GString buf)
	304	{
	305	struct context *inc;
	306	int ret;
	307	char channelname[8];
	308
	309	g_string_append_len(in->buf, buf->str, buf->len);
	310
	311	if (in->buf->len < MIN_DATA_CHUNK_OFFSET) {
	312	/*
	313	* Don't even try until there's enough room
	314	* for the data segment to start.
	315	*/
	316	return SR_OK;
	317	}
	318
	319	inc = in->priv;
	320	if (!in->sdi_ready) {
	321	if ((ret = parse_wav_header(in->buf, inc)) == SR_ERR_NA)
	322	/* Not enough data yet. */
	323	return SR_OK;
	324	else if (ret != SR_OK)
	325	return ret;
	326
	327	for (int i = 0; i < inc->num_channels; i++) {
	328	snprintf(channelname, 8, "CH%d", i + 1);
	329	sr_channel_new(in->sdi, i, SR_CHANNEL_ANALOG, TRUE, channelname);
	330	}
	331
	332	/* sdi is ready, notify frontend. */
	333	in->sdi_ready = TRUE;
	334	return SR_OK;
	335	}
	336
	337	ret = process_buffer(in);
	338
	339	return ret;
	340	}
	341
	342	static int end(struct sr_input *in)
	343	{
	344	struct context *inc;
	345	int ret;
	346
	347	if (in->sdi_ready)
	348	ret = process_buffer(in);
	349	else
	350	ret = SR_OK;
	351
	352	inc = in->priv;
	353	if (inc->started)
	354	std_session_send_df_end(in->sdi, LOG_PREFIX);
	355
	356	return ret;
	357	}
	358
	359	SR_PRIV struct sr_input_module input_wav = {
	360	.id = "wav",
	361	.name = "WAV",
	362	.desc = "WAV file",
	363	.exts = (const char*[]){"wav", NULL},
	364	.metadata = { SR_INPUT_META_HEADER \| SR_INPUT_META_REQUIRED },
	365	.format_match = format_match,
	366	.init = init,
	367	.receive = receive,
	368	.end = end,
	369	};