Freesurfer

freesurfer@nmr.mgh.harvard.edu

21 participants
30070 discussions

talairach registration
by Cameron Ellis 08 Jan '07

08 Jan '07

Hi, Hopefully this is something simple. When I try to check the talairach transform for my subjects I type: tkregister2 --mgz --s SM1C2_autoreconall --fstal --surf orig and then I get this error: subject = SM1C2_autoreconall Reading lh surface orig MRISread(/space/troy/3/users/retrievalnewest/SM1C2_autoreconall/surf/lh.orig): could not open file No such file or directory /usr/local/freesurfer/stable/bin/tkregister2.bin: could not read surface orig No such file or directory orig.mgz is in the mri folder. why is it looking for lh.orig? Without the surf tag, tkregister loads the transform fine but I do not get the green surface outlines to assist with the registration. A more detailed transcription of the error is pasted below if needed. Thanks, Cameron [troy:SM1C2_autoreconall] (nmr-std-env) cd .. [troy:retrievalnewest] (nmr-std-env) pwd /autofs/space/troy_003/users/retrievalnewest [troy:retrievalnewest] (nmr-std-env) tkregister2 --mgz --s SM1C2_autoreconall --fstal --surf orig surfname set to orig INFO: no target volume specified, assuming FreeSurfer orig volume. target volume orig movable volume (null) reg file (null) LoadVol 1 $Id: tkregister2.c,v 1.49.2.4 2006/03/27 21:48:16 greve Exp $ Diagnostic Level -1 %Mon Jan 8 14:55:02 2007>>> /usr/local/freesurfer/dev/mni/bin/minctracc -clobber -debug /tmp/mritotal_2618/src_8_dxyz.mnc /usr/local/freesurfer/dev/mni/bin/../share/mni_autoreg/average_305_8_dxyz.mnc transforms/talairach.auto.xfm -transformation /tmp/mritotal_2618/src_8tmp2c.xfm -lsq9 -xcorr -model_mask /usr/local/freesurfer/dev/mni/bin/../share/mni_autoreg/average_305_8_mask.mnc -center -2.114464 54.148323 -58.055290 -step 4 4 4 -tol 0.004 -simplex 2 talairach.xfm --------------------- 0.790 -0.163 0.022 8.207; 0.189 0.845 -0.535 -58.555; 0.093 0.573 0.939 13.334; 0.000 0.000 0.000 1.000; INFO: loading target /space/troy/3/users/retrievalnewest/SM1C2_autoreconall/mri/orig.mgz INFO: changing target type to float Ttarg: -------------------- -1.000 -0.000 0.000 125.117; 0.000 -0.000 1.000 -84.596; -0.000 -1.000 0.000 84.914; 0.000 0.000 0.000 1.000; INFO: loading movable /usr/local/freesurfer/stable/average/mni305.cor.mgz INFO: changing move type to float Tmov: -------------------- -1.000 0.000 0.000 128.000; 0.000 0.000 1.000 -128.000; 0.000 -1.000 0.000 128.000; 0.000 0.000 0.000 1.000; mkheaderreg = 0, float2int = 0 ---- Input registration matrix -------- 0.790 -0.163 0.022 8.207; 0.189 0.845 -0.535 -58.555; 0.093 0.573 0.939 13.334; 0.000 0.000 0.000 1.000; subject = SM1C2_autoreconall Reading lh surface orig MRISread(/space/troy/3/users/retrievalnewest/SM1C2_autoreconall/surf/lh.orig): could not open file No such file or directory /usr/local/freesurfer/stable/bin/tkregister2.bin: could not read surface orig No such file or directory

4 4

ca normalize
by Don Hagler 08 Jan '07

08 Jan '07

I have a subject that has gotten stuck on the CA Normalize step. I have tried running it several times and it always gets stuck in the same place, running away for days if I let it, using 99.9% of the CPU. We are running RH4-x86_64-R302. Is this something that has been observed / dealt with before? This is the screen output for the CA Normalize step up to the point of getting stuck: #-------------------------------------- #@# CA Normalize Fri Jan 5 11:02:19 PST 2007 /space/emc1/4/data/MMILDB/ADNI/Containers/FREESURFERRECON_005_S_0448_20060504.085355_1/mri mri_ca_normalize -mask brainmask.mgz nu.mgz /space/monkeys/1/pubsw/packages/freesurfer/RH4-x86_64-R302/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 '/space/monkeys/1/pubsw/packages/freesurfer/RH4-x86_64-R302/average/RB_all_2006-02-15.gca'... setting gca type = Normal gca type reading transform from 'transforms/talairach.lta'... INFO: compute sample coordinates transform 1.405 -0.358 0.038 35.157; 0.193 0.687 -0.595 80.138; 0.176 0.895 1.062 -269.933; 0.000 0.000 0.000 1.000; INFO: transform used _________________________________________________________________ Fixing up the home? Live Search can help http://imagine-windowslive.com/search/kits/default.aspx?kit=improve&locale=…

3 3

Flag: -subjectid unrecognized
by Chunchun Ni 08 Jan '07

08 Jan '07

