Re: [caret-users] projecting functional MRI to gii surfaces

2014-08-05 Thread Tang, Yan
Maybe I do some wrong things.  I used  the Freesurfer_to_fs_LR Pipeline to get  
164k fs_LR surface. Then, I want to get the points in 164k fs_LR surface and 
use these points as seed to analysis the resting state functional connectivity. 
So I do project functional MRI to surface. But after that, I don't know how to 
do it in next step. I get metric file, but these files cannot be read in 
Matlab. And I also don't get the coordinates of the points in 164k fs_LR 
surface.

Just you mention the method of Alex Cohen. I am  beginner of Caret. Could you 
tell me the method in detail?

thank you

From: caret-users-boun...@brainvis.wustl.edu 
[caret-users-boun...@brainvis.wustl.edu] on behalf of Donna Dierker 
[do...@brainvis.wustl.edu]
Sent: Tuesday, August 05, 2014 10:26 AM
To: Caret, SureFit, and SuMS software users
Cc: Tang, Yiyuan
Subject: Re: [caret-users] projecting functional MRI to gii surfaces

You can do that, but I am not used to seeing this done with the mean 
midthickness surface.  Alex Cohen did something like this when he was using 
resting state functional connectivity to find gradients in a subject's cortical 
networks.  He used the PALS flat or spherical map to get a common grid onto a 
standard mesh.  Once there, he projected that grid onto the individual's 
midthickness surface on the standard mesh.  (Like the grid is getting folded 
back up into the individual's anatomical pattern.)  Then he unprojected the 
points to use as seeds for his analysis.

Do you mind if I ask how the mean midthickness surface comes into play?

Caret's Layers: Borders has options for making grids on the flat map.  People 
make grids on the sphere in matlab.

There are caret_command tools for unprojecting borders.

  caret_command -surface-border-unprojection

Border files differ from border projection files in that they are points not 
tied to a particular mesh.  The advantage is that you can open the same border 
point on, say, a native and 164k mesh, and it will align with both, if it 
aligns with one (and they are identical except for mesh).  The advantage of 
borderproj files is that they open on multiple configurations - flat, 
midthickness, inflated, etc. - but the price is that you're tied to a mesh.


On Aug 5, 2014, at 10:06 AM, "Tang, Yan"  wrote:

> Thank you for your help. Another problem is how to use the Caret software to 
> generate  regularly spaced Cartesian grids  on the flattened PALS-B12 average 
> surface of the left and right hemispheres.  Can I use this 3-dimensional (3D) 
> stereotactic coordinates from the PALS-B12 average fiducial (midthickness) 
> surface for each grid location to obtain the voxel coordinates (3 × 3 × 3 mm 
> resolution) containing that point in fMRI?
>
> 
> From: caret-users-boun...@brainvis.wustl.edu 
> [caret-users-boun...@brainvis.wustl.edu] on behalf of Donna Dierker 
> [donna.dier...@sbcglobal.net]
> Sent: Friday, August 01, 2014 5:34 PM
> To: Caret, SureFit, and SuMS software users
> Cc: Tang, Yiyuan
> Subject: Re: [caret-users] projecting functional MRI to gii surfaces
>
> Push Toolbar: D/C and make sure the primary overlay is Metric.
>
> Make sure the right column is selected.
>
> If that check out okay, then I would do:
>
> File: Open Data File: Volume Functional File
> Load the volume you just mapped
> Switch to volume view and select view All (as opposed to H (horizontal or 
> axial).
> Select D/C and on the page selection drop-down menu, scroll all the way to 
> the bottom
>something like volume surface outline
> Toggle on the fiducial surface used for the mapping, so that you can see how 
> the surface aligns with the volume.
>
> Sometimes there are header issues, and the origin is not set correctly, 
> resulting in faulty volume-surface alignment.
>
>
> On Aug 1, 2014, at 4:29 PM, "Tang, Yan"  wrote:
>
>> When I open the spec file and mapped the Metric, only the surface was 
>> displayed. The result is in the attachment. How should I do?
>> 
>> From: caret-users-boun...@brainvis.wustl.edu 
>> [caret-users-boun...@brainvis.wustl.edu] on behalf of Donna Dierker 
>> [donna.dier...@sbcglobal.net]
>> Sent: Friday, August 01, 2014 4:02 PM
>> To: Caret, SureFit, and SuMS software users
>> Cc: Tang, Yiyuan
>> Subject: Re: [caret-users] projecting functional MRI to gii surfaces
>>
>> Hmmm.  Sounds like more than a header to me.
>>
>> When you open the spec file you selected when you mapped the data, and 
>> select the output file that is 446kb, what happens?
>>
>> You must make sure you select Metric on the D/C: Overlay/Underlay Surface 
>> menu (primary or secondary, typically).  Else it won't display.
>>
>>
>> On Aug 1, 2014, at 2:47 PM, "Tang, Yan"  wrote:
>>
>>> I am sure that the file exist and the size of file is 446KB. Is It correct?
>>> 
>>> From: caret-users-boun...@brainvis.wustl.edu 
>

Re: [caret-users] projecting functional MRI to gii surfaces

2014-08-05 Thread Donna Dierker
You can do that, but I am not used to seeing this done with the mean 
midthickness surface.  Alex Cohen did something like this when he was using 
resting state functional connectivity to find gradients in a subject's cortical 
networks.  He used the PALS flat or spherical map to get a common grid onto a 
standard mesh.  Once there, he projected that grid onto the individual's 
midthickness surface on the standard mesh.  (Like the grid is getting folded 
back up into the individual's anatomical pattern.)  Then he unprojected the 
points to use as seeds for his analysis.

Do you mind if I ask how the mean midthickness surface comes into play?

Caret's Layers: Borders has options for making grids on the flat map.  People 
make grids on the sphere in matlab.

There are caret_command tools for unprojecting borders.

  caret_command -surface-border-unprojection  

Border files differ from border projection files in that they are points not 
tied to a particular mesh.  The advantage is that you can open the same border 
point on, say, a native and 164k mesh, and it will align with both, if it 
aligns with one (and they are identical except for mesh).  The advantage of 
borderproj files is that they open on multiple configurations - flat, 
midthickness, inflated, etc. - but the price is that you're tied to a mesh.


On Aug 5, 2014, at 10:06 AM, "Tang, Yan"  wrote:

> Thank you for your help. Another problem is how to use the Caret software to 
> generate  regularly spaced Cartesian grids  on the flattened PALS-B12 average 
> surface of the left and right hemispheres.  Can I use this 3-dimensional (3D) 
> stereotactic coordinates from the PALS-B12 average fiducial (midthickness) 
> surface for each grid location to obtain the voxel coordinates (3 × 3 × 3 mm 
> resolution) containing that point in fMRI?
> 
> 
> From: caret-users-boun...@brainvis.wustl.edu 
> [caret-users-boun...@brainvis.wustl.edu] on behalf of Donna Dierker 
> [donna.dier...@sbcglobal.net]
> Sent: Friday, August 01, 2014 5:34 PM
> To: Caret, SureFit, and SuMS software users
> Cc: Tang, Yiyuan
> Subject: Re: [caret-users] projecting functional MRI to gii surfaces
> 
> Push Toolbar: D/C and make sure the primary overlay is Metric.
> 
> Make sure the right column is selected.
> 
> If that check out okay, then I would do:
> 
> File: Open Data File: Volume Functional File
> Load the volume you just mapped
> Switch to volume view and select view All (as opposed to H (horizontal or 
> axial).
> Select D/C and on the page selection drop-down menu, scroll all the way to 
> the bottom
>something like volume surface outline
> Toggle on the fiducial surface used for the mapping, so that you can see how 
> the surface aligns with the volume.
> 
> Sometimes there are header issues, and the origin is not set correctly, 
> resulting in faulty volume-surface alignment.
> 
> 
> On Aug 1, 2014, at 4:29 PM, "Tang, Yan"  wrote:
> 
>> When I open the spec file and mapped the Metric, only the surface was 
>> displayed. The result is in the attachment. How should I do?
>> 
>> From: caret-users-boun...@brainvis.wustl.edu 
>> [caret-users-boun...@brainvis.wustl.edu] on behalf of Donna Dierker 
>> [donna.dier...@sbcglobal.net]
>> Sent: Friday, August 01, 2014 4:02 PM
>> To: Caret, SureFit, and SuMS software users
>> Cc: Tang, Yiyuan
>> Subject: Re: [caret-users] projecting functional MRI to gii surfaces
>> 
>> Hmmm.  Sounds like more than a header to me.
>> 
>> When you open the spec file you selected when you mapped the data, and 
>> select the output file that is 446kb, what happens?
>> 
>> You must make sure you select Metric on the D/C: Overlay/Underlay Surface 
>> menu (primary or secondary, typically).  Else it won't display.
>> 
>> 
>> On Aug 1, 2014, at 2:47 PM, "Tang, Yan"  wrote:
>> 
>>> I am sure that the file exist and the size of file is 446KB. Is It correct?
>>> 
>>> From: caret-users-boun...@brainvis.wustl.edu 
>>> [caret-users-boun...@brainvis.wustl.edu] on behalf of Donna Dierker 
>>> [donna.dier...@sbcglobal.net]
>>> Sent: Friday, August 01, 2014 10:10 AM
>>> To: Caret, SureFit, and SuMS software users
>>> Cc: Tang, Yiyuan
>>> Subject: Re: [caret-users] projecting functional MRI to gii surfaces
>>> 
>>> Hi Yan,
>>> 
>>> Could you use a terminal window or file manager to check whether the file 
>>> exists, and if so, what its size is.
>>> 
>>> We have seen cases before where the file was just a header -- no data.  
>>> Inexplicably, the presence of a non-english character set on the system 
>>> used has caused this sort of trouble.  If there is a system nearby that 
>>> does not have a non-english character set installed, you might see if Caret 
>>> works there.  Or remove any non-english character sets on your system and 
>>> see if it helps.
>>> 
>>> Donna
>>> 
>>> 
>>> On Jul 31, 2014, at 3:52 PM, "Tang, Yan"  wrote:
>>> 
>>

Re: [caret-users] projecting functional MRI to gii surfaces

2014-08-05 Thread Tang, Yan
Thank you for your help. Another problem is how to use the Caret software to 
generate  regularly spaced Cartesian grids  on the flattened PALS-B12 average 
surface of the left and right hemispheres.  Can I use this 3-dimensional (3D) 
stereotactic coordinates from the PALS-B12 average fiducial (midthickness) 
surface for each grid location to obtain the voxel coordinates (3 × 3 × 3 mm 
resolution) containing that point in fMRI?
 

From: caret-users-boun...@brainvis.wustl.edu 
[caret-users-boun...@brainvis.wustl.edu] on behalf of Donna Dierker 
[donna.dier...@sbcglobal.net]
Sent: Friday, August 01, 2014 5:34 PM
To: Caret, SureFit, and SuMS software users
Cc: Tang, Yiyuan
Subject: Re: [caret-users] projecting functional MRI to gii surfaces

Push Toolbar: D/C and make sure the primary overlay is Metric.

Make sure the right column is selected.

If that check out okay, then I would do:

File: Open Data File: Volume Functional File
Load the volume you just mapped
Switch to volume view and select view All (as opposed to H (horizontal or 
axial).
Select D/C and on the page selection drop-down menu, scroll all the way to the 
bottom
something like volume surface outline
Toggle on the fiducial surface used for the mapping, so that you can see how 
the surface aligns with the volume.

Sometimes there are header issues, and the origin is not set correctly, 
resulting in faulty volume-surface alignment.


On Aug 1, 2014, at 4:29 PM, "Tang, Yan"  wrote:

> When I open the spec file and mapped the Metric, only the surface was 
> displayed. The result is in the attachment. How should I do?
> 
> From: caret-users-boun...@brainvis.wustl.edu 
> [caret-users-boun...@brainvis.wustl.edu] on behalf of Donna Dierker 
> [donna.dier...@sbcglobal.net]
> Sent: Friday, August 01, 2014 4:02 PM
> To: Caret, SureFit, and SuMS software users
> Cc: Tang, Yiyuan
> Subject: Re: [caret-users] projecting functional MRI to gii surfaces
>
> Hmmm.  Sounds like more than a header to me.
>
> When you open the spec file you selected when you mapped the data, and select 
> the output file that is 446kb, what happens?
>
> You must make sure you select Metric on the D/C: Overlay/Underlay Surface 
> menu (primary or secondary, typically).  Else it won't display.
>
>
> On Aug 1, 2014, at 2:47 PM, "Tang, Yan"  wrote:
>
>> I am sure that the file exist and the size of file is 446KB. Is It correct?
>> 
>> From: caret-users-boun...@brainvis.wustl.edu 
>> [caret-users-boun...@brainvis.wustl.edu] on behalf of Donna Dierker 
>> [donna.dier...@sbcglobal.net]
>> Sent: Friday, August 01, 2014 10:10 AM
>> To: Caret, SureFit, and SuMS software users
>> Cc: Tang, Yiyuan
>> Subject: Re: [caret-users] projecting functional MRI to gii surfaces
>>
>> Hi Yan,
>>
>> Could you use a terminal window or file manager to check whether the file 
>> exists, and if so, what its size is.
>>
>> We have seen cases before where the file was just a header -- no data.  
>> Inexplicably, the presence of a non-english character set on the system used 
>> has caused this sort of trouble.  If there is a system nearby that does not 
>> have a non-english character set installed, you might see if Caret works 
>> there.  Or remove any non-english character sets on your system and see if 
>> it helps.
>>
>> Donna
>>
>>
>> On Jul 31, 2014, at 3:52 PM, "Tang, Yan"  wrote:
>>
>>> Dear all,
>>>
>>> I used  the Freesurfer_to_fs_LR Pipeline to get  164k fs_LR surface. Now I 
>>> want to map functional volumes to surfaces.
>>> In volume selection page, I choose my file 'ff001_010.nii'; In spec file 
>>> and surface selection page, I choose the file 'study1.L.orig.164k_fs_LR'. I 
>>> get a file 'map_data_0_31_Jul_2014_15_10_13.metric'. But I find I couldn't 
>>> open this file. Which step is wrong? How can you do it?
>>> ___
>>> caret-users mailing list
>>> caret-users@brainvis.wustl.edu
>>> http://brainvis.wustl.edu/mailman/listinfo/caret-users
>>
>>
>> ___
>> caret-users mailing list
>> caret-users@brainvis.wustl.edu
>> http://brainvis.wustl.edu/mailman/listinfo/caret-users
>>
>> ___
>> caret-users mailing list
>> caret-users@brainvis.wustl.edu
>> http://brainvis.wustl.edu/mailman/listinfo/caret-users
>
>
> ___
> caret-users mailing list
> caret-users@brainvis.wustl.edu
> http://brainvis.wustl.edu/mailman/listinfo/caret-users
> ___
> caret-users mailing list
> caret-users@brainvis.wustl.edu
> http://brainvis.wustl.edu/mailman/listinfo/caret-users


___
caret-users mailing list
caret-users@brainvis.wustl.edu
http://brainvis.wustl.edu/mailman/listinfo/caret-users

___
care

[caret-users] Metric file must be provided for the metric mapping

2014-08-05 Thread Donna Dierker
Yesterday, I ran into two arcane issues worth passing on, one of which was 
covered in an earlier thread:

1.  First I tried mapping a paint volume onto an atlas surface, but it turned 
out the volume was scalar values -- not discrete integers as I'd expected for a 
label/ROI volume.  Caret generated paint files that were just headers -- no 
data, but didn't complain.

2.  After realizing these were scalars, I mapped as functional, but the mapping 
algorithm stayed on the paint-related algorithm, as Tony explains here:

http://www.mail-archive.com/caret-users%40brainvis.wustl.edu/msg02887.html

When I tried to map, I got the "Metric file must be provided for the metric 
mapping" error.  Switching to a metric algorithm resolved the error.

No action required; this is in case someone runs into the same issue and 
searches caret-users for answers.
___
caret-users mailing list
caret-users@brainvis.wustl.edu
http://brainvis.wustl.edu/mailman/listinfo/caret-users


Re: [caret-users] projecting functional MRI to gii surfaces

2014-08-05 Thread Donna Dierker
I'm not clear on what you mean by "I want get fMRI time course for surface 
vertices of every subject."

If you just mean how do you scale up -- map that many volumes to all your 
subjects -- then I recommend scripting it and using caret_command.  (Note that 
Workbench, the software that is superseding Caret 5.*, has more robust mapping 
features than caret_command, but I am going to provide the caret_command usage, 
since this is the caret-users list.  There is also a hcp-users list that covers 
workbench.)

Here is the usage for the command that maps volumes onto surfaces:

  caret_command -volume-map-to-surface  
 
 
 
 
 
 
 [-av  average-voxel-neighbor-cube-size (mm)]
 [-bf  brain-fish-max-distance
   brain-fish-splat factor]
 [-g   gaussian-neighbor-cube-size (mm)
   sigma-norm
   sigma-tang
   norm below cutoff (mm)
   norm above cutoff (mm)
   tang-cutoff (mm)]
 [-mv  maximum-voxel-neighbor-cube-size (mm)]
 [-sv  strongest-voxel-neighbor-cube-size (mm)]
 
 Map volume(s) to a surface metric or paint file.
 
 For successful mapping, both the surface and the volume
 must be in the same stereotaxic space.
 
 "algorithm" is one of:
METRIC_AVERAGE_NODES
METRIC_AVERAGE_VOXEL
METRIC_ENCLOSING_VOXEL
METRIC_GAUSSIAN
METRIC_INTERPOLATED_VOXEL
METRIC_MAXIMUM_VOXEL
METRIC_MCW_BRAIN_FISH
METRIC_STRONGEST_VOXEL
PAINT_ENCLOSING_VOXEL
 
 If the input metric or paint file name is not an empty string
  (""), the newly create metric or paint columns will be 
 appended to the file and then written with the output file 
 name.



On Aug 4, 2014, at 3:49 PM, "Tang, Yan"  wrote:

> Thank you for your help.  Now, I meet another problem. I project all 
> functional MRI to 164k fs_LR surface. Every subject have 135 volumes. I want 
> get fMRI time course for surface vertices of every subject. How should I do?
> 
> From: caret-users-boun...@brainvis.wustl.edu 
> [caret-users-boun...@brainvis.wustl.edu] on behalf of Donna Dierker 
> [donna.dier...@sbcglobal.net]
> Sent: Friday, August 01, 2014 5:34 PM
> To: Caret, SureFit, and SuMS software users
> Cc: Tang, Yiyuan
> Subject: Re: [caret-users] projecting functional MRI to gii surfaces
> 
> Push Toolbar: D/C and make sure the primary overlay is Metric.
> 
> Make sure the right column is selected.
> 
> If that check out okay, then I would do:
> 
> File: Open Data File: Volume Functional File
> Load the volume you just mapped
> Switch to volume view and select view All (as opposed to H (horizontal or 
> axial).
> Select D/C and on the page selection drop-down menu, scroll all the way to 
> the bottom
>something like volume surface outline
> Toggle on the fiducial surface used for the mapping, so that you can see how 
> the surface aligns with the volume.
> 
> Sometimes there are header issues, and the origin is not set correctly, 
> resulting in faulty volume-surface alignment.
> 
> 
> On Aug 1, 2014, at 4:29 PM, "Tang, Yan"  wrote:
> 
>> When I open the spec file and mapped the Metric, only the surface was 
>> displayed. The result is in the attachment. How should I do?
>> 
>> From: caret-users-boun...@brainvis.wustl.edu 
>> [caret-users-boun...@brainvis.wustl.edu] on behalf of Donna Dierker 
>> [donna.dier...@sbcglobal.net]
>> Sent: Friday, August 01, 2014 4:02 PM
>> To: Caret, SureFit, and SuMS software users
>> Cc: Tang, Yiyuan
>> Subject: Re: [caret-users] projecting functional MRI to gii surfaces
>> 
>> Hmmm.  Sounds like more than a header to me.
>> 
>> When you open the spec file you selected when you mapped the data, and 
>> select the output file that is 446kb, what happens?
>> 
>> You must make sure you select Metric on the D/C: Overlay/Underlay Surface 
>> menu (primary or secondary, typically).  Else it won't display.
>> 
>> 
>> On Aug 1, 2014, at 2:47 PM, "Tang, Yan"  wrote:
>> 
>>> I am sure that the file exist and the size of file is 446KB. Is It correct?
>>> 
>>> From: caret-users-boun...@brainvis.wustl.edu 
>>> [caret-users-boun...@brainvis.wustl.edu] on behalf of Donna Dierker 
>>> [donna.dier...@sbcglobal.net]
>>> Sent: Friday, August 01, 2014 10:10 AM
>>> To: Caret, SureFit, and SuMS software users
>>> Cc: Tang, Yiyuan
>>> Subject: Re: [caret-users] projecting functional MRI to gii surfaces
>>> 
>>> Hi Yan,
>>> 
>>> Could you use a terminal window or file manager to check whether the file 
>>> exists, and if so, what its size is.
>>> 
>>> We have seen cases before where the file was just a header -- no data.  
>>> Inexplicably, the presence of a non-english character set o