docs: def _get_indices(self, data):
          """ Compute indices along temporal dimension corresponding to the sought percentile
  
              :param data: Input 3D array holding the reference band
              :type data: numpy array
              :return: 2D array holding the temporal index corresponding to percentile
          """
          indices = self._index_by_percentile(data, self.percentile)
          return indices