Freesurfer

freesurfer@nmr.mgh.harvard.edu

4 participants
30090 discussions

content of "subjects_dir"
by Schumman Resonance 05 Jul '17

05 Jul '17

Dear free surfer community, Silly question here. I note the default "subjects_dir" in freesurfer contains a bunch of directories and files. If I set the environment so that my "subjects_dir" is not the default location, do I need to move any of the files in the original location to the new one? I hope this makes sense. Thanks, Kev

2 1

Question mri_glmfit fwhm smoothing
by Van Der Deijl, Rosanne 05 Jul '17

05 Jul '17

Hi, I have a question about mri_glmfit. As an input, you give e.g. the stacked thickness data from all subjects, smoothed with a certain fwhm. But then as an output file, it gives an fwhm.dat file, which, at least in my case, contains another smoothing level. I noticed this because I ran mri_glmfit-sim forclusterwise multiple comparison correction, which gives me the error: ERROR: cannot find /usr/local/freesurfer/stable6/average/mult-comp-cor/fsaverage/rh/cortex/fwhm33/abs/th13/mc-z.csd As input I gave data smoothed to 25fwhm, and not fwhm33, as it seems to think. But indeed, in the fwhm.dat file it gives fwhm = 32.940499. I read in the tutorial that the fwhm.dat file contains the average FWHM of the residual. But still my questions are: - what this means - why mri_glmfit-sim would need this number, and whether/how I can run the mult. comparison correction on the 25fwhm-smoothed data. Thank you in advance! - Rose Commands: mri_glmfit: mri_glmfit --fsgd $CALMGREY/CALM_FSGD.txt dods --C CALMlong_contrast1.mat --C CALMlong_contrast2.mat --C CALMlong_contrast3.mat --C CALMlong_contrast4.mat --glmdir rh.CALMlong.thickness-spc.fwhm25 --y rh.CALMlong.thickness-spc.stack.fwhm25.mgh --cortex --surf fsaverage rh mri_glmfit-sim: mri_glmfit-sim --glmdir rh.CALMlong.thickness-spc.fwhm0 --cache 1.3 neg --cwp --2spaces

2 1

mri_glmfit error (mri_reshape number of elements)
by Mcnorgan, Christopher 05 Jul '17

05 Jul '17

Hi, I’m running a group-level analysis and have encountered an error that results in a core dump. The first-level analyses were run using the fsaverage surfaces, rather than self surfaces. I’ve already inspected the individual contrasts, loading the fsaverage surface in tksurfer and rendering each subject’s contrast as an overlay, and would now like to compute the group-level results. I first ran isxconcat-sess, and saw no errors indicated in the program output. mri_glmfit was then called as followed, and produced the following output (same result for the rh): $ mri_glmfit --y ces.nii.gz --wls cesvar.nii.gz --osgm --surf fsaverage lh --glmdir glm.wls --nii.gz Reading source surface /home/chris/LDT/fsaverage/surf/lh.white Number of vertices 163842 Number of faces 327680 Total area 65416.648438 AvgVtxArea 0.399267 AvgVtxDist 0.721953 StdVtxDist 0.195470 $Id: mri_glmfit.c,v 1.196.2.8 2012/11/01 18:51:41 greve Exp $ cwd /home/chris/LDT/RFX/LDT_fsaverage.sm4.rh/hi-v-lo cmdline mri_glmfit --y ces.nii.gz --wls cesvar.nii.gz --osgm --surf fsaverage lh --glmdir glm.wls --nii.gz sysname Linux hostname brocasarea machine x86_64 user chris FixVertexAreaFlag = 1 UseMaskWithSmoothing 1 OneSampleGroupMean 1 y /home/chris/LDT/RFX/LDT_fsaverage.sm4.rh/hi-v-lo/ces.nii.gz logyflag 0 usedti 0 labelmask /home/chris/LDT/fsaverage/label/lh.cortex.label maskinv 0 glmdir glm.wls IllCondOK 0 ReScaleX 1 DoFFx 0 wFile cesvar.nii.gz weightinv 1 weightsqrt 1 Creating output directory glm.wls Loading y from /home/chris/LDT/RFX/LDT_fsaverage.sm4.rh/hi-v-lo/ces.nii.gz Saving design matrix to glm.wls/Xg.dat Normalized matrix condition is 1 Matrix condition is 1 Found 149955 points in label. ERROR: mri_reshape: number of elements cannot change nv1 = 163842, nv1 = 537168 Pruning voxels by thr: 0.000000 Segmentation fault (core dumped) This seems to have its root in a mismatch between the vertices in the fsaverage lh.cortex label and those in the contrast data. However, I carried out subject-level analyses in the fsaverage surface space specifically to facilitate a group-level analysis in fsaverage surface space, so I can’t see where I went wrong. The only other post I could find that was similar was related to an attempt to carry out a group analysis in native surface space, which is not the case here. Thanks for your time

