Only compare the first two dimensions.

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

geoarray/baseclasses.py

@@ -1609,12 +1609,12 @@ class GeoArray(object):
         self.arr, self.gt, self.projection = self.get_mapPos(mapBounds=poly.bounds)
         self.mask_nodata.arr, self.mask_nodata.gt, self.mask_nodata.projection = \
             self.mask_nodata.get_mapPos(mapBounds=poly.bounds, mapBounds_prj=self.prj)
-        assert self.shape == self.mask_nodata.shape
+        assert self.shape[:2] == self.mask_nodata.shape
 
         if self._mask_baddata is not None:
             self.mask_baddata.arr, self.mask_baddata.gt, self.mask_baddata.projection = \
                 self.mask_baddata.get_mapPos(mapBounds=poly.bounds)
-            assert self.shape == self.mask_baddata.shape
+            assert self.shape[:2] == self.mask_baddata.shape
 
         # update footprint polygon
         if self._footprint_poly: