Freesurfer

freesurfer@nmr.mgh.harvard.edu

16 participants
30085 discussions

talairach2ras
by Ben Hamida Cammoun Leila Faiza 08 Feb '10

08 Feb '10

Hi Can you please help me to recover RAS voxel coordinates (I need to use them in the label file) from the talairach ones. I think that the transformation I should use is the one in Talairach.xfm (T), but If I copy one vertex coordinates P (MNI Tal) from tksurfer, apply the transformation: I don't have the coordinate of the same RAS vertex. Ps:I tried to add 0 0 0 1 to the transform to make it 4x4, and 1 to the point to make it 4x1 Please correct me if the transformation and the displayed values in Tksurfer are not related to the Subject transformation in Talaraich space. Thanks a lot Leila

1 0

Cannot Configure Segmentation Brush WHen Editing Subcortical Parcellations
by Jennifer Bramen 08 Feb '10

08 Feb '10

Dear Freesurfer Support Cut-and-paste the following info into your FreeSurfer problem report: --------------------------------------------------------------------- FREESURFER_HOME: /Applications/freesurfer Build stamp: freesurfer-Darwin-leopard-i686-stable-pub-v4.5.0 Kernel info: Darwin 9.8.0 i386 --------------------------------------------------------------------- Please include the following additional information in your report: 1) subject name 6002 2) the entire command-line executed tkmedit 6002 brainmask.mgz lh.white -aux T1.mgz -aux-surface rh.white - segmentation aseg.mgz $FREESURFER_HOME/FreeSurferColorLUT.txt & 3) the error message generated n/a I am trying to edit the subcortical parcellations. I found the wiki instructions, and have followed (I think) to a t. I am unable to effect the Segmentation Brush label by using the Configure Segmentation Brush menu. But regardless of my selection, my drawing is labeled "out of bounds". I have updated to the most recent version of Freesurfer. I also attempted to use an older version, in case it was a recent bug, incompatible with recon-all done in an older version, etc. However, no version changes have helped. Can someone tell me what I am missing? FYI - When I select a label type, I have tried both having the "Segmenmtation Brush Info" window open, with my slection highlighted, and highlighting and closing the window. I am hoping I am simply missing a step though obvious to the wiki creator. Here is the link to the "Editing" instructions I have been following: http://surfer.nmr.mgh.harvard.edu/fswiki/TkMeditGuide/TkMeditWorkingWithDat… 4) optionally include the subject's /script/recon-all.log (n/a) Thanks for your assistance. Best regards, Jennifer Bramen, Ph.D. Assistant Research Neuroscientist UCLA Center for Cognitive Neuroscience Developmental Cognitive Neuroimaging Suite 17-369 740 Westwood Plaza Los Angeles, CA 90095 Phone: (310)525-0865 Fax: (310)206-5518 Email: jbramen(a)mednet.ucla.edu

1 0

Re: [Freesurfer] Error mris_thickness_diff
by Martin Reuter 08 Feb '10

08 Feb '10

