fixing MultiRaster reprojection

pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "rasters"
-version = "1.4.2"
+version = "1.4.3"
 description = "raster processing toolkit"
 readme = "README.md"
 authors = [

rasters/multi_raster.py

@@ -73,6 +73,27 @@ def __init__(
         self.cmap = cmap
         self._source_metadata = {}
 
+    def contain(self, array=None, geometry=None, nodata=None) -> Raster:
+        if array is None:
+            array = self.array
+
+        if geometry is None:
+            geometry = self.geometry
+
+        if np.size(array) == 1:
+            return array
+
+        if nodata is None:
+            nodata = self.nodata
+
+        return MultiRaster(
+            array,
+            nodata=nodata,
+            metadata=self.metadata,
+            cmap=self.cmap,
+            geometry=geometry
+        )
+
     @classmethod
     def stack(cls, rasters: List[Raster], *args, **kwargs) -> MultiRaster:
         geometry = rasters[0].geometry

rasters/raster.py

@@ -532,25 +532,33 @@ def contain(self, array=None, geometry=None, nodata=None) -> Raster:
         if nodata is None:
             nodata = self.nodata
 
-        # try:
-        if len(array.shape) == 2:
-            return Raster(
-                array,
-                nodata=nodata,
-                metadata=self.metadata,
-                cmap=self.cmap,
-                geometry=geometry
-            )
-        elif len(array.shape) == 3:
-            return MultiRaster(
-                array,
-                nodata=nodata,
-                metadata=self.metadata,
-                cmap=self.cmap,
-                geometry=geometry
-            )
-        else:
-            raise ValueError(f"invalid raster array with shape {array.shape}")
+        return Raster(
+            array,
+            nodata=nodata,
+            metadata=self.metadata,
+            cmap=self.cmap,
+            geometry=geometry
+        )
+
+        # # try:
+        # if len(array.shape) == 2:
+        #     return Raster(
+        #         array,
+        #         nodata=nodata,
+        #         metadata=self.metadata,
+        #         cmap=self.cmap,
+        #         geometry=geometry
+        #     )
+        # elif len(array.shape) == 3:
+        #     return MultiRaster(
+        #         array,
+        #         nodata=nodata,
+        #         metadata=self.metadata,
+        #         cmap=self.cmap,
+        #         geometry=geometry
+        #     )
+        # else:
+        #     raise ValueError(f"invalid raster array with shape {array.shape}")
         # except Exception as e:
         #     return array
 

rasters/version.txt

@@ -1 +1 @@
-1.4.2
+1.4.3