Yet in fact, the 'nearest' interpolation is not used and instead, the default 'trilinear' is used (this has important effects especially at GM-WM and WM-CSF interface).
The reason is that the 4th argument of dtiGetValFromFibers() is not the interpMethod, but rather the valName. It should therefore read:
Dear all,
I think that dtiCreateQuenchStats.m has a small bug that need fixing. Line 197 reads:
vals = dtiGetValFromFibers(data, fg, mat,'nearest');
Yet in fact, the 'nearest' interpolation is not used and instead, the default 'trilinear' is used (this has important effects especially at GM-WM and WM-CSF interface). The reason is that the 4th argument of dtiGetValFromFibers() is not the interpMethod, but rather the valName. It should therefore read:
vals = dtiGetValFromFibers(data, fg, mat, [], 'nearest');
Best, Roey