Adaptive Forward-Backward Diffusion Flows in Image Processing
Spatial regularization in the diffusion coefficient alters discontinuities in a synthetic corner image. (a) Original synthetic image of size 31 x 31, a square (gray value = 160) at the bottom right corner with uniform background (gray value = 219). (b) Input image obtained by adding Gaussian std = 30 to the original image. This noisy image is used as the initial value u0 for the nonlinear PDEs with C1 diffusion coefficient and K = 20. Results of PMADE (1) with 20 iterations in (left) image (right) surface format (c), and GRADE (3) with 20 iterations in (left) image (right) surface format (d). The intersection of red dotted lines indicate the exact corner location of the square.
Diffusion process for a simple synthetic image. (a) Original synthetic image of size 31 × 31, a square (2 × 2, gray value = 1) at the center with uniform background(gray value = 0). (b) Input image obtained by adding Gaussian noise !n = 30 to the original image. This noisy image is used as the initial value u0. (c) Diffusion coefficient C1 in (2), with K = 20. This acts as a discontinuity detector and stops the diffusion spread across edges. (d) Flux function C1(|grad u| ) · |grad u| . (e) Result of heat equation with 20 iterations in (left) image (right) surface format. (f) Result of PMADE equation (1) with 20 iterations in (left) image (right) surface format. The white dotted lines indicate the influence region at the center.
Inverse mollification function when combined with power growth numerator can stop diffusion across edges. Solution of the PDE (*) on noisy synthetic Brain image (noise standard deviation std = 30) with power growth phi(|grad u|)= |grad u|^p, p = 1, 2, 3, 4, 5 (left to right) without inverse mollification (a) g = 0, and with (b) g (x) = x^2 . In both cases we used K = 10^−4 and terminal time 100. It is clear visually that the inverse mollification has an effect in keeping homogenous regions separated by strong edges and avoids leakage.