This is a quirk of the wfile format (just one of the reasons to move away from it:). The wfile is "compressed" in that values that are 0 are not represented. The indices of the values that are there are given in the 2nd output arugment to read_wfile. But it is much easier to just:
mri_surf2surf --srcsubject sp01t7 --hemi lh --srcsurfval thickness --src_type curv --trgsubject sp01MULTt1 --trgsurfval ../../sp01MULTt1/surf/t7tot1-lh.thickness.mgh
doug
On Mon, 7 May 2007, David Perlman wrote:
I think I must be misunderstanding something about the way mri_surf2surf works. I thought it was supposed to resample onto the target mesh, which I would think would mean that resampling two different data files onto the same target mesh would give outputs with the same numbers of nodes. But when I do this and then read the resampled data into Matlab, the nodes do not match up. Here are the commands I used:
mri_surf2surf --srcsubject sp01t7 --hemi lh --srcsurfval thickness --src_type curv --trgsubject sp01MULTt1 --trgsurfval ../../sp01MULTt1/ surf/t7tot1-lh.thickness.w --trg_type w
mri_surf2surf --srcsubject sp01t6 --hemi lh --srcsurfval thickness --src_type curv --trgsubject sp01MULTt1 --trgsurfval ../../sp01MULTt1/ surf/t6tot1-lh.thickness.w --trg_type w
And here's what I see in Matlab:
[vt6to1,ft6to1]=read_wfile('t6tot1-lh.thickness.w'); [vt7to1,ft7to1]=read_wfile('t7tot1-lh.thickness.w'); size(vt6to1)
ans =
124463 1size(vt7to1)
ans =
124134 1What am I missing here?
-- -dave---------------------------------------------------------------- Science arose from poetry... when times change the two can meet again on a higher level as friends. -Göthe
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer