[4cc9fe5] | 1 | #! /usr/bin/python |
---|
| 2 | |
---|
| 3 | from aubio.bench.node import * |
---|
[0029638] | 4 | from aubio.tasks import * |
---|
| 5 | |
---|
| 6 | class benchpitch(bench): |
---|
| 7 | |
---|
[ee7b9da] | 8 | """ list of values to store per file """ |
---|
| 9 | valuenames = ['mode'] |
---|
| 10 | """ list of lists to store per file """ |
---|
| 11 | valuelists = ['orig', 'mean', 'med'] |
---|
| 12 | """ list of values to print per dir """ |
---|
| 13 | printnames = [ 'mode'] |
---|
| 14 | |
---|
| 15 | """ per dir """ |
---|
| 16 | formats = {'mode': "%12s" , 'thres': "%5.4s", |
---|
| 17 | 'dist': "%5.4s", 'prec': "%5.4s", 'recl': "%5.4s", |
---|
| 18 | 'Ttrue': "%5.4s", 'Tfp': "%5.4s", 'Tfn': "%5.4s", |
---|
| 19 | 'Tm': "%5.4s", 'Td': "%5.4s", |
---|
| 20 | 'aTtrue':"%5.4s", 'aTfp': "%5.4s", 'aTfn': "%5.4s", |
---|
| 21 | 'aTm': "%5.4s", 'aTd': "%5.4s", |
---|
| 22 | 'mean': "%5.40s", 'smean': "%5.40s", |
---|
| 23 | 'amean': "%5.40s", 'samean': "%5.40s"} |
---|
| 24 | |
---|
| 25 | def dir_eval(self): |
---|
| 26 | """ evaluate statistical data over the directory """ |
---|
| 27 | v = self.v |
---|
| 28 | |
---|
| 29 | v['mode'] = self.params.pitchmode |
---|
| 30 | |
---|
[4f4a8a4] | 31 | def file_exec(self,input,output): |
---|
[0029638] | 32 | filetask = self.task(input,params=self.params) |
---|
| 33 | computed_data = filetask.compute_all() |
---|
[ee7b9da] | 34 | orig,mean,med = filetask.eval(computed_data) |
---|
| 35 | |
---|
| 36 | self.v['orig'].append(orig) |
---|
| 37 | self.v['mean'].append(mean) |
---|
| 38 | self.v['med'].append(med) |
---|
| 39 | #print results#, computed_data |
---|
[0029638] | 40 | #print input, results, results - float(input.split('.')[-2]) |
---|
| 41 | |
---|
[4f4a8a4] | 42 | def run_bench(self,modes=['schmitt']): |
---|
[ee7b9da] | 43 | from os.path import basename |
---|
| 44 | d = [] |
---|
[0029638] | 45 | self.modes = modes |
---|
[ee7b9da] | 46 | self.pretty_titles() |
---|
[0029638] | 47 | for mode in self.modes: |
---|
[4f4a8a4] | 48 | self.params.pitchmode = mode |
---|
[ee7b9da] | 49 | self.dir_eval_print() |
---|
| 50 | self.plotpitchtessiture(d, |
---|
| 51 | self.v['orig'], |
---|
| 52 | self.v['med'], |
---|
| 53 | plottitle=self.v['mode'], |
---|
| 54 | plotmode='points') |
---|
| 55 | #d.append('beta = .25,orig(x) title \"-2 octave\"') |
---|
| 56 | d.append('beta = .50,orig(x) title \"-1 octave\"') |
---|
| 57 | d.append('beta = 1.0,orig(x) title \"original\"') |
---|
| 58 | d.append('beta = 2.0,orig(x) title \"+1 octave\"') |
---|
| 59 | title = basename(self.datadir) |
---|
| 60 | outplot = "_-_".join(('pitchtessiture',title)) |
---|
| 61 | for ext in ('ps','png','svg',''): |
---|
| 62 | self.plotplotpitchtessiture(d, |
---|
| 63 | plottitle=title, |
---|
| 64 | outplot=outplot, |
---|
| 65 | extension=ext) |
---|
| 66 | |
---|
| 67 | """ |
---|
| 68 | Plot functions |
---|
| 69 | """ |
---|
| 70 | |
---|
| 71 | def plotpitchtessiture(self,d,lx,ly,plottitle="",plotmode='linespoints'): |
---|
| 72 | import Gnuplot, Gnuplot.funcutils |
---|
| 73 | d.append(Gnuplot.Data(lx, ly, with=plotmode, title="%s" % (plottitle) )) |
---|
| 74 | |
---|
| 75 | def plotplotpitchtessiture(self,d,plottitle='',outplot=0,extension=''): |
---|
| 76 | from aubio.gnuplot import gnuplot_create |
---|
| 77 | g = gnuplot_create(outplot=outplot,extension=extension) |
---|
| 78 | g.title(plottitle) |
---|
| 79 | g('orig(x) = beta*x') |
---|
| 80 | g.xlabel('original pitch (Hz)') |
---|
| 81 | g.ylabel('detected pitch (Hz)') |
---|
| 82 | g('set log xy') |
---|
| 83 | g.plot(*d) |
---|
| 84 | |
---|
[0029638] | 85 | |
---|
| 86 | if __name__ == "__main__": |
---|
| 87 | import sys |
---|
| 88 | if len(sys.argv) > 1: datapath = sys.argv[1] |
---|
| 89 | else: print "error: a path is required"; sys.exit(1) |
---|
[4f4a8a4] | 90 | if len(sys.argv) > 2: |
---|
| 91 | for each in sys.argv[3:-1]: print each |
---|
[ee7b9da] | 92 | modes = ['schmitt', 'yin', 'mcomb', 'fcomb'] |
---|
| 93 | #modes = ['fcomb'] |
---|
[0029638] | 94 | |
---|
[ee7b9da] | 95 | params = taskparams() |
---|
| 96 | params.bufsize = 2048 |
---|
| 97 | params.hopsize = params.bufsize/2 |
---|
| 98 | benchpitch = benchpitch(datapath,params=params) |
---|
[0029638] | 99 | benchpitch.task = taskpitch |
---|
| 100 | |
---|
[ee7b9da] | 101 | #benchpitch.titles = [ 'mode', 'thres', 'avg', 'avgdist' ] |
---|
| 102 | #benchpitch.formats = ["%12s" , "| %6s", "| %6s", "| %6s", "| %6s", "| %6s" ] |
---|
[0029638] | 103 | try: |
---|
| 104 | benchpitch.run_bench(modes=modes) |
---|
| 105 | except KeyboardInterrupt: |
---|
| 106 | print "Interrupted by user" |
---|
| 107 | sys.exit(1) |
---|
| 108 | |
---|
| 109 | sys.exit(0) |
---|