diff options
Diffstat (limited to 'examples/plotpartials.py')
| -rw-r--r-- | examples/plotpartials.py | 38 |
1 files changed, 38 insertions, 0 deletions
diff --git a/examples/plotpartials.py b/examples/plotpartials.py new file mode 100644 index 0000000..64974d5 --- /dev/null +++ b/examples/plotpartials.py @@ -0,0 +1,38 @@ +# Copyright (c) 2009 John Glover, National University of Ireland, Maynooth +# +# 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 2 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, write to the Free Software +# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA + +import simpl +import simpl.plot +import numpy as np +from scipy.io.wavfile import read + +input_file = '../test/audio/flute.wav' +audio_in_data = read(input_file) +audio_in = np.asarray(audio_in_data[1], np.float32) / 32768.0 # values between -1 and 1 +sample_rate = audio_in_data[0] + +# take just the first few frames +audio = audio_in[0:4096] +# Peak detection and partial tracking using SMS +pd = simpl.SMSPeakDetection() +pd.max_peaks = 20 +pd.hop_size = 147 +peaks = pd.find_peaks(audio) +pt = simpl.MQPartialTracking() +pt.max_partials = 20 +partials = pt.find_partials(peaks) +simpl.plot.plot_partials(partials) +simpl.plot.show() |