# Extract power bands from epochs

**URL:** <https://mne.discourse.group/t/extract-power-bands-from-epochs/5743>\
**Category:** Support & Discussions\
**Tags:** eeg, epochs\
**Created:** [October 16, 2022, 8:28pm UTC](https://mne.discourse.group/t/extract-power-bands-from-epochs/5743 "2022-10-16T20:28:15Z")\
**Posts on this page:** 3\
**Page:** 1

<div class="post-metadata">

**Author:** ![MNE](https://avatars.discourse-cdn.com/v4/letter/m/8edcca/32.png) [@MNE](https://mne.discourse.group/u/MNE)\
**Post date:** [October 16, 2022, 8:28pm UTC](https://mne.discourse.group/t/extract-power-bands-from-epochs/5743/1 "2022-10-16T20:28:15Z")

</div>

**❓ If you have a question or issue with MNE-Python, please include the following info:**

- MNE version: e.g. 0.24.0
- operating system: e.g. macOS 12 / Windows 10 / Ubuntu 18.04

Is there any function to extract frequency bands from epochs not the raw data?

---

<div class="post-metadata">

**Author:** ![mscheltienne](https://yyz2.discourse-cdn.com/free1/user_avatar/mne.discourse.group/mscheltienne/32/827_2.png) [@mscheltienne](https://mne.discourse.group/u/mscheltienne)\
**Post date:** [October 17, 2022, 11:02am UTC](https://mne.discourse.group/t/extract-power-bands-from-epochs/5743/2 "2022-10-17T11:02:21Z")

</div>

Hello,

The PSD functions, `multitaper` or `welch` can both be run on `Raw`,`Epochs` or `Evoked` object.  
See for instance the argument for `multitaper` in the version 0.24

[https://mne.tools/0.24/generated/mne.time\_frequency.psd\_multitaper.html?highlight=psd#mne.time\_frequency.psd\_multitaper](https://mne.tools/0.24/generated/mne.time_frequency.psd_multitaper.html?highlight=psd#mne.time_frequency.psd_multitaper)

Mathieu

---

<div class="post-metadata">

**Author:** ![drammock](https://yyz2.discourse-cdn.com/free1/user_avatar/mne.discourse.group/drammock/32/4_2.png) [@drammock](https://mne.discourse.group/u/drammock)\
**Post date:** [October 17, 2022, 4:24pm UTC](https://mne.discourse.group/t/extract-power-bands-from-epochs/5743/3 "2022-10-17T16:24:59Z")

</div>

Are you actually using version 0.24, or is that just the placeholder text? If you’re on the current stable version (1.2) you can use `epochs.compute_psd()` to get a `Spectrum` object, and then `psd, freqs = spectrum.get_data(return_freqs=True)` . From this you can figure out which bin(s) contain the frequencies of interest, and then average those bins together to get average power in the band.