SlideShare a Scribd company logo
Normalized averaging
      using adaptive applicability functions

                            with applications in image reconstruction
                            from sparsely and randomly sampled data




      Presented at SCIA 2003

      Tuan Q. Pham and Lucas J. van Vliet

      July 17, 2009

                                                                        1




Pattern Recognition Group
Overview


• Normalized averaging
• Local structure adaptive filtering
• Experimental results
• Comparison with diffusion-based image inpainting
• Directions for further research




July 17, 2009                                        2
Normalized averaging (Knutsson’93)
•Weighted average filtering:                                                 r = s*a

•Normalized averaging = weighted average + signal/certainty principle:
        •each signal s is associated with a certainty c
        •s & c have to be processed separately                               ( s . c) ∗ a
                                                                          r=
                                                                                c∗a
                where s :signal, c :certainty, a :filter, r :result, * :convolution




                input with 10%          Gaussian smoothing          NA with Gaussian
                 original pixels              (σ = 1)              applicability (σ = 1)
July 17, 2009                                                                               3
Normalized averaging: An example
Reconstruction from 10% random pixels




Nearest neighbor interpolation   NA with adaptive applicability
July 17, 2009                                                     4
Image reconstruction using
Adaptive Normalized Averaging


   Input Image   σ=1     Normalized    Output Image
   (sparsely &           Averaging     (with local
   randomly                            structure
   sampled)            Adaptive        extended into
                       applicability   missing regions)

                          Structure
                          Analysis


July 17, 2009                                             5
Local structure adaptive filtering
•Local structure from the structure tensor                                          r
                                                                                    v
                                                                            r
                      r r T                                                 u ϕ
                                  rr T  rr T
                T = ∇I ∇I = λu uu + λv vv
                                      r
                •orientation φ = arg(u)
                                                                       y = 1 κ x2
                •anisotropy A = (λu - λv)/(λu + λv)                        2
                                        r
                •curvature    κ = ∂φ /∂ v
                •scale    rdensity = sample density

•Scale-adaptive curvature-bent anisotropic Gaussian
kernel with scales in 2 orthogonal directions:
     σ u = C (1 − A)α rdensity σ v = C (1 + A)α rdensity
where           C ~ SNR       α   ~ degree of structure enhancement   kernel aligns with
                                                                       local structure
July 17, 2009                                                                           6
Sample Density Transform
     •Definition: Smallest radius of a pillbox, centered at each pixel, that
     encompasses total certainty of at least 1

     •Role: Automatic scale selection of the applicability in the NA equation
     to avoid unnecessary smoothing




Lena with missing hole   Density transform   NA with Gaussian(σ=1)   Adap. Norm. Avg.

     July 17, 2009                                                              7
4x4 super-resolution from 4 noisy frames
•     4 input LowRes captured with fill-factor = 25%, intensity noise
      (σ=10), registration noise (σ=0.2 LR pitch)




1 of 4 input 64x64 LR         SR using triangulation   SR using adaptive NA
• 16 times upsampling from only 4 frames. How is it possible: along
  linear structures, only 4 samples are enough for 4x super-resolution
 July 17, 2009                                                          8
Orientation            Anisotropy                     Curvature




Sample density   Scale along linear structures Scale in perpendicular direction
July 17, 2009                                                         9
Comparison with image inpainting

• Image inpainting (Sapiro) = diffusion with level line evolution
        • also extending orientation into the missing regions
        • slow due to iterative nature
        • poor result for large holes




     input                inpainting          inpainting +      Adapt. Norm. Avg.
                       110 iters (6 min)    texture synthesis     0 iters (6 sec)
July 17, 2009                                                                10
Directions for Further Research

• Applications
        • Image filtering (noise/watermark removal, edge enhancement...)
        • Image interpolation from sparsely and randomly sampled data
          (image inpainting, image fusion, super-resolution...)


• Further improvements
        • Scale-space local structure analysis.
        • Detect multiple orientations using orientation space.
        • Robust neighborhood operator than the weighted mean.




July 17, 2009                                                      11
Image inpainting of thin scribbles




       input




  inpainting     Adaptive Normalized Averaging (10 sec)
July 17, 2009                                         12
Simultaneous geometry/texture inpainting




                                        texture
  input




geometry                               Adaptive NA
                                         (1 min)

    July 17, 2009                            13
Inpainting of ambiguous discontinuity




      original   input   inpainting   Adaptive NA (1 sec)




July 17, 2009                                        14

More Related Content

PPTX
14.20 o1 m miller
PDF
57969246 finite-volume
PPTX
3.13 usefulness and limitation of the cie system
PPTX
3.8 the 1931 cie s ystem
PDF
Biao Hou--SAR IMAGE DESPECKLING BASED ON IMPROVED DIRECTIONLET DOMAIN GAUSSIA...
PPTX
3.9 addition to the cie system
PDF
High-Precision Optical Interferometry of Faint and Low-Contrast Targets Requi...
PPTX
Standard illuminants
14.20 o1 m miller
57969246 finite-volume
3.13 usefulness and limitation of the cie system
3.8 the 1931 cie s ystem
Biao Hou--SAR IMAGE DESPECKLING BASED ON IMPROVED DIRECTIONLET DOMAIN GAUSSIA...
3.9 addition to the cie system
High-Precision Optical Interferometry of Faint and Low-Contrast Targets Requi...
Standard illuminants

What's hot (6)

PDF
Robust Shape and Topology Optimization - Northwestern
PPTX
study Image and video abstraction by multi scale anisotropic kuwahara
PDF
Learning Moving Cast Shadows for Foreground Detection (VS 2008)
PPT
Pisek_20110726_IGARSS2011_v2.ppt
PDF
Topological Inference via Meshing
PDF
Ostst 2014 - scagliola
Robust Shape and Topology Optimization - Northwestern
study Image and video abstraction by multi scale anisotropic kuwahara
Learning Moving Cast Shadows for Foreground Detection (VS 2008)
Pisek_20110726_IGARSS2011_v2.ppt
Topological Inference via Meshing
Ostst 2014 - scagliola
Ad

Viewers also liked (15)

PDF
Super resolution microscopy & astronomy
PPS
Icdecs 2011
PDF
Assignment2 analog to digital conversion soumit_mukherjee
PDF
Kalman filter - Applications in Image processing
PPT
Retina preliminary
PPTX
Single Image Super-Resolution from Transformed Self-Exemplars (CVPR 2015)
PPTX
Data compression
PDF
Image processing fundamentals
PPT
The Low Vision Examination
PPTX
Matlab Image Restoration Techniques
PPT
Sensors On 3d Digitization
PPT
Image processing
PPT
Bio-metrics Technology
PPTX
Data Compression (Lossy and Lossless)
PPTX
Digital image processing img smoothning
Super resolution microscopy & astronomy
Icdecs 2011
Assignment2 analog to digital conversion soumit_mukherjee
Kalman filter - Applications in Image processing
Retina preliminary
Single Image Super-Resolution from Transformed Self-Exemplars (CVPR 2015)
Data compression
Image processing fundamentals
The Low Vision Examination
Matlab Image Restoration Techniques
Sensors On 3d Digitization
Image processing
Bio-metrics Technology
Data Compression (Lossy and Lossless)
Digital image processing img smoothning
Ad

Similar to Normalized averaging using adaptive applicability functions with applications in image reconstruction from sparsely and randomly sampled data (20)

PDF
Normalized averaging using adaptive applicability functions with applications...
PDF
Image Smoothing for Structure Extraction
PDF
Modified adaptive bilateral filter for image contrast enhancement
PDF
Andrew Zisserman Talk - Part 1a
PPT
Fingerprint High Level Classification
PPT
Edges and lines
PPT
SIFT.ppt
PPTX
study Domain Transform for Edge-Aware Image and Video Processing
PPT
M.sc. m hassan
PDF
Semantic Filtering (An Image Processing Method)
PDF
Image Restoration for 3D Computer Vision
PDF
SinogramReconstruction
PDF
Miller - Remote Sensing and Imaging Physics - Spring Review 2012
PDF
Image processing
PPTX
chapter 4 computervision PART1.pcomputerptx
PPTX
chapter 4 computervision.PPT.pptx ABOUT COMPUTER VISION
PDF
Pixel Matching from Stereo Images (Callan seminar)
PDF
Lecture 02 internet video search
PDF
Robust Super-Resolution by minimizing a Gaussian-weighted L2 error norm
Normalized averaging using adaptive applicability functions with applications...
Image Smoothing for Structure Extraction
Modified adaptive bilateral filter for image contrast enhancement
Andrew Zisserman Talk - Part 1a
Fingerprint High Level Classification
Edges and lines
SIFT.ppt
study Domain Transform for Edge-Aware Image and Video Processing
M.sc. m hassan
Semantic Filtering (An Image Processing Method)
Image Restoration for 3D Computer Vision
SinogramReconstruction
Miller - Remote Sensing and Imaging Physics - Spring Review 2012
Image processing
chapter 4 computervision PART1.pcomputerptx
chapter 4 computervision.PPT.pptx ABOUT COMPUTER VISION
Pixel Matching from Stereo Images (Callan seminar)
Lecture 02 internet video search
Robust Super-Resolution by minimizing a Gaussian-weighted L2 error norm