2 1

Building your own fsaverage
by Silas 05 Jul '17

05 Jul '17

Dear FreeSurfer team, I'm making a surface based group analysis of 40 multiple sclerosis subjects using FreeSurfer tools. What are the advantages of building your own fsaverage (make_average_subject) instead of using the premade fsaverage? Would you in general recommend this? What is the pro's and con's of building your own fsaverage? All the best, Silas I searched the mailing stream and found this answer: 'make_average_subject script, but as Doug implies you probably don't need to. Our warp is high-enough dimensional that you rarely need to build average subjects from your own population'

2 1

freesurfer@nmr.mgh.harvard.edu
by 糜超 05 Jul '17

05 Jul '17

Hi, I am preprocessing my fMRI data acquired with Philips. when I using the command 'dcmunpack -src . -targ sess01 -run 101 bold nii f.nii .....' The error occurred likes below: Found 31 total files. Found 12 unique series: 0 101 301 401 501 601 701 801 901 1001 1101 1201 Subject xxx Date xxx Time 110223 Institution xxxxx. 0 Survey unknown unknown unknown unknown unknown unknown ./PS_0002 101 Survey unknown unknown unknown unknown unknown 140 ./XX_0003 301 T1W_3D_TFE_ref unknown unknown unknown unknown unknown 190 ./IM_0004 401 FE_EPI unknown unknown unknown unknown unknown 2835 ./XX_0007 501 FE_EPI_EX1 unknown unknown unknown unknown unknown 2835 ./IM_0008 601 FE_EPI_EX2 unknown unknown unknown unknown unknown 2835 ./IM_0011 701 FE_EPI_EX3 unknown unknown unknown unknown unknown 2835 ./IM_0013 801 FE_EPI_EX4 unknown unknown unknown unknown unknown 2835 ./XX_0016 901 FE_EPI_EX5 unknown unknown unknown unknown unknown 2835 ./XX_0018 1001 FE_EPI_EX6 unknown unknown unknown unknown unknown 2835 ./IM_0019 1101 FE_EPI_EX7 unknown unknown unknown unknown unknown 2835 ./XX_0022 1201 FE_EPI_EX8 unknown unknown unknown unknown unknown 2835 ./XX_0024 mri_convert ./XX_0003 ./bold/101/f.nii --nskip 0 --ndrop 0 nskip = 0 ndrop = 0 $Id: mri_convert.c,v 1.179.2.7 2012/09/05 21:55:16 mreuter Exp $ reading from ./XX_0003... Starting DICOMRead2() dcmfile = ./XX_0003 dcmdir = . WARNING: tag image orientation not found in ./XX_0003 Ref Series No = 101 ERROR: bits = 0 not supported. My os is Linux 16.04 and the Freesurfer's version is 5.3.

2 1

GLM fitting for the cortical volume
by Ali Radaideh 05 Jul '17

05 Jul '17

Dear FS experts, I have managed to run the *glm_fit* on the cortical thickness and area on a group of MS patients. However, i am not sure how to do this on the cortical volume. which file represents the cortical volume in subject's directory? Many thanks, -- *Ali M. Al-Radaideh. PhDHead of Department of Medical ImagingVice Dean of the Faculty of Allied Health SciencesThe Hashemite University, Zarqa, JordanW.phone +962 5 3903333 ext.5422, 5355, 5364Email: ali.radaideh(a)hu.edu.jo <ali.radaideh(a)hu.edu.jo>* *webpage: http://staff.hu.edu.jo/Default.aspx?id=ABZUPDzYeqM= <http://staff.hu.edu.jo/Default.aspx?id=ABZUPDzYeqM=>*

