Freesurfer October 2008

freesurfer@nmr.mgh.harvard.edu

71 participants
110 discussions

Start a nNew thread

Permissions Denied (Control Points)
by Kellen Mobilia 09 Oct '08

09 Oct '08

Help-

1 0

FW: [Freesurfer] Talairach_with_skull.lta accuracy
by Jeff Sadino 09 Oct '08

09 Oct '08

Hello, I am trying to validate the ICV value by lining up the nu_noneck image with the RB image. These are the commands I am using: mri_convert nu_noneck.mgz --apply_transform transforms/talairach_with_skull.lta nu_noneck5.mgz tkmedit -f nu_noneck5.mgz & tkmedit -f /apps/freesurfer/average/RB_all_withskull_2006-02-15.gca & When I compare these two, they don't line up. The RB image is all blurry (Maybe it's a combination of several subjects and is supposed to be like that). Any suggestions as to how to line up the noneck transformed imaged to the RB image is greatly appreciated! I also attached three files to show what's happening in my transform process. Thanks, Jeff ---------------------------------------- > Date: Fri, 5 Sep 2008 20:31:55 -0400 > From: fischl(a)nmr.mgh.harvard.edu > To: jsadino(a)hotmail.com > Subject: RE: [Freesurfer] Talairach_with_skull.lta accuracy > > can you send this to the list so others can respond? Where did you get the > tmp.reg file? If you have transformed them already you shouldn't need > tkregister2 (or a registration file in general) to compare them. > > On Fri, 5 Sep 2008, > Jeff Sadino wrote: > >> >> Thank you very much for your help. Your knowledge and quick responses are very much appreciated. >> >> I converted my image into transformed space, but it still does not match the RB image. These are the commands I used: >> >> mri_convert nu_noneck.mgz --apply_transform > transforms/talairach_with_skull.lta nu_noneckt3.mgz >> >> tkregister2 --targ > /apps/freesurfer/average/RB_all_withskull_2006-02-15.gca --reg tmp.reg > --lta /data/040002_S04_bakj/mri/transforms/talairach_with_skull.lta --mov > /data/040002_S04_bakj/mri/nu_noneckt3.mgz >> >> Again, my main goal is to verify the ICV value, so to do that I need to verify that the lta file is transforming the nu_noneck image onto the RB image...I think. >> >> What step am I missing? >> >> Thank you >> Jeff >> >> >> ---------------------------------------- >>> Date: Fri, 5 Sep 2008 13:09:21 -0400 >>> From: fischl(a)nmr.mgh.harvard.edu >>> To: jsadino(a)hotmail.com >>> CC: greve(a)nmr.mgh.harvard.edu >>> Subject: RE: [Freesurfer] Talairach_with_skull.lta accuracy >>> >>> that won't generate a transformed volume, just display it in the >>> transformed space. That is, the volumes you sent were in their native >>> space and would have no reason to be aligned. The neck removal seems to >>> have been fine in both cases. >>> On Thu, 4 Sep 2008, Jeff Sadino wrote: >>> >>>> >>>> ps, sorry, I forgot to include my command: >>>> >>>> >>>> tkregister2 --targ /apps/freesurfer/average/RB_all_withskull_2006-02-15.gca --mov /data/${file}/mri/nu_noneck.mgz --lta /data/${file}/mri/transforms/talairach_with_skull.lta --reg tmp.reg >>>> >>>> I believe this will transform the nu_neck.mgz image via the talairach matrix, so that it matches the RB image. >>>> >>>> btw, we processed the scans using freesurfer 3, but the tkregister2 that came with freesurfer 3 did not include the --lta tag, so I am using the tkregister2 from freesurfer4. >>>> >>>> Thanks! >>>> >>>> ---------------------------------------- >>>>> Date: Thu, 4 Sep 2008 21:49:27 -0400 >>>>> From: fischl(a)nmr.mgh.harvard.edu >>>>> To: jsadino(a)hotmail.com >>>>> CC: freesurfer(a)nmr.mgh.harvard.edu >>>>> Subject: Re: [Freesurfer] Talairach_with_skull.lta accuracy >>>>> >>>>> Hi Jeff, >>>>> >>>>> why do you say that one is bad? We don't resample, so neither one will be >>>>> aligned with the atlas unless you apply the transform. >>>>> >>>>> We are aware that the eTIV measures have more variability than we would >>>>> like, and we're looking into it, but don't have a good solution from only >>>>> T1 images. >>>>> >>>>> cheers, >>>>> Bruce >>>>> On Thu, 4 Sep 2008, >>>>> Jeff Sadino wrote: >>>>> >>>>>> >>>>>> Hello all, >>>>>> >>>>>> According to aseg.stats, our ICV values fluctuate on the same subject across scans, sometimes by up to 10% I have been validating the talairach_with_skull.lta accuracy by comparing the nu_neck.mgz to the RB_all_withskull_2006-02-15.gca template. I read on the boards that it is not all that necessary that the brains match, that it is more important for the skulls to match. The RB is blurred (I don't fully understand where the RB comes from, I think from 22 subjects that the freesurfers combined into one?), so am I matching up the whitest areas of the RB to the skull on the nu_neck image? I noticed also that the bottom of the brains rarely match at all to the RB image...is this normal? Are there any other parts of the image that need to be checked for validity? >>>>>> >>>>>> I attached two nu_noneck.mgz files, one was good and the other bad. Are these images typical? >>>>>> >>>>>> Thank you for your help, >>>>>> Jeff >>>>>> >>>>>> _________________________________________________________________ >>>>>> See how Windows Mobile brings your life together?at home, work, or on the go. >>>>>> http://clk.atdmt.com/MRT/go/msnnkwxp1020093182mrt/direct/01/ >>>> >>>> _________________________________________________________________ >>>> Get more out of the Web. Learn 10 hidden secrets of Windows Live. >>>> http://windowslive.com/connect/post/jamiethomson.spaces.live.com-Blog-cns!5… >>>> >>>> >>>> >> >> _________________________________________________________________ >> See how Windows connects the people, information, and fun that are part of your life. >> http://clk.atdmt.com/MRT/go/msnnkwxp1020093175mrt/direct/01/ >> >> >> _________________________________________________________________ See how Windows Mobile brings your life together—at home, work, or on the go. http://clk.atdmt.com/MRT/go/msnnkwxp1020093182mrt/direct/01/

3 2

Re: [Freesurfer] Fwd: glmfit and surfcluster
by Doug Greve 09 Oct '08

09 Oct '08

Hi Joost, can you look at the original sig map thresholded at 1.3 to see if there is a cluster of that size? Hi, version 4.0.5 after: mri_glmfit --y qdec/N83_only_boys.mgh --fsgd qdec/qdec.fsgd doss --surf fsaverage lh --C qdec/diff-1-2.txt --sim perm 1000 1.3 stats/newperm_lh --glmdir stats/newperm_lh/ i got an csd file (attached), "stats/newperm_lh and --glmdir stats/newperm_lh/" were ignored (as expected). i want to use this .csd file to inspect clusters: mri_surfcluster --in qdec/N83_only_boys.mgh --csd stats/perm_lh-diff-1-2.txt.csd --cwsig stats/cluster_map.mgh --sum qdec/summary.txt which gives the following output (summary.txt is attached): ******************** # Cluster Growing Summary (mri_surfcluster) # $Id: mri_surfcluster.c,v 1.39 2007/07/31 00:34:19 greve Exp $ # $Id: mrisurf.c,v 1.557.2.13 2008/05/23 00:06:05 nicks Exp $ # CreationTime 2008/10/08-08:52:02-GMT # cmdline mri_surfcluster --in qdec/N83_only_boys.mgh --csd stats/perm_lh-diff-1-2.txt.csd --cwsig stats/cluster_map.mgh --sum q dec/summary.txt # cwd /home/neuro/joost/redes/maranon/freesurfer # sysname Linux # hostname fedoraneuro.lim.local # machine x86_64 # FixVertexAreaFlag = 1 # # Input qdec/N83_only_boys.mgh # Frame Number 0 # srcsubj fsaverage # hemi lh # surface white # SUBJECTS_DIR /home/neuro/joost/redes/maranon/freesurfer # Minimum Threshold 1.3 # Maximum Threshold infinity # Threshold Sign abs # AdjustThreshWhenOneTail 1 # Area Threshold 0 mm^2 # CSD thresh 1.300000 # CSD nreps 1000 # CSD simtype perm # CSD contrast diff-1-2.txt # CSD confint 90.000000 # Overall max 4.37947 at vertex 3484 # Overall min 1.58841e-06 at vertex 136280 # NClusters 1 # Total Cortical Surface Area 65416.6 (mm^2) # FixMNI = 1 # # ClusterNo Max VtxMax Size(mm^2) TalX TalY TalZ CWP CWPLow CWPHi NVtxs 1 4.379 3484 98860.55 -34.7 11.5 -6.9 0.00100 0.00000 0.00200 155267 [neuro@fedoraneuro qdec]$ cd .. [neuro@fedoraneuro freesurfer]$ mri_surfcluster --in qdec/N83_only_boys.mgh --csd stats/perm_lh-diff-1-2.txt.csd --cwsig stats/cluster_map.mgh --sum qdec/summary.txt thsign = abs, id = 0 version $Id: mri_surfcluster.c,v 1.39 2007/07/31 00:34:19 greve Exp $ hemi = lh srcid = qdec/N83_only_boys.mgh paint srcsubjid = fsaverage srcsurf = white srcframe = 0 thsign = abs thmin = 1.3 thmax = -1 fdr = -1 minarea = 0 xfmfile = talairach.xfm nth = -1 sumfile = qdec/summary.txt subjectsdir = /home/neuro/joost/redes/maranon/freesurfer FixMNI = 1 ------------- XFM matrix (RAS2RAS) --------------- /home/neuro/joost/redes/maranon/freesurfer/fsaverage/mri/transforms/talairach.xfm 1.000 0.000 0.000 0.000; 0.000 1.000 0.000 0.000; 0.000 0.000 1.000 0.000; 0.000 0.000 0.000 1.000; ---------------------------------------------------- Reading source surface /home/neuro/joost/redes/maranon/freesurfer/fsaverage/surf/lh.white reading group avg surface area 822 cm^2 from file Reading in average area /home/neuro/joost/redes/maranon/freesurfer/fsaverage/surf/lh.white.avg.area.mgh Done reading source surface Computing metric properties Loading source values number of voxels in search space = 163842 Done loading source values (nvtxs = 163842) overall max = 4.37947 at vertex 3484 overall min = 1.58841e-06 at vertex 136280 surface nvertices 163842 surface area 65417.097656 surface area 65416.648438 NOT Adjusting threshold for 1-tailed test Searching for Clusters ... thmin=1.300000 (1.300000), thmax=-1.000000 (-1), thsignid=0, minarea=0.000000 Found 1 clusters Max cluster size 98860.554688 INFO: fixing MNI talairach coordinates Saving cluster pval stats/cluster_map.mgh ************************ - the cluster size seems way to big. your thoughts? - for visualizing the cluster i thought: loading the --cwsig file (cluster_map.mgh) as an overlay on the (inflated) surface of fsaverage and then configure the overlay to set the minimal threshold to 1.3. your thoughts? thanks! -joost On Tue, Oct 7, 2008 at 7:37 PM, Pratap Kunwar <pkunwar(a)nmr.mgh.harvard.edu <mailto:pkunwar@nmr.mgh.harvard.edu>> wrote: Hi, You have a continuous variable "age" in your fsgd file. As far as i know, permutation can't have continuous variables. Try mc-z that should work. pratap > Hi, > > using version 4.0.5 > > i ran the command: > mri_glmfit --y qdec/N83_only_boys.mgh --fsgd qdec/qdec.fsgd doss --surf > fsaverage lh --C qdec/diff-1-2.txt --glmdir stats/ > > it worked. > > then i did: > > mri_glmfit --y qdec/N83_only_boys.mgh --fsgd qdec/qdec.fsgd doss --surf > fsaverage lh --C qdec/diff-1-2.txt --sim perm 1000 1.3 stats/newperm_lh > --glmdir stats/newperm_lh/ > > output: > > gdfReadHeader: reading qdec/qdec.fsgd > INFO: DeMeanFlag keyword not found, DeMeaning will NOT be done. > Continuous Variable Means (all subjects) > 0 Age 15.4819 > Class Means of each Continuous Variable > 1 control 15.0882 > 2 patient 15.7551 > INFO: gd2mtx_method is doss > Reading source surface > /home/neuro/joost/redes/maranon/freesurfer/fsaverage/surf/lh.white > reading group avg surface area 822 cm^2 from file > Reading in average area > /home/neuro/joost/redes/maranon/freesurfer/fsaverage/surf/lh.white.avg.area.mgh > simbase stats/perm_lh > 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.138.2.1 <http://1.138.2.1> 2007/09/12 15:38:19 nicks Exp $ > cwd /home/neuro/joost/redes/maranon/freesurfer > cmdline mri_glmfit --y qdec/N83_only_boys.mgh --fsgd qdec/qdec.fsgd doss > --surf fsaverage lh --C qdec/diff-1-2.txt --sim perm 1000 1.3 > stats/perm_lh > --glmdir stats/ > sysname Linux > hostname fedoraneuro.lim.local > machine x86_64 > user neuro > FixVertexAreaFlag = 1 > UseMaskWithSmoothing 1 > OneSampleGroupMean 0 > y /home/neuro/joost/redes/maranon/freesurfer/qdec/N83_only_boys.mgh > logyflag 0 > usedti 0 > FSGD qdec/qdec.fsgd > glmdir stats/ > DoFFx 0 > Loading y from > /home/neuro/joost/redes/maranon/freesurfer/qdec/N83_only_boys.mgh > Matrix condition is 404.504 > search space = 82220 > ERROR: design matrix is not orthogonal, cannot be used with permutation. > If this something you really want to do, run with --perm-force > > mri_info N83_only_boys.mgh gives 83 frames. > > what´s wrong with my design matrix? > > thanks! > -joost > _______________________________________________ > Freesurfer mailing list > Freesurfer(a)nmr.mgh.harvard.edu <mailto:Freesurfer@nmr.mgh.harvard.edu> > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer version 4.0.5. <http://4.0.5.> ------------------------------------------------------------------------ # Cluster Growing Summary (mri_surfcluster) # $Id: mri_surfcluster.c,v 1.39 2007/07/31 00:34:19 greve Exp $ # $Id: mrisurf.c,v 1.557.2.13 2008/05/23 00:06:05 nicks Exp $ # CreationTime 2008/10/08-08:52:02-GMT # cmdline mri_surfcluster --in qdec/N83_only_boys.mgh --csd stats/perm_lh-diff-1-2.txt.csd --cwsig stats/cluster_map.mgh --sum qdec/summary # cwd /home/neuro/joost/redes/maranon/freesurfer # sysname Linux # hostname fedoraneuro.lim.local # machine x86_64 # FixVertexAreaFlag = 1 # # Input qdec/N83_only_boys.mgh # Frame Number 0 # srcsubj fsaverage # hemi lh # surface white # SUBJECTS_DIR /home/neuro/joost/redes/maranon/freesurfer # Minimum Threshold 1.3 # Maximum Threshold infinity # Threshold Sign abs # AdjustThreshWhenOneTail 1 # Area Threshold 0 mm2 # CSD thresh 1.300000 # CSD nreps 1000 # CSD simtype perm # CSD contrast diff-1-2.txt # CSD confint 90.000000 # Overall max 4.37947 at vertex 3484 # Overall min 1.58841e-06 at vertex 136280 # NClusters 1 # Total Cortical Surface Area 65416.6 (mm2) # FixMNI = 1 # # ClusterNo Max VtxMax Size(mm2) TalX TalY TalZ CWP CWPLow CWPHi NVtxs 1 4.379 3484 98860.55 -34.7 11.5 -6.9 0.00100 0.00000 0.00200 155267 ------------------------------------------------------------------------ 1 -1 ------------------------------------------------------------------------ _______________________________________________ Freesurfer mailing list Freesurfer(a)nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer j janssen wrote: -- Douglas N. Greve, Ph.D. MGH-NMR Center greve(a)nmr.mgh.harvard.edu Phone Number: 617-724-2358 Fax: 617-726-7422 In order to help us help you, please follow the steps in: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting

1 0

glmfit and surfcluster
by j janssen 09 Oct '08

09 Oct '08

Hi, version 4.0.5. after: mri_glmfit --y qdec/N83_only_boys.mgh --fsgd qdec/qdec.fsgd doss --surf fsaverage lh --C qdec/diff-1-2.txt --sim perm 1000 1.3 stats/newperm_lh --glmdir stats/newperm_lh/ i got an csd file (attached), "stats/newperm_lh and --glmdir stats/newperm_lh/" were ignored (as expected). i want to use this .csd file to inspect clusters: mri_surfcluster --in qdec/N83_only_boys.mgh --csd stats/perm_lh-diff-1-2.txt.csd --cwsig stats/cluster_map.mgh --sum qdec/summary.txt which gives the following output (summary.txt is attached): ******************** # Cluster Growing Summary (mri_surfcluster) # $Id: mri_surfcluster.c,v 1.39 2007/07/31 00:34:19 greve Exp $ # $Id: mrisurf.c,v 1.557.2.13 2008/05/23 00:06:05 nicks Exp $ # CreationTime 2008/10/08-08:52:02-GMT # cmdline mri_surfcluster --in qdec/N83_only_boys.mgh --csd stats/perm_lh-diff-1-2.txt.csd --cwsig stats/cluster_map.mgh --sum q dec/summary.txt # cwd /home/neuro/joost/redes/maranon/freesurfer # sysname Linux # hostname fedoraneuro.lim.local # machine x86_64 # FixVertexAreaFlag = 1 # # Input qdec/N83_only_boys.mgh # Frame Number 0 # srcsubj fsaverage # hemi lh # surface white # SUBJECTS_DIR /home/neuro/joost/redes/maranon/freesurfer # Minimum Threshold 1.3 # Maximum Threshold infinity # Threshold Sign abs # AdjustThreshWhenOneTail 1 # Area Threshold 0 mm^2 # CSD thresh 1.300000 # CSD nreps 1000 # CSD simtype perm # CSD contrast diff-1-2.txt # CSD confint 90.000000 # Overall max 4.37947 at vertex 3484 # Overall min 1.58841e-06 at vertex 136280 # NClusters 1 # Total Cortical Surface Area 65416.6 (mm^2) # FixMNI = 1 # # ClusterNo Max VtxMax Size(mm^2) TalX TalY TalZ CWP CWPLow CWPHi NVtxs 1 4.379 3484 98860.55 -34.7 11.5 -6.9 0.00100 0.00000 0.00200 155267 [neuro@fedoraneuro qdec]$ cd .. [neuro@fedoraneuro freesurfer]$ mri_surfcluster --in qdec/N83_only_boys.mgh --csd stats/perm_lh-diff-1-2.txt.csd --cwsig stats/cluster_map.mgh --sum qdec/summary.txt thsign = abs, id = 0 version $Id: mri_surfcluster.c,v 1.39 2007/07/31 00:34:19 greve Exp $ hemi = lh srcid = qdec/N83_only_boys.mgh paint srcsubjid = fsaverage srcsurf = white srcframe = 0 thsign = abs thmin = 1.3 thmax = -1 fdr = -1 minarea = 0 xfmfile = talairach.xfm nth = -1 sumfile = qdec/summary.txt subjectsdir = /home/neuro/joost/redes/maranon/freesurfer FixMNI = 1 ------------- XFM matrix (RAS2RAS) --------------- /home/neuro/joost/redes/maranon/freesurfer/fsaverage/mri/transforms/talairach.xfm 1.000 0.000 0.000 0.000; 0.000 1.000 0.000 0.000; 0.000 0.000 1.000 0.000; 0.000 0.000 0.000 1.000; ---------------------------------------------------- Reading source surface /home/neuro/joost/redes/maranon/freesurfer/fsaverage/surf/lh.white reading group avg surface area 822 cm^2 from file Reading in average area /home/neuro/joost/redes/maranon/freesurfer/fsaverage/surf/lh.white.avg.area.mgh Done reading source surface Computing metric properties Loading source values number of voxels in search space = 163842 Done loading source values (nvtxs = 163842) overall max = 4.37947 at vertex 3484 overall min = 1.58841e-06 at vertex 136280 surface nvertices 163842 surface area 65417.097656 surface area 65416.648438 NOT Adjusting threshold for 1-tailed test Searching for Clusters ... thmin=1.300000 (1.300000), thmax=-1.000000 (-1), thsignid=0, minarea=0.000000 Found 1 clusters Max cluster size 98860.554688 INFO: fixing MNI talairach coordinates Saving cluster pval stats/cluster_map.mgh ************************ - the cluster size seems way to big. your thoughts? - for visualizing the cluster i thought: loading the --cwsig file (cluster_map.mgh) as an overlay on the (inflated) surface of fsaverage and then configure the overlay to set the minimal threshold to 1.3. your thoughts? thanks! -joost On Tue, Oct 7, 2008 at 7:37 PM, Pratap Kunwar <pkunwar(a)nmr.mgh.harvard.edu>wrote: > > Hi, > > You have a continuous variable "age" in your fsgd file. As far as i know, > permutation can't have continuous variables. > > Try mc-z that should work. > > pratap > > > Hi, > > > > using version 4.0.5 > > > > i ran the command: > > mri_glmfit --y qdec/N83_only_boys.mgh --fsgd qdec/qdec.fsgd doss --surf > > fsaverage lh --C qdec/diff-1-2.txt --glmdir stats/ > > > > it worked. > > > > then i did: > > > > mri_glmfit --y qdec/N83_only_boys.mgh --fsgd qdec/qdec.fsgd doss --surf > > fsaverage lh --C qdec/diff-1-2.txt --sim perm 1000 1.3 stats/newperm_lh > > --glmdir stats/newperm_lh/ > > > > output: > > > > gdfReadHeader: reading qdec/qdec.fsgd > > INFO: DeMeanFlag keyword not found, DeMeaning will NOT be done. > > Continuous Variable Means (all subjects) > > 0 Age 15.4819 > > Class Means of each Continuous Variable > > 1 control 15.0882 > > 2 patient 15.7551 > > INFO: gd2mtx_method is doss > > Reading source surface > > /home/neuro/joost/redes/maranon/freesurfer/fsaverage/surf/lh.white > > reading group avg surface area 822 cm^2 from file > > Reading in average area > > > /home/neuro/joost/redes/maranon/freesurfer/fsaverage/surf/lh.white.avg.area.mgh > > simbase stats/perm_lh > > 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.138.2.1 2007/09/12 15:38:19 nicks Exp $ > > cwd /home/neuro/joost/redes/maranon/freesurfer > > cmdline mri_glmfit --y qdec/N83_only_boys.mgh --fsgd qdec/qdec.fsgd doss > > --surf fsaverage lh --C qdec/diff-1-2.txt --sim perm 1000 1.3 > > stats/perm_lh > > --glmdir stats/ > > sysname Linux > > hostname fedoraneuro.lim.local > > machine x86_64 > > user neuro > > FixVertexAreaFlag = 1 > > UseMaskWithSmoothing 1 > > OneSampleGroupMean 0 > > y /home/neuro/joost/redes/maranon/freesurfer/qdec/N83_only_boys.mgh > > logyflag 0 > > usedti 0 > > FSGD qdec/qdec.fsgd > > glmdir stats/ > > DoFFx 0 > > Loading y from > > /home/neuro/joost/redes/maranon/freesurfer/qdec/N83_only_boys.mgh > > Matrix condition is 404.504 > > search space = 82220 > > ERROR: design matrix is not orthogonal, cannot be used with permutation. > > If this something you really want to do, run with --perm-force > > > > mri_info N83_only_boys.mgh gives 83 frames. > > > > what´s wrong with my design matrix? > > > > thanks! > > -joost > > _______________________________________________ > > Freesurfer mailing list > > Freesurfer(a)nmr.mgh.harvard.edu > > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > >

1 1

[Fwd: Re: [Freesurfer] A methodological question]
by Michael Harms 09 Oct '08

09 Oct '08

I'm trying to post this again to the list, since the papers attached (in response to a post from Pedro) may be of interest to others, but they appear to have not been distributed to the overall list in my previous reply to Pedro. Does the list server perhaps reject certain posts based on number or size of attachments? cheers, Mike H. -------- Forwarded Message -------- From: Michael Harms <mharms(a)conte.wustl.edu> To: Pedro Paulo de Magalhães Oliveira Junior <ppj(a)netfilter.com.br> Cc: Freesurfer Mailing List <freesurfer(a)nmr.mgh.harvard.edu> Subject: Re: [Freesurfer] A methodological question Date: Fri, 03 Oct 2008 13:29:29 -0500 Hi Pedro, Attached are 4 papers that I'm aware of that explicitly compare aspects of FS to other approaches (Lee 2006; Tae 2008; Pengas 2008; Ju 2005) Those particularly interested in TICV should take a look at the Pengas paper (which is timely, given the recent discussion on the list regarding the reliability of the "eTIV" measure). cheers, Mike H. On Fri, 2008-10-03 at 13:50 -0200, Pedro Paulo de Magalhães Oliveira Junior wrote: > Some weeks ago,I have posted here a compilation of the main papers > that describe FreeSurfer implementation. > > > I was wondering if someone has researched in comparison between > FreeSurfer methodology and other ones. (For instance: > doi:10.1016/j.neuroimage.2004.12.056) > > > Best Regards > > > PPJ > > -- > ----------------------------------------------------------- > Pedro Paulo de M. Oliveira Junior > Diretor de Operações > Netfilter & SpeedComm Telecom > > > _______________________________________________ > Freesurfer mailing list > Freesurfer(a)nmr.mgh.harvard.edu > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer _______________________________________________ Freesurfer mailing list Freesurfer(a)nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

1 0

Voxel-wise analysis of subcortical strucutures
by Juergen Haenggi 09 Oct '08

09 Oct '08

Dear FS experts I am wondering whether it is possible to do an voxel-wise analysis on the subcortical segmentations similar to a voxel-based morphometry (VBM) in SPM or FSL. Because there is no .mgz file that encodes the probability that a particular voxel belongs to a particular tissue type (GM or WM) in the mri folder of the subject, I guess that it is not possible (yet), right? Any advise is appreciated Regards Juergen ----------------------------------------------------------- Juergen Haenggi Ph.D. (Dr. des.) Division of 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(a)psychologie.uzh.ch (email) http://www.psychologie.uzh.ch/neuropsy/ (website) http://www.juergenhaenggi.ch (private website) -----------------------------------------------------------

2 1

Problem tkmedit after recon-all -long
by Juergen Haenggi 09 Oct '08

09 Oct '08

Dear FS experts I have processed longitudinal data with recon-all -long of three patients with atrophy in the caudate nucleus. When I produce a figure from time point 1 and one figure from time point 2 in order to visualize the atrophy, the figures look different due the different orientations of the two images in native space. How can I reorient the two T1.mgz images and the aseg and the surfaces in a way that the AC-PC plane is horizontal in the native space in order to produce figures that show both images at the same plane. Is there any option that can be used with tkmedit to accomplish this? Thanks in advance Best regards Juergen ----------------------------------------------------------- Juergen Haenggi Ph.D. (Dr. des.) Division of 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(a)psychologie.uzh.ch (email) http://www.psychologie.uzh.ch/neuropsy/ (website) http://www.juergenhaenggi.ch (private website) -----------------------------------------------------------

2 1

tksurfer display tcl error & FS training
by Derin Cobia 08 Oct '08

08 Oct '08

Two questions: 1. We installed the latest version of 32-bit FS on several workstations running CentOS 5.1 in our lab. However, once we ran the CentOS updates (~280 of them through up2date), tkmedit now fails to open and tksurfer exhibits behavior exactly as described in this message from the mail archive: http://www.mail-archive.com/freesurfer@nmr.mgh.harvard.edu/msg07570.html Having attempted to troubleshoot it I discovered that a blt lib was missing, but even with re-installation of this tksurfer continues to not display correctly. I believe it's a tcl problem, but I don't understand why performing incremental updates to CentOS would break FS. More importantly, I'm not sure which of the updates does "the breaking" so I can just avoid it. Any ideas? 2. We are in the process of training new research assistants in how to use FS. Typically we have approached this in a more qualitative way by working side-by-side with them to guide their edits and processing. I was wondering if there is a quantitative way of comparing their work that you (or others) have used when training "new recruits" in FS. Maybe something like was done in the Han et al (2006) paper, such as compare their thickness (or other) maps with someone we know does well with manual interventions? In essence, we're trying to think of a rigorous way to train new people in FS, thanks! -Derin

2 1

qcache commands
by Mishkin Derakhshan 08 Oct '08

08 Oct '08

Hi, I am trying to understand what happens when you run the recon-all -qcache command I noticed these three commands get run: mris_preproc --s tester2 --hemi lh --meas thickness --target fsaverage --out lh.thickness.fsaverage.mgh mri_surf2surf --hemi lh --srcsubject tester2 --surfreg sphere.reg --sval /lab1/FreeSurfer/4.0.5/ms-mni/testing/tester2/surf/lh.thickness --sfmt curv --trgsubject fsaverage --tval ./tmp.mris_preproc.12176/tester2.1.mgh --noreshape mri_concat ./tmp.mris_preproc.12176/tester2.1.mgh --o lh.thickness.fsaverage.mgh My questions are: 1. Am I right to assume that the first command (mris_preproc) is just a wrapper script that calls the 2nd and 3rd commands? 2a. In the case above, does the 3rd command (mri_concat) simply copy the contents of tester2.1.mgh into lh.thickness.fsaverage.mgh. 2b. Is there a case when mri_concat will actually concat more than one thing into lh.thickness.fsaverage.mgh? thanks, mishkin

2 3

Re: [Freesurfer] qdec.table.dat error
by Mishkin Derakhshan 08 Oct '08

08 Oct '08

Hi Nick, I tried unalias -a and then runnig qdec but it didn't work either. I think the problem is with bash vs. sh in ubuntu (possible all debian systems). ie. /bin/sh is not the same as /bin/bash so not all bash commands will work. I did this: sudo mv /bin/sh /bin/sh.bak sudo ln -s /bin/bash /bin/sh and then it worked. I'm hoping this doesn't have repercussion for my other programs though. thanks, mishkin On Tue, Oct 7, 2008 at 11:20 PM, Nick Schmansky <nicks(a)nmr.mgh.harvard.edu> wrote: > yes, this sounds like a tough one. the VerifySubjects routine that > produces that error message just runs the 'ls' command on that subject > directory name, so that is why i asked. if you run the 'ls' command on > '/lab1/FreeSurfer/4.0.5/ms-mni/testing/tester', does it return a correct > response? ie a non-zero status? > > what does 'which ls' say? possibly it is aliased to something that is > causing a problem. > > does qdec work despite the error? the code is setup to just print the > error, but continue anyway. if the subject doesnt exist, it will fail > downstream (during mri_glmfit). > > i'll try to think of other things. > > nick > > > On Tue, 2008-10-07 at 21:26 -0400, Mishkin Derakhshan wrote: >> Hi Nick, >> Yes I've checked most of the obvious things I could think of. >> The directories all exist, permissions on all files are 775, i will >> try 777 but don't suspect that will make a difference. >> My env variables seem to be set properly as well. >> >> opus[/lab1/FreeSurfer/4.0.5/ms-mni/testing/tester]$ . >> /lab1/FreeSurfer/4.0.5/quarantine/SetUpFreeSurfer.sh >> -------- freesurfer-Linux-centos4_x86_64-stable-pub-v4.0.5 -------- >> Setting up environment for FreeSurfer/FS-FAST (and FSL) >> FREESURFER_HOME /lab1/FreeSurfer/4.0.5/quarantine >> FSFAST_HOME /lab1/FreeSurfer/4.0.5/quarantine/fsfast >> FSF_OUTPUT_FORMAT nii >> SUBJECTS_DIR /lab1/FreeSurfer/4.0.5/ms-mni/testing/ >> MNI_DIR /lab1/FreeSurfer/4.0.5/quarantine/mni >> >> I get the same error when trying to load the tutorial qdec.table.dat as well. >> >> My analysis is stuck at this point so any help is greatly appreciated. >> >> thanks, >> mishkin >> >> >> On Tue, Oct 7, 2008 at 7:27 PM, Nick Schmansky >> <nicks(a)nmr.mgh.harvard.edu> wrote: >> > Mishkin, >> > >> > Does this directory exist? >> > >> > /lab1/FreeSurfer/4.0.5/ms-mni/testing/tester >> > >> > The error message you are seeing is checking for the existence of that >> > directory. >> > >> > Nick >> > >> > >> > On Tue, 2008-10-07 at 14:30 -0400, Mishkin Derakhshan wrote: >> >> Hi, >> >> When I try to load a qdec.table.dat I always get the same error: >> >> >> >> sh: Syntax error: Bad fd number >> >> ERROR: QdecProject::VerifySubjects: Couldn't find subject 'tester' in >> >> SUBJECTS_DIR >> >> >> >> I've tried digging around through the archives and still can't seem to >> >> find a solution for this problem. >> >> I've tried tcsh and the bash shell. >> >> Here is a bit more information and my qdec.table.dat is attached. Any >> >> ideas on how to get qdec up and running would be very useful. >> >> thanks, >> >> mishkin >> >> >> >> >> >> FREESURFER_HOME: /lab1/FreeSurfer/4.0.5/quarantine >> >> Build stamp: freesurfer-Linux-centos4_x86_64-stable-pub-v4.0.5 >> >> Debian version: lenny/sid >> >> Kernel info: Linux 2.6.24-18-generic x86_64 >> >> >> >> opus:/lab1/FreeSurfer/4.0.5/ms-mni/testing> qdec & >> >> opus:/lab1/FreeSurfer/4.0.5/ms-mni/testing> >> >> Reading /lab1/FreeSurfer/4.0.5/quarantine/lib/tcl/tkUtils.tcl >> >> >> >> Using /lab1/FreeSurfer/4.0.5/quarantine/lib/tcl/fsgdfPlot.tcl >> >> >> >> Loading data table /lab1/FreeSurfer/4.0.5/ms-mni/testing/qdec/qdec.table.dat... >> >> Setting SUBJECTS_DIR to '/lab1/FreeSurfer/4.0.5/ms-mni/testing' >> >> Number of columns: 4 >> >> fsid column: 1 >> >> Number of factors: 3 >> >> Number of subjects: 2 >> >> Reading discrete factor levels from config file >> >> /lab1/FreeSurfer/4.0.5/ms-mni/testing/qdec/gender.levels...'Male','Female', >> >> done. >> >> >> >> Data table /lab1/FreeSurfer/4.0.5/ms-mni/testing/qdec/qdec.table.dat loaded. >> >> SUBJECTS_DIR is '/lab1/FreeSurfer/4.0.5/ms-mni/testing' >> >> sh: Syntax error: Bad fd number >> >> ERROR: QdecProject::VerifySubjects: Couldn't find subject 'tester' in >> >> SUBJECTS_DIR >> >> sh: Syntax error: Bad fd number >> >> ERROR: QdecProject::VerifySubjects: Couldn't find subject 'tester2' in >> >> SUBJECTS_DIR >> >> Error loading the data table. >> >> _______________________________________________ >> >> Freesurfer mailing list >> >> Freesurfer(a)nmr.mgh.harvard.edu >> >> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer >> > >> > >> >> > >

2 1

← Newer
1
...
5
6
7
8
9
10
11
Older →

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

2006

2005

2004

Freesurfer October 2008