Freesurfer

freesurfer@nmr.mgh.harvard.edu

18 participants
30067 discussions

FSGD and design matrix
by Reschke, Cole 30 Nov '06

30 Nov '06

All- I am looking to run freesurfer group analysis on a longitudinal study, examining the changes between baseline and follow-up. I would like to set up a paired t-test, yet I am still in need of some assistance creating my FSGD file and my design matrix. I have read the GroupAnalysis wiki page, and did not see anything regarding t-tests. Also, is "Class" is equivalent to "Group" in the FSGD file? I have been using your sample GroupDescriptorFile to create my own, and pasted an excerpt below (my edits are in red and underlined). GroupDescriptorFile 1 Title MyTitle Class bsl plus blue Class flp circle green Variables Age Gender Input subjid1 bsl 20 M Input subjid1 flp 20 M DefaultVariable Age Cole Reschke Data Manager, PET Center Banner Alzheimer Institute Banner Good Samaritan Medical Center Phone: (602) 239-4868

2 1

segmentation fault during mri_glmfit
by Bum Seok Jeong 30 Nov '06

30 Nov '06

Dear freesurfer users I got an error message during mri_glmfit. (please see below my error message) I was able to find 'mri_glmfit.log' file in './glm/lh.dx_age.glmdir' directory Thank you for advance --------------------------------------------------- These are my constrasts. : to construct the contrast vector for testing volume change according to age = age.mat 0 0 0 1 : for effect of interest (ANOVA: f test) -> Am I right? 1 0 0 0 0 1 0 0 0 0 1 0 : for difference btw 3 group independent of age -> Am I right? -1 0 1 0 0 -1 1 0 There are three class and one variable (see below my FSGD file) # FSGD file GroupDescriptorFile 1 Title Thickness vs Dx and Age Class A circle blue Class B circle green Class C plus red Variables Age Input 61 A 71 Input 54 A 72 Input 12 B 74 Input 33 B 77 Input 63 C 80 Input 01 C 87 ...... DefaultVariable Age # Error message mri_glmfit --y lh.dx_age.thickness.10.mgh \--fsgd my_dx_age_1_fsgd.txt doss \--glmdir lh.dx_age.glmdir \--surf average lh \--C age.mat gdfReadHeader: reading my_dx_age_1_fsgd.txt INFO: ignoring tag vs INFO: gd2mtx_method is doss Reading source surface /GR/average/surf/lh.white reading group avg surface area 936 cm^2 from file Reading in average area /GR/average/surf/lh.white.avg.area.mgh Number of vertices 163842 Number of faces 327680 Total area 62167.195312 AvgVtxArea 0.379434 AvgVtxDist 0.717000 StdVtxDist 0.244431 $Id: mri_glmfit.c,v 1.83 2006/03/10 23:26:51 greve Exp $ cwd /GR/group_analysis/glm cmdline mri_glmfit --y lh.dx_age.thickness.10.mgh --fsgd my_dx_age_1_fsgd.txt doss --glmdir lh.dx_age.glmdir --surf average lh --C age.mat sysname Linux hostname *** machine x86_64 user *** FixVertexAreaFlag = 1 OneSampleGroupMean 0 y /GR/group_analysis/glm/lh.dx_age.thickness.10.mgh FSGD my_dx_age_1_fsgd.txt glmdir lh.dx_age.glmdir Creating output directory lh.dx_age.glmdir Loading y from /GR/group_analysis/glm/lh.dx_age.thickness.10.mgh Matrix condition is 177.843 search space = 93594.7 Segmentation fault

2 1

error message during CA normalize (3T image)
by Bum Seok Jeong 30 Nov '06

30 Nov '06

