390 likes | 417 Views
Lecture 7. Patch based methods: nonlocal means, BM3D, K-SVD, data-driven (tight) frame. Outline. Nonlocal methods Nonlocal means BM3D Dictionary learning methods K-SVD Data-driven tight frame Low rank models. Nonlocal Methods. Local means and BM3D.
E N D
Lecture 7 Patch based methods: nonlocal means, BM3D, K-SVD, data-driven (tight) frame
Outline • Nonlocal methods • Nonlocal means • BM3D • Dictionary learning methods • K-SVD • Data-driven tight frame • Low rank models
Nonlocal Methods Local means and BM3D
Denoising: Gaussian Smoothing Revisited *h = * *h *h h
Perona-Malik: Anisotropic Filtering • Edges smooth only along edges • “Smooth” regions smooth isotropically gradient
Ideal: Spatially Adaptive Smoothing • Non uniform smoothingDepending on image content: • Smooth where possible • Preserve fine details How? *h3 *h2 *h1
Slides taken from Sylvain Paris, Siggraph 2007 Gaussian Smoothing * output input * * Same Gaussian kernel everywhere Averages across edges blur
Slides taken from Sylvain Paris, Siggraph 2007 Bilateral Filtering * output input * * Kernel shape depends on image content Avoids averaging across edges
Nonlocal Means (NLM): Motivation • Assume a static scene • Consider multiple images at different times • The signal remains constant • varies over time with zero mean
Nonlocal Means (NLM): Motivation Average multiple images over time
Nonlocal Means (NLM): Motivation Average multiple images over time
Nonlocal Means (NLM): Motivation Average multiple images over time
Redundancy in natural images Glasneret al. (2009)
Single image “time-like” denoising Unfortunately, patches are not exactly the same simple averaging just won’t work
NLM • Buades, Cole and Morel (2005) Use a weighted average based on similarity
From Bilateral Filters to NLM Patch similarity
Performance Evaluation Gaussian Smoothing Anisotropic Filtering Bilateral Filtering NLM Windowed Weiner Hard WT Soft WT Buadeset al. (2005)
Variational Formulation of NLM • Gilboa and Osher(2008) • Nonlocal (partial) derivative • Nonlocal gradient: • Inner product: • With such inner product, we can define divergence • And Laplacian
Variational Formulation of NLM • Gilboa and Osher (2008) • Nonlocal TV • Nonlocal ROF • Nonlocal TV-
Variational Formulation of NLM • Gilboa and Osher (2008)
Application in Surface Denoising • Dong et al. (2008)
Application in Surface Denoising • Dong et al. (2008)
Application in 4D CT Reconstruction • Tian et al. (2011)
What’s Next? • The idea of grouping sounds good reduces mixing • Denoise = “extract the common (the signal)” • NLM: common = weighted average • Can a sparser representation do better?
BM3D: Dabov et al. (2007) • Block Matching 3D collaborative filtering • Group patches with similar local structure (BM) • Jointly denoise each group (3D) • Smart Fusion of multiple estimates
BM3D: Dabov et al. (2007) Single BM3D Estimate Block matching Inverse 3D transform Filter / thresholding R R R 3D grouping Denoised 3D group 3D transform
BM3D: Dabov et al. (2007) • For every noisy reference block: • Calculate SSD between noisy blocks • If SSD<thr add to group
BM3D: Dabov et al. (2007) • 3D transform 3D transform Image Patch Domain Sparse Domain
BM3D: Dabov et al. (2007) • Collaborative filtering • Use hard thresholding or Wiener filter • Each patch in the group gets a denoised estimate • Unlike NLM – where only central pixel in reference patch got an estimate Filter / Thresholding Noisy patches Denoised Patches R R
BM3D: Dabov et al. (2007) Collaborative filtering R R thr thr Multiple BM3D Estimate R R R R R Collaborative filtering R
BM3D: Dabov et al. (2007) R R t t ? Fusion R R R R R R
BM3D: Dabov et al. (2007) • Each pixel gets multiple estimates from different groups • Naive approach Average all estimates of each pixel …. not all estimates are as good • SuggestionGive higher weight to more reliable estimates
BM3D: Dabov et al. (2007) • Give each estimate a weight according to denoising quality of its group • Quality = Sparsity induced by the denoising Hard thresholding Weiner filtering
BM3D: Dabov et al. (2007) • Noise may result in poor matching Degrades de-noising performance • Improvements: • Match using a smoothed version of the image • Perform BM3D in 2 phases: • Basic BM3D estimate improved 3D groups • Final BM3D • Variational formulation: BM3D frame Danielyan, Katkovnikand Egiazarian. BM3D frames and variational image deblurring. IEEE TIP. 21(4):1715-28. 2012.
BM3D: Dabov et al. (2007) • Results
Dictionary Learning K-SVD, Data-Driven Tight Frame
Dictionary Selection • Which D to use? • A fixed set of basis: • Steerable wavelet • Contourlet • DCT Basis • …… • Data adaptive dictionary – learn from data • K-SVD (-norm) http://bicmr.pku.edu.cn/~dongbin/Teaching_files/%E5%9B%BE%E5%83%8F%E5%A4%84%E7%90%86%E4%B8%AD%E7%9A%84%E6%95%B0%E5%AD%A6%E6%96%B9%E6%B3%95-18-19/ClassMaterials/K-SVD.pptx • Data-Driven Tight Frame (DDTF) http://bicmr.pku.edu.cn/~dongbin/Teaching_files/%E5%9B%BE%E5%83%8F%E5%A4%84%E7%90%86%E4%B8%AD%E7%9A%84%E6%95%B0%E5%AD%A6%E6%96%B9%E6%B3%95-18-19/ClassMaterials/DDTF.pdf