1 / 32

330 likes | 614 Views

Markov random field: A brief introduction (2). Tzu-Cheng Jen Institute of Electronics, NCTU 2007-07-25. Outline. Markov random field: Review Edge-preserving regularization in image processing. Markov random field: Review. Prior knowledge.

Download Presentation
## Markov random field: A brief introduction (2)

**An Image/Link below is provided (as is) to download presentation**
Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.
Content is provided to you AS IS for your information and personal use only.
Download presentation by click this link.
While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.
During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

**Markov random field: A brief introduction (2)**Tzu-Cheng Jen Institute of Electronics, NCTU 2007-07-25**Outline**• Markov random field: Review • Edge-preserving regularization in image processing**Prior knowledge**• In order to explain the concept of the MRF, we first introduce following definition: 1. i: Site (Pixel) 2. fi: The value at site i (Intensity) 3. S: Set of sites (Image) 4. Ni: The neighboring site of i (1st order neighborhood of i: f2, f4, f5, f7 ) 5. Ci: Clique, the subset of S and the element in this subset must be neighboring A 3x3 imagined image**Markov random field (MRF)**• View the 2D image f as the collection of the random variables (Random field) • A random field is said to be Markov random field if it satisfies following properties Red: Neighboring site**Gibbs random field (GRF) and Gibbs**distribution • A random field is said to be a Gibbs random field if and only if its configuration f obeys Gibbs distribution, that is: A 3x3 imagined image U(f): Energy function; T: Temperature Vi(f): Clique potential**Markov-Gibbs equivalence**• Hammersley-Clifford theorem: A random field F is an MRF if and only if F is a GRF Red: Neighboring site**MAP formulation for denoising problem**Noisy signal d denoised signal f d = f + N(0, σ)**MAP formulation for denoising problem**• A signal denoising problem could be modeled as the MAP estimation problem, that is, (Observation model) (Prior model)**MAP formulation for denoising problem**• Assume the observation is the true signal plus the independent Gaussian noise, that is • Assume the unknown data f is MRF, the prior model is:**MAP formulation for denoising problem**• Substitute above information into the MAP estimator, we could get: Observation model (Similarity measure) Prior model (Reconstruction constrain, Regularization)**The solver of the optimization problem: Gradient descent**algorithm**Simulation results for denoising problem**Simulation result Simulation result Processed profiles are blurred !**Discussion for the phenomenon of blur (1)**• From the potential function point of view: Energy 1st derivative Quadratic function is used as potential function g=x2 Simulation result**Discussion for the phenomenon of blur (2)**• From the optimization process point of view (gradient descent algorithm): Update equation of gradient descent:**Edge-preserving regularization**• S. Geman and D. Geman, "Stochastic relaxation, Gibbs distributions, and the Bayesian restoration of images," IEEE Trans. Pattern Anal. Mach. Intell, 6, 721-741, 1984. • S.Z. Li, “On Discontinuity-adaptive smoothness priors in computer vision,” IEEE Trans. Pattern Anal. Mach. Intell, June, 1995. • Pierre Charbonnier et al, “Deterministic edge preserving regularization in computed imaging,” IEEE Trans. Image Processing, Feb, 1997. • S.Z. Li, “Markov random field modeling in computer vision,“ Springer, 1995**MRF with pixel process and line process (Geman and Geman)**Lattice of pixel site: SP Labeling value: fip (real value) Lattice of line site: SE Labeling value: fii’E (only 0 or 1)**MRF with pixel process and line process (Geman and Geman)**• Based on the concept of line process, we could modify the original restoration problem as follows: ? Goal: Find realizations fp and fE such that edge-preserving regularization could be achieved**MRF with pixel process and line process (Geman and Geman)**• Define the prior: • Substitute above information into the MAP estimator, we could get: The above optimization problem is a combination of real and combinatorial problem**MRF with pixel process and line process (Geman and Geman)**• Blake and Zisserman covert previous restoration problem into real minimization problem by introducing truncated quadratic function as potential function Energy Truncated quadratic function**MRF with pixel process and line process (Geman and Geman)**• Simulation results Original image Degraded image Restoration result (1000 iterations)**MRF with pixel process and line process (Geman and Geman)**• Simulation results Original image Degraded image Restoration result (100 iterations) Restoration result (1000 iterations)**Discontinuity-adaptive regularization (S. Z. Li)**• Revisit the gradient descent algorithm Adjust it adaptively ! Derivative or compensator Weight or interaction function**Discontinuity-adaptive regularization (S. Z. Li)**• For edge-preserving regularization, interaction function hr should satisfy following property:**Discontinuity-adaptive regularization (S. Z. Li)**• Possible choices for interaction function hr**Discontinuity-adaptive regularization (S. Z. Li)**• Simulation results (1D)**Discontinuity-adaptive regularization (S. Z. Li)**• Simulation results (2D) Original image Edge-preserving restoration Restoration without preserving edge**Discontinuity-adaptive regularization (Pierre Charbonnier**et al ) • Pierre Charbonnier et al impose following conditions on potential function φfor edge preserving regularization**Discontinuity-adaptive regularization (Pierre Charbonnier**et al ) • Possible choices for potential function φ**Other related techniques for**edge-preserving regularization • P. Perona and J. Malik, “Scale-space and edge detection using anisotropic diffusion,” IEEE Trans. Pattern Anal. Mach. Intell, July, 1990. Dropping observation model (w=1) when evaluating f**Other related techniques for**edge-preserving regularization • L.I. Rudin, S. Osher, E. Fatemi (1992): Nonlinear Total Variation Based Noise Removal Algorithms, Physica D, 60(1-4), 259-268. Energy Replace the quadratic potential function with absolute value function Quadratic function versus Absolute value function

More Related