Next: Conclusions
Up: Regularization Schemes
Previous: Gradient magnitude and Laplacian
Until now we have used a roughener as the regularization operator to take into account the prior knowledge that letters are piecewise constant. Turning it off at the edges and on away from the edges.
To do the opposite also makes sense, using a smoothing operator at the edges and turning it off away from the edges doesn't make the model smoother away from the edges but increases the sharpness of the model at the edges.
The gradient magnitude with triangular smoothing edgepreserving regularization fitting goal was set
following the nonlinear iterations:
starting with , at the k^{th} iteration the algorithm solves
 

 (5) 
where
 
(6) 
is a nonstationary convolution matrix, is the result of the k^{th} nonlinear iteration, is the (k1)^{th} diagonal weighting operators, is the identity matrix, is triangular smoother operator,
and the scalar balances the relative importance of the data and model residuals.
comp_images_smooth_2d
Figure 8 A) Original image, B) Deblurred image using LS with the edgepreserving regularization gradient magnitude with triangular smoothing
comp_graph_2d_smooth
Figure 9 Comparison between Figures A and B; A) Slice y=229 and B) Slice x=229of.
Figures and show what we were expecting, i.e. the edges look sharp, but away from the edges the noise dominates. However, we can see in figure some details that we are not visible in Figures , , and . Looking that a combination of the styling goals in equations (3) (smoothing away from the edges) and (5) (roughening at the edges) should give a good result.
Next: Conclusions
Up: Regularization Schemes
Previous: Gradient magnitude and Laplacian
Stanford Exploration Project
10/14/2003