I'm trying to analyze 3T spgr image with freesurfer. I faced error message during autorecon2 step. And typed 'recon-all -autorecon2 -conorm -subjid case**** to figure out which process it occured in. This is error message after putting '-conorm' tag. ------------------------------------------------------------------------------------ ..... registration took 7 minutes and 3 seconds. #-------------------------------------- #@# CA Normalize Thu Nov 30 11:12:00 EST 2006 /gr/caseD00916/mri mri_ca_normalize -mask brainmask.mgz nu.mgz /freesurfer/average/RB_all_2006-02-15.gca transforms/talairach.lta norm.mgz using MR volume brainmask.mgz to mask input volume... reading 1 input volumes reading atlas from '/freesurfer/average/RB_all_2006-02-15.gca'... setting gca type = Normal gca type reading transform from 'transforms/talairach.lta'... GCAhistoScaleImageIntensities: could not find wm peak reading input volume from nu.mgz... resetting wm mean[0]: 109 --> 116 resetting gm mean[0]: 68 --> 68 input volume #1 is the most T1-like using real data threshold=39.0 using (126, 155, 152) as brain centroid... mean wm in atlas = 116, using box (126,155,136) --> (125, 154,168) to find MRI wm before smoothing, mri peak at 0 after smoothing, mri peak at 0, scaling input intensities by inf Linux 2.6.17-1.2157_FC5 #1 SMP Tue Jul 11 22:53:56 EDT 2006 x86_64 x86_64 x86_64 GNU/Linux recon-all exited with ERRORS at Thu Nov 30 11:12:16 EST 2006

2 1

ASL perfusion maps
by Abi Gopal 30 Nov '06

30 Nov '06

Hi all I have collected some raw Arterial Spin Labeling (ASL) data. I would like to use freesurfer to create perfusion maps of this data. Is it possible to take the raw data, motion correct it, perform a pairwise subtraction, and then average all the subtractions to create a mean perfusion map? Please let me know if you have ideas on how to do this sort of thing. Thanks Abi -- Abilash A. Gopal Doris Duke Research Fellow, Harvard Medical School MD Candidate 2008, Tufts Medical School

3 3

Extracting surface normals using mri_surf2surf
by XJ Kang 29 Nov '06

29 Nov '06

Hi, All, I am trying to extract the surface normal values using mri_surf2surf: %mri_surf2surf --s subject1 --hemi lh --sval-nxyz white --tval lh.white.surfn.img --trg_type analyze4d The error mesage is : ...... Loading source data Extracting surface normals INFO: surfcluster: NOT fixing group surface area INFO: trgsubject = srcsubject Saving target data Reshaping 131441 (nvertices = 131441) ANALYZE FORMAT ERROR: nslices 131441 in volume exceeds 32768 I am using FreeSurfer v303. I tested the above command on the subjects inflated by old version of FreeSurfer, and I could get the surface normals. Anything wrong with the new inflation? Thank you for your time. XJ Kang

3 3

mri_segstats - missing data in FreeSurferColorLUT.txt
by Darren Weber 29 Nov '06

29 Nov '06

Why do we see the error below from CTABread? Is it important to have 6 parameters in FreeSurferColorLUT.txt? Thanks, Darren mri_segstats --seg mri/aseg.mgz --sum stats/aseg.stats --pv mri/norm.mgz --ctab-default --excludeid 0 --brain-vol-from-seg --brainmask mri/brainmask.mgz --in mri/norm.mgz --in-intensity-name norm --in-intensity-units MR --etiv --subject ucsf_disc08 Using defalt ctab /home/dweber/freesurfer/home/FreeSurferColorLUT.txt atlas_icv = 1.3356e+06 Loading mri/aseg.mgz Loading mri/norm.mgz Loading mri/norm.mgz Loading mri/brainmask.mgz # nbrainmaskvoxels 1453903 # brainmaskvolume 1453903.0 # nbrainsegvoxels 1047830 # brainsegvolume 1047830.0 Voxel Volume is 1 mm^3 Generating list of segmentation ids CTABread(/home/dweber/freesurfer/home/FreeSurferColorLUT.txt): could not scan 6 parms from line 1100 ctx-lh-Unknown 0 0 0 (403) CTABread(/home/dweber/freesurfer/home/FreeSurferColorLUT.txt): could not scan 6 parms from line 1101 ctx-lh-Corpus_callosum 50 50 50 (403) ....

