(active_win, baseline_win, epochs)
| 103 | |
| 104 | |
| 105 | def _gen_dics(active_win, baseline_win, epochs): |
| 106 | freqs = np.logspace(np.log10(12), np.log10(30), 9) |
| 107 | csd = csd_morlet(epochs, freqs, tmin=-1, tmax=1.5, decim=20) |
| 108 | csd_baseline = csd_morlet( |
| 109 | epochs, freqs, tmin=baseline_win[0], tmax=baseline_win[1], decim=20 |
| 110 | ) |
| 111 | csd_ers = csd_morlet( |
| 112 | epochs, freqs, tmin=active_win[0], tmax=active_win[1], decim=20 |
| 113 | ) |
| 114 | filters = make_dics( |
| 115 | epochs.info, |
| 116 | fwd, |
| 117 | csd.mean(), |
| 118 | pick_ori="max-power", |
| 119 | reduce_rank=True, |
| 120 | real_filter=True, |
| 121 | rank=rank, |
| 122 | ) |
| 123 | stc_base, freqs = apply_dics_csd(csd_baseline.mean(), filters) |
| 124 | stc_act, freqs = apply_dics_csd(csd_ers.mean(), filters) |
| 125 | stc_act /= stc_base |
| 126 | return stc_act |
| 127 | |
| 128 | |
| 129 | # generate lcmv source estimate |
no test coverage detected