Freesurfer

freesurfer@nmr.mgh.harvard.edu

30132 discussions

Co-ordinate systems in FS
by Martin Juneja 04 Aug '17

04 Aug '17

Hi FreeSurfer experts, I am trying to correlate my behavioral data with whole brain cortical measures using cortical measures as IVs and behavioral data DV. Somehow the significant clusters I get from cortical thickness, area and volume are in MNI co-ordinates but significant clusters from local gyrification-behavioral data are in Talairach co-ordinates. I used identical procedures to preprocess my data but still I am getting my results in different co-ordinate systems (thickness, volume, area in MNI and gyrification in Tal). I would really appreciate any help. MJ

2 3

label2annot issue
by Redwan Maatoug 04 Aug '17

04 Aug '17

Hi all, So my issue is the following : I have 112 mgh files. (They come from the Yeo 17 Network atlas (we split the atlas in smaller regions). I have extracted the coordinates with : mri_surfcluster (playing with the thmin option) mri_vol2label and then I have used mri_label2annot with my own color lookup table to combine the 112 files into 1 file. (I have played with --noverbose, --maxstatwinner) When I try to display, I have the feeling that there is some overlapping because what tksufer dispays, is not very relevant and it is spotty. I have attached some screenshots Thank you very much for your help, Redwan

3 9

Optseq2 onset times
by Ruthie Shaffer 04 Aug '17

04 Aug '17

Hi, I’m using Optseq2 to generate stimulus / jitter schedules, and, depending on the TR used, in the output timing files I see several event onset times very close to the correct time, but slightly off (for example, something like 124.7999, instead of 124.8). I’ve used the following command to generate timing files: optseq2 --ntp 75 --tr 3.2 --psdwin 0 22.4 3.2 --ev wordPair 3.2 45 --nkeep 6 --o CRtiming.2017 --nsearch 100000 --tnullmin 0 --tnullmax 6.4 Should I be concerned about optimal onset times that are slightly off from the actual possible onset times based on the TR, or is it correct to assume it should be 124.4, instead of 124.4001 or 125.8, instead of 125.7999, for example? Thank you in advance for any help on this! Best wishes, Ruth Shaffer

2 1

Pineal gland volume?
by Rotem Saar 04 Aug '17

04 Aug '17

Dear Freesurfer experts, I wondered if its possible to extract the volume of the *pineal gland* using the basic cortical + sub-cortical segmentation (recon-all)? I could not find the volume of this ROI in the output. If the volume output for this ROI is available - please guide me where to look (maybe I missed it?) If not - does manual ROI drawing and following,volume calculation is possible using freesurfer? Thanks ! Rotem Saar-Ashkenazy Department of Brain and Cognitive Science Ben Gurion University of the Negev, Beer-Sheva, 84105 School of Social Work Ashkelon Academic College, Ashkelon, 78211 Israel

3 2

Re: [Freesurfer] output with input matrix size and referential
by pierre deman-kibleur 04 Aug '17

04 Aug '17

yes, it's what i wanted to say. On Aug 2, 2017 17:46, "Douglas Greve" <greve(a)nmr.mgh.harvard.edu> wrote: what do you mean by referential? the coordinate system of the vertex coords? On 8/2/17 4:49 AM, pierre deman-kibleur wrote: Hi, Is it possible to have the gifti output in the same referential as the orig nii file as well ? I tried mris_convert lh.white lh.white.gii -rl mri/orig/001.mgz but mris_convert doesn't have the -rl option. Regards, Pierre On Wed, Aug 10, 2016 at 2:40 PM, Iglesias, Eugenio <e.iglesias(a)ucl.ac.uk> wrote: > You can always run mri_convert with the -rl (“reslice like”) flag to > achieve this. Something like: mri_convert input.mgz output.mgz -rl > reference.mgz -odt float > If it’s a segmentation volume, you should add -rt nearest, in order to use > nearest neighbor interpolation. > Cheers, > Eugenio > > Juan Eugenio Iglesias > Translational Imaging Group > University College London > http://www.jeiglesias.com > http://cmictig.cs.ucl.ac.uk/ > > > On 10 Aug 2016, at 13:31, pierre deman <deman.pierre(a)gmail.com> wrote: > > Hello, > > I am trying to have the atlases aparc.a2009+aseg.mgz, hippocampal subfield > etc ... > in a matrix with the same size as the input of recon-all (and not the > cubic matrix 256x256x256) and with the same orientation. > > what is the best way to do it ? > > Cheers, > Pierre > > -- > DEMAN Pierre > Mobile : +33 7 82 57 80 94 <+33%207%2082%2057%2080%2094> > _______________________________________________ > Freesurfer mailing list > Freesurfer(a)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. > > > > _______________________________________________ > Freesurfer mailing list > Freesurfer(a)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. > > -- DEMAN Pierre Mobile : +33 7 82 57 80 94 <+33%207%2082%2057%2080%2094> _______________________________________________ Freesurfer mailing listFreesurfer@nmr.mgh.harvard.eduhttps://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer _______________________________________________ Freesurfer mailing list Freesurfer(a)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.

2 2

Re: [Freesurfer] [Announcement] Freesurfer v6.0 beta release and testing
by Elijah Mak 03 Aug '17

03 Aug '17

Hi ZK, Thanks. It is great news! I would like to report the following. asegstats2table seems to be acting strangely. File "/Applications/freesurfer_b6/bin/asegstats2table", line 195 print 'ERROR: subjects are not specified (use --subjects SUBJECTS)' Even typing asegstats2table --help leads to the same error message. It was not found in a prev dev built or in 5.3. Thanks. Also, not sure if it's placebo, but freeview seems to be running quicker now. Best Wishes, Elijah

4 3

Medial Surface Problem
by Bronwyn Overs 03 Aug '17

03 Aug '17

Hi mailing list, I am currently processing a longitudinal image and have come across a problem at the medial surface in the long image (step 3). As you can see in the attached image, the rh medial surface has been significantly overestimated and extends into the lh surface. This error showed up only in only 1 of the 2 long images for this subject, and both the time-point images and the base are perfectly fine. Do you know why this is happening? Kind regards, Bronwyn Overs Research Assistant Neuroscience Research Australia Margarete Ainsworth Building Barker Street Randwick Sydney NSW 2031 Australia M 0411 308 769 T +61 2 9399 1883 F +61 2 9399 1265 neura.edu.au <http://neura.edu.au/> <https://twitter.com/neuraustralia> <https://www.facebook.com/NeuroscienceResearchAustralia> <http://www.neura.edu.au/help-research/subscribe>

4 19

Univariate Linear Mixed Effects
by Worker, Amanda 03 Aug '17

03 Aug '17

Hi there, I am currently working through the linear mixed effects example for hippocampal volume on the wiki and I can't quite get my head around how the contrast matrix is set up. In the example, the design matrix has 14 columns but the contrast matrix only seems to have 5 (below), could you explain why that is please? And also what hypotheses these three different contrasts are testing? C = [zeros(3,3) [1 0 0 0 0; -1 0 1 0 0; 0 0 -1 0 1] zeros(3,6)]; At the moment I have a simple design with two groups, 4-5 time points so I have a design matrix with 4 columns (intercept, weeks, group, group*time) and the contrast I am using is [0 0 0 1] to test the interaction term. If I were to add covariates such as age and gender would I simply add two more columns to my design matrix and contrast to look like this...[0 0 0 1 0 0]? Thanks in advance for your help. Amanda

2 2

reducing a matrix
by Linda Eva Heidinger 02 Aug '17

02 Aug '17

Dear experts in freesurfer, I myself am a very beginner in freesurfer but I come with a difficult question to answer. As you can see from the script below, I am trying to obtain the Bold in CVS35MNI152 template space which hopefully has better properties for accuracy and determination of location of a voxel for targeting in use with TMS. One can see that bbregister and mri_cvs_register were used followed by mri_vol2vol to register the Bold and transform it from subject space to the above template. Unfortunately, the template space is in freesurfer space (1mm^3,256^3) . When we try to run the result that we obtain after transformation in our AFNI programs, they couldn't run on our computers because of the size of the Bold matrix being too large. I am wondering if it is possible for mri_cvs_register to use a template that is of the same size as the T1? or really we should just somehow reformulate the transformation matrix resulting from mri_cvs_register which I believe is: final_CVSmorph_tocvs_avg35_MNI152.m3z to allow for a reduced final Bold matrix from mri_vol2vol which will have the same dimensions as the subject space which the original Bold had? I appreciate the time and attention you give this question. Thanks, -Linda ? # FreeSurfer commands nmr-std-env setenv SUBJECTS_DIR <path to the directory where all MRI data for this study reside> setenv SUBJECT <my subject code> cd $SUBJECTS_DIR # First, run FS recon-all on the T1 (MEMPRAGE) anatomical as usual. recon-all -subjid $SUBJECT -i temp_${SUBJECT}/001.mgz -all -openmp 4 # Second, do co-registration at the individual level between the hires anatomical and the bold (functional) image. # In the below example the already motion corrected (and otherwise preprocessed) bold data for run 005 are in # $SUBJECTS_DIR/$SUBJECT/bold/005/005.mc.nii # Use boundary-based registration (bbregister) - this by far the most accurate tool. # Note: This is rigid 6 DOF. It is almost always best to use rigid DOF 6 for within-subject registration. bbregister --mov $SUBJECT/bold/005/005.mc.nii --frame 0 --bold --s $SUBJECT --init-fsl --reg boldregister.dat --fslmat ${SUBJECT}_anat2bold.mat # Just in case, let us make FSL type coregistration matrices as well. convert_xfm -omat ${SUBJECT}_bold2anat.mat -inverse ${SUBJECT}_anat2bold.mat # Now, we will compute a group-level coregistration from the individual brain to a template brain for both volumetric non-linear morphing. # We will use the standard brain cvs_avg35_inMNI152 as the target/template. mri_cvs_register --mov $SUBJECT --mni --outdir $SUBJECT/cvs_CVS35MNI152 --openmp 8 # Finally, we can move the bold image from the individual spac into the standard space: mri_vol2vol \ --mov $SUBJECT/bold/005/005.mc.nii \ --targ $FREESURFER_HOME/subjects/cvs_avg35_inMNI152/mri/norm.mgz \ --o Grand_Average/${SUBJECT}_CVS35MNI152_005.mc.nii \ --m3z $SUBJECT/cvs_CVS35MNI152/final_CVSmorph_tocvs_avg35_inMNI152.m3z \ --noDefM3zPath \ --interp nearest \ --reg $SUBJECT/bold/boldregister.dat \ --no-save-reg \

3 5

error with fcseed-config & fcseed-sess
by Jessica Huang 02 Aug '17

02 Aug '17

Hi freesurfer experts, I’m using the fcseed-config and fcseed-sess commands to create a seed at V1 in accordance with the FsFastFunctionalConnectivity Walkthrough <http://surfer.nmr.mgh.harvard.edu/fswiki/FsFastFunctionalConnectivityWalkth…> Steps: *STEP 4: Use fcseed-config to configure the parameters you wish to pass to your connectivity analysis. Sample command: fcseed-config -segid 1010 -fcname L_Posteriorcingulate.dat -fsd bold -mean -cfg mean.L_Posteriorcingulate.config *STEP 5: Create the FC seed for an individual fcseed-sess -s sessionid -cfg L_Posteriorcingulate.config V1 corresponds to the number 400 in $FREESRUFER_HOME/FreeSurferColorLUT.txt. The subject session I’m running is B101. The data is rest data but I named the file bold. Overall, with this information specific to my data, I ran the following commands: INPUT $ fcseed-config -segid 400 -fcname V1.dat -fsd bold -mean -cfg mean.V1.config OUTPUT fcname is V1.dat fcseed-config done INPUT $ fcseed-sess -s B101 -cfg mean.V1.config OUTPUT Logfile is /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/log/fcseed-sess.V1.dat.log.B101 Wed Jul 26 12:27:17 EDT 2017 Darwin Caramazza-Labs-Mac-Pro.local 16.6.0 Darwin Kernel Version 16.6.0: Fri Apr 14 16:21:16 PDT 2017; root:xnu-3789.60.24~6/RELEASE_X86_64 x86_64 setenv SUBJECTS_DIR /Users/caramazzalab/Desktop/Jess/Subjects cd /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir /Applications/freesurfer/fsfast/bin/fcseed-sess -s B101 -cfg mean.V1.config freesurfer-Darwin-OSX-stable-pub-v6.0.0-2beb96c Wed Jul 26 12:27:17 EDT 2017 Darwin Caramazza-Labs-Mac-Pro.local 16.6.0 Darwin Kernel Version 16.6.0: Fri Apr 14 16:21:16 PDT 2017; root:xnu-3789.60.24~6/RELEASE_X86_64 x86_64 setenv SUBJECTS_DIR /Users/caramazzalab/Desktop/Jess/Subjects cd /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir /Applications/freesurfer/fsfast/bin/fcseed-sess -s B101 -cfg mean.V1.config freesurfer-Darwin-OSX-stable-pub-v6.0.0-2beb96c mri_label2vol --seg /Users/caramazzalab/Desktop/Jess/Subjects/B201_mprage/mri/aparc+aseg.mgz --reg /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/register.dof6.lta --temp /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/template.nii.gz --fillthresh .5 --o /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/tmp.fcseed-sess.47801/seg.nii.gz --pvf /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/tmp.fcseed-sess.47801/pvf.nii.gz PVF /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/tmp.fcseed-sess.47801/pvf.nii.gz Number of labels: 0 Annot File: (null) Template Volume: /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/template.nii.gz Outut Volume: /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/tmp.fcseed-sess.47801/seg.nii.gz Registration File: /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/register.dof6.lta Fill Threshold: 0.5 Label Vox Vol: 1 ProjType: (null) ProjTypeId: 0 ProjStart: 0 ProjStop: 0 ProjDelta: 0.1 Subject: (null) Hemi: (null) UseNewASeg2Vol: 0 DoLabelStatVol 0 LabelCodeOffset 0 setenv SUBJECTS_DIR /Users/caramazzalab/Desktop/Jess/Subjects $Id: mri_label2vol.c,v 1.46 2014/12/08 21:11:54 greve Exp $ Template RAS-to-Vox: -------- -0.50000 -0.00000 -0.00000 56.00000; -0.00000 -0.00000 -0.50000 56.00000; -0.00000 0.45455 -0.00000 32.00000; -0.00000 -0.00000 -0.00000 1.00000; Template Voxel Volume: 8.8 nHits Thresh: 4.4 Loading registration from /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/register.dof6.lta regio_read_register: loading lta RegMat: -------- 0.99785 0.02664 0.05987 0.19449; -0.06077 0.03415 0.99757 -39.76084; -0.02453 0.99906 -0.03570 2.27075; 0.00000 0.00000 0.00000 1.00000; Label RAS-to-Vox: -------- -0.49893 -0.01332 -0.02994 55.90276; 0.01227 -0.49953 0.01785 54.86462; -0.02762 0.01552 0.45344 13.92689; 0.00000 0.00000 0.00000 1.00000; PVF /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/tmp.fcseed-sess.47801/pvf.nii.gz Computing PVF 8.8 mri_label2vol done mri_binarize --i /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/tmp.fcseed-sess.47801/seg.nii.gz --o /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/tmp.fcseed-sess.47801/mask.nii.gz --match 400 $Id: mri_binarize.c,v 1.43 2016/06/09 20:46:21 greve Exp $ cwd /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir cmdline mri_binarize.bin --i /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/tmp.fcseed-sess.47801/seg.nii.gz --o /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/tmp.fcseed-sess.47801/mask.nii.gz --match 400 sysname Darwin hostname Caramazza-Labs-Mac-Pro.local machine x86_64 user caramazzalab input /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/tmp.fcseed-sess.47801/seg.nii.gz frame 0 nErode3d 0 nErode2d 0 output /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/tmp.fcseed-sess.47801/mask.nii.gz Binarizing based on matching values nMatch 1 0 400 binval 1 binvalnot 0 fstart = 0, fend = 0, nframes = 1 Found 0 values in range Counting number of voxels in first frame Found 0 voxels in final mask Count: 0 0.000000 802816 0.000000 mri_binarize done mri_segstats --i /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/fmcpr.nii.gz --seg /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/tmp.fcseed-sess.47801/mask.nii.gz --id 1 --sum /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/tmp.fcseed-sess.47801/junk.sum --avgwfvol /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/tmp.fcseed-sess.47801/avgwf.mgh $Id: mri_segstats.c,v 1.121 2016/05/31 17:27:11 greve Exp $ cwd cmdline mri_segstats --i /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/fmcpr.nii.gz --seg /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/tmp.fcseed-sess.47801/mask.nii.gz --id 1 --sum /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/tmp.fcseed-sess.47801/junk.sum --avgwfvol /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/tmp.fcseed-sess.47801/avgwf.mgh sysname Darwin hostname Caramazza-Labs-Mac-Pro.local machine x86_64 user caramazzalab UseRobust 0 Loading /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/tmp.fcseed-sess.47801/mask.nii.gz Loading /Users/caramazzalab/Desktop/Jess/V1_REST_STUDY/dicomdir/B101/bold/001/fmcpr.nii.gz Voxel Volume is 8.8 mm^3 Generating list of segmentation ids Found 1 segmentations Computing statistics for each segmentation Reporting on 0 segmentations Using PrintSegStat Computing spatial average of each frame ERROR: no voxels found in segmentation I attached the error log below - any idea what happened? Is there something I can change in these commands to yield the timecourse of V1? These are all subjects 18+ years of age. What do you think? Thank you in advance, any and all feedback would be greatly appreciated. Thanks, Jess

3 9

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

2006

2005

2004

Freesurfer