Bugfixes.

geoarray/baseclasses.py

@@ -4,6 +4,7 @@ import os
 import warnings
 from importlib import util
 from collections import OrderedDict
+from copy import deepcopy
 
 import numpy as np
 from matplotlib import pyplot as plt
@@ -1476,7 +1477,11 @@ class GeoArray(object):
         # apply zslice to bandnames and metadata
         if zslicing:
             bNs_out = list(np.array(list(self.bandnames))[zslice])
-            _meta_out = None if self._metadata is None else self.meta[list(np.array(range(self.bands))[zslice])]
+
+            _meta_out = deepcopy(self.meta)
+            if self._metadata:
+                for k, v in _meta_out.band_meta.items():
+                    _meta_out.band_meta[k] = list(np.array(v)[zslice])
         else:
             bNs_out = list(self.bandnames)
             _meta_out = self.meta

tests/test_geoarray.py

@@ -24,7 +24,6 @@ Program started in Mai 2017.
 from __future__ import print_function
 from collections import OrderedDict
 import dill
-import geopandas
 import numpy as np
 import os
 from os import path
@@ -43,6 +42,7 @@ from py_tools_ds.geo.vector import geometry
 from py_tools_ds.compatibility.python.os import makedirs
 matplotlib.use('Template')  # disables matplotlib figure popups
 from geoarray import GeoArray, masks, __file__  # noqa E402 module level import not at top of file
+from geoarray.metadata import GDAL_Metadata  # noqa E402 module level import not at top of file
 
 __author__ = 'Jessica Palka'
 
@@ -453,7 +453,7 @@ class Test_GeoarrayFunctions(unittest.TestCase):
         Test, if the output of the metadata-function is an instance of GeoDataFrame.
         """
 
-        self.assertIsInstance(self.testtiff.metadata, geopandas.GeoDataFrame)
+        self.assertIsInstance(self.testtiff.metadata, GDAL_Metadata)
 
     def test_tiles(self):
         test_gAs = [self.testtiff,  # 3D