<div>Hi Robert,</div>
<div> </div>
<div>Thanks very much, it looks like you're right. While plotting the src data just came up with an empty image, BESA is supposed to output coordinates according to the Talairach atlas, so [0 0 0] should be on the anterior commisure. Plotting the structural MRI shows that [0 0 0] is clearly
<em>not</em> on the AC, so it needs to be realigned. I'll try volumerealign now.</div>
<div> </div>
<div>-Jared<br><br> </div>
<div><span class="gmail_quote">On 10/11/06, <b class="gmail_sendername">Robert Oostenveld</b> <<a href="mailto:r.oostenveld@fcdonders.ru.nl">r.oostenveld@fcdonders.ru.nl</a>> wrote:</span>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Hi Jared,<br><br>On 10 Oct 2006, at 22:38, Jared Van Snellenberg wrote:<br>> I've been analyzing the output of BESA's multiple source beamformer
<br>> (MEG<br>> data) using Fieldtrip, and it appears that there is a problem with the<br>> coordinate system. The results of volumewrite, when overlayed with<br>> the spm<br>> T1.img, appear to be substantially outside of the brain, and in oddly
<br>> shaped patterns that suggest an error in assigning statical values<br>> to the<br>> correct coordinates.<br><br>I suspect a misalignment between your functional and anatomical data.<br><br>> mri=read_fcdc_mri('
mri.img');<br><br>read_fcdc_mri does not do automatic alignment, but just returns the<br>data with the coordinate system as defined in the mri file.<br><br>Do you know what the coordinate system is to which the anatomical MRI
<br>is aligned? I.e. where is the origin [0,0,0] in anatomical terms?<br><br>> for k=1:length(files)<br>>     src(k)=besa2fieldtrip([fdir filesep files(k).name]);<br>>     interp(k)=sourceinterpolate(sinterp_cfg,src(k),mri);
<br>>     volume(k)=volumenormalise(norm_cfg,interp(k));<br>> end<br><br>And what is the coorrdinate system in which the functional source<br>data is described? I.e. at which location in the head would a dipole<br>at the coordinate [0 0 0] be?
<br><br>For sourceinterpolate the two input arguments (functional +<br>anatomical data) should be defined w.r.t. the same coordinate system.<br><br>I suggest that you use sourceplot for debugging, e.g. something like<br><br>
cfg =[]<br>cfg.funparameyter = 'pow'          % don't know how the functional<br>data in src is called, probably 'pow'<br>cfg.location      = 'interactive'  % or try [0 0 0]<br>sourceplot(cfg, src(1));           % there is no anatomy, only
<br>functional data<br><br>and<br><br>cfg =[]<br>cfg.funparameyter = []             % no functional data, only anatomy<br>cfg.location      = 'interactive'  % or try [0 0 0]<br>sourceplot(cfg, mri);<br><br>that should give you both volumes (in seperate windows) and you can
<br>compare whether they seem to be aligned. You could use the<br>volumerealign(*) function for manual realignment of the anatomical<br>MRI. Or you can use MRIcro. Also have a look at<br><a href="http://www2.ru.nl/fcdonders/fieldtrip/doku.php">
http://www2.ru.nl/fcdonders/fieldtrip/doku.php</a>?<br>id=fieldtrip:documentation:frequently_asked_questions#how_are_the_differ<br>ent_head_and_mri_coordinate_systems_defined<br><br>I hope this helps.<br>Robert<br><br>*) that is a completely new function that was added yesterdays to
<br>fieldtrip, hence you should fetch the latest version from ftp<br></blockquote></div><br><br clear="all"><br>-- <br>Jared Van Snellenberg<br>Social Cognitive Affective Neuroscience Unit<br><a href="http://scan.psych.columbia.edu">
http://scan.psych.columbia.edu</a><br>(212) 854-7858 p<br>(212) 854-3609 f<br>Department of Psychology, Columbia University<br>406 Schermerhorn Hall<br>1190 Amsterdam Avenue, Mail Code 5501<br>New York, NY 10027<br>