Freesurfer

freesurfer@nmr.mgh.harvard.edu

17 participants
30099 discussions

Start a nNew thread

White matter completely filled in
by Jonathan Holt 21 Jul '14

21 Jul '14

The problem is detailed in pictures, I'm emailing to find out how to avoid this problem in the future.

1 0

Edits to Pial and CP at the same time?
by Mark Fletcher 21 Jul '14

21 Jul '14

Dear FreeSurfers, In the training that I attended in May, I believe it was demonstrated by Allison that you could fix multiple problems at the same time. I understand how to make the edits at the same time, but am unsure what command to run to integrate all of these changes concurrently. Since I currently do not know how to do this, it is taking me multiple days to fix some problems on subjects with multiple types of errors. Currently I first manually delete some dura on the brainmask.mgz (thus I would run recon-all -autorecon-pial -subjid xxxx). After examining the output the next day, I then add control points in the anterior temporal lobe or other places if needed on the same subject ( recon-all -autorecon2-cp -autorecon3 -subjid xxxx ). I was wondering how to both manually delete the dura, and add cp at the same time for one processing event (thus saving me a lot of time)? Would the correct command line be: recon-all -autorecon-pial autorecon2-cp -autorecon3 -subjid xxxx ? Also, is it possible to edit the pial, add white matter, and control points all at the same time for one recon-all processing event? Would the correct command line be recon-all -autorecon-pial autorecon2*-cp-wm* -autorecon3 -subjid xxxx or recon-all -autorecon-pial *autorecon2-cp* *autorecon2-wm* -autorecon3 -subjid xxxx Thanks for your help :) Sincerely, Mark Fletcher

1 0

freeview brush value problem
by Tracy Riggins 21 Jul '14

21 Jul '14

Hello Freesurfer experts I am trying to edit an aseg file. I want to set my brush value to a value that does not exist in the default lookup table (e.g,. 570). So, I edited the lookup table, saved this , and loaded it in freeview. My new label and brush value is now visible. However when I set the brush value to this value/label (both the main GUI window and edit voxel window say "Brush value = 570") and draw on the aseg file it does not use this brush value. It uses a different value (eg 48). I can't seem to figure out why. Any thoughts? Thanks in advance, Tracy Tracy Riggins, Ph.D. Assistant Professor Department of Psychology BPS 2147J University of Maryland College Park, MD 20742 Office: (301) 405-5905 Fax: (301) 314-9566 http://ncdl.umd.edu<http://ncdl.umd.edu/> riggins(a)umd.edu<mailto:riggins@umd.edu>

2 1

Tools for FreeSurfer group analysis in R
by Andreas Berger 21 Jul '14

21 Jul '14

Hello FreeSurfers, during the course of my diploma thesis i wrote a few R functions around FreeSurfer, mainly for group analysis of FreeSurfer surface data in R. I thought i'd just contribute the code here, on the off chance that someone may find it useful, and despite the code quality (i wrote this as part of me learning R, so it may actually set your house on fire and so on). It could be of use for someone already working in R and looking to bring FreeSurfer data into it, or someone looking for an alternative to QDEC and the plots created by QDEC. Hardware requirement is not entirely modest unfortunately; it eats up about 8G of memory and takes about 15 minutes around here. Here's an example of how it is used: ############################ BEGIN R CODE ############################## # Presupposing that you have a dataframe 'data' containing your group # description and that all subjects have been processed in ~/subjects ### Group analysis with lm ### # Run the model. This takes some time and memory. lm <- freesurfer.group.analysis.lm(data = data, formula = thickness ~ Age + Gender + Diagnosis, hemi = "lh", fwhm = 10, template = "fsaverage", freesurfer_subjects_dir = "~/subjects") # Extract t values from the model, do FDR correction t <- freesurfer.group.analysis.lm.stats(lm, type.of.value = "t", correction="fdr") # Write the t values to mgh files and open them in FreeSurfer's freeview freesurfer.group.analysis.lm.stats.to.mgh(t, hemi = "lh", template = "fsaverage", freesurfer_subjects_dir = "~/subjects", launch.freeview = TRUE) # Write the t values to .dset ASCII files and open them in AFNI's SUMA freesurfer.group.analysis.lm.stats.to.suma(t, hemi = "lh", template = "fsaverage", freesurfer_subjects_dir = "~/subjects", launch.suma = TRUE) # Make nice plots for vertex #1000 freesurfer.group.analysis.lm.plot(lm, vertex = 1000, interactive = TRUE) ### Group analysis with QDEC ### # Call QDEC freesurfer.group.analysis.QDEC(data = data, freesurfer_subjects_dir = "~/subjects") ### Extract stats ### # Get segmentation stats segstats <- freesurfer.get.segstats(subjects = data$ID, freesurfer_subjects_dir = "~/subjects") ############################## END R CODE ############################## i uploaded the R package here, i called it freesurfR: https://www.nitrc.org/projects/freesurf_r/ thanks also goes to Heath Pardoe, whose code for reading and writing MGH files (also posted on this list) i was allowed to include. regards, Andreas

