Freesurfer March 2004

freesurfer@nmr.mgh.harvard.edu

21 participants
23 discussions

Looking for correct coordinates
by Daniel Goldenholz 31 Mar '04

31 Mar '04

Hi freesurfer folks. I have a pair of inflated brains that I'd like to work with in matlab. I found a nice function called "read_surf.m" which works very well. But I noticed that the two hemispheres cannot be viewed at the same time, because their locations overlap! I pushed one to the left and one to the right by an arbitrary amount, and this looked a bit better. But I was hoping to put the position of the inflated brain in the correct spot so that I could correlate MEG info with these surface files. What else do I need to know about positioning these brain surfaces? Thanks -- Daniel Goldenholz -----------------------------------------------------

2 1

remote desktop turns blue?
by Nguyen, Thang Q 31 Mar '04

31 Mar '04

Hello all, I don't know if anyone has had this problem before but when I remote desktop into a PC running freesurfer (csurf), the brain image in tkmedit turns blue. In addition, tksurfer doesn't load. When I relog on to the PC, I can reload tkmedit and it works fine. Here are some specifications: Running freesurfer-20030724 Remote desktop with windows XP Pro to another windows XP Pro PC Running Reflections X on the PC desktop Running RedHat v9.0 on the compute nodes In addition to all this trouble, after editing in tkmedit and creating a new surface, the surface doesn't reflect the new changes. This hasn't happened before. Does anybody have any ideas on what contributes to these problems? Thanks, Thang ------------------------ Thang Nguyen Psychiatry and Human Behavior - Brain Imaging Center University of California, Irvine This e-mail/fax 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/fax and destroy all copies of the original message.

3 2

Pial Surface
by Keyma Prince 30 Mar '04

30 Mar '04

When you look at the values that are produced by mris_anatomical_stats for surface area, are these values the same as pial surface area? Is the pial surface the same as the inflated surface? -Keyma

2 5

having problem with motion correction.
by Tamara Knutsen 26 Mar '04

26 Mar '04

error is: >> >> >> >> >> >> ??? Undefined function or variable 'afni2bvol'. >> ------------------------------------------ ERROR: creating 011/fmc i'm just running: mc-sess -sf sf -df df at the command line (the files are correct, and the motion correction seems to start fine, but then it can't create these fmc files. any ideas? thanks tamara

1 0

Snapshot of version
by Jonathan Wisco 25 Mar '04

25 Mar '04

Hi, I was wondering if I someone could send me a snapshot or pointer to the freesurfer recon tools for cortical inflation that was being used between Sept 2003 and Mar 2004. I suspect that the version changed somewhere in the middle of this time period and would like to analyze my data using the same version of freesurfer. Thanks, Jon --- Jonathan J. Wisco, Ph.D. Massachusetts General Hospital-NMR Center Building 149, 13th St., Mailcode 149-2301 Charlestown, MA 02129-2060 e-mail: jjwisco(a)nmr.mgh.harvard.edu tel: 617-851-8492 fax: 617-726-7422

6 9

Re: [Freesurfer] is there any more documentation for mri_segment?
by Tamara Knutsen 25 Mar '04

25 Mar '04

