From 928d4448f307d50a278ccd5f090a7e926e06809a Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Sun, 30 May 2021 09:57:09 -0500
Subject: [PATCH] remove old code, add some new code
---
modules/icing/pirep_goes.py | 33 ++++++++++++++++++++++++++-------
1 file changed, 26 insertions(+), 7 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 7d01553f..93f222d5 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1326,12 +1326,31 @@ def normalize(data, param, mean_std_dict):
return data
-def test(filename, skip=1):
- h5f = h5py.File(filename, 'r')
- time = h5f['time']
- intsty = h5f['icing_intensity']
+lon_space = np.linspace(-180, 180, 361)
+lat_space = np.linspace(-90, 90, 181)
- trn_idxs, tst_idxs = split_data(time.shape[0], skip=skip)
- print(np.histogram(intsty[trn_idxs], bins=7))
- print(np.histogram(intsty[tst_idxs], bins=7))
+def spatial_filter(icing_dict):
+ keys = icing_dict.keys()
+ grd_x_hi = lon_space.shape[0] - 1
+ grd_y_hi = lat_space.shape[0] - 1
+
+ for key in keys:
+
+ grd_bins = np.full((lat_space.shape[0], lon_space.shape[0]), 0)
+
+ tup = icing_dict.get(key)
+ lat = tup[0]
+ lon = tup[1]
+
+ lon_idx = np.searchsorted(lon_space, lon)
+ lat_idx = np.searchsorted(lat_space, lat)
+
+ if lon_idx < 0 or lon_idx > grd_x_hi:
+ continue
+ if lat_idx < 0 or lat_idx > grd_y_hi:
+ continue
+
+ grd_bins[lat_idx, lon_idx] += 1
+
+ return grd_bins
--
GitLab