If I follow this, you're eroding one 1x1x1 voxel in hires space and then doing a nearest-neighbor interpolation to EPI space? That's probably going to end up with some "WM" voxels in the grey matter.