docs: def execute(self, eopatch):
          """ Compute argmax/argmin of specified `data_feature` and `data_index`
  
          :param eopatch: Input eopatch
          :return: eopatch with added argmax/argmin features
          """
          if self.mask_data:
              valid_data_mask = eopatch.mask['VALID_DATA']
          else:
              valid_data_mask = eopatch.mask['IS_DATA']
  
          if self.data_index is None:
              data = eopatch.data[self.data_feature]
          else:
              data = eopatch.data[self.data_feature][..., self.data_index]
  
          madata = np.ma.array(data,
                               dtype=np.float32,
                               mask=~valid_data_mask.astype(np.bool))
  
          argmax_data = np.ma.MaskedArray.argmax(madata, axis=0)
          argmin_data = np.ma.MaskedArray.argmin(madata, axis=0)
  
          if argmax_data.ndim == 2:
              argmax_data = argmax_data.reshape(argmax_data.shape + (1,))
  
          if argmin_data.ndim == 2:
              argmin_data = argmin_data.reshape(argmin_data.shape + (1,))
  
          eopatch.data_timeless[self.amax_feature] = argmax_data
          eopatch.data_timeless[self.amin_feature] = argmin_data
  
          return eopatch