code: def read(self, indexes=None, **kwargs):
          
          return read_raster_window(
              self.raster_file.path,
              self.tile,
              indexes=self._get_band_indexes(indexes),
              resampling=self.resampling,
              gdal_opts=self.gdal_opts
          )
docs: Read reprojected & resampled input data.
  
          Returns
          -------
          data : array