you have to switch environments type use_dev (brings you to nmr-std-dev environment) and then type mri_segment -u and you will get the appropriate output. yes...documentation in general for all these things are indeed lacking. On Thu, 25 Mar 2004, Ned T. Sahin wrote: > Hi All, > > This thread brings up a more general point. It starts from the idea that > probably (just my guess) many people were surprised to find out that the > "-u" argument would yield more information about mri_segment. > > Let me make the strong assertion that well-behaved code should spit out > in the short version of feedback/documentation the method for getting > the long version. As in, when you just type in mri_segment, you should > get back a simple usage file that includes the flag you would need to > supply in order to get more usage information. > > How would a user know to try -u? It is only listed in the output you get > when you supply the flag. > > In fact, I just tried it: the -u flag does not give the output Brian > posted when most users give it. It only works if you have access to the > development environment. > > Furthermore, this is in the context of there being many standards for > how to evince extra online help from the various utilities used as part > of FreeSurfer and FS-FAST. > > In some cases, the flag is "-help" in others it is "--help" and in > others, apparently, it is "-u". When it is one, it is not the other - > for instance neither 'help' flag works with mri_segment. > > Of course I realize that much of this is historical - different > standards at different times, when various utilities were written. > Perhaps one is the chosen standard now. > > Any one person having to check and re-work all ~832 utilities, well, > would be quite annoyed. Especially because it is a mindless task. But, > could we utilize how mindless it is, and pawn it off in chunks to do in > spare moments to a fairly large net of developers at the Center? 5 > programs a week times 10 programmers would mean roughly 1 minute per > day, and yet the task would be done this summer. And of course it could > be much faster than that. The task, again, would be standardizing all > programs to give full info when one flag (say "-help") is supplied, and > making the simple usage prompt reflect this fact. > > The payoff is all users would be able to get the information below > immediately, rather than after 5 emails to the whole list, and in many > cases would be able to answer their own questions. > > Sorry for the rant, here. I hope it was useful, and the suggestion can > be implemented. > > > Thanks! > > > -ned t. sahin > > > > Brian T. Quinn wrote: > > > try the -u option for the following: > > > > -slope <float s> set the curvature slope (both n and p) > > -pslope <float p> set the curvature pslope (default=1.0) > > -nslope <float n> set the curvature nslope (default=1.0) > > -debug_voxel <int x y z> set voxel for debugging > > -auto automatically detect class statistics (default) > > -noauto don't automatically detect class statistics > > -log log to ./segment.dat > > -keep keep wm edits. maintains all values of 0 and 255 > > -ghi, -gray_hi <int h> set the gray matter high limit > > (default=100.000) > > -wlo, -wm_low <int l> set the white matter low limit > > (default=90.000) > > -whi, -wm_hi <int h> set the white matter high limit > > (default=125.000) > > -nseg <int n> thicken the n largest thin strands (default=20) > > -thicken toggle thickening step (default=ON) > > -fillbg toggle filling of the basal ganglia > > (default=OFF) > > -fillv toggle filling of the ventricles (default=OFF) > > -b <float s> set blur sigma (default=0.25) > > -n <int i> set # iterations of border classification > > (default=1) > > -t <int t> set limit to thin strands in mm (default=4) > > -v verbose > > -p <float p> set % threshold (default=0.80) > > -x <filename> extract options from filename > > -w <int w> set wsize (default=11) > > -u usage > > > > good luck, > > -- > > brian t. quinn > > > > > > On Thu, 25 Mar 2004, Tamara Knutsen wrote: > > > > > >>this is all i find: > >> > >>[whatever:ben1] (nmr-std-env) mri_segment > >>usage: mri_segment <input volume> <output volume> > >> > >>thanks! > >> > >> > >>_______________________________________________ > >>Freesurfer mailing list > >>Freesurfer(a)nmr.mgh.harvard.edu > >>http://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > >> > > > > > > _______________________________________________ > > Freesurfer mailing list > > Freesurfer(a)nmr.mgh.harvard.edu > > http://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > >

1 0

stage1 error
by Lukas Pezawas 24 Mar '04

24 Mar '04

