Changeset a0fd4e4 for examples/aubionotes.c

Timestamp:

Apr 2, 2005, 3:29:23 PM (19 years ago)

Author:

Paul Brossier <piem@altern.org>

Branches:

feature/autosink, feature/cnn, feature/cnn_org, feature/constantq, feature/crepe, feature/crepe_org, feature/pitchshift, feature/pydocstrings, feature/timestretch, fix/ffmpeg5, master, pitchshift, sampler, timestretch, yinfft+

Children:

cf7c76a

Parents:

e50b695

Message:

merge aubionotes and aubionotesmedian

add make_audio_plot to gnuplot.py
adds autocorrelation to mathutils
adds Makefile.in to junk files

File:

• examples/aubionotes.c (modified) (10 diffs)

examples/aubionotes.c

@@ -22 +22 @@
 #include <getopt.h>
 #include <unistd.h>
-#include <math.h> // how do i do a floorf with a mask again ?
+#include <math.h>
 #include "aubio.h"
 #include "aubioext.h"
@@ -35 +35 @@
 int usedoubled = 1;
 int usemidi = 1;
+int median = 0;
 
 /* energy,specdiff,hfc,complexdomain,phase */
 aubio_onsetdetection_type type_onset  = hfc;
-aubio_onsetdetection_type type_onset2 = hfc;
+aubio_onsetdetection_type type_onset2 = complexdomain;
 smpl_t threshold                      = 0.3;
 smpl_t threshold2                     = -90.;
@@ -61 +62 @@
 
 /* pitch objects */
-int curnote                = 0;
 smpl_t pitch               = 0.;
 aubio_pitchdetection_t * pitchdet;
-int bufpos = 0;
-smpl_t curlevel = 0;
+
+fvec_t * note_buffer = NULL;
+fvec_t * note_buffer2 = NULL;
+uint_t medianfiltlen = 6;
+smpl_t curlevel = 0.;
 smpl_t maxonset = 0.;
 
@@ -101 +104 @@
 }
 
+
+/** append new note candidate to the note_buffer and return filtered value. we
+ * need to copy the input array as vec_median destroy its input data.*/
+
+void note_append(fvec_t * note_buffer, smpl_t curnote); 
+void note_append(fvec_t * note_buffer, smpl_t curnote) {
+  uint_t i = 0;
+  for (i = 0; i < note_buffer->length - 1; i++) { 
+      note_buffer->data[0][i] = note_buffer->data[0][i+1];
+  }
+  note_buffer->data[0][note_buffer->length - 1] = curnote;
+  return;
+}
+
+uint_t get_note(fvec_t *note_buffer, fvec_t *note_buffer2);
+uint_t get_note(fvec_t *note_buffer, fvec_t *note_buffer2){
+  uint_t i = 0;
+  for (i = 0; i < note_buffer->length; i++) { 
+      note_buffer2->data[0][i] = note_buffer->data[0][i];
+  }
+  return vec_median(note_buffer2);
+}
+
+
+smpl_t curnote = 0.;
+smpl_t newnote = 0.;
+uint_t isready = 0;
+
 int aubio_process(float **input, float **output, int nframes);
 int aubio_process(float **input, float **output, int nframes) {
@@ -126 +157 @@
       
       pitch = aubio_pitchdetection(pitchdet,ibuf);
+      if(median){
+      newnote = (FLOOR)(bintomidi(pitch,samplerate,buffer_size*4)+.5);
+      note_append(note_buffer, newnote);
+      }
 
       /* curlevel is negatif or 1 if silence */
@@ -133 +168 @@
         if (curlevel == 1.) {
           isonset=0;
+          if (median) isready = 0;
           /* send note off */
           send_noteon(mdriver,curnote,0);
         } else {
+          if (median) {
+                  isready = 1;
+          } else {
           /* kill old note */
           send_noteon(mdriver,curnote,0);
@@ -147 +186 @@
             send_noteon(mdriver,curnote,127+(int)FLOOR(curlevel));
           }
+          }
+          
           for (pos = 0; pos < overlap_size; pos++){
             obuf->data[0][pos] = woodblock->data[0][pos];
@@ -152 +193 @@
         }
       } else {
+              if (median) {
+//        if (curlevel != 1) {
+//          if (bufpos == note_buffer->length)
+//            curnote = aubio_getnote(note_buffer);
+//        }
+//
+        if (isready > 0)
+            isready++;
+        if (isready == note_buffer->length)
+        {
+            //outmsg("%d, %f\n", isready, curnote);
+            send_noteon(mdriver,curnote,0);
+        /* kill old note */
+            newnote = get_note(note_buffer, note_buffer2);
+            curnote = newnote;
+        /* get and send new one */
+        /*if (curnote<45){
+          send_noteon(mdriver,curnote,0);
+        } else {*/
+            if (curnote>45){
+                send_noteon(mdriver,curnote,127+(int)FLOOR(curlevel));
+            }
+        }
+        } // if median
         for (pos = 0; pos < overlap_size; pos++)
           obuf->data[0][pos] = 0.;
@@ -188 +253 @@
   fftgrain    = new_cvec(buffer_size, channels);
 
-  pitchdet = new_aubio_pitchdetection(buffer_size*4, overlap_size, channels, samplerate, mcomb, freq);
-  //pitchdet = new_aubio_pitchdetection(buffer_size*2, overlap_size, channels, samplerate, yin, freq);
+  aubio_pitchdetection_type mode = yin; // mcomb
+  pitchdet = new_aubio_pitchdetection(buffer_size*4, overlap_size, channels, samplerate, mode, freq);
   
+if (median) {
+  note_buffer = new_fvec(medianfiltlen, 1);
+  note_buffer2= new_fvec(medianfiltlen, 1);
+}
   /* read the output sound once */
   file_read(onsetfile, overlap_size, woodblock);
@@ -262 +331 @@
   del_aubio_pitchdetection(pitchdet);
   del_fvec(onset);
+  if (median) {
+  del_fvec(note_buffer);
+  del_fvec(note_buffer2);
+  }
 
   debug("End of program.\n");