hi, I converted one analyze file to MGZ format, then try use the recon-all -autorecon-all -subjid bert (replace subjectid by myid ), but the command run nothing but Flag: -subjectid unrecognized. Anybody has some idea how to fix it? Thanks for any help! Chunchun

2 1

FS-FAST: problem with extracting BOLD time courses with a tksurfer tcl script
by raij＠nmr.mgh.harvard.edu 05 Jan '07

05 Jan '07

Hi, I am trying to extract numerical values for BOLD time courses of (a large number of) ROIs. For this purpose I use a tksurfer tcl script. Everything works fine when I only have a small number (~6) of contrasts/stimulus categories (as defined in the paradigm files), but if I increase the number of stimulus categories to 14-18, the script dies without producing anything. I am able to produce the numerical tables semi-manually with the GUI/tcl shell even for the large number of stimulus categories so I do not quite figure what gives. However, doing this manually for all our data would be impossible. Any suggestions how I could circumvent this problem? (The visualization part in the script below, that draws the actual time course display, is not necessary for producing the numerical tables, but removing that part of the script does not help - it still dies.) Any suggestions much appreciated! Also if anyone needs the full script I am happy to send you a copy - however notice that currently one of the command lines only works in the developmental environment (nmr-dev-env). Cheers, Tommi Specs: on machine ai nmr-dev-env cd /space/cuzco/12/users/raij/avml_fmri_mctogether/ setenv SUBJECTS_DIR /space/cognito/5/users/raij/subjects_mri setenv SUBJECT fsaverage #Here is the part of the subroutine that does the extraction (to run the whole thing from the master script that sets the variable values, launch run-label-FIR-h.csh in the same directory): # Set curvature for lh set curv lh.curv read_binary_curv set curvflag 1 # Clear earlier labels labl_remove_all # Clear earlier timecourses func_clear_selection # Load new label (use actual path to label file instead of $lname if running the script by itself) labl_load $lpath/$lname.label # Load new timecourse. 3 = No registration needed. # (use actual path to time course file [e.g., bold/<analysisname>/h-lh.bhdr] instead of $timecourse if running the script by itself) func_load_timecourse $timecourse 3 # Set timecourse to show average value across all vertices included in the label set gaLinkedVars(func_graph_avg_mode) 2 # Some draw command, dunno what it does really SendLinkedVarGroup graph # Select only some of the contrasts for display. Specifically, here turn off contrasts 4-6 # General format: set gGraphSetting(<color>,visible) <0|1> # where <color> is: Red/Green/Blue/Purple/Brown/Pink/Gray/LightBlue/Yellow/Orange # "Red" is condition zero (NULL), "Green" is condition one, etc. set gGraphSetting(Brown,visible) 0 set gGraphSetting(Pink,visible) 0 set gGraphSetting(Gray,visible) 0 # Adjust amplitude (y-axis, % BOLD signal change) scale from -1 to +2 .wwGraph.gwGraph.gwGraph axis config y -min -1 -max 2 # Pre-stimulus baseline correction (to turn OFF, use 0) set gbPreStimOffset 1 # Draw time course display (this command works only in nmr-DEV-env version!) Graph_SelectVerticesFromMode # Save text file listing the numerical time course values for all categories func_print_timecourse_selection $lpath/$lname-$aname-h.xls # Save time course display as postscript Graph_SaveToPS $lpath/$lname-$aname-h.ps

1 0

control points and white matter edits
by Paul Greenberg 05 Jan '07

05 Jan '07

Hi all, I'm having trouble getting accurate white and pial surfaces in inferior temporal regions on a few lower quality scans with significant intensity variation. I've added hundreds of control points and drawn in white matter in these regions when control points failed to recover unlabeled cortex, but when re-running the scripts with either -autorecon2-wm, or -autorecon2-cp, followed by -autorecon3 the surfaces are still not accurate. Are there other interventions I can to to get more accurate surfaces? If you manually add white matter why will the regenerated white matter surfaces not always extend around the hand drawn parts? Thanks, Paul

2 1

Coloring problem after autorecon2
by Niels Bergsland 05 Jan '07

05 Jan '07

Hi list, After completing autorecon1, I found the cerebellum missing and I manually editted for this case. But I found after autorecon2, the cerebellum were not filled with appropriate color at all. I have uploaded a picture to my webspace. Please see: http://www.lucidthoughts.org/download/color_failure.png I can also upload the dicoms to the drop box if this would help. Thanks! Niels Bergsland

2 1

label
by Ran Tao 05 Jan '07

05 Jan '07

Hi All, When I ran mri_annotation2label --subject UNC_4 --hemi rh --labelbase ./labels/aparc-rh --table $FREESURFER_HOME/surface_labels.txt Some labels are always skipped such as 6,7... I also tried on tutorial data, the same thing happened. Any ideas? Is surface_labels.txt a complete label list I should use? I also want to add some own labels(not extract them from annot), and cannot figure it out. How can I do this? Thank you, Ran

1 0

ages of subjecs in buckner_data-raw_dicoms.tar.gz
by Jose Luis Cantero 05 Jan '07

05 Jan '07

