From 8972d98a32993ccfe70133d364e0e894cd81aa8e Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Thu, 12 Jan 2023 16:23:55 -0600
Subject: [PATCH] snapshot...

---
 modules/util/plot.py | 9 ++++++---
 1 file changed, 6 insertions(+), 3 deletions(-)

diff --git a/modules/util/plot.py b/modules/util/plot.py
index dfe6a0b4..442d2882 100644
--- a/modules/util/plot.py
+++ b/modules/util/plot.py
@@ -242,8 +242,9 @@ def make_icing_image(h5f, probs, ice_lons, ice_lats, clvrx_str_time, satellite,
     if bg_image is not None:
         bg_im = ax.imshow(bg_image, origin='upper', extent=(xmin, xmax, ymin, ymax), cmap='Greys', transform=geos)
     if probs is not None:
-        p_im = ax.imshow(probs, cmap='Purples', origin='upper', extent=(xmin, xmax, ymin, ymax),
-                         vmin=0.20, vmax=1.0, alpha=0.8, transform=geos)
+        # 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),
+                         vmin=0.55, vmax=1.0, alpha=0.8, transform=geos)
     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.BORDERS, linewidth=1.50, edgecolor='green')
@@ -251,7 +252,9 @@ def make_icing_image(h5f, probs, ice_lons, ice_lats, clvrx_str_time, satellite,
     if ice_lons is not None:
         ax.scatter(ice_lons, ice_lats, s=40.0, marker='o', color='blue', transform=ccrs.PlateCarree())
     if ice_lats_vld is not None:
-        ax.scatter(ice_lons_vld, ice_lats_vld, s=180.0, marker='o', color='red', transform=ccrs.PlateCarree())
+        # ax.scatter(ice_lons_vld, ice_lats_vld, s=180.0, marker='o', color='red', transform=ccrs.PlateCarree())
+        ax.scatter(ice_lons_vld, ice_lats_vld, s=500.0, marker='o', color='white', transform=ccrs.PlateCarree())
+        ax.scatter(ice_lons_vld, ice_lats_vld, s=200.0, marker='o', color='black', transform=ccrs.PlateCarree())
 
     title = satellite+', '+clvrx_str_time+'  >50% Probability Icing'
     plt.title(title, loc='left', fontweight='bold', fontsize=24)
-- 
GitLab