Freesurfer

freesurfer@nmr.mgh.harvard.edu

4 participants
30076 discussions

Desikan Atlas in MNI Space
by ece ulug 11 Oct '18

11 Oct '18

External Email - Use Caution Dear All, I am a master student in Turkey. For my thesis; I have to do exactly this; Desikan atlas (which I have in the tailarach space) in the MNI space, as a file (Desikan_mni.nii). I am aware that transforming from one space to the other is not a very good idea but nevertheless I would like to test some computations. Anyone can attach/send a link/know of if exists such a file? the aparc+aseg.mgz in fsaverage is in mni305 space. will that work? This is a mail dated 2014, freesurfer maillist, and I find this via internet. Desikan Atlas 'aparc+aseg.mgz' or .nii format where should I look for these file? And what is your opinion, is it good idea? (To transform Desikan atlas to mni space for visualization VBM results?) Thanks in advance E. Ece Uluğ

3 6

Hippocampal/Amygdala subfield segmentation error - Freesurfer devel-20180612
by Sally Grace 11 Oct '18

11 Oct '18

External Email - Use Caution Dear experts, I am running the hippocampal/amygdala segmentation on T1's preprocessed in -reconall in FS 5.3 using FS devel-20180612 using the command segmentHA_T1.sh bert [SUBJECTS_DIR] This has worked successfully for roughly 1000 participants but for one participant I keep getting the error below, despite several troubleshooting attempts. The error appears to be "Fitting mesh to synthetic image resulted in no deformation" - I have had a look at the hippocampus masks and they appear to be registered incorrectly, despite the grey matter volume output from -reconall having no errors. Do you have any idea what I can do to troubleshoot or fix this? Thanks in advance, Sally. Error Log: ------------------------------ USER valentil HOST m3a007 PROCESSID 29440 PROCESSOR x86_64 OS Linux Linux m3a007 3.10.0-514.26.1.el7.x86_64 #1 SMP Thu Jun 29 16:05:25 UTC 2017 x86_64 x86_64 x86_64 GNU/Linux ------------------------------ #-------------------------------------------- #@# Hippocampal Subfields processing (T1) left Tue Sep 25 10:18:56 AEST 2018 ------------------------------------------ Setting up environment variables --- LD_LIBRARY_PATH is .:/usr/local/freesurfer/devel-20180612/MCRv84//runtime/glnxa64:/usr/local/freesurfer/devel-20180612/MCRv84//bin/glnxa64:/usr/local/freesurfer/devel-20180612/MCRv84//sys/os/glnxa64:/usr/local/freesurfer/devel-20180612/lib/tcltktixblt:/usr/local/freesurfer/devel-20180612/lib/qt/lib:/usr/local/freesurfer/devel-20180612/lib/petsc/lib:/usr/local/freesurfer/devel-20180612/lib/KWWidgets/lib/KWWidgets:/usr/local/freesurfer/devel-20180612/lib/cuda/lib64:/usr/local/freesurfer/devel-20180612/lib/cuda/lib:/usr/local/freesurfer/devel-20180612/lib/bem:/usr/local/freesurfer/devel-20180612/lib:/usr/local/freesurfer/devel-20180612/lib/vtk/lib/vtk-5.6:/usr/local/fsl/5.0.11/fsl/lib:/usr/local/libjpeg-turbo/1.4.2/lib64:/usr/local/cuda/7.5/extras/CUPTI/lib64:/usr/local/cuda/7.5/lib64:/usr/local/cuda/7.5/lib:/usr/local/cuda/7.5/lib64/stubs:/opt/munge-0.5.11/lib:/opt/slurm-17.11.4/lib:/opt/slurm-17.11.4/lib/slurm:/usr/local/tigervnc/1.8.0/lib64:/usr/local/tigervnc/1.8.0/lib:/opt/munge-0.5.11/lib:/opt/slurm-17.11.4/lib:/opt/slurm-17.11.4/lib/slurm::/usr/local/freesurfer/devel-20180612/MCRv84//sys/opengl/lib/glnxa64:/usr/local/freesurfer/devel-20180612/lib/tcltktixblt:/usr/local/freesurfer/devel-20180612/lib/qt/lib:/usr/local/freesurfer/devel-20180612/lib/petsc/lib:/usr/local/freesurfer/devel-20180612/lib/KWWidgets/lib/KWWidgets:/usr/local/freesurfer/devel-20180612/lib/cuda/lib64:/usr/local/freesurfer/devel-20180612/lib/cuda/lib:/usr/local/freesurfer/devel-20180612/lib/bem:/usr/local/freesurfer/devel-20180612/lib:/usr/local/freesurfer/devel-20180612/lib/vtk/lib/vtk-5.6:/usr/local/fsl/5.0.11/fsl/lib:/usr/local/libjpeg-turbo/1.4.2/lib64:/usr/local/cuda/7.5/extras/CUPTI/lib64:/usr/local/cuda/7.5/lib64:/usr/local/cuda/7.5/lib:/usr/local/cuda/7.5/lib64/stubs:/opt/munge-0.5.11/lib:/opt/slurm-17.11.4/lib:/opt/slurm-17.11.4/lib/slurm:/usr/local/tigervnc/1.8.0/lib64:/usr/local/tigervnc/1.8.0/lib:/opt/munge-0.5.11/lib:/opt/slurm-17.11.4/lib:/opt/slurm-17.11.4/lib/slurm: Registering imageDump.mgz to hippocampal mask from ASEG $Id: mri_robust_register.cpp,v 1.77 2016/01/20 23:36:17 greve Exp $ --mov: Using imageDump.mgz as movable/source volume. --dst: Using /home/valentil/kg98/Valentina/ENIGMA_sex_differences/Scripts/rerun_reconall/NIAAA_MRI00032_1-out/tmp/hippoSF_T1_v21_left//hippoAmygBinaryMask_autoCropped.mgz as target volume. --lta: Output transform as trash.lta . --mapmovhdr: Will save header adjusted movable as imageDump_coregistered.mgz ! --sat: Using saturation 50 in M-estimator! reading source 'imageDump.mgz'... reading target '/home/valentil/kg98/Valentina/ENIGMA_sex_differences/Scripts/rerun_reconall/NIAAA_MRI00032_1-out/tmp/hippoSF_T1_v21_left//hippoAmygBinaryMask_autoCropped.mgz'... Registration::setSourceAndTarget(MRI s, MRI t, keeptype = TRUE ) Type Source : 0 Type Target : 3 ensure both FLOAT (3) Reordering axes in mov to better fit dst... ( -1 3 -2 ) MRIreorder() ----------- xdim=-1 ydim=3 zdim=-2 src 131 241 99, 0.250000 0.250000 0.250000 dst 131 99 241, 0.250000 0.250000 0.250000 Determinant after swap : 0.015625 Mov: (0.25, 0.25, 0.25)mm and dim (131, 99, 241) Dst: (1, 1, 1)mm and dim (36, 31, 52) Asserting both images: 1mm isotropic - reslicing Mov ... -- changing data type from 0 to 3 (noscale = 0)... -- Original : (0.25, 0.25, 0.25)mm and (131, 99, 241) voxels. -- Resampled: (1, 1, 1)mm and (36, 31, 61) voxels. -- Reslicing using cubic bspline MRItoBSpline degree 3 - reslicing Dst ... -- Original : (1, 1, 1)mm and (36, 31, 52) voxels. -- Resampled: (1, 1, 1)mm and (36, 31, 61) voxels. -- Reslicing using cubic bspline MRItoBSpline degree 3 Registration::computeMultiresRegistration - computing centroids - computing initial transform -- using translation info - Get Gaussian Pyramid Limits ( min size: 16 max size: -1 ) - Build Gaussian Pyramid ( Limits min steps: 0 max steps: 0 ) - Build Gaussian Pyramid ( Limits min steps: 0 max steps: 0 ) - initial transform: Ti = [ ... 1.0000000000000 0 0 0.8149113483238 0 1.0000000000000 0 0.6091444604181 0 0 1.0000000000000 -10.4099368843102 0 0 0 1.0000000000000 ] - initial iscale: Ii =1 Resolution: 0 S( 36 31 61 ) T( 36 31 61 ) Iteration(f): 1 -- diff. to prev. transform: 8.43588 Iteration(f): 2 -- diff. to prev. transform: 6.13157 Iteration(f): 3 -- diff. to prev. transform: 5.20397 Iteration(f): 4 -- diff. to prev. transform: 5.47621 Iteration(f): 5 -- diff. to prev. transform: 8.03452 max it: 5 reached! - final transform: Tf = [ ... 0.9476591624290 -0.1857780011274 0.2596702642998 -6.9596494990127 0.2461981066435 0.9430226622820 -0.2238185664945 2.3191313534194 -0.2032943780569 0.2760340426825 0.9394022584233 -11.1737673430040 0 0 0 1.0000000000000 ] - final iscale: If = 1 ********************************************************** * * WARNING: Registration did not converge in 5 steps! * Problem might be ill posed. * Please inspect output manually! * ********************************************************** Final Transform: Adjusting final transform due to initial resampling (voxel or size changes) ... M = [ ... -0.2369147906072 0.0649175660749 0.0464445002818 20.3444784868448 -0.0615495266609 -0.0559546416236 -0.2357556655705 36.6877108257354 0.0508235945142 0.2348505646058 -0.0690085106706 -14.6334717181848 0 0 0 1.0000000000000 ] Determinant : -0.015625 writing output transformation to trash.lta ... converting VOX to RAS and saving RAS2RAS... mapmovhdr: Changing vox2ras MOV header (to map to DST) ... To check aligned result, run: freeview -v /home/valentil/kg98/Valentina/ENIGMA_sex_differences/Scripts/rerun_reconall/NIAAA_MRI00032_1-out/tmp/hippoSF_T1_v21_left//hippoAmygBinaryMask_autoCropped.mgz imageDump_coregistered.mgz Registration took 0 minutes and 1 seconds. Thank you for using RobustRegister! If you find it useful and use it for a publication, please cite: Highly Accurate Inverse Consistent Registration: A Robust Approach M. Reuter, H.D. Rosas, B. Fischl. NeuroImage 53(4):1181-1196, 2010. http://dx.doi.org/10.1016/j.neuroimage.2010.07.020 http://reuter.mit.edu/papers/reuter-robreg10.pdf Highly Accurate Inverse Consistent Registration: A Robust ...<http://reuter.mit.edu/papers/reuter-robreg10.pdf> reuter.mit.edu Highly Accurate Inverse Consistent Registration: A Robust Approach Martin Reutera,b,c,∗, H. Diana Rosasa,b, Bruce Fischla,b,c aMassachusetts General Hospital / Harvard Medical School, Boston, MA, USA bMartinos Center for Biomedical Imaging, 143 13th Street, Charlestown, MA, USA cMIT Computer Science and AI Lab, Cambridge, MA, USA Abstract ... $Id: mri_robust_register.cpp,v 1.77 2016/01/20 23:36:17 greve Exp $ --mov: Using imageDump.mgz as movable/source volume. --dst: Using /home/valentil/kg98/Valentina/ENIGMA_sex_differences/Scripts/rerun_reconall/NIAAA_MRI00032_1-out/tmp/hippoSF_T1_v21_left//hippoAmygBinaryMask_autoCropped.mgz as target volume. --lta: Output transform as trash.lta . --mapmovhdr: Will save header adjusted movable as imageDump_coregistered.mgz ! --affine: Enableing affine transform! --sat: Using saturation 50 in M-estimator! reading source 'imageDump.mgz'... reading target '/home/valentil/kg98/Valentina/ENIGMA_sex_differences/Scripts/rerun_reconall/NIAAA_MRI00032_1-out/tmp/hippoSF_T1_v21_left//hippoAmygBinaryMask_autoCropped.mgz'... Registration::setSourceAndTarget(MRI s, MRI t, keeptype = TRUE ) Type Source : 0 Type Target : 3 ensure both FLOAT (3) Reordering axes in mov to better fit dst... ( -1 3 -2 ) MRIreorder() ----------- xdim=-1 ydim=3 zdim=-2 src 131 241 99, 0.250000 0.250000 0.250000 dst 131 99 241, 0.250000 0.250000 0.250000 Determinant after swap : 0.015625 Mov: (0.25, 0.25, 0.25)mm and dim (131, 99, 241) Dst: (1, 1, 1)mm and dim (36, 31, 52) Asserting both images: 1mm isotropic - reslicing Mov ... -- changing data type from 0 to 3 (noscale = 0)... -- Original : (0.25, 0.25, 0.25)mm and (131, 99, 241) voxels. -- Resampled: (1, 1, 1)mm and (36, 31, 61) voxels. -- Reslicing using cubic bspline MRItoBSpline degree 3 - reslicing Dst ... -- Original : (1, 1, 1)mm and (36, 31, 52) voxels. -- Resampled: (1, 1, 1)mm and (36, 31, 61) voxels. -- Reslicing using cubic bspline MRItoBSpline degree 3 Registration::computeMultiresRegistration - computing centroids - computing initial transform -- using translation info - Get Gaussian Pyramid Limits ( min size: 16 max size: -1 ) - Build Gaussian Pyramid ( Limits min steps: 0 max steps: 0 ) - Build Gaussian Pyramid ( Limits min steps: 0 max steps: 0 ) - initial transform: Ti = [ ... 1.0000000000000 0 0 0.8149148880414 0 1.0000000000000 0 0.6090566088288 0 0 1.0000000000000 -10.4098892270167 0 0 0 1.0000000000000 ] - initial iscale: Ii =1 Resolution: 0 S( 36 31 61 ) T( 36 31 61 ) Iteration(f): 1 -- diff. to prev. transform: 41.9486 Iteration(f): 2 -- diff. to prev. transform: 11.6817 Iteration(f): 3 -- diff. to prev. transform: 6.67838 Iteration(f): 4 -- diff. to prev. transform: 7.29474 Iteration(f): 5 -- diff. to prev. transform: 6.81314 max it: 5 reached! - final transform: Tf = [ ... 0.7187784317010 -0.7039049689666 0.2089453505114 9.6512053486057 0.0746036190898 0.1293914974162 0.0117417560783 14.1524707759350 -0.2465560963278 -0.6917283813412 1.0760848577522 2.1884563417970 0 0 0 1.0000000000000 ] - final iscale: If = 1 ********************************************************** * * WARNING: Registration did not converge in 5 steps! * Problem might be ill posed. * Please inspect output manually! * ********************************************************** Final Transform: Adjusting final transform due to initial resampling (voxel or size changes) ... M = [ ... -0.1796946067577 0.0522363297057 0.1759762437984 14.8124989197474 -0.0186509047068 0.0029354397057 -0.0323478742666 20.2771618184274 0.0616390300951 0.2690212102028 0.1729321033528 -29.4307376768930 0 0 0 1.0000000000000 ] Determinant : -0.00250544 Decompose into Rot * Shear * Scale : Rot = [ ... -0.8007108168802 0.0158873366257 0.5988403629232 -0.5106339926748 0.5045973184270 -0.6961569304125 0.3132333207947 0.8632086299216 0.3959239168991 ] Shear = [ ... 1.0000000000000 0.1745654451980 -0.3575927026144 0.2370458412539 1.0000000000000 0.6912992608189 -0.4065670242970 0.5788093472423 1.0000000000000 ] Scale = diag([ 0.1727145993906 0.2345325414370 0.1963688302739 ]) writing output transformation to trash.lta ... converting VOX to RAS and saving RAS2RAS... mapmovhdr: Changing vox2ras MOV header (to map to DST) ... To check aligned result, run: freeview -v /home/valentil/kg98/Valentina/ENIGMA_sex_differences/Scripts/rerun_reconall/NIAAA_MRI00032_1-out/tmp/hippoSF_T1_v21_left//hippoAmygBinaryMask_autoCropped.mgz imageDump_coregistered.mgz Registration took 0 minutes and 1 seconds. Thank you for using RobustRegister! If you find it useful and use it for a publication, please cite: Highly Accurate Inverse Consistent Registration: A Robust Approach M. Reuter, H.D. Rosas, B. Fischl. NeuroImage 53(4):1181-1196, 2010. http://dx.doi.org/10.1016/j.neuroimage.2010.07.020 http://reuter.mit.edu/papers/reuter-robreg10.pdf Highly Accurate Inverse Consistent Registration: A Robust ...<http://reuter.mit.edu/papers/reuter-robreg10.pdf> reuter.mit.edu Highly Accurate Inverse Consistent Registration: A Robust Approach Martin Reutera,b,c,∗, H. Diana Rosasa,b, Bruce Fischla,b,c aMassachusetts General Hospital / Harvard Medical School, Boston, MA, USA bMartinos Center for Biomedical Imaging, 143 13th Street, Charlestown, MA, USA cMIT Computer Science and AI Lab, Cambridge, MA, USA Abstract ... Reading contexts of file /usr/local/freesurfer/devel-20180612/average/HippoSF/atlas/compressionLookupTable.txt Constructing image-to-world transform from header information (asegModCHA.mgz) Constructing image-to-world transform from header information (/home/valentil/kg98/Valentina/ENIGMA_sex_differences/Scripts/rerun_reconall/NIAAA_MRI00032_1-out/tmp/hippoSF_T1_v21_left/imageDump.mgz) Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Transforming points Reading contexts of file /usr/local/freesurfer/devel-20180612/average/HippoSF/atlas/compressionLookupTable.txt -------------- Making Left-Cerebral-Cortex map to reduced label 1 Making alveus map to reduced label 1 Making subiculum-body map to reduced label 1 Making subiculum-head map to reduced label 1 Making Hippocampal_tail map to reduced label 1 Making molecular_layer_HP-body map to reduced label 1 Making molecular_layer_HP-head map to reduced label 1 Making GC-ML-DG-body map to reduced label 1 Making GC-ML-DG-head map to reduced label 1 Making CA4-body map to reduced label 1 Making CA4-head map to reduced label 1 Making CA1-body map to reduced label 1 Making CA1-head map to reduced label 1 Making CA3-body map to reduced label 1 Making CA3-head map to reduced label 1 Making HATA map to reduced label 1 Making fimbria map to reduced label 1 Making presubiculum-body map to reduced label 1 Making presubiculum-head map to reduced label 1 Making parasubiculum map to reduced label 1 Making Lateral-nucleus map to reduced label 1 Making Paralaminar-nucleus map to reduced label 1 Making Basal-nucleus map to reduced label 1 Making Accessory-Basal-nucleus map to reduced label 1 Making Corticoamygdaloid-transitio map to reduced label 1 Making Central-nucleus map to reduced label 1 Making Cortical-nucleus map to reduced label 1 Making Medial-nucleus map to reduced label 1 Making Anterior-amygdaloid-area-AAA map to reduced label 1 -------------- Making Left-Cerebral-White-Matter map to reduced label 2 -------------- Making Left-Lateral-Ventricle map to reduced label 3 -------------- Making Left-choroid-plexus map to reduced label 4 -------------- Making hippocampal-fissure map to reduced label 5 Making Unknown map to reduced label 5 -------------- Making Left-VentralDC map to reduced label 6 -------------- Making Left-Putamen map to reduced label 7 -------------- Making Left-Pallidum map to reduced label 8 -------------- Making Left-Thalamus-Proper map to reduced label 9 -------------- Making Left-Accumbens-area map to reduced label 10 -------------- Making Left-Caudate map to reduced label 11 Smoothing mesh collection with kernel size 3.000000 ...Smoothing class: 0 Smoothing class: 1 Smoothing class: 2 Smoothing class: 3 Smoothing class: 4 Smoothing class: 5 Smoothing class: 6 Smoothing class: 7 Smoothing class: 8 Smoothing class: 9 Smoothing class: 10 done numberOfContrasts = 1 Image: 4 AtlasMeshToIntensityImage 1 0.398942 SLIDING L-BFGS Verbose: 0 MaximalDeformationStopCriterion: 1e-10 LineSearchMaximalDeformationIntervalStopCriterion: 1e-10 MaximumNumberOfIterations: 1000 BFGS-MaximumMemoryLength: 12 Resolution 1, iteration 1 Did one deformation step of max. 0 voxels in 4.2457 seconds minLogLikelihoodTimesPrior = 2.4171e-16 Smoothing mesh collection with kernel size 2.000000 ...Smoothing class: 0 Smoothing class: 1 Smoothing class: 2 Smoothing class: 3 Smoothing class: 4 Smoothing class: 5 Smoothing class: 6 Smoothing class: 7 Smoothing class: 8 Smoothing class: 9 Smoothing class: 10 done numberOfContrasts = 1 Image: 4 AtlasMeshToIntensityImage 1 0.398942 SLIDING L-BFGS Verbose: 0 MaximalDeformationStopCriterion: 1e-10 LineSearchMaximalDeformationIntervalStopCriterion: 1e-10 MaximumNumberOfIterations: 1000 BFGS-MaximumMemoryLength: 12 Resolution 2, iteration 1 Did one deformation step of max. 0 voxels in 4.2086 seconds minLogLikelihoodTimesPrior = 2.4171e-16 Fitting mesh to synthetic image from ASEG took 21.6669 seconds Error using segmentSubjectT1_autoEstimateAlveusML (line 533) Fitting mesh to synthetic image resulted in no deformation @#@FSTIME 2018:09:25:10:18:56 run_segmentSubjectT1_autoEstimateAlveusML.sh N 13 e 65.81 S 1.73 U 34.89 P 55% M 1012356 F 271 R 219755 W 0 c 3155 w 11451 I 146312 O 477024 L 11.09 10.51 11.31 @#@FSLOADPOST 2018:09:25:10:20:02 run_segmentSubjectT1_autoEstimateAlveusML.sh N 13 35.39 16.50 13.25 Linux m3a007 3.10.0-514.26.1.el7.x86_64 #1 SMP Thu Jun 29 16:05:25 UTC 2017 x86_64 x86_64 x86_64 GNU/Linux T1 hippocampal subfields exited with ERRORS at Tue Sep 25 10:20:02 AEST 2018 -- Sally Grace, BSc (Hons), PhD Candidate. Centre for Mental Health, Faculty of Health, Arts and Design Swinburne University of Technology

