docs: def execute(self, eopatch):
          """ Compute spatio-temporal features for input eopatch
  
          :param eopatch: Input eopatch
          :return: eopatch with computed spatio-temporal features
          """
          # pylint: disable=invalid-name
          amax_ndvi, amin_ndvi = eopatch.data_timeless[self.argmax_ndvi], eopatch.data_timeless[self.argmin_ndvi]
          amax_ndvi_slope, amin_ndvi_slope = eopatch.data_timeless[self.argmax_ndvi_slope], \
                                             eopatch.data_timeless[self.argmin_ndvi_slope]
          amax_red = eopatch.data_timeless[self.argmax_red]
  
          stf_idx = [amax_ndvi, amin_ndvi, amax_ndvi_slope, amin_ndvi_slope, amax_red]
  
          bands = eopatch.data[self.data_feature][..., self.indices]
  
          _, h, w, _ = bands.shape
          hh, ww = np.ogrid[:h, :w]
          stf = np.concatenate([bands[ii.squeeze(), hh, ww] for ii in stf_idx if ii is not None], axis=-1)
  
          eopatch.data_timeless[self.feats_feature] = stf
  
          return eopatch