2 1

Strange Discrepancy Between Activation Location on Surface vs. Volume
by Taylor, Johnmark 05 Jul '17

05 Jul '17

Hello, I am running into an odd issue. In particular, I have run a GLM analysis in *volume* space, and I am visualizing the results both on the surface using tksurfer, and on the volume using tkmedit. In both cases, I am using the register.dof6.dat file to do the registration. Oddly, when I visualize the results on the *surface*, the activation patch falls on the medial surface of the hemisphere, in the *precuneus*, whereas when I visualize the results on the *volume*, the patch is in the *intraparietal sulcus* (which is where we expect to see this particular activation). I have made a label out of the activation in both cases and viewed it in the other space just to make sure that it is in fact the same activation, and it is. Any idea what might be going on? It seems to me that it can't be an issue with the BOLD volume --> structural volume registration, since the activations look normal on tkmedit. It seems to me, then, that it must have to do with the relationship between the anatomical volume and the anatomical surface--could this be where the problem lies? Maybe an imperfect surface reconstruction? I've attached pictures of what I mean. The activation in question is in the right hemisphere. On the surface it can be seen on the medial part of the hemisphere , and in the volume it's where the red cursor is (as I mentioned, I verified by making labels that these are in fact the same activation). Thank you very much, JohnMark

2 1

Using multimodal slab input for hippocampus subfield segmentation in FS6
by Frauke Beyer 05 Jul '17

05 Jul '17

Dear experts, I am trying to use T2-weighted images (0.5x0.5x1mm) as an additional scan to a T1-weighted image (0.7x0.7x0.7mm) for hippocampal subfield segmentation in Freesurfer 6. The T2-weighted data is a slab which includes the whole hippocampus and its y-axis is parallel to the hippocampal long axis (image dimensions [384,384,50]). When I try to use recon-all -s <subject_name> -hippocampal-subfields-T2 <file name of additional scan> <analysisID>I always get a faulty alignment (through visual inspection by overlaying T2.FSspace.mgz on T1.mgz in freeview or by looking at the quality control image T1_to_T2.v10.QC.gif), even when I manually align the images beforehand. Do you have any idea or hint how to improve the registration and thus start the segmentation correctly? I also tried to find out which registration tool is used for this (mri_robust_register?) in order to play around with it and would also be thankful for any hint on this. Thank you for your advice and time, Best, Frauke

2 1

Recon-all: too many ('s
by M del Mar Velasco 05 Jul '17

05 Jul '17

Dear all, I am running a recon-all with Freesurfer-Linux-centos6_x86_64-dev-20170330-de24995 and Ubuntu 16.04: recon-all -subjid subject_name -all I am using 0.65m3 images from RM. I get the message ' Too many ('s ' and the recon-all stops at Aparc-to-Aseg. The recon.all.log doesn't state anything about this message. The analysis just stops. I have attached the log file and the screenshot about the message. Any idea about the reason about this issue? And how can I solve it? Thanks you in advance M del Mar

2 1

error message while using autorecon1
by Buckless, Colleen 05 Jul '17

05 Jul '17

Hi All, I am using the newest version of freesurfer, and I keep getting errors when I run recon-all -autorecon1. It appears that the error might have to do with nu_correct, and I have attached the recon-all.log. Any suggestions on how to fix this? Thanks, Colleen Buckless Research Assistant Center for Neurodevelopmental and Imaging Research Kennedy Krieger Institute 443-923-2747 Buckless(a)kennedykrieger.org ---------------------------------------------------------------------- Disclaimer: The materials in this e-mail are private and may contain Protected Information. Please note that e-mail communication is not encrypted by default. You have the right to request further emails be encrypted by notifying the sender. Your continued use of e-mail constitutes your acknowledgment of these confidentiality and security limitations. If you are not the intended recipient, be advised that any unauthorized use, disclosure, copying, distribution, or the taking of any action in reliance on the contents of this information is strictly prohibited. If you have received this e-mail in error, please immediately notify the sender via telephone or return e-mail.

2 5

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