def load(filename):
"""
Load a wave file and return the signal, sample rate and number of channels.
Can be any format supported by the underlying library (libsndfile or SciPy)
"""
if wav_loader == 'pysoundfile':
sf = SoundFile(filename)
signal = sf.read()
channels = sf.channels
sample_rate = sf.samplerate
samples = len(sf)
file_format = sf.format_info + ' ' + sf.subtype_info
sf.close()
elif wav_loader == 'scikits.audiolab':
sf = Sndfile(filename, 'r')
signal = sf.read_frames(sf.nframes)
channels = sf.channels
sample_rate = sf.samplerate
samples = sf.nframes
file_format = sf.format
sf.close()
elif wav_loader == 'scipy.io.wavfile':
sample_rate, signal = read(filename)
try:
channels = signal.shape[1]
except IndexError:
channels = 1
samples = signal.shape[0]
file_format = str(signal.dtype)
return signal, sample_rate, channels
def load_dict(filename):
"""
Load a wave file and return the signal, sample rate and number of channels.
Can be any format supported by the underlying library (libsndfile or SciPy)
"""
soundfile = {}
if wav_loader == 'pysoundfile':
sf = SoundFile(filename)
soundfile['signal'] = sf.read()
soundfile['channels'] = sf.channels
soundfile['fs'] = sf.samplerate
soundfile['samples'] = len(sf)
soundfile['format'] = sf.format_info + ' ' + sf.subtype_info
sf.close()
elif wav_loader == 'scikits.audiolab':
sf = Sndfile(filename, 'r')
soundfile['signal'] = sf.read_frames(sf.nframes)
soundfile['channels'] = sf.channels
soundfile['fs'] = sf.samplerate
soundfile['samples'] = sf.nframes
soundfile['format'] = sf.format
sf.close()
elif wav_loader == 'scipy.io.wavfile':
soundfile['fs'], soundfile['signal'] = read(filename)
try:
soundfile['channels'] = soundfile['signal'].shape[1]
except IndexError:
soundfile['channels'] = 1
soundfile['samples'] = soundfile['signal'].shape[0]
soundfile['format'] = str(soundfile['signal'].dtype)
return soundfile
def load(filename):
"""
Load a wave file and return the signal, sample rate and number of channels.
Can be any format supported by the underlying library (libsndfile or SciPy)
"""
if wav_loader == 'pysoundfile':
sf = SoundFile(filename)
signal = sf.read()
channels = sf.channels
sample_rate = sf.samplerate
sf.close()
elif wav_loader == 'scikits.audiolab':
sf = Sndfile(filename, 'r')
signal = sf.read_frames(sf.nframes)
channels = sf.channels
sample_rate = sf.samplerate
sf.close()
elif wav_loader == 'scipy.io.wavfile':
sample_rate, signal = read(filename)
try:
channels = signal.shape[1]
except IndexError:
channels = 1
return signal, sample_rate, channels
def analyze(filename):
if wav_loader == 'pysoundfile':
sf = SoundFile(filename)
signal = sf.read()
channels = sf.channels
sample_rate = sf.samplerate
samples = len(sf)
file_format = sf.format_info + ' ' + sf.subtype_info
sf.close()
elif wav_loader == 'audiolab':
sf = Sndfile(filename, 'r')
signal = sf.read_frames(sf.nframes)
channels = sf.channels
sample_rate = sf.samplerate
samples = sf.nframes
file_format = sf.format
sf.close()
elif wav_loader == 'scipy':
sample_rate, signal = read(filename)
try:
channels = signal.shape[1]
except IndexError:
channels = 1
samples = signal.shape[0]
file_format = str(signal.dtype)
# Scale common formats
# Other bit depths (24, 20) are not handled by SciPy correctly.
if file_format == 'int16':
signal = signal.astype(float) / (2**15)
elif file_format == 'uint8':
signal = (signal.astype(float) - 128) / (2**7)
elif file_format == 'int32':
signal = signal.astype(float) / (2**31)
elif file_format == 'float32':
pass
else:
raise Exception("Don't know how to handle file "
"format {}".format(file_format))
else:
raise Exception("wav_loader has failed")
header = 'dBFS values are relative to a full-scale square wave'
if samples / sample_rate >= 1:
length = str(samples / sample_rate) + ' seconds'
else:
length = str(samples / sample_rate * 1000) + ' milliseconds'
results = [
'Properties for "' + filename + '"',
str(file_format),
'Channels:\t%d' % channels,
'Sampling rate:\t%d Hz' % sample_rate,
'Samples:\t%d' % samples,
'Length: \t' + length,
'-----------------',
]
if channels == 1:
# Monaural
results += properties(signal, sample_rate)
elif channels == 2:
# Stereo
if array_equal(signal[:, 0], signal[:, 1]):
results += ['Left and Right channels are identical:']
results += properties(signal[:, 0], sample_rate)
else:
results += ['Left channel:']
results += properties(signal[:, 0], sample_rate)
results += ['Right channel:']
results += properties(signal[:, 1], sample_rate)
else:
# Multi-channel
for ch_no, channel in enumerate(signal.transpose()):
results += ['Channel %d:' % (ch_no + 1)]
results += properties(channel, sample_rate)
display(header, results)
plot_histogram = False
if plot_histogram:
histogram(signal)
def analyze(filename):
if wav_loader == 'pysoundfile':
sf = SoundFile(filename)
signal = sf.read()
channels = sf.channels
sample_rate = sf.samplerate
samples = len(sf)
file_format = sf.format_info + ' ' + sf.subtype_info
sf.close()
elif wav_loader == 'scikits.audiolab':
sf = Sndfile(filename, 'r')
signal = sf.read_frames(sf.nframes)
channels = sf.channels
sample_rate = sf.samplerate
samples = sf.nframes
file_format = sf.format
sf.close()
elif wav_loader == 'scipy.io.wavfile':
sample_rate, signal = read(filename)
try:
channels = signal.shape[1]
except IndexError:
channels = 1
samples = signal.shape[0]
file_format = str(signal.dtype)
# Scale common formats
# Other bit depths (24, 20) are not handled by SciPy correctly.
if file_format == 'int16':
signal = signal.astype(float) / (2**15)
elif file_format == 'uint8':
signal = (signal.astype(float) - 128) / (2**7)
elif file_format == 'int32':
signal = signal.astype(float) / (2**31)
elif file_format == 'float32':
pass
else:
raise Exception("Don't know how to handle file "
"format {}".format(file_format))
else:
raise Exception("wav_loader has failed")
header = 'dBFS values are relative to a full-scale square wave'
if samples/sample_rate >= 1:
length = str(samples/sample_rate) + ' seconds'
else:
length = str(samples/sample_rate*1000) + ' milliseconds'
results = [
"Using sound file backend '" + wav_loader + "'",
'Properties for "' + filename + '"',
str(file_format),
'Channels:\t%d' % channels,
'Sampling rate:\t%d Hz' % sample_rate,
'Samples:\t%d' % samples,
'Length: \t' + length,
'-----------------',
]
if channels == 1:
# Monaural
results += properties(signal, sample_rate)
elif channels == 2:
# Stereo
if array_equal(signal[:, 0], signal[:, 1]):
results += ['Left and Right channels are identical:']
results += properties(signal[:, 0], sample_rate)
else:
results += ['Left channel:']
results += properties(signal[:, 0], sample_rate)
results += ['Right channel:']
results += properties(signal[:, 1], sample_rate)
else:
# Multi-channel
for ch_no, channel in enumerate(signal.transpose()):
results += ['Channel %d:' % (ch_no + 1)]
results += properties(channel, sample_rate)
display(header, results)
plot_histogram = False
if plot_histogram:
histogram(signal)