Added more documentation

testing/correlation_test.py

@@ -14,17 +14,21 @@ def main(path):
     dmv = radiance.get_radiance_from_rnc(path)
     hk = housekeeping.get_all_housekeeping(path)
 
+    #only use dmv data from scene H
     data = pd.DataFrame(index=hk.index, columns=('dmv', 'hk'))
     data['dmv'] = dmv.xs(str.encode('H'), level='scene').mean(axis=1)
     data['hk'] = hk['SCEtemp']
 
+    #get rid of hk rows where scene isn't H
     data.dropna(inplace=True)
 
     try:
+        #subtract the rolling average
         for x in data.columns:
             tmp = scipy.signal.savgol_filter(data[x].values[:], 201, 3)
             data[x] = data[x] - tmp
 
+        #correlate the data
         correlation = np.correlate(data['dmv'].values[:], data['hk'].values[:], mode='same')
 
         print('max = ', np.amax(correlation), ' : min = ', np.amin(correlation))
@@ -44,6 +48,7 @@ def main(path):
         plt.savefig(name)
         plt.clf()
 
+    #probably the scipy.signal.savgol_filter failing
     except:
         print('FAIL')