support multiple bands in rasterio read

[gidden]

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

[EdgardBontempo]

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