From c1ab71748e0c20b956e28bb9ea955ad092eafc04 Mon Sep 17 00:00:00 2001
From: Ben Jolly <jollyb@landcareresearch.co.nz>
Date: Wed, 22 Jan 2025 14:16:10 +1300
Subject: [PATCH] fix for issue where np.nan was not ignored for stats

---
 rios/calcstats.py | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/rios/calcstats.py b/rios/calcstats.py
index 16da717..6809667 100644
--- a/rios/calcstats.py
+++ b/rios/calcstats.py
@@ -639,6 +639,8 @@ def doStatsAccum(self, arr):
         """
         if self.nullval is None:
             values = arr.flatten()
+        elif numpy.isnan(self.nullval):
+            values = arr[~numpy.isnan(arr)]
         else:
             values = arr[arr != self.nullval]
         if len(values) > 0: