<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Dear Ravi,<div class=""><br class=""></div><div class="">as Carsten already said, calculating FEM with anisotropic conductivities is not directly supported by the FieldTrip-SimBio implementation. However, if you are willing to invest a bit of time it is possible to work around this.</div><div class=""><br class=""></div><div class="">The „only“ thing that needs to be changed is the calculation of the FEM stiffness matrix, which is performed by the routine „calc_stiffness_matrix_val“ in the function sb_calc_stiff (usually called from ft_prepare_headmodel). The problem is that FieldTrip does not support anisotropic conductivities, so that you would have to call calc_stiffness_matrix_val directly. You can see the correct call in sb_calc_stiff. For anisotropic conductivities you have to replace the input „cond“ by a #elements x 6 matrix containing your anisotropic conductivities in the format "xx yy zz xy yz zx“. If you now follow the normal FieldTrip-SimBio workflow using the resulting stiffness matrix, you will get results for anisotropic conductivities.</div><div class=""><br class=""></div><div class="">Best,</div><div class=""><span class="Apple-tab-span" style="white-space: pre;"> </span>Johannes</div><div class=""><br class=""></div><div><blockquote type="cite" class=""><div class="">Am 29.10.2018 um 12:31 schrieb Carsten Wolters <<a href="mailto:carsten.wolters@uni-muenster.de" class="">carsten.wolters@uni-muenster.de</a>>:</div><br class="Apple-interchange-newline"><div class=""><div class="moz-cite-prefix" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; background-color: rgb(255, 255, 255);">Dear Ravi,<br class=""><br class="">1) You can use the pure SimBio-code from<span class="Apple-converted-space"> </span><br class=""><a class="moz-txt-link-freetext" href="https://www.mrt.uni-jena.de/simbio/index.php/Main_Page">https://www.mrt.uni-jena.de/simbio/index.php/Main_Page</a><br class="">to treat WM anisotropy.<br class="">While it would in principle also be possible to use anisotropic conductivities with FieldTrip-SimBio,<span class="Apple-converted-space"> </span><br class="">this is currently not implemented using ft_prepare_headmodel. Johannes (in CC), who implemented<span class="Apple-converted-space"> </span><br class="">Fieldtrip-SimBio, answered a same question by Junjie Wu in March 2018:<span class="Apple-converted-space"> </span><br class="">"Depending on your matlab skills and your available time, I could help you to give it a<span class="Apple-converted-space"> </span><br class="">try though. It should be possible with using some direct function calls instead of the high-level fieldtrip-functions."<br class=""><br class="">2) We recommend<span class="Apple-converted-space"> </span><br class=""><a class="moz-txt-link-freetext" href="http://www.sci.utah.edu/~wolters/PaperWolters/2012/RuthottoEtAl_PhysMedBiol_2012.pdf">http://www.sci.utah.edu/~wolters/PaperWolters/2012/RuthottoEtAl_PhysMedBiol_2012.pdf</a><br class="">on individual data. I could imagine that an atlas does a reasonable job w.r.t. the main<br class="">bigger fiber tracts such as corpus callosum or pyramidal tracts, but that the finer details<br class="">in the cortices are individual. We always measure T1, T2 and DTI from each subject<br class="">and I personally do not have experience with such a group-level anisotropy compared<span class="Apple-converted-space"> </span><br class="">to the individual one. Might be interesting to hear from others what they think!?<br class=""><br class="">BR<br class=""> Carsten<br class=""><br class=""><br class=""><br class="">Am 25.10.18 um 23:05 schrieb Ravi Mill:<br class=""></div><blockquote type="cite" cite="mid:BN4PR14MB059603BCC3645A4E117876C7B2F70@BN4PR14MB0596.namprd14.prod.outlook.com" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; background-color: rgb(255, 255, 255);" class=""><div id="divtagdefaultwrapper" dir="ltr" class=""><div style="margin-top: 0px; margin-bottom: 0px;" class="">Dear Fieldtrippers</div><div style="margin-top: 0px; margin-bottom: 0px;" class=""><br class=""></div><div style="margin-top: 0px; margin-bottom: 0px;" class="">I have applied the FEM simbio head modeling pipeline implemented in Fieldtrip to my EEG data. <span class="">My understanding is that this pipeline assumes isotropic conductivities for 5 head compartments (as specified by cfg.conductivity in ft_prepare_headmodel). </span>After reading some papers (e.g. Vorwerk et al 2014 <a href="https://doi.org/10.1016/j.neuroimage.2014.06.040" class="OWAAutoLink" id="LPlnk216235" previewremoved="true" moz-do-not-send="true">https://doi.org/10.1016/j.neuroimage.2014.06.040</a>), it seems like incorporating white matter conductivity anisotropy has a relatively small albeit significant effect on the source solution. I am interested in comparing FEM results when treating white matter as anisotropic. <span style="font-size: 12pt;" class="">My questions are as follows:</span></div><div style="margin-top: 0px; margin-bottom: 0px;" class=""><span style="font-size: 12pt;" class=""><br class=""></span></div><ol style="margin-bottom: 0px; margin-top: 0px;" class=""><li class="">Is there a way to implement the FEM simbio head model whilst treating WM as anisotropic within Fieldtrip? If so, how would one do this (or are there any resources available that demonstrate this)?</li><li class="">From previous papers and some simbio documentation (<a href="https://www.mrt.uni-jena.de/simbio/index.php/SIMBIO/Releasenotes/Examples" class="OWAAutoLink" id="LPlnk493580" previewremoved="true" moz-do-not-send="true">https://www.mrt.uni-jena.de/simbio/index.php/SIMBIO/Releasenotes/Examples</a>) it seems like diffusion MRI data is required to calculate the WM conductivity for each individual subject. I only have T1 and T2 scans for my subjects. So would it be possible to use WM anisotropic information obtained from some kind of diffusion MRI group average/atlas instead (accepting some loss in subject-level precision)? If so, does such a group average/atlas exist?</li></ol><br class=""><div style="margin-top: 0px; margin-bottom: 0px;" class="">Any help would be greatly appreciated!</div><div style="margin-top: 0px; margin-bottom: 0px;" class=""><br class=""></div><div style="margin-top: 0px; margin-bottom: 0px;" class="">Thanks</div><div style="margin-top: 0px; margin-bottom: 0px;" class="">Ravi</div><br class=""><p style="margin-top: 0px; margin-bottom: 0px;" class=""><span style="font-size: 12pt;" class=""></span></p><div class=""> </div></div><br class=""><fieldset class="mimeAttachmentHeader"></fieldset><pre class="moz-quote-pre" wrap="">_______________________________________________
fieldtrip mailing list
<a class="moz-txt-link-freetext" href="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip">https://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a>
<a class="moz-txt-link-freetext" href="https://doi.org/10.1371/journal.pcbi.1002202">https://doi.org/10.1371/journal.pcbi.1002202</a>
</pre></blockquote><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; background-color: rgb(255, 255, 255);" class=""><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; background-color: rgb(255, 255, 255);" class=""><pre class="moz-signature" cols="72" style="font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; word-spacing: 0px; -webkit-text-stroke-width: 0px; background-color: rgb(255, 255, 255);">--
Prof. Dr.rer.nat. Carsten H. Wolters
University of Münster
Institute for Biomagnetism and Biosignalanalysis
Malmedyweg 15
48149 Münster, Germany
Phone:
+49 (0)251 83 56904
+49 (0)251 83 56865 (secr.)
Fax:
+49 (0)251 83 56874
Email: <a class="moz-txt-link-abbreviated" href="mailto:carsten.wolters@uni-muenster.de">carsten.wolters@uni-muenster.de</a>
Web: <a class="moz-txt-link-freetext" href="https://campus.uni-muenster.de/biomag/das-institut/mitarbeiter/carsten-wolters/">https://campus.uni-muenster.de/biomag/das-institut/mitarbeiter/carsten-wolters/</a></pre></div></blockquote></div><br class=""></body></html>