Freesurfer

freesurfer@nmr.mgh.harvard.edu

30132 discussions

corpus callosum volume in aseg.mgz...
by Gonzalo Rojas Costa 16 Mar '11

16 Mar '11

Hi: The CC volume in aseg.mgz is computed using two slices to the left and two sagital slices to the right the midline ?... Sincerely, -- Gonzalo Rojas Costa Laboratory for Advanced Medical Image Processing Department of Radiology Las Condes Clinic Lo Fontecilla 441, Las Condes, Santiago, Chile. Tel: 56-2-2105170 Cel: 56-9-97771785 www.clc.cl

2 1

MCSRCH error
by Juranek, Jenifer 16 Mar '11

16 Mar '11

While processing a normal brain through fsv5.0.0, an unfamiliar message popped up in the screen output (see below). Any ideas what I should do about this message (anything in particular I should check)? Many Thanks, Jenifer IFLAG= -1 LINE SEARCH FAILED. SEE DOCUMENTATION OF ROUTINE MCSRCH ERROR RETURN OF LINE SEARCH: INFO= 3 POSSIBLE CAUSES: FUNCTION OR GRADIENT ARE INCORRECT OR INCORRECT TOLERANCESpass 2 through quasi-newton minimization... IFLAG= -1 LINE SEARCH FAILED. SEE DOCUMENTATION OF ROUTINE MCSRCH ERROR RETURN OF LINE SEARCH: INFO= 3 POSSIBLE CAUSES: FUNCTION OR GRADIENT ARE INCORRECT OR INCORRECT TOLERANCESoutof QuasiNewtonEMA: 009: -log(p) = 1145193.0 tol 0.000000 final transform: 1.386 -0.067 0.094 -63.518; 0.088 1.120 -0.108 -0.481; -0.131 0.093 0.972 -2.704; 0.000 0.000 0.000 1.000; writing output transformation to transforms/talairach_with_skull.lta... registration took 34 minutes and 38 seconds. mri_watershed -T1 -brain_atlas /usr/local/freesurfer/average/RB_all_withskull_2008-03-26.gca transforms/talairach_with_skull.lta T1.mgz brainmask.auto.mgz Mode: T1 normalized volume Mode: Use the information of atlas (default parms, --help for details)

2 1

Is there a way to obtain the thickness of every vertex of the whole brain?
by 汪贵宏 16 Mar '11

16 Mar '11

Hi,all: Just as the title says,I wonder if there is a method to obtain the thickness of each vertex in the whole brain,and so,if we can,how do we kown which of the area is the vertex belong to? I am using the freesurfer 5.0. thanks in advance! wang

2 1

Re: [Freesurfer] corpus calosum volume ?...
by Gonzalo Rojas Costa 15 Mar '11

15 Mar '11

Hi Bruce: Thanks you very much for your answer... And, how can I get the volume of the new segmentation of the CC that I got with mri_cc command ?... Sincerely, Gonzalo Rojas Costa El 15-03-2011 21:27, Bruce Fischl escribió: > I mean that if the volume is rotated in any 2d slice it can look as if > only some of the cc is properly labeled, even if the full 3d labeling is > correct > On Tue, 15 Mar 2011, Gonzalo Rojas Costa wrote: > >> Hi Bruce; >> >> I don't understand why you said that in 3D the labeling is usually >> correct... ???... >> >> Sincerely, >> >> >> Gonzalo Rojas Costa >> >> >> >> El 15-03-2011 15:33, Bruce Fischl escribió: >>> Hi Gonzalo >>> >>> mri_cc computes the optimal rotation. If the patient's head is not >>> straight this can make it look like an inconsistent labeling, but in 3d >>> it usually is correct. >>> >>> Bruce >>> >>> On Tue, 15 Mar 2011, Gonzalo Rojas Costa wrote: >>> >>>> Hi: >>>> >>>> The CC_Posterior, CC_Central, etc, volumes in "aseg.mgz" are the >>>> complete CC volume, or it includes only some sagital slices ?... >>>> >>>> How can I meassure the CC volume for some sagital slices only (5 in >>>> total for example) ?... When I tried using "mri_cc" I got an error >>>> (for example) in the case that the patient is slightly rotated >>>> (freesurfer segments only one section of the CC as CC and the other as >>>> WM)... >>>> >>>> Sincerely, >>>> >>>> >>>> >>>> Gonzalo Rojas Costa >>>> Laboratory for Advanced Medical Image Processing >>>> Department of Radiology >>>> Las Condes Clinic >>>> Lo Fontecilla 441, Las Condes, Santiago, Chile. >>>> Tel: 56-2-2105170 >>>> Cel: 56-9-97771785 >>>> www.clc.cl >>>> >>>> >>>> >>> >>> >>> 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.

