Tweaks to resample()

Author: eggsyntax

datamapper.py

@@ -98,10 +98,9 @@ def sample_rate(self):
     def sample_rate(self, r):
         self._sample_rate = r
 
-	def resample(self, factor):
-		for dob in self:
-			print 'DOB:', dob
-			dob.resample(factor)
+    def resample(self, factor):
+        for dob in self:
+            dob.resample(factor)
 
     def append(self, dob):
         assert isinstance(dob, DataObject), "You've got a " + str(type(dob)) + ", not a DataObject"
@@ -256,9 +255,9 @@ def trim_to(self, length):
         for ts in self.values():
             ts.replace_data(ts[:length])
 
-	def resample(self, factor):
-		for ts in self.values():
-			ts.resample(factor)
+    def resample(self, factor):
+        for ts in self.values():
+            ts.resample(factor)
 
     def __setitem__(self, key, ts):
          assert(isinstance(ts, TimeSeries))

geothermophone.py

@@ -0,0 +1,51 @@
+''' This is a parser for a particular project (the geothermophone). The data
+is already in optimized form (a list of dicts from varname to a list of values).
+'''
+from datamapper import DataObjectCollection, DataObject, TimeSeries
+from dataparser import DataParser
+from renderers.midirenderers import MidiCCRenderer
+from renderers.visual_renderers import LineGraphRenderer
+
+
+class GriddedDataParser(DataParser):
+    ''' Convert the input data into a DataObjectCollection and return it. '''
+
+    def __init__(self):
+
+        self.file = open(
+            '/Users/egg/Documents/Work In Progress/geothermophone/final_values/final_values.txt')
+
+    def parse(self):
+        d = eval(self.file.readline().strip())
+        doc = DataObjectCollection(sample_rate=1 / 3.0)
+        for i, octant in enumerate(d):
+            do = DataObject()
+            for varname, values in octant.items():
+                ts = TimeSeries(values)
+                do[varname] = ts
+            doc.append(do)
+        return doc
+
+
+def build_data():
+    parser = GriddedDataParser()
+    doc = parser.parse()
+    doc.resample(.1403) # Trial and error to get 30 min of data
+
+    # renderer = LineGraphRenderer()
+    # plot = renderer.render(doc, showplot=True, render_separate=True)
+
+    # Produce MIDI data
+    renderer = MidiCCRenderer()
+    # divide 20 by 1.782
+    renderer.tempo = 20
+    var_to_midi_map = {'air': 74,
+                       'prate': 75,
+                       'rhum': 76,
+                       'wspd': 77}
+    transformed_doc = doc.transform(var_to_midi_map, renderer)
+    renderer.render(transformed_doc,
+                    output_file='/Users/egg/Documents/Work In Progress/geothermophone/midi/geothermophone.mid')
+
+build_data()
+

renderers/midirenderers.py

@@ -56,8 +56,6 @@ def render(self, doc, output_file='output.mid'):
                 for i, val in enumerate(time_series):
                     time = float(i) / time_series.sample_rate
                     logging.debug(str(time) + ' ' + str(val))
-                    #MyMIDI.addNote(track,channel,pitch,time,duration,volume)
-                    #print 'adding: {} {} {} {} {}'.format(track, channel+1, time*8, cc_number, val)
                     MyMIDI.addControllerEvent(track, channel, time, cc_number, int(val))
 
         # And write it to disk.