Skip to content
Snippets Groups Projects
Commit 700d357b authored by tomrink's avatar tomrink
Browse files

snapshot...

parent cc4263fc
No related branches found
No related tags found
No related merge requests found
...@@ -1581,7 +1581,8 @@ def run_icing_predict_image_1x1(clvrx_dir='/Users/tomrink/data/clavrx/RadC/', ou ...@@ -1581,7 +1581,8 @@ def run_icing_predict_image_1x1(clvrx_dir='/Users/tomrink/data/clavrx/RadC/', ou
varX_std = np.where(np.isnan(varX_std), 0, varX_std) varX_std = np.where(np.isnan(varX_std), 0, varX_std)
probs = nght_model.predict_proba(varX_std)[:, 1] probs = nght_model.predict_proba(varX_std)[:, 1]
fd_probs[np.invert(day_idxs)] = probs[np.invert(day_idxs)] fd_probs[np.invert(day_idxs)] = probs[np.invert(day_idxs)]
fd_probs[clr_idxs] = 0.0
fd_probs[clr_idxs] = -1.0
max_prob = fd_probs.reshape((num_lines, num_elems)) max_prob = fd_probs.reshape((num_lines, num_elems))
... ...
......
...@@ -277,7 +277,7 @@ def make_icing_image(h5f, probs, ice_lons, ice_lats, clvrx_str_time, satellite, ...@@ -277,7 +277,7 @@ def make_icing_image(h5f, probs, ice_lons, ice_lats, clvrx_str_time, satellite,
# p_im = ax.imshow(probs, cmap='Purples', origin='upper', extent=(xmin, xmax, ymin, ymax), # p_im = ax.imshow(probs, cmap='Purples', origin='upper', extent=(xmin, xmax, ymin, ymax),
p_im = ax.imshow(probs, cmap='jet', origin='upper', extent=(xmin, xmax, ymin, ymax), p_im = ax.imshow(probs, cmap='jet', origin='upper', extent=(xmin, xmax, ymin, ymax),
# p_im = ax.imshow(probs, cmap=icing_cmap, origin='upper', extent=(xmin, xmax, ymin, ymax), # p_im = ax.imshow(probs, cmap=icing_cmap, origin='upper', extent=(xmin, xmax, ymin, ymax),
vmin=0.50, vmax=1.0, alpha=0.8, transform=geos) vmin=0.55, vmax=1.0, alpha=0.8, transform=geos)
ax.coastlines(resolution='50m', color='green', linewidth=1.50) ax.coastlines(resolution='50m', color='green', linewidth=1.50)
ax.add_feature(ccrs.cartopy.feature.STATES, linewidth=1.50, edgecolor='green') ax.add_feature(ccrs.cartopy.feature.STATES, linewidth=1.50, edgecolor='green')
ax.add_feature(ccrs.cartopy.feature.BORDERS, linewidth=1.50, edgecolor='green') ax.add_feature(ccrs.cartopy.feature.BORDERS, linewidth=1.50, edgecolor='green')
... ...
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please to comment