Dear Fieldtrip users,<div><br></div><div>Thank you Patrick, Nathan & Gio for your inputs. My MRI anatomical images need extensive cleaning (enhance contrasts, get rid of non-zero voxels outside the brain, etc...) before segmentation, as well as adjusting cfg.threshold for each tissue type during segmentation and each subject. At least testing with SPM standard MRI, segmentation and forward modeling works OK now.</div>
<div><br></div><div>Individual subject's MRI, though, is still problematic, and I have a question for Fieldtrip developers: could you let me know whether "Create BEM headmodel for EEG (  <a href="http://fieldtrip.fcdonders.nl/example/create_bem_headmodel_for_eeg">http://fieldtrip.fcdonders.nl/example/create_bem_headmodel_for_eeg</a> )" routine is what you recommend, or it's integrated into a newer function with newer cfg options? Is there a page I missed on the web?</div>
<div><br></div><div>Thank you so much for your time!</div><div>Akiko</div><div><br><div class="gmail_quote">On Fri, Feb 10, 2012 at 4:36 AM, Gio Piantoni <span dir="ltr"><<a href="mailto:g.piantoni@nin.knaw.nl">g.piantoni@nin.knaw.nl</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Akiki and Patrick,<br>
<br>
I agree with Nathan about segmentation problems. I had the same<br>
problem when two of the three meshes intersected each other. Before<br>
you run the forward model, you should check the segmentation (bnd)<br>
using ft_plot_mesh.<br>
ft_plot_mesh(bnd(1), 'facealpha', .5)<br>
hold on<br>
ft_plot_mesh(bnd(2), 'facecolor', 'r')<br>
You will see that the two meshes intersect with each other. Depending<br>
on how you created the meshes, you should fix the parameters in such a<br>
way that they don't intersect. For example, ft_volumesegment has a<br>
cfg.threshold option (you can check how the segmentation went with<br>
ft_sourceplot) and ft_prepare_mesh_new has cfg.thresholdseg and<br>
cfg.smoothseg.<br>
<br>
I think that the meshes in SPM are computed on an (affine-transformed)<br>
template and not on the anatomically derived meshes. The template is<br>
more robust and that's why the two surfaces do not intersect.<br>
<br>
In general, the segmentation depends on the quality of your MRIs and<br>
you should really visually check the segmentation and meshes. When<br>
those look ok, the BEM forward model usually works. Hope this helps.<br>
<br>
Cheers,<br>
<br>
Gio<br>
<div class="HOEnZb"><div class="h5"><br>
On Fri, Feb 10, 2012 at 09:00, Nathan Weisz <<a href="mailto:nathanweisz@mac.com">nathanweisz@mac.com</a>> wrote:<br>
> Hi Akiki,<br>
><br>
> I am not a headmodel expert.<br>
> but looking at the output message, it seems your segementation did not go<br>
> too well (see the "intersect" message). my intuition is then then the next<br>
> temporary file *.ama is not created and ft_prepare_bemmodel then simply says<br>
> ciao.<br>
><br>
> perhaps you want to check / redo the segmentation?<br>
><br>
> good luck,<br>
> nathan<br>
><br>
><br>
> On 09.02.2012, at 23:04, Akiko Ikkai wrote:<br>
><br>
> Hi,<br>
><br>
> I'm trying to create a volume conduction model for my EEG study based on<br>
> anatomical MRI, and having trouble with ft_prepare_bemmodel<br>
><br>
> when I run<br>
> cfg                = [];<br>
> cfg.tissue         = [7 4 1]; % value for brain, skull, and scalp<br>
> cfg.conductivity   = [1 1/80 1]*.33; % after standard_BEM.zip data<br>
> cfg.isolatedsource = true;<br>
> cfg.method         = 'dipoli';<br>
> cfg.sourceunits = 'mm';<br>
> cfg.mriunits = 'mm';<br>
> vol = ft_prepare_bemmodel(cfg, seg4bem);<br>
><br>
> I get the following error message while ft_prepare_bemmodel is running:<br>
><br>
> Fatal error in dipoli:  interface<br>
> /private/tmp/tpd8a04422_e19e_4978_9b9b_56d1f2696492.tri and<br>
> /private/tmp/tp313c9a32_1991_4df6_a404_8efc41b6302e.tri intersect<br>
>  at vertex 1402 of /private/tmp/tpd8a04422_e19e_4978_9b9b_56d1f2696492.tri<br>
><br>
> Warning: an error ocurred while running dipoli<br>
>> In dipoli at 94<br>
>   In ft_prepare_bemmodel at 112<br>
> Error using ==> fread<br>
> Invalid file identifier.  Use fopen to generate a valid file identifier.<br>
> Warning: File '/private/tmp/tp0c0779e3_a657_4ca5_b9e1_845f41e9574a.ama' not<br>
> found.<br>
>> In dipoli at 102<br>
>   In ft_prepare_bemmodel at 112<br>
><br>
><br>
> and the output vol is missing .mat, which causes an error<br>
> during ft_sourceanalysis.<br>
>>> vol<br>
> vol =<br>
>             cond: [0.3300 0.0041 0.3300]<br>
>              bnd: [1x3 struct]<br>
>           source: 3<br>
>     skin_surface: 3<br>
>             skin: 1<br>
>             type: 'dipoli'<br>
><br>
> Does anyone know how to correct for "fatal error" so<br>
> that ft_prepare_bemmodel creates the correct volume conduction model? I<br>
> opened .tri file and removed line 1402, but got the same results.<br>
><br>
> Thanks in advance! Akiko<br>
><br>
> --<br>
> Akiko Ikkai, Ph.D.<br>
> Postdoctoral Fellow<br>
> Department of Psychological and Brain Sciences<br>
> Johns Hopkins University<br>
> Ames Hall, 3400 N. Charles St.<br>
> Baltimore, MD 21218<br>
><br>
><br>
> _______________________________________________<br>
> fieldtrip mailing list<br>
> <a href="mailto:fieldtrip@donders.ru.nl">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>
><br>
><br>
><br>
> _______________________________________________<br>
> fieldtrip mailing list<br>
> <a href="mailto:fieldtrip@donders.ru.nl">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>
<br>
_______________________________________________<br>
fieldtrip mailing list<br>
<a href="mailto:fieldtrip@donders.ru.nl">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></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><font><span style="font-family:arial,helvetica,sans-serif">Akiko Ikkai, Ph.D. <br>
Postdoctoral Fellow<br style="font-family:arial,helvetica,sans-serif"></span></font><font style="font-family:arial,helvetica,sans-serif" face="'PrimaSans BT,Verdana,sans-serif'" size="2">Department of 
Psychological and Brain Sciences<br>Johns Hopkins University<br>Ames 
Hall, 3400 N. Charles St.<br>Baltimore, MD 21218</font><br style="font-family:arial,helvetica,sans-serif"><br><br>
</div>