Freesurfer

freesurfer@nmr.mgh.harvard.edu

3 participants
30135 discussions

Re: [Freesurfer] trac-all -prep error field map error
by Eileen Moore 05 Sep '14

05 Sep '14

Hi - I still have not been able to apply the field map correction within tracula. Tracula does not appear to recognize the format of my field maps in either dicom or nifti formats. However, as a work around I first did the field map correction outside tracula (using our wrapper script for epidewarp). Next I used this field map corrected image as the input diffusion image for tracula and set dob0 = 0 in my dmrirc file (so as to skip tracula's field map correction). This appears to be working. Eileen. > > ------------------------------ > > Message: 5 > Date: Wed, 27 Aug 2014 13:51:06 -0700 > From: Eileen Moore <eimoore(a)mail.sdsu.edu> > Subject: Re: [Freesurfer] trac-all -prep error field map error > To: "freesurfer(a)nmr.mgh.harvard.edu" <freesurfer(a)nmr.mgh.harvard.edu> > Message-ID: > <CA++dajQo_2EF= > LvUFCKbs+_061v7dcTFC5T11PN0K4G9d-CuAw(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Anastasia, I'm not sure how I can give tracula the phase difference maps, > can you provide an example? > > As a reminder, I've run a wrapper script around epidewarp which takes our > dicom images and converts them into a phase difference map (it is based on > the script that was authored by Doug Greve for the fBIRN consortium). > > In the tutorial dmrirc example I see only the 'set b0mlist' and 'set > b0plist' options, and am unsure where I could give it the phase difference > map. > > I've also tried converting the dicom magnitude and phase maps to nifti > format and inputting these into the dmrirc file but I am still getting an > error. (See attached log file). > > Thanks, > Eileen. > > > > > ------------------------------ > > > > Message: 10 > > Date: Thu, 17 Jul 2014 04:14:16 -0400 (EDT) > > From: Anastasia Yendiki <ayendiki(a)nmr.mgh.harvard.edu> > > Subject: Re: [Freesurfer] trac-all -prep error field map error > > To: Freesurfer support list <freesurfer(a)nmr.mgh.harvard.edu> > > Message-ID: <alpine.LRH.2.03.1407170411110.12767(a)nmr.mgh.harvard.edu> > > Content-Type: TEXT/PLAIN; charset=US-ASCII; format=flowed > > > > > > Hi Eileen - All image formats that mri_convert (and thus any freesurfer > > program) can read will do, you don't need dicom. Try giving trac-all the > > phase difference maps you have and let's see what happens. > > > > Best, > > a.y > > > > On Wed, 16 Jul 2014, Eileen Moore wrote: > > > > > Hi Anastasia, > > > > > > OK - our imaging physicists had developed a wrapper script around > > epidewarp > > > which takes our dicom images and converts them into the phase > difference > > map > > > (it is based on the script that was authored by Doug Greve for the > fBIRN > > > consortium). It looks like the tracula script requires dicom images. Is > > > there a way for me to do epidewarp manually and then follow through > with > > the > > > remaining steps using trac-all? > > > > > > Thanks, > > > Eileen. > > > > > > > > > Date: Mon, 14 Jul 2014 11:23:47 -0400 (EDT) > > > From: Anastasia Yendiki <ayendiki(a)nmr.mgh.harvard.edu> > > > Subject: Re: [Freesurfer] trac-all -prep error field map error > > > To: Freesurfer support list <freesurfer(a)nmr.mgh.harvard.edu> > > > Message-ID: > > > <alpine.LRH.2.03.1407141119560.23459(a)nmr.mgh.harvard.edu> > > > Content-Type: text/plain; charset="utf-8" > > > > > > > > > Hi Eileen - That step is run with the program epidewarp.fsl, > > > which you can > > > run on the command line to see the available options. The input > > > phase map > > > volume is expected to have either 2 frames (2 separate phase > > > maps) or 1 > > > frame (the difference phase map). This error from trac-all means > > > that it > > > found that your input phase map has something other than 1 or 2 > > > frames, so > > > it doesn't know how to handle it. > > > > > > Hope this helps, > > > a.y > > > > > > On Sat, 12 Jul 2014, Eileen Moore wrote: > > > > > > > Hi - I'm using freesurfer 5.3 on an Ubuntu operating system. > > > I'm tyring to > > > > run tracula on some? images acquired from a GE scanner. I am > > > running into > > > > problems when I try to perform the registration-based > > > B0-inhomogeneity > > > > compensation, getting an error in the terminal output: "ERROR: > > > Unrecognized > > > > format of phase map". However, trac-all -prep runs fine when I > > > skip the > > > > B0-inhomogeneity compensation step. Is there something else I > > > need to do to > > > > use field map correction? > > > > > > > > I've attached the trac-all.log and a copy of the terminal > > > output. > > > > > > > > Thanks, > > > > Eileen. > > > > > > > > > > > > > > > > > > > > > > > > > > > >

1 0

Re: [Freesurfer] longitudinal QDEC help
by Martin Reuter 05 Sep '14

05 Sep '14

Hi Taylor, You can do the glmfit osgm on whatever you like (rate, percentchange pc1, or symm. percent change spc). I recommend SPC. The 2stage model page shows an example of glmfit with osgm flag and explains how to stack your data first. Best Martin Sent via my smartphone, please excuse brevity. <div>-------- Original message --------</div><div>From: Taylor North <taylornorth189(a)hotmail.ca> </div><div>Date:09/04/2014 11:14 PM (GMT-05:00) </div><div>To: Freesurfer support list <freesurfer(a)nmr.mgh.harvard.edu> </div><div>Subject: Re: [Freesurfer] longitudinal QDEC help </div><div> </div>Hi Martin Yes I that definitely clears things up, thanks very much for the help. As for the QCAHCE we did use the long_mris_slopes described in the tutorial you linked in your email. Just to confirm I believe I should be using the mri_glmfit -osgm model for my one sample group mean analysis and then when I want to compare against my continuous/dependent variables I should be using QDEC? Also it looks like mri_glmfit -osgm calculates % change, so should I be analyzing % change, not symmetrized percent change on my QDEC analysis as well, even though it say symmetrized percent change is the more robust result? Thanks again for your help Taylor From: mreuter(a)nmr.mgh.harvard.edu Date: Thu, 4 Sep 2014 22:44:16 -0400 To: freesurfer(a)nmr.mgh.harvard.edu Subject: Re: [Freesurfer] longitudinal QDEC help Hi Taylor, first, in your description you mention 'base and 6 month'. Better call it 'baseline' or 'first and second time points at 0 and 6 months' in order to not confuse the first time point with the 'base' = within-subject template created as the second step. I think you did not confuse it, but just want to make sure. From your description it is not clear how you computed the symmetrized percent change. You say you ran QCACHE, but how? Did you use long_mris_slopes as described here: http://freesurfer.net/fswiki/LongitudinalTwoStageModel http://freesurfer.net/fswiki/FsTutorial/LongitudinalTutorial That step computes slopes for each subject and (instead of having 2 thickness maps) reduces the information to a single map (here the symmetrized percent change is simply: 100 * tp2-tp1/(0.5*(tp1+tp2)) ) . These maps are stored in the base directories (one for each subject and of corse each hemisphere that you run). Now to your questions: 1) I think a 'one sample group mean' analysis is not possible in QDEC. For that you have to use mri_glmfit -osgm …. I wish we had this functionality and I'll see how hard it would be to implement it in future versions of QDEC (this came up frequently and it should not be too hard). 2) a) the longitudinal qdec table use the cross id's in the fsid column for simplicity. When it is used in many programs (e.g. long_mris_slopes, long_stats_slopes and other scripts with the long-qdec option) the scripts knows that it is a longitudinal qdec table (because of the existing fsid-base column) and internally use the longitudinal directories for processing. The cross sectional directories are never used. b) as mentioned above, the results from the first stage analysis (creation of the within subject slopes and percent changes) are stored in the base directories. They are used in this step and found there. (and you are right, the recon-all results from the base should not be analyzed directly, but here the data that is used is derived from the longitudinal directories and just placed in the base so that it can be found for each subject). Hope that explains it. Best, Martin On Sep 4, 2014, at 9:33 PM, Taylor North <taylornorth189(a)hotmail.ca> wrote: Hello, I am having some troubles using QDEC to analyze the cortical thickness changes in my group study and I was hoping someone may be able to give me a hand. The steps I have done are as follows: I am working on a study looking at 40 subjects scanned at two time points (base and 6 months after). I have followed all the steps for longitudinal processing for each subject (cross,base, and longitudinals) and made edits where they were necessary to the base. I then entered all 40 subjects into a QDEC table as follows fsid fsid-base time gender age crosstpN1 base1 0 female 20 crosstpN2 base1 6 female 20 .......... continues for rest of the subjects then I ran the QCHACHE on the data, made a Qdecrc. file in the subject DIR, and then transformed the longitudinal QDEC table into a QDEC table which looks like: fsid time gender age Base1 3 female 20 .......... continues for the rest of the subjects. I then followed the QDEC steps outlined In the group tutorial. (import cross sectional table, generate stats data table, import aseg volume, set longitudinal dependent variables and pick continuous variable. 1) the first problem I am having is that when I go to process the longitudinal group results through QDEC it makes me pick a continuous variable when I just want to analyze the dependent variable The symmetrized percent change. So I have tried picking different continuous variables ( such as age or importing rhcortex volume etc...) with choosing symmetrized percent change as the dependent variable. No matter what continuous variable I pick the question: ``does the average symmetrized percent change differ from zero"? always comes up, the problem is depending on what continuous variable I pick the results change. Is there a way to just look at the longitudinal thickness changes across a group alone, because I first want to look at the atrophy across the whole group over 6 months and then compare the differences in atrophy against the diagnosed and non diagnosed group. 2) another thing I don't quite understand is it seems like the longitudinal qdec tables are made with the cross of the subjects under fsid, but why wouldn't you use the longitudinals under the fsid for the longitudinal table instead? Also it seems like in the cross sectional qdec table the base is the subject being analyzed under fsid, which also doesn't make sense to me because in the tutorial it says base should not be used directly for analysis. thanks for the help! Taylor _______________________________________________ Freesurfer mailing list Freesurfer(a)nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer --------------------------------- Dr. Martin Reuter Assistant in Neuroscience - Massachusetts General Hospital Instructor in Neurology - Harvard Medical School MGH / HMS / MIT A.A.Martinos Center for Biomedical Imaging 149 Thirteenth Street, Suite 2301 Charlestown, MA 02129 Phone: +1-617-724-5652 Email: mreuter(a)nmr.mgh.harvard.edu reuter(a)mit.edu Web : http://reuter.mit.edu _______________________________________________ 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.

1 0

longitudinal QDEC help
by Taylor North 04 Sep '14

04 Sep '14

Hello, I am having some troubles using QDEC to analyze the cortical thickness changes in my group study and I was hoping someone may be able to give me a hand. The steps I have done are as follows: I am working on a study looking at 40 subjects scanned at two time points (base and 6 months after). I have followed all the steps for longitudinal processing for each subject (cross,base, and longitudinals) and made edits where they were necessary to the base. I then entered all 40 subjects into a QDEC table as follows fsid fsid-base time gender age crosstpN1 base1 0 female 20 crosstpN2 base1 6 female 20 .......... continues for rest of the subjects then I ran the QCHACHE on the data, made a Qdecrc. file in the subject DIR, and then transformed the longitudinal QDEC table into a QDEC table which looks like: fsid time gender age Base1 3 female 20 .......... continues for the rest of the subjects. I then followed the QDEC steps outlined In the group tutorial. (import cross sectional table, generate stats data table, import aseg volume, set longitudinal dependent variables and pick continuous variable. 1) the first problem I am having is that when I go to process the longitudinal group results through QDEC it makes me pick a continuous variable when I just want to analyze the dependent variable The symmetrized percent change. So I have tried picking different continuous variables ( such as age or importing rhcortex volume etc...) with choosing symmetrized percent change as the dependent variable. No matter what continuous variable I pick the question: ``does the average symmetrized percent change differ from zero"? always comes up, the problem is depending on what continuous variable I pick the results change. Is there a way to just look at the longitudinal thickness changes across a group alone, because I first want to look at the atrophy across the whole group over 6 months and then compare the differences in atrophy against the diagnosed and non diagnosed group. 2) another thing I don't quite understand is it seems like the longitudinal qdec tables are made with the cross of the subjects under fsid, but why wouldn't you use the longitudinals under the fsid for the longitudinal table instead? Also it seems like in the cross sectional qdec table the base is the subject being analyzed under fsid, which also doesn't make sense to me because in the tutorial it says base should not be used directly for analysis. thanks for the help! Taylor

2 2

recon-all input
by Xinfa Shi 04 Sep '14

04 Sep '14

Hello, Now I have a subject data,including 156 dicom images(MR0001~MR0156),amd they were converted into a .nii.gz file(eg.CON_T1_01.nii.gz) .I am wondering whether they are the same : recon-all -i MR00001.dcm -s con01 -all recon-all -i CON_T1_01.nii.gz -s con01 -all Many Regards. Xin-Fa Shi.

2 1

invalid=96777 in recon-all
by Pom Sailasuta 04 Sep '14

04 Sep '14

Hello, I am using FreeSurfer version 5.3.0 in 64bit desktop, I got these errors at the begining of recon-all but the reconstruction is not stop, keep on going with lots of invalid=96777. **************** pass 1 of 1 ************************ setting smoothness coefficient to 0.039 blurring input image with Gaussian with sigma=2.000... 0000: dt=0.000, rms=0.849, neg=0, invalid=96777 0001: dt=129.472000, rms=0.789 (7.080%), neg=0, invalid=96777 0002: dt=295.936000, rms=0.765 (3.083%), neg=0, invalid=96777 0003: dt=129.472000, rms=0.759 (0.766%), neg=0, invalid=96777 0004: dt=517.888000, rms=0.749 (1.353%), neg=0, invalid=96777 0005: dt=92.480000, rms=0.744 (0.567%), neg=0, invalid=96777 0006: dt=369.920000, rms=0.740 (0.600%), neg=0, invalid=96777 0007: dt=73.984000, rms=0.738 (0.229%), neg=0, invalid=96777 0008: dt=73.984000, rms=0.737 (0.148%), neg=0, invalid=96777 0009: dt=73.984000, rms=0.735 (0.250%), neg=0, invalid=96777 0010: dt=73.984000, rms=0.733 (0.313%), neg=0, invalid=96777 0011: dt=73.984000, rms=0.730 (0.371%), neg=0, invalid=96777 0012: dt=73.984000, rms=0.727 (0.408%), neg=0, invalid=96777 0013: dt=73.984000, rms=0.724 (0.430%), neg=0, invalid=96777 0014: dt=73.984000, rms=0.721 (0.450%), neg=0, invalid=96777 0015: dt=73.984000, rms=0.718 (0.450%), neg=0, invalid=96777 0016: dt=73.984000, rms=0.714 (0.455%), neg=0, invalid=96777 0017: dt=73.984000, rms=0.711 (0.442%), neg=0, invalid=96777 0018: dt=73.984000, rms=0.708 (0.424%), neg=0, invalid=96777 0019: dt=73.984000, rms=0.705 (0.403%), neg=0, invalid=96777 0020: dt=73.984000, rms=0.703 (0.388%), neg=0, invalid=96777 0021: dt=73.984000, rms=0.700 (0.347%), neg=0, invalid=96777 0022: dt=73.984000, rms=0.698 (0.329%), neg=0, invalid=96777 thanks for the comments. Pom -- Pom

3 2

Re: [Freesurfer] problems with tracula pre-processing (bvecs)
by Michele Cavallari 04 Sep '14

04 Sep '14

Hi Barbara, thanks for your reply. I checked the dwi series on a viewer and it looks fine: I have 1960 dicoms, corresponding to 56 slices x 35 directions specified in the bval/bvec. Also, I was able to obtain the FA and MD maps through fsl using the very same dicoms and bval/bvec. Any suggestion? Thanks On Thu, Aug 28, 2014 at 7:30 PM, Barbara Kreilkamp <bakk.hbg(a)googlemail.com> wrote: Dear Michele, I cannot find anything wrong with your bvec and bval files. Seeing as you had an errormessage related to the dwi (the one about the dwi_frame.nii.gz): Did you check that your dwi data have the same amount of volumes as entries in bvec and bval? Good luck, Barbara On 28/08/2014 22:47, Michele Cavallari wrote: Hi, I'm having some problems with the first command of the tracula pipeline. I am using FreeSurfer 5.3 on a mac (OS 10.9). The problem seems to be related to the bvecs file. I received the following error message: niiRead(): error opening file /path/to/dmri/dwi_frame.nii.gz I read some threads available on your website, but couldn't figured out a solution yet. I tried to organize the bvecs in columns, as opposed to the original configuration in rows. By doing that and re-launching the pre-processing command I obtained a differente error message: Error: bvecs and bvals don't have the same number of entries The original bvals and bvecs files seem to have the same numbers of entries to me. But the bvecs file produced by the processing - both the bvecs.norot and the bvecs files in the dmri folder - don't have all the information of the vectors. I am enclosing a zip folder with attachments: 1) original bvecs and bvals files organized in rows and columns 2) bvecs and bvals files generated by tracula 3) error logs Let me know if you need any other information. Thanks in advance for your help. _______________________________________________ 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 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.

