2009/8/25 Feng-Xian Yan <mimigirln@gmail.com>

Hi,

   First, I use mri_surf2surf (as follows) to map our subjects to fsaverage.

mri_surf2surf –srcsubject subject --sval-tal-xyz white --hemi lh --trgsubject fsaverage --tval-xyz --tval lh.subject.tal

 

Next step, I use mris_convert to show the cortical thickness value of this subject on this MNI coordinates. But, this has an error as follows.
mris_convert -c subject/surf/lh.thickness fsaverage /surf/lh.white  lh.subject.thickness.asc
 

ERROR: number of vertices in subject/surf/lh.thickness does not match surface  (154962,163842)

I know what’s happen, but I don’t know how to solve.

I want to obtain the asc file containing the same vertices as fsaverage and with the cortical thickness value and their coordinates in the file. Could you teach me how to obtain the match fsaverage surface file?

Thank you in advance!

Feng-Xian

 

 



2009/8/20 Bruce Fischl <fischl@nmr.mgh.harvard.edu>

you can use mri_surf2surf to map them all the fsaverage, then mris_convert to convert each one to asciii.

cheers
Bruce

On Thu, 20 Aug 2009, Feng-Xian Yan wrote:

Hi,

 But I want to obtain the asc file or other files containing the same
vertices and with the cortical thickness value and their coordinates in the
file.

Could some ways to obtain?



Thank you in advance!



Feng-Xian




2009/8/20 Doug Greve <greve@nmr.mgh.harvard.edu>


Yes, this is what our group analysis tools do. You can find tutorials
on the wiki. You'll probably start with mris_preproc.

doug




On Thu, 20 Aug 2009, Feng-Xian Yan wrote:

Hi,

 Thank you for your responses before.



 I have another problem want to solve. I check the number of vertex for
all
subjects, but I found that the number of vertex for each subjects and each
hemispheres is different. Because I want to compare cortical thickness
using
vertex-by-vertex, I would to obtain the same vertex number.

Therefore, have something to do to obtain the same vertex number, i.e.
template the cortical thickness at each vertex for each subject on the
fsaverage surface, or other methods?



How can I do? Would you get me a hand?



Thank you in advance!



Feng-Xian




2009/8/18 Feng-Xian Yan <mimigirln@gmail.com>


Thank s for your response.

So, the coordinates are MNI coordinates, not Talairach coordinates.



Thank you in advance.



Feng-Xian



2009/8/18 Douglas N Greve <greve@nmr.mgh.harvard.edu>

The values in the asc file are in MNI305 space. You can see this in
tksurfer with View->Configure->Information and select the mni305 button.

Feng-Xian Yan wrote:


Thank you for your response.

So, you mean that I have to do "mris_convert" two times. One for
lh.white.tal.asc, another for lh.thickness.asc. And the Tal coordinates
in
lh.white.tal.asc file corresponds to the cortical thickness value in
lh.thickness.asc file.


But, the Tal coordinates in lh.white.tal.asc is difference from the
result in the tksurfer tool box. That is to say, when I open this
subject by
tksurfer, I look at the coordinate of the vertex Tal in the tksurfer
tool
box, and the coordinate is difference from the coordinates that obtain
from
the command ?mris_convert ?c lh.thickness lh.white lh.thickness.asc?.


What do have some problems?


Thank you in advance!


Feng-Xian


 2009/8/18 Douglas N Greve <greve@nmr.mgh.harvard.edu <mailto:
greve@nmr.mgh.harvard.edu>>

 You have convert the thickness file (lh.thickness) to ascii like
 you did
 with the surface (with mris_convert). You will then have two
 files, and
 the second file will have the thickness values.

 doug

 Feng-Xian Yan wrote:
>
> Hi,
>
> Sorry, I don't know what you mean. I know each lines is a vertex and
> represent a Tal coordinate, but each lines only contain xyz, but not
> contain cortical thickness value of these coordinates. That is
 to say,
> the following is
>
> Tal x Tal y Tal z
>
> -13.424789 -125.068245 12.671550 0
>
> What?s wrong do I think?
>
> Thank you in advance.
>
> Feng-Xian
>
> 2009/8/18 Douglas N Greve <greve@nmr.mgh.harvard.edu
 <mailto:greve@nmr.mgh.harvard.edu>
 > <mailto:greve@nmr.mgh.harvard.edu
 <mailto:greve@nmr.mgh.harvard.edu>>>
>
>     Just convert the lh.thickness to ascii as well. Each line of the
>     xyz coords corresponds to a line in the thickness ascii
 (each line
>     is a vertex).
>
>
>     doug
>
>     Feng-Xian Yan wrote:
>
>
>         Hi,
>
>         I try it again with you recommend me, but the produced asc
>         file (lh.white.tal.asc) only contain x, y, z
 coordinates. (The
>         following shows that.)
>
>
>         #!ascii version of lh.white.tal
>
>         155006 310008
>
>         -13.424789 -125.068245 12.671550 0
>
>         -13.939213 -125.063232 12.714314 0
>
>         -14.861158 -125.268639 12.404987 0
>
>
>
>         But, I want to have cortical thickness value with Tal
>         coordinates. This file doesn?t satisfy our demand.
>
>
>         Our demand should contain two points.
>
>         1) Each vertex across the mantle was represented with an
>         estimated value of cortical thickness and these estimated
>         values of cortical thickness should have their coordinates.
>
>         2) We want to analysis the cortical thickness difference of
>         three groups, so we want to the cortical thickness value of
>         each vertex across each subjects should be matched. That
 is to
>         say, the cortical thickness should superimposed on a
 template
