code: def _compute_pow_plv(
      data,
      K,
      sel,
      Ws,
      source_ori,
      use_fft,
      Vh,
      with_power,
      with_plv,
      pick_ori,
      decim,
      noise_norm=None,
      verbose=None,
  ):
      """Aux function for induced power and PLV."""
      shape, is_free_ori = _prepare_tfr(
          data, decim, pick_ori, Ws, K, source_ori
      )
      power = np.zeros(shape, dtype=np.float64)  # power or raw TFR
      # phase lock
      plv = np.zeros(shape, dtype=np.complex128) if with_plv else None
  
      for epoch in data:
          epoch = epoch[sel]  # keep only selected channels
  
          if Vh is not None:
              epoch = np.dot(Vh, epoch)  # reducing data rank
  
          power_e, plv_e = _single_epoch_tfr(
              data=epoch,
              is_free_ori=is_free_ori,
              K=K,
              Ws=Ws,
              use_fft=use_fft,
              decim=decim,
              shape=shape,
              with_plv=with_plv,
              with_power=with_power,
          )
  
          power += power_e
          if with_plv:
              plv += plv_e
  
      if noise_norm is not None:
          power *= noise_norm[:, :, np.newaxis] ** 2
  
      return power, plv
  
imports: import numpy
  from .._fiff.constants import FIFF
  from .._fiff.pick import pick_info
  from ..baseline import _log_rescale, rescale
  from ..epochs import Epochs
  from ..event import make_fixed_length_events
  from ..evoked import EvokedArray
  from ..fixes import _safe_svd
  from ..label import BiHemiLabel, Label
  from ..parallel import parallel_func
  from ..source_estimate import _make_stc
  from ..time_frequency.multitaper import _compute_mt_params, _mt_spectra, _psd_from_mt, _psd_from_mt_adaptive
  from ..time_frequency.tfr import cwt, morlet
  from ..utils import ProgressBar, _check_option, _pl, _validate_type, logger, verbose
  from .inverse import INVERSE_METHODS, _assemble_kernel, _check_or_prepare, _check_ori, _pick_channels_inverse_operator, _subject_from_inverse, combine_xyz
label_desc: Use image node_id 5 for calling a global function with example usage: _compute_pow_plv(data, K, sel, Ws, source_ori, use_fft, Vh, with_power, with_plv, pick_ori, decim, noise_norm, verbose) and returns: power,  plv