2 2

denoising and recon-all
by Lisa Crystal Krishnamurthy 11 Oct '18

11 Oct '18

External Email - Use Caution Hi, I have found that denoising the MPRAGE prior to FS recon-all seems to improve the segmentation (see attached pdf). However, it is not clear if the denoising algorithm may cause some signal intensity changes (especially in voxels with partial voluming at the edge of the brain) that violate assumptions of recon-all. Could you help me understand what the assumptions of your algorithm are, and what I need to do to make sure my images conform to those assumptions? Your help is greatly appreciated. Best, -Lisa

1 0

Recon-all with/without T2 give the same stats
by Jianzhong Chen 11 Oct '18

11 Oct '18

External Email - Use Caution Dear All, I tried recon-all with/without T2 image using FS 5.3.0. Then I found that the aseg stats and aparc stats with/without T2 are exactly the same. When I check the pial surfaces, I can see that they are very similar but not the same. So is it reasonable to see exactly the same stats? An example of the pial surfaces with/without T2 is attached here. The red line is lh.pial with T2 and the yellow line is the lh.pial without T2 My command for recon-all with T2 is: recon-all -s ${subject} -i T1.nii -all -T2 T2.nii -T2pial -sd ${T2_dir} My command for recon-all without T2 is: recon-all -s ${subject} -i T1.nii -all -sd ${T1_dir} Thanks, Jianzhong