1 0

Philips DICOM format problem
by leonardo kay 15 Mar '11

15 Mar '11

Hi All I'm new in freesurfer. I have a question regards Philips DICOM format. I tried to convert from dicom to .mgz using mri_convert but i get the following error message: MRIresample(): error inverting matrix; determinant is nan, matrix is: -0.000 -0.000 nan nan; 0.000 0.000 nan nan; 0.000 -0.000 nan nan; 0.000 0.000 0.000 1.000; Have anyone experienced this before? Any help converting from philips dicom format? Best regards Leo.

3 2

Re: [Freesurfer] mri_glmfit allocate memory error
by Douglas N Greve 15 Mar '11

15 Mar '11

So you re-ran mris_preproc and you still get the error? If so, try running mris_preproc --debug <Your Options> |& tee mris_preprog.err.log Send me the mris_preprog.err.log, but don't email it as it will likely be huge. Instead drop it to me a the file drop listed at the end of this email. doug Andrew Mendlowitz wrote: > The mri_prepoc runs fine still but the glmfit still says can not > allocate memory, there is 8gb of ram and very little disk space is > used up, there is 876.1 gb still free. > > On Fri, Mar 11, 2011 at 12:49 PM, Douglas N Greve > <greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu>> wrote: > > If you re-run mris_preproc does the same thing happen? Are you > running out of disk space? > doug > > Andrew Mendlowitz wrote: > > The code for the lh.paired-diff.thickness.mgh file is > mris_preproc --target fsaverage --hemi lh \ > > --meas thickness --out lh.paired-diff.thickness.mgh \ > > --fsgd pairss.fsgd --paired-diff > > > On Thu, Mar 10, 2011 at 10:13 AM, Douglas N Greve > <greve(a)nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu> > <mailto:greve@nmr.mgh.harvard.edu > <mailto:greve@nmr.mgh.harvard.edu>>> wrote: > > There appears to be something wrong with the > lh.paired-diff.thickness.mgh file. How was it created? > doug > > Andrew Mendlowitz wrote: > > Hello All, > I get the following error while trying to run > mri_glmfit. I am > running a 2.93 ghz imac with 8gb of ram. I am also running > freesurfer version 5.0.0 on Snow > Leopard OS version 10.6. > > > /Applications/freesurfer/subjects/qdec$ mri_glmfit > --glmdir > lh.paired-diff --fsgd paired-diff.fsgd --C mean.mtx --y > lh.paired-diff.thickness.mghgdfReadHeader: reading > paired-diff.fsgd > INFO: DeMeanFlag keyword not found, DeMeaning will NOT > be done. > Continuous Variable Means (all subjects) > 0 Age 1.5 > Class Means of each Continuous Variable > 1 Main 1.5000 > INFO: gd2mtx_method is dods > > $Id: mri_glmfit.c,v 1.187.2.1 2010/07/26 15:54:39 greve > Exp $ > cwd /Applications/freesurfer/subjects/qdec > cmdline mri_glmfit --glmdir lh.paired-diff --fsgd > paired-diff.fsgd --C mean.mtx > --y lh.paired-diff.thickness.mgh > sysname Darwin > hostname larry > machine i386 > user amend > FixVertexAreaFlag = 1 > UseMaskWithSmoothing 1 > OneSampleGroupMean 0 > y > /Applications/freesurfer/subjects/qdec/lh.paired-diff.thickness.mgh > logyflag 0 > usedti 0 > FSGD paired-diff.fsgd > glmdir lh.paired-diff > IllCondOK 0 > DoFFx 0 > Creating output directory lh.paired-diff > Loading y from > > /Applications/freesurfer/subjects/qdec/lh.paired-diff.thickness.mgh > mri_glmfit(79576) malloc: *** mmap(size=2563764224) failed > (error code=12) > *** error: can't allocate region > *** set a breakpoint in malloc_error_break to debug > MRIallocIndices: could not allocate 1714511945 elt > index array > Cannot allocate memory > > The mri_info for the file is : > > > mghRead(/Applications/freesurfer/subjects/qdec/fsgd_input.mgh, > -1): could not > open file > amend@larry /Applications/freesurfer/subjects/qdec$ > mri_info > > /Applications/freesurfer/subjects/qdec/lh.paired-diff.thickness.mgh > Volume information for > > /Applications/freesurfer/subjects/qdec/lh.paired-diff.thickness.mgh > type: MGH > dimensions: 1714510945 x 1853057372 x 1634628457 x > 1668310833 > voxel sizes: 72950800688805540228332034981888.0000, > 0.0000, > 70387746064971392794503413760.0000 > type: UNKNOWN (842347100) > fov: inf > dof: 1668244335 > xstart: -inf, xend: inf > ystart: -39.8, yend: 39.8 > zstart: > -57528904905339976080461439864148066304.0, zend: > 57528904905339976080461439864148066304.0 > TR: 71446969504754118412748034408448.00 > msec, TE: > 17862977898376175754203558313984.00 msec, TI: > 4464583905807213894543171649536.00 msec, flip angle: > 4128358109894647822155776.00 degrees > nframes: 1668310833 > PhEncDir: UNKNOWN > ras xform present > xform info: x_r = > 19441367164500690594879212879872.0000, y_r = > 282670772483715414097920.0000, z_r = > 4630517930517356511053742080.0000, c_r = > 18179207505345150839094025125888.0000 > : x_a = > 72953300747400503281465327353856.0000, y_a = > 71756401321839400437789128392704.0000, z_a = > 4721224652433019196440490016768.0000, c_a = 0.0000 > : x_s = 0.0000, y_s = > 208210237828607676252160.0000, z_s = > 4286666299869639999488.0000, c_s = > 4622889417339745644740714954752.0000 > > talairach xfm : > Orientation : AAA > Primary Slice Direction: coronal > > voxel to ras transform: > inf 12133604733747200.0000 inf > -inf > inf 3080133850722982464323584.0000 > inf > -inf > 36752183205510726549504.0000 > 8937395835633664.0000 inf -inf > 0.0000 0.0000 0.0000 1.0000 > > voxel-to-ras determinant nan > > ras to voxel transform: > mat = NULL! > Invalid argument > MatrixFree: NULL mat POINTER! > > Invalid argument > > -Andrew > > ------------------------------------------------------------------------ > > _______________________________________________ > Freesurfer mailing list > Freesurfer(a)nmr.mgh.harvard.edu > <mailto:Freesurfer@nmr.mgh.harvard.edu> > <mailto:Freesurfer@nmr.mgh.harvard.edu > <mailto:Freesurfer@nmr.mgh.harvard.edu>> > > > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer > > > -- Douglas N. Greve, Ph.D. > MGH-NMR Center > greve(a)nmr.mgh.harvard.edu > <mailto:greve@nmr.mgh.harvard.edu> > <mailto:greve@nmr.mgh.harvard.edu > <mailto:greve@nmr.mgh.harvard.edu>> > Phone Number: 617-724-2358 <tel:617-724-2358> > <tel:617-724-2358 <tel:617-724-2358>> Fax: 617-726-7422 > <tel:617-726-7422> > <tel:617-726-7422 <tel:617-726-7422>> > > > Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting > <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> > <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting> > > FileDrop: > www.nmr.mgh.harvard.edu/facility/filedrop/index.html > <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> > <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> > > > > > 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. > > > > -- > 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: www.nmr.mgh.harvard.edu/facility/filedrop/index.html > <http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html> > > -- 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: www.nmr.mgh.harvard.edu/facility/filedrop/index.html

