Freesurfer

freesurfer@nmr.mgh.harvard.edu

16 participants
30085 discussions

ASCII to SURFACE
by João Ricardo Sato 02 Dec '09

02 Dec '09

Dear Freesurfers, Is it possible to convert ascii files in surface files? I found some posts in this mailing list only about converting surfaces to ascii files using mris_convert. Actually, I would like to use the thickness ascii files (converted using mris_convert) to process the data in another software (such as R or Matlab), and then retrieve the results (in ascii) back to freesurfer, for visualisation in tksurfer. Could anyone help me, please? Thank you in advance! All the best João Ricardo Sato Universidade Federal do ABC Brazil

2 1

"dummy variable" in mri-glmfit
by Stefan Brauns 02 Dec '09

02 Dec '09

Hi there, we would like to test the effect of a binary variable (genotype = carrier vs. homozygous) on cortical thicknes in mri-glmfit. Since we are also controlling for gender and aquisition site (4 sites) we already have 16 groups. In order to control for age as a covariate we need at least 2 subjects per group to be able to estimate an age slope. If we include the aforementioned binary variable (genotype) as a factor (two different "groups"), we would have 32 groups and unfortunately not enough subjects per group. Is it possible to include binary variables ("dummy variable" coded as 0 and 1) such as genotype or gender as covariates (slope), in order to reduce the number of groups and examine the effect on thickness? In simple regression this would not affect the results - what would we expect here? Many thanks, Stefan

3 3

Flattening hemisphere: Lateral extents of cuts
by Dickinson, Annelise (NIH/NIMH) [F] 02 Dec '09

02 Dec '09

Hello, I'm making whole-hemisphere flat maps, and I've been referring to the 1999 paper as well as the FreeSurfer manual for where to place the cuts. However, my flat maps are not turning out like the examples and I suspect that has something to do with how far AROUND the brain the cuts should extend. I've just been making them go out to the edge, as shown in Picture 1 attached (this is Figure 4 from the 1999 Fischl et al. paper). To my knowledge, the only extant example views of the cuts are medial views, like this one. If anyone has a lateral view of the cuts (showing where they should terminate on the outer surface of the hemisphere), or an explanation of where they should go, I would really really appreciate it!! Please let me know, and thanks! -Elise ***************************************************** Annelise Dickinson Post-Bac IRTA Laboratory of Brain and Cognition, NIMH/NIH/DHHS Building 10, Room 4C212 10 CENTER DR MSC 1366 BETHESDA MD 20892-1366 ***************************************************** Phone: (301) 435-4941 FAX: (301) 402-0921 E-Mail : dickinsona(a)mail.nih.gov

3 2

Neonatal brain
by Schneider, Jacques 02 Dec '09

02 Dec '09

Hello to the freesurfer group, I would like to analyze neonatal cortices with freesurfer. I ran a test data with a 10 yo child which came through wihtout problem. But a neonatal data set failed the recon at the level of talairach transformation. I suspect that the adult template is inconsistent with the neonatal immature folding and makes the transformation impossible. Now, Kazemi et al in NeuroImage 2007 published a neonatal template, and gave me the files on a research basis. My question is: how to replace the talairach coordinates in freesurfer with these neonatal data to make it work for this population ? I suspect the answer is not easy. But we are willing to give it a try, or perhaps are you already working on it ? Cheers Jack Dr. J.F. Schneider Head of Pediatric Radiology University Children's Hospital UKBB CH-4005 Basel Switzerland Phone +41 61 685 6282 Fax +41 61 685 6015 e-mail: jacques.schneider(a)ukbb.ch Diese Mitteilung ist vertraulich. Sie kann durch das Berufsgeheimnis geschuetzte Informationen enthalten und dient nur zum Gebrauch durch den Adressaten. Sollte diese Mitteilung irrtuemlich an Sie erfolgt sein, informieren Sie uns bitte und loeschen Sie diese Mitteilung sowie allfaellige Anhaenge unverzueglich. Jedes Verbreiten, Weiterleiten oder Vervielfaeltigen dieser Mitteilung ist verboten. This message is confidential and may contain confidential information which is covered by legal or professional privilege. If you are not the intended addressee and have received this e-mail by error, please advise the sender and immediately delete this message and any attachments. Any dissemination, distribution or copying of this communication is prohibited.

2 1

Incorrect hemispheres
by Mayuresh Korgaonkar 02 Dec '09

