code: def _open_xarray_dataset(self, val, chunks=CHUNK_SIZE):
          
          dask_arr = from_sds(val, chunks=chunks)
          attrs = val.attributes()
          return xr.DataArray(dask_arr, dims=('y', 'x'),
                              attrs=attrs)
docs: Read the band in blocks.