The central differences method does only take 6 of the surrounding 26 voxels into account.
The so called Sobel operator applies the central differences method to the entire neighborhood, resulting in a smoother change of the normal across voxel boundaries:
Left: Central Differences, Right: Sobel Operator