02 Dec '09

Hello Freesurfer authors, I have processed a few datasets without realising that the hemispheric orientation was provided incorrectly i.e the labelled left hemisphere is actually the right hemisphere. Does this mean all the analysis is incorrect and unusable? What should I do to correct for this hemispheric flip? Thanks in advance, Mayuresh The INTERNET now has a personality. YOURS! See your Yahoo! Homepage. http://in.yahoo.com/

2 1

(no subject)
by Guang Zeng 01 Dec '09

01 Dec '09

Hi, there, I have made our own template using 20 subjects follwoing the instruction in FreeSurfer WIKI. (http://surfer.nmr.mgh.harvard.edu/fswiki/SurfaceRegAndTemplates) I have a quesition now, Does the quality of the new template depends on the first subject I picked to built the tempalte? Maybe I do not interpret my question clearly, I think I'd like to ask is: If we pick different subject to built a template at the fist step of the all procedure, will the output templates have big difference? If yes, is there a way in FS to pick one subject which is most similar to others, and use it at the first s tep of the build own template procedure? Thanks a lot! Guang _________________________________________________________________ Chat with Messenger straight from your Hotmail inbox. http://www.microsoft.com/windows/windowslive/hotmail_bl1/hotmail_bl1.aspx?o…

2 1

Problem with mri_aparc2aseg
by Jürgen Hänggi 01 Dec '09

01 Dec '09

Dear FS experts I run recon-all -long on 12 subjects (3-4 time points each) and there is a problem with mri_aparc2aseg. The processor is running, but the the number of labeling slices does not increase (hangs by 4 slices since hours). The other time points of this subject as well as the cross-sectional analysis worked fine. The .log file shows following: #@# WMParc Sat Nov 28 02:34:54 CET 2009 /Applications/freesurfer/subjects/thomas_sub04_tp3.long.thomas_sub04 \n mri_aparc2aseg --s thomas_sub04_tp3.long.thomas_sub04 --labelwm --hypo-as-wm --rip-unknown --volmask --o mri/wmparc.mgz --ctxseg aparc+aseg.mgz \n reading colortable from annotation file... colortable with 36 entries read (originally /autofs/space/amaebi_026/users/buckner_cortical_atlas/scripts/colortable_fin al.txt) reading colortable from annotation file... colortable with 36 entries read (originally /autofs/space/amaebi_026/users/buckner_cortical_atlas/scripts/colortable_fin al.txt) SUBJECTS_DIR /Applications/freesurfer/subjects subject thomas_sub04_tp3.long.thomas_sub04 outvol mri/wmparc.mgz useribbon 0 baseoffset 0 labeling wm labeling hypo-intensities as wm dmaxctx 5.000000 RipUnknown 1 CtxSeg /Applications/freesurfer/subjects/thomas_sub04_tp3.long.thomas_sub04/mri/apa rc+aseg.mgz Reading lh white surface /Applications/freesurfer/subjects/thomas_sub04_tp3.long.thomas_sub04/surf/l h.white Reading lh pial surface /Applications/freesurfer/subjects/thomas_sub04_tp3.long.thomas_sub04/surf/l h.pial Loading lh annotations from /Applications/freesurfer/subjects/thomas_sub04_tp3.long.thomas_sub04/label/l h.aparc.annot Reading rh white surface /Applications/freesurfer/subjects/thomas_sub04_tp3.long.thomas_sub04/surf/r h.white Reading rh pial surface /Applications/freesurfer/subjects/thomas_sub04_tp3.long.thomas_sub04/surf/r h.pial Loading rh annotations from /Applications/freesurfer/subjects/thomas_sub04_tp3.long.thomas_sub04/label/r h.aparc.annot Have color table for lh white annotation Have color table for rh white annotation Loading ribbon segmentation from /Applications/freesurfer/subjects/thomas_sub04_tp3.long.thomas_sub04/mri/rib bon.mgz Loading filled from /Applications/freesurfer/subjects/thomas_sub04_tp3.long.thomas_sub04/mri/fil led.mgz Ripping vertices labeled as unkown Ripped 6730 vertices from left hemi Ripped 6931 vertices from right hemi Building hash of lh white Building hash of lh pial Building hash of rh white Building hash of rh pial Loading aseg from /Applications/freesurfer/subjects/thomas_sub04_tp3.long.thomas_sub04/mri/ase g.mgz Loading Ctx Seg File /Applications/freesurfer/subjects/thomas_sub04_tp3.long.thomas_sub04/mri/apa rc+aseg.mgz ASeg Vox2RAS: ----------- -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; ------------------------- Labeling Slice 0 1 2 3 4 Any idea? I repeated the analysis several times and always the same problem occured. Thanks in advance Regards Jürgen ---------------------------------------------------------------------------- Jürgen Hänggi, Ph.D. Division Neuropsychology Institute of Psychology University of Zurich Binzmuehlestrasse 14, PO Box 25 8050 Zurich, Switzerland 0041 44 635 73 97 (phone office) 0041 76 445 86 84 (phone mobile) 0041 44 635 74 09 (fax office) BIN 4.D.04 (office room number) j.haenggi[at]psychologie.uzh.ch (email) http://www.psychologie.uzh.ch/neuropsy/ (website) http://www.juergenhaenggi.ch (private website) This e-mail (and any attachment/s) contains confidential and/or privileged information. If you are not the intended recipient (or have received this e-mail in error) please notify the sender immediately and destroy this e-mail. Any unauthorised copying, disclosure or distribution of the material in this e-mail is strictly forbidden. ----------------------------------------------------------------------------

2 1

colin27 skullstrip problem
by Mishkin Derakhshan 01 Dec '09

01 Dec '09

Hi, I'm trying to get the surfaces for the colin27 brain that is included with the mni packages. Has anyone managed to do this before? I thought I would have an easy time given that this linear average of 27 brains has very good signal to noise and nice contrasts, and even though recon-all completes without errors, the surfaces are way off. The main problem appears to be the skull stripping. My input volume is here: http://opus.bic.mni.mcgill.ca/public/colin27_t1_tal_lin.mnc.gz These are the steps I took and where I think the problem might be: 1. gunzip colin27_t1_tal_lin.mnc.gz 2. recon-all -i colin27_t1_tal_lin.mnc -subjid colin27 3. tkmedit conlin27 orig/001.mgz At this point the image the brightness and contrast levels of the image look very washed out. I am hoping this is only a display issue and that the underlying intensities of the image are still intact and is not a result of something going on in mri_convert, but maybe this is what is causing the downstream skull strip to fail. Can someone take a look? http://opus.bic.mni.mcgill.ca/public/001.mgz 4. recon-all -autorecon1 -subjid colin27 5. tkregister2 --mgz --s colin27 --fstal Results are not perfect but reasonable ie. about the same as many other runs I've had that have passed so I don't think it is the registration. 6. tkmedit colin27 brainmask.mgz -aux T1.mgz It is clear this is where the problem is. http://opus.bic.mni.mcgill.ca/public/brainmask.mgz http://opus.bic.mni.mcgill.ca/public/T1.mgz 7. I have subsequently tried: recon-all -skullstrip -wsthresh 10 -clean-bm -no-wsgcaatlas -subjid colin27 (basically extract the entire head) recon-all -skullstrip -wsthresh 5 -clean-bm -no-wsgcaatlas -subjid colin27 (basically extracted the entire head) recon-all -skullstrip -wsthresh 5 -clean-bm -subjid colin27 (extracted the brain, but left large chunks of skull) recon-all -skullstrip -wsthresh 2 -clean-bm -subjid colin27 (extracted the brain, but left quite a few chunks of skull) I'm concerned that the -wsthresh is so low in order to get a reasonable result, that's why I think the intensities might be off (see step 3). I'm re-running -autorecon2/3 now to see if the surfaces are better, but if someone could take a look at my input file and the converted 001.mgz to see if it isn't a problem there it would be much appreciated. thanks, mishkin

4 5

recon-all -qcache question
by Guang Zeng 01 Dec '09

01 Dec '09

Hi, there, I have 60 subjects which have been FS analyzed without the flag -qcache. Now, I wan to do group analysis using QDEC. I need run recon-all -s subjectid -qcache for each of my subjects. Just want to know how it will take to finish the recon-all -s subjectid -qcache command for one subject? Will it go through the whole -autorecon1, 2, 3 process again? Thanks a lot! Guang _________________________________________________________________ Hotmail: Trusted email with powerful SPAM protection. http://clk.atdmt.com/GBL/go/177141665/direct/01/

3 2

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

2006

2005

2004

Freesurfer