2 6

Difference between 2 mgh files in surf dir (or: effects of 'mri_surf2surf --prune' without smoothing)
by ts+ml＠rcmd.org 11 Oct '18

11 Oct '18

External Email - Use Caution Dear FreeSurfer experts, what is the difference between the following 2 files found in a subject's surf/ directory: * file #1: lh.area.fsaverage.mgh * file #2: lh.area.fwhm0.fsaverage.mgh I see in the logs (<SUBJECT>/scripts/recon-all.log), that file #2 is generated from file #1 by running the following command: mri_surf2surf --prune --s fsaverage --hemi lh --fwhm 0 --sval lh.area.fsaverage.mgh --tval lh.area.fwhm0.fsaverage.mgh --cortex I know that if fwhm > 0, spatial smoothing will be applied. But what happens if fwhm is 0? From a quick comparison of the data in both files, they seem to be identical, but the hashes of the files differ. Is the difference only about the header? Thanks for your help, -- Dr. Tim Schäfer Postdoc Computational Neuroimaging Department of Child and Adolescent Psychiatry, Psychosomatics and Psychotherapy University Hospital Frankfurt, Goethe University Frankfurt am Main, Germany

2 2

masking cortical thickness
by N Saf 11 Oct '18

11 Oct '18

External Email - Use Caution Dear Freesurfer Experts, I wonder how can I show the cortical thickness surface on just one region of my interest in Freeview. assume that I want to see the thickness surface of one of the regions in DKTatlas. does it work if I make a binary mask of that specific region on aparcDKTatlas.annot and then multiply it to the ?h.thickness surface ?? any help will be appreciated. Best Regards, Nazanin