2 1

re: group analysis
by Emily Cooper 27 Nov '06

27 Nov '06

Hi Doug, Thanks! I was wondering: when I have selected "Calculate Identity Matrix" to overly functional data on the fly in the tksurfer gui before, it seems to align it perfectly to the freesurfer anatomical. Is there a way to output the matrix that it's using as a register.dat file? I was trying last week to make the file using tkregister2 and for some reason it always came out misaligned (that's why I started just using the gui). mri_vol2surf is working to make surface overlays, but they don't line up right, presumably because it's using the old register.dat file from tkregister2. mri_vol2surf --src 003/bold/03.trim_time.tval+orig.BRIK --src_type afni --srcreg register.dat --hemi lh --projfrac .3 --out trim_time.tval.surf.5 --out_type paint Emily

1 0

group level analysis
by Emily Cooper 27 Nov '06

27 Nov '06

Hi, I have functional datasets (.HEAD and .BRIK) from Afni that I'm trying to analyze in the surface domain using Freesurfer. I have Freesurfer anatomical surfaces for all my subjects, and have managed to use the tkmedit and tksurfer gui's to map my functional overlays (using the "calculate identity matrix" option) onto the volume and surface. What I'd like to do now is map the functional data onto a surface with a common number of nodes so that I can extract a matrix of beta values and do group level analysis using R. Tksurfer seems to be able to do this with the "Save plotted data to table" option (greyed over right now), however, I've read in other messages on the mail archive that the sphere.reg surface is not in common space. I've been trying to use func2sph-sess, but I seem to be missing something, getting "ERROR: cannot find any sessions." Maybe that's the wrong way to go about it? Thanks, Emily

2 1

Partial brain statistics registered to a full brain surface-based map?
by emoulton＠mclean.harvard.edu 27 Nov '06

27 Nov '06

Can a partial brain statistical image volume be accurately registered to a whole brain surface-based map? When I try to convert a partial brain statistical file (zstat1.nii.gz, created using FSL) into an overlay for a surface-based full brain (created using FreeSurfer), the result is a statistical overlay that spreads across the entire surface of the brain. My command line looks like this: mris_preproc --out lh-zstat1.mgz target SUBJ1 --hemi lh --mean --iv zstat1.nii.gz ../../../anat/Dummy.feat/reg/freesurfer/anat2exf.register.dat My register.dat transformation matrices work fine when I use tkmedit to view the cope1.nii.gz image as an overlay of the orig.mgz file (partial brain zstats where I would expect them on the image of the full brain). Does mris_preproc adher strictly to the register.dat parameters, or is it possibly over-ridden by the target surface? Thanks, Eric Any information, including protected health information (PHI), transmitted in this email is intended only for the person or entity to which it is addressed and may contain information that is privileged, confidential and or exempt from disclosure under applicable Federal or State law. Any review, retransmission, dissemination or other use of or taking of any action in reliance upon, protected health information (PHI) by persons or entities other than the intended recipient is prohibited. If you received this email in error, please contact the sender and delete the material from any computer.

3 5

vertex coordinates / volume coordinates
by Roberto Toro 27 Nov '06

27 Nov '06

Hello, I have a C code that reads FS surfaces and does some processing on them. I would like to map this back to the original MRI volume. Googling around I found that vertex coordinates may correspond to coordinates in the COR volume.... but COR volumes are now deprecated....so two questions: 1. what is the relationship between vertex coordinates and volume coordinates? 2. is there a special volume where vertex coordinates just map volume coordinates? thank you very much! roberto This message has been checked for viruses but the contents of an attachment may still contain software viruses, which could damage your computer system: you are advised to perform your own checks. Email communications with the University of Nottingham may be monitored as permitted by UK legislation.

3 2

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