code: def _read_as_tiledir(
          self,
          out_tile=None,
          td_crs=None,
          tiles_paths=None,
          profile=None,
          validity_check=False,
          indexes=None,
          resampling=None,
          dst_nodata=None,
          gdal_opts=None,
          **kwargs
      ):
          
          return _read_as_tiledir(
              data_type=self.METADATA["data_type"],
              out_tile=out_tile,
              td_crs=td_crs,
              tiles_paths=tiles_paths,
              profile=profile,
              validity_check=validity_check,
              indexes=indexes,
              resampling=resampling,
              dst_nodata=dst_nodata,
              gdal_opts=gdal_opts,
              **{k: v for k, v in kwargs.items() if k != "data_type"}
          )
docs: Read reprojected & resampled input data.
  
          Parameters
          ----------
          validity_check : bool
              vector file: also run checks if reprojected geometry is valid,
              otherwise throw RuntimeError (default: True)
  
          indexes : list or int
              raster file: a list of band numbers; None will read all.
          dst_nodata : int or float, optional
              raster file: if not set, the nodata value from the source dataset
              will be used
          gdal_opts : dict
              raster file: GDAL options passed on to rasterio.Env()
  
          Returns
          -------
          data : list for vector files or numpy array for raster files