def METHOD_NAME(self):
    import dask.array as da
    self._get_schema()
    arrs = [NumpyAccess(f, self.shape, self.dtype) for f in self.files]
    if len(arrs) > 1:
        chunks = (self._chunks or -1,) + (-1,) * (len(self.shape) - 1)
    else:
        chunks = (self._chunks or -1,) + (-1,) * (len(self.shape) - 2)
    self._arrs = [da.from_array(arr, chunks) for arr in arrs]
    self._arr = da.stack(self._arrs)
    return self._arr