Skip to content
Snippets Groups Projects
Commit 665157a1 authored by tomrink's avatar tomrink
Browse files

fix stupid bug

parent e8798f2f
No related branches found
No related tags found
No related merge requests found
......@@ -1088,7 +1088,7 @@ def fov_extract(outfile='/home/rink/fovs_l1b_out.h5', train_params=l1b_ds_list,
def tile_extract(trnfile='/home/rink/tiles_l1b_train.h5', tstfile='/home/rink/tiles_l1b_test.h5',
train_params=l1b_ds_list, ds_types=l1b_ds_types, cld_temp_name='temp_11_0um_nom', augment=False, split=0.2):
train_params=l1b_ds_list, ds_types=l1b_ds_types, cld_mask_name='cloud_mask', augment=False, split=0.2):
icing_int_s = []
ice_time_s = []
no_ice_time_s = []
......@@ -1116,10 +1116,12 @@ def tile_extract(trnfile='/home/rink/tiles_l1b_train.h5', tstfile='/home/rink/ti
icing_int = f['icing_intensity'][:]
for i in range(num_obs):
cld_tmp = f[cld_temp_name][i, 12:28, 12:28]
cld_msk = f[cld_mask_name][i, 12:28, 12:28]
for ds_name in train_params:
dat = f[ds_name][i, 12:28, 12:28]
#np.where(cld_tmp < 273.0, dat, np.nan)
# np.where(cld_tmp < 273.0, dat, np.nan)
#keep = np.logical_or(cld_msk[i,] == 2, cld_msk[i,] == 3) # cloudy
#np.where(keep, dat, np.nan)
icing_data_dct[ds_name].append(dat)
icing_int_s.append(icing_int[i])
......@@ -1151,10 +1153,12 @@ def tile_extract(trnfile='/home/rink/tiles_l1b_train.h5', tstfile='/home/rink/ti
lats = f['latitude']
for i in range(num_obs):
cld_tmp = f[cld_temp_name][i, 12:28, 12:28]
cld_msk = f[cld_mask_name][i, 12:28, 12:28]
for ds_name in train_params:
dat = f[ds_name][i, 12:28, 12:28]
#np.where(cld_tmp < 273.0, dat, np.nan)
# np.where(cld_tmp < 273.0, dat, np.nan)
#keep = np.logical_or(cld_msk[i,] == 2, cld_msk[i,] == 3) # cloudy
#np.where(keep, dat, np.nan)
no_icing_data_dct[ds_name].append(dat)
num_no_ice += 1
no_ice_time_s.append(times[i])
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment