Fixed ValueError in case NaN is used as nodata value for float data type. Bumped version.

Signed-off-by: Daniel Scheffler <danschef@gfz-potsdam.de>

HISTORY.rst

@@ -3,10 +3,11 @@ History
 =======
 
 
-0.15.2 (coming soon)
---------------------
+0.15.2 (29.10.2021)
+-------------------
 
 * Added compatibility to the 'HDF:"/path/file.hdf":subdataset' syntax when opening HDF sub-datasets directly.
+* Fixed ValueError in case NaN is used as nodata value for float data type.
 
 
 0.15.1 (22.10.2021)

geoarray/baseclasses.py

@@ -156,8 +156,11 @@ class GeoArray(object):
             return
 
         if not dt_min <= _nodata <= dt_max:
-            raise ValueError("The given no-data value (%s) is out range for data type %s."
-                             % (self._initParams['nodata'], str(np.dtype(self.dtype))))
+            if np.issubdtype(self.dtype, np.floating) and np.isnan(_nodata):
+                pass
+            else:
+                raise ValueError("The given no-data value (%s) is out range for data type %s."
+                                 % (self._initParams['nodata'], str(np.dtype(self.dtype))))
 
     @property
     def arr(self):

geoarray/version.py

@@ -22,5 +22,5 @@
 # limitations under the License.
 
 
-__version__ = '0.15.1'
-__versionalias__ = '20211022.01'
+__version__ = '0.15.2'
+__versionalias__ = '20211029.01'

tests/test_geoarray.py

@@ -185,6 +185,9 @@ class Test_GeoArray(TestCase):
         # test float nodata value together with integer array data type
         GeoArray(np.array([1, 2]).astype(np.uint8), nodata=0.)
 
+        # test NaN nodata value together with float array data type
+        GeoArray(np.array([1, 2]).astype(np.float), nodata=np.nan)
+
     @parameterized.expand(_get_gA_inMem_notInMem())
     def test_calc_mask_nodata(self, _, gA):
         def test_with_flag_all(_gA):