Thanks that seemed to help the scaling issue. However, my beamforming is still not producing results. There is a huge bias towards the center of the head, even though I accounted for the neural activity index. It seems like the beamforming is just projecting everything inwards. I'm not sure where to go from here. Do you have any suggestions on different configuration settings I could try for the beamforming? <br>

<br>-Allison<br>
<br><div class="gmail_quote">On Fri, Feb 11, 2011 at 4:02 PM, Luisa Frei <span dir="ltr"><<a href="mailto:l.frei@psy.gla.ac.uk" target="_blank">l.frei@psy.gla.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">


Hi Alison,<br>
how do you change the units? I convert them using:<br>
<br>
data.grad=ft_convert_units(data.grad,'cm')<br>
<br>
Maybe this helps.<br>
<br>
Luisa<div><div></div><div><br>
<br>
On 11 Feb 2011, at 21:19, Allison Connolly wrote:<br>
<br>
</div></div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div><div></div><div>
Hello,<br>
<br>
I am trying to perform beamforming on some MEG data (148 channels). I have the dicom files for the MRI, and I have converted them to .img and .hdr files using Curry.<br>
I am importing them into fieldtrip via<br>
mri=ft_read_mri('MRI.img').<br>
<br>
Then I manually declaring the fiducials with<br>
cfg.method                  = 'interactive';<br>
mrialign                    = ft_volumerealign(cfg, mri);<br>
<br>
Then segmenting with<br>
cfg = [];<br>
cfg.write        = 'no';<br>
[segment] = ft_volumesegment(cfg, mri);<br>
<br>
Then creating the head model with<br>
cfg = [];<br>
cfg.smooth                  = 'no';<br>
vol = ft_prepare_singleshell(cfg, segment);<br>
<br>
Then creating the grid with<br>
cfg                 = [];<br>
cfg.grad            = freq.grad;<br>
cfg.vol             = vol;<br>
cfg.reducerank      = 'no';<br>
cfg.channel                 = {'MEG','-TRIGGER','-RESPONSE',...<br>
     '-UACurrent','-MxA','-MyA','-MzA','-GxxA','-GyxA','-GyyA',...<br>
     '-GzyA','-GzxA','-MxaA','-MyaA','-MzaA','-E1','-E2','-E3','-E4'};<br>
cfg.grid.xgrid              = [-10:0.5:15]; % adjust according to vol<br>
cfg.grid.ygrid              = [-15:0.5:15]; % adjust according to vol<br>
cfg.grid.zgrid              = [-10:0.5:15]; % adjust according to vol<br>
[grid] = ft_prepare_leadfield(cfg);<br>
<br>
I am importing my MEG data with<br>
cfg = [];  % empty configuration<br>
cfg.dataset                 = 'C:/MEG Study/Subjectd001/d001/Billy_180s/06@-09@-10@_19_41/4/c,rfhp0.1Hz';<br>
cfg.trialfun                = 'trialfun_general';<br>
cfg.trialdef.triallength    = 2;  % duration in seconds<br>
cfg.trialdef.ntrials        = inf;  % number of trials, inf results in as many as possible<br>
cfg                         = ft_definetrial(cfg);<br>
<br>
% preprocess the data<br>
cfg.channel                 = {'MEG','-TRIGGER','-RESPONSE',...<br>
    '-UACurrent','-MxA','-MyA','-MzA','-GxxA','-GyxA','-GyyA',...<br>
    '-GzyA','-GzxA','-MxaA','-MyaA','-MzaA','-E1','-E2','-E3','-E4'}; %read all MEG channels except those with '-'<br>
cfg.lpfilter                = 'yes';<br>
cfg.lpfreq                  = 200;<br>
cfg.bsfilter                = 'yes'; %60 Hz noise<br>
cfg.bsfreq                  = [58 62];<br>
<br>
data = ft_preprocessing(cfg);<br>
<br>
In the end, when I go to plot them all on the same graph, the head model seems to be 100 times larger than the gradiometer sensors (they appear as a small ball in the middle of the head). I've tried adjusting the grad.unit parameter from 'm' to 'cm', but it doesn't help the downstream frequency, leadfield, and source analysis.<br>



<br>
Can someone tell me how to ensure the gradiometer and mri/grid are on the same scale so I can be sure I have the right model for my beamforming?<br>
<br>
Thank you,<br>
-Allison<br></div></div>
_______________________________________________<br>
fieldtrip mailing list<br>
<a href="mailto:fieldtrip@donders.ru.nl" target="_blank">fieldtrip@donders.ru.nl</a><br>
<a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip" target="_blank">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
</blockquote>
<br>
_______________________________________________<br>
fieldtrip mailing list<br>
<a href="mailto:fieldtrip@donders.ru.nl" target="_blank">fieldtrip@donders.ru.nl</a><br>
<a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip" target="_blank">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
</blockquote></div><br>