Hi Jürgen, the transforms found in the base are for each timepoint to the base and vice-versa. You'd need the registrations between the time points, so for example: tp2_to_tp1.long.base.lta from the tp1.long.base/mri/transform directory (registering tp2 to tp1). >From my test-retest analysis I see that the lta is more accurate than the -register (internal surface registration) of mris_thickness_diff. But I only tested this on a 15 subjects test population. Anyway, there is a third way of doing the thickness comparison for longitudinally processed data (and this is what I would recommend, also based on my test with the test-retest data set): As the surfaces in the longitudinal stream are in vertex correspondence across time, you do not need to spatially register at all. Just smooth, subtract, average etc the time points surface measures directly. So you can e.g. modify Han's method for step 2 similar to this: mris_thickness_diff -src_type curv \ -out $SUBJECTS_DIR/1001/surf/lh.thickness_diff.mgh \ -nsmooth 50 \ -abs -S 1001 \ -L $SUBJECTS_DIR/1001/surf/lh_thickness_diff.log \ $SUBJECTS_DIR/1001/surf/lh.white \ $SUBJECTS_DIR/1001/surf/lh.thickness \ $SUBJECTS_DIR/1001/surf/lh.white \ <------ ! $SUBJECTS_DIR/1001a/surf/lh.thickness Note that here I specify the same lh.white surface for both time points. Therefore the registration is skipped and only smoothing etc. applied. Alternatively you can do all this directly with other freesurfer commands (mris_calc ...). Han's steps 3 and 4 deal with Group analysis and should stay the same independently of cross sectional or longitudinal processing. Cheers, Martin On Mon, 2010-02-08 at 18:41 +0100, Jürgen Hänggi wrote: > Hi Martin > > Thanks a lot. In the past, I used the -xform /../mri/transforms/.. .lta, but > then I saw an email from Doug where he recommended to use -register for > longitudinal data. > However if I use lta, which transform do I have to use? > > Base_to_tp1.lta > Base_to_tp2.lta > tp1_to_base.lta > tp2_to_base-lta > > From Base/mri/transform > > Or any of the numerous transforms found in the tp1.long or tp2.long folder? > > In older versions, I think, there was a kind of combined transform that have > to be used. > > Do I'am right that when using longitudinal processed data, only step2 of > Han's four steps described under > > http://surfer.nmr.mgh.harvard.edu/fswiki/ThicknessRepeatibility > > have to be run? > > Thanks in advance > Cheers > Jürgen > > > > On [DATE], "Martin Reuter" <[ADDRESS]> wrote: > > > Hi Jürgen, > > > > I have seen this error before. Not sure exactly what is going on. There > > could be a bug in Han's code, but I could not find it quickly. > > I don't recommend the rigid alignment in mris_thickness_diff anyway. I > > have not seen any evaluation on the accuracy of the registration and do > > not trust it (it is a rigid surface based alignment using the nearest > > neighbor approach). > > > > Instead try to specify the lta (transform) for the two corresponding > > volumes. If you use the longitudinal stream, these transforms should be > > available (either in the <base>/mri/transform or in the > > <longitudinal-tp>/mri/transform dir). If you are processing your data > > cross sectionally: > > I would recommend to do it longitudinally. > > If that is not an option use mri_robust_register to register the > > norm.mgz and use that transform for the mris_thickness_diff. > > > > Cheers, Martin > > > > > > > > On Mon, 2010-02-08 at 08:30 +0100, Jürgen Hänggi wrote: > >> Dear FS experts > >> > >> when running mris_thickness_diff the following error occurred: > >> > >> mris_thickness_diff \ > >> ? -src_type curv \ > >> ? -out amblyo_pat01_tp1.long.amblyo_pat01_base1_3/surf/lh.thickness_diff.mgh > >> \ > >> ? -register \ > >> ? -abs -S amblyo_pat01_tp3.long.amblyo_pat01_base1_3 \ > >> ? -L amblyo_pat01_tp1.long.amblyo_pat01_base1_3/surf/lh_thickness_diff.log \ > >> ? amblyo_pat01_tp3.long.amblyo_pat01_base1_3/surf/lh.white \ > >> ? amblyo_pat01_tp3.long.amblyo_pat01_base1_3/surf/lh.thickness \ > >> ? amblyo_pat01_tp1.long.amblyo_pat01_base1_3/surf/lh.white \ > >> ? amblyo_pat01_tp1.long.amblyo_pat01_base1_3/surf/lh.thickness > >> Output differences to file > >> amblyo_pat01_tp1.long.amblyo_pat01_base1_3/surf/lh.thickness_diff.mgh > >> Do a rigid alignment of two surfaces before mapping to each other > >> logging results to > >> amblyo_pat01_tp1.long.amblyo_pat01_base1_3/surf/lh_thickness_diff.log > >> Reading first surface file > >> amblyo_pat01_tp3.long.amblyo_pat01_base1_3/surf/lh.white > >> Surface1 has 137991 vertices > >> Reading in first data file > >> amblyo_pat01_tp3.long.amblyo_pat01_base1_3/surf/lh.thickness > >> Reading second surface file > >> amblyo_pat01_tp1.long.amblyo_pat01_base1_3/surf/lh.white > >> Surface2 has 137991 vertices > >> Reading in second data file > >> amblyo_pat01_tp1.long.amblyo_pat01_base1_3/surf/lh.thickness > >> Register surface 2 to surface 1 (rigid mapping using ICP) > >> > >> ERROR: FindClosest: MHTfindClosestVertexNo returned annIndex=-1, which is > >> invalid! > >> > >> What does it mean? > >> > >> > >> 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. > >> ---------------------------------------------------------------------------- > >> > >> > >> > >> > >> _______________________________________________ > >> Freesurfer mailing list > >> Freesurfer(a)nmr.mgh.harvard.edu > >> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > >> > >> > > > > > ---------------------------------------------------------------------------- > 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. > ---------------------------------------------------------------------------- > > > > >

