diff --git a/modules/util/viirs_surfrad.py b/modules/util/viirs_surfrad.py index c1ebe11ec49ffa9a8edab8f963be05c1e3a1bfd2..78c522eb80df3d4dfe54b2681389266971f9d6bd 100644 --- a/modules/util/viirs_surfrad.py +++ b/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