Hi everybody, I've been having a strange problem running stage1 of recon-all. The process finishes successfully when I select "Process Volume" from the GUI, but fails when I start the stage via terminal command (using the same brain). I would prefer using the terminal command because I've been told it is a more up-to-date version than the GUI command. The command I entered was "recon-all -stage1 -nomorphrgb -subjid test0317bTerm" and it appeared to terminate during intensity normalization. Below is the error log. Thanks. -Aaron- on Mar 22 15:09:11 EST 2004 /usr/local/freesurfer_dev/bin/noarch/recon-all -stage1 -nomorphrgb -subjid test0317bTerm /j/freesurfsubj/test0317bTerm subjid test0317bTerm SUBJECTS_DIR /j/freesurfsubj FREESURFER_HOME /usr/local/freesurfer_dev Linux Filesystem 1K-blocks Used Available Use% Mounted on freud:/export/freesurfsubj 74697640 36308304 34594784 52% /j/freesurfsubj --------------------------------------------- MotionCor Mon Mar 22 15:09:12 EST 2004 INFO: no runs were found in /j/freesurfsubj/test0317bTerm/mri/orig/XXX, however, there is a volume in /j/freesurfsubj/test0317bTerm/mri/orig, so I am just going to proceed with that as the input. --------------------------------------------- Talairach Mon Mar 22 15:09:14 EST 2004 /j/freesurfsubj/test0317bTerm/mri talairach2 test0317bTerm INFO: registering test0317bTerm Mon Mar 22 15:09:14 EST 2004 -------------------------------------------- mri_convert orig orig/orig.mnc mri_convert orig orig/orig.mnc reading from orig... TR=0.00, TE=0.00, TI=0.00, flip angle=0.00 i_ras = (-1, 0, 0) j_ras = (0, 0, -1) k_ras = (0, 1, 0) writing to orig/orig.mnc... -------------------------------------------- mritotal -verbose -debug -clobber -protocol icbm orig/orig.mnc transforms/talairach.xfm -------------------------------------------- mri_add_xform_to_header /j/freesurfsubj/test0317bTerm/mri/transforms/talairach.xfm /j/freesurfsubj/test0317bTerm/mri/orig /j/freesurfsubj/test0317bTerm/mri/orig Mon Mar 22 15:11:14 EST 2004 --------------------------------------------- Nu Intensity Correction Mon Mar 22 15:11:14 EST 2004 mri_convert orig nu/nu0.mnc mri_convert orig nu/nu0.mnc reading from orig... TR=0.00, TE=0.00, TI=0.00, flip angle=0.00 i_ras = (-1, 0, 0) j_ras = (0, 0, -1) k_ras = (0, 1, 0) writing to nu/nu0.mnc... nu_correct -clobber nu/nu0.mnc nu/nu1.mnc [goldmana@feynman:/j/freesurfsubj/test0317bTerm/mri/] [2004-03-22 15:11:32] running: /usr/local/bin/nu_estimate_np_and_em -parzen -log -sharpen 0.15 0.01 -iterations 50 -stop 0.001 -shrink 4 -auto_mask -nonotify -b_spline 1 -distance 200 -quiet -execute -clobber -nokeeptmp -tmpdir /usr/tmp/nu_correct_21409/ nu/nu0.mnc nu/nu1.imp Processing:............................................................. .... Done Not implemented yet in cache_volume_range_has_changed() Not implemented yet in cache_volume_range_has_changed() Number of iterations: 1 CV of field change: 2.16148e-26 [goldmana@feynman:/j/freesurfsubj/test0317bTerm/mri/] [2004-03-22 15:11:45] running: /usr/local/bin/make_template -quiet -shrink 3 nu/nu0.mnc /usr/tmp/nu_correct_21409//template.mnc Transforming slices:................................................................. .... .................Done 256+0 records in 256+0 records out Not implemented yet in cache_volume_range_has_changed() Not implemented yet in cache_volume_range_has_changed() Transforming slices:................................................................. .... ........................................................................ .... ........................................................................ .... ...................................Done nu_correct -clobber nu/nu1.mnc nu/nu2.mnc [goldmana@feynman:/j/freesurfsubj/test0317bTerm/mri/] [2004-03-22 15:12:16] running: /usr/local/bin/nu_estimate_np_and_em -parzen -log -sharpen 0.15 0.01 -iterations 50 -stop 0.001 -shrink 4 -auto_mask -nonotify -b_spline 1 -distance 200 -quiet -execute -clobber -nokeeptmp -tmpdir /usr/tmp/nu_correct_21464/ nu/nu1.mnc nu/nu2.imp Processing:............................................................. .... Done Not implemented yet in cache_volume_range_has_changed() Not implemented yet in cache_volume_range_has_changed() Number of iterations: 1 CV of field change: 3.04238e-26 [goldmana@feynman:/j/freesurfsubj/test0317bTerm/mri/] [2004-03-22 15:12:28] running: /usr/local/bin/make_template -quiet -shrink 3 nu/nu1.mnc /usr/tmp/nu_correct_21464//template.mnc Transforming slices:................................................................. .... .................Done 256+0 records in 256+0 records out Not implemented yet in cache_volume_range_has_changed() Not implemented yet in cache_volume_range_has_changed() Transforming slices:................................................................. .... ........................................................................ .... ........................................................................ .... ...................................Done nu_correct -clobber nu/nu2.mnc nu/nu3.mnc [goldmana@feynman:/j/freesurfsubj/test0317bTerm/mri/] [2004-03-22 15:12:59] running: /usr/local/bin/nu_estimate_np_and_em -parzen -log -sharpen 0.15 0.01 -iterations 50 -stop 0.001 -shrink 4 -auto_mask -nonotify -b_spline 1 -distance 200 -quiet -execute -clobber -nokeeptmp -tmpdir /usr/tmp/nu_correct_21517/ nu/nu2.mnc nu/nu3.imp Processing:............................................................. .... Done Not implemented yet in cache_volume_range_has_changed() Not implemented yet in cache_volume_range_has_changed() Number of iterations: 1 CV of field change: 2.0113e-26 [goldmana@feynman:/j/freesurfsubj/test0317bTerm/mri/] [2004-03-22 15:13:11] running: /usr/local/bin/make_template -quiet -shrink 3 nu/nu2.mnc /usr/tmp/nu_correct_21517//template.mnc Transforming slices:................................................................. .... .................Done 256+0 records in 256+0 records out Not implemented yet in cache_volume_range_has_changed() Not implemented yet in cache_volume_range_has_changed() Transforming slices:................................................................. .... ........................................................................ .... ........................................................................ .... ...................................Done nu_correct -clobber nu/nu3.mnc nu/nu4.mnc [goldmana@feynman:/j/freesurfsubj/test0317bTerm/mri/] [2004-03-22 15:13:43] running: /usr/local/bin/nu_estimate_np_and_em -parzen -log -sharpen 0.15 0.01 -iterations 50 -stop 0.001 -shrink 4 -auto_mask -nonotify -b_spline 1 -distance 200 -quiet -execute -clobber -nokeeptmp -tmpdir /usr/tmp/nu_correct_21575/ nu/nu3.mnc nu/nu4.imp Processing:............................................................. .... Done Not implemented yet in cache_volume_range_has_changed() Not implemented yet in cache_volume_range_has_changed() Number of iterations: 1 CV of field change: 2.08277e-26 [goldmana@feynman:/j/freesurfsubj/test0317bTerm/mri/] [2004-03-22 15:13:57] running: /usr/local/bin/make_template -quiet -shrink 3 nu/nu3.mnc /usr/tmp/nu_correct_21575//template.mnc Transforming slices:................................................................. .... .................Done 256+0 records in 256+0 records out Not implemented yet in cache_volume_range_has_changed() Not implemented yet in cache_volume_range_has_changed() Transforming slices:................................................................. .... ........................................................................ .... ........................................................................ .... ...................................Done mri_convert nu/nu4.mnc nu mri_convert nu/nu4.mnc nu reading from nu/nu4.mnc... DimSizes: 3, 256, 256, 256, 1 DataType: 1 MINC Linear Transform -1.0000 0.0000 0.0000 123.0000 0.0000 0.0000 1.0000 -122.3080 0.0000 -1.0000 0.0000 139.8577 0.0000 0.0000 0.0000 1.0000 Done reading minc TR=0.00, TE=0.00, TI=0.00, flip angle=0.00 i_ras = (-1, 0, 0) j_ras = (0, 0, -1) k_ras = (0, 1, 0) Original Data has (1, 1, 1) mm size and (256, 256, 256) voxels. Data is conformed to 1 mm size and 256 voxels for all directions writing to nu... --------------------------------------------- Intensity Normalization Mon Mar 22 15:14:37 EST 2004 /j/freesurfsubj/test0317bTerm/mri mri_normalize orig T1 reading from orig... normalizing image... MRInormalize: could not find any valid peaks No such file or directory mri_normalize: normalization failed No such file or directory ERROR: mri_normalize exited with non-zero status ------------------------------------------ Filesystem 1K-blocks Used Available Use% Mounted on freud:/export/freesurfsubj 74697640 36325824 34577264 52% /j/freesurfsubj ------------------------------------------ recon-all exited with errors at Mon Mar 22 15:14:50 EST 2004

