From 9bb2a68ceb5179ec1aa8a07e30ef2ecd4dea41be Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Sun, 12 Dec 2021 14:50:08 -0600
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 12 ++++++++----
1 file changed, 8 insertions(+), 4 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 25b24d37..a0c32e75 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -1795,6 +1795,8 @@ def analyze(preds_file, test_file='/Users/tomrink/data/icing_ml/tiles_2021092400
print('Icing: ', np.histogram(nda[iint == 1], bins=5)[0])
print('---------------------------')
+ print('No Icing(Negative): mean cld_dz, cld_hgt')
+ print('Icing(Positive): ", "')
print('level 0: ')
print(np.nanmean(cld_dz[(nda == 0) & (iint == 0)]), np.nanmean(cld_hgt[(nda == 0) & (iint == 0)]))
print(np.nanmean(cld_dz[(nda == 0) & (iint == 1)]), np.nanmean(cld_hgt[(nda == 0) & (iint == 1)]))
@@ -1832,10 +1834,12 @@ def analyze(preds_file, test_file='/Users/tomrink/data/icing_ml/tiles_2021092400
false_ice = (labels == 0) & (preds == 1)
print('Total (Positive/Icing Prediction: ')
- print('True icing: ', np.histogram(nda[true_ice], bins=5)[0])
+ print('True icing: ', np.histogram(nda[true_ice], bins=5)[0])
print('-------------------------')
- print('False icing (False Positive/Alarm): ', np.histogram(nda[false_ice], bins=5)[0])
+ print('False icing (False Positive/Alarm): ', np.histogram(nda[false_ice], bins=5)[0])
print('By flight level:')
+ print('No Icing(Negative): mean cld_dz, cld_hgt')
+ print('Icing(Positive): ", "')
print('level 0: ')
print(np.nanmean(cld_dz[(nda == 0) & false_ice]), np.nanmean(cld_hgt[(nda == 0) & false_ice]))
print(np.nanmean(cld_dz[(nda == 0) & true_ice]), np.nanmean(cld_hgt[(nda == 0) & true_ice]))
@@ -1866,9 +1870,9 @@ def analyze(preds_file, test_file='/Users/tomrink/data/icing_ml/tiles_2021092400
false_no_ice = (labels == 1) & (preds == 0)
print('Total (Negative/No Icing Prediction: ')
- print('True no icing: ', np.histogram(nda[true_no_ice], bins=5)[0])
+ print('True no icing: ', np.histogram(nda[true_no_ice], bins=5)[0])
print('-------------------------')
- print('* False no icing (False Negative/Miss) *: ', np.histogram(nda[false_no_ice], bins=5)[0])
+ print('* False no icing (False Negative/Miss) *: ', np.histogram(nda[false_no_ice], bins=5)[0])
print('-------------------------')
print('level 0: ')
--
GitLab