Hi Doug
Thank you very much for your reply.
I am doing these analysis in each subject’s self space, and my ROI is a
handmade cortical label in ascci format.
In this case, would the indices be the vertex Id in the label file?
also, instead of f=MRIread('f.nii.gz’); can I use “sig.nii.gz” to get the
threshold of each voxel (sig.nii.gz is the result of selxavg for a specific
contrast)?
Thanks
Mona
> On Mar 20, 2020, at 10:49 AM, Douglas N. Greve <dgr...@mgh.harvard.edu> wrote:
>
> Do you mean all the voxels separately? I don't think I have anything,
> but you can do it in matlab with something like
> seg = MRIread('seg.nii.gz');
> indices= find(seg.vol == IndexOfYourROI);
> f = MRIread('f.nii.gz');
> fmat = fast_vol2mat(f);
> froi = fmat(:,indices);
>
> this assumes taht f.nii.gz and seg.nii.gz are in vox-for-vox registration
>
>
>
> On 3/18/2020 5:13 PM, Nasiriavanaki, Zahra wrote:
>> Hi all
>>
>> Is there any command for extracting the fmri activation from all voxels in
>> an ROI?
>> I know I can use funcroi-config, funcroi-sess and funcroi-table-sess to
>> extract the average activation in a specific ROI from a specific contrast.
>> But what if I want to extract the values of each voxel within that ROI?
>>
>> I appreciate any comments on this.
>>
>> Thanks
>> Mona
>> _______________________________________________
>> 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
_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