1 0

display multiple functional overlays simultaneously
by Scott Hayes 15 Mar '11

15 Mar '11

Hi, I have 2 functional maps (Patient > Controls and Controls > Patients) and I would like to view both maps simultaneously on the surface (one with "normal" values and the other using the "inverse" values so that I get a "warm" and "cool" map simultaneously displayed on the surface). I'm using FS 5.0 and mac os 10.6.6 command from the mac terminal: > tksurfer fsaverage/ lh inflated -overlay PCC_nosmooth_level3_FE_8matched_vs_MA016.gfeat/cope1.feat/thresh_zstat1.nii.gz -overlay PCC_nosmooth_level3_FE_8matched_vs_MA016.gfeat/cope1.feat/thresh_zstat2.nii.gz Both maps load, but I have to click on one or the other (rather than both being overlaid simultaneously), and the options (selecting inverse) seem to be applied to both layers (even though I have not selected "apply changes to all layers"). Is it possible to view these maps at the same time using tksurfer and apply normal values to one map and inverse values to the other? It was also suggested that freeview *might* be able to do this, but this does not launch from the command line (FS, tkmedit and tksurfer all work fine). I've copied the first few lines from the freeview crash if that helps. Any suggestions/guidance would be appreciated. Thanks, Scott Process: freeview.bin [36175] Path: /Applications/freesurfer/bin/Freeview.app/Contents/MacOS/freeview.bin Identifier: freeview.bin Version: 'latest trunk' ('latest trunk') Code Type: X86 (Native) Parent Process: bash [30033] Date/Time: 2011-03-15 08:57:43.954 -0400 OS Version: Mac OS X 10.6.6 (10J567) Report Version: 6 Interval Since Last Report: 1791542 sec Crashes Since Last Report: 39 Per-App Crashes Since Last Report: 3 Anonymous UUID: E5903541-C168-4ED7-92B6-6741CD098D1A Exception Type: EXC_BAD_ACCESS (SIGBUS) Exception Codes: KERN_PROTECTION_FAILURE at 0x0000000000000000 Crashed Thread: 0 Dispatch queue: com.apple.main-thread Thread 0 Crashed: Dispatch queue: com.apple.main-thread 0 ??? 0000000000 0 + 0 1 libvtkRendering.5.2.dylib 0x02b3d1a9 vtkCarbonRenderWindow::CreateAWindow() + 281 2 libvtkRendering.5.2.dylib 0x02b3bfa4 vtkCarbonRenderWindow::WindowInitialize() + 52 3 libvtkRendering.5.2.dylib 0x02b3a435 vtkCarbonRenderWindow::Start() + 21 4 libvtkRendering.5.2.dylib 0x02aa11db vtkRenderWindow::DoStereoRender() + 27 5 libvtkRendering.5.2.dylib 0x02aa13a9 vtkRenderWindow::DoFDRender() + 153 6 libvtkRendering.5.2.dylib 0x02aa1832 vtkRenderWindow::DoAARender() + 98 7 libvtkRendering.5.2.dylib 0x02a9fc54 vtkRenderWindow::Render() + 644 8 freeview.bin 0x0011d547 wxVTKRenderWindowInteractor::OnPaint(wxPaintEvent&) + 151 9 freeview.bin 0x00407b53 wxAppConsole::HandleEvent(wxEvtHandler*, void (wxEvtHandler::*)(wxEvent&), wxEvent&) const + 65 10 freeview.bin 0x00452f30 wxEvtHandler::ProcessEventIfMatches(wxEventTableEntryBase const&, wxEvtHandler*, wxEvent&) + 212 11 freeview.bin 0x004532d4 wxEventHashTable::HandleEvent(wxEvent&, wxEvtHandler*) + 204 12 freeview.bin 0x00453f5e wxEvtHandler::ProcessEvent(wxEvent&) + 240 13 freeview.bin 0x003525e3 wxWindow::MacDoRedraw(void*, long) + 975 14 freeview.bin 0x00353985 wxMacWindowControlEventHandler(OpaqueEventHandlerCallRef*, OpaqueEventRef*, void*) + 815 15 freeview.bin 0x00354135 wxMacWindowEventHandler(OpaqueEventHandlerCallRef*, OpaqueEventRef*, void*) + 181 16 com.apple.HIToolbox 0x95953ecf DispatchEventToHandlers(EventTargetRec*, OpaqueEventRef*, HandlerCallRec*) + 1567 17 com.apple.HIToolbox 0x95953196 SendEventToEventTargetInternal(OpaqueEventRef*, OpaqueEventTargetRef*, HandlerCallRec*) + 411