1 0

qdec problem
by Hannah 08 Feb '10

08 Feb '10

Hi, I have a problem using qdec, which I don't know how to solve- I get the following error message: Error in Analyze: command failed: mri_glmfit --y /Users/Hannah/ Borderline/subjects/qdec/Untitled/y.mgh --fsgd /Users/Hannah/ Borderline/subjects/qdec/Untitled/qdec.fsgd dods --glmdir /Users/ Hannah/Borderline/subjects/qdec/Untitled --surf fsaverage lh --C / Users/Hannah/Borderline/subjects/qdec/Untitled/contrasts/lh-Avg- Intercept-thickness.mat --C /Users/Hannah/Borderline/subjects/qdec/ Untitled/contrasts/lh-Avg-thickness-alter-Cor.mat --C /Users/Hannah/ Borderline/subjects/qdec/Untitled/contrasts/lh-Diff-BPD-KG-Intercept- thickness.mat --C /Users/Hannah/Borderline/subjects/qdec/Untitled/ contrasts/lh-Diff-BPD-KG-Cor-thickness-alter.mat --C /Users/Hannah/ Borderline/subjects/qdec/Untitled/contrasts/lh-Diff-NonPTSD-PTSD- Intercept-thickness.mat --C /Users/Hannah/Borderline/subjects/qdec/ Untitled/contrasts/lh-Diff-NonPTSD-PTSD-Cor-thickness-alter.mat --C / Users/Hannah/Borderline/subjects/qdec/Untitled/contrasts/lh-X- diagnostik-PTSD-Intercept-thickness.mat --C /Users/Hannah/Borderline/ subjects/qdec/Untitled/contrasts/lh-X-diagnostik-PTSD-Cor-thickness- alter.mat Do you have any suggestions? Thanks a lot! hannah ----------------------------------------------------------------------------------------------- Dr. Hannah Bruehl Freie Universität Berlin Emotionspsychologie und affektive Neurowissenschaft Habelschwerdter Allee 45 - Raum JK 25/222h 14195 Berlin +49 30 838 564 72 -----------------------------------------------------------------------------------------------

2 1

Error mris_thickness_diff
by Jürgen Hänggi 08 Feb '10

08 Feb '10

Dear FS experts when running mris_thickness_diff the following error occurred: mris_thickness_diff \ ? -src_type curv \ ? -out amblyo_pat01_tp1.long.amblyo_pat01_base1_3/surf/lh.thickness_diff.mgh \ ? -register \ ? -abs -S amblyo_pat01_tp3.long.amblyo_pat01_base1_3 \ ? -L amblyo_pat01_tp1.long.amblyo_pat01_base1_3/surf/lh_thickness_diff.log \ ? amblyo_pat01_tp3.long.amblyo_pat01_base1_3/surf/lh.white \ ? amblyo_pat01_tp3.long.amblyo_pat01_base1_3/surf/lh.thickness \ ? amblyo_pat01_tp1.long.amblyo_pat01_base1_3/surf/lh.white \ ? amblyo_pat01_tp1.long.amblyo_pat01_base1_3/surf/lh.thickness Output differences to file amblyo_pat01_tp1.long.amblyo_pat01_base1_3/surf/lh.thickness_diff.mgh Do a rigid alignment of two surfaces before mapping to each other logging results to amblyo_pat01_tp1.long.amblyo_pat01_base1_3/surf/lh_thickness_diff.log Reading first surface file amblyo_pat01_tp3.long.amblyo_pat01_base1_3/surf/lh.white Surface1 has 137991 vertices Reading in first data file amblyo_pat01_tp3.long.amblyo_pat01_base1_3/surf/lh.thickness Reading second surface file amblyo_pat01_tp1.long.amblyo_pat01_base1_3/surf/lh.white Surface2 has 137991 vertices Reading in second data file amblyo_pat01_tp1.long.amblyo_pat01_base1_3/surf/lh.thickness Register surface 2 to surface 1 (rigid mapping using ICP) ERROR: FindClosest: MHTfindClosestVertexNo returned annIndex=-1, which is invalid! What does it mean? 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

Segmentation fault?
by Raghav Puranmalka 08 Feb '10

08 Feb '10

