SlideShare a Scribd company logo
Conf. 5672: Image Processing
         Algorithms and Systems IV




Influence of Signal-to-Noise Ratio
  and Point Spread Function on
    Limits of Super-Resolution

               Tuan Pham

         Quantitative Imaging Group
        Delft University of Technology
               The Netherlands
Super-Resolution: an example




128x128x100 infra-red sequence   4x super-resolution
  © 2004 Tuan Pham                                     2
Super-Resolution: an example




         Low resolution   4x super-resolution
© 2004 Tuan Pham                                3
Overview and Goal
                                   System inputs
                   No. of inputs       SNR         PSF




          Positioning                  SNR           Resolving
             limit                     limit           limit


  GOAL: Derive the Limits of Super-Resolution given system inputs
© 2004 Tuan Pham                                                    4
Limit of registration
• Cramer-Rao Lower Bound for 2D shift: I2(x, y) = I1(x+vx, y+vy) :
                      var(v x ) ≥ F111 = σ n ∑ I2 Det (F)
                                    −      2
                                                y
                                               S

                      var(v y ) ≥ F221 = σ n ∑ I2 Det (F)
                                   −       2
                                                x
                                               S

where I x = ∂I / ∂x , I y = ∂I / ∂y , σ n is noise variance, and F is the
                                        2


 Fisher Information Matrix:
                                     ⎡ ∑ I2        ∑I I ⎤
                                   1 ⎢ S x              x y
                                                        ⎥
                          F( v ) = 2 ⎢             S

                                  σ n ∑ IxIy
                                     ⎢S
                                     ⎣
                                                   ∑I ⎥
                                                    S
                                                        ⎥
                                                        ⎦
                                                         2
                                                         y




• Optimal registration is achievable by iterative optimization

• CRLB also exists for more complicated motion models:
        - 2D projective                                - optic flow

© 2004 Tuan Pham                                                            5
Noise of HR image after fusion

 • Total noise = Intensity noise + Noise due to registration error
                                                                             μ : zoom factor
                                  μ2                      2                  N : # of LR images
               σ   2
                   n     =             σ   I
                                            2
                                                +   ∇I σ      2
                                                              reg
                                                                                2
                                                                             ∇I : gradient energy
                                  N
                                                               position error distribution




                                                                                                Intensity error distribution
                                                                              σreg
                                                                         I
                                                                                         σI
                                                                                     x
                                                               local
                                                               signal
                                                                                ∂I
                                                                        σI =       σ
                                                                                ∂x
                                                                                         re g

  Blurred & mis-registered              Noise due to
5x5 box blur, σ reg = 0.2 pixel        mis-registration         mis-registration → noise
  © 2004 Tuan Pham                                                                              6
The need for deconvolution

• After fusion, the High-Resolution image is still blurry due to:
   – Sensor integration blur (severe if high fill-factor)
   – Optical blur (severe if high sampling factor)




                                                       On-chip microlens of
                                                      Sony Super HAD CCD
© 2004 Tuan Pham                                                              7
The necessity of aliasing

                                         • Spectrum is cut off beyond fc due to optics → data forever lost

                                           1                                                                                                    1
                                                                      OTF (sampling factor = 0.25)                                                                              OTF (sampling factor = 1)
frequency spectra / transfer functions




                                                                                                     frequency spectra / transfer functions
                                                                      STF (fill factor = 1)                                                                                     STF (fill factor = 1)
                                          0.8                         Original scene spectrum                                                  0.8                              Original scene spectrum
                                                                      Band−limited spectrum                                                                                     Band−limited spectrum
                                                                      Aliased image spectrum                                                                                    Sampled image spectrum
                                          0.6                                                                                                  0.6


                                          0.4                                                                                                  0.4


                                          0.2                                                                                                  0.2


                                           0                                                                                                    0


                                         −0.2                                                                                                 −0.2


                                         −0.4                                                                                                 −0.4
                                             0         0.5             1             1.5                                                          0         0.5             1             1.5               2
                                                 frequency in unit of sampling frequency (f/fs)                                                       frequency in unit of sampling frequency (f/fs)


                                                     Aliasing due to                                                                                         No aliasing at
                                                 under-sampling (fs < 2fc)                                                                            critical sampling (fs = 2fc)

                                © 2004 Tuan Pham                                                                                                                                                       8
Limit of deconvolution
• Blur = attenuation of HF spectrum                                      recoverable

• Deconvolution = amplify HF spectrum:
   – noise is also amplified → limit the deconvolution             PS>PN

• Deconvolution can only recover:
                                                               Not
   – Spectrum whose signal power > noise power             recoverable
                                                           resolution factor = 0.44




        fusion result
© 2004 Tuan Pham               after deconvolution       simulated at resolution = 0.44
                                                                                    9
SR reconstruction experiment
  • Aim: show that the attainable SR factor agrees with the prediction
  • Experiment:
     – Inputs: sufficient shifted LR images of the Pentagon
     – Output: SR image and a measure of SR factor from edge width




    64x64 LR input               4xHR after fusion        4xSR after deconvolution
sampling=1/4, fill = 100%         BSNR = 20 dB                SR factor = 3.4
  © 2004 Tuan Pham                                                           10
SR reconstruction experiment
  • Aim: show that the attainable SR factor agrees with the prediction
  • Experiment:
     – Inputs: sufficient shifted LR images of the Pentagon
     – Output: SR image and a measure of SR factor from edge width




    64x64 LR input               4xHR after fusion        4xSR after deconvolution
sampling=1/4, fill = 100%         BSNR = 20 dB                SR factor = 3.4
  © 2004 Tuan Pham                                                           11
SR factor at BSNR=20dB
            • Consistent results between prediction and measurement:
                – Attainable SR factor depends mainly on sampling factor (i.e. level of aliasing)




            6                                         3.2
                                                                           6                                             3.0

                                                     3.4                                                             2.5
                                                                           4




                                                                SR limit
SR factor




            4
                                                                                                                 1.7
                                               1.9                         2
            2                                                                                              1.0
                                         1.0
                                                                           0                        0.6
            0                                                              0                                                   0
            0                                               0
                                0.6
                                                                                      0.5                     1
                       0.5                 1                                                                sampling factor
                                       sampling factor                         fill factor          1 2       (f /2f )
                fill factor    1 2       (fs/2fc)                                                                s   c



                     Measured SR factor                                                      Predicted SR factor

            © 2004 Tuan Pham                                                                                              12
Summary


 • Limit of Super-Resolution depends on:
    – input Signal-to-Noise Ratio

    – System’s Point Spread Function and how well it can be estimated



 • Procedure for estimating SR factor directly from inputs:
    – Measure noise variance from LR images   σ I2
    – Derive registration error   σ reg
                                    2


    – Determine SR factor from the Power Spectrum Density (PS > PN)




© 2004 Tuan Pham                                                        13

More Related Content

PPTX
Facial Expression Recognition / Removal
PPT
Nanotechnology
PDF
Cd Simon
PDF
CVPR2010: Advanced ITinCVPR in a Nutshell: part 4: Isocontours, Registration
PDF
Mesh Processing Course : Differential Calculus
PDF
Elementary Landscape Decomposition of Combinatorial Optimization Problems
PDF
Identity Based Encryption
PDF
Modern features-part-2-descriptors
Facial Expression Recognition / Removal
Nanotechnology
Cd Simon
CVPR2010: Advanced ITinCVPR in a Nutshell: part 4: Isocontours, Registration
Mesh Processing Course : Differential Calculus
Elementary Landscape Decomposition of Combinatorial Optimization Problems
Identity Based Encryption
Modern features-part-2-descriptors

What's hot (20)

PDF
Color Img at Prisma Network meeting 2009
PDF
Modern features-part-1-detectors
PPTX
Kccsi 2012 a real-time robust object tracking-v2
PDF
Exponentials integrals
PDF
23 industrial engineering
PDF
Learning Spectral Graph Transformations for Link Prediction
PDF
Tro07 sparse-solutions-talk
PPTX
14.20 o1 m miller
PDF
Biao Hou--SAR IMAGE DESPECKLING BASED ON IMPROVED DIRECTIONLET DOMAIN GAUSSIA...
PDF
Visual Odomtery(2)
PPT
Chapter14
PPSX
Indices and standard form
PDF
Csr2011 june14 17_00_pospelov
PDF
Bouguet's MatLab Camera Calibration Toolbox
PDF
A Physical Approach to Moving Cast Shadow Detection (ICASSP 2009)
PDF
Modern features-part-3-software
PDF
Fatigue damage in solder joint interconnects - presentation
PDF
Form 5 formulae and note
PDF
Learning Moving Cast Shadows for Foreground Detection (VS 2008)
PPT
17.04.2012 seminar trions_kochereshko
Color Img at Prisma Network meeting 2009
Modern features-part-1-detectors
Kccsi 2012 a real-time robust object tracking-v2
Exponentials integrals
23 industrial engineering
Learning Spectral Graph Transformations for Link Prediction
Tro07 sparse-solutions-talk
14.20 o1 m miller
Biao Hou--SAR IMAGE DESPECKLING BASED ON IMPROVED DIRECTIONLET DOMAIN GAUSSIA...
Visual Odomtery(2)
Chapter14
Indices and standard form
Csr2011 june14 17_00_pospelov
Bouguet's MatLab Camera Calibration Toolbox
A Physical Approach to Moving Cast Shadow Detection (ICASSP 2009)
Modern features-part-3-software
Fatigue damage in solder joint interconnects - presentation
Form 5 formulae and note
Learning Moving Cast Shadows for Foreground Detection (VS 2008)
17.04.2012 seminar trions_kochereshko
Ad

Similar to Influence of Signal-to-Noise Ratio and Point Spread Function on Limits of Super-Resolution (20)

PDF
Mapping the discrete logarithm
PPT
fghdfh
PDF
Image formation
PDF
RSA Cryptosystem
PDF
RSA Cryptosystem
PDF
RSA Cryptosystem
PDF
Topological Inference via Meshing
PDF
Second Order Perturbations During Inflation Beyond Slow-roll
PDF
IGARSS2011 FR3.T08.3 BenDavid.pdf
PDF
A Decomposition-based Approach to Modeling and Understanding Arbitrary Shapes
PDF
Lecture 02 internet video search
PDF
Reconstruction (of micro-objects) based on focus-sets using blind deconvoluti...
PDF
CVPR2010: Advanced ITinCVPR in a Nutshell: part 4: additional slides
PDF
Bubbles of True Vacuum and Duality in M-theory
PDF
Matrix Models of 2D String Theory in Non-trivial Backgrounds
PDF
Dsp U Lec06 The Z Transform And Its Application
PDF
Signal Processing Course : Compressed Sensing
PDF
Mapping Ash Tree Colonization in an Agricultural Moutain Landscape_ Investiga...
PPT
Crystallite size nanomaterials
PDF
Solar Cells Lecture 3: Modeling and Simulation of Photovoltaic Devices and Sy...
Mapping the discrete logarithm
fghdfh
Image formation
RSA Cryptosystem
RSA Cryptosystem
RSA Cryptosystem
Topological Inference via Meshing
Second Order Perturbations During Inflation Beyond Slow-roll
IGARSS2011 FR3.T08.3 BenDavid.pdf
A Decomposition-based Approach to Modeling and Understanding Arbitrary Shapes
Lecture 02 internet video search
Reconstruction (of micro-objects) based on focus-sets using blind deconvoluti...
CVPR2010: Advanced ITinCVPR in a Nutshell: part 4: additional slides
Bubbles of True Vacuum and Duality in M-theory
Matrix Models of 2D String Theory in Non-trivial Backgrounds
Dsp U Lec06 The Z Transform And Its Application
Signal Processing Course : Compressed Sensing
Mapping Ash Tree Colonization in an Agricultural Moutain Landscape_ Investiga...
Crystallite size nanomaterials
Solar Cells Lecture 3: Modeling and Simulation of Photovoltaic Devices and Sy...
Ad

More from Tuan Q. Pham (15)

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
PDF
Robust Super-Resolution by minimizing a Gaussian-weighted L2 error norm
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
PPT
Normalized averaging using adaptive applicability functions with applications...
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
Robust Super-Resolution by minimizing a Gaussian-weighted L2 error norm
Resolution enhancement of low-quality videos using a high-resolution frame
Separable bilateral filtering for fast video preprocessing
Performance of Optimal Registration Estimator
Normalized averaging using adaptive applicability functions with applications...