1 0

create surface
by Bob Thoma 24 Mar '04

24 Mar '04

Hello all, I've had a problem with a few subjects. When I run Create Surface it fails at fill white matter. In the log file, it states, could not find corpus callosum. Then an alert states /home dir/<file name>/mri/filled/COR-001 .... not found ... wmfill first.. when you say ok, it says process interrupted, and closes. Anyone have any idea as to why it is failing, and what I can do to fix it?? Thanks Natalie _________________________________________________________________ Check out MSN PC Safety & Security to help ensure your PC is protected and safe. http://specials.msn.com/msn/security.asp

1 0

Map_all_labels
by Keyma Prince 24 Mar '04

24 Mar '04

When map_all_labels is run is it using the average7 brain to know where to map the labels onto the surfaces? Also why are the labels (such as the central sulcus) generated by map_all_labels and parcellate_subject different? -Keyma Keyma Prince Massachusetts Institute of Technology Department of Brain and Cognitive Sciences NE20-392 Cambridge, MA 02139

2 1

Detecting errors
by Satrajit Ghosh 22 Mar '04

22 Mar '04

Hi, Here is the output of lh.euler.orig for a particular subject: euler # = v-e+f = 2g-2: 149150 - 447444 + 298297 = 3 --> 0 holes F =2V-4: 298297 != 298300-4 (-1) 2E=3F: 894888 != 894891 (-3) total defect index = 4 Is there a way of determining which parts of the surface and therefore which parts of the volume are giving rising to the defects? It's likely I can write some surface analysis code to do this, but I'm hoping freesurfer already has something of this sort available. Thanks, Satra -- Satrajit S. Ghosh Department of Cognitive and Neural Systems, Boston University (I): http://www.cns.bu.edu/~satra

3 3

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

2006

2005

2004

Freesurfer March 2004