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
1268956107,I_kwDOAMm_X85LorvL,6693,"When a raster is modified, it cannot be read correctly.",27909361,closed,0,,,4,2022-06-13T06:15:44Z,2022-06-13T20:44:00Z,2022-06-13T20:44:00Z,NONE,,,,"### What happened?

As the following code shows, when I read this raster for the first time, xarray can read its content correctly. When this raster is modified, it can't be read correctly.

### What did you expect to happen?

_No response_

### Minimal Complete Verifiable Example

```Python
import rasterio
import xarray as xr
import numpy as np
import os

data = np.arange(10*10.).reshape(10, 10)
inp_rs = ""demo.tif""

def generate_raster(data):
    with rasterio.open(
        inp_rs,
        'w',
        driver='GTiff',
        height=data.shape[0],
        width=data.shape[1],
        count=1,
        dtype=data.dtype,
        crs='+proj=latlong',
    ) as dst:
        dst.write(data, 1)

generate_raster(data)

rs1 = rasterio.open(inp_rs)
print(""rasterio :"", rs1.shape)
rs2 = xr.open_dataset(inp_rs)
print(""xarray: "", rs2.band_data.shape)

# modify it
os.remove(inp_rs)
data_new = np.arange(9*9.).reshape(9, 9)
generate_raster(data_new)

rs1 = rasterio.open(inp_rs)
print(""rasterio :"", rs1.shape)
rs2 = xr.open_dataset(inp_rs)
print(""xarray: "", rs2.band_data.shape)
```


### MVCE confirmation

- [X] Minimal example — the example is as focused as reasonably possible to demonstrate the underlying issue in xarray.
- [X] Complete example — the example is self-contained, including all data and the text of any traceback.
- [X] 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.
- [X] New issue — a search of GitHub Issues suggests this is not a duplicate.

### Relevant log output

```Python
Read correctly the first time：
rasterio : (10, 10)
xarray:  (1, 10, 10)

Incorrectly read the second time：
rasterio : (9, 9)
xarray:  (1, 10, 10)
```


### Anything else we need to know?

_No response_

### Environment

<details>

INSTALLED VERSIONS
------------------
commit: None
python: 3.8.12 | packaged by conda-forge | (default, Oct 12 2021, 21:57:06) 
[GCC 9.4.0]
python-bits: 64
OS: Linux
OS-release: 5.3.0-7625-generic
machine: x86_64
processor: x86_64
byteorder: little
LC_ALL: None
LANG: en_US.UTF-8
LOCALE: ('en_US', 'UTF-8')
libhdf5: 1.12.0
libnetcdf: 4.7.4
xarray: 2022.3.0
pandas: 1.4.2
numpy: 1.22.3
scipy: 1.8.0
netCDF4: 1.5.8
pydap: None
h5netcdf: None
h5py: 3.6.0
Nio: None
zarr: None
cftime: 1.6.0
nc_time_axis: None
PseudoNetCDF: None
rasterio: 1.2.10
cfgrib: None
iris: None
bottleneck: None
dask: 2022.02.1
distributed: 2022.2.1
matplotlib: 3.3.2
cartopy: None
seaborn: None
numbagg: None
fsspec: 2022.3.0
cupy: None
pint: None
sparse: None
setuptools: 59.8.0
pip: 22.0.4
conda: None
pytest: None
IPython: 8.1.1
sphinx: None

</details>
","{""url"": ""https://api.github.com/repos/pydata/xarray/issues/6693/reactions"", ""total_count"": 0, ""+1"": 0, ""-1"": 0, ""laugh"": 0, ""hooray"": 0, ""confused"": 0, ""heart"": 0, ""rocket"": 0, ""eyes"": 0}",,completed,13221727,issue