More from Tuan Q. Pham (14)

PDF
Oral presentation on Asymmetric recursive Gaussian filtering for space-varia...
PDF
Asymmetric recursive Gaussian filtering for space-variant artificial bokeh
PPTX
Parallel implementation of geodesic distance transform with application in su...
PDF
Parallel implementation of geodesic distance transform with application in su...
PDF
Multi-hypothesis projection-based shift estimation for sweeping panorama reco...
PDF
Multi-hypothesis projection-based shift estimation for sweeping panorama reco...
PPT
Non-maximum suppression using fewer than two comparison per pixels
PDF
Paper fingerprinting using alpha-masked image matching
PDF
Paper fingerprinting using alpha-masked image matching
PDF
Bidirectional bias correction for gradient-based shift estimation
PPT
Resolution enhancement of low-quality videos using a high-resolution frame
PPT
Separable bilateral filtering for fast video preprocessing
PDF
Performance of Optimal Registration Estimator
PDF
Influence of Signal-to-Noise Ratio and Point Spread Function on Limits of Sup...
Oral presentation on Asymmetric recursive Gaussian filtering for space-varia...
Asymmetric recursive Gaussian filtering for space-variant artificial bokeh
Parallel implementation of geodesic distance transform with application in su...
Parallel implementation of geodesic distance transform with application in su...
Multi-hypothesis projection-based shift estimation for sweeping panorama reco...
Multi-hypothesis projection-based shift estimation for sweeping panorama reco...
Non-maximum suppression using fewer than two comparison per pixels
Paper fingerprinting using alpha-masked image matching
Paper fingerprinting using alpha-masked image matching
Bidirectional bias correction for gradient-based shift estimation
Resolution enhancement of low-quality videos using a high-resolution frame
Separable bilateral filtering for fast video preprocessing
Performance of Optimal Registration Estimator
Influence of Signal-to-Noise Ratio and Point Spread Function on Limits of Sup...

Recently uploaded (20)

PDF
NewMind AI Weekly Chronicles - August'25-Week II
PDF
Profit Center Accounting in SAP S/4HANA, S4F28 Col11
PDF
Agricultural_Statistics_at_a_Glance_2022_0.pdf
DOCX
The AUB Centre for AI in Media Proposal.docx
PDF
Building Integrated photovoltaic BIPV_UPV.pdf
PDF
Blue Purple Modern Animated Computer Science Presentation.pdf.pdf
PPTX
KOM of Painting work and Equipment Insulation REV00 update 25-dec.pptx
PDF
TokAI - TikTok AI Agent : The First AI Application That Analyzes 10,000+ Vira...
PDF
Empathic Computing: Creating Shared Understanding
PPTX
20250228 LYD VKU AI Blended-Learning.pptx
PDF
Peak of Data & AI Encore- AI for Metadata and Smarter Workflows
PDF
Mobile App Security Testing_ A Comprehensive Guide.pdf
PPT
Teaching material agriculture food technology
PDF
Review of recent advances in non-invasive hemoglobin estimation
PPTX
ACSFv1EN-58255 AWS Academy Cloud Security Foundations.pptx
PDF
cuic standard and advanced reporting.pdf
PDF
Reach Out and Touch Someone: Haptics and Empathic Computing
PDF
A comparative analysis of optical character recognition models for extracting...
PDF
Diabetes mellitus diagnosis method based random forest with bat algorithm
PPTX
Machine Learning_overview_presentation.pptx
NewMind AI Weekly Chronicles - August'25-Week II
Profit Center Accounting in SAP S/4HANA, S4F28 Col11
Agricultural_Statistics_at_a_Glance_2022_0.pdf
The AUB Centre for AI in Media Proposal.docx
Building Integrated photovoltaic BIPV_UPV.pdf
Blue Purple Modern Animated Computer Science Presentation.pdf.pdf
KOM of Painting work and Equipment Insulation REV00 update 25-dec.pptx
TokAI - TikTok AI Agent : The First AI Application That Analyzes 10,000+ Vira...
Empathic Computing: Creating Shared Understanding
20250228 LYD VKU AI Blended-Learning.pptx
Peak of Data & AI Encore- AI for Metadata and Smarter Workflows
Mobile App Security Testing_ A Comprehensive Guide.pdf
Teaching material agriculture food technology
Review of recent advances in non-invasive hemoglobin estimation
ACSFv1EN-58255 AWS Academy Cloud Security Foundations.pptx
cuic standard and advanced reporting.pdf
Reach Out and Touch Someone: Haptics and Empathic Computing
A comparative analysis of optical character recognition models for extracting...
Diabetes mellitus diagnosis method based random forest with bat algorithm
Machine Learning_overview_presentation.pptx

