python/aubionotes

   1 #!/usr/bin/python
   2
   3 def do(filein,threshold):
   4
   5     import aubio.aubioclass
   6     import aubio.median
   7     #from numarray import around
   8     from math import floor
   9     hopsize   = 512
  10     bufsize   = 4096
  11     channels  = 1
  12     frameread = 0
  13     silthres  = -80.
  14     filei     = aubio.aubioclass.sndfile(filein)
  15     srate     = filei.samplerate()
  16     myvec     = aubio.aubioclass.fvec(hopsize,channels)
  17     readsize  = filei.read(hopsize,myvec)
  18     ppick     = aubio.aubioclass.pitchpick(bufsize,hopsize,channels,myvec,srate)
  19     opick     = aubio.aubioclass.onsetpick(bufsize,hopsize,channels,myvec,threshold)
  20     mylist    = list()
  21
  22     wassilence = 0
  23     lastpitch = 0
  24     starttime = 0
  25     while(readsize==hopsize):
  26         readsize = filei.read(hopsize,myvec)
  27         val = ppick.do(myvec)
  28         midival = aubio.aubioclass.bintomidi(val,srate,bufsize)
  29         isonset,onset = opick.do(myvec)
  30         now = (frameread)*hopsize/(srate+0.)
  31         issilence = aubio.aubioclass.aubio_silence_detection(myvec.vec,silthres)
  32
  33         estmidival = 0
  34         if (issilence == 1):
  35             if (wassilence == 0):
  36                 #outputnow
  37                 endtime = (frameread-3)*hopsize/(srate+0.)
  38                 if len(mylist) > 5 :
  39                     estmidival = aubio.median.percental(mylist,len(mylist)/2)
  40                     print "sil", starttime, endtime, estmidival
  41                 #resetnow
  42                 mylist = list()
  43             else:
  44                 wassilence = 1
  45         else:
  46             if isonset == 1:
  47                 if (wassilence == 0):
  48                     #outputnow
  49                     endtime = (frameread-3)*hopsize/(srate+0.)
  50                     #estmidival = aubio.median.percental(around(array(mylist)),len(mylist)//2)
  51                     if len(mylist) > 5 :
  52                         estmidival = aubio.median.percental(mylist,len(mylist)/2)
  53                         print starttime, endtime, estmidival
  54                 #resetnow
  55                 mylist = list()
  56                 #store start time
  57                 starttime = (frameread-3)*hopsize/(srate+0.)
  58             else:
  59                 """
  60                 if(listfull):
  61                     #outputnow
  62                     endtime = (frameread-3)*hopsize/(srate+0.)
  63                     print starttime, endtime, estimmidival
  64                 else:
  65                 """
  66                 #bufferize
  67                 if midival > 50 and midival < 75:
  68                     mylist.append(floor(midival))
  69             wassilence = 0
  70
  71
  72         #elif( midival > 45 ):
  73         #    mylist.append(( now , midival+12 ))
  74         #mylist.append(toappend)
  75         frameread += 1
  76
  77
  78 if __name__ == "__main__":
  79     import sys
  80     do(sys.argv[1],sys.argv[2])