Freesurfer

freesurfer@nmr.mgh.harvard.edu

19 participants
30088 discussions

Help ! Error for generating *.mgz files
by xiaotian yin 23 Feb '11

23 Feb '11

Dear Bruce, Thank you so much for taking care of my stupid questions :) I searched online regarding the ANALYZE vs Dicom, and I understand your point now ... The problem is, for this scan, I am afraid the only data format that I can play with is ANALYZE ... (although there might be Dicom format available on later scans) So ... I tried some other tools and got the orientation information out of the ANALYZE data (not sure if it makes sense ...) 0.7000 0 0 0 0 0.7000 0 0 0 0 1.5000 0 0 0 0 1.0000 I save them in a *.mat file, and seems that mri_convert did the job successfully, I got the *.mgz file. Then I setup the necessary directories ( i.e. my_new_subjects_dir/subj001/mri/orig ), put the *.mgz file there, and launched " recon-all -all -s subj001 ". Then I got some warnings and errors (see below), and here comes my new questions, if you do not mind ... ** Question-1 : After launching recon-all, got a warning immediately (see below). Sorry for my stupid questions ... but what is a typical input I should feed to recon-all ? Two scans of the same whole brain? Or one scan for the left hemisphere scan and another one for the right hemisphere? Or something else ? "WARNING: only one run found. This is OK, but motion correction cannot be performed on one run, so I'll copy the run to rawavg and continue." ** Question-2 : The program seems running well after the warning, until it terminates a couple of minutes later, throwing out some error like below. I just have no idea what to do now ... Could you please give me some hints? "ERROR: talairach_afd: Talairach Transform: transforms/talairach.xfm ***FAILED*** (p=0.0000, pval=0.0000 < threshold=0.0050)" The complete output from recon-all is attached, for your information. If you need the data, is there any way I can send it to you without bothering other members in the mailing list? Thank you ! Best Regards Tim On Sun, Feb 20, 2011 at 9:47 AM, Bruce Fischl <fischl(a)nmr.mgh.harvard.edu> wrote: > Hi Tim, > > it doesn't help to convert from analyze to something else. Once you have > gone through analyze all direction information is lost. Can you get access > to the data before it was converted to analyze in whatever format the > scanner produced? > Bruce > > On Sun, 20 Feb 2011, xiaotian yin wrote: > >> Dear Bruce, >> >> >> THanks for the quick response! >> >> I searched online for Analyze --> Dicom converter ... no luck ... >> >> Also tried to read the Analyze file (*.hdr) into Matlab, ending up >> with the following info : >> >> Filename: 'IA_1_orig.hdr' >> FileModDate: '04-ÿÿÿÿ-2010 01:08:20' >> HdrFileSize: 348 >> ImgFileSize: 8126464 >> Format: 'Analyze' >> FormatVersion: '7.5' >> Width: 256 >> Height: 256 >> BitDepth: 8 >> ColorType: 'grayscale' >> ByteOrder: 'ieee-le' >> HdrDataType: 'volume ' >> DatabaseName: 'IA_1_%%%.pic.' >> Extents: 16384 >> SessionError: -28067 >> Regular: 1 >> Dimensions: [256 256 124 1 1] >> VoxelUnits: '' >> CalibrationUnits: 'Y??·2w?·' >> ImgDataType: 'DT_UNSIGNED_CHAR' >> PixelDimensions: [1x5 single] >> VoxelOffset: 0 >> CalibrationMax: 255 >> CalibrationMin: 0 >> Compressed: -2.7107e-005 >> Verified: -2.7172e-005 >> GlobalMax: 255 >> GlobalMin: 0 >> Descriptor: [1x80 char] >> AuxFile: '?J?· ?ù·? ?oü·????°???' >> Orientation: 'Orientation unavailable' >> Originator: ' ·???' >> Generated: ' p' >> Scannumber: '??· ' >> PatientID: '' >> ExposureDate: '' >> ExposureTime: '' >> Views: 0 >> VolumesAdded: 0 >> StartField: -1209816480 >> FieldSkip: -1208305456 >> OMax: -1208193036 >> OMin: -1208460092 >> SMax: -1208306728 >> SMin: -1077308224 >> >> >> >> Did not see any "hist.orient" field or affine transformation field ... >> >> Do you have any recommendations ? >> >> >> Best >> Tim >> >> >> >> >> On Sat, Feb 19, 2011 at 9:46 PM, Bruce Fischl >> <fischl(a)nmr.mgh.harvard.edu> wrote: >>> >>> Hi Tim, >>> >>> the .hdr/.img format is called Analyze, and it doesn't come with >>> orientation informaton. This means that it is easy for example to get a >>> left/right swap that you will never know was in your data. Can you get it >>> in any other format, like dicom, or whatever format it came off the >>> scanner >>> in? >>> >>> If not, you can send us one of the pairs of files that gives the segfault >>> and we'll take a look >>> >>> cheers >>> Bruce >>> >>> >>> >>> On Sat, 19 Feb 2011, xiaotian yin wrote: >>> >>>> Hi, >>>> >>>> Thank you for reading this email ! >>>> >>>> I am a newbie to Medical Imaging, and am learning to use FreeSurfer >>>> for surface construction. >>>> >>>> Already installed FreeSurfer in Ubuntu, and the installation has been >>>> verified to be 100% successful. >>>> >>>> I got a bunch of data (coming in pairs of *.hdr and *.img), and I am >>>> trying to convert them to *.mgz files, so that I can use recon-all >>>> then. >>>> >>>> but when I call mri_convert as follows: >>>> >>>> mri_convert Rcing.img Rcing.mgz >>>> >>>> I got the following error : >>>> >>>> >>>> mri_convert Rcing.img Rcing.mgz >>>> $Id: mri_convert.c,v 1.166.2.2 2010/08/10 19:11:50 greve Exp $ >>>> reading from Rcing.img... >>>> ----------------------------------------------------------------- >>>> INFO: could not find /mnt/Data/Rcing.mat file for direction cosine info. >>>> INFO: use Analyze 7.5 hdr->hist.orient value: -40, . >>>> INFO: if not valid, please provide the information in >>>> /mnt/Data/Rcing.mat file >>>> ----------------------------------------------------------------- >>>> Segmentation fault >>>> >>>> >>>> I have searched the freesurfer website ... no luck .... >>>> Would you mind to give me some advice or hints ? Thank you so much ! >>>> >>>> >>>> Best >>>> Tim >>>> _______________________________________________ >>>> 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 >>> >>> >>> 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. >>> >>> >> >> >

