Freesurfer

freesurfer@nmr.mgh.harvard.edu

19 participants
30088 discussions

affine registration with mri_robust_register
by Nathaniel Reynolds 22 Nov '13

22 Nov '13

Hi all, I can't find any explanation of how to perform an affine registration with mdi_robust_register. Could someone please explain? Thanks, Nathaniel

1 0

Re: [Freesurfer] scale bar question
by Douglas N Greve 22 Nov '13

22 Nov '13

Just click on a point and look at the vertex RAS. At zoom=1, you have to translate the surface over to the scale bar. At zoom=1 I got about 10mm, but it did seem to change when I zoomed in. doug On 11/22/2013 12:23 PM, Cesar Echavarria wrote: > Actually, Doug can you tell me how to get the xyz coordinates on > Tksurfer so I can figure out the length of the bar as you did? > > Thanks > > Cesar > > > On Fri, Nov 22, 2013 at 12:19 PM, Cesar Echavarria > <cechavarria60(a)gmail.com <mailto:cechavarria60@gmail.com>> wrote: > > I meant zoom in on the brain. I guess I'll take Doug's approach > and estimate the scale of the bar with the xyz coordinates. > > Cesar > > > On Fri, Nov 22, 2013 at 12:16 PM, Bruce Fischl > <fischl(a)nmr.mgh.harvard.edu <mailto:fischl@nmr.mgh.harvard.edu>> > wrote: > > Hi Cesar > > what do you mean by "scale the brain up"? > > Bruce > > On Fri, 22 Nov 2013, Cesar Echavarria wrote: > > Hi freesurfers, > Anyone know what the scale bar in TkSurfer represents? > >From past threads I > gathered that it represent 1 cm on the human brain. > However, the size of the > bar displayed gets smaller as I scale the brain up which > does not make sense > if the bar is supposed to represent a constant length. > > Cesar > > > > > The information in this e-mail is intended only for the person > to whom it is > addressed. If you believe this e-mail was sent to you in error > and the e-mail > contains patient information, please contact the Partners > Compliance HelpLine at > http://www.partners.org/complianceline . If the e-mail was > sent to you in error > but does not contain patient information, please contact the > sender and properly > dispose of the e-mail. > > > > > -- > Cesar Echavarria > Research Assistant > Martinos Center for Biomedical Imaging, Harvard-MGH > Brain & Cognitive Sciences | MIT 2012 > > > > > -- > Cesar Echavarria > Research Assistant > Martinos Center for Biomedical Imaging, Harvard-MGH > Brain & Cognitive Sciences | MIT 2012 -- Douglas N. Greve, Ph.D. MGH-NMR Center greve(a)nmr.mgh.harvard.edu Phone Number: 617-724-2358 Fax: 617-726-7422 Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting FileDrop: https://gate.nmr.mgh.harvard.edu/filedrop2 www.nmr.mgh.harvard.edu/facility/filedrop/index.html Outgoing: ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/

1 0

Re: [Freesurfer] scale bar question
by Cesar Echavarria 22 Nov '13

22 Nov '13

I meant zoom in on the brain. I guess I'll take Doug's approach and estimate the scale of the bar with the xyz coordinates. Cesar On Fri, Nov 22, 2013 at 12:16 PM, Bruce Fischl <fischl(a)nmr.mgh.harvard.edu>wrote: > Hi Cesar > > what do you mean by "scale the brain up"? > > Bruce > > On Fri, 22 Nov 2013, Cesar Echavarria wrote: > > Hi freesurfers, >> Anyone know what the scale bar in TkSurfer represents? >From past threads >> I >> gathered that it represent 1 cm on the human brain. However, the size of >> the >> bar displayed gets smaller as I scale the brain up which does not make >> sense >> if the bar is supposed to represent a constant length. >> >> Cesar >> >> >> > > The information in this e-mail is intended only for the person to whom it > is > addressed. If you believe this e-mail was sent to you in error and the > e-mail > contains patient information, please contact the Partners Compliance > HelpLine at > http://www.partners.org/complianceline . If the e-mail was sent to you in > error > but does not contain patient information, please contact the sender and > properly > dispose of the e-mail. > > -- Cesar Echavarria Research Assistant Martinos Center for Biomedical Imaging, Harvard-MGH Brain & Cognitive Sciences | MIT 2012

1 0

scale bar question
by Cesar Echavarria 22 Nov '13

22 Nov '13