Hi all, I am running recon-all as follows: recon-all -subjid subj1 -all and (after about 4 hours of running,) get this error: 0516: dt=0.000000, rms=0.507 (0.000%), neg=0, invalid=766 Segmentation fault ERROR: mri_ca_register with non-zero status 139 but continuing despite the error #-------------------------------------- #@# CA Reg Inv Sat Feb 6 02:56:03 EST 2010 /usr/local/freesurfer/subjects/be496/mri mri_ca_register -invert-and-save transforms/talairach.m3z Loading, Inverting, Saving, Exiting ... Reading transforms/talairach.m3z ERROR: cannot find or read transforms/talairach.m3z ERROR: mri_ca_register with non-zero status 0 It seems that perhaps this is because it cannot find "transforms/talairach.m3z" - indeed, this file doesn't exist. But should this be something generated by recon-all? Should I be running recon-all differently? Thanks for your help, Raghav

2 1

longitudinal design
by Casper van Oers 07 Feb '10

07 Feb '10

Hi, I am trying to specify a longitudinal design in an fsgd file. there is 2 groups and each subject has 3 timepoints. can I use a similar approach as paired-diff in the tutorial, but now using 3 timepoints like this: GroupDescriptorFile 1 Title ALL_subjects Class Stroke Class Control Input subject1 Control Input subject1_time2 Control Input subject1_time3 Control Input subject2 Stroke Input subject2_time2 Stroke Input subject2_time3 Stroke ... as I fear the approach above will not work, is there another way? thanks, casper -- Casper A.M.M. van Oers, MD Department of Neurology, C03.236 University Medical Centre PO Box 85500 3508 GA Utrecht the Netherlands T: +31 30 2503184 F: +31 30 2542100 (Dept. of neurology)

1 0

question, how to copy the liscence.
by jijun wang 06 Feb '10

06 Feb '10

Dear expert, I can not find the file of ".license" in the Freesurfer folder when I finishing the installation. who could give me some hints? Thanks a lot. -- Best regards, Jijun Wang, MD & PhD Shanghai Mental Health Center 600 Wanping Road Shanghai 20030, China. E-mail: jijunwang27(a)gmail.com

2 1

user specified output orientation using mri_convert?
by Michael Harms 05 Feb '10

05 Feb '10

Hello Doug, Is it possible using mri_convert to just reorder the way in which the data is stored, without any resampling of the actual voxel values? i.e., I have Siemen's MPRAGE DICOMs, and if run mri_convert, I get a default "PIL" orientation. I would like to just reorder the way in which the data is stored, to get a "LAS" orientation instead. I tried using "--out_orientation LAS", but that does a resampling (interpolation) such that the resulting xform has strictly 0's on the off-diagonal. And reading the mri_convert help, it clearly states that the --out_orientation flag is to be used when the input orientation is WRONG, which is not the case here (since I presume that the input orientation is correctly read out of the DICOMs). I tried "--reorder 2 1 3", to see what effect that has on a file already in nii format, but I got the message: in nii format, but I got the message: reordering axes... MRIreorder: incorrect dst width error reordering axes I've checked the archives, and read through the mri_convert help, but just can't seem to figure out if reordering from one orientation to another (without resampling) is possible using mri_convert. thanks, -MH -- Michael Harms, Ph.D. -------------------------------------------------------------------- Conte Center for the Neuroscience of Mental Disorders Washington University School of Medicine Department of Psychiatry, Box 8134 Renard Hospital, Room 6604 Tel: 314-747-6173 660 South Euclid Ave. Fax: 314-747-2182 St. Louis, MO 63110 Email: mharms(a)wustl.edu --------------------------------------------------------------------

2 5

Difference in sub-cortical volumes
by Walimuni, Indika S 05 Feb '10

05 Feb '10

Hi I noticed in my FreeSurfer runs, that there is a significant volume difference for sub-cortical structures obtained using the aseg.stat file and the mri_binaraize command. For example I got 14454mm3 for Thalamus-Proper (obtained by summing up the values for Left_Thalamus-Proper and Right-Thalamus-Proper) using aseg.stats file. For the same subject I got 15091mm3 using the command " mri_binarize --i aseg.mgz --match 10 --match 49 --o Thalamus-Proper_mask.nii.gz". Is this normal ? If so what is the most accurate method ? Thanks in advance, Indika Walimuni Research Associate University of Texas Health Science Center Houston, TX.

2 1

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