External Email - Use Caution
Dear Freesurfer Expert: First, I got the lh.parsorbitalis.label according to the command:mri_annotation2label --subject fsaverage --hemi lh --outdir /Users/qianqianyu/Desktop/fsl_bold/label_fsaverageNow, I want to make the lh.parsorbitalis.label to fsaverage surface file (nii.gz). How?
I did the conversion myself using the following method, but an error occurred. The structure of lh.parsorbitalis.label is as follows:
I think the first column is coordinate information. So, I use the following commands to convert in Matlab.
a=[116,255,256,735,...3806,...];%All values in the first column are assigned to the matrix, a.
n=1:163842; %cause the vol matrix of fsaverage lh is1*163842
t=zeros(1,length(n));
t(1,a)=1;%value of lh.parsorbitalis area were 1
MRIread('/Users/qianqianyu/Desktop/Project/Sess06/bold/my-syn.sm04.lh/sem-v-fj/sig.nii');%read a fsaverage image
ans.vol=t;
err = MRIwrite(ans,'lh.parsorbitalis.mask.nii.gz')
But, when I check the lh.parsorbitalis.mask.nii.gz in Freeview, I find it is not consistent with the lh.parsorbitalis.label
I am now very very anxious to know how to get the right lh.parsorbitalis.mask.nii.gz!Looking forward to your reply!
Best wishes!
Qianqian Yu
Use mri_label2label. Run it with --help to get examples
On 12/4/19 1:58 AM, 于倩倩 wrote:
External Email - Use Caution
Dear Freesurfer Expert: First, I got the lh.parsorbitalis.label according to the command: mri_annotation2label --subject fsaverage --hemi lh --outdir /Users/qianqianyu/Desktop/fsl_bold/label_fsaverage Now, I want to make the lh.parsorbitalis.label to fsaverage surface file (nii.gz). How?
I did the conversion myself using the following method, but an error occurred. The structure of lh.parsorbitalis.label is as follows: I think the first column is coordinate information. So, I use the following commands to convert in Matlab.
a=[116,255,256,735,...3806,...];%All values in the first column are assigned to the matrix, a.
n=1:163842; %cause the vol matrix of fsaverage lh is1*163842
t=zeros(1,length(n));
t(1,a)=1;%value of lh.parsorbitalis area were 1
MRIread('/Users/qianqianyu/Desktop/Project/Sess06/bold/my-syn.sm04.lh/sem-v-fj/sig.nii');%read a fsaverage image
ans.vol=t;
err = MRIwrite(ans,'lh.parsorbitalis.mask.nii.gz')
But, when I check the lh.parsorbitalis.mask.nii.gz in Freeview, I find it is not consistent with the lh.parsorbitalis.label
I am now very very anxious to know how to get the right lh.parsorbitalis.mask.nii.gz!Looking forward to your reply!
Best wishes!
Qianqian Yu
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
freesurfer@nmr.mgh.harvard.edu