robintw
support multiple bands in rasterio read
It looks like this only supports a single band for now [1]. This could pretty easily be extended to 3 dims (most common for geoscience cases) by assuming a 3D array (with appropriate 2D support). I'd be happy to look into this in the near term as this is something that would help me.
@robintw, do you have any interest in connecting this with pydata/xarray#1070 ?
[1] https://github.com/robintw/XArrayAndRasterio/blob/master/rasterio_to_xarray.py#L26
You guys been working on this or has some better way to do it come up since then?
I tried using it - as I appreciate xarray's functionality but need to get geotiffs from my data - but had problems I feel may have been caused by own (almost complete) inexperience with python.
Details HERE, at stackoverflow.
Anyway, I really appreciate the code. Kudos, Robin Wilson