Recently uploaded (20)

PDF
Profit Center Accounting in SAP S/4HANA, S4F28 Col11
PDF
Build a system with the filesystem maintained by OSTree @ COSCUP 2025
PPTX
KOM of Painting work and Equipment Insulation REV00 update 25-dec.pptx
PDF
cuic standard and advanced reporting.pdf
PPTX
Detection-First SIEM: Rule Types, Dashboards, and Threat-Informed Strategy
PPTX
ACSFv1EN-58255 AWS Academy Cloud Security Foundations.pptx
PDF
How UI/UX Design Impacts User Retention in Mobile Apps.pdf
PDF
Advanced methodologies resolving dimensionality complications for autism neur...
PDF
Encapsulation theory and applications.pdf
PDF
Chapter 3 Spatial Domain Image Processing.pdf
PDF
Blue Purple Modern Animated Computer Science Presentation.pdf.pdf
PDF
Machine learning based COVID-19 study performance prediction
PDF
Reach Out and Touch Someone: Haptics and Empathic Computing
PDF
Electronic commerce courselecture one. Pdf
PPT
“AI and Expert System Decision Support & Business Intelligence Systems”
PDF
MIND Revenue Release Quarter 2 2025 Press Release
PPTX
Big Data Technologies - Introduction.pptx
PDF
Peak of Data & AI Encore- AI for Metadata and Smarter Workflows
PPTX
Understanding_Digital_Forensics_Presentation.pptx
PDF
Review of recent advances in non-invasive hemoglobin estimation
Profit Center Accounting in SAP S/4HANA, S4F28 Col11
Build a system with the filesystem maintained by OSTree @ COSCUP 2025
KOM of Painting work and Equipment Insulation REV00 update 25-dec.pptx
cuic standard and advanced reporting.pdf
Detection-First SIEM: Rule Types, Dashboards, and Threat-Informed Strategy
ACSFv1EN-58255 AWS Academy Cloud Security Foundations.pptx
How UI/UX Design Impacts User Retention in Mobile Apps.pdf
Advanced methodologies resolving dimensionality complications for autism neur...
Encapsulation theory and applications.pdf
Chapter 3 Spatial Domain Image Processing.pdf
Blue Purple Modern Animated Computer Science Presentation.pdf.pdf
Machine learning based COVID-19 study performance prediction
Reach Out and Touch Someone: Haptics and Empathic Computing
Electronic commerce courselecture one. Pdf
“AI and Expert System Decision Support & Business Intelligence Systems”
MIND Revenue Release Quarter 2 2025 Press Release
Big Data Technologies - Introduction.pptx
Peak of Data & AI Encore- AI for Metadata and Smarter Workflows
Understanding_Digital_Forensics_Presentation.pptx
Review of recent advances in non-invasive hemoglobin estimation