4 9

problems with tracula pre-processing (bvecs)
by Michele Cavallari 04 Sep '14

04 Sep '14

Hi, I'm having some problems with the first command of the tracula pipeline. I am using FreeSurfer 5.3 on a mac (OS 10.9). The problem seems to be related to the bvecs file. I received the following error message: niiRead(): error opening file /path/to/dmri/dwi_frame.nii.gz I read some threads available on your website, but couldn't figured out a solution yet. I tried to organize the bvecs in columns, as opposed to the original configuration in rows. By doing that and re-launching the pre-processing command I obtained a differente error message: Error: bvecs and bvals don't have the same number of entries The original bvals and bvecs files seem to have the same numbers of entries to me. But the bvecs file produced by the processing - both the bvecs.norot and the bvecs files in the dmri folder - don't have all the information of the vectors. I am enclosing a zip folder with attachments: 1) original bvecs and bvals files organized in rows and columns 2) bvecs and bvals files generated by tracula 3) error logs Let me know if you need any other information. Thanks in advance for your help.

4 3

Re: [Freesurfer] tracula: bvec-error
by Barbara Kreilkamp 04 Sep '14

04 Sep '14

Dear all, Robby, are you still having the problem? Did you solve it? I am really having the same problem, just that I am using nifti files instead of dicoms. I am using the same version as Robby. Please find my command+output, trac-all.log and configuration file in the attachment. I would also appreciate any help. Thank you, Barbara On Tue, Aug 26, 2014 at 8:21 AM, Robby De Pauw <Robby.DePauw(a)ugent.be> > wrote: > >> Hi Anastasia, >> >> I’m using this version of tracula. >> >> trac-all,v 1.22.2.12 2013/02/23 02:01:02 ayendiki Exp >> >> I’ve attached the log-file you requested. >> >> Robby De Pauw, dra. >> *Ghent University* >> Department of Physiotherapy and Rehabilitation Sciences >> 3B3 >> De Pintelaan 185 >> B-9000 Ghent >> >> Robby.DePauw(a)Ugent.be >> >> >> >> >> >> >> >> >> >> _______________________________________________ >> 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. >> >> >

