Computation of normal vectors of discrete 3D objects: application to natural snow images from X-ray tomography