diffusion_tensors function

Compute field of diffusion tensors for edge-preserving smoothing.