3 3

watershed surface to volume rendering repost
by Sherwood, Matt 04 Sep '14

04 Sep '14

Freesurfer, Sorry for the repost, but I did not get any further details: We have scans that have recon-all completed. Following recon-all, we completed watershed (under the MNE toolbox) to extract skull and scalp from the T1 scans. The outputs of this algorithm are outer brain, inner skull, outer skull and outer scalp surfaces. We are particularly interested in these volumes for our study. Is there a way to convert these surface to volumes representative of CSF (layer between brain/skull), skull (inner to outer skull surface) and skin (outer skull to outer scalp)? Matthew Sherwood Medical Research Engineer Wright State Research Institute 4035 Colonel Glenn Hwy. Dayton, OH 45431 PH: 937.705.1025 FAX: 937.705.1095 EMAIL: matt.sherwood(a)wright.edu<mailto:matt.sherwood@wright.edu>

3 3

FreeSurfer NVIDIA card drivers with CentOS 6.5
by raij＠nmr.mgh.harvard.edu 04 Sep '14

04 Sep '14

Hi, Is there a FreeSurfer/MNE (maybe even FSL, TrackVis) recommendation to use the "noveau" drivers that come with CentOS 6.5 or would it be better to go with the proprietary drivers from NVIDIA? I would imagine that since FS 5.3 does not use GPU for analysis, this is mainly relevant for visualization (freeview etc) requiring 3D rendering of volumes, surfaces, and tractography results simultaneously. In case it matters, Dell Precision T7500 + NVIDIA Quadra FX 5800 here. Thanks! Tommi --- Tommi Raij, MD, PhD

3 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