Re: [Freesurfer] convert DK atlas ROI to .nii.gz in MNI space

2013-07-17 Thread Yang, Daniel
Hi Doug,

Thanks very much!!

Best,
Daniel

On Jul 16, 2013, at 9:14 PM, Douglas Greve 
gr...@nmr.mgh.harvard.edumailto:gr...@nmr.mgh.harvard.edu wrote:


use

$FREESURFER_HOME/average/mni152.register.dat

doug

On 7/16/13 8:36 PM, Yang, Daniel wrote:
Dear FreeSurfer Experts,

I am interested in a specific ROI in the DK atlas, and I know that I can use 
tkregister2 to obtain a register.dat, and then mri_annotation2label + 
mri_label2vol to convert the annotation of the ROI to an .nii.gz.

However, it seems that the end product is in the TAL space? If yes, how can I 
convert it to the MNI 152 space?

Many thanks!!

Best,
Daniel
--
Yung-Jui Daniel Yang, PhD
Postdoctoral Researcher
Yale Child Study Center
New Haven, CT
(203) 737-5454



___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edumailto:Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edumailto:Freesurfer@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@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.


Re: [Freesurfer] bvecs error from Tracula

2013-07-17 Thread amirhossein manzouri
Thanks Anastasia,
IT HELPED!