2 1

Why conn't find subject in Segmentation of hippocampal subfields
by wang kangcheng 20 Jul '14

20 Jul '14

Dear experts When I do Segmentation of hippocampal subfields, it report that ERROR:connot find /nfs/preproc/muti_first/s.... What should I do next? THe scripts was following: SUBJECTS_DIR=/nfs/preproc/multi_first subject_list=${SUBJECTS_DIR}/scripts/check/hipp.list ## Get subjects to run subjects=$( cat ${subject_list} ) ## A. SUBJECT LOOP for subject in $subjects do echo -------------------------- echo running subject ${subject} echo -------------------------- recon-all -s ${subject} -hippo-subfields done echo All have finished Wangkangcheng

1 0

R: Re: Yeo parcellation
by stdp82＠virgilio.it 20 Jul '14

20 Jul '14

Hi list and Thomas, I'm running the first step to use Yeo's parcellation. The running is wrong. mri_surf2surf --s Control01 --hemi lh --sval-tal-xyz white --tval lh.white.Control01 --tval-xyz --trgsubject fsaverage5/label/rh.Yeo2011_7NetworksConfidence_N1000.mgz Reading target surface reg $SUBJECTS/fsaverage5/label/rh.Yeo2011_7NetworksConfidence_N1000.mgz/surf/lh.sphere.regMRISread( $SUBJECTS/fsaverage5/label/rh.Yeo2011_7NetworksConfidence_N1000.mgz/surf/lh.sphere.reg): could not open fileNot a directorymri_surf2surf: could not read surface $SUBJECTS/fsaverage5/label/rh.Yeo2011_7NetworksConfidence_N1000.mgz/surf/lh.sphere.reg Thanks, Stefano ----Messaggio originale---- Da: ythomas(a)csail.mit.edu Data: 7-feb-2014 9.49 A: <stdp82(a)virgilio.it> Cc: "freesurfer"<freesurfer(a)nmr.mgh.harvard.edu> Ogg: Re: [Freesurfer] Yeo parcellation Hi Stefano, A) From FreeSurfer 5.2 onwards, the Yeo parcellation files can be found inside $FREESURFER_HOME/subjects/fsaverage/label/, so you can simply use those files there instead of downloading from the wiki. B) You can use mri_surf2surf to transform the annot files from the fsaverage "subject" to the individual subjects' surface. There is an example of how to do it if you run "mri_surf2surf --help" C) Once you have created the subject's specific annot file from part B, I believe you must first run mris_anatomical_stats, and then run aparcstats2table. However, I am not sure about this, so someone on the list should correct me if I am wrong! Regards, Thomas On Fri, Feb 7, 2014 at 3:58 PM, <stdp82(a)virgilio.it> wrote: > Hi list, > > I have download the files of Yeo parcellation from > ftp://surfer.nmr.mgh.harvard.edu/pub/data/Yeo_JNeurophysiol11_MNI152.zip. > > I have two questions, please. > > A- Which is the directory where I should put these files? > B-I have already run by recon-all -all all my subjects. How can I added to > them Yeo parcellation? > C-I'd like also calculate the cortical thickness for each Yeo's ROI. Is > possible with aparcstats2table? > > Thanks. > > Regards, > > > Stefano > > _______________________________________________ > Freesurfer mailing list > Freesurfer(a)nmr.mgh.harvard.edu > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > > > 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. > _______________________________________________ Freesurfer mailing list Freesurfer(a)nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

