[Simnibs-discuss] Focality with a specific threshold

[LAURENT TORLAY]

I answer to myself if anyone search about the same question : 

You have to modify mesh_io.py in this way : line 2673 
focality[i] = np.sum(weights[norm > c * peak_value]) -> replace by focality[i] = np.sum(weights[norm > c * your_threshold]) 
additionnaly, replace cuttofs by [100,100] in 'def get_focality' (line ) and 'def summary' (line 2677) 

L.T. 



[ https://www.univ-grenoble-alpes.fr/ | 
                            
                         ] 	Laurent Torlay Ingénieur - CNRS 


Bureau E 110 
33 (0)4 76 74 81 38 [ http://www.univ-grenoble-alpes.fr/ | www.univ-grenoble-alpes.fr ] Rejoignez-nous sur [ https://www.facebook.com/UGrenobleAlpes/ | Facebook ] - Suivez-nous sur [ https://twitter.com/UGrenobleAlpes | Twitter ] 
[ https://zimbra.univ-grenoble-alpes.fr/zimbra/VARIABLE-URL | VARIABLE-URL ] 


De: "LAURENT TORLAY" <laurent.torlay at univ-grenoble-alpes.fr> 
À: "discuss" <discuss at simnibs.org> 
Envoyé: Mardi 15 Novembre 2022 10:49:25 
Objet: Focality with a specific threshold 


Hello, 

First, I use Simnibs for several years and thank you for developping such a software ! 
My question : in the file fields_summary.txt, focality is measured by 'area with a field >= X% of the 99.9th percentile', how can we get a similar measurement with a specific threshold in V/m ? (moreover, why the result is a volume ? The area is multiplied by a mean cortical thickness (of the area, the whole brain ?) ) 
Thanks in advance, 
Best, 

L. Torlay 


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.drcmr.dk/pipermail/simnibs-discuss/attachments/20221121/a37b8487/attachment.html>