>         surface like fsaverage.
>
>         Therefore, how do we do to obtain the asc file matching our
>         demands?
>
>
>         Thank you in advance.
>
>
>         Feng-Xian
>
>
>
>
>         2009/8/14 Douglas N Greve <greve@nmr.mgh.harvard.edu
 <mailto:greve@nmr.mgh.harvard.edu>
>         <mailto:greve@nmr.mgh.harvard.edu
 <mailto:greve@nmr.mgh.harvard.edu>>
>         <mailto:greve@nmr.mgh.harvard.edu
 <mailto:greve@nmr.mgh.harvard.edu>
>         <mailto:greve@nmr.mgh.harvard.edu
 <mailto:greve@nmr.mgh.harvard.edu>>>>
>
>
>         Try these commands:
>
>         mri_surf2surf --s subject --hemi lh --sval-tal-xyz white
>         --tval-xyz --tval lh.white.tal
>         mris_convert lh.white.tal lh.white.tal.asc
>
>         doug
>
>         Feng-Xian Yan wrote:
>
>
>         Hi,
>
>         I have some problems want to understand.
>
>         (1) I want to analysis the cortical thickness difference of
>         three groups. But, three groups (one factor/three levels) of
>         Freesurfer is not satisfied with our demand. Because we want
>         to get a thickness difference map that satisfied with
 the many
>         voxels containing A group > B group > C group and A
 group < B
>         group < C group, Freesurfer can?t do that. So, we want
 to get
>         this map by myself and use Matlab or other tool to obtain
>         that. And, now I want to obtain a file that containing
 the Tal
>         coordinate and the thickness value of this coordinate
 for all
>         subjects.
>
>         Would you tell me how to get the file that containing
 the Tal
>         coordinate and the thickness value of this coordinate?
>
>         (2) Before, I want to use the following command to get the
>         thickness value.
>
>         The command line is "mri_surf2surf --srcsubject noise0.5_0
>         --srcsurfval thickness --sfmt surface --trgsubject
 noise0.5_0
>         --trgsurfval noise0.5_0_tal.xfm --tfmt surface --hemi lh
>         --sval-tal-xyz orig --tval-xyz", and then i type
 ?mris_convert
>         lh.noise0.5_0_tal.xfm lh.white lh.thickness.asc?.
>
>         But, the result file(lh.thickness.asc) does not contain the
>         thickness value of this coordinate.
>
>         So, I think these commands have some error.
>
>         (3) Use the command ?mris_convert ?c lh.thickness lh.white
>         lh.thickness.asc? to obtain the thickness value of the RAS
>         coordinates. But, I?m confused about that if I change
 lh.white
>         to lh.pail or lh.inflated.
>         What?s difference for the three?
>
>         (4) And, I ask some problem before.
>
>         Everyone said the result file(lh.thickness.asc) is using RAS
>         coordinates through ?mris_convert ?c lh.thickness lh.white
>         lh.thickness.asc?. But, I have confused about this. When I
>         open this subject by tksurfer, I look at the coordinate
 of the
>         vertex RAS in the tksurfer tool box, and the coordinate is
>         difference from the coordinates that obtain from the command
>         ?mris_convert ?c lh.thickness lh.white lh.thickness.asc?.
>         So, which one is right?
>
>         Thank you in advance!
>
>         Feng-Xian
>

 ------------------------------------------------------------------------
>
>         _______________________________________________
>         Freesurfer mailing list
>         Freesurfer@nmr.mgh.harvard.edu
 <mailto:Freesurfer@nmr.mgh.harvard.edu>
>         <mailto:Freesurfer@nmr.mgh.harvard.edu
 <mailto:Freesurfer@nmr.mgh.harvard.edu>>
>         <mailto:Freesurfer@nmr.mgh.harvard.edu
 <mailto:Freesurfer@nmr.mgh.harvard.edu>
>         <mailto:Freesurfer@nmr.mgh.harvard.edu
 <mailto:Freesurfer@nmr.mgh.harvard.edu>>>
>
>
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
>
>
>         -- Douglas N. Greve, Ph.D.
>         MGH-NMR Center
>         greve@nmr.mgh.harvard.edu
 <mailto:greve@nmr.mgh.harvard.edu>
 <mailto:greve@nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu

>         <mailto:greve@nmr.mgh.harvard.edu
 <mailto:greve@nmr.mgh.harvard.edu>
>         <mailto:greve@nmr.mgh.harvard.edu
  <mailto:greve@nmr.mgh.harvard.edu>>>
>
>         Phone Number: 617-724-2358 Fax: 617-726-7422
>
>         In order to help us help you, please follow the steps in:
>         surfer.nmr.mgh.harvard.edu/fswiki/BugReporting
 <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting>
>         <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting>
>         <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting>
>
>
>
>
>     --
>     Douglas N. Greve, Ph.D.
>     MGH-NMR Center
>     greve@nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu>
 <mailto:greve@nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu

>     Phone Number: 617-724-2358 Fax: 617-726-7422
>
>     In order to help us help you, please follow the steps in:
>     surfer.nmr.mgh.harvard.edu/fswiki/BugReporting
 <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting>
>     <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting>
>
>
>

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

 In order to help us help you, please follow the steps in:
 surfer.nmr.mgh.harvard.edu/fswiki/BugReporting
 <http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting>




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

In order to help us help you, please follow the steps in:
surfer.nmr.mgh.harvard.edu/fswiki/BugReporting






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

In order to help us help you, please follow the steps in:
surfer.nmr.mgh.harvard.edu/fswiki/BugReporting