Edge Aware Anisotropic Diffusion for 3D Scalar Data

Fig. 1: The left half of the figure demonstrates the consistency in smoothing of our method compared to the existing method. The right half of the figure demonstrates the de-noising capabilities of our method. All the images from (a-c) were obtained byrenderingan iso-surface of 153. (a) Diffused with an existing diffusion model proposed by Krissian et al. [20] with k = 40, and100 iterations (b) The original Sheep’s heart data. (c) Diffused with our method with σ = 1 and the same number of iterations. The yellow circle indicates aregionwheretheiso-surfacehasbothhighandmediumrangegradient magnitude, and the blue circle marks a region where the gradient magnitude is much lower. Note the inconsistent smoothing in (a) inside the yellow circle. (d) The tooth data contaminated with Poisson noise (SNR=12.8) (e)Theoriginaltoothdata(f)Diffusedwithourmethod(SNR=25.76) withσ = 1 and 25 iterations. We used the exact same transfer function to render all the images in(d-f). Abstract—Inthispaperwepresentanovelanisotropicdiffusionmodel targeted for 3D scalar field data. Our model preserves material boundaries as well as fine tubular structures while noise is smoothed out. One of the major novelties is the use of the directional second derivative to define material boundaries instead of the gradient magnitude for thresholding. This results in a diffusion model that has much lower sensitivity to the diffusion parameter and smoothes material boundaries consistently compared to gradient magnitude based techniques. We empirically analyze the stability and convergence of the proposed diffusion and demonstrate its de-noising capabilities for both analytic and real data. We also discuss applications in the context of volume rendering