Influence of Signal-to-Noise Ratio and Point Spread Function on Limits of Super-Resolution

  • 1. Conf. 5672: Image Processing Algorithms and Systems IV Influence of Signal-to-Noise Ratio and Point Spread Function on Limits of Super-Resolution Tuan Pham Quantitative Imaging Group Delft University of Technology The Netherlands
  • 2. Super-Resolution: an example 128x128x100 infra-red sequence 4x super-resolution © 2004 Tuan Pham 2
  • 3. Super-Resolution: an example Low resolution 4x super-resolution © 2004 Tuan Pham 3
  • 4. Overview and Goal System inputs No. of inputs SNR PSF Positioning SNR Resolving limit limit limit GOAL: Derive the Limits of Super-Resolution given system inputs © 2004 Tuan Pham 4
  • 5. Limit of registration • Cramer-Rao Lower Bound for 2D shift: I2(x, y) = I1(x+vx, y+vy) : var(v x ) ≥ F111 = σ n ∑ I2 Det (F) − 2 y S var(v y ) ≥ F221 = σ n ∑ I2 Det (F) − 2 x S where I x = ∂I / ∂x , I y = ∂I / ∂y , σ n is noise variance, and F is the 2 Fisher Information Matrix: ⎡ ∑ I2 ∑I I ⎤ 1 ⎢ S x x y ⎥ F( v ) = 2 ⎢ S σ n ∑ IxIy ⎢S ⎣ ∑I ⎥ S ⎥ ⎦ 2 y • Optimal registration is achievable by iterative optimization • CRLB also exists for more complicated motion models: - 2D projective - optic flow © 2004 Tuan Pham 5
  • 6. Noise of HR image after fusion • Total noise = Intensity noise + Noise due to registration error μ : zoom factor μ2 2 N : # of LR images σ 2 n = σ I 2 + ∇I σ 2 reg 2 ∇I : gradient energy N position error distribution Intensity error distribution σreg I σI x local signal ∂I σI = σ ∂x re g Blurred & mis-registered Noise due to 5x5 box blur, σ reg = 0.2 pixel mis-registration mis-registration → noise © 2004 Tuan Pham 6
  • 7. The need for deconvolution • After fusion, the High-Resolution image is still blurry due to: – Sensor integration blur (severe if high fill-factor) – Optical blur (severe if high sampling factor) On-chip microlens of Sony Super HAD CCD © 2004 Tuan Pham 7
  • 8. The necessity of aliasing • Spectrum is cut off beyond fc due to optics → data forever lost 1 1 OTF (sampling factor = 0.25) OTF (sampling factor = 1) frequency spectra / transfer functions frequency spectra / transfer functions STF (fill factor = 1) STF (fill factor = 1) 0.8 Original scene spectrum 0.8 Original scene spectrum Band−limited spectrum Band−limited spectrum Aliased image spectrum Sampled image spectrum 0.6 0.6 0.4 0.4 0.2 0.2 0 0 −0.2 −0.2 −0.4 −0.4 0 0.5 1 1.5 0 0.5 1 1.5 2 frequency in unit of sampling frequency (f/fs) frequency in unit of sampling frequency (f/fs) Aliasing due to No aliasing at under-sampling (fs < 2fc) critical sampling (fs = 2fc) © 2004 Tuan Pham 8
  • 9. Limit of deconvolution • Blur = attenuation of HF spectrum recoverable • Deconvolution = amplify HF spectrum: – noise is also amplified → limit the deconvolution PS>PN • Deconvolution can only recover: Not – Spectrum whose signal power > noise power recoverable resolution factor = 0.44 fusion result © 2004 Tuan Pham after deconvolution simulated at resolution = 0.44 9
  • 10. SR reconstruction experiment • Aim: show that the attainable SR factor agrees with the prediction • Experiment: – Inputs: sufficient shifted LR images of the Pentagon – Output: SR image and a measure of SR factor from edge width 64x64 LR input 4xHR after fusion 4xSR after deconvolution sampling=1/4, fill = 100% BSNR = 20 dB SR factor = 3.4 © 2004 Tuan Pham 10
  • 11. SR reconstruction experiment • Aim: show that the attainable SR factor agrees with the prediction • Experiment: – Inputs: sufficient shifted LR images of the Pentagon – Output: SR image and a measure of SR factor from edge width 64x64 LR input 4xHR after fusion 4xSR after deconvolution sampling=1/4, fill = 100% BSNR = 20 dB SR factor = 3.4 © 2004 Tuan Pham 11
  • 12. SR factor at BSNR=20dB • Consistent results between prediction and measurement: – Attainable SR factor depends mainly on sampling factor (i.e. level of aliasing) 6 3.2 6 3.0 3.4 2.5 4 SR limit SR factor 4 1.7 1.9 2 2 1.0 1.0 0 0.6 0 0 0 0 0 0.6 0.5 1 0.5 1 sampling factor sampling factor fill factor 1 2 (f /2f ) fill factor 1 2 (fs/2fc) s c Measured SR factor Predicted SR factor © 2004 Tuan Pham 12
  • 13. Summary • Limit of Super-Resolution depends on: – input Signal-to-Noise Ratio – System’s Point Spread Function and how well it can be estimated • Procedure for estimating SR factor directly from inputs: – Measure noise variance from LR images σ I2 – Derive registration error σ reg 2 – Determine SR factor from the Power Spectrum Density (PS > PN) © 2004 Tuan Pham 13