[Freesurfer] new version of mri_segstats

2013-11-04 Thread Matura, Silke
Dear Freesurfer community, we are doing cortical thickness analyses with Freesurfer version 5.3.0 in a Linux environment. In order to do the autorecon-1 and autorecon-2 steps, we need to download the new version of mri_segstats. Unfortunately, the link on the Freesurfer Release Notes Webpage

[Freesurfer] Whole brain tractography using Freesurfer surfaces as seed and target

2013-11-04 Thread Karen Marie Sandø Ambrosen
Hi Freesurfer experts, I am trying to perform whole brain tractography using FSL's Probtrackx2 with Freesurfer white surfaces as seed and pial surfaces as targets. To obtain a whole brain connectivity graph I need the streamlines between the hemispheres. This means that I need to remove the

Re: [Freesurfer] minimal freesurfer pieces for recon-all

2013-11-04 Thread Satrajit Ghosh
hi, an update on the files needed for a vanilla recon-all run (evaluation done using cde: http://www.pgbovine.net/cde.html): freesurfer files: https://www.dropbox.com/s/a56nzr3b7b8uqz6/fslist.txt all files (fs + system+ output): https://www.dropbox.com/s/1mhcuugb20cebod/allfiles.txt size of fs

Re: [Freesurfer] A question as to subcortical volumetric measurement

2013-11-04 Thread Bruce Fischl
Hi Kim what structure is it? Our definitions are different - for example I believe we include lateral thalamic nuclei in the thalamus and they do not. Have you visually inspected both sets of segmentations? cheers Bruce On Mon, 4 Nov 2013, jh kim wrote: Dear Doug and experts I've

Re: [Freesurfer] new version of mri_segstats

2013-11-04 Thread Satrajit Ghosh
hi doug, could you please highlight the changes - especially if there are any critical bug fixes? and should we rerun the new mri_segstats on our current 5.3 processed data? cheers, satra On Mon, Nov 4, 2013 at 9:26 AM, Douglas Greve gr...@nmr.mgh.harvard.eduwrote: I've put new versions

[Freesurfer] mri_surf2vol questions

2013-11-04 Thread Govindarajan, Koushik Athreya
Hello Freesurfers, I am trying to use mri_surf2vol to convert my thickness files into analyze/nifty format to use with a custom code. Is there a way I can write all output files from mri_surf2vol in axial format? Also, is there a way in FreeSurfer to combine these thickness files

[Freesurfer] Research Felow/Postdoc: Neuroimaging in Cerebrovascular Disease

2013-11-04 Thread Direktion, Institut f. Schlaganfall- und Demenzforschung
The Hospital of the University of Munich, Germany, is one of the largest and most competitive university hospitals in Germany and Europe. 45 specialized hospitals, departments and institutions harboring excellent research and education provide patient care at the highest medical level with

Re: [Freesurfer] mri_segment WM/GM limits

2013-11-04 Thread Louis Nicholas Vinke
Hi Jon, The limit numbers correspond to voxel intensities in the brainmask.mgz. I believe you will want a lower number for the -wlo/-wm_low flag to make the wm surf extend further out, assuming that the wm voxels that are not being captured by the surface is hypo-intense ( 110). -Louis On

[Freesurfer] expand GM label into white matter

2013-11-04 Thread Corinna Bauer
Hello all, Is there a way to expand a GM label into the white matter by n voxels? Corinna ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer The information in this e-mail is

Re: [Freesurfer] expand GM label into white matter

2013-11-04 Thread Douglas N Greve
There is not a single command to do it, but you could do something like mri_binarize --i aparc+aseg.mgz --wm --o wmmask.mgz # mask of white matter mri_binarize --i seg.mgz --match 17 --o mask17.mgz # mask of seg 17 mri_binarize --i seg.mgz --match 17 --dilate 2 --o mask17.dil2.mgz # mask of

Re: [Freesurfer] Inaccurate white matter and pial surfaces

2013-11-04 Thread Bruce Fischl
Hi Linda try using the expert options to specify the gray and white matter intensity ranges cheers Bruce On Mon, 4 Nov 2013, Linda Zhang wrote: Dear all, A few of my subjects end up with surfaces where the white matter surface is actually the pial surface, probably due to less than ideal

Re: [Freesurfer] expand GM label into white matter

2013-11-04 Thread Bruce Fischl
Hi Corinna I think this what Doug's WM parcellation tool does cheers Bruce On Mon, 4 Nov 2013, Corinna Bauer wrote: Hello all, Is there a way to expand a GM label into the white matter by n voxels? Corinna ___ Freesurfer mailing list

Re: [Freesurfer] mri_segment WM/GM limits

2013-11-04 Thread Jonathan Holt
Hey Louis, thanks for your response. It’s actually a bit different, in my case an entire sulcus, including the empty space within it has been segmented and surfaced as white matter. I’m trying to reduce the WM surface such that the empty space is not included. I’ve tried removing segmentation

Re: [Freesurfer] expand GM label into white matter

2013-11-04 Thread Corinna Bauer
ok thanks, I'll give that a try. I would have to do that for each of the segmentations of interest and I get the number from the LUT? On Mon, Nov 4, 2013 at 11:47 AM, Douglas N Greve gr...@nmr.mgh.harvard.eduwrote: There is not a single command to do it, but you could do something like

Re: [Freesurfer] mri_segment WM/GM limits

2013-11-04 Thread Bruce Fischl
Hi Jon maybe you should upload your data and we'll take a look cheers Bruce On Mon, 4 Nov 2013, Jonathan Holt wrote: Bruce, If I understand you correctly, you?re saying I should tweak my expert options? Should I further change the limits, how so? You previously advised that I not adjust

Re: [Freesurfer] mri_segment WM/GM limits

2013-11-04 Thread Jonathan Holt
OK, been uploaded as JONbrain.tar Please let me know! On Nov 4, 2013, at 1:15 PM, Bruce Fischl fis...@nmr.mgh.harvard.edu wrote: Hi Jon maybe you should upload your data and we'll take a look cheers Bruce On Mon, 4 Nov 2013, Jonathan Holt wrote: Bruce, If I understand you

Re: [Freesurfer] mri_segment WM/GM limits

2013-11-04 Thread Louis Nicholas Vinke
Got it, I'll take a look. On Mon, 4 Nov 2013, Jonathan Holt wrote: OK, been uploaded as JONbrain.tar Please let me know! On Nov 4, 2013, at 1:15 PM, Bruce Fischl fis...@nmr.mgh.harvard.edu wrote: Hi Jon maybe you should upload your data and we'll take a look cheers Bruce On Mon,

[Freesurfer] convert (resample) fsaverage/mri/aseg.mgz into $FSLDIR Nifti with MNI152_T1_1mm target

2013-11-04 Thread Roderick McColl
I am obviously using mri_vol2vol incorrectly but I'm stumped. Goal is to get the $SUBJECTS/fsaverage/mri/aseg.mgz file into the FSL MNI152 space to do some fslmaths. From looking around it seemed the command to accomplish this is (inside $SUBJECTS/fsaverage/mri) mri_vol2vol --targ

Re: [Freesurfer] convert (resample) fsaverage/mri/aseg.mgz into $FSLDIR Nifti with MNI152_T1_1mm target

2013-11-04 Thread Bruce Fischl
Hi Roddy can you check whether fslview shows the directions as correct even if it displays it upside down? That is, do they still know which direction is superior, anterior, etc...? cheers Bruce On Mon, 4 Nov 2013, Roderick McColl wrote: I am obviously using mri_vol2vol incorrectly but

[Freesurfer] Which boost version is used in CVS, mris_resample

2013-11-04 Thread Martin Kavec
Hi guys, When running mri_cvs_register on linux (5.3.0), I get into problem with mris_resample, which cannot find libboost_programs_options.so.5. I installed the latest version of boost-1.52.0 for my system, but there are unresolved symbols. Which version is suitable for mris_resample?

Re: [Freesurfer] expand GM label into white matter

2013-11-04 Thread Douglas N Greve
Yes. As Bruce mentioned you can also use the program that creates the wmparc (mri_aparc2aseg). You can look in the recon-all.log file for the command line (there will be two in there, one to create aparc+aseg and one to create wmparc). You can run this command with --wmparc-dmax to set the

Re: [Freesurfer] new version of mri_segstats

2013-11-04 Thread Satrajit Ghosh
thanks doug. cheers, satra On Mon, Nov 4, 2013 at 5:33 PM, Douglas N Greve gr...@nmr.mgh.harvard.eduwrote: It gives the same results as the 5.3 version. The difference is that it does not require the ribbon.mgz file to be there in order to run. If you run recon-all from start to finish,

Re: [Freesurfer] convert (resample) fsaverage/mri/aseg.mgz into $FSLDIR Nifti with MNI152_T1_1mm target

2013-11-04 Thread Douglas N Greve
try this command instead mri_vol2vol --mov $FSLDIR/data/standard/MNI152_T1_1mm.nii.gz \ --targ aseg.mgz --inv --reg transforms/reg.mni152.1mm.dat --o asegmni152.nii.gz On 11/04/2013 04:37 PM, Roderick McColl wrote: I am obviously using mri_vol2vol incorrectly but I'm stumped. Goal is

Re: [Freesurfer] mri_surf2vol questions

2013-11-04 Thread Douglas N Greve
To get them both in one file, run mri_surf2vol to get the left hemisphere. When you run it a 2nd time for the right hemisphere add --merge vol.lh.nii. I think FSL has software to swap dimensions doug On 11/04/2013 09:53 AM, Govindarajan, Koushik Athreya wrote: Hello Freesurfers,

Re: [Freesurfer] tkregister2 and mri_label2vol

2013-11-04 Thread Douglas N Greve
* try running this command mri_label2vol --label lh.words10_10.label --temp rawavg.mgz --reg rawreg.dat --fillthresh .5 --o lh.words10_10.nii * On 11/04/2013 07:55 AM, Zhao Jingjing wrote: Hi Doug, I asked a question a couple days ago for how to convert .label files into .nii files in

Re: [Freesurfer] convert (resample) fsaverage/mri/aseg.mgz into $FSLDIR Nifti with MNI152_T1_1mm target - aparc+aseg.mgz problem

2013-11-04 Thread Roderick McColl
Hmmm, Actually using aseg.mgz works just fine (just tried this now), displays correctly on fslview. What I was really trying to do was convert aparc+aseg.mgz which does not provide the expected orientation in fslview. Does this mean aparc+aseg.mgz and aseg.mgz are in different spaces somehow?

Re: [Freesurfer] A question as to subcortical volumetric measurement

2013-11-04 Thread Douglas N Greve
Hi Kim, the answer is both easy and difficult. The easy answer is that they are completely different methods and will often give inconsistent results. I'm not sure what answer the reviewer will accept, and you might need to post the same question to the VBM and DARTEL crews. In your VBM

Re: [Freesurfer] Error when using --mask flag in mri_glmfit

2013-11-04 Thread Douglas N Greve
Hi Suzanne, there is probably an easier way to do this. If you create contrasts of each condition vs baseline then run isxconcat-sess on each, you can then run something like mri_segstats --seg mc-z.pos.23.lh.sig.ocn.mgh --i condition1.nii.gz --avgwf condition1.table.dat --excludeid 0 where

Re: [Freesurfer] mri_segstats url in 5.3.0 release notes does not exist

2013-11-04 Thread Douglas N Greve
Hi Jonathan, see my email from earlier about the new permanent location for mri_segstats doug On 10/30/2013 12:26 PM, Jonathan Harlap wrote: The version of mri_segstats that fixes the issue where mri_segstats fails unless -autorecon3 is run returns a file not found at the url

Re: [Freesurfer] Fwd: Qdec Questions

2013-11-04 Thread Douglas N Greve
Yes, true for all models. I'm not sure what I'm looking at in this scatter plot. The scatter plot does not have the nuisance effects removed, so it might be that when you remove the nuisance effect, then HCSAD doug On 10/30/2013 01:32 PM, Ashley Shurick wrote: Great - and this would be

Re: [Freesurfer] How do you create labels from binary masks?

2013-11-04 Thread Douglas N Greve
That vol2surf command might miss the values you want to sample. Try adding --projdist-max -2 0 .1 This will sample from 2mm into the hippo out to the surface in .1mm increments and then take the maximum value it finds doug On 10/30/2013 03:08 PM, Marcos Martins da Silva wrote: Hi, Doug.

Re: [Freesurfer] glmfit for paired-diff on 2 groups

2013-11-04 Thread Douglas N Greve
Yes, that looks correct, though I could be more confident if you send the command lines used and your fsgd file(s). doug On 10/31/2013 07:32 AM, Vivian R. Steiger wrote: Dear FS-Experts, I would like to compare 2 Groups (Responder (N=19) vs. Non-Responder (N=11)) with longitudinal data.

Re: [Freesurfer] mri_segstats url in 5.3.0 release notes does not exist

2013-11-04 Thread Jonathan Harlap
Thanks Doug, I did see that. Cheers, Jonathan On Nov 4, 2013, at 6:23 PM, Douglas N Greve gr...@nmr.mgh.harvard.edu wrote: Hi Jonathan, see my email from earlier about the new permanent location for mri_segstats doug On 10/30/2013 12:26 PM, Jonathan Harlap wrote: The version of

Re: [Freesurfer] convert (resample) fsaverage/mri/aseg.mgz into $FSLDIR Nifti with MNI152_T1_1mm target - aparc+aseg.mgz problem

2013-11-04 Thread Douglas N Greve
Hi Roddy, aparc+aseg.mgz and aseg.mgz are in the exact same space and the same command line should generate the same result. Does it display properly in tkmedit? If so and if you're having problems visualizing with fslview, contact the FSL crew. doug On 11/04/2013 05:52 PM, Roderick McColl

[Freesurfer] Seg fault with mri_segstats?

2013-11-04 Thread Michael Waskom
Hi Doug, I'm running into a strange seg fault when trying to use mri_segstats (it's actually ocurring within mri_glmfit-sim, but I can reproduce it manually) This is the line that fails (stripping the leading paths for readability): mri_segstats --seg cache.th20.pos.sig.ocn.mgh --exclude 0 --i

[Freesurfer] Creating group average ROI

2013-11-04 Thread Miggy Chuapoco
Hi all, We have ROIs (defined in mrVista) that we converted to label files and would like to now make an average ROI. Does a pipeline exist that will help us create an average ROI composed of these label files. We were thinking of using mri_label2label to register the ROIs to a common surface but

Re: [Freesurfer] Creating group average ROI

2013-11-04 Thread Bruce Fischl
Hi Miggy Are the ROIs on different subjects? Do you have surfaces reconstructed for each of them? cheers Bruce On Mon, 4 Nov 2013, Miggy Chuapoco wrote: Hi all,  We have ROIs (defined in mrVista) that we converted to label files and would like to now make an average ROI. Does a pipeline

Re: [Freesurfer] Creating group average ROI

2013-11-04 Thread Douglas Greve
I would convert them to a binary mask (0,1) then average the masks together. Each vertex would then be the probability that a given vertex was in the ROI. You can convert the label to a binary mask with a second call to mri_label2label (--outmask option) doug On 11/4/13 8:47 PM, Miggy