<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Dear Hanneke,<div><br></div><div>I assume that t(1:3,4) is the center of the figure of 8 coil (the crossing point). Adding or subtracting 0.008 along the correct dimension seems to be necessary to arrive at the center of each of the two loops that make up the planar gradiometer.</div><div><br></div><div>Joachim</div><div><br></div><div><br><div><div>On 10 Sep 2009, at 10:37, Hanneke van Dijk wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">Dear all,<br><br>At the moment I am using the MNE toolbox and Fieldtrip to localize an alpha source (1.5 seconds spontaneous data recorded with a Neuromag 122 system). <br><br>With data from the new 306 channel system it seems to work out well. However with data from the old system the source ends up in the cerebellum (attached a picture of Neural Activity index, cfg.lambda = '10%'). Of course it should be in the parieto-occipital cortex or the sensorimotor cortex. I have the feeling there is problem with the position of the gradiometers relative to the headmodel or the other way around. The head model seems to be too far away from the sensors (see attached headmodel1), if my hunch is correct I don't know how to fix this or where exactly I can find the problem. <br> <br>How I made the headmodel you can find on: <br><a href="http://fieldtrip.fcdonders.nl/example/read_neuromag_mri_and_create_single-subject_grids_in_individual_head_space_that_are_all_aligned_in_mni_space">http://fieldtrip.fcdonders.nl/example/read_neuromag_mri_and_create_single-subject_grids_in_individual_head_space_that_are_all_aligned_in_mni_space</a><br> <br>I'm constructing the grad files (to do sourceanalysis with) seperately using read_header and mne2grad (with orig.chs(n).coil_type=2, because the coil_type is 2 for the old system).<br><br>Another related question: in mne2grad on line 104 and 108 it says <br> <br> 104 >> grad.pnt(kCoil,:)=100*(t(1:3,4)-0.008*t(1:3,1)); % multiply with 100 to get cm<br> .....<br> 108 >> grad.pnt(kCoil,:)=100*(t(1:3,4)+0.008*t(1:3,1));<br><br>Where does this 0.008 come from? I have changed this number to try it out but did't see any clear differences in the gradiometer positions.<br> <br>Thanks in advance for your help,<br>Wishing you all the best,<br><br>Hanneke<br><br><br><p>----------------------------------</p><p>The aim of this list is to facilitate the discussion between users of the FieldTrip toolbox, to share experiences and to discuss new ideas for MEG and EEG analysis.</p><p> <a href="http://listserv.surfnet.nl/archives/fieldtrip.html">http://listserv.surfnet.nl/archives/fieldtrip.html</a></p><p> <a href="http://www.ru.nl/fcdonders/fieldtrip/">http://www.ru.nl/fcdonders/fieldtrip/</a></p> <span><headmodel1.jpeg></span><span><alphasource.jpeg></span></blockquote></div><br></div></body></html><p>----------------------------------</p>
<p>The aim of this list is to facilitate the discussion between users of the FieldTrip toolbox, to share experiences and to discuss new ideas for MEG and EEG analysis.</p>
<p> http://listserv.surfnet.nl/archives/fieldtrip.html</p>
<p> http://www.ru.nl/fcdonders/fieldtrip/</p>