Hi Sara, it looks like you are doing a functional constraint. See these
slides
http://surfer.nmr.mgh.harvard.edu/pub/docs/fs.multimodal-integration.Part2.pptx
(tutorial here
https://surfer.nmr.mgh.harvard.edu/fswiki/MultiModalTutorialV6.0)
and let me know if you have any questions
doug
On 1/29/18 11:33 AM, Makaretz, Sara Johanna wrote:
Hi FS group,
My goal is to calculate the % of voxels in a given region that are
higher than a region-specific threshold. This is what I've done so far
for each subject:
* recon T1
* coregister vol.nii to T1
* move aparc+aseg.mgz into functional space
* mri_segstats --i vol.nii --seg aparc+aseg.functionalspace.nii
--sum vol.aparc+aseg.stats to get NVoxels & mean stats
For each subject, for each cortical region (Seg1001-1035 &
Seg2001-2035), I now want to count the number of voxels that are above
region-specific thresholds as found in thresh.lst, ie.
ctx-lh-bankssts Seg1001 1.383697722
ctx-lh-caudalanteriorcingulate Seg1002 1.272858762
ctx-lh-caudalmiddlefrontal Seg1003 1.167540497
...
Do you have any advice on the best way to do this? I think I can use
mri_segstats, but am not sure if I'm looking at the right arguments
and/or if there is a more straightforward way to do this. This is what
I have so far:
foreach s (`cat s.lst`) / foreach seg (`cat thresh.lst | awk
'{print $2}'`)
set thresh = `cat thresh.lst | grep $seg | awk '{print $3}'`
mri_segstats --i vol.nii --seg aparc+aseg.functionalspace.nii --id
$seg --maskthresh $thresh --masksign pos --sum tmp.$seg.stats
...and then grep NVoxels from tmp.$seg.stats, grep total NVoxels
from vol.aparc+aseg.stats, divide/echo into summary stats file?
Thank you in advance for any advice!
Best,
Sara
___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
The information in this e-mail is intended only for the person to whom it is
addressed. If you believe this e-mail was sent to you in error and the e-mail
contains patient information, please contact the Partners Compliance HelpLine at
http://www.partners.org/complianceline . If the e-mail was sent to you in error
but does not contain patient information, please contact the sender and properly
dispose of the e-mail.