On Tue, Jul 16, 2013 at 5:25 PM, Anastasia Yendiki 
ayend...@nmr.mgh.harvard.edu wrote:


 Hi Amirhossein - I'm attaching sample files.

 a.y


 On Tue, 16 Jul 2013, amirhossein manzouri wrote:

  Hi Anastasia,
 I have tried both , 3-row and 3-column format and still wrong bvecc and
 bval
 from flip4fsl. Could you please send me a sample of your original files
 so I
 can compare?


 On Mon, Jul 15, 2013 at 7:44 PM, Anastasia Yendiki
 ayend...@nmr.mgh.harvard.edu wrote:

   Hi Amirhossein - Are your original bvecs/bvals in 3-row format
   instead of 3-column format?
   
 http://surfer.nmr.mgh.harvard.**edu/fswiki/FsTutorial/Traculahttp://surfer.nmr.mgh.harvard.edu/fswiki/FsTutorial/Tracula

   Thanks,
   a.y

   On Mon, 15 Jul 2013, amirhossein manzouri wrote:

 Dear Experts,
 I am running trac-all -prep -c dmrirc on my DWI
 data. After flip4fsl step I get the attached bval
 and bvec file which the bvec one is wrong so the
 process exits
 with error in dtifit. I have also attached  original
 bvec and bval.

 --
 Best regards,
 Amirhossein Manzouri







 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/**compliancelinehttp://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.




 --
 Best regards,
 Amirhossein Manzouri







-- 
Best regards,
Amirhossein Manzouri
___
Freesurfer mailing list
Freesurfer@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] Missing flirt.fsl file in bin folder

2013-07-17 Thread Gundran, Andrew
Hello freesurfers,

I’ve been trying to run tracula on freesurfer and it seems that when I get to 
the point where it uses FSL’s flirt I receive an error. I looked into it 
further and it seems I am missing a file to run flirt:

$FREESURFER_HOME/bin/flirt.fsl

It seems that I somehow do not have that file to complete tracula.  I was 
hoping someone could help with this or perhaps send me the file?


I am using freesurfer-
Darwin-snowleopard-i686-stable-pub-v5.3.0-
MacBook-Pro.local 10.8.0 Darwin Kernel Version 10.8.0:

Hope to hear from someone soon,
Andrew 
___
Freesurfer mailing list
Freesurfer@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] Propagating changes to segmentation

2013-07-17 Thread Liane Hunter
Hello,

I have a question about propagating changes to segmentation. Once I have
saved my edits what is the appropriate command to propagate these changes
for both cortical and subcortical changes that I have made.

The workflow instructs that after segmentation edits are made to run

 -normalization -maskbfs -segmentation -autorecon2-wm -autorecon3 to the wm.
mgz file. However the instructions indicate to open brainmask.mgz to make
edits.  Should I make changes on the wm.mgz file?

Any clarification on this topic would be much appreciated.

-- 
Liane Hunter
MD/Ph.D Candidate
Albert Einstein College of Medicine
___
Freesurfer mailing list
Freesurfer@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] Freesurfer installation

2013-07-17 Thread pablo najt
Dear Freesurfer experts,
 I have a question regarding the installation/running of FreeSurfer.
I have been trying to install freesurfer in a macbook pro (snow leopard) but no 
success. 
I followed the webpage from freesurfer (link below). 


https://surfer.nmr.mgh.harvard.edu/fswiki/FsTutorial/PersonalLaptopSetup


But after the installations when I open the terminal nothing happens.


I also tried to launch freesurfer this way:
export FREESURFER_HOME=/Applications/freesurfer
source $FREESURFER_HOME/SetUpFreeSurfer.sh
If I do this I get:
command not found
command not found
command not found


Any advice on how to finally set up the software would be great help.
Thank you in anticipation.Pablo   ___
Freesurfer mailing list
Freesurfer@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] kvlQuantifyHippocampalSubfieldSegmentations.sh empty text files

2013-07-17 Thread Jasmeet Hayes
Dear Freesurfer experts,

I tried running the kvlQuantifyHippocampalSubfieldSegmentations.sh script
and get two text files (left  right stats) but the only information in the
text files is the subject ID number.

However, on my terminal window, I can see the volume output for each
subfield and for each subject (see below for an example of 1 subject). And
there were no errors after running the script.

Every subject in the folder went through the hippocampal subfield recon
procedure, and I have verified that there were no errors (each subject's
folder has the left and right posterior_* files).

We made a slight modification to the kvl script, only to change the Subject
directory. Any ideas about what might be the problem? Thanks! -Jasmeet

cd .. Quantifying subject er002_recon/ right Doing right side cd mri THIS
IS THE COMPRESSION LOOKUP TABLE FILE
/Applications/freesurfer/data/GEMS/compressionLookupTable.txt PWD before
running KVLQuantifyPosteriorProb IS:
/Volumes/VA_Imaging/Projects/jpannuhayes/emoreg/HIPPOCAMPAL_SUBFIELDS/er002_recon/mri
kvlQuantifyPosteriorProbabilityImages
/Applications/freesurfer/data/GEMS/compressionLookupTable.txt
posterior_Right-Hippocampus.mgz posterior_right_presubiculum.mgz
posterior_right_CA1.mgz posterior_right_CA2-3.mgz
posterior_right_fimbria.mgz posterior_right_subiculum.mgz
posterior_right_CA4-DG.mgz posterior_right_hippocampal_fissure.mgz 
volumeStats_right.txt volumeInVoxels: Right-Hippocampus: 3094.49
right_presubiculum: 4373.15 right_CA1: 2483.23 right_CA2-3: 8192.11
right_fimbria: 740.934 right_subiculum: 5402.79 right_CA4-DG: 4688.93
right_hippocampal_fissure: 541.872 cd .. cd ..
___
Freesurfer mailing list
Freesurfer@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.


Re: [Freesurfer] Freesurfer installation

2013-07-17 Thread Z K
Pablo,

That page describes how to setup a personal laptop for usage in the 
Freesurfer course. It does not apply to a standard freesurfer 
installation.

Try opening a new terminal and typing bash before the export and 
source commands you give below.

-Zeke


On 07/17/2013 10:27 AM, pablo najt wrote:
 Dear Freesurfer experts,

   I have a question regarding the installation/running of FreeSurfer.

 I have been trying to install freesurfer in a macbook pro (snow leopard)
 but no success.

 I followed the webpage from freesurfer (link below).


 https://surfer.nmr.mgh.harvard.edu/fswiki/FsTutorial/PersonalLaptopSetup


 But after the installations when I open the terminal nothing happens.


 I also tried to launch freesurfer this way:

 export FREESURFER_HOME=/Applications/freesurfer

 source $FREESURFER_HOME/SetUpFreeSurfer.sh

 If I do this I get:

 command not found

 command not found

 command not found


 Any advice on how to finally set up the software would be great help.

 Thank you in anticipation.

 Pablo



 ___
 Freesurfer mailing list
 Freesurfer@nmr.mgh.harvard.edu
 https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

___
Freesurfer mailing list
Freesurfer@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.



Re: [Freesurfer] kvlQuantifyHippocampalSubfieldSegmentations.sh empty text files

2013-07-17 Thread Juan Eugenio Iglesias

Hi Jasmeet,
what modification did you exactly make? Whih version of FreeSurfer is this?
Cheers,
/Eugenio

On 07/17/2013 10:22 AM, Jasmeet Hayes wrote:

Dear Freesurfer experts,

I tried running the kvlQuantifyHippocampalSubfieldSegmentations.sh 
script and get two text files (left  right stats) but the only 
information in the text files is the subject ID number.


However, on my terminal window, I can see the volume output for each 
subfield and for each subject (see below for an example of 1 subject). 
And there were no errors after running the script.


Every subject in the folder went through the hippocampal subfield 
recon procedure, and I have verified that there were no errors (each 
subject's folder has the left and right posterior_* files).


We made a slight modification to the kvl script, only to change the 
Subject directory. Any ideas about what might be the problem? Thanks! 
-Jasmeet


cd .. Quantifying subject er002_recon/ right Doing right side cd mri 
THIS IS THE COMPRESSION LOOKUP TABLE FILE 
/Applications/freesurfer/data/GEMS/compressionLookupTable.txt PWD 
before running KVLQuantifyPosteriorProb IS: 
/Volumes/VA_Imaging/Projects/jpannuhayes/emoreg/HIPPOCAMPAL_SUBFIELDS/er002_recon/mri 
kvlQuantifyPosteriorProbabilityImages 
/Applications/freesurfer/data/GEMS/compressionLookupTable.txt 
posterior_Right-Hippocampus.mgz posterior_right_presubiculum.mgz 
posterior_right_CA1.mgz posterior_right_CA2-3.mgz 
posterior_right_fimbria.mgz posterior_right_subiculum.mgz 
posterior_right_CA4-DG.mgz posterior_right_hippocampal_fissure.mgz  
volumeStats_right.txt volumeInVoxels: Right-Hippocampus: 3094.49 
right_presubiculum: 4373.15 right_CA1: 2483.23 right_CA2-3: 8192.11 
right_fimbria: 740.934 right_subiculum: 5402.79 right_CA4-DG: 4688.93 
right_hippocampal_fissure: 541.872 cd .. cd ..



___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer


--
-
Juan Eugenio Iglesias, PhD
http://www.jeiglesias.com
igles...@nmr.mgh.harvard.edu
Athinoula A. Martinos Center for Biomedical Imaging
149 Thirteenth Street, Suite 2301
Charlestown, Massachusetts 2129
U.S.A.

___
Freesurfer mailing list
Freesurfer@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.


Re: [Freesurfer] bvecs error from Tracula

2013-07-17 Thread Anastasia Yendiki


Good to hear! Can you please share with all of us what was wrong with your 
original files? It may help others who have the same problem in the 
future. Thanks!


On Wed, 17 Jul 2013, amirhossein manzouri wrote:


Thanks Anastasia, 
IT HELPED!


On Tue, Jul 16, 2013 at 5:25 PM, Anastasia Yendiki 
ayend...@nmr.mgh.harvard.edu
wrote:

  Hi Amirhossein - I'm attaching sample files.

  a.y

  On Tue, 16 Jul 2013, amirhossein manzouri wrote:

Hi Anastasia, 
I have tried both , 3-row and 3-column format and still
wrong bvecc and bval
from flip4fsl. Could you please send me a sample of your
original files so I
can compare?


On Mon, Jul 15, 2013 at 7:44 PM, Anastasia Yendiki
ayend...@nmr.mgh.harvard.edu wrote:

      Hi Amirhossein - Are your original bvecs/bvals in
3-row format
      instead of 3-column format?
     
http://surfer.nmr.mgh.harvard.edu/fswiki/FsTutorial/Tracula

      Thanks,
      a.y

      On Mon, 15 Jul 2013, amirhossein manzouri wrote:

            Dear Experts, 
            I am running trac-all -prep -c dmrirc on my
DWI
            data. After flip4fsl step I get the attached
bval
            and bvec file which the bvec one is wrong so
the
            process exits
            with error in dtifit. I have also attached
 original
            bvec and bval.

            --
            Best regards,
            Amirhossein Manzouri 

             





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.




--
Best regards,
Amirhossein Manzouri 

 





--
Best regards,
Amirhossein Manzouri 

 


___
Freesurfer mailing list
Freesurfer@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.


Re: [Freesurfer] kvlQuantifyHippocampalSubfieldSegmentations.sh empty text files

2013-07-17 Thread Jasmeet Hayes
Hi Eugenio,

We are using version 5.1. I attached our script. We added to line 60,
setting the subject directory. We also changed the statFileName (see line
291 on the attached script).  I should also note that the volumeStats_left
and right text files within each subject's mri folder is empty. Thanks for
your help.


On Wed, Jul 17, 2013 at 11:16 AM, Juan Eugenio Iglesias 
igles...@nmr.mgh.harvard.edu wrote:

  Hi Jasmeet,
 what modification did you exactly make? Whih version of FreeSurfer is this?
 Cheers,
 /Eugenio


 On 07/17/2013 10:22 AM, Jasmeet Hayes wrote:

 Dear Freesurfer experts,

  I tried running the kvlQuantifyHippocampalSubfieldSegmentations.sh
 script and get two text files (left  right stats) but the only information
 in the text files is the subject ID number.

  However, on my terminal window, I can see the volume output for each
 subfield and for each subject (see below for an example of 1 subject). And
 there were no errors after running the script.

  Every subject in the folder went through the hippocampal subfield recon
 procedure, and I have verified that there were no errors (each subject's
 folder has the left and right posterior_* files).

  We made a slight modification to the kvl script, only to change the
 Subject directory. Any ideas about what might be the problem? Thanks!
 -Jasmeet

  cd .. Quantifying subject er002_recon/ right Doing right side cd mri
 THIS IS THE COMPRESSION LOOKUP TABLE FILE
 /Applications/freesurfer/data/GEMS/compressionLookupTable.txt PWD before
 running KVLQuantifyPosteriorProb IS:
 /Volumes/VA_Imaging/Projects/jpannuhayes/emoreg/HIPPOCAMPAL_SUBFIELDS/er002_recon/mri
 kvlQuantifyPosteriorProbabilityImages
 /Applications/freesurfer/data/GEMS/compressionLookupTable.txt
 posterior_Right-Hippocampus.mgz posterior_right_presubiculum.mgz
 posterior_right_CA1.mgz posterior_right_CA2-3.mgz
 posterior_right_fimbria.mgz posterior_right_subiculum.mgz
 posterior_right_CA4-DG.mgz posterior_right_hippocampal_fissure.mgz 
 volumeStats_right.txt volumeInVoxels: Right-Hippocampus: 3094.49
 right_presubiculum: 4373.15 right_CA1: 2483.23 right_CA2-3: 8192.11
 right_fimbria: 740.934 right_subiculum: 5402.79 right_CA4-DG: 4688.93
 right_hippocampal_fissure: 541.872 cd .. cd ..


 ___
 Freesurfer mailing 
 listfreesur...@nmr.mgh.harvard.eduhttps://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer


 --
 -
 Juan Eugenio Iglesias, 
 PhDhttp://www.jeiglesias.comigles...@nmr.mgh.harvard.edu
 Athinoula A. Martinos Center for Biomedical Imaging
 149 Thirteenth Street, Suite 2301
 Charlestown, Massachusetts 2129
 U.S.A.

  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.



kvlQuantifyHippocampalSubfieldSegmentations.sh
Description: Bourne shell script
___
Freesurfer mailing list
Freesurfer@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.


Re: [Freesurfer] SWI acquisition question

2013-07-17 Thread Bruce Fischl
Hi Panos

sorry, I think this is pretty far outside our expertise, and it would be 
hard in any case to try to figure out what happened to give you these 
strange looking results. I certainly don't have any insight

sorry,
Bruce
On Wed, 17 Jul 
2013, Fotiadis, Panagiotis wrote:

 Hi FS Community,

 I was wondering whether anyone had a chance to look into this. Thanks again 
 for your time!

 Best,
 Panos
 
 From: freesurfer-boun...@nmr.mgh.harvard.edu 
 [freesurfer-boun...@nmr.mgh.harvard.edu] on behalf of Fotiadis, Panagiotis
 Sent: Tuesday, July 16, 2013 2:00 PM
 To: freesurfer@nmr.mgh.harvard.edu
 Subject: [Freesurfer] SWI acquisition question

 Hi FreeSurfer community,

 I am analyzing some SWI scans of some old subjects. For most of my subjects 
 when a SWI was acquired, the resulting files that I would get would be a 
 Magnitude image, a Phase image, a mIP image, and the SWI image. However, I 
 have a couple of subjects for whom I get only the mIP and the SWI images, and 
 the SWI image looks more blurry than it's supposed to be. I was wondering:

 1) Why do I get only the mIP and the SWI images for those subjects and not 
 the magnitude and phase images as well,
 2) Is it normal for the mIP to look like that? Usually it looks different, and
 3) Why is the final SWI image more blurry than usual? By looking at the 
 image, I believe that maybe the head coil was not properly connected to the 
 scanner at the time of the acquisition but I don't know whether that would be 
 a sufficient explanation.

 I have attached two attachments: One of the resulting mIP and one of the 
 equivalent blurry SWI of the same subject.

 Thank you in advance for your help and time!

 Best,
 Panos

___
Freesurfer mailing list
Freesurfer@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.



Re: [Freesurfer] kvlQuantifyHippocampalSubfieldSegmentations.sh empty text files

2013-07-17 Thread Juan Eugenio Iglesias
OK, then the question is why kvlQuantifyPosteriorProbabilityImages is 
producing an empty file
I think the problem is that 
kvlQuantifyHippocampalSubfieldSegmentations.sh behaves in a different 
way depending on the number arguments, etc. I would suggest that you go 
back to the original version of the script, and use it the way it's 
described in 
http://ftp.nmr.mgh.harvard.edu/fswiki/HippocampalSubfieldSegmentation, 
that is, setting the environment variable SUBJECTS_DIR and calling the 
script without arguments.

Kind regards,
/Eugenio



On 07/17/2013 11:46 AM, Jasmeet Hayes wrote:

Hi Eugenio,

We are using version 5.1. I attached our script. We added to line 60, 
setting the subject directory. We also changed the statFileName (see 
line 291 on the attached script).  I should also note that the 
volumeStats_left and right text files within each subject's mri folder 
is empty. Thanks for your help.



On Wed, Jul 17, 2013 at 11:16 AM, Juan Eugenio Iglesias 
igles...@nmr.mgh.harvard.edu mailto:igles...@nmr.mgh.harvard.edu 
wrote:


Hi Jasmeet,
what modification did you exactly make? Whih version of FreeSurfer
is this?
Cheers,
/Eugenio


On 07/17/2013 10:22 AM, Jasmeet Hayes wrote:

Dear Freesurfer experts,

I tried running the
kvlQuantifyHippocampalSubfieldSegmentations.sh script and get two
text files (left  right stats) but the only information in the
text files is the subject ID number.

However, on my terminal window, I can see the volume output for
each subfield and for each subject (see below for an example of 1
subject). And there were no errors after running the script.

Every subject in the folder went through the hippocampal subfield
recon procedure, and I have verified that there were no errors
(each subject's folder has the left and right posterior_* files).

We made a slight modification to the kvl script, only to change
the Subject directory. Any ideas about what might be the problem?
Thanks! -Jasmeet

cd .. Quantifying subject er002_recon/ right Doing right side cd
mri THIS IS THE COMPRESSION LOOKUP TABLE FILE
/Applications/freesurfer/data/GEMS/compressionLookupTable.txt PWD
before running KVLQuantifyPosteriorProb IS:

/Volumes/VA_Imaging/Projects/jpannuhayes/emoreg/HIPPOCAMPAL_SUBFIELDS/er002_recon/mri
kvlQuantifyPosteriorProbabilityImages
/Applications/freesurfer/data/GEMS/compressionLookupTable.txt
posterior_Right-Hippocampus.mgz posterior_right_presubiculum.mgz
posterior_right_CA1.mgz posterior_right_CA2-3.mgz
posterior_right_fimbria.mgz posterior_right_subiculum.mgz
posterior_right_CA4-DG.mgz
posterior_right_hippocampal_fissure.mgz  volumeStats_right.txt
volumeInVoxels: Right-Hippocampus: 3094.49 right_presubiculum:
4373.15 right_CA1: 2483.23 right_CA2-3: 8192.11 right_fimbria:
740.934 right_subiculum: 5402.79 right_CA4-DG: 4688.93
right_hippocampal_fissure: 541.872 cd .. cd ..


___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu  mailto:Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer


-- 
-

Juan Eugenio Iglesias, PhD
http://www.jeiglesias.com
igles...@nmr.mgh.harvard.edu  mailto:igles...@nmr.mgh.harvard.edu
Athinoula A. Martinos Center for Biomedical Imaging
149 Thirteenth Street, Suite 2301
Charlestown, Massachusetts 2129
U.S.A.

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.




--
-
Juan Eugenio Iglesias, PhD
http://www.jeiglesias.com
igles...@nmr.mgh.harvard.edu
Athinoula A. Martinos Center for Biomedical Imaging
149 Thirteenth Street, Suite 2301
Charlestown, Massachusetts 2129
U.S.A.

___
Freesurfer mailing list
Freesurfer@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.


Re: [Freesurfer] convert DK atlas ROI to .nii.gz in MNI space

2013-07-17 Thread Yang, Daniel
Hi Doug,

Can I ask a basic question about mri_label2vol? in mri_label2vol, what should 
be the input for --temp? I understand it's supposed to be the template image. 
If it's the native space, should it be rawavg.mgz or orig.mgz for that 
individual?

Great thanks!

Best,
Daniel

--
Yung-Jui Daniel Yang, PhD
Postdoctoral Researcher
Yale Child Study Center
New Haven, CT
(203) 737-5454

On 7/16/13 9:13 PM, Douglas Greve 
gr...@nmr.mgh.harvard.edumailto:gr...@nmr.mgh.harvard.edu wrote:


use

$FREESURFER_HOME/average/mni152.register.dat

doug

On 7/16/13 8:36 PM, Yang, Daniel wrote:
Dear FreeSurfer Experts,

I am interested in a specific ROI in the DK atlas, and I know that I can use 
tkregister2 to obtain a register.dat, and then mri_annotation2label + 
mri_label2vol to convert the annotation of the ROI to an .nii.gz.

However, it seems that the end product is in the TAL space? If yes, how can I 
convert it to the MNI 152 space?

Many thanks!!

Best,
Daniel
--
Yung-Jui Daniel Yang, PhD
Postdoctoral Researcher
Yale Child Study Center
New Haven, CT
(203) 737-5454



___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edumailto:Freesurfer@nmr.mgh.harvard.eduhttps://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

___
Freesurfer mailing list
Freesurfer@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.


Re: [Freesurfer] convert DK atlas ROI to .nii.gz in MNI space

2013-07-17 Thread Douglas N Greve
does htis answer your question? It is from running mri_label2vol with --help
doug

--temp tempvolid

Template volume. The output volume will have the same size and geometry
as the template. Template must have geometry information (ie, direction
cosines and voxel sizes). Required.



On 07/17/2013 12:11 PM, Yang, Daniel wrote:
 Hi Doug,

 Can I ask a basic question about mri_label2vol? in mri_label2vol, what 
 should be the input for --temp? I understand it's supposed to be the 
 template image. If it's the native space, should it be rawavg.mgz or 
 orig.mgz for that individual?

 Great thanks!

 Best,
 Daniel

 -- 
 Yung-Jui Daniel Yang, PhD
 Postdoctoral Researcher
 Yale Child Study Center
 New Haven, CT
 (203) 737-5454

 On 7/16/13 9:13 PM, Douglas Greve gr...@nmr.mgh.harvard.edu 
 mailto:gr...@nmr.mgh.harvard.edu wrote:


 use

 $FREESURFER_HOME/average/mni152.register.dat

 doug

 On 7/16/13 8:36 PM, Yang, Daniel wrote:
 Dear FreeSurfer Experts,

 I am interested in a specific ROI in the DK atlas, and I know
 that I can use tkregister2 to obtain a register.dat, and then
 mri_annotation2label + mri_label2vol to convert the annotation of
 the ROI to an .nii.gz.

 However, it seems that the end product is in the TAL space? If
 yes, how can I convert it to the MNI 152 space?

 Many thanks!!

 Best,
 Daniel
 -- 
 Yung-Jui Daniel Yang, PhD
 Postdoctoral Researcher
 Yale Child Study Center
 New Haven, CT
 (203) 737-5454


 ___
 Freesurfer mailing list
 
 freesur...@nmr.mgh.harvard.eduhttps://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer


-- 
Douglas N. Greve, Ph.D.
MGH-NMR Center
gr...@nmr.mgh.harvard.edu
Phone Number: 617-724-2358
Fax: 617-726-7422

Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting
FileDrop: https://gate.nmr.mgh.harvard.edu/filedrop2
www.nmr.mgh.harvard.edu/facility/filedrop/index.html
Outgoing: ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/

___
Freesurfer mailing list
Freesurfer@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.



Re: [Freesurfer] Missing flirt.fsl file in bin folder

2013-07-17 Thread Z K
Hello,

Individuals who downloaded Freesurfer v5.3 for snow leopard (32bit) will 
be missing six fsl routines. A patch has been created and can be 
downloaded from this link:

ftp://surfer.nmr.mgh.harvard.edu/pub/dist/freesurfer/patches/freesurfer_fsl.tar.gz

To apply the patch, simply extract the contents of the above file into 
the /Applications/freesurfer/bin directory. This should result in six 
new files in the /Applications/freesurfer/bin directory all with the 
.fsl extension.

-Zeke



On 07/17/2013 09:33 AM, Gundran, Andrew wrote:
 Hello freesurfers,

 I’ve been trying to run tracula on freesurfer and it seems that when I get to 
 the point where it uses FSL’s flirt I receive an error. I looked into it 
 further and it seems I am missing a file to run flirt:

 $FREESURFER_HOME/bin/flirt.fsl

 It seems that I somehow do not have that file to complete tracula.  I was 
 hoping someone could help with this or perhaps send me the file?


 I am using freesurfer-
 Darwin-snowleopard-i686-stable-pub-v5.3.0-
 MacBook-Pro.local 10.8.0 Darwin Kernel Version 10.8.0:

 Hope to hear from someone soon,
 Andrew
 ___
 Freesurfer mailing list
 Freesurfer@nmr.mgh.harvard.edu
 https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer


___
Freesurfer mailing list
Freesurfer@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.



Re: [Freesurfer] kvlQuantifyHippocampalSubfieldSegmentations.sh empty text files

2013-07-17 Thread Jasmeet Hayes
That worked, thanks for your help!


On Wed, Jul 17, 2013 at 12:00 PM, Juan Eugenio Iglesias 
igles...@nmr.mgh.harvard.edu wrote:

  OK, then the question is why kvlQuantifyPosteriorProbabilityImages is
 producing an empty file
 I think the problem is that kvlQuantifyHippocampalSubfieldSegmentations.sh
 behaves in a different way depending on the number arguments, etc. I would
 suggest that you go back to the original version of the script, and use it
 the way it's described in
 http://ftp.nmr.mgh.harvard.edu/fswiki/HippocampalSubfieldSegmentation,
 that is, setting the environment variable SUBJECTS_DIR and calling the
 script without arguments.
 Kind regards,
 /Eugenio




 On 07/17/2013 11:46 AM, Jasmeet Hayes wrote:

 Hi Eugenio,

  We are using version 5.1. I attached our script. We added to line 60,
 setting the subject directory. We also changed the statFileName (see line
 291 on the attached script).  I should also note that the volumeStats_left
 and right text files within each subject's mri folder is empty. Thanks for
 your help.


 On Wed, Jul 17, 2013 at 11:16 AM, Juan Eugenio Iglesias 
 igles...@nmr.mgh.harvard.edu wrote:

  Hi Jasmeet,
 what modification did you exactly make? Whih version of FreeSurfer is
 this?
 Cheers,
 /Eugenio


 On 07/17/2013 10:22 AM, Jasmeet Hayes wrote:

  Dear Freesurfer experts,

  I tried running the kvlQuantifyHippocampalSubfieldSegmentations.sh
 script and get two text files (left  right stats) but the only information
 in the text files is the subject ID number.

  However, on my terminal window, I can see the volume output for each
 subfield and for each subject (see below for an example of 1 subject). And
 there were no errors after running the script.

  Every subject in the folder went through the hippocampal subfield recon
 procedure, and I have verified that there were no errors (each subject's
 folder has the left and right posterior_* files).

  We made a slight modification to the kvl script, only to change the
 Subject directory. Any ideas about what might be the problem? Thanks!
 -Jasmeet

  cd .. Quantifying subject er002_recon/ right Doing right side cd mri
 THIS IS THE COMPRESSION LOOKUP TABLE FILE
 /Applications/freesurfer/data/GEMS/compressionLookupTable.txt PWD before
 running KVLQuantifyPosteriorProb IS:
 /Volumes/VA_Imaging/Projects/jpannuhayes/emoreg/HIPPOCAMPAL_SUBFIELDS/er002_recon/mri
 kvlQuantifyPosteriorProbabilityImages
 /Applications/freesurfer/data/GEMS/compressionLookupTable.txt
 posterior_Right-Hippocampus.mgz posterior_right_presubiculum.mgz
 posterior_right_CA1.mgz posterior_right_CA2-3.mgz
 posterior_right_fimbria.mgz posterior_right_subiculum.mgz
 posterior_right_CA4-DG.mgz posterior_right_hippocampal_fissure.mgz 
 volumeStats_right.txt volumeInVoxels: Right-Hippocampus: 3094.49
 right_presubiculum: 4373.15 right_CA1: 2483.23 right_CA2-3: 8192.11
 right_fimbria: 740.934 right_subiculum: 5402.79 right_CA4-DG: 4688.93
 right_hippocampal_fissure: 541.872 cd .. cd ..


  ___
 Freesurfer mailing 
 listfreesur...@nmr.mgh.harvard.eduhttps://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer


 --
 -
 Juan Eugenio Iglesias, 
 PhDhttp://www.jeiglesias.comigles...@nmr.mgh.harvard.edu
 Athinoula A. Martinos Center for Biomedical Imaging
 149 Thirteenth Street, Suite 2301
 Charlestown, Massachusetts 2129
 U.S.A.

 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.



 --
 -
 Juan Eugenio Iglesias, 
 PhDhttp://www.jeiglesias.comigles...@nmr.mgh.harvard.edu
 Athinoula A. Martinos Center for Biomedical Imaging
 149 Thirteenth Street, Suite 2301
 Charlestown, Massachusetts 2129
 U.S.A.


___
Freesurfer mailing list
Freesurfer@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.


Re: [Freesurfer] SWI acquisition question

2013-07-17 Thread Fotiadis, Panagiotis
Hi Bruce,

Don't worry about it!
In addition, I was wondering about something else as well. I've been using 
MRIcron to analyze some of my Susceptibility Clinical Scans and even though 
they used to open fine, some of them (like the one I have attached in this 
email) appear distorted in a way. In the images I attached the Original one 
is how the Susceptibility clinical image should appear like, and the 
Distorted is how it appears for some of my subjects. Would you happen to know 
why is this happening?
Thanks again for your time and help!

Best,
Panos

From: Bruce Fischl [fis...@nmr.mgh.harvard.edu]
Sent: Wednesday, July 17, 2013 11:54 AM
To: Fotiadis, Panagiotis
Cc: freesurfer@nmr.mgh.harvard.edu
Subject: Re: [Freesurfer] SWI acquisition question

Hi Panos

sorry, I think this is pretty far outside our expertise, and it would be
hard in any case to try to figure out what happened to give you these
strange looking results. I certainly don't have any insight

sorry,
Bruce
On Wed, 17 Jul
2013, Fotiadis, Panagiotis wrote:

 Hi FS Community,

 I was wondering whether anyone had a chance to look into this. Thanks again 
 for your time!

 Best,
 Panos
 
 From: freesurfer-boun...@nmr.mgh.harvard.edu 
 [freesurfer-boun...@nmr.mgh.harvard.edu] on behalf of Fotiadis, Panagiotis
 Sent: Tuesday, July 16, 2013 2:00 PM
 To: freesurfer@nmr.mgh.harvard.edu
 Subject: [Freesurfer] SWI acquisition question

 Hi FreeSurfer community,

 I am analyzing some SWI scans of some old subjects. For most of my subjects 
 when a SWI was acquired, the resulting files that I would get would be a 
 Magnitude image, a Phase image, a mIP image, and the SWI image. However, I 
 have a couple of subjects for whom I get only the mIP and the SWI images, and 
 the SWI image looks more blurry than it's supposed to be. I was wondering:

 1) Why do I get only the mIP and the SWI images for those subjects and not 
 the magnitude and phase images as well,
 2) Is it normal for the mIP to look like that? Usually it looks different, and
 3) Why is the final SWI image more blurry than usual? By looking at the 
 image, I believe that maybe the head coil was not properly connected to the 
 scanner at the time of the acquisition but I don't know whether that would be 
 a sufficient explanation.

 I have attached two attachments: One of the resulting mIP and one of the 
 equivalent blurry SWI of the same subject.

 Thank you in advance for your help and time!

 Best,
 Panos

attachment: Distorted.jpegattachment: Original.jpeg___
Freesurfer mailing list
Freesurfer@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] Gray matter missed in the pial layer after mri_segment

2013-07-17 Thread Erin Browning
Hello--

We have a set of scans with bad artifact throughout the temporal lobe and
the superior cortex. We've had to do the normalization steps (n3
correction, normalization) manually on most of these scans to preserve the
gm/wm boundary. Now, when we run mri_segment, it leaves out some gray
matter along the pial border. Is there a correction for this? Mri_segment
has wlo and whi and ghi, but no glo, and it doesn't look like there's a
method to add gray matter to the pial boundary like there is to remove it.

Thank you,
Erin Browning
___
Freesurfer mailing list
Freesurfer@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] Crash in mri_register_ca during recon-all -base id -tp id1 -tp id2 -all

2013-07-17 Thread Knut J Bjuland
Hi,
I have some trouble running recon-all longitudinal pipelines. When 
recon-all comes to mri_register_ca it crashes without giving a message, 
even tough 8g ram was available. I have run the program with -debug. Is 
there any other option I may use? The system is running Rocks 6.0 
(Mamba), and it uses Slurm queue system.

Can I use the FreeSurfer 5.3 longitudinal process with data processed 
with FreeSurfer 5.1? I wish to use the new Longitudinal linear mixed 
effects model  MATLAB tools.

Knut J
___
Freesurfer mailing list
Freesurfer@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.



Re: [Freesurfer] beta weights from FS-Fast analysis

2013-07-17 Thread Joseph Dien
I implemented the ROI percent signal change formula following the MarsBaR FAQ 
(http://marsbar.sourceforge.net/faq.html) but the values I'm getting seem too 
small (on the order of .0002%).  Basically the formula is the (beta * peak 
absolute value of the canonical HRF regressor * 100)/(run mean).  No 
derivatives in this case as it is a boxcar design.

I took the mean across all the runs since FSFAST uses the same regressor across 
the entire experiment (unlike SPM).
I used the X.runflac(1).flac.ev(m).Xirf values for the canonical HRF as you 
suggested (where m equals the condition+1). 

Is it possible that I'm missing something in the scaling here?  Especially with 
a boxcar design, the signal change should be much larger than this for a 
significant cluster, I think.  For example, the peak HRF value for one of the 
conditions is 0.0092.  If the betas are already scaled according to the peak 
value, then it would come out as .02%, which is more reasonable, although still 
too small.

Thanks for your help with this!

Joe



On May 31, 2013, at 5:02 PM, Douglas N Greve gr...@nmr.mgh.harvard.edu wrote:

 
 Oh, right, it is probably not there for subcortical. I don't know what I 
 would have to do to write it out. It won't be something that happens 
 before I get back from HBM. Can you remind me after HBM?
 doug
 
 On 05/31/2013 04:44 PM, Joseph Dien wrote:
 It looks like the corrected vertex p-values 
 (ex: cache.th13.abs.sig.voxel.nii.gz) are only available for the 
 surface-based lh and rh spaces.  For the subcortical volume-based 
 analysis I don't see the corresponding corrected voxel p-values being 
 available?
 
 On May 31, 2013, at 2:46 PM, Joseph Dien jdie...@mac.com 
 mailto:jdie...@mac.com wrote:
 
 
 On May 31, 2013, at 12:11 PM, Douglas N Greve 
 gr...@nmr.mgh.harvard.edu mailto:gr...@nmr.mgh.harvard.edu wrote:
 
 
 On 05/31/2013 01:49 AM, Joseph Dien wrote:
 I was able to make more progress so I'm mostly good at this point but
 I have a remaining question:
 
 I assume the contents of sig.nii.gz (which I assume are the vertex
 p-values) are not FWE corrected.  Is it possible to get FWE-corrected
 vertex p-values?  Or are only clusterwise corrections available?
 There should be something like cache.th13.abs.sig.voxel.mgh which is
 corrected on a voxelwise basis (the th13 is just part of the name 
 but it
 should be the same regardless of the threshold you choose)
 doug
 
 Excellent!  Thanks!  :)
 
 
 Thanks again for your patience!
 
 Joe
 
 On May 30, 2013, at 4:37 PM, Joseph Dien jdie...@mac.com 
 mailto:jdie...@mac.com
 mailto:jdie...@mac.com wrote:
 
 Just to make sure I'm doing this right, I'm going to summarize what
 I've taken away from your answers and to ask some new questions. In
 order to present the results, I need two things:
 
 1) A set of histograms (with error bars) for each cluster figure to
 show the % signal change for each of the four contrasts of interest.
 The cache.th20.pos.y.ocn.dat file only gives it for the condition
 where the cluster was significant so I can't use that.
 So I could use mri_label2vol to convert cache.th20.neg.sig.ocn.annot
 from the group level analysis to generate a mask for each cluster of
 interest.
 Then I could extract the value of the voxels from each
 subject's cespct file for each contrast, average them across the
 cluster ROI, then average them across each subject, to generate the
 histogram?
 This would suffice to give me the %age signal change?
 I would be doing these computations in Matlab using MRIread.
 
 2) A results table with the headings:
 
 Cluster p (FWE corrected)
 Cluster size
 Peak Voxel p (FWE corrected)
 Peak Voxel T
 Peak Voxel Coords
 BA
 Anatomical Landmark
 
 I can get the first two from
 the cache.th20.pos/neg.sig.cluster.summary files from the group level
 analysis.
 I can get the peak voxel coordinates from the summary files as well.
 I can use this to get the peak voxel p from the group
 level sig.nii.gz file.  Is this FWE corrected?  If not, how can I get
 this information?
 I can use these coordinates to get the peak voxel T by getting the
 value from the group level F.nii.gz file and taking its square root.
 How can I get the sign of the T statistic?
 I can use the Lancaster transform to convert the MNI305 peak voxel
 coordinates into the Atlas coordinates to look up the putative BA and
 landmarks (unless there is a better way with Freesurfer?  I'm seeing
 some references to some BA labels in the forum but it doesn't look
 like this is a complete set yet?).
 
 Sorry for all these questions!  I got some nice results from FSFAST
 and would like to get them written up.
 
 Cheers!
 
 Joe
 
 
 
 
 On May 29, 2013, at 10:53 PM, Douglas Greve
 gr...@nmr.mgh.harvard.edu mailto:gr...@nmr.mgh.harvard.edu 
 mailto:gr...@nmr.mgh.harvard.edu wrote:
 
 
 On 5/29/13 10:42 PM, Joseph Dien wrote:
 
 On May 29, 2013, at 11:40 AM, Douglas N Greve
 gr...@nmr.mgh.harvard.edu mailto:gr...@nmr.mgh.harvard.edu 
 

Re: [Freesurfer] beta weights from FS-Fast analysis

2013-07-17 Thread Joseph Dien
then I get on the order of .02% difference between the contrasted conditions.
The run mean values are in my expected ballpark of about 100 or so.
The condition betas are just very very small.
Or perhaps this is typical of FSFAST analyses?

On Jul 17, 2013, at 2:00 PM, Douglas N Greve gr...@nmr.mgh.harvard.edu wrote:

 
 The beta's have already been scaled. What do you get if you just 
 beta/runmean ?
 
 
 
 On 07/17/2013 01:45 PM, Joseph Dien wrote:
 I implemented the ROI percent signal change formula following the 
 MarsBaR FAQ (http://marsbar.sourceforge.net/faq.html) but the values 
 I'm getting seem too small (on the order of .0002%).  Basically the 
 formula is the (beta * peak absolute value of the canonical HRF 
 regressor * 100)/(run mean).  No derivatives in this case as it is a 
 boxcar design.
 
 I took the mean across all the runs since FSFAST uses the same 
 regressor across the entire experiment (unlike SPM).
 I used the X.runflac(1).flac.ev(m).Xirf values for the canonical HRF 
 as you suggested (where m equals the condition+1).
 
 Is it possible that I'm missing something in the scaling here? 
 Especially with a boxcar design, the signal change should be much 
 larger than this for a significant cluster, I think.  For example, the 
 peak HRF value for one of the conditions is 0.0092.  If the betas are 
 already scaled according to the peak value, then it would come out as 
 .02%, which is more reasonable, although still too small.
 
 Thanks for your help with this!
 
 Joe
 
 
 
 On May 31, 2013, at 5:02 PM, Douglas N Greve 
 gr...@nmr.mgh.harvard.edu mailto:gr...@nmr.mgh.harvard.edu wrote:
 
 
 Oh, right, it is probably not there for subcortical. I don't know what I
 would have to do to write it out. It won't be something that happens
 before I get back from HBM. Can you remind me after HBM?
 doug
 
 On 05/31/2013 04:44 PM, Joseph Dien wrote:
 It looks like the corrected vertex p-values
 (ex: cache.th13.abs.sig.voxel.nii.gz) are only available for the
 surface-based lh and rh spaces.  For the subcortical volume-based
 analysis I don't see the corresponding corrected voxel p-values being
 available?
 
 On May 31, 2013, at 2:46 PM, Joseph Dien jdie...@mac.com 
 mailto:jdie...@mac.com
 mailto:jdie...@mac.com wrote:
 
 
 On May 31, 2013, at 12:11 PM, Douglas N Greve
 gr...@nmr.mgh.harvard.edu mailto:gr...@nmr.mgh.harvard.edu 
 mailto:gr...@nmr.mgh.harvard.edu wrote:
 
 
 On 05/31/2013 01:49 AM, Joseph Dien wrote:
 I was able to make more progress so I'm mostly good at this point but
 I have a remaining question:
 
 I assume the contents of sig.nii.gz (which I assume are the vertex
 p-values) are not FWE corrected.  Is it possible to get FWE-corrected
 vertex p-values?  Or are only clusterwise corrections available?
 There should be something like cache.th13.abs.sig.voxel.mgh which is
 corrected on a voxelwise basis (the th13 is just part of the name
 but it
 should be the same regardless of the threshold you choose)
 doug
 
 Excellent!  Thanks!  :)
 
 
 Thanks again for your patience!
 
 Joe
 
 On May 30, 2013, at 4:37 PM, Joseph Dien jdie...@mac.com 
 mailto:jdie...@mac.com
 mailto:jdie...@mac.com
 mailto:jdie...@mac.com wrote:
 
 Just to make sure I'm doing this right, I'm going to summarize what
 I've taken away from your answers and to ask some new questions. In
 order to present the results, I need two things:
 
 1) A set of histograms (with error bars) for each cluster figure to
 show the % signal change for each of the four contrasts of interest.
 The cache.th20.pos.y.ocn.dat file only gives it for the condition
 where the cluster was significant so I can't use that.
 So I could use mri_label2vol to convert cache.th20.neg.sig.ocn.annot
 from the group level analysis to generate a mask for each cluster of
 interest.
 Then I could extract the value of the voxels from each
 subject's cespct file for each contrast, average them across the
 cluster ROI, then average them across each subject, to generate the
 histogram?
 This would suffice to give me the %age signal change?
 I would be doing these computations in Matlab using MRIread.
 
 2) A results table with the headings:
 
 Cluster p (FWE corrected)
 Cluster size
 Peak Voxel p (FWE corrected)
 Peak Voxel T
 Peak Voxel Coords
 BA
 Anatomical Landmark
 
 I can get the first two from
 the cache.th20.pos/neg.sig.cluster.summary files from the group 
 level
 analysis.
 I can get the peak voxel coordinates from the summary files as well.
 I can use this to get the peak voxel p from the group
 level sig.nii.gz file.  Is this FWE corrected?  If not, how can 
 I get
 this information?
 I can use these coordinates to get the peak voxel T by getting the
 value from the group level F.nii.gz file and taking its square root.
 How can I get the sign of the T statistic?
 I can use the Lancaster transform to convert the MNI305 peak voxel
 coordinates into the Atlas coordinates to look up the putative 
 BA and
 landmarks (unless there is a better way with 

Re: [Freesurfer] Q) Medial_wall.label

2013-07-17 Thread Bruce Fischl

Hi Yune

you are better off using the ?h.cortex.label file to mask in what is 
needed for each subject


cheers
Bruce

On Wed, 17 Jul 2013, Yune Lee wrote:


Hello freesurfer experts, 
For visualization purpose (i.e., displaying group maps), I applied a mask
for medial wall regions using the files (e.g., rh.Medial_wall.label or
lh.Medial.wall.label), which I found in fsaverage/label/ directory. 

However, I noticed a somewhat imperfect medial wall masking in both
hemispheres, such that 
transversing lines are seen in unthresholded t maps for both
hemispheres (see attached ). 

Those lines eventually disappear when I increased threshold, but I wonder
why this happened. 

Any ideas? 

best,
yune 




___
Freesurfer mailing list
Freesurfer@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.


Re: [Freesurfer] Gray matter missed in the pial layer after mri_segment

2013-07-17 Thread ye tian
Hello Erin,

Thank you so much for asking. I have been curious to know this, too.

Sincerely,
Ye


On Wed, Jul 17, 2013 at 12:41 PM, Erin Browning brown...@uwm.edu wrote:

 Hello--

 We have a set of scans with bad artifact throughout the temporal lobe and
 the superior cortex. We've had to do the normalization steps (n3
 correction, normalization) manually on most of these scans to preserve the
 gm/wm boundary. Now, when we run mri_segment, it leaves out some gray
 matter along the pial border. Is there a correction for this? Mri_segment
 has wlo and whi and ghi, but no glo, and it doesn't look like there's a
 method to add gray matter to the pial boundary like there is to remove it.

 Thank you,
 Erin Browning

 ___
 Freesurfer mailing list
 Freesurfer@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@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.


Re: [Freesurfer] beta weights from FS-Fast analysis

2013-07-17 Thread Joseph Dien
It's a boxcar design so 20.265.

  mkanalysis-sess -fsd bold -analysis CPA.sm05.lh -surface fsaverage lh 
-fwhm 5 -event-related  -paradigm CPA1.par -nconditions 20 -spmhrf 0 -TR 2 
-refeventdur 20.265 -polyfit 2 -per-run -force -nuisreg nuisreg2.dat 7 
-tpexclude tpexclude.dat

On Jul 17, 2013, at 3:50 PM, Douglas N Greve gr...@nmr.mgh.harvard.edu wrote:

 when you ran mkanalysis-sess, what did you set --refeventdur to?
 On 07/17/2013 02:50 PM, Joseph Dien wrote:
 then I get on the order of .02% difference between the contrasted conditions.
 The run mean values are in my expected ballpark of about 100 or so.
 The condition betas are just very very small.
 Or perhaps this is typical of FSFAST analyses?
 
 On Jul 17, 2013, at 2:00 PM, Douglas N Greve gr...@nmr.mgh.harvard.edu 
 mailto:gr...@nmr.mgh.harvard.edu wrote:
 
 
 The beta's have already been scaled. What do you get if you just
 beta/runmean ?
 
 
 
 On 07/17/2013 01:45 PM, Joseph Dien wrote:
 I implemented the ROI percent signal change formula following the
 MarsBaR FAQ (http://marsbar.sourceforge.net/faq.html) but the values
 I'm getting seem too small (on the order of .0002%).  Basically the
 formula is the (beta * peak absolute value of the canonical HRF
 regressor * 100)/(run mean).  No derivatives in this case as it is a
 boxcar design.
 
 I took the mean across all the runs since FSFAST uses the same
 regressor across the entire experiment (unlike SPM).
 I used the X.runflac(1).flac.ev(m).Xirf values for the canonical HRF
 as you suggested (where m equals the condition+1).
 
 Is it possible that I'm missing something in the scaling here?
 Especially with a boxcar design, the signal change should be much
 larger than this for a significant cluster, I think.  For example, the
 peak HRF value for one of the conditions is 0.0092.  If the betas are
 already scaled according to the peak value, then it would come out as
 .02%, which is more reasonable, although still too small.
 
 Thanks for your help with this!
 
 Joe
 
 
 
 On May 31, 2013, at 5:02 PM, Douglas N Greve
 gr...@nmr.mgh.harvard.edu mailto:gr...@nmr.mgh.harvard.edu 
 mailto:gr...@nmr.mgh.harvard.edu wrote:
 
 
 Oh, right, it is probably not there for subcortical. I don't know what I
 would have to do to write it out. It won't be something that happens
 before I get back from HBM. Can you remind me after HBM?
 doug
 
 On 05/31/2013 04:44 PM, Joseph Dien wrote:
 It looks like the corrected vertex p-values
 (ex: cache.th13.abs.sig.voxel.nii.gz) are only available for the
 surface-based lh and rh spaces.  For the subcortical volume-based
 analysis I don't see the corresponding corrected voxel p-values being
 available?
 
 On May 31, 2013, at 2:46 PM, Joseph Dien jdie...@mac.com 
 mailto:jdie...@mac.com
 mailto:jdie...@mac.com
 mailto:jdie...@mac.com wrote:
 
 
 On May 31, 2013, at 12:11 PM, Douglas N Greve
 gr...@nmr.mgh.harvard.edu mailto:gr...@nmr.mgh.harvard.edu 
 mailto:gr...@nmr.mgh.harvard.edu
 mailto:gr...@nmr.mgh.harvard.edu wrote:
 
 
 On 05/31/2013 01:49 AM, Joseph Dien wrote:
 I was able to make more progress so I'm mostly good at this point but
 I have a remaining question:
 
 I assume the contents of sig.nii.gz (which I assume are the vertex
 p-values) are not FWE corrected.  Is it possible to get FWE-corrected
 vertex p-values?  Or are only clusterwise corrections available?
 There should be something like cache.th13.abs.sig.voxel.mgh which is
 corrected on a voxelwise basis (the th13 is just part of the name
 but it
 should be the same regardless of the threshold you choose)
 doug
 
 Excellent!  Thanks!  :)
 
 
 Thanks again for your patience!
 
 Joe
 
 On May 30, 2013, at 4:37 PM, Joseph Dien jdie...@mac.com 
 mailto:jdie...@mac.com
 mailto:jdie...@mac.com
 mailto:jdie...@mac.com
 mailto:jdie...@mac.com wrote:
 
 Just to make sure I'm doing this right, I'm going to summarize what
 I've taken away from your answers and to ask some new questions. In
 order to present the results, I need two things:
 
 1) A set of histograms (with error bars) for each cluster figure to
 show the % signal change for each of the four contrasts of interest.
 The cache.th20.pos.y.ocn.dat file only gives it for the condition
 where the cluster was significant so I can't use that.
 So I could use mri_label2vol to convert cache.th20.neg.sig.ocn.annot
 from the group level analysis to generate a mask for each cluster of
 interest.
 Then I could extract the value of the voxels from each
 subject's cespct file for each contrast, average them across the
 cluster ROI, then average them across each subject, to generate the
 histogram?
 This would suffice to give me the %age signal change?
 I would be doing these computations in Matlab using MRIread.
 
 2) A results table with the headings:
 
 Cluster p (FWE corrected)
 Cluster size
 Peak Voxel p (FWE corrected)
 Peak Voxel T
 Peak Voxel Coords
 BA
 Anatomical Landmark
 
 I can get the first two from
 the cache.th20.pos/neg.sig.cluster.summary files 

Re: [Freesurfer] beta weights from FS-Fast analysis

2013-07-17 Thread Joseph Dien
It's a little complicated.  Basically there were eight runs, comprising four 
conditions (me, we, you1, you2) each with two adjoining runs.  For the 
analysis, I merged each of the pairs into a single run and added a nuisance 
regressor to account for the difference in run means.  There were a total of 
four different kinds of boxcars (AR, CS, EM, MP).  So 4x4=16 conditions.  There 
was also a covariate of non-interest to mark the switch point for each boxcar, 
one for each run, so 20 total.

The 7 nuisance regressors are six movement covariates plus one to account for 
merging eight runs into four (it consists of 1 for the first half and -1 for 
the second, so the difference in the run means).  I'm using the movement 
covariates from a prior SPM run since ARTdetect (for detecting bad volumes) 
isn't set up for AFNI style data.  From all published accounts the different 
movement detection routines yield similar enough results that it shouldn't be a 
problem (consistent with what I found when I compared them for this dataset).

You're thinking that collinearity could have reduced the effect sizes?  When I 
correlate the X.X regressor matrix, the 20 predictors don't correlate by more 
than about .2 at worst.  I do see greater correlations with some of the 
nuisance regressors (as high as the .4 range).  Are my betas unusually small 
for FSFAST analyses?  They did come up clusterwise significant at least. Or 
should I not worry?  I'm not sure what to expect from FSFAST analyses.

Thanks!

Joe


On Jul 17, 2013, at 3:58 PM, Douglas N Greve gr...@nmr.mgh.harvard.edu wrote:

 why do you have 20 conditions? And what are the 7 nuisance regressors?
 
 On 07/17/2013 03:54 PM, Joseph Dien wrote:
 It's a boxcar design so 20.265.
 
 mkanalysis-sess -fsd bold -analysis CPA.sm05.lh -surface fsaverage lh 
 -fwhm 5 -event-related-paradigm CPA1.par -nconditions 20 -spmhrf 0 -TR 
 2 -refeventdur 20.265 -polyfit 2 -per-run -force -nuisreg nuisreg2.dat 
 7 -tpexclude tpexclude.dat
 
 On Jul 17, 2013, at 3:50 PM, Douglas N Greve 
 gr...@nmr.mgh.harvard.edu mailto:gr...@nmr.mgh.harvard.edu wrote:
 
 when you ran mkanalysis-sess, what did you set --refeventdur to?
 On 07/17/2013 02:50 PM, Joseph Dien wrote:
 then I get on the order of .02% difference between the contrasted 
 conditions.
 The run mean values are in my expected ballpark of about 100 or so.
 The condition betas are just very very small.
 Or perhaps this is typical of FSFAST analyses?
 
 On Jul 17, 2013, at 2:00 PM, Douglas N Greve 
 gr...@nmr.mgh.harvard.edu mailto:gr...@nmr.mgh.harvard.edu 
 mailto:gr...@nmr.mgh.harvard.edu wrote:
 
 
 The beta's have already been scaled. What do you get if you just
 beta/runmean ?
 
 
 
 On 07/17/2013 01:45 PM, Joseph Dien wrote:
 I implemented the ROI percent signal change formula following the
 MarsBaR FAQ (http://marsbar.sourceforge.net/faq.html) but the values
 I'm getting seem too small (on the order of .0002%).  Basically the
 formula is the (beta * peak absolute value of the canonical HRF
 regressor * 100)/(run mean).  No derivatives in this case as it is a
 boxcar design.
 
 I took the mean across all the runs since FSFAST uses the same
 regressor across the entire experiment (unlike SPM).
 I used the X.runflac(1).flac.ev(m).Xirf values for the canonical HRF
 as you suggested (where m equals the condition+1).
 
 Is it possible that I'm missing something in the scaling here?
 Especially with a boxcar design, the signal change should be much
 larger than this for a significant cluster, I think.  For example, the
 peak HRF value for one of the conditions is 0.0092.  If the betas are
 already scaled according to the peak value, then it would come out as
 .02%, which is more reasonable, although still too small.
 
 Thanks for your help with this!
 
 Joe
 
 
 
 On May 31, 2013, at 5:02 PM, Douglas N Greve
 gr...@nmr.mgh.harvard.edu mailto:gr...@nmr.mgh.harvard.edu 
 mailto:gr...@nmr.mgh.harvard.edu 
 mailto:gr...@nmr.mgh.harvard.edu wrote:
 
 
 Oh, right, it is probably not there for subcortical. I don't know 
 what I
 would have to do to write it out. It won't be something that happens
 before I get back from HBM. Can you remind me after HBM?
 doug
 
 On 05/31/2013 04:44 PM, Joseph Dien wrote:
 It looks like the corrected vertex p-values
 (ex: cache.th13.abs.sig.voxel.nii.gz) are only available for the
 surface-based lh and rh spaces.  For the subcortical volume-based
 analysis I don't see the corresponding corrected voxel p-values 
 being
 available?
 
 On May 31, 2013, at 2:46 PM, Joseph Dien jdie...@mac.com 
 mailto:jdie...@mac.com mailto:jdie...@mac.com
 mailto:jdie...@mac.com
 mailto:jdie...@mac.com wrote:
 
 
 On May 31, 2013, at 12:11 PM, Douglas N Greve
 gr...@nmr.mgh.harvard.edu mailto:gr...@nmr.mgh.harvard.edu 
 mailto:gr...@nmr.mgh.harvard.edu 
 mailto:gr...@nmr.mgh.harvard.edu
 mailto:gr...@nmr.mgh.harvard.edu wrote:
 
 
 On 05/31/2013 01:49 AM, Joseph Dien wrote:
 I was able to make more progress so I'm mostly 

[Freesurfer] Cortical surface

2013-07-17 Thread pablo najt
Dear Freesurfer experts,  I am 
interested in running cortical surface area on longitudinal data from MRI scans 
and run multisubject comparisons (patients vs controls). I would like to ask 
you where could I find either a tutorial or a guide for performing such 
analysis. Also I am not familiar with freesurfer, so I am not sure which type 
of preprocessing should be undertaken.I would greatly appreciate your valuable 
help on this topic.Pablo 
  ___
Freesurfer mailing list
Freesurfer@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] recon-all steps specified into volume and surface based stream?

2013-07-17 Thread Martina Papmeyer
Dear FreeSurfer experts,

I was wondering if there is anywhere a table or hint, suggesting which  
particular steps of all the recon-all steps (i.e.  
https://surfer.nmr.mgh.harvard.edu/fswiki/ReconAllDevTable) are part  
of the volume-based or the surface-based stream or both? I also found  
this diagram  
https://surfer.nmr.mgh.harvard.edu/fswiki/ReconAllBlockDiagram but I'm  
unsure if the colours are supposed to give me any clues.

I'm familiar with the broad processing stages of the two streams off  
course but would very much love to understand the details better :)

Thank you very much for your help, Martina

-- 
The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.


___
Freesurfer mailing list
Freesurfer@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.



Re: [Freesurfer] beta weights from FS-Fast analysis

2013-07-17 Thread Douglas N Greve
Oh, right, the nuisance regressors are given a separate regressor for 
each run. I had forgotten!

On 07/17/2013 04:55 PM, Joseph Dien wrote:

 On Jul 17, 2013, at 4:50 PM, Douglas N Greve 
 gr...@nmr.mgh.harvard.edu mailto:gr...@nmr.mgh.harvard.edu wrote:


 On 07/17/2013 04:48 PM, Joseph Dien wrote:
 The single nuisance regressor models the difference between the merged
 runs.  So if the first run had a mean of 90 and the second run had a
 mean of 110 then the merged run mean would be 100.  The nuisance
 regressor has 1s for the volumes of the first run and -1s for the
 volumes of the second run so it ends up with a beta value of 10, thus
 accounting for the difference between these two sets of volumes.  I
 think it does make sense.
 Do you do this in a single regressor? So you would have a +1 -1 pattern
 repeated 4 times? I think to make it work, you would need 4 regressors.
 In any event, FSFAST will do the right thing, so maybe it is not 
 important.


 Looking at the X.X file, it did create four nuisance regressors.  I'm 
 really impressed with how well FSFAST handles all this!  :)


 In any case, while it was necessary to do so for my original SPM
 analyses since it uses a separate covariates for each run, after
 working through the FSFAST procedures with your help I see that is not
 the case for FSFAST (a single regressor models a given condition in
 all the runs).  I'll try doing as you suggest to see what difference
 it makes.

 It is indeed cognitive areas and the manipulations are subtle social
 cognition manipulations so perhaps not surprising after all.

 I'll send you the paradigm file separately.

 Thanks for taking the time to look into this!

 Joe


 On Jul 17, 2013, at 4:38 PM, Douglas N Greve
 gr...@nmr.mgh.harvard.edu mailto:gr...@nmr.mgh.harvard.edu 
 mailto:gr...@nmr.mgh.harvard.edu wrote:


 It is not necessary or beneficial to combine the runs in this way.
 FSFAST will do all this for you and keep account of all the runs and
 transitions. FSFAST will put in regressors to fit each of the run 
 means.
 The single regressor you have is not the right way to hand this (at
 least I don't understand how it works). It could be that the low %
 signal change is related to the colinearity between the task waveform
 and the mean regressors. Can you set things up in the way that FSFAST
 expects them and don't use any nuisance regressors?

 Also, in what area are you looking at the percent change? .02% sounds
 very small, but may be if it is in some cognitive area, maybe it is ok.
 If it is in visual cortex, then it looks way too low.

 Also, can you send the paradigm file?

 doug




 On 07/17/2013 04:27 PM, Joseph Dien wrote:
 It's a little complicated.  Basically there were eight runs,
 comprising four conditions (me, we, you1, you2) each with two
 adjoining runs.  For the analysis, I merged each of the pairs into a
 single run and added a nuisance regressor to account for the
 difference in run means.  There were a total of four different kinds
 of boxcars (AR, CS, EM, MP).  So 4x4=16 conditions.  There was also a
 covariate of non-interest to mark the switch point for each boxcar,
 one for each run, so 20 total.

 The 7 nuisance regressors are six movement covariates plus one to
 account for merging eight runs into four (it consists of 1 for the
 first half and -1 for the second, so the difference in the run means).
 I'm using the movement covariates from a prior SPM run since
 ARTdetect (for detecting bad volumes) isn't set up for AFNI style
 data.  From all published accounts the different movement detection
 routines yield similar enough results that it shouldn't be a problem
 (consistent with what I found when I compared them for this dataset).

 You're thinking that collinearity could have reduced the effect sizes?
 When I correlate the X.X regressor matrix, the 20 predictors don't
 correlate by more than about .2 at worst.  I do see greater
 correlations with some of the nuisance regressors (as high as the .4
 range).  Are my betas unusually small for FSFAST analyses?  They did
 come up clusterwise significant at least. Or should I not worry?  I'm
 not sure what to expect from FSFAST analyses.

 Thanks!

 Joe


 On Jul 17, 2013, at 3:58 PM, Douglas N Greve
 gr...@nmr.mgh.harvard.edu mailto:gr...@nmr.mgh.harvard.edu 
 mailto:gr...@nmr.mgh.harvard.edu
 mailto:gr...@nmr.mgh.harvard.edu wrote:

 why do you have 20 conditions? And what are the 7 nuisance 
 regressors?

 On 07/17/2013 03:54 PM, Joseph Dien wrote:
 It's a boxcar design so 20.265.

 mkanalysis-sess -fsd bold -analysis CPA.sm05.lh -surface 
 fsaverage lh
 -fwhm 5 -event-related-paradigm CPA1.par -nconditions 20 -spmhrf 
 0 -TR
 2 -refeventdur 20.265 -polyfit 2 -per-run -force -nuisreg 
 nuisreg2.dat
 7 -tpexclude tpexclude.dat

 On Jul 17, 2013, at 3:50 PM, Douglas N Greve s.
 gr...@nmr.mgh.harvard.edu mailto:gr...@nmr.mgh.harvard.edu 
 mailto:gr...@nmr.mgh.harvard.edu
 mailto:gr...@nmr.mgh.harvard.edu
 

Re: [Freesurfer] Artefactual skull stripping for some subjects

2013-07-17 Thread Tudor Popescu
Dear Bruce,

Could you please let me know whether my message of last week, in which I
uploaded the subject's volssurfs, have reached the list? I haven't
received a reply about those questions and was wondering whether that email
maybe got rejected due to the large size of the attachments. The questions
should be visible in plaintext below, please let me know if I should upload
the files on a sharing server

Thanks!
Tudor

On 15 July 2013 19:51, Tudor Popescu tud...@gmail.com wrote:

 Hi again,

 I was wondering whether my previous message was received, since it had
 rather large attachments. I look forward to receiving a reply with regards
 to my questions!

 Thank you,
 Tudor


 On 12 July 2013 16:03, Tudor Popescu tud...@gmail.com wrote:

 Thanks Bruce. So the massive chunks that appear as missing in freeview,
 e.g. at coronal slice 153, show correctly in tkmedit, although the surfaces
 don't exist there. If I mouse above that area, brainmask.mgz intensities
 (which seem to be always identical to T1.mgz ones) seem quite normal to me,
 i.e. they're not all 0 but vary in a range centred around ~100.

 As for the segmentation (aseg) - not quite sure where I need to look for
 that. I attach here the main volumes  surfaces of this subject. Could my
 -wsthresh 35 recon not have replaced the existing one?

 Many thanks again!
 Tudor

 PS: why does freeview show the coordinates in the order [x z y] instead
 of [x y z]? I.e., as I move from anterior to posterior, I'd expect the
 middle coordinate to decrease, rather than the third one.

 PPS: what is the difference between brain.mgz and brainmask.mgz? Is it to
 do with conformity and intensity correction, just like the difference
 between 001.mgz and T1.mgz? I would have thought that brainmask.mgz would
 be a binary image (a mask consisting of 0s and 1s), rather than a regular
 (skull-stripped) structural.


 On 11 July 2013 21:40, Bruce Fischl fis...@nmr.mgh.harvard.edu wrote:

 Hi Tudor

 can you cc the list so others can answer? What are the intensities in
 the excluded right temporal wm in the brainmask.mgz? What is the aseg like
 there? Looks like an intensity normalization failure, but I'm not sure why.
 If you want to upload the subject I'll take a look
 Bruce

 On Thu, 11 Jul 2013, Tudor Popescu wrote:

  Hi Bruce,

 Sorry, wasn't sure whether you'd seen my previous message. So in the
 meantime, after the following command for subject C06 (whose artefactual
 skull strip I previously included as a screenshot) completed without
 error,

 recon-all -skullstrip -wsthresh 35 -clean-bm -subjid C06

 ..the volumessurfaces still look bad (see screenshot and recon log
 attached).

 C06 is just one of several subjects for which this has happened. Is
 there
 anything else I can try? Thanks!

 Tudor

 On 9 July 2013 13:33, Tudor Popescu tud...@gmail.com wrote:
   Thanks for clarifying this Doug!

   Please let me know if anyone has suggestions about my other
   questions below, re skull stripping.

   Cheers!
   T

   On 8 July 2013 18:01, Douglas N Greve
   gr...@nmr.mgh.harvard.edu wrote:

 001.mgz - raw.mgz - orig.mgz - nu.mgz - T1.mgz

 If you only have one run, then 001 and raw are the
 same
 orig.mgz is the first conformed volume

 doug



 On 07/08/2013 12:38 PM, Tudor Popescu wrote:
  Thanks Bruce. I have looked at that part of the
 wiki but I was
  wondering whether the screenshots at all suggest
 that all artefacts
  are likely to be the result of the same problem,
 and that e.g.
  re-running recon for all subjects with a
 watershed25 and with the
  -no-wsgcaatlas flag might be the thing to do in
 this case, rather than
  manually editing slices in tkmedit.
 
  So after making these corrections to all affected
 subjects, before
  proceeding to analyses, is a visual inspection in
 freeview of the
  brain.mgz of all subjects enough? It seems like
 some of these
  skull-stripping errors are quite subtle and can be
 easily missed (as I
  have, in fact, when I checked the output after my
 first recon-all).
 
  the T1.mgz is intensity corrected and
 conformed (the orig.mgz is
  just conformed).
 
  So are both derived from 001.mgz? Otherwise, what
 is the difference
  between 001.mgz and orig.mgz?
 
  Thanks!
  Tudor
 
 
  __**_
  Freesurfer mailing list
  Freesurfer@nmr.mgh.harvard.edu
  

Re: [Freesurfer] Artefactual skull stripping for some subjects

2013-07-17 Thread Bruce Fischl

Hi Tudor,

yes, it got through. It's on our list of things to get to, we just 
haven't gotten there yet

Bruce
On Wed, 17 Jul 2013, Tudor Popescu wrote:


Dear Bruce,

Could you please let me know whether my message of last week, in which I
uploaded the subject's volssurfs, have reached the list? I haven't received
a reply about those questions and was wondering whether that email maybe got
rejected due to the large size of the attachments. The questions should be
visible in plaintext below, please let me know if I should upload the files
on a sharing server

Thanks!
Tudor

On 15 July 2013 19:51, Tudor Popescu tud...@gmail.com wrote:
  Hi again,

  I was wondering whether my previous message was received, since
  it had rather large attachments. I look forward to receiving a
  reply with regards to my questions!

  Thank you,
  Tudor

  On 12 July 2013 16:03, Tudor Popescu tud...@gmail.com wrote:
Thanks Bruce. So the massive chunks that appear as
missing in freeview, e.g. at coronal slice 153, show
correctly in tkmedit, although the surfaces don't
exist there. If I mouse above that area,
brainmask.mgz intensities (which seem to be always
identical to T1.mgz ones) seem quite normal to me,
i.e. they're not all 0 but vary in a range centred
around ~100.

As for the segmentation (aseg) - not quite sure
where I need to look for that. I attach here the
main volumes  surfaces of this subject. Could my
-wsthresh 35 recon not have replaced the existing
one?

Many thanks again!
Tudor

PS: why does freeview show the coordinates in the
order [x z y] instead of [x y z]? I.e., as I move
from anterior to posterior, I'd expect the middle
coordinate to decrease, rather than the third one.

PPS: what is the difference between brain.mgz and
brainmask.mgz? Is it to do with conformity and
intensity correction, just like the difference
between 001.mgz and T1.mgz? I would have thought
that brainmask.mgz would be a binary image (a mask
consisting of 0s and 1s), rather than a regular
(skull-stripped) structural.

On 11 July 2013 21:40, Bruce Fischl
fis...@nmr.mgh.harvard.edu wrote:
  Hi Tudor

  can you cc the list so others can
  answer? What are the intensities in the
  excluded right temporal wm in the
  brainmask.mgz? What is the aseg like
  there? Looks like an intensity
  normalization failure, but I'm not sure
  why. If you want to upload the subject
  I'll take a look
  Bruce
  On Thu, 11 Jul 2013, Tudor Popescu
  wrote:

Hi Bruce,

Sorry, wasn't sure whether
you'd seen my previous
message. So in the
meantime, after the
following command for
subject C06 (whose
artefactual
skull strip I previously
included as a screenshot)
completed without error,

recon-all -skullstrip
-wsthresh 35 -clean-bm
-subjid C06

..the volumessurfaces still
look bad (see screenshot and
recon log
attached).

C06 is just one of several
subjects for which this has
happened. Is there
anything else I can try?
Thanks!

Tudor

On 9 July 2013 13:33, Tudor
Popescu tud...@gmail.com
wrote:
      Thanks for clarifying
this Doug!

      Please let me know if
anyone has suggestions about
my other
      questions below, re
skull stripping.

      Cheers!
      T

      On 8 July 2013 18:01,
Douglas N Greve
     
gr...@nmr.mgh.harvard.edu
wrote:

            001.mgz -
raw.mgz - orig.mgz -
nu.mgz - T1.mgz

            If you only have