Print

Print


Dear List,

I've noticed that the same image smoothed with the same FWHM kernel yields
subtly different values.

Comparing the code for spm_smooth.m between the two versions yields the
following difference around line 95:


in SPM5:
   x  = exp(-(x).^2/(2*(s(1)).^2));
   y  = exp(-(y).^2/(2*(s(2)).^2));
   z  = exp(-(z).^2/(2*(s(3)).^2));


in SPM8:
   x = spm_smoothkern(s(1),x,1);
   y = spm_smoothkern(s(2),y,1);
   z = spm_smoothkern(s(3),z,1);

I'm assuming that the SPM8 is the favored code at this point; just curious
why it changed (i.e., although spm_smoothkern is part of the SPM5 package,
it's not referenced in the SPM5 version of spm_smooth.  Couldn't find an
explanation in the SPM8 manual or on the list previously.

Cheers!

rob

--
Robert J Ellis, PhD
  Postdoctoral Research Fellow in Neurology
  Beth Israel Deaconess Medical Center
  Harvard Medical School