Skip to content
Snippets Groups Projects
Commit 34f42b0b authored by tomrink's avatar tomrink
Browse files

snapshot...

parent ecc1a315
No related branches found
No related tags found
No related merge requests found
...@@ -14,6 +14,7 @@ import h5py ...@@ -14,6 +14,7 @@ import h5py
from util.util import get_grid_values_all, get_cartopy_crs, homedir from util.util import get_grid_values_all, get_cartopy_crs, homedir
import csv import csv
from matplotlib.colors import LinearSegmentedColormap from matplotlib.colors import LinearSegmentedColormap
from matplotlib.markers import MarkerStyle
from util.setup import home_dir from util.setup import home_dir
import mpl_scatter_density import mpl_scatter_density
...@@ -255,6 +256,7 @@ def plot_image3(image, bt, cmap='Greys'): ...@@ -255,6 +256,7 @@ def plot_image3(image, bt, cmap='Greys'):
def make_icing_image(h5f, probs, ice_lons, ice_lats, clvrx_str_time, satellite, domain, ice_lons_vld=None, ice_lats_vld=None, def make_icing_image(h5f, probs, ice_lons, ice_lats, clvrx_str_time, satellite, domain, ice_lons_vld=None, ice_lats_vld=None,
ice_intensity=None, imagename='icing', extent=[-110, -60, 10, 55]): ice_intensity=None, imagename='icing', extent=[-110, -60, 10, 55]):
mrkr_styl = MarkerStyle('o', fillstyle=None)
geos, xlen, xmin, xmax, ylen, ymin, ymax = get_cartopy_crs(satellite, domain) geos, xlen, xmin, xmax, ylen, ymin, ymax = get_cartopy_crs(satellite, domain)
bg_image = None bg_image = None
if h5f is not None: if h5f is not None:
...@@ -292,12 +294,12 @@ def make_icing_image(h5f, probs, ice_lons, ice_lats, clvrx_str_time, satellite, ...@@ -292,12 +294,12 @@ def make_icing_image(h5f, probs, ice_lons, ice_lats, clvrx_str_time, satellite,
lats = ice_lats_vld[ice_intensity == 1] lats = ice_lats_vld[ice_intensity == 1]
ax.scatter(lons, lats, s=500.0, marker='o', color='white', transform=ccrs.PlateCarree()) ax.scatter(lons, lats, s=500.0, marker='o', color='white', transform=ccrs.PlateCarree())
ax.scatter(lons, lats, s=200.0, marker='o', color='black', transform=ccrs.PlateCarree()) ax.scatter(lons, lats, s=200.0, marker='o', color='black', transform=ccrs.PlateCarree())
ax.scatter(lons, lats, s=1000.0, fillstyle=None, marker='o', color='black', transform=ccrs.PlateCarree()) ax.scatter(lons, lats, s=1000.0, marker=mrkr_styl, color='black', transform=ccrs.PlateCarree())
lons = ice_lons_vld[ice_intensity > 1] lons = ice_lons_vld[ice_intensity > 1]
lats = ice_lats_vld[ice_intensity > 1] lats = ice_lats_vld[ice_intensity > 1]
ax.scatter(lons, lats, s=500.0, marker='^', color='white', transform=ccrs.PlateCarree()) ax.scatter(lons, lats, s=500.0, marker='^', color='white', transform=ccrs.PlateCarree())
ax.scatter(lons, lats, s=200.0, marker='^', color='black', transform=ccrs.PlateCarree()) ax.scatter(lons, lats, s=200.0, marker='^', color='black', transform=ccrs.PlateCarree())
ax.scatter(lons, lats, s=1000.0, fillstyle=None, marker='o', color='black', transform=ccrs.PlateCarree()) ax.scatter(lons, lats, s=1000.0, marker=mrkr_styl, color='black', transform=ccrs.PlateCarree())
lons = ice_lons_vld[ice_intensity == -1] lons = ice_lons_vld[ice_intensity == -1]
lats = ice_lats_vld[ice_intensity == -1] lats = ice_lats_vld[ice_intensity == -1]
# ax.scatter(lons, lats, s=500.0, marker='o', color='white', transform=ccrs.PlateCarree()) # ax.scatter(lons, lats, s=500.0, marker='o', color='white', transform=ccrs.PlateCarree())
......
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