From 789f30ed3dd49927206ff9811e1e824e34b4ba7b Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Fri, 13 May 2022 10:59:41 -0500
Subject: [PATCH] snapshot...
---
modules/util/viirs.py | 33 +++++++++++++++++++++++++++++++++
1 file changed, 33 insertions(+)
diff --git a/modules/util/viirs.py b/modules/util/viirs.py
index 4e1e67f3..b4d66cbc 100644
--- a/modules/util/viirs.py
+++ b/modules/util/viirs.py
@@ -95,6 +95,39 @@ def run(mod_res_filename, img_res_filename, mod_tiles, img_tiles):
img_h5f.close()
+def run_mean_std(directory):
+
+ data_dct = {name: [] for name in mod_res_params}
+ mean_dct = {name: 0 for name in mod_res_params}
+ std_dct = {name: 0 for name in mod_res_params}
+
+ cnt = 0
+ for p in os.scandir(directory):
+ if not p.is_dir():
+ continue
+ cnt += 1
+ mod_files = glob.glob(directory+p.name+'/'+'*.uwssec.nc')
+
+ for idx, mfile in enumerate(mod_files):
+ if idx % 8 == 0:
+ h5f = h5py.File(mfile, 'r')
+ for param in mod_res_params:
+ gvals = get_grid_values_all(h5f, param, range_name=None, stride=10)
+ data_dct[param].append(gvals.flatten())
+ print(mfile)
+ h5f.close()
+
+ for param in mod_res_params:
+ data = data_dct[param]
+ data = np.concatenate(data)
+
+ mean_dct[param] = np.nanmean(data)
+ std_dct[param] = np.nanstd(data)
+
+ print(mean_dct)
+ print(std_dct)
+
+
def run_test(directory):
# files = glob.glob(directory + 'clavrx_snpp_viirs*.h5')
files = Path(directory).rglob('clavrx_snpp_viirs*.h5')
--
GitLab