Hi Ariel,
> Our original raw data files have a voxel size of 3.12 x 3.12 x 6, and we
> are setting the write normalize defaults to a voxel size of 3 x 3 x 3.
> When we run through our normalization process the w* files come out with
> voxel sizes of 3x3x3 as we would like them to. However, then when we
> smooth the data, the smoothed images (sw*) come out with a voxel size of
> 3.12x3.12x6, returning to the original raw data size! This is curious,
This is very curious indeed... looking at spm_smooth_ui.m and
spm_smooth.m, it appears to me that the smoothed output gets its
volume info from the input image and hence should have exactly the
same matrix (and vox dims) as the input.
Can you post the results of something like:
in = 'wExample.img'; %%% change this
[pth fnm ext] = fileparts(in);
out = fullfile(pth, ['s' fnm ext]);
spm_smooth(in, out, 12);
Vin = spm_vol(in)
Vout = spm_vol(out)
Vin.mat
Vout.mat
Cheers,
Ged.
|