added air_interferometer_check to quick_vis.py and fixed an issue so...

added air_interferometer_check to quick_vis.py and fixed an issue so qc_percent is 1 when missing_data_flag_check is 1; make minor adjustments to detectory_temp.py and testing_bst.py

testing/detector_temp_flag.py

@@ -159,7 +159,7 @@ def main(path):
     plt.plot(dmv['Time'], std_data['std_mean_std'])
 
     print('title: ', title)
-    plt.savefig("sgp_detect_pngs/detector_temp_flag_test_%s.png"%title,bbox_inches='tight')
+    # plt.savefig("sgp_detect_pngs/detector_temp_flag_test_%s.png"%title,bbox_inches='tight')
     # plt.savefig("detector_temp_flag_test_%s_2.png"%title,bbox_inches='tight')
     # plt.show()
     plt.clf()

testing/quick_vis.py

@@ -20,7 +20,8 @@ def main(filename):
         'safing_check', 'encoder_check', 'detector_check',
         'hbb_thermistor_check', 'abb_thermistor_check', 'spike_check',
         'hbb_temp_outlier_check', 'abb_temp_outlier_check',
-        'bst_temp_outlier_check', 'sce_temp_deviation_check',
+        'bst_temp_outlier_check', 'air_interferometer_outlier_check',
+        'sce_temp_deviation_check',
         'hbb_stable_check', 'hbb_covariance_check', 'imaginary_radiance_check',
         'detector_temp_check',
         'sky_brightness_temp_spectral_averages_ch1_check',
@@ -56,6 +57,7 @@ def main(filename):
     #infers times for NaT values
     data.loc[data['time'] < 0] = 0
     data['missing_data_flag_check'].loc[data['time'] == 0] = 1
+    data['qc_percent'].loc[data['time'] == 0] = 1
     for t in range(len(data['time'])):
         if data['time'].iloc[t] == 0 and t != 0:
             data['time'].iloc[t] = data['time'].iloc[t-1] + time_increment
@@ -90,6 +92,8 @@ def main(filename):
                 var_name = 'surface_bt_ch1' + ' '*30
             elif 'ch2' in var_name:
                 var_name = 'surface_bt_ch2' + ' '*30
+        if 'air_interferometer_outlier' in var_name:
+            var_name = 'air_inter_outlier' + ' '*30
         if 'qc_percent' in var_name:
             var_name = 'qc_percent ({:3.2f}%)'.format(qc_percent_num) + ' '*30
 
@@ -127,10 +131,10 @@ def main(filename):
 
     print('saving...')
     end_of_name = filename.split('/')[-1].split('.')[0] + '.png'
-    # plt.savefig('/Users/adiebold/aeri_quality_control/testing/pngs/awr/' + end_of_name)
+    plt.savefig('/Users/adiebold/aeri_quality_control/testing/pngs/sgp/' + end_of_name)
     # plt.savefig('/Users/adiebold/aeri_quality_control/testing/' + end_of_name)
     #comment out plt.show() when doing a directory
-    plt.show()
+    # plt.show()
     plt.clf()
 
     print('finished')
@@ -147,7 +151,7 @@ if __name__ == '__main__':
     print(args.filepath)
 
     #amount of files to skip
-    skip_num = 0
+    skip_num = 1095    #sgp at 245
     curr_num = 0
     print('skip_num = ', skip_num, '\n')
     if os.path.isdir(args.filepath):

testing/testing_bst.py

@@ -13,7 +13,9 @@ def main(path, bounds = None):
     data = housekeeping.get_all_housekeeping(path)
     if bounds != None:
         data = data[(data['Time'] > bounds[0]) & (data['Time'] < bounds[1])]
-    variables = ['BBsupportStructureTemp', 'HBBapexTemp', 'HBBbottomTemp', 'HBBtopTemp', 'ABBapexTemp', 'ABBbottomTemp', 'ABBtopTemp', 'airNearBBsTemp', 'SCEtemp', 'sceneMirrorPosition']
+    variables = ['BBsupportStructureTemp', 'HBBapexTemp', 'HBBbottomTemp', 
+                'HBBtopTemp', 'ABBapexTemp', 'ABBbottomTemp', 'ABBtopTemp',
+                'airNearBBsTemp', 'SCEtemp', 'sceneMirrorPosition']
     fig, ax = plt.subplots(len(variables),figsize=(25,15),sharex=True)
     for x,v in enumerate(variables):
         name = v + ' '*30