3 6

Pet Surfer
by David.Mackarthy 10 Oct '18

10 Oct '18

External Email - Use Caution Dear Dr Greeve, I am following petsurfer pipeline to do surface based analysis of PET-FDG data and I have interest in doing partial volume correction. I would like to inquire about: 1) Can partial volume correction be applied on white mater, or just on gray matter? 2) Can partial volume correction be applied between gray matter and CSF? According to Muller-Gartner (MG) Method we segment the brain to gray/white and csf, then we apply partial volume correction only on gray matter. why not on the other parts. I apologies if my question is basics and I appreciate your input. Thanks for publishing this valuable manuscript https://www.ncbi.nlm.nih.gov/pubmed/26915497 With all respect! Dev

2 3

Moving VOI to FreeSurfer Space
by Sarosh, Cyrus 10 Oct '18

10 Oct '18

External Email - Use Caution Hello FreeSurfer Experts, I'm trying to find the best method and command to move a hand drawn VOI in original MRI space to FreeSurfer space. These scans have already been ran through the entire FreeSurfer process and the hand drawn VOIs have already been converted to .mgz files. Thanks for your help Thanks, Cyrus Cyrus Sarosh [Description: Description: Description: Description: Description: Description: Description: cid:image002.png@01CC3BB7.F8D9A990] Research Area Specialist Associate Department of Radiology Functional Neuroimaging, Cognitive and Mobility Laboratory Domino's Farms Lobby B Suite 1200 24 Frank Lloyd Wright Dr PO BOX 362 Ann Arbor, MI 48106 t: 734-998-4790 f: 734-998-8403 ********************************************************** Electronic Mail is not secure, may not be read every day, and should not be used for urgent or sensitive issues

