id,node_id,number,title,user,state,locked,assignee,milestone,comments,created_at,updated_at,closed_at,author_association,active_lock_reason,draft,pull_request,body,reactions,performed_via_github_app,state_reason,repo,type 1441345452,I_kwDOAMm_X85V6S-s,7273,concat breaks raster stack,630436,closed,1,,,5,2022-11-09T04:15:07Z,2023-08-01T16:29:10Z,2023-08-01T16:29:09Z,NONE,,,,"### What happened? Loading multiple tiffs into a raster stack via `concat`, `open_mfdataset`, `open_dataset`, `combine_by_coords` corrupts the data, causing them to not display correctly, or write out to disk correctly. The rasters show as banded outputs missing data. Plotting individual pre-stacked looks as expected. If a slice of the stacked DA is written to a file, that file is also corrupted. This is confirmed by loading the file in QGIS. ### What did you expect to happen? The stack to not mangle the data ### Minimal Complete Verifiable Example ```Python import xarray as xr import rioxarray as rio import glob files = sorted(glob.glob('nomax_tiff/*tiff')) print(files) # open the files d0=rio.open_rasterio(files[0]).to_dataset(name='vel').assign_coords(time=0).expand_dims(dim=""time"") d1=rio.open_rasterio(files[1]).to_dataset(name='vel').assign_coords(time=1).expand_dims(dim=""time"") # should look fine d0.vel.plot() d1.vel.plot() ds = xr.concat([d0,d1],dim='time') #mangled ds.vel.plot(col='time') #this is also mangled ds.drop('band').squeeze().isel(time=0).rio.to_raster('mangled.tiff') #without rioxarray/rasterio it is differently mangled xr.open_mfdataset(files,concat_dim='time',combine='nested').isel(time=0).band_data.plot() ``` ### MVCE confirmation - [ ] Minimal example — the example is as focused as reasonably possible to demonstrate the underlying issue in xarray. - [ ] Complete example — the example is self-contained, including all data and the text of any traceback. - [ ] Verifiable example — the example copy & pastes into an IPython prompt or [Binder notebook](https://mybinder.org/v2/gh/pydata/xarray/main?urlpath=lab/tree/doc/examples/blank_template.ipynb), returning the result. - [ ] New issue — a search of GitHub Issues suggests this is not a duplicate. ### Relevant log output _No response_ ### Anything else we need to know? [nomax_tiff.zip](https://github.com/pydata/xarray/files/9967216/nomax_tiff.zip) ### Environment On a linux cluster:
xr.show_versions() INSTALLED VERSIONS ------------------ commit: None python: 3.8.10 (default, Jun 16 2021, 14:20:20) [GCC 9.3.0] python-bits: 64 OS: Linux OS-release: 3.10.0-1160.36.2.el7.x86_64 machine: x86_64 processor: byteorder: little LC_ALL: None LANG: en_CA.UTF-8 LOCALE: ('en_CA', 'UTF-8') libhdf5: 1.10.6 libnetcdf: 4.7.4 xarray: 2022.6.0+computecanada pandas: 1.4.0 numpy: 1.22.2 scipy: 1.8.0 netCDF4: 1.5.8 pydap: None h5netcdf: None h5py: 3.1.0 Nio: None zarr: None cftime: 1.6.0 nc_time_axis: None PseudoNetCDF: None rasterio: 1.3.0 cfgrib: 0.9.10.1 iris: None bottleneck: None dask: 2022.7.0 distributed: None matplotlib: 3.5.1 cartopy: 0.20.3 seaborn: 0.11.2 numbagg: None fsspec: 2022.5.0 cupy: None pint: None sparse: None flox: None numpy_groupies: None setuptools: 46.1.3 pip: 20.0.2 conda: None pytest: None IPython: 7.31.1 sphinx: None
This also reproduces on a macos install. But xr show_versions segfaults. No idea, seems weird.
Python 3.10.4 (main, Jun 14 2022, 14:00:56) [Clang 13.0.0 (clang-1300.0.27.3)] on darwin rioxarray 0.12.4 xarray 2022.11.0 rasterio 1.3.3
","{""url"": ""https://api.github.com/repos/pydata/xarray/issues/7273/reactions"", ""total_count"": 0, ""+1"": 0, ""-1"": 0, ""laugh"": 0, ""hooray"": 0, ""confused"": 0, ""heart"": 0, ""rocket"": 0, ""eyes"": 0}",,completed,13221727,issue