2 1

Re: [Freesurfer] mri_cor2label revisited
by Douglas N Greve 20 Jul '14

20 Jul '14

Are you trying to view it in the volume or the surface? When I run freeview viewing the surface and load the label using the "Load Label" button (not File->LoadLabel), it works for me. Can you try the tksurfer command below? It runs for me. You can ignore the warning about 0 nonzero vertices since they are all set to 1. tksurfer fsaverage lh inflated -label ~/tmp/lh_binarymap_area_covar8.label subject is fsaverage hemi is lh surface is inflated surfer: current subjects dir: /autofs/cluster/con_009/users/greve/ixi-subjects surfer: not in "scripts" dir ==> using cwd for session root surfer: session root data dir ($session) set to: surfer: /autofs/cluster/con_009/users/greve/ixi-subjects/ixi320/mri checking for nofix files in 'inflated' Reading image info (/autofs/cluster/con_009/users/greve/ixi-subjects/fsaverage) Reading /autofs/cluster/con_009/users/greve/ixi-subjects/fsaverage/mri/orig.mgz surfer: Reading header info from /autofs/cluster/con_009/users/greve/ixi-subjects/fsaverage/mri/orig.mgz surfer: vertices=163842, faces=327680 surfer: curvature read: min=-0.673989 max=0.540227 surfer: single buffered window surfer: tkoInitWindow(fsaverage) setting percentile thresholds (-1.00, 0.00, 0.00) loading label /homes/4/greve/tmp/lh_binarymap_area_covar8.label reading white matter vertex locations... ******************** 0 nonzero vertices found ******************** label stat field identically zero - setting to 1 surfer: using interface /autofs/space/tanha_002/users/greve/fsdev.build/tktools/tksurfer.tcl Reading /autofs/space/tanha_002/users/greve/fsdev.build/tktools/tkm_common.tcl Reading /autofs/space/tanha_002/users/greve/fsdev.build/tktools/tkm_wrappers.tcl Reading /autofs/space/tanha_002/users/greve/fsdev.build/lib/tcl/fsgdfPlot.tcl Reading /autofs/space/tanha_002/users/greve/fsdev.build/tktools/tkUtils.tcl Successfully parsed tksurfer.tcl reading white matter vertex locations... On 07/15/2014 09:04 AM, Lars M. Rimol wrote: > Hi, > > I go to File -> Label -> load label and there I pick > lh_binarymap_area_covar8.label and then I get the error message: > > > % ******************** 0 nonzero vertices found ******************** > label stat field identically zero - setting to 1 > > > That's it! > > I don't know, but it seems to be saying that all vertices are zero, > doesn't it? > > When I cat the file, I see this: > > 163839 -34.195 -22.897 -24.093 0.0000000000 > > 163840 -34.295 -23.062 -23.449 0.0000000000 > > 163841 -35.093 -23.358 -22.788 0.0000000000 > > The leftmost column surely is vertex number, I don't know what the > next two are (if this were a volume, I'd say x,y,z coordinates) and if > the final column is the value of each vertex, I would have to agree > with tksurfer that there are no nonzero vertices. > > When I load the mgh file, thogh, it has some zero vertices and 98877 > that are 100, that is they are larger than 99.999999999999992. > > >> size(find(gt(testing,99.999999999999992))) > > ans = > > 98877 1 > > > Here's the code I used to make the binary map and save it as an mgh file: > > log10pID = significance threshold > > binaryVLBWlhcovar8(find(gt(binaryVLBWlhcovar8,log10pID))) = 100; > binaryVLBWlhcovar8(find(lt(binaryVLBWlhcovar8,log10pID))) = 0; > > fs_save_mgh(binaryVLBWlhcovar8,'lh_binarymap_area_covar8.mgh',mri); > > > So mri_cor2label must be changing some of these vertices from > 99.99999etc. to 0, right? > > > LMR > > > > On Mon, Jul 14, 2014 at 7:33 PM, Douglas N Greve > <greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu>> wrote: > > > It loads ok for me both in 5.1 and 5.3. Can you give more info on > what is going wrong with tksurfer? > > > On 07/13/2014 05:35 AM, Lars M. Rimol wrote: > > Hi Doug, > > I'm having some problems with the label creation again. When I > run mri_cor2label, it seems OK (please see below), but when I > try to open it in tksurfer, it gives me an error message (bottom). > > Would you mind taking a look at this label to see what's wrong > with it? I'm attaching the binary mgh file and the label. NB! > The binary .mgh file is from matlab, where I have set all > suprathreshold vertices to 100 and all others to 0. Of course, > Matlab doesn't really do this, it sets them to somewhere close > to 99.999999999999992. I don't suppose I could use "greater > than 99" in mri_cor2label? > > ######### > creating label: > > [LMR in ~ ] mri_cor2label --i > /home/lmr/ANALYSIS/STUDIER/VMI_Kam_103/OUTDIR/binary_pmaps/lh_binarymap_area_covar8.mgh > --id 100 --l > /home/lmr/ANALYSIS/STUDIER/VMI_Kam_103/OUTDIR/binary_pmaps/lh_binarymap_area_covar8 > --surf fsaverage lh > $Id: mri_cor2label.c,v 1.12 2011/03/02 00:04:14 nicks Exp $ > Loading mri > /home/lmr/ANALYSIS/STUDIER/VMI_Kam_103/OUTDIR/binary_pmaps/lh_binarymap_area_covar8.mgh > Loading /home/lmr/subjects/fsaverage/surf/lh.white > Scanning the volume > Found 98877 label voxels > Writing label file > /home/lmr/ANALYSIS/STUDIER/VMI_Kam_103/OUTDIR/binary_pmaps/lh_binarymap_area_covar8 > Centroid: -28.93 -15.97 10.14 > mri_cor2label completed SUCCESSFULLY > > ####### > loading it in tksurfer: > > % ******************** 0 nonzero vertices found > ******************** > label stat field identically zero - setting to 1 > > > > Thank you! > > -- > yours, > > Lars M. Rimol, PhD > St. Olavs Hospital > Trondheim, > Norway > > > -- > Douglas N. Greve, Ph.D. > MGH-NMR Center > greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu> > Phone Number: 617-724-2358 <tel:617-724-2358> > Fax: 617-726-7422 <tel:617-726-7422> > > Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting > <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> > FileDrop: https://gate.nmr.mgh.harvard.edu/filedrop2 > www.nmr.mgh.harvard.edu/facility/filedrop/index.html > <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> > Outgoing: > ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/ > > > > 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. > > > > > -- > yours, > > Lars M. Rimol, PhD > St. Olavs Hospital > Trondheim, > Norway -- 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/

2 1

Group analysis tutorial
by Vilde 19 Jul '14

19 Jul '14

Dear mailing list, I am doing the group analysis freesurfer tutorial. When I try to open the significance map, I get the error message "unrecognized sub-opton flag 'overlay_threshold'. " The command is freeview -f $SUBJECTS_DIR/fsaverage/surf/lh.inflated:annot=aparc.annot:annot_outline=1:overlay=lh.gender_age.glmdir/lh-Avg-thickness-age-Cor/sig.mgh:overlay_threshold=4,5 -viewport 3d I get the same error message doing the tutorial on my mac os x and on a linux, so I doubt there is a problem with the tutorial data. All files needed seem to be in the directory. I very much appreciate any help Best, Vilde

2 1

Tracula AD weighted average
by Cat Chong 18 Jul '14

18 Jul '14

Dear Anastasia, Another quick question: We found significant differences in two subject groups for the AD_AVG as well as the AD_AVG_Weight, which is stated as the "weighted average over the entire path distribution'. Could you tell me how the average weight is calculated and why it would be better to use than the AD_AVG perhaps? Thank you very much, Cat

2 1

Re: [Freesurfer] recon-all ERROR, no idea why
by Remy Wahnoun 18 Jul '14

18 Jul '14

Thanks Douglas, I definitely appreciate your help. Here is what it Spit out, *Rémy Wahnoun, PhD* Biomedical Research Group Manager Barrow Neurological Institute at Phoenix Children's Hospital Adjunct Faculty, Arizona State University *www.NeuralControl.net <http://www.neuralcontrol.net/>* set echo = 1 ; set debug = 1 ; set debug = 1 breaksw breaksw end end while ( $#argv != 0 ) while ( 11 != 0 ) set flag = $argv[1] ; shift ; set flag = --n shift switch ( $flag ) switch ( --n ) if ( $#argv == 0 ) goto arg1err ; if ( 10 == 0 ) goto arg1err set nIters = $argv[1] ; shift ; set nIters = 1 shift breaksw breaksw end end while ( $#argv != 0 ) while ( 9 != 0 ) set flag = $argv[1] ; shift ; set flag = --proto-iters shift switch ( $flag ) switch ( --proto-iters ) if ( $#argv == 0 ) goto arg1err ; if ( 8 == 0 ) goto arg1err set nProtoIters = $argv[1] ; shift ; set nProtoIters = 1000 shift breaksw breaksw end end while ( $#argv != 0 ) while ( 7 != 0 ) set flag = $argv[1] ; shift ; set flag = --distance shift switch ( $flag ) switch ( --distance ) if ( $#argv == 0 ) goto arg1err ; if ( 6 == 0 ) goto arg1err set Distance = $argv[1] ; shift ; set Distance = 50 shift breaksw breaksw end end while ( $#argv != 0 ) while ( 5 != 0 ) set flag = $argv[1] ; shift ; set flag = --no-rescale shift switch ( $flag ) switch ( --no-rescale ) set DoRescale = 0 ; set DoRescale = 0 breaksw breaksw end end while ( $#argv != 0 ) while ( 4 != 0 ) set flag = $argv[1] ; shift ; set flag = --i shift switch ( $flag ) switch ( --i ) if ( $#argv == 0 ) goto arg1err ; if ( 3 == 0 ) goto arg1err set InVol = $argv[1] ; shift ; set InVol = orig.mgz shift breaksw breaksw end end while ( $#argv != 0 ) while ( 2 != 0 ) set flag = $argv[1] ; shift ; set flag = --o shift switch ( $flag ) switch ( --o ) if ( $#argv == 0 ) goto arg1err ; if ( 1 == 0 ) goto arg1err set OutVol = $argv[1] ; shift ; set OutVol = orig_nu.mgz shift breaksw breaksw end end while ( $#argv != 0 ) while ( 0 != 0 ) goto parse_args_return ; goto parse_args_return goto check_params ; goto check_params if ( $#InVol == 0 ) then if ( 1 == 0 ) then if ( ! -e $InVol ) then if ( ! -e orig.mgz ) then if ( $#OutVol == 0 ) then if ( 1 == 0 ) then goto check_params_return ; goto check_params_return set OutDir = `dirname $OutVol` ; set OutDir = `dirname $OutVol` dirname orig_nu.mgz mkdir -p $OutDir ; mkdir -p . if ( $#LF == 0 ) set LF = $OutDir/mri_nu_correct.mni.log if ( 0 == 0 ) set LF = ./mri_nu_correct.mni.log set LF = ./mri_nu_correct.mni.log if ( -e $LF ) mv $LF $LF.bak if ( -e ./mri_nu_correct.mni.log ) mv ./mri_nu_correct.mni.log ./mri_nu_correct.mni.log.bak mv ./mri_nu_correct.mni.log ./mri_nu_correct.mni.log.bak pwd | tee -a $LF pwd tee -a ./mri_nu_correct.mni.log /media/sf_Imagery/FreeSurfer/xx/mri which mri_nu_correct.mni | tee -a $LF which mri_nu_correct.mni tee -a ./mri_nu_correct.mni.log /usr/local/freesurfer/bin/mri_nu_correct.mni echo $cmdargs | tee -a $LF echo --debug --n 1 --proto-iters 1000 --distance 50 --no-rescale --i orig.mgz --o orig_nu.mgz tee -a ./mri_nu_correct.mni.log --debug --n 1 --proto-iters 1000 --distance 50 --no-rescale --i orig.mgz --o orig_nu.mgz echo "nIters $nIters" | tee -a $LF echo nIters 1 tee -a ./mri_nu_correct.mni.log nIters 1 echo $VERSION | tee -a $LF tee -a ./mri_nu_correct.mni.log echo $Id: mri_nu_correct.mni,v 1.18.2.1 2013/01/09 21:23:42 nicks Exp $ $Id: mri_nu_correct.mni,v 1.18.2.1 2013/01/09 21:23:42 nicks Exp $ uname -a | tee -a $LF uname -a tee -a ./mri_nu_correct.mni.log Linux xubuntu-VirtualBox 3.2.0-64-generic #97-Ubuntu SMP Wed Jun 4 22:03:48 UTC 2014 i686 i686 i386 GNU/Linux date | tee -a $LF date tee -a ./mri_nu_correct.mni.log Mon Jun 23 15:32:45 EDT 2014 nu_correct -version | tee -a $LF nu_correct -version tee -a ./mri_nu_correct.mni.log Program nu_correct, built from: Package MNI N3, version 1.10, compiled by nicks@fishie (i686-pc-linux-gnu) on 2010-02-20 at 20:19:13 if ( ! $#tmpdir ) set tmpdir = $OutDir/tmp.mri_nu_correct.mni.$$ if ( ! 0 ) set tmpdir = ./tmp.mri_nu_correct.mni.1745 set tmpdir = ./tmp.mri_nu_correct.mni.1745 mkdir -p $tmpdir mkdir -p ./tmp.mri_nu_correct.mni.1745 echo "tmpdir is $tmpdir" | tee -a $LF echo tmpdir is ./tmp.mri_nu_correct.mni.1745 tee -a ./mri_nu_correct.mni.log tmpdir is ./tmp.mri_nu_correct.mni.1745 set cmd = ( mri_convert ${HiRes} $InVol $tmpdir/nu0.mnc ) set cmd = ( mri_convert orig.mgz ./tmp.mri_nu_correct.mni.1745/nu0.mnc ) if ( $UseFloat ) set cmd = ( $cmd -odt float ) if ( 1 ) set cmd = ( mri_convert orig.mgz ./tmp.mri_nu_correct.mni.1745/nu0.mnc -odt float ) set cmd = ( mri_convert orig.mgz ./tmp.mri_nu_correct.mni.1745/nu0.mnc -odt float ) pwd | & tee -a $LF pwd tee -a ./mri_nu_correct.mni.log /media/sf_Imagery/FreeSurfer/xx/mri echo $cmd | & tee -a $LF echo mri_convert orig.mgz ./tmp.mri_nu_correct.mni.1745/nu0.mnc -odt float tee -a ./mri_nu_correct.mni.log mri_convert orig.mgz ./tmp.mri_nu_correct.mni.1745/nu0.mnc -odt float $cmd | & tee -a $LF mri_convert orig.mgz ./tmp.mri_nu_correct.mni.1745/nu0.mnc -odt float tee -a ./mri_nu_correct.mni.log mri_convert orig.mgz ./tmp.mri_nu_correct.mni.1745/nu0.mnc -odt float $Id: mri_convert.c,v 1.179.2.7 2012/09/05 21:55:16 mreuter Exp $ reading from orig.mgz... TR=25.00, TE=5.36, TI=0.00, flip angle=30.00 i_ras = (-1, 1.49012e-08, 0) j_ras = (0, 0, -1) k_ras = (-1.30385e-08, 1, -5.96046e-08) changing data type from uchar to float (noscale = 0)... writing to ./tmp.mri_nu_correct.mni.1745/nu0.mnc... if ( $status ) then if ( 0 ) then @ nthIter = 1 ; @ nthIter = 1 while ( $nthIter < = $nIters ) while ( 1 < = 1 ) echo " " | & tee -a $LF echo tee -a ./mri_nu_correct.mni.log echo "--------------------------------------------------------" | & tee -a $LF echo -------------------------------------------------------- tee -a ./mri_nu_correct.mni.log -------------------------------------------------------- echo "Iteration $nthIter `date`" | & tee -a $LF echo Iteration 1 `date` tee -a ./mri_nu_correct.mni.log date Iteration 1 Mon Jun 23 15:32:49 EDT 2014 @ m = $nthIter - 1 @ m = 1 - 1 set cmd = ( nu_correct -clobber $tmpdir/nu${m}.mnc $tmpdir/nu${nthIter}.mnc ) set cmd = ( nu_correct -clobber ./tmp.mri_nu_correct.mni.1745/nu0.mnc ./tmp.mri_nu_correct.mni.1745/nu1.mnc ) set cmd = ( $cmd -tmpdir ${tmpdir}/${m}/ ) ; set cmd = ( nu_correct -clobber ./tmp.mri_nu_correct.mni.1745/nu0.mnc ./tmp.mri_nu_correct.mni.1745/nu1.mnc -tmpdir ./tmp.mri_nu_correct.mni.1745/0/ ) if ( $#nProtoIters ) set cmd = ( $cmd -iterations $nProtoIters ) ; if ( 1 ) set cmd = ( nu_correct -clobber ./tmp.mri_nu_correct.mni.1745/nu0.mnc ./tmp.mri_nu_correct.mni.1745/nu1.mnc -tmpdir ./tmp.mri_nu_correct.mni.1745/0/ -iterations 1000 ) set cmd = ( nu_correct -clobber ./tmp.mri_nu_correct.mni.1745/nu0.mnc ./tmp.mri_nu_correct.mni.1745/nu1.mnc -tmpdir ./tmp.mri_nu_correct.mni.1745/0/ -iterations 1000 ) if ( $#StopThresh ) set cmd = ( $cmd -stop $StopThresh ) ; if ( 0 ) set cmd = ( nu_correct -clobber ./tmp.mri_nu_correct.mni.1745/nu0.mnc ./tmp.mri_nu_correct.mni.1745/nu1.mnc -tmpdir ./tmp.mri_nu_correct.mni.1745/0/ -iterations 1000 -stop ) if ( $#Distance ) set cmd = ( $cmd -distance $Distance ) ; if ( 1 ) set cmd = ( nu_correct -clobber ./tmp.mri_nu_correct.mni.1745/nu0.mnc ./tmp.mri_nu_correct.mni.1745/nu1.mnc -tmpdir ./tmp.mri_nu_correct.mni.1745/0/ -iterations 1000 -distance 50 ) set cmd = ( nu_correct -clobber ./tmp.mri_nu_correct.mni.1745/nu0.mnc ./tmp.mri_nu_correct.mni.1745/nu1.mnc -tmpdir ./tmp.mri_nu_correct.mni.1745/0/ -iterations 1000 -distance 50 ) if ( $#FWHM ) set cmd = ( $cmd -fwhm $FWHM ) ; if ( 0 ) set cmd = ( nu_correct -clobber ./tmp.mri_nu_correct.mni.1745/nu0.mnc ./tmp.mri_nu_correct.mni.1745/nu1.mnc -tmpdir ./tmp.mri_nu_correct.mni.1745/0/ -iterations 1000 -distance 50 -fwhm ) echo $cmd | & tee -a $LF echo nu_correct -clobber ./tmp.mri_nu_correct.mni.1745/nu0.mnc ./tmp.mri_nu_correct.mni.1745/nu1.mnc -tmpdir ./tmp.mri_nu_correct.mni.1745/0/ -iterations 1000 -distance 50 tee -a ./mri_nu_correct.mni.log nu_correct -clobber ./tmp.mri_nu_correct.mni.1745/nu0.mnc ./tmp.mri_nu_correct.mni.1745/nu1.mnc -tmpdir ./tmp.mri_nu_correct.mni.1745/0/ -iterations 1000 -distance 50 $cmd | & tee -a $LF nu_correct -clobber ./tmp.mri_nu_correct.mni.1745/nu0.mnc ./tmp.mri_nu_correct.mni.1745/nu1.mnc -tmpdir ./tmp.mri_nu_correct.mni.1745/0/ -iterations 1000 -distance 50 tee -a ./mri_nu_correct.mni.log [fsuser@xubuntu-VirtualBox:/media/sf_Imagery/FreeSurfer/xx/mri/] [2014-06-23 15:32:49] running: /usr/local/freesurfer/mni/bin/nu_estimate_np_and_em -parzen -log -sharpen 0.15 0.01 -iterations 1000 -stop 0.001 -shrink 4 -auto_mask -nonotify -b_spline 1 -distance 50 -quiet -execute -clobber -nokeeptmp -tmpdir ./tmp.mri_nu_correct.mni.1745/0/ ./tmp.mri_nu_correct.mni.1745/nu0.mnc ./tmp.mri_nu_correct.mni.1745/nu1.imp Assertion failed at line 827 in file templates/CachedArray.cc nu_estimate_np_and_em: crashed while running volume_stats (termination status=256) nu_correct: crashed while running nu_estimate_np_and_em (termination status=256) if ( $status ) then if ( 0 ) then rm -f $tmpdir/nu${m}.mnc $tmpdir/nu${m}.imp rm -f ./tmp.mri_nu_correct.mni.1745/nu0.mnc ./tmp.mri_nu_correct.mni.1745/nu0.imp @ nthIter = $nthIter + 1 ; @ nthIter = 1 + 1 echo " " | & tee -a $LF echo tee -a ./mri_nu_correct.mni.log end end while ( $nthIter < = $nIters ) while ( 2 < = 1 ) echo " " | & tee -a $LF echo tee -a ./mri_nu_correct.mni.log echo " " | & tee -a $LF echo tee -a ./mri_nu_correct.mni.log set numnc = $tmpdir/nu$nIters.mnc set numnc = ./tmp.mri_nu_correct.mni.1745/nu1.mnc if ( $DoRescale ) then if ( 0 ) then if ( ! -e $numnc ) then if ( ! -e ./tmp.mri_nu_correct.mni.1745/nu1.mnc ) then echo "ERROR: file $numnc does not exist!" | & tee -a $LF echo ERROR: file ./tmp.mri_nu_correct.mni.1745/nu1.mnc does not exist! tee -a ./mri_nu_correct.mni.log ERROR: file ./tmp.mri_nu_correct.mni.1745/nu1.mnc does not exist! exit 1 ; exit 1 ------------------------------ This transmission, including any attachments, is for the sole use of the intended recipient (s) and may contain information that is confidential, proprietary, legally privileged, or otherwise protected by law from disclosure. Any unauthorized review, use, copying, disclosure, or distribution is prohibited. If you are not the intended recipient, or the person responsible for delivering this to an addressee, you should notify the sender immediately by telephone or by reply e-mail, and destroy all copies of the original message.

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