Hi freesurfers, Anyone know what the scale bar in TkSurfer represents? From past threads I gathered that it represent 1 cm on the human brain. However, the size of the bar displayed gets smaller as I scale the brain up which does not make sense if the bar is supposed to represent a constant length. Cesar

1 0

Order of edits for longitudinal data
by Rongxiang Tang 22 Nov '13

22 Nov '13

Dear All, I was would like to double check the order of edits for longitudinal data. I did edits for my cross-sectional data, and for my base data, I did the following: 1. skull strips edits (if necessary) 2. control points edits 3. white matter edits 4. pial edits 5.brain final surfs edits. I was also wondering if correction for topological defects is needed for the base data. Thanks, Rongxiang

2 1

Re: [Freesurfer] help with save_nifti
by Leila Reddy 22 Nov '13

22 Nov '13

Hi Doug, Thanks for your reply. I have the problem even if I use MRIread and MRIwrite. But I have done a little more digging following your suggestions. Basically, there are differences in the output of mri_vol2surf (that I should have spotted earlier), as well as differences in the header information. I highlight the differences below, and I have put the detailed outputs of mri_vol2surf at the end in case you still need to look at it. I think the clearest way of explaining the different outcomes is as follows.: 1) mri_vol2surf performed on a nifti volume generated by Nibabel. Nothing done in Matlab. mri_vol2surf ouput: "Number of source voxels hit = 22587". "niiRead(): NIFTI_UNITS_UNKNOWN, assuming mm" Results looks fine in tksurfer (i.e., I have a distribution of voxel values in the Configure Overlay window). If I look at the nifti header, hdr.volres= [2.0000e+10 2.0000e+10 2.0000e+10] and several other fields (e.g., vox2ras etc) also have the same magnitude. In addition, hdr.niftihdr.dim'=[3 76 76 93 1 1 1 1] 2) The Nibabel nifti is opened in matlab with MRIread and immediately saved with MRIwrite. Nothing else is done to the data in Matlab. Output of MRIread=WARNING: xyz units code 0 is unrecognized Rescaling NIFTI: slope = 1, intercept = 0 mri_vol2surf output: "Number of source voxels hit = 1". Nothing about NIFTI_UNITS (note the header is the same as in 1 above). Rest of mri_vol2surf output is exactly as for #1 above. The header values are as above except hdr.niftihdr.dim'=[4 76 76 93 1 1 1 1] 3) The Nibabel nifti is opened in matlab with MRIread. The 3d volume is saved with MRIwrite but with the header information from a FSFAST created nifti. mri_vol2surf says "Number of source voxels hit = 22587". hdr.volres=[2 2 2] and other fields have values of the same order. hdr.niftihdr.dim'=[4 76 76 93 1 1 1 1] I don't know if I have focused on the relevant header information here so I also attach the headers from steps 1-3 above. Thanks a lot for your help and time on this. I will often have to go back and forth between python functions and Matlab/FS-FAST scripts so this will be very useful. Leila Detailed output of mri_vol2surf with fsfast header (#3 above): mri_vol2surf --src fs_hdr.nii --hemi rh --out_type mgh --float2int round --o fshdr-rh.mgh --projfrac 0.5 --trgsubject fsaverage --regheader fsaverage srcvol = fs_hdr.nii srcreg unspecified srcregold = 0 srcwarp unspecified surf = white hemi = rh trgsubject = fsaverage surfreg = sphere.reg ProjFrac = 0.5 thickness = thickness reshape = 0 interp = nearest float2int = round GetProjMax = 0 INFO: float2int code = 0 Done loading volume Computing registration from header. Using /APPLICATIONS/freesurfer/subjects/fsaverage/mri/orig.mgz as target reference. Reading surface /APPLICATIONS/freesurfer/subjects/fsaverage/surf/rh.white Done reading source surface Reading thickness /APPLICATIONS/freesurfer/subjects/fsaverage/surf/rh.thickness Done Mapping Source Volume onto Source Subject Surface 1 0.5 0.5 0.5 using old Done mapping volume to surface Number of source voxels hit = 22587 Writing to fshdr-rh.mgh Dim: 163842 1 1 Detailed output of mri_vol2surf with nibabel header (after opening and saving in matlab) (#2 above): mri_vol2surf --src test.nii --hemi rh --out_type mgh --float2int round --o test-rh.mgh --projfrac 0.5 --trgsubject fsaverage --regheader fsaverage srcvol = test.nii srcreg unspecified srcregold = 0 srcwarp unspecified surf = white hemi = rh trgsubject = fsaverage surfreg = sphere.reg ProjFrac = 0.5 thickness = thickness reshape = 0 interp = nearest float2int = round GetProjMax = 0 INFO: float2int code = 0 Done loading volume Computing registration from header. Using /APPLICATIONS/freesurfer/subjects/fsaverage/mri/orig.mgz as target reference. Reading surface /APPLICATIONS/freesurfer/subjects/fsaverage/surf/rh.white Done reading source surface Reading thickness /APPLICATIONS/freesurfer/subjects/fsaverage/surf/rh.thickness Done Mapping Source Volume onto Source Subject Surface 1 0.5 0.5 0.5 using old Done mapping volume to surface Number of source voxels hit = 1 Writing to test-rh.mgh Dim: 163842 1 1

1 0

comparing left and right cortical thickness
by Jacobs H (NP) 22 Nov '13

22 Nov '13

Dear FreeSurfers, Is it possible to compare left and right cortical thickness patterns, vertex-wise? Thanks Heidi ================================================ Dr. Heidi Jacobs Postdoc researcher Faculty of Health, Medicine and Life Sciences School for Mental Health and Neurosciences Division Cognitive Neuropsychiatry and Clinical Neurosciences Alzheimer Center Limburg h.jacobs(a)maastrichtuniversity.nl www.maastrichtuniversity.nl www.heidijacobs.nl Dr. Tanslaan 12, 6229 ET Maastricht P.O. Box 616, 6200 MD Maastricht, The Netherlands T +31 43 38 84 090 F +31 43 38 84 092 ================================================

1 0

problem feeding design matrix to mri_glmfit
by Laura M. Tully 21 Nov '13

21 Nov '13

Hello freesurfers, I'm trying to run mri_glmfit between groups analysis (patients vs. controls) with gender, age, medication (continuous), and duration of illness as covariates. Because the control group have values of zeros for medication and duration of illness (which results in an illconditioned matrix), I have created my own design matrix so that controls do not have a column for medication or duration of illness and am feeding this matrix into mri_glmfit instead of the fsgd file. However, this seems to cause a problem with the --y input; when I try to run mri_glmfit, i get an error message "error in uncompressing" - I'm not quite sure what this means, perhaps I did something wrong in creating the design matrix? See the full output below. I've also attached the design matrix, contrast file, and the original fsgd (which results in the ill conditioned matrix). Best, Laura. *Command line: * mri_glmfit --y lh_grpDiffs_thickness10.mgh --X new_mat_files/X_lh_grp_thick_med_dur.mat --C Contrast1B_grp_Diffs_HCvsSZ_Med_Dur.mtx --surf fsaverage lh --cortex --glmdir lh.grp_diffs_thickness_med_dur_112013.glmdir --illcond *Output:* Reading source surface /ncf/snp/04/SCORE/freesurfer_analysis/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.6 2011/05/05 20:54:25 greve Exp $ cwd /ncf/snp/04/SCORE/freesurfer_analysis/glm_analysis cmdline mri_glmfit --y lh_grpDiffs_thickness10.mgh --X new_mat_files/X_lh_grp_thick_med_dur.mat --C Contrast1B_grp_Diffs_HCvsSZ_Med_Dur.mtx --surf fsaverage lh --cortex --glmdir lh.grp_diffs_thickness_med_dur_112013.glmdir --illcond sysname Linux hostname ncfws13.rc.fas.harvard.edu machine x86_64 user ltully FixVertexAreaFlag = 1 UseMaskWithSmoothing 1 OneSampleGroupMean 0 y /ncf/snp/04/SCORE/freesurfer_analysis/glm_analysis/lh_grpDiffs_thickness10.mgh logyflag 0 X new_mat_files/X_lh_grp_thick_med_dur.mat usedti 0 labelmask /ncf/snp/04/SCORE/freesurfer_analysis/fsaverage/label/lh.cortex.label maskinv 0 glmdir lh.grp_diffs_thickness_med_dur_112013.glmdir IllCondOK 1 ReScaleX 1 DoFFx 0 Creating output directory lh.grp_diffs_thickness_med_dur_112013.glmdir Loading y from /ncf/snp/04/SCORE/freesurfer_analysis/glm_analysis/lh_grpDiffs_thickness10.mgh MatrixReadTxT: could not scan value [1][1] File exists ERROR: in uncompressing MATFILE5: 0 x 0, type 0, imagf 0, name '' File exists unsupported matlab format 0 (unknown) File exists Saving design matrix to lh.grp_diffs_thickness_med_dur_112013.glmdir/Xg.dat *** glibc detected *** mri_glmfit: double free or corruption (out): 0x000000003b1279c0 *** ======= Backtrace: ========= /lib64/libc.so.6[0x3229470d7f] /lib64/libc.so.6(cfree+0x4b)[0x32294711db] mri_glmfit[0x99fea2] mri_glmfit[0x5bcbe3] mri_glmfit[0x433608] mri_glmfit[0x433802] mri_glmfit[0x40f2f7] /lib64/libc.so.6(__libc_start_main+0xf4)[0x322941d994] mri_glmfit(__gxx_personality_v0+0x172)[0x40b2ba] ======= Memory map: ======== 00400000-00b86000 r-xp 00400000 00:00 0 00b86000-00c86000 ---p 00b86000 00:00 0 00c86000-26a08000 rwxp 00c86000 00:00 0 35b4b000-3d4fa000 rwxp 35b4b000 00:00 0 [heap] 3229000000-322901c000 r-xp 00000000 fd:00 557272 /lib64/ld-2.5.so 322901c000-322921c000 ---p 322901c000 00:00 0 322921c000-322921d000 r-xp 0001c000 fd:00 557272 /lib64/ld-2.5.so 322921d000-322921e000 rwxp 0001d000 fd:00 557272 /lib64/ld-2.5.so 3229400000-322954d000 r-xp 00000000 fd:00 557376 /lib64/libc-2.5.so 322954d000-322974d000 ---p 0014d000 fd:00 557376 /lib64/libc-2.5.so 322974d000-3229751000 r-xp 0014d000 fd:00 557376 /lib64/libc-2.5.so 3229751000-3229752000 rwxp 00151000 fd:00 557376 /lib64/libc-2.5.so 3229752000-3229757000 rwxp 3229752000 00:00 0 3229800000-3229882000 r-xp 00000000 fd:00 557399 /lib64/libm-2.5.so 3229882000-3229a81000 ---p 00082000 fd:00 557399 /lib64/libm-2.5.so 3229a81000-3229a82000 r-xp 00081000 fd:00 557399 /lib64/libm-2.5.so 3229a82000-3229a83000 rwxp 00082000 fd:00 557399 /lib64/libm-2.5.so 3229c00000-3229c02000 r-xp 00000000 fd:00 557398 /lib64/libdl-2.5.so 3229c02000-3229e02000 ---p 00002000 fd:00 557398 /lib64/libdl-2.5.so 3229e02000-3229e03000 r-xp 00002000 fd:00 557398 /lib64/libdl-2.5.so 3229e03000-3229e04000 rwxp 00003000 fd:00 557398 /lib64/libdl-2.5.so 322a000000-322a014000 r-xp 00000000 fd:00 557401 /lib64/libz.so.1.2.3 322a014000-322a213000 ---p 00014000 fd:00 557401 /lib64/libz.so.1.2.3 322a213000-322a214000 rwxp 00013000 fd:00 557401 /lib64/libz.so.1.2.3 322a400000-322a416000 r-xp 00000000 fd:00 557391 /lib64/libpthread-2.5.so 322a416000-322a615000 ---p 00016000 fd:00 557391 /lib64/libpthread-2.5.so 322a615000-322a616000 r-xp 00015000 fd:00 557391 /lib64/libpthread-2.5.so 322a616000-322a617000 rwxp 00016000 fd:00 557391 /lib64/libpthread-2.5.so 322a617000-322a61b000 rwxp 322a617000 00:00 0 322f000000-322f00d000 r-xp 00000000 fd:00 557370 /lib64/libgcc_s-4.1.2-20080825.so.1 322f00d000-322f20d000 ---p 0000d000 fd:00 557370 /lib64/libgcc_s-4.1.2-20080825.so.1 322f20d000-322f20e000 rwxp 0000d000 fd:00 557370 /lib64/libgcc_s-4.1.2-20080825.so.1 3237200000-32372e6000 r-xp 00000000 fd:00 3309621 /usr/lib64/libstdc++.so.6.0.8 32372e6000-32374e5000 ---p 000e6000 fd:00 3309621 /usr/lib64/libstdc++.so.6.0.8 32374e5000-32374eb000 r-xp 000e5000 fd:00 3309621 /usr/lib64/libstdc++.so.6.0.8 32374eb000-32374ee000 rwxp 000eb000 fd:00 3309621 /usr/lib64/libstdc++.so.6.0.8 32374ee000-3237500000 rwxp 32374ee000 00:00 0 3239600000-3239609000 r-xp 00000000 fd:00 557404 /lib64/libcrypt-2.5.so 3239609000-3239808000 ---p 00009000 fd:00 557404 /lib64/libcrypt-2.5.so 3239808000-3239809000 r-xp 00008000 fd:00 557404 /lib64/libcrypt-2.5.so 3239809000-323980a000 rwxp 00009000 fd:00 557404 /lib64/libcrypt-2.5.so 323980a000-3239838000 rwxp 323980a000 00:00 0 2b1225c45000-2b1225c48000 rwxp 2b1225c45000 00:00 0 2b1225c69000-2b1225c6d000 rwxp 2b1225c69000 00:00 0 2b1225c6d000-2b1225ca2000 r-xs 00000000 fd:00 3703836 /var/run/nscd/dbC3ggnz (deleted) 2b1225d6d000-2b12952bb000 rwxp 2b1225d6d000 00:00 0 7fff44349000-7fff44362000 rwxp 7ffffffe5000 00:00 0 [stack] 7fff443cb000-7fff443ce000 r-xp 7fff443cb000 00:00 0 [vdso] ffffffffff600000-ffffffffffe00000 ---p 00000000 00:00 0 [vsyscall] Aborted -- -- Laura M. Tully, PhD Post-Doctoral Fellow in Psychiatry UC Davis Imaging Research Center 4701 X Street, Sacramento, CA 95817 Phone: (916) 734-7927 Fax: (916) 734-8750 Alumnus of Social Neuroscience & Psychopathology Lab, Harvard University ltully(a)fas.harvard.edu Follow me on twitter: @tully_laura -- *Confidentiality Notice:* This e-mail message, including any attachments, is for the sole use of the intended recipient(s) and may contain confidential and privileged information. Any unauthorized review, use, disclosure or distribution is prohibited. If you are not the intended recipient, please contact the sender by reply e-mail and destroy all copies of the original message. -- -- Laura M. Tully, PhD Post-Doctoral Fellow in Psychiatry UC Davis Imaging Research Center 4701 X Street, Sacramento, CA 95817 Phone: (916) 734-7927 Fax: (916) 734-8750 Alumnus of Social Neuroscience & Psychopathology Lab, Harvard University ltully(a)fas.harvard.edu Follow me on twitter: @tully_laura -- *Confidentiality Notice:* This e-mail message, including any attachments, is for the sole use of the intended recipient(s) and may contain confidential and privileged information. Any unauthorized review, use, disclosure or distribution is prohibited. If you are not the intended recipient, please contact the sender by reply e-mail and destroy all copies of the original message.

2 2

Destrieux atlas in FSL
by Kirstie Whitaker 21 Nov '13

21 Nov '13

Hi Freesurfer and FSL communities, Firstly, I apologise for the double posting. I suspect there are experts in both groups that can help me. I'd like to be able to visualise the Destrieux atlas that freesurfer uses at the same time as my MNI standard space results in FSL. Is that possible? I'm afraid I've not been successful in my search of the mailing list archives. Thank you for your help. Kirstie -- Kirstie Whitaker, PhD Research Associate tel: +44 7583 535 307 Department of Psychiatry University of Cambridge Herchel Smith Building Robinson Way Cambridge CB2 0SZ

2 1

optseq2 for a (sparse) 1-back task
by Chris McNorgan 21 Nov '13

21 Nov '13

Hi all, I figured I would have some luck tapping expertise on using optseq2 in this mailing list. There are 8 trials times 8 conditions in an event-related design. We are using a 1-back task just to keep participants paying attention. Thus, we are not attempting to match the number of 1-back-yes and 1-back-no trials. In fact, ideally we would have exactly 1 1-back-yes trial for each of the 8 conditions. The plan would be to arrive at a paradigm schedule where, for each of the 8 conditions, there is at least 1 set of consecutive trials for that condition. In this way, if we use only 7 stimuli, then the consecutive trials can use the same item, and that would constitute the 1-back-yes event for that condition as in: 1 a1.bmp 2 b1.bmp 0 NULL 3 c1.bmp 2 b2.bmp <-- 0 NULL 2 b2.bmp <--first instance of consecutive events from same condition assigned same item to create a 1-back-yes event The problem is that I have not discovered any means within optseq to enforce the criterion that at two consecutive events are from the same condition for all 8 conditions. The best I can do is create a large enough pool of paradigm files, and sort through them in the hopes of finding several that meet this requirement. Can optseq2 be forced to meet this requirement, either as a normal command line switch, or else through some clever condition labeling trickery? Thanks, Chris

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