2 5

corpus calosum volume ?...
by Gonzalo Rojas Costa 15 Mar '11

15 Mar '11

Hi: The CC_Posterior, CC_Central, etc, volumes in "aseg.mgz" are the complete CC volume, or it includes only some sagital slices ?... How can I meassure the CC volume for some sagital slices only (5 in total for example) ?... When I tried using "mri_cc" I got an error (for example) in the case that the patient is slightly rotated (freesurfer segments only one section of the CC as CC and the other as WM)... Sincerely, Gonzalo Rojas Costa Laboratory for Advanced Medical Image Processing Department of Radiology Las Condes Clinic Lo Fontecilla 441, Las Condes, Santiago, Chile. Tel: 56-2-2105170 Cel: 56-9-97771785 www.clc.cl

2 1

V1 prediction
by Thesen, Thomas 15 Mar '11

15 Mar '11

Dear Freesurfer List, I ran the script for predicting V1 location on the surface: http://surfer.nmr.mgh.harvard.edu/fswiki/V1 According to the paper, the output should be a probability map, but the script only seems to produce a label file: $SUBJECT_DIR/label/?h.v1.prob.label . When I load this file (it only loads as 'label') it is a binary label, with no probability values associated with it. That means it is a pretty generous description of predicted V1. How can I get probability values and trim it to, let's say, >80% probability? Thanks for your input. Thomas </PRE> <html> <body> ------------------------------------------------------------<br /> This email message, including any attachments, is for the sole use of the intended recipient(s) and may contain information that is proprietary, confidential, and exempt from disclosure under applicable law. Any unauthorized review, use, disclosure, or distribution is prohibited. If you have received this email in error please notify the sender by return email and delete the original message. Please note, the recipient should check this email and any attachments for the presence of viruses. The organization accepts no liability for any damage caused by any virus transmitted by this email.<br /> ================================= </body> </html> <PRE>

3 2

tksurfer error
by Don Hagler 15 Mar '11

15 Mar '11

Hi, In response to the message below, I have been getting the exact same error when running tksurfer on a particular workstation. tksurfer does display correctly, but the drawing is very slow, as if running remotely (or not using the video card). This happens with versions 3.0.5 and 4.5.0. I tried uninstalling and reinstalling the nvidia driver, but that didn't help. glxgears and glxinfo work fine, no warnings or error messages. In Xorg.0.log, there is one warning: (WW) NVIDIA(0): UBB is incompatible with the Composite extension. Disabling (WW) NVIDIA(0): UBB. After googling, I found that UBB stands for Unified Back Buffer. http://http.download.nvidia.com/XFree86/Linux-x86/1.0-8178/README/appendix-… This only applies to Quadro GPUs. This particular workstation has this type: NVIDIA GPU Quadro FX 580 I haven't yet found a solution. Evan, do you have a Quadro GPU? Did you find a solution? Don _________________________________________________________________________ Dear freesurfer community, I am currently trying to view one of the surfaces generated by recon-all command in tksurfer but when I enter the command: tksurfer LuceHR_FS_SURF lh inflated.nofix I receive a series of errors that look like this: Received X error! Error code : 1 Request code : 147 Minor code : 181 Error text : 'BadRequest (invalid request code or no such operation)' Subsequently the surface window does not display what I would believe is the surface I am trying to view and then the only way I can exit the tksurfer environment is by exiting the terminal that I invoked it in. Any ideas or solutions? Thank you very much, Evan Luther

1 1

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