Hello,
I am new freesurfer user, I need to have the labels of parcellation on
analyze format. I did recon all command on my subject, I see the very
beautiful parcellation on freesurfer. After reading all the
question/response on the freesurfer website, I write my own txt file
describing the 35 labels: I just copy paste the data obtained from the
variable colortable.struct_name and colortable.table: the output of
read_annotation('rh.aparc.annot'). It seems working correctly, because the
mri_annotation2label creates 35 files of all the labels. To check some of
them I used tkmedit chuv44 brain.mgz and I loaded one label, it seems good.
Is this correct?
But when I tried to export this labels to analyze format (readable with
mricro for example), I use mri_label2vol:
Mri_label2vol -label ./subects/chuv44/mylab/rh.superiorfrontal.label -temp
$SUBJECTS_DIR/chuv44/mri/orig/chuv44.img -o ./subjects/chuv44_rh_one.img
I had a warning saying that the mat file is missing, and the result
orientation is not good.
I then open my analyze file with v=spm_vol('chuv44.img') and copy the v.mat
as the mat file freesurfer asked for:
v.mat
ans =
-2 0 0 129
0 2 0 -129
0 0 2 -69
0 0 0 1
The command line Mri_label2vol -label
./subects/chuv44/mylab/rh.superiorfrontal.label -temp
$SUBJECTS_DIR/chuv44/mri/orig/chuv44.img -o ./subjects/chuv44_rh_one.img
works then correctly and write in the end:
Analyze output matrix
-2 0 0 129
0 2 0 -129
0 0 2 -69
0 0 0 1
---------------------
INFO: set hdr.hist.orient to -1
But when I open the output file chuv44_rh_one.img (analyze) the format is
not the same as my image: the axial slice is coronal and vice versa.
First question: Is the staff I did with my home-made txt labels_file is
correct?
And do you know what I should do to obtain my label mask in correct
orientation?
Thanks in advance
Leila Cammoun