2 1

Re: [Freesurfer] error glm-fit-sim (message 11)
by Caroline Beelen 10 Oct '18

10 Oct '18

External Email - Use Caution Few clusters survive glm fit. No clusters seem to survive glm-fit-sim (when visualizing). However, I'm not sure if I should also see this in a text file. Opening cache.th40(or 30).abs.sig.cluster.summary it seems to have found no clusters. At least, under the line ClusterNo Max VtxMax Size MNIX MNIY etc. nothing is written. But does a lack of significant clusters also explain the weird messages of "colortable with 1 entries read (originally none)" (& "CTABisEntryValid: index -1 was OOB" & "Resource temporarily unavailable") that show up in the terminal? Thanks very much, Caroline -----Oorspronkelijk bericht----- Van: freesurfer-bounces(a)nmr.mgh.harvard.edu [mailto:freesurfer-bounces@nmr.mgh.harvard.edu] Namens freesurfer-request(a)nmr.mgh.harvard.edu Verzonden: dinsdag 9 oktober 2018 18:00 Aan: freesurfer(a)nmr.mgh.harvard.edu Onderwerp: Freesurfer Digest, Vol 176, Issue 9 Send Freesurfer mailing list submissions to freesurfer(a)nmr.mgh.harvard.edu To subscribe or unsubscribe via the World Wide Web, visit https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer or, via email, send a message with subject or body 'help' to freesurfer-request(a)nmr.mgh.harvard.edu You can reach the person managing the list at freesurfer-owner(a)nmr.mgh.harvard.edu When replying, please edit your Subject line so it is more specific than "Re: Contents of Freesurfer digest..." Today's Topics: 1. FW: RE: error glm-fit-sim (Caroline Beelen) 2. Re: output values of vertices on surface into a textfile (Sims, Sara A) 3. Re: output values of vertices on surface into a textfile (Bruce Fischl) 4. Re: output values of vertices on surface into a textfile (Sims, Sara A) 5. Re: Hippocampal segmentation with an additional scan ERROR kvlGEMSMatlab (Pradeep) 6. Re: Hippocampal segmentation with an additional scan ERROR kvlGEMSMatlab (Iglesias Gonzalez, Eugenio) 7. hippocampal subfields (Marcel Heers) 8. Re: hippocampal subfields (Iglesias Gonzalez, Eugenio) 9. Re: Difference between 6.0.1 and 6.0.0 (Dicamillo, Robert) 10. Talairach failed in monkey brain (Frehiwot Woldeyes) 11. Re: FW: RE: error glm-fit-sim (Greve, Douglas N.,Ph.D.) 12. Mri Convert Help (Justin Arnett) 13. Re: Mri Convert Help (Greve, Douglas N.,Ph.D.) 14. Re: Recon-all with/without T2 give the same stats (Greve, Douglas N.,Ph.D.) 15. Re: Recon-all with/without T2 give the same stats (Jianzhong Chen) 16. m.khalilian(a)sutech.ac.ir sent you files via WeTransfer (WeTransfer) 17. very bad result using "mris_ca_train" and "mris_ca_label" (Maedeh Khalilian) 18. Re: Talairach failed in monkey brain (Greve, Douglas N.,Ph.D.) 19. Re: Recon-all with/without T2 give the same stats (Greve, Douglas N.,Ph.D.) 20. Re: very bad result using "mris_ca_train" and "mris_ca_label" (Bruce Fischl) 21. Re: Extract individual vertex thickness values from cortical parcellations (Greve, Douglas N.,Ph.D.) 22. Re: nu_correct disk i/o issues (Greve, Douglas N.,Ph.D.) 23. Re: masking cortical thickness (Greve, Douglas N.,Ph.D.) ---------------------------------------------------------------------- Message: 1 Date: Mon, 8 Oct 2018 16:12:06 +0000 From: Caroline Beelen <caroline.beelen(a)kuleuven.be> Subject: [Freesurfer] FW: RE: error glm-fit-sim To: "freesurfer-bounces(a)nmr.mgh.harvard.edu" <freesurfer-bounces(a)nmr.mgh.harvard.edu>, "freesurfer(a)nmr.mgh.harvard.edu" <freesurfer(a)nmr.mgh.harvard.edu> Message-ID: <4051d561c94b40e0b641202b7bad7192(a)ICTS-S-EXMBX25.luna.kuleuven.be> Content-Type: text/plain; charset="us-ascii" External Email - Use Caution The command: gop@gop-linux:~/freesurfer$ freeview -f $SUBJECTS_DIR/avgsubject/surf/lh.inflated:overlay=lh.area.glmdir/groupeffect/cache.th40.abs.sig.cluster.mgh:overlay_threshold=2,5:annot=lh.area.glmdir/groupeffect/cache.th40.abs.sig.ocn.annot -viewport 3d The output: reading colortable from annotation file... colortable with 1 entries read (originally none) colortable with 1 entries read (originally none) CTABisEntryValid: index -1 was OOB Resource temporarily unavailable Resource temporarily unavailable B ginally none) CTABisEntryValid: index -1 was OOB Resource temporarily unavailable Resource temporarily unavailable B ginally none) CTABisEntryValid: index -1 was OOB Resource temporarily unavailable Resource temporarily unavailable B ginally none) The command with a bit less strict significance level: gop@gop-linux:~/freesurfer$ freeview -f $SUBJECTS_DIR/avgsubject/surf/lh.inflated:overlay=lh.area.glmdir/groupeffect/cache.th40.abs.sig.cluster.mgh:overlay_threshold=3,5:annot=lh.area.glmdir/groupeffect/cache.th40.abs.sig.ocn.annot -viewport 3d The output: reading colortable from annotation file... colortable with 1 entries read (originally none) colortable with 1 entries read (originally none) (and then it says no more). Both commands visualize the data, but no clusters are found. Is this how it should be? Does it mean that there is a problem with reading colortables or not? Best, Caroline Message: 2 Date: Mon, 1 Oct 2018 17:10:28 +0000 From: "Greve, Douglas N.,Ph.D." <DGREVE(a)mgh.harvard.edu<mailto:DGREVE@mgh.harvard.edu>> Subject: Re: [Freesurfer] error glm-fit-sim To: "freesurfer(a)nmr.mgh.harvard.edu<mailto:freesurfer@nmr.mgh.harvard.edu>" <freesurfer(a)nmr.mgh.harvard.edu<mailto:freesurfer@nmr.mgh.harvard.edu>> Message-ID: <504d6069-399f-163d-39a9-dcc6c02d32f9(a)mgh.harvard.edu<mailto:504d6069-399f-163d-39a9-dcc6c02d32f9@mgh.harvard.edu>> Content-Type: text/plain; charset="Windows-1252" What says that? What command are you running? What is the other terminal output? On 10/01/2018 04:54 AM, Caroline Beelen wrote: > > ????????External Email - Use Caution > > Dear FS, > > I followed the group analysis page for glm analysis and it worked fine > until the final step. I?m using version 5.3. After performing > glm-fit-sim it says ?glm-fit done? (for all files), but when trying to > open these it says: > > ?Reading colortable from annotation file? > > Colortable with 1 entries read (originally none) > > Colortable with 1 entries read (originally none) > > CTABisentryvalid: index -1 was 00B > > Resource temporarily unavailable > > CTABisentryvalid: index -1 was 00B > > Resource temporarily unavailable > > Resource temporarily unavailable? > > (etc..) > > It opens, but nothing is visible and I have the feeling it cannot open > color files?? > > MRI glmfit gives positive clusters for some values. I therefore used > the abs sign for glm-fit-sim. > > What might be wrong? > > Thanks. Caroline