Hi there, we are trying to determine small atrophy changes based on novel analytic approaches. At this point of our project MPRAGE sequences were adquired in a group of MCI patients (65-85 yrs.) but we do not have healthy controls of the same age interval. For this purpose, I decided to download the file buckner_data-raw_dicoms.tar.gz, but I am unable to figure out the ages of each subject within the file. I was wondering if you wouldn´t mind to provide me with the ages for the 40 subjects contained in this file in order to test my algorithms by combining Freesurfer and Matlab functions. Thanks in advance. Best, Jose L. Cantero -- Jose Luis Cantero (Lorente), Ph.D. Laboratory of Functional Neuroscience Department of Physiology, Anatomy and Cellular Biology University Pablo de Olavide Ctra. de Utrera, Km.1 41013 Seville (Spain) Phone: +34 95 4348366 Fax: +34 95 4349151 email: jlcanlor(a)upo.es

2 1

Cortical thickness group analyses
by Dziobek, Isabel 04 Jan '07

04 Jan '07

Dear All, I have some questions in regard to the FreeSurfer group analyses (I am using the 3.0 Version), for which I would appreciate your help. I am trying to compare a patient to a control group with respect to cortical thickness (regressing out variables such as age and IQ) and also to see if associations exist in specific ROIs between cortical thickness and cognitive variables (for the groups separately). I started constructing a group descriptor file containing the two classes (patient and control) and the variable age. First, I wanted to compare the groups disregarding any other variable and so I constructed the following contrast vector: (1 -1 0), assuming doss and thinking that the ³0² meant that this column of the design matrix would not be taken into account in the model. However, this yielded the same result as running (1 -1 1). Do I have to construct a separate fsgd file only containing the classes and then a contrast vector (1 -1)? What about regressing out more than 1 variable (e.g. age and IQ), how does the contrast vector have to look like then? Also, how do I look at correlations between those variables and cortical thickness for the groups separately? When loading the scatter plots of the fsgd, the regression line for the separate groups seems to run always parallel, which does not seem very likely (see attached jpg). Are these regression lines at all? Where can I obtain the Rs in numbers? Can I, alternatively, stream out the numbers for the cortical thickness of a defined ROI to then run correlations in SPSS. Where would I find those numbers? I was also wondering about corrections for intracranial vault. Does the registration to the group¹s average brain takes care of differences in overall head size or is it advisable to put icv as a separate regressor (as is done in very few studies)? Thanks a lot for your help and happy holidays, Isabel ************************************** Isabel Dziobek, PhD Postdoctoral Research Fellow Independent Junior Research Group Neurocognition of Decision-Making Max-Planck Institute for Human Development Lentzeallee 94 14195 Berlin Phone: +49 30 82406 619 Fax: +49 30 82406 616 dziobek(a)mpib-berlin.mpg.de

2 1

mri_surf2surf resampling error
by Julien 04 Jan '07

04 Jan '07

Hello there, I'm trying to resample some thickness difference maps using mri_surf2surf. Unfortunately, I keep getting the following error: ------------------------------------------------------------------- [engram:Multispect] (nmr-std-env) mri_surf2surf --srcsubject 002_MEF_A \ --sval $SUBJECTS_DIR/thickness/002_MEF_A-002_MPRAGE_A/surf/lh.thickness_diff.w \ --src_type curv \ --trgsubject 001_MEF_A \ --tval $SUBJECTS_DIR/thickness/002_MEF_A-002_MPRAGE_A/surf/lh.thickness_diff_resampled_to_001_MEF_A.w \ --trg_type paint \ --hemi lh \ --surfreg sphere.reg Registration surface changed to sphere.reg srcsubject = 002_MEF_A srcval = /space/engram/1/users/Multispect/thickness/002_MEF_A-002_MPRAGE_A/surf/lh.thickness_diff.w srctype = curv trgsubject = 001_MEF_A trgval = /space/engram/1/users/Multispect/thickness/002_MEF_A-002_MPRAGE_A/surf/lh.thickness_diff_resampled_to_001_MEF_A.w trgtype = paint surfreg = sphere.reg srchemi = lh trghemi = lh frame = 0 fwhm-in = 0 fwhm-out = 0 Reading source surface reg /space/engram/1/users/Multispect/002_MEF_A/surf/lh.sphere.reg Loading source data Reading curvature file /space/engram/1/users/Multispect/thickness/002_MEF_A-002_MPRAGE_A/surf/lh.thickness_diff.ÔC MRISreadCurvature: could not open /space/engram/1/users/Multispect/thickness/002_MEF_A-002_MPRAGE_A/surf/lh.thickness_diff.ÔC No such file or directory ERROR: reading curvature file ------------------------------------------------------------------- From the initial command output, it appears as though the paths are input correctly, but upon reading the file it appears to be trying to open "lh.thickness_diff.ÔC" instead of "lh.thickness_diff.w" and subsequently failing. Am I doing anything incorrectly here? Thanks, Julien -- Julien Wonderlick Research Assistant, Behavioral Neuroscience Laboratory Department of Brain and Cognitive Sciences, MIT Building 46-5127 617.253.7437 jswonder(a)mit.edu

7 13

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