Dear FreeSurfers,

Please can you explain how mri_surf2surf's --cortex option for
smoothing is implemented?

I assumed it would be along the following lines:
 - mask the data to zero-out the medial wall
 - smooth this masked data
 - smooth the mask itself
 - divide the smoothed data by the smoothed mask
 - (perhaps) re-mask the data

But I've tried doing that long-hand (i.e. calling mri_surf2surf
without the --cortex flag for each smoothing, and using l
load_mgh and save_mgh in MATLAB for the rest), and it doesn't seem to
quite match the results of the --cortex flag.

Now I'm guessing that it instead might do something like the
following. For each connected-neighbour iteration of the smoothing,
only include the neighbour in the average if it is inside the cortex
(not medial wall). However, I can't easily/lazily reimplement this to
check, so could someone reveal the truth? ;-)

Many thanks,
Ged
_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

Reply via email to