2 1

Re: [Freesurfer] Disjointed thalamic nuclei segmentation
by Iglesias Gonzalez, Eugenio 10 Oct '18

10 Oct '18

External Email - Use Caution Thanks, Niels. It is definitely a matter of hyperintensity of voxels close to the lateral boundary, particularly in the anterior region (which has a fuzzy boundary with the white matter in first place). Whether these hyperintensities are due to MS pathology, that I do not know (I’m really not the person to tell…). These could be cleaned up with postprocessing (e.g., filling holes, or using a Markov Random Field), but the impact on the volumes would be negligible. Cheers, /Eugenio -- Juan Eugenio Iglesias ERC Senior Research Fellow Centre for Medical Image Computing (CMIC) University College London http://www.jeiglesias.com http://cmictig.cs.ucl.ac.uk/ From: Niels Bergsland <theniels(a)gmail.com> Date: Wednesday, 10 October 2018 at 15:42 To: "Iglesias Gonzalez, Eugenio" <e.iglesias(a)ucl.ac.uk> Subject: Re: [Freesurfer] Disjointed thalamic nuclei segmentation Hi Eugenio, Thanks for the lightning fast response! :) I've attached them here. It looks like it might actually be the opposite in this case in that the voxels here are darker. Some of the cases are MS patients and so it's conceivable that they have thalamic lesions. Admittedly though, I have not really spent time in the past looking for thalamic lesions on T1-weighted images and this hasn't been an issue in the past since the thalamic segmentation in the aseg doesn't end up being affected by these voxels with lower intensities. If you would like to look at some other cases or need other outputs, just let me know and I'm happy to send them along! Thanks again, Niels On Wed, Oct 10, 2018 at 4:30 PM Iglesias Gonzalez, Eugenio <e.iglesias(a)ucl.ac.uk<mailto:e.iglesias@ucl.ac.uk>> wrote: External Email - Use Caution Thanks, Niels. I suspect the problem is that, for some reason, such voxels have much brighter intensity than their surroundings, and are classified as white matter / reticular nucleus. But that shouldn’t be happening that far from the lateral boundary of the thalamus… Can you please send us norm.mgz and the thalamic segmentation file, so I can take a look? Cheers, /Eugenio -- Juan Eugenio Iglesias ERC Senior Research Fellow Centre for Medical Image Computing (CMIC) University College London http://www.jeiglesias.com http://cmictig.cs.ucl.ac.uk/ From: <freesurfer-bounces(a)nmr.mgh.harvard.edu<mailto:freesurfer-bounces@nmr.mgh.harvard.edu>> on behalf of Niels Bergsland <theniels(a)gmail.com<mailto:theniels@gmail.com>> Reply-To: Freesurfer support list <freesurfer(a)nmr.mgh.harvard.edu<mailto:freesurfer@nmr.mgh.harvard.edu>> Date: Wednesday, 10 October 2018 at 15:22 To: Freesurfer Mailing List <freesurfer(a)nmr.mgh.harvard.edu<mailto:freesurfer@nmr.mgh.harvard.edu>> Subject: [Freesurfer] Disjointed thalamic nuclei segmentation External Email - Use Caution Hi Eugenio, Thank you again for the fantastic thalamic nuclei segmentation tool! I'm going through and QC'ing a batch of subjects that were processed through the pipeline. Data was originally processed with freesurfer-6.0.0 and then processed through the thalamic nuclei stream (dev-20180818-e30e6f9). I have found that for some cases, there are some isolated patches of voxels that are disconnected and seemingly in the wrong place. I've noticed it primarily for the MDm and my impression is that it tends to happen more in the right hemisphere. I've attached a sample image to show you what I mean. The red cross has been placed on a set of these voxels. It is not the case that they are connected in 3D. Any input is appreciated and thanks again! If it can be useful, I'm happy to upload the data -Niels _______________________________________________ Freesurfer mailing list Freesurfer(a)nmr.mgh.harvard.edu<mailto:Freesurfer@nmr.mgh.harvard.edu> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

1 0

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