Hello Freesurfer Experts,
Long story short--I would like to extract BOLD values from each TR across all vertices for one subjects flattened surface Following is a brief overview of my steps and at the end you can see where I am stuck. First, after *recon-all I* followed the steps to cut the lh inflated surface, saved as a patch, ran *mris_flatten, *converted patch into asc file. Second, I used read_patch.m to extract all spatial information and a net loss of approximately 10k vertices (127k to 116k)
We had all functional images processed in FSL, the 4D file has 64x64x36x240 dimensions (voxels are 3.4x3.4x3.0). Next, co-registered the functional and anatomicals together via the following cmd: * bbregister --s cbp001_v1 --mov filtered_func_data.nii.gz --bold --init-fsl --reg dummy1.da*t. Afterwards, converted the volume to a surface using the following cmd*: mri_vol2surf --mov filtered_func_data.nii.gz --reg dummy1.dat --projfrac 0.5 --interp trilinear --hemi lh --o ./lh.func.vol2surf.mgh. *THe dimensions are 127027 x 1 x 1 x 240. Visually no problem when using* tksurfer cbp001_v1 lh inflated -patch /path to flattened patch/ -overlay /lh.func.vol2surf.mgh/ -timecourse lh.func.vol2surf.mgh;* click on the patch and shows the timecourse for that selected vertex. Using the View>Configure>overlay I can shuffle through the TRs to inspect the change in raw BOLD signal per vertex.
I have been perusing the email web server searching for how to extract the hemodynamic waveform for each vertex across the flattened surface and ultimately will be using matlab to understand how the spatial transformation is happening. As well I have all the matlab files that seemed relevant to my query (read_surf.m, read_patch.m, and readMRI.m). I was hoping that I would be able to have a text file with all the vertices (127K not the flattened 116k) in rows and each column would have the TRs; I ran this command;* mri_segstats --slabel cbp001_v1 lh /home/share/freesurfer/subjects/cbp001_v1/label/lh.cortex.label --avgwf mri_seg_stats.dat --i lh.func.vol2surf.mgh*, output was the 240 TRs as rows and seems like the average global BOLD signal in the single column corresponding with each TR. Excuse my naiveness, I just recently (1 month ago) started using freesurfer and I feel like I have exhausted as much of the information available on the FS wiki and email server. Any information or advice would be great! I put the commands just to give an idea of my workflow (most of the command lines are from the email server or the FS Wiki) and if there are any issues with my steps please let me know so I can correct them before starting the group analysis.
Best, Taha