snapshot...

modules/util/viirs_surfrad.py

@@ -39,6 +39,8 @@ solzen_name = group_name_m + 'solar_zenith'
 
 params_i = [group_name_i+'temp_11_0um', group_name_i+'refl_0_65um', group_name_i+target_param]
 params_m = [group_name_m+'temp_11_0um', group_name_m+'refl_0_65um', group_name_m+target_param]
+# params_i = [group_name_i+'temp_ch38', group_name_i+'refl_ch01', group_name_i+target_param]
+# params_m = [group_name_m+'temp_ch38', group_name_m+'refl_ch01', group_name_m+target_param]
 
 param_idx_m = params_m.index(group_name_m + target_param)
 param_idx_i = params_i.index(group_name_i + target_param)
@@ -123,7 +125,7 @@ def run_all(directory, out_directory, day_night='ANY', pattern='clavrx_*.nc', st
             try:
                 num_not_missing = run(h5f, params_m, train_tiles_m, valid_tiles_m,
                                       params_i, train_tiles_i, valid_tiles_i,
-                                      num_keep_x_tiles=num_keep_x_tiles, tile_width=64, kernel_size=7, day_night=day_night)
+                                      num_keep_x_tiles=num_keep_x_tiles, tile_width=64, kernel_size=7, factor=2, day_night=day_night)
             except Exception as e:
                 print(e)
                 h5f.close()
@@ -199,7 +201,7 @@ def run_all(directory, out_directory, day_night='ANY', pattern='clavrx_*.nc', st
 #  tile_width: Must be even!
 #  kernel_size: Must be odd!
 def run(h5f, params_m, train_tiles_m, valid_tiles_m, params_i, train_tiles_i, valid_tiles_i,
-        num_keep_x_tiles=8, tile_width=64, kernel_size=3, day_night='ANY'):
+        num_keep_x_tiles=8, tile_width=64, kernel_size=3, factor=2, day_night='ANY'):
 
     border = int((kernel_size - 1)/2) + 1  # Need to add for interpolation with no edge effects
 
@@ -224,7 +226,7 @@ def run(h5f, params_m, train_tiles_m, valid_tiles_m, params_i, train_tiles_i, va
     grd_s = []
     for param in params_i:
         try:
-            grd = get_grid_values(h5f, param, 0, 0, None, num_lines*2, num_pixels*2)
+            grd = get_grid_values(h5f, param, 0, 0, None, num_lines*factor, num_pixels*factor)
             grd_s.append(grd)
         except Exception as e:
             print(e)
@@ -258,7 +260,7 @@ def run(h5f, params_m, train_tiles_m, valid_tiles_m, params_i, train_tiles_i, va
                 continue
 
             nda_m = data_m[:, j_a:j_b, i_a:i_b]
-            nda_i = data_i[:, j_a*2:j_b*2, i_a*2:i_b*2]
+            nda_i = data_i[:, j_a*factor:j_b*factor, i_a*factor:i_b*factor]
             if is_missing(param_idx_i, nda_i):
                 continue
             num_not_missing += 1