<strong>This message contained a link to a web site whose IP address has recently been associated with malicious activity. Due to that activity, the site is not accessible from within the Partners Healthcare network. The link has been modified in the message below. If you feel this is an error, please contact the IS Service Desk and open a ticket with the network security - phs assignment group.<br><span style="font-size: 13px;"><br> IS Service Desk<br></span><span style="font-size: 13px;"><span style="font-size: 13px;">BWH 617-732-5927 BWFH 617-983-7454 BWH-RICS 617-525-0848<br>DFCI 617-632-3399 LCC 857-307-4150 MCL 781-416-8940<br>MGH 617-726-5085 NHP 617-643-2020 NSMC 978-354-2014<br>NWH 617-243-6001 PCPO 781-433-3757 PHH 617-726-0790</span><br>PHS 857-282-4357 SRN 617-952-5555 CDHC 413-582-5005</span></strong><br>*********************************<br><br><span style="font-size: 13px;"> </span><br>
<p><span style="padding: 3px 10px; border-radius: 5px; color: #ffffff; font-weight: bold; display: inline-block; background-color: #ff0000;"> External Email - Use Caution </span></p><p></p><div dir="ltr"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div lang="EN-GB"><div class="gmail-m_3637811604928529703WordSection1"><p class="MsoNormal">I’m currently running source recon MEG data using individual MRIs. I found that when running make_inverse_operator(<a href="http://raw.info" target="_blank">raw.info</a>, fwd, noise_cov) I was getting an error saying ‘no channels match the selection’. I checked in <a href="http://raw.info" target="_blank">raw.info</a> and no channels
were marked as bad, so I tried again using <a href="BLOCKEDepoch[.]infoBLOCKED" target="_blank"><font color="red"><b>MailScanner has detected a possible fraud attempt from "blockedepoch[.]infoblocked" claiming to be</b></font> epoch.info</a> (also no bad channels) instead of <a href="http://raw.info" target="_blank">raw.info</a> and got the same error.</p></div></div></blockquote><div><br></div><div>raw.ch_names, fwd['ch_names'], and/or noise_cov['names'] must not match. This suggests that something is going wrong in your pipeline as these should all be consistent.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div lang="EN-GB"><div class="gmail-m_3637811604928529703WordSection1"><p class="MsoNormal"><u></u><u></u></p>
<p class="MsoNormal">I have diagnosed this further up the pipeline – I was making the forward solution with:<u></u><u></u></p>
<p class="MsoNormal">make_forward_solution(<a href="http://raw.info" target="_blank">raw.info</a>, trans, src, bem, meg=True, eeg=False)
<u></u><u></u></p>
<p class="MsoNormal">Thinking that this would pick only the MEG channels.</p></div></div></blockquote><div><br></div><div>Indeed it should. You can check with fwd['sol']['row_names'].</div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div lang="EN-GB"><div class="gmail-m_3637811604928529703WordSection1"><p class="MsoNormal">However, my noise covariance is calculated from epochs, not the raw data, and I think that creating a forward solution using <a href="http://raw.info" target="_blank">raw.info</a> and the covariance using <a href="BLOCKEDepoch[.]infoBLOCKED" target="_blank"><font color="red"><b>MailScanner has detected a possible fraud attempt from "blockedepoch[.]infoblocked" claiming to be</b></font> epoch.info</a> resulted in a
clash of channels, even though both of them have 270 MEG channels.</p></div></div></blockquote><div><br></div><div>It shouldn't. Creating epochs should not change the channel names. I would try making your script as short as possible to reproduce these steps and see if they fail. Something like:\</div><div></div><div class="markdown-here-wrapper" style=""><pre style="font-size:0.85em;font-family:Consolas,Inconsolata,Courier,monospace;font-size:1em;line-height:1.2em;margin:1.2em 0px"><code style="font-size:0.85em;font-family:Consolas,Inconsolata,Courier,monospace;margin:0px 0.15em;padding:0px 0.3em;white-space:pre-wrap;border:1px solid rgb(234,234,234);background-color:rgb(248,248,248);border-radius:3px;display:inline;white-space:pre;overflow:auto;border-radius:3px;border:1px solid rgb(204,204,204);padding:0.5em 0.7em;display:block!important">raw = mne.io.read_raw_something(...)
events = mne.make_fixed_length_events(raw)
epochs = mne.Epochs(raw, events)
assert raw.ch_names == epochs.ch_names
cov = mne.compute_covariance(epochs)
assert cov['names'] == epochs.ch_names
fwd = make_forward_solution(<a href="http://raw.info">raw.info</a>, ...)
inv = make_inverse_operator(<a href="http://epochs.info">epochs.info</a>, fwd, noise_cov)
</code></pre><div title="MDH:PGRpdj5gYGA8L2Rpdj48ZGl2PnJhdyA9IG1uZS5pby5yZWFkX3Jhd19zb21ldGhpbmcoLi4uKTwv
ZGl2PjxkaXY+ZXZlbnRzID0gbW5lLm1ha2VfZml4ZWRfbGVuZ3RoX2V2ZW50cyguLi4pPC9kaXY+
PGRpdj5lcG9jaHMgPSBtbmUuRXBvY2hzKHJhdywgZXZlbnRzKTwvZGl2PjxkaXY+YXNzZXJ0IHJh
dy5jaF9uYW1lcyA9PSBlcG9jaHMuY2hfbmFtZXM8L2Rpdj48ZGl2PjwvZGl2PjxkaXY+Y292ID0g
bW5lLmNvbXB1dGVfY292YXJpYW5jZShlcG9jaHMpPC9kaXY+PGRpdj5hc3NlcnQgY292WyduYW1l
cyddID09IGVwb2Nocy5jaF9uYW1lczwvZGl2PjxkaXY+ZndkID0gbWFrZV9mb3J3YXJkX3NvbHV0
aW9uKHJhdy5pbmZvLCAuLi4pPC9kaXY+PGRpdj5pbnYgPSBtYWtlX2ludmVyc2Vfb3BlcmF0b3Io
ZXBvY2hzLmluZm8sIGZ3ZCwgbm9pc2VfY292KSZuYnNwOzxicj48L2Rpdj48ZGl2PmBgYDwvZGl2
PjxkaXY+PC9kaXY+" style="height:0;width:0;max-height:0;max-width:0;overflow:hidden;font-size:0em;padding:0;margin:0"></div></div><div>If the filled-out version of this fails on your data, there is probably some bug in MNE with these functions. If it works, then you can gradually add lines to it (or remove lines from your other script(s)) to figure out what step breaks things. It could be some bug in MNE or some bug in your script.<br></div><div><br></div><div>Eric</div><div><br></div></div></div>