diff --git a/MANIFEST b/MANIFEST new file mode 100644 index 0000000..bfed780 --- /dev/null +++ b/MANIFEST @@ -0,0 +1,2 @@ +README.markdown +setup.py diff --git a/README.markdown b/README.markdown new file mode 100644 index 0000000..287fd55 --- /dev/null +++ b/README.markdown @@ -0,0 +1,9 @@ +## Introduction + +`wavebender` is an audio synthesis library for Python. + +## Usage + +``` +from wavebender import * +``` diff --git a/wavebender.py b/build/lib.linux-x86_64-2.6/wavebender/__init__.py old mode 100755 new mode 100644 similarity index 90% rename from wavebender.py rename to build/lib.linux-x86_64-2.6/wavebender/__init__.py index 840efc2..73f088a --- a/wavebender.py +++ b/build/lib.linux-x86_64-2.6/wavebender/__init__.py @@ -1,4 +1,10 @@ #!/usr/bin/env python +""" +An audio synthesis library for Python. + +It makes heavy use of the `itertools` module. +Good luck! (This is a work in progress.) +""" import sys import wave import math @@ -7,6 +13,18 @@ import argparse from itertools import * +# metadata +__author__ = 'Zach Denton' +__author_email__ = 'zacharydenton@gmail.com' +__version__ = '0.2' +__url__ = 'http://github.com/zacharydenton/wavebender' +__longdescr__ = ''' +An audio synthesis library for Python. +''' +__classifiers__ = [ + 'Topic :: Multimedia :: Sound/Audio :: Sound Synthesis' +] + def grouper(n, iterable, fillvalue=None): "grouper(3, 'ABCDEFG', 'x') --> ABC DEF Gxx" args = [iter(iterable)] * n diff --git a/wave.py b/build/lib.linux-x86_64-2.6/wavebender/wave.py similarity index 100% rename from wave.py rename to build/lib.linux-x86_64-2.6/wavebender/wave.py diff --git a/binaural.py b/examples/binaural.py similarity index 100% rename from binaural.py rename to examples/binaural.py diff --git a/damped.py b/examples/damped.py similarity index 100% rename from damped.py rename to examples/damped.py diff --git a/focus3.py b/examples/focus3.py similarity index 100% rename from focus3.py rename to examples/focus3.py diff --git a/sbagen.py b/examples/sbagen.py similarity index 100% rename from sbagen.py rename to examples/sbagen.py diff --git a/square.py b/examples/square.py similarity index 100% rename from square.py rename to examples/square.py diff --git a/violin.py b/examples/violin.py similarity index 100% rename from violin.py rename to examples/violin.py diff --git a/whitenoise.py b/examples/whitenoise.py similarity index 100% rename from whitenoise.py rename to examples/whitenoise.py diff --git a/setup.py b/setup.py new file mode 100755 index 0000000..7eb4009 --- /dev/null +++ b/setup.py @@ -0,0 +1,20 @@ +#!/usr/bin/env python +''' +Installer script for the wavebender module. +''' + +from distutils.core import setup +import wavebender + +setup ( + name = "wavebender", + description = "An audio synthesis library for Python.", + + author = wavebender.__author__, + author_email = wavebender.__author_email__, + version = wavebender.__version__, + url = wavebender.__url__, + long_description = wavebender.__longdescr__, + classifiers = wavebender.__classifiers__, + packages = ['wavebender',], +) diff --git a/wavebender/__init__.py b/wavebender/__init__.py new file mode 100644 index 0000000..73f088a --- /dev/null +++ b/wavebender/__init__.py @@ -0,0 +1,129 @@ +#!/usr/bin/env python +""" +An audio synthesis library for Python. + +It makes heavy use of the `itertools` module. +Good luck! (This is a work in progress.) +""" +import sys +import wave +import math +import struct +import random +import argparse +from itertools import * + +# metadata +__author__ = 'Zach Denton' +__author_email__ = 'zacharydenton@gmail.com' +__version__ = '0.2' +__url__ = 'http://github.com/zacharydenton/wavebender' +__longdescr__ = ''' +An audio synthesis library for Python. +''' +__classifiers__ = [ + 'Topic :: Multimedia :: Sound/Audio :: Sound Synthesis' +] + +def grouper(n, iterable, fillvalue=None): + "grouper(3, 'ABCDEFG', 'x') --> ABC DEF Gxx" + args = [iter(iterable)] * n + return izip_longest(fillvalue=fillvalue, *args) + +def sine_wave(frequency=440.0, framerate=44100, amplitude=0.5): + ''' + Generate a sine wave at a given frequency of infinite length. + ''' + period = int(framerate / frequency) + if amplitude > 1.0: amplitude = 1.0 + if amplitude < 0.0: amplitude = 0.0 + lookup_table = [float(amplitude) * math.sin(2.0*math.pi*float(frequency)*(float(i%period)/float(framerate))) for i in xrange(period)] + return (lookup_table[i%period] for i in count(0)) + +def square_wave(frequency=440.0, framerate=44100, amplitude=0.5): + for s in sine_wave(frequency, framerate, amplitude): + if s > 0: + yield amplitude + elif s < 0: + yield -amplitude + else: + yield 0.0 + +def damped_wave(frequency=440.0, framerate=44100, amplitude=0.5, length=44100): + if amplitude > 1.0: amplitude = 1.0 + if amplitude < 0.0: amplitude = 0.0 + return (math.exp(-(float(i%length)/float(framerate))) * s for i, s in enumerate(sine_wave(frequency, framerate, amplitude))) + +def white_noise(amplitude=0.5): + ''' + Generate random samples. + ''' + return (float(amplitude) * random.uniform(-1, 1) for i in count(0)) + +def compute_samples(channels, nsamples=None): + ''' + create a generator which computes the samples. + + essentially it creates a sequence of the sum of each function in the channel + at each sample in the file for each channel. + ''' + return islice(izip(*(imap(sum, izip(*channel)) for channel in channels)), nsamples) + +def write_wavefile(filename, samples, nframes=None, nchannels=2, sampwidth=2, framerate=44100, bufsize=2048): + "Write samples to a wavefile." + if nframes is None: + nframes = -1 + + w = wave.open(filename, 'w') + w.setparams((nchannels, sampwidth, framerate, nframes, 'NONE', 'not compressed')) + + max_amplitude = float(int((2 ** (sampwidth * 8)) / 2) - 1) + + # split the samples into chunks (to reduce memory consumption and improve performance) + for chunk in grouper(bufsize, samples): + frames = ''.join(''.join(struct.pack('h', int(max_amplitude * sample)) for sample in channels) for channels in chunk if channels is not None) + w.writeframesraw(frames) + + w.close() + + return filename + +def write_pcm(f, samples, sampwidth=2, framerate=44100, bufsize=2048): + "Write samples as raw PCM data." + max_amplitude = float(int((2 ** (sampwidth * 8)) / 2) - 1) + + # split the samples into chunks (to reduce memory consumption and improve performance) + for chunk in grouper(bufsize, samples): + frames = ''.join(''.join(struct.pack('h', int(max_amplitude * sample)) for sample in channels) for channels in chunk if channels is not None) + f.write(frames) + + f.close() + + return filename + +def main(): + parser = argparse.ArgumentParser() + parser.add_argument('-c', '--channels', help="Number of channels to produce", default=2, type=int) + parser.add_argument('-b', '--bits', help="Number of bits in each sample", choices=(16,), default=16, type=int) + parser.add_argument('-r', '--rate', help="Sample rate in Hz", default=44100, type=int) + parser.add_argument('-t', '--time', help="Duration of the wave in seconds.", default=60, type=int) + parser.add_argument('-a', '--amplitude', help="Amplitude of the wave on a scale of 0.0-1.0.", default=0.5, type=float) + parser.add_argument('-f', '--frequency', help="Frequency of the wave in Hz", default=440.0, type=float) + parser.add_argument('filename', help="The file to generate.") + args = parser.parse_args() + + # each channel is defined by infinite functions which are added to produce a sample. + channels = ((sine_wave(args.frequency, args.rate, args.amplitude),) for i in range(args.channels)) + + # convert the channel functions into waveforms + samples = compute_samples(channels, args.rate * args.time) + + # write the samples to a file + if args.filename == '-': + filename = sys.stdout + else: + filename = args.filename + write_wavefile(filename, samples, args.rate * args.time, args.channels, args.bits / 8, args.rate) + +if __name__ == "__main__": + main() diff --git a/wavebender/wave.py b/wavebender/wave.py new file mode 100644 index 0000000..78f54b3 --- /dev/null +++ b/wavebender/wave.py @@ -0,0 +1,503 @@ +"""Stuff to parse WAVE files. + +Usage. + +Reading WAVE files: + f = wave.open(file, 'r') +where file is either the name of a file or an open file pointer. +The open file pointer must have methods read(), seek(), and close(). +When the setpos() and rewind() methods are not used, the seek() +method is not necessary. + +This returns an instance of a class with the following public methods: + getnchannels() -- returns number of audio channels (1 for + mono, 2 for stereo) + getsampwidth() -- returns sample width in bytes + getframerate() -- returns sampling frequency + getnframes() -- returns number of audio frames + getcomptype() -- returns compression type ('NONE' for linear samples) + getcompname() -- returns human-readable version of + compression type ('not compressed' linear samples) + getparams() -- returns a tuple consisting of all of the + above in the above order + getmarkers() -- returns None (for compatibility with the + aifc module) + getmark(id) -- raises an error since the mark does not + exist (for compatibility with the aifc module) + readframes(n) -- returns at most n frames of audio + rewind() -- rewind to the beginning of the audio stream + setpos(pos) -- seek to the specified position + tell() -- return the current position + close() -- close the instance (make it unusable) +The position returned by tell() and the position given to setpos() +are compatible and have nothing to do with the actual position in the +file. +The close() method is called automatically when the class instance +is destroyed. + +Writing WAVE files: + f = wave.open(file, 'w') +where file is either the name of a file or an open file pointer. +The open file pointer must have methods write(), tell(), seek(), and +close(). + +This returns an instance of a class with the following public methods: + setnchannels(n) -- set the number of channels + setsampwidth(n) -- set the sample width + setframerate(n) -- set the frame rate + setnframes(n) -- set the number of frames + setcomptype(type, name) + -- set the compression type and the + human-readable compression type + setparams(tuple) + -- set all parameters at once + tell() -- return current position in output file + writeframesraw(data) + -- write audio frames without pathing up the + file header + writeframes(data) + -- write audio frames and patch up the file header + close() -- patch up the file header and close the + output file +You should set the parameters before the first writeframesraw or +writeframes. The total number of frames does not need to be set, +but when it is set to the correct value, the header does not have to +be patched up. +It is best to first set all parameters, perhaps possibly the +compression type, and then write audio frames using writeframesraw. +When all frames have been written, either call writeframes('') or +close() to patch up the sizes in the header. +The close() method is called automatically when the class instance +is destroyed. +""" + +import __builtin__ + +__all__ = ["open", "openfp", "Error"] + +class Error(Exception): + pass + +WAVE_FORMAT_PCM = 0x0001 + +_array_fmts = None, 'b', 'h', None, 'l' + +# Determine endian-ness +import struct +if struct.pack("h", 1) == "\000\001": + big_endian = 1 +else: + big_endian = 0 + +from chunk import Chunk + +class Wave_read: + """Variables used in this class: + + These variables are available to the user though appropriate + methods of this class: + _file -- the open file with methods read(), close(), and seek() + set through the __init__() method + _nchannels -- the number of audio channels + available through the getnchannels() method + _nframes -- the number of audio frames + available through the getnframes() method + _sampwidth -- the number of bytes per audio sample + available through the getsampwidth() method + _framerate -- the sampling frequency + available through the getframerate() method + _comptype -- the AIFF-C compression type ('NONE' if AIFF) + available through the getcomptype() method + _compname -- the human-readable AIFF-C compression type + available through the getcomptype() method + _soundpos -- the position in the audio stream + available through the tell() method, set through the + setpos() method + + These variables are used internally only: + _fmt_chunk_read -- 1 iff the FMT chunk has been read + _data_seek_needed -- 1 iff positioned correctly in audio + file for readframes() + _data_chunk -- instantiation of a chunk class for the DATA chunk + _framesize -- size of one frame in the file + """ + + def initfp(self, file): + self._convert = None + self._soundpos = 0 + self._file = Chunk(file, bigendian = 0) + if self._file.getname() != 'RIFF': + raise Error, 'file does not start with RIFF id' + if self._file.read(4) != 'WAVE': + raise Error, 'not a WAVE file' + self._fmt_chunk_read = 0 + self._data_chunk = None + while 1: + self._data_seek_needed = 1 + try: + chunk = Chunk(self._file, bigendian = 0) + except EOFError: + break + chunkname = chunk.getname() + if chunkname == 'fmt ': + self._read_fmt_chunk(chunk) + self._fmt_chunk_read = 1 + elif chunkname == 'data': + if not self._fmt_chunk_read: + raise Error, 'data chunk before fmt chunk' + self._data_chunk = chunk + self._nframes = chunk.chunksize // self._framesize + self._data_seek_needed = 0 + break + chunk.skip() + if not self._fmt_chunk_read or not self._data_chunk: + raise Error, 'fmt chunk and/or data chunk missing' + + def __init__(self, f): + self._i_opened_the_file = None + if isinstance(f, basestring): + f = __builtin__.open(f, 'rb') + self._i_opened_the_file = f + # else, assume it is an open file object already + try: + self.initfp(f) + except: + if self._i_opened_the_file: + f.close() + raise + + def __del__(self): + self.close() + # + # User visible methods. + # + def getfp(self): + return self._file + + def rewind(self): + self._data_seek_needed = 1 + self._soundpos = 0 + + def close(self): + if self._i_opened_the_file: + self._i_opened_the_file.close() + self._i_opened_the_file = None + self._file = None + + def tell(self): + return self._soundpos + + def getnchannels(self): + return self._nchannels + + def getnframes(self): + return self._nframes + + def getsampwidth(self): + return self._sampwidth + + def getframerate(self): + return self._framerate + + def getcomptype(self): + return self._comptype + + def getcompname(self): + return self._compname + + def getparams(self): + return self.getnchannels(), self.getsampwidth(), \ + self.getframerate(), self.getnframes(), \ + self.getcomptype(), self.getcompname() + + def getmarkers(self): + return None + + def getmark(self, id): + raise Error, 'no marks' + + def setpos(self, pos): + if pos < 0 or pos > self._nframes: + raise Error, 'position not in range' + self._soundpos = pos + self._data_seek_needed = 1 + + def readframes(self, nframes): + if self._data_seek_needed: + self._data_chunk.seek(0, 0) + pos = self._soundpos * self._framesize + if pos: + self._data_chunk.seek(pos, 0) + self._data_seek_needed = 0 + if nframes == 0: + return '' + if self._sampwidth > 1 and big_endian: + # unfortunately the fromfile() method does not take + # something that only looks like a file object, so + # we have to reach into the innards of the chunk object + import array + chunk = self._data_chunk + data = array.array(_array_fmts[self._sampwidth]) + nitems = nframes * self._nchannels + if nitems * self._sampwidth > chunk.chunksize - chunk.size_read: + nitems = (chunk.chunksize - chunk.size_read) / self._sampwidth + data.fromfile(chunk.file.file, nitems) + # "tell" data chunk how much was read + chunk.size_read = chunk.size_read + nitems * self._sampwidth + # do the same for the outermost chunk + chunk = chunk.file + chunk.size_read = chunk.size_read + nitems * self._sampwidth + data.byteswap() + data = data.tostring() + else: + data = self._data_chunk.read(nframes * self._framesize) + if self._convert and data: + data = self._convert(data) + self._soundpos = self._soundpos + len(data) // (self._nchannels * self._sampwidth) + return data + + # + # Internal methods. + # + + def _read_fmt_chunk(self, chunk): + wFormatTag, self._nchannels, self._framerate, dwAvgBytesPerSec, wBlockAlign = struct.unpack(' 4: + raise Error, 'bad sample width' + self._sampwidth = sampwidth + + def getsampwidth(self): + if not self._sampwidth: + raise Error, 'sample width not set' + return self._sampwidth + + def setframerate(self, framerate): + if self._datawritten: + raise Error, 'cannot change parameters after starting to write' + if framerate <= 0: + raise Error, 'bad frame rate' + self._framerate = framerate + + def getframerate(self): + if not self._framerate: + raise Error, 'frame rate not set' + return self._framerate + + def setnframes(self, nframes): + if self._datawritten: + raise Error, 'cannot change parameters after starting to write' + self._nframes = nframes + + def getnframes(self): + return self._nframeswritten + + def setcomptype(self, comptype, compname): + if self._datawritten: + raise Error, 'cannot change parameters after starting to write' + if comptype not in ('NONE',): + raise Error, 'unsupported compression type' + self._comptype = comptype + self._compname = compname + + def getcomptype(self): + return self._comptype + + def getcompname(self): + return self._compname + + def setparams(self, params): + nchannels, sampwidth, framerate, nframes, comptype, compname = params + if self._datawritten: + raise Error, 'cannot change parameters after starting to write' + self.setnchannels(nchannels) + self.setsampwidth(sampwidth) + self.setframerate(framerate) + self.setnframes(nframes) + self.setcomptype(comptype, compname) + + def getparams(self): + if not self._nchannels or not self._sampwidth or not self._framerate: + raise Error, 'not all parameters set' + return self._nchannels, self._sampwidth, self._framerate, \ + self._nframes, self._comptype, self._compname + + def setmark(self, id, pos, name): + raise Error, 'setmark() not supported' + + def getmark(self, id): + raise Error, 'no marks' + + def getmarkers(self): + return None + + def tell(self): + return self._nframeswritten + + def writeframesraw(self, data): + self._ensure_header_written(len(data)) + nframes = len(data) // (self._sampwidth * self._nchannels) + if self._convert: + data = self._convert(data) + if self._sampwidth > 1 and big_endian: + import array + data = array.array(_array_fmts[self._sampwidth], data) + data.byteswap() + data.tofile(self._file) + self._datawritten = self._datawritten + len(data) * self._sampwidth + else: + self._file.write(data) + self._datawritten = self._datawritten + len(data) + self._nframeswritten = self._nframeswritten + nframes + + def writeframes(self, data): + self.writeframesraw(data) + if self._datalength != self._datawritten: + self._patchheader() + + def close(self): + if self._file: + self._ensure_header_written(0) + if self._datalength != self._datawritten: + self._patchheader() + self._file.flush() + self._file = None + if self._i_opened_the_file: + self._i_opened_the_file.close() + self._i_opened_the_file = None + + # + # Internal methods. + # + + def _ensure_header_written(self, datasize): + if not self._datawritten: + if not self._nchannels: + raise Error, '# channels not specified' + if not self._sampwidth: + raise Error, 'sample width not specified' + if not self._framerate: + raise Error, 'sampling rate not specified' + self._write_header(datasize) + + def _write_header(self, initlength): + self._file.write('RIFF') + if not self._nframes: + self._nframes = initlength / (self._nchannels * self._sampwidth) + self._datalength = self._nframes * self._nchannels * self._sampwidth + self._form_length_pos = 4 + wave_header_format = '