Consider array band count when writing Image about ukis-pysat

[MWieland]

I totally agree with that. So from a user's perspective we would call something like this:

img = Image(array=myarray, crs=mycrs, transform=mytransform, dimorder="first")

This would be more meaningful, less confusing and more flexible (I like the link to easily go from GDAL to ukis-pysat).

[fwfichtner]

OT: Please be aware that you are quoting an old version of write_to_file().

I also don't quite understand what you're doing, img is not a dataset. This method is also not meant to write a new array to file, at least that's not how I intended it, for such purposes I would suggest using rasterio directly with something like this:

        with rasterio.open(
            destination,
            "w",
            driver="GTiff",
            width=ndvi_arr.shape[-1],
            height=ndvi_arr.shape[-2],
            count=1,
            dtype=np.float32,
            transform=self.transform,
        ) as dst:
            dst.write(ndvi_arr, indexes=1)

[MWieland]

Just edited the above comment. Should be img.dataset.

[fwfichtner]

The dataset is not the dataset you're looking for and this all becomes very confusing then. But I agree that it is very inviting to do something like that.
As an alternative we could consider something else:
Instead of initializing Image with a dataset (itself confusing, GDAL datasets will not work I guess) and array we could initalize it with crs, affine_transform and array. Then we could use an adaption of __update_dataset() to create the correct dataset internally. This would also be nice to go from GDAL to ukis-pysat in memory because you just have to convert the transform.