diff trunk/src/Modules/Output/FileOutputHTK.cc @ 277:6b4921704eb1

- Ported over HTK file output
- Added some more meat to the Slaney IIR gammatone implementation
- Ported over the AIM-MAT sf2003 parabola strobe algorithm
- Finished making the SAI implementation compile
- Ported over the strobe list class (now uses STL deques internally)

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/trunk/src/Modules/Output/FileOutputHTK.cc	Thu Feb 18 16:55:40 2010 +0000
@@ -0,0 +1,194 @@
+// Copyright 2006-2010, Thomas Walters
+//
+// AIM-C: A C++ implementation of the Auditory Image Model
+// http://www.acousticscale.org/AIMC
+//
+// 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/>.
+
+/*!
+ * \file
+ * \brief File output in the HTK format. 
+ *
+ * \author Tom Walters <[e-mail]> 
+ * \author Willem van Engen <[e-mail]>
+ * \date created 2006/10/30
+ * \version \$Id$
+ */
+
+#ifdef _WINDOWS
+#	include <direct.h> // for _mkdir&_rmdir
+#else
+#	include <sys/types.h>
+#	include <dirent.h> // for opendir&friends
+#endif
+#include <stdio.h>
+#include <string.h>
+#include <cmath>
+
+#include "Modules/Output/FileOutputHTK.h"
+
+namespace aimc {
+FileOutputHTK::FileOutputHTK(Parameters *params) : Module(params) {
+	file_handle_ = NULL;
+	header_written_ = false;
+  filename_[0] = '\0';
+  frame_period_ms_ = 0.0f;
+}
+
+FileOutputHTK::~FileOutputHTK() {
+	if (file_handle_ != NULL)
+		CloseFile();
+}
+
+bool FileOutputHTK::OpenFile(const char* filename, float frame_period_ms) {
+	if (file_handle_ != NULL) {
+	  LOG_ERROR(_T("Couldn't open output file. A file is already open."));
+  	return false;
+	}
+
+	// Check that the output file exists and is writeable
+	if ((file_handle_ = fopen(filename, "wb"))==NULL ) {
+		LOG_ERROR(_T("Couldn't open output file '%s' for writing."), filename);
+		return false;
+	}
+	strcpy(filename_, filename);
+	sample_count_ = 0;
+  frame_period_ms_ = frame_period_ms;
+  header_written_ = false;
+	return true;
+}
+
+bool FileOutputHTK::InitializeInternal(const SignalBank &input) {
+  if (file_handle_ == NULL) {
+    LOG_ERROR(_T("Couldn't initialize file output. "
+                 "Please call FileOutputHTK::OpenFile first"));
+		return false;
+  }
+  if (header_written_) {
+    LOG_ERROR(_T("A header has already been written on the output file."
+                 "Please call FileOutputHTK::CloseFile to close that file, "
+                 "and FileOutputHTK::OpenFile to open an new one before "
+                 "calling FileOutputHTK::Initialize again."));
+		return false;
+  }
+  channel_count_ = input.channel_count();
+  buffer_length_ = input.buffer_length();
+  WriteHeader(channel_count_ * buffer_length_, frame_period_ms_);
+  return true;
+}
+
+void FileOutputHTK::ResetInternal() {
+  if (file_handle_ != NULL && !header_written_) {
+    WriteHeader(channel_count_ * buffer_length_, frame_period_ms_);
+  }
+}
+
+void FileOutputHTK::WriteHeader(int num_elements, float period_ms) {
+	if (header_written_)
+		return;
+
+	/* HTK format file: (taken from the HTK book - section 5.10.1)
+	 * Header: 12 bytes in total, contains:
+	 * sample_count - number of samples in file (4-byte integer)(long)
+	 * sample_period - sample period in 100ns units (4-byte integer)(long)
+	 * sample_size - number of bytes per sample (2-byte integer)(short)
+	 * parameter_kind - a code indicating the sample kind (2-byte integer)(short)
+	 */
+
+	 // To be filled in when the file is done
+	int32_t sample_count = 0;
+
+  int32_t sample_period = floor(1e4 * period_ms);
+  int16_t sample_size = num_elements * sizeof(float);
+
+  // User-defined coefficients with energy term
+	int16_t parameter_kind = H_USER + H_E;
+
+	// Fix endianness
+	sample_count = ByteSwap32(sample_count);
+	sample_period = ByteSwap32(sample_period);
+	sample_size = ByteSwap16(sample_size);
+	parameter_kind = ByteSwap16(parameter_kind);
+
+  // Enter header values. sample_count is a dummy value which is filled in on
+  // file close
+	fwrite(&sample_count, sizeof(sample_count), 1, file_handle_);
+	fwrite(&sample_period, sizeof(sample_period), 1, file_handle_);
+	fwrite(&sample_size, sizeof(sample_size), 1, file_handle_);
+	fwrite(&parameter_kind, sizeof(parameter_kind), 1, file_handle_);
+	fflush(file_handle_);
+
+	header_written_ = true;
+}
+
+
+void FileOutputHTK::Process(const SignalBank &input) {
+  if (file_handle_ == NULL) {
+    LOG_ERROR(_T("Couldn't process file output. No file is open."
+                 "Please call FileOutputHTK::OpenFile first"));
+		return;
+  }
+
+  if (!header_written_) {
+    LOG_ERROR(_T("No header has been written on the output file yet. Please"
+                 "call FileOutputHTK::Initialize() before calling "
+                 "FileOutputHTK::Process()"));
+		return;
+  }
+	float s;
+
+	for (int ch = 0; ch < input.channel_count(); ch++) {
+    for (int i = 0; i < input.buffer_length(); i++) {
+      s = input.sample(ch, i);
+      s = ByteSwapFloat(s);
+      fwrite(&s, sizeof(float), 1, file_handle_);
+    }
+  }
+	sample_count_++;
+}
+
+bool FileOutputHTK::CloseFile() {
+	if (file_handle_ == NULL)
+    return false;
+
+	// Write the first 4 bytes of the file
+	// with how many samples there are in the file
+	fflush(file_handle_);
+	rewind(file_handle_);
+	fflush(file_handle_);
+  int32_t samples = sample_count_;
+	samples = ByteSwap32(samples);
+	fwrite(&samples, sizeof(samples), 1, file_handle_);
+
+	// And close the file
+	fclose(file_handle_);
+	file_handle_ = NULL;
+	return true;
+}
+
+float FileOutputHTK::ByteSwapFloat(float d) {
+  // Endianness fix
+  float a;
+  unsigned char *dst = (unsigned char *)&a;
+  unsigned char *src = (unsigned char *)&d;
+
+  dst[0] = src[3];
+  dst[1] = src[2];
+  dst[2] = src[1];
+  dst[3] = src[0];
+
+  return a;
+}
+}  //namespace aimc
+