2 1

Neuroimaging Data Processing Workshop 23-24 March 2011, Edinburgh, UK
by Chris Filo Gorgolewski 23 Feb '11

23 Feb '11

Dear all, I am delighted to invite you to a Neuroimaging Data Processing Workshop which will be held on 23-24th of March 2011 at University of Edinburgh, UK. During the workshop we will cover neuroimaging data analysis and management with reusable, scalable and reproducible processing workflows that combine different software packages (such as SPM, FSL, and FreeSurfer) using NiPyPE framework. The workshop will also feature talks and demos covering data management, encapsulation and visualization using PyXNAT and the Connectome Library and Viewer. More details can be found at the workshop website: http://nipype.blogspot.com Best regards, Chris Gorgolewski

1 1

longitudinal analysis
by nadia merhoum 23 Feb '11

23 Feb '11

Hello I have just started to use FreeSurfer (version 4.5) and have some questions on the longitudinal stream: I have started to process cross-sectionally all time points (7 tp) of my patients with the following script: recon-all –i tpidN.nii –subjid tpidN* –all* Then, created an unbiased template from all time points for each subject and process it with recon-all with the following script: recon-all -base <templateid> -tp <tp1id> -tp <tp2id> ... -all and finally did the longitudinal process to all time points with the following script: recon-all -long <tpNid> <templateid> -all I would like to know if what I have been doing till this final step is correct or not? I would like also to know how we could calculate the differences between tpid1.long.template and tpid2.long.template for example? And is there any automatic way to do so? Thank you for your help Sincerely Nadia

2 1

mkanalysis-sess says ERROR: bad format for seq.info file
by Pablo Polosecki 23 Feb '11

23 Feb '11

Hi all, I am trying to define an analysis for a block-based functional scan using mkanalysis-sess. When I run it, I get the following output with an error message: Setting TER to .10000000000000000000 /home/ppolosecki/space/data/ppolosecki/cooked/2011/110211Quincy ERROR: bad format for seq.info file, /home/ppolosecki/space/data/ppolosecki/cooked/2011/110211Quincy The seq.info file (attached) was produced by the unpacksdcmdir command and it looks no different from other ones I have seen. The information there seems correct to me. Could someone perhaps point out possible reasons why this could happen? The script I use to call mkanalysis-sess containing all the parameters used is attached. Thank you in advance. Pablo -- Pablo Polosecki Graduate Fellow Laboratory of Neural Systems The Rockefeller University 1230 York Avenue, Box 9 New York, NY 10065 Tel: 212-327-7699 (Lab) 917-842-2510 (Cell) email: ppolosecki(a)mail.rockefeller.edu

2 10

recon-all on data with anatomical abnormality or strong artifact
by Eveline Geiser 22 Feb '11

22 Feb '11

Hi, I run into the following error when running recon-all. ERROR: _FindFacePath: could not find path! in the process: Correcting Topology of defect 2 with euler number -9 (5 loops) computing statistics for defect 2: 1934 vertices location: [ (93,119,57) - average intensity = 85.136 ].... The subject also has a huge anatomical abnormality or possibly an artifact ventral of the cerebellum (between cerebellum and pons). Brainmask.mgz also shows striping. When running the old topology fixer the following error occurs: mris_fix_topology: could not read brain volume from /mindhive/gablab/rhythm/cacophonix_July10/freesurfer_data/cacix7/mri/brain however, brain.mgz is there and loadable in tkmedit. Could the artifact be the reason for all this or is there another issue? Thanks for your help Eveline -- --------------------- Eveline Geiser PhD Massachusetts Institute of Technology Department of Brain and Cognitive Sciences phone: +617 324 63 71 e-mail: egeiser(a)mit.edu http://web.mit.edu/gabrieli-lab/

2 1

question about weighted regression analysis
by Katie Bettencourt 22 Feb '11

22 Feb '11

So I created a weighted regression analysis to look at the effect of memory load in a particular brain region. Basically, I weighted the paradigms by a behavioral measure that reflected the number of items actually remembered (as set size was increased). As far as Doug told me there are basically 2 ways to weight your paradigm files. Version 1: Have 2 conditions, baseline (condition 0) and all the set sizes (condition 1). Condition 1 would then be weighted by the behavioral measure. Version 2: Have 3 conditions, baseline (condition 0), and then I represented each presentation as two different conditions, one with a weight that is always 1 (condition 1), the other weighted according to the behavioral measure (condition 2). The difference, as far as I understand it, in version 1, it is assumed that the response amplitude is ) when the weight is 0 (ie. that when you are attending to 0 items, brain activity = 0). Whereas, version 2, tests the slope of the HRF amplitude vs weight without the assumption above. However, I'm a bit confused as to the results I got. When I looked at the data from both versions, version 1 provided a much higher amount of activation and more areas activated than version 2. However, I believe version 2 better fits with the multiple regression analysis that is done in Brain Voyager. Can anyone give me a better explanation of what the difference between these analysis models is? Katie

2 3

mri_glmfit-sim doubt
by Rafa x 22 Feb '11

22 Feb '11

Dear Freesurfers, I have performed a permutation test with the command: mri_glmfit-sim --glmdir rh.fsgd_subj_gauss005 --sim mc-z 10000 THRESHOLD csd --sim-sign pos --seed 1297708255 And after that I have obtained the clusters with: mri_surfcluster --subject CNTRL_MCI_RAFA_average –mask rh.fsgd_subj_gauss005/mask.mgh --sum rh.fsgd_subj_gauss005/contrast/clustsum.txt --in rh.fsgd_subj_gauss005/contrast/sig.mgh --cwsig rh.fsgd_subj_gauss005/contrast/cwsig.mgh --vwsig rh.fsgd_subj_gauss005/contrast/vwsig.mgh --sign pos --csd rh.fsgd_subj_gauss005/csd/csd.j001-contrast.csd I’ve did it over the same dataset but with two different thresholds, THRESHOLD: 2.301 (pvalue=0.005) and 1.301 (pvalue=0.05). As far I know the cwsig.mgh resulting using the 1.3 threshold should be similar than the cwsig.mgh from the second thr, but with a more restrictive significative clusters (some coloured areas should disappear, in other words, clusters returned in the first case should be included in clusters returned by the 1.3 threshold). My results didn’t show that behaviour, some of the cwsig clusters can be contained in the other, but others don’t (see attached Figure). It is my reasoning correct? If so, why the images mismatching? Thank you all in advance. Rafa.

3 2

Gammafit parameters
by Mandy Nagy 22 Feb '11

22 Feb '11

Hi all, We are trying to use -gammafit when running the mkanalysis-sess command. In the FsFast Tutorial, the flag -gammafit 2.25 1.25 is used. Are those numbers the default? If not, what is the default? Thanks in advance. Best, Mandy Nagy

2 1

overlay sig.mgh onto T1 vol
by Fernanda Palhano 22 Feb '11

22 Feb '11

Hello, Does anyone knows if is it possible to overlay the surface file sig.mgh onto a T1 vol? I tried the command: tkmedit fsaverage T1.mgz -fthresh 2 -fmax 5 -overlay qdec/rh-Diff-grupo1-grupo2-Intercept-thickness/sig.mgh (where 2 and 5 were the thresholds I used in qdec), but it seems like I have only one point of the sig.mgh (in the slice 128) which is outside the volume. There's another way to do this? Thanks for the help! Fernanda

4 5

FW: Surface Generation Problem Pertaining to FreeSurfer
by Saad Shakeel 22 Feb '11

22 Feb '11

From: saad.akhun(a)hotmail.com To: fsregister(a)nmr.mgh.harvard.edu; freesurfer(a)nmr.mgh.harvard.edu Subject: Surface Generation Problem Pertaining to FreeSurfer Date: Tue, 22 Feb 2011 11:43:36 +0500 Hello, We are graduate students from Giki university, Pakistan. We are doing our Final Year Project on shape Analysis Techniques.Our Project steps are a. Segmentation b. Surface Generation c. Spherical Harmonics for shape Analysis. We need your help pertaining to Surface Generation module. I have certain problems:1. How to create : surface generation of hippo-campus ( Exists In Sub-cortical region of Human Brain)2. Can I use mri_vol2surf tool in FreeSurfer for this purpose. If so, then what are the advantages and dis-advantages Thankyou for your time. Kindly send me the supporter's email address if you don't deal with this kind of situation.

2 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