Normalized averaging using adaptive applicability functions with applications in image reconstruction from sparsely and randomly sampled data

  • 1. Normalized averaging using adaptive applicability functions with applications in image reconstruction from sparsely and randomly sampled data Presented at SCIA 2003 Tuan Q. Pham and Lucas J. van Vliet July 17, 2009 1 Pattern Recognition Group
  • 2. Overview • Normalized averaging • Local structure adaptive filtering • Experimental results • Comparison with diffusion-based image inpainting • Directions for further research July 17, 2009 2
  • 3. Normalized averaging (Knutsson’93) •Weighted average filtering: r = s*a •Normalized averaging = weighted average + signal/certainty principle: •each signal s is associated with a certainty c •s & c have to be processed separately ( s . c) ∗ a r= c∗a where s :signal, c :certainty, a :filter, r :result, * :convolution input with 10% Gaussian smoothing NA with Gaussian original pixels (σ = 1) applicability (σ = 1) July 17, 2009 3
  • 4. Normalized averaging: An example Reconstruction from 10% random pixels Nearest neighbor interpolation NA with adaptive applicability July 17, 2009 4
  • 5. Image reconstruction using Adaptive Normalized Averaging Input Image σ=1 Normalized Output Image (sparsely & Averaging (with local randomly structure sampled) Adaptive extended into applicability missing regions) Structure Analysis July 17, 2009 5
  • 6. Local structure adaptive filtering •Local structure from the structure tensor r v r r r T u ϕ rr T rr T T = ∇I ∇I = λu uu + λv vv r •orientation φ = arg(u) y = 1 κ x2 •anisotropy A = (λu - λv)/(λu + λv) 2 r •curvature κ = ∂φ /∂ v •scale rdensity = sample density •Scale-adaptive curvature-bent anisotropic Gaussian kernel with scales in 2 orthogonal directions: σ u = C (1 − A)α rdensity σ v = C (1 + A)α rdensity where C ~ SNR α ~ degree of structure enhancement kernel aligns with local structure July 17, 2009 6
  • 7. Sample Density Transform •Definition: Smallest radius of a pillbox, centered at each pixel, that encompasses total certainty of at least 1 •Role: Automatic scale selection of the applicability in the NA equation to avoid unnecessary smoothing Lena with missing hole Density transform NA with Gaussian(σ=1) Adap. Norm. Avg. July 17, 2009 7
  • 8. 4x4 super-resolution from 4 noisy frames • 4 input LowRes captured with fill-factor = 25%, intensity noise (σ=10), registration noise (σ=0.2 LR pitch) 1 of 4 input 64x64 LR SR using triangulation SR using adaptive NA • 16 times upsampling from only 4 frames. How is it possible: along linear structures, only 4 samples are enough for 4x super-resolution July 17, 2009 8
  • 9. Orientation Anisotropy Curvature Sample density Scale along linear structures Scale in perpendicular direction July 17, 2009 9
  • 10. Comparison with image inpainting • Image inpainting (Sapiro) = diffusion with level line evolution • also extending orientation into the missing regions • slow due to iterative nature • poor result for large holes input inpainting inpainting + Adapt. Norm. Avg. 110 iters (6 min) texture synthesis 0 iters (6 sec) July 17, 2009 10
  • 11. Directions for Further Research • Applications • Image filtering (noise/watermark removal, edge enhancement...) • Image interpolation from sparsely and randomly sampled data (image inpainting, image fusion, super-resolution...) • Further improvements • Scale-space local structure analysis. • Detect multiple orientations using orientation space. • Robust neighborhood operator than the weighted mean. July 17, 2009 11
  • 12. Image inpainting of thin scribbles input inpainting Adaptive Normalized Averaging (10 sec) July 17, 2009 12
  • 13. Simultaneous geometry/texture inpainting texture input geometry Adaptive NA (1 min) July 17, 2009 13
  • 14. Inpainting of ambiguous discontinuity original input inpainting Adaptive NA (1 sec) July 17, 2009 14