SlideShare a Scribd company logo
ICME 2005
IEEE International Conference on Multimedia & Expo
July 6-8, 2005, Amsterdam, The Netherlands




       Separable bilateral filtering
      for fast video preprocessing

                       Tuan Q. Pham      &    Lucas J. van Vliet



                      Delft University of Technology
                             The Netherlands                       1



        Quantitative Imaging Group
        Faculty of Applied Sciences
Contents

1. Bilateral filtering:
    Edge-preserving filter
    High computational complexity


2. Separable implementation:
    Good approximation of the original filter
    Linear computational complexity


3. Application to video preprocessing:
    Noise reduction
    Better compressed video
Gaussian filtering revisited

 Gaussian filtering: weights depend on distance to the center pixel



                         *                             =
     Noisy step edge              Gaussian weights         Gaussian filtered result


 Adaptive filtering: avoid filtering across edges



                         .                             =
     Noisy step edge         Edge-preserving weights   Edge-presered filtered result
Gaussian vs bilateral filtering

 Gaussian filtering: weights depend on distance to the center pixel



                        *                              =
     Noisy step edge              Gaussian weights         Gaussian filtered result


 Bilateral filtering: weights depend on both spatial closeness and
  photometric similarity


                         .                             =
     Noisy step edge         Edge-preserving weights   Edge-presered filtered result
How bilateral filtering works?




 Every sample is a weighted average of its neighbors:
                             1
                   O ( s0 ) = ∑ w ( s, s0 ) . I ( s )
                             K s
                                                                          Local modes
 The weight w = w s . wt is product of two Gaussian                 weights: after
                                                                            bilateral
                                   − ( s − s0 ) 2                          filtering
    Spatial proximity: w s = exp                 
                                       2σ s2 
                                     − ( I ( s ) − I ( s0 ) ) 2 
    Tonal similarity:     wt = exp                             
                                              2σ t 2
                                                                 
                                                                
Example: Gaussian filtering




 Noisy input: PSNR = 39.1 dB   Gaussian filtered: PSNR = 67.9 dB
Example: Bilateral filtering




 Noisy input: PSNR = 39.1 dB   Bilateral filtered: PSNR = 41.6 dB
Computational complexity

 Bilateral filtering kernel is space-variant → complexity is:

                         O ( Nm d )
      where N: number of pixels in the image
            m: size of filtering kernel (m ≈ 7 is good enough)
             d: image dimensionality

 Previous attempt for fast bilateral filtering:
   Piecewise-linear: approximate bilateral filtering with M Gaussian
     filtering - Durand & Dorsey (SIGGRAPH 2002) → complexity is:

                    O ( N log( N ) . M )       M ≈ 17 for 8-bit images

 Our approach: separable bilateral filtering
                          O ( Nmd )
Is bilateral filter separable?

 Gaussian filter is space-invariant and separable:
                                                           g(x)
                                       g(y)
         ∑ w ( s − s ) . I ( s)
                  s       0
O( s ) =  s
                                                  *               =
          ∑ w (s − s )
    0
                      s       0
              s
                                              y
                                                                  x

                                                  w s ( s ) = gauss( x ) . gauss( y )
                                                                                  Kernel center
 Bilateral filter is NOT separable:                                  I(s0)
     Space-variant kernel due to local                                   I(s)
      intensity dependency



 However, even a highly non-linear filter like median filter is
  approximately separable (Narendra – PAMI 1981)
Separable bilateral filtering result

 Separable bilateral filtering is a good approximation of full kernel filtering:




 noisy Erika σnoise = 10   bilateral filtered in x-dimension followed by y-dimension filtering
Separable bilateral filtering result

 Separable bilateral filtering is a good approximation of full kernel filtering:




  noisy Erika σnoise = 10            bilateral filtered in x-dimension followed by y-dimension filtering

Image size             Brute-force           Durrand 2002        Separable        Aniso. diffusion
256x256                      4.46                 0.37              0.21                2.76
512x512                     17.88                 1.59              0.89               26.02
61x61x61                     5.29                 4.20              0.45                5.25
256x256x212                 56 min          Out of memory           50.3          Out of memory
How separable bilateral filtering works?




   Pros: extremely fast (fixed spatial weight + LUT for tonal weight)
   Cons: effective filtering kernel is a slightly distorted
Performance of separable bilateral
filtering




 Very good approximation of    Almost linear execution time
  the full-kernel filter         per pixel
MPEG-1 Foreman with bilateral preprocessing

               16        without preprocessing
                         with 3x3x3 full-kernel
                         with 9x9x5 separable

               15
        RMSE




               14



               13
                0         400              800         1200
                           bit-rate (K bits/s)


   Better RMSE is achieved with separable bilateral filtering compared to
    full-kernel bilateral filtering with the same computation requirement
Less artifact with Bilateral preprocessing
Conclusions


 Separable implementation of bilateral filtering:
    Very good approximation of the original filter
    Much faster than the original or other approximations


 Applications in video preprocessing:
    Improved quality of compressed video
    Reduced processing time → real-time possibility
Literature

   C. Tomasi and R. Manduchi, Bilateral fitering for gray and color
    images, Proc. of ICCV, USA, 1998, 839-846.

   T.Q. Pham and L.J. van Vliet, Separable bilateral filtering for fast
    video preprocessing, Proc. of ICME’05.

   F. Durrand and J. Dorsey, Fast bilateral filtering for the display of
    high dynamic range images, Proc. of SIGGRAPH’02, 2002, 844-847.

   P. Perona and J. Malik, Scale-space filtering and edge detection
    using anisotropic diffusion, PAMI, vol. 12, no. 7, 1990, 629-639.

   R. v.d. Boomgaard and J. v.d. Weijer. On the equivalence of local-
    mode finding, robust estimation and mean-shift analysis as used in
    early vision tasks. In Proc. of ICPR, pages 927-930, 2002.

More Related Content

PDF
Cisco Security Presentation
PPTX
Blueprint Tips 30選
PDF
Need Of Security Operations Over SIEM
PPTX
Unreal Engine (For Creating Games) Presentation
PPTX
Siber Tehdit Gözetleme ve SIEM Olarak Açık Kaynak Sistemlerin Kullanımı
PDF
Get Real-Time Cyber Threat Protection with Risk Management and SIEM
PPSX
Game Design Document
PPT
Slide Tesi informatica
Cisco Security Presentation
Blueprint Tips 30選
Need Of Security Operations Over SIEM
Unreal Engine (For Creating Games) Presentation
Siber Tehdit Gözetleme ve SIEM Olarak Açık Kaynak Sistemlerin Kullanımı
Get Real-Time Cyber Threat Protection with Risk Management and SIEM
Game Design Document
Slide Tesi informatica

What's hot (20)

PDF
【UE4.25 新機能】新しいシリアライゼーション機能「Unversioned Property Serialization」について
PPTX
Python / BlueprintによるUnreal Engineの自動化 / GTMF2019
PDF
Securing a Web App with Security Keys
PPTX
Taking Splunk to the Next Level - Architecture
PDF
Cybersecurity in Industrial Control Systems (ICS)
PPTX
CCSK Certificate of Cloud Computing Knowledge - overview
PPTX
Drone tech hand gesture controlled drones
PPS
God Of War : post mortem
PPT
Design Document
PPT
Extend Your Market Reach with IBM Security QRadar for MSPs
PDF
遅延の少ないLivePreview方法
PDF
第5回ue4ハンズオンセミナー
PPTX
Adopting A Zero-Trust Model. Google Did It, Can You?
PPTX
Nessus-Vulnerability Tester
PPTX
Indie Game Development Intro
PPT
Level Design Challenges & Solutions - Mirror's Edge
PDF
IBM QRadar Security Intelligence Overview
PPTX
Yapım ekibindeki kişiler ve görevleri
PDF
Physically Based Lighting in Unreal Engine 4
PDF
個人製作インディーゲーム”ジラフとアンニカ” のUE4 制作事例紹介 | UNREAL FEST EXTREME 2020 WINTER
【UE4.25 新機能】新しいシリアライゼーション機能「Unversioned Property Serialization」について
Python / BlueprintによるUnreal Engineの自動化 / GTMF2019
Securing a Web App with Security Keys
Taking Splunk to the Next Level - Architecture
Cybersecurity in Industrial Control Systems (ICS)
CCSK Certificate of Cloud Computing Knowledge - overview
Drone tech hand gesture controlled drones
God Of War : post mortem
Design Document
Extend Your Market Reach with IBM Security QRadar for MSPs
遅延の少ないLivePreview方法
第5回ue4ハンズオンセミナー
Adopting A Zero-Trust Model. Google Did It, Can You?
Nessus-Vulnerability Tester
Indie Game Development Intro
Level Design Challenges & Solutions - Mirror's Edge
IBM QRadar Security Intelligence Overview
Yapım ekibindeki kişiler ve görevleri
Physically Based Lighting in Unreal Engine 4
個人製作インディーゲーム”ジラフとアンニカ” のUE4 制作事例紹介 | UNREAL FEST EXTREME 2020 WINTER
Ad

Similar to Separable bilateral filtering for fast video preprocessing (20)

PPTX
Summarized notes on Interference and Diffraction for JEE Main
PDF
Lecture 02 internet video search
PDF
Random Valued Impulse Noise Removal in Colour Images using Adaptive Threshold...
PPT
Spatial Filtering in intro image processingr
PDF
channel_mzhazbay.pdf
PDF
Goddard-DR-2010
PPTX
DONY Simple and Practical Algorithm sft.pptx
PPTX
Pres Simple and Practical Algorithm sft.pptx
PDF
Biao Hou--SAR IMAGE DESPECKLING BASED ON IMPROVED DIRECTIONLET DOMAIN GAUSSIA...
PDF
Biao Hou--SAR IMAGE DESPECKLING BASED ON IMPROVED DIRECTIONLET DOMAIN GAUSSIA...
PDF
EDGE DETECTION IN RADAR IMAGES USING WEIBULL DISTRIBUTION
PDF
EDGE DETECTION IN RADAR IMAGES USING WEIBULL DISTRIBUTION
PPTX
Image processing 2
PDF
Cb25464467
PDF
Xu_IGARSS11_snow_f.pdf
PPT
cos323_s06_lecture13_sigproc bio signal processing
PDF
Robust Super-Resolution by minimizing a Gaussian-weighted L2 error norm
PDF
OE Instrumentation_03_Interferometry_2.pdf
PDF
Presentation of Understanding Sharpness Dynamics in NN Training with a Minima...
PDF
Edge Detection 2 digital image processing
Summarized notes on Interference and Diffraction for JEE Main
Lecture 02 internet video search
Random Valued Impulse Noise Removal in Colour Images using Adaptive Threshold...
Spatial Filtering in intro image processingr
channel_mzhazbay.pdf
Goddard-DR-2010
DONY Simple and Practical Algorithm sft.pptx
Pres Simple and Practical Algorithm sft.pptx
Biao Hou--SAR IMAGE DESPECKLING BASED ON IMPROVED DIRECTIONLET DOMAIN GAUSSIA...
Biao Hou--SAR IMAGE DESPECKLING BASED ON IMPROVED DIRECTIONLET DOMAIN GAUSSIA...
EDGE DETECTION IN RADAR IMAGES USING WEIBULL DISTRIBUTION
EDGE DETECTION IN RADAR IMAGES USING WEIBULL DISTRIBUTION
Image processing 2
Cb25464467
Xu_IGARSS11_snow_f.pdf
cos323_s06_lecture13_sigproc bio signal processing
Robust Super-Resolution by minimizing a Gaussian-weighted L2 error norm
OE Instrumentation_03_Interferometry_2.pdf
Presentation of Understanding Sharpness Dynamics in NN Training with a Minima...
Edge Detection 2 digital image processing
Ad

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
PDF
Performance of Optimal Registration Estimator
PDF
Influence of Signal-to-Noise Ratio and Point Spread Function on Limits of Sup...
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
Resolution enhancement of low-quality videos using a high-resolution frame
Performance of Optimal Registration Estimator
Influence of Signal-to-Noise Ratio and Point Spread Function on Limits of Sup...
Normalized averaging using adaptive applicability functions with applications...

Recently uploaded (20)

PDF
Chapter 3 Spatial Domain Image Processing.pdf
PDF
Reach Out and Touch Someone: Haptics and Empathic Computing
PDF
Profit Center Accounting in SAP S/4HANA, S4F28 Col11
PDF
Optimiser vos workloads AI/ML sur Amazon EC2 et AWS Graviton
PDF
Approach and Philosophy of On baking technology
PDF
Network Security Unit 5.pdf for BCA BBA.
PDF
NewMind AI Weekly Chronicles - August'25 Week I
DOCX
The AUB Centre for AI in Media Proposal.docx
PDF
Electronic commerce courselecture one. Pdf
PDF
Empathic Computing: Creating Shared Understanding
PDF
Advanced methodologies resolving dimensionality complications for autism neur...
PDF
MIND Revenue Release Quarter 2 2025 Press Release
PDF
TokAI - TikTok AI Agent : The First AI Application That Analyzes 10,000+ Vira...
PPTX
ACSFv1EN-58255 AWS Academy Cloud Security Foundations.pptx
PPTX
Effective Security Operations Center (SOC) A Modern, Strategic, and Threat-In...
PPTX
Spectroscopy.pptx food analysis technology
PPT
“AI and Expert System Decision Support & Business Intelligence Systems”
PDF
The Rise and Fall of 3GPP – Time for a Sabbatical?
PDF
Dropbox Q2 2025 Financial Results & Investor Presentation
PPTX
KOM of Painting work and Equipment Insulation REV00 update 25-dec.pptx
Chapter 3 Spatial Domain Image Processing.pdf
Reach Out and Touch Someone: Haptics and Empathic Computing
Profit Center Accounting in SAP S/4HANA, S4F28 Col11
Optimiser vos workloads AI/ML sur Amazon EC2 et AWS Graviton
Approach and Philosophy of On baking technology
Network Security Unit 5.pdf for BCA BBA.
NewMind AI Weekly Chronicles - August'25 Week I
The AUB Centre for AI in Media Proposal.docx
Electronic commerce courselecture one. Pdf
Empathic Computing: Creating Shared Understanding
Advanced methodologies resolving dimensionality complications for autism neur...
MIND Revenue Release Quarter 2 2025 Press Release
TokAI - TikTok AI Agent : The First AI Application That Analyzes 10,000+ Vira...
ACSFv1EN-58255 AWS Academy Cloud Security Foundations.pptx
Effective Security Operations Center (SOC) A Modern, Strategic, and Threat-In...
Spectroscopy.pptx food analysis technology
“AI and Expert System Decision Support & Business Intelligence Systems”
The Rise and Fall of 3GPP – Time for a Sabbatical?
Dropbox Q2 2025 Financial Results & Investor Presentation
KOM of Painting work and Equipment Insulation REV00 update 25-dec.pptx

Separable bilateral filtering for fast video preprocessing

  • 1. ICME 2005 IEEE International Conference on Multimedia & Expo July 6-8, 2005, Amsterdam, The Netherlands Separable bilateral filtering for fast video preprocessing Tuan Q. Pham & Lucas J. van Vliet Delft University of Technology The Netherlands 1 Quantitative Imaging Group Faculty of Applied Sciences
  • 2. Contents 1. Bilateral filtering:  Edge-preserving filter  High computational complexity 2. Separable implementation:  Good approximation of the original filter  Linear computational complexity 3. Application to video preprocessing:  Noise reduction  Better compressed video
  • 3. Gaussian filtering revisited  Gaussian filtering: weights depend on distance to the center pixel * = Noisy step edge Gaussian weights Gaussian filtered result  Adaptive filtering: avoid filtering across edges . = Noisy step edge Edge-preserving weights Edge-presered filtered result
  • 4. Gaussian vs bilateral filtering  Gaussian filtering: weights depend on distance to the center pixel * = Noisy step edge Gaussian weights Gaussian filtered result  Bilateral filtering: weights depend on both spatial closeness and photometric similarity . = Noisy step edge Edge-preserving weights Edge-presered filtered result
  • 5. How bilateral filtering works?  Every sample is a weighted average of its neighbors: 1 O ( s0 ) = ∑ w ( s, s0 ) . I ( s ) K s Local modes  The weight w = w s . wt is product of two Gaussian weights: after bilateral  − ( s − s0 ) 2  filtering  Spatial proximity: w s = exp    2σ s2   − ( I ( s ) − I ( s0 ) ) 2   Tonal similarity: wt = exp    2σ t 2   
  • 6. Example: Gaussian filtering Noisy input: PSNR = 39.1 dB Gaussian filtered: PSNR = 67.9 dB
  • 7. Example: Bilateral filtering Noisy input: PSNR = 39.1 dB Bilateral filtered: PSNR = 41.6 dB
  • 8. Computational complexity  Bilateral filtering kernel is space-variant → complexity is: O ( Nm d ) where N: number of pixels in the image m: size of filtering kernel (m ≈ 7 is good enough) d: image dimensionality  Previous attempt for fast bilateral filtering: Piecewise-linear: approximate bilateral filtering with M Gaussian filtering - Durand & Dorsey (SIGGRAPH 2002) → complexity is: O ( N log( N ) . M ) M ≈ 17 for 8-bit images  Our approach: separable bilateral filtering O ( Nmd )
  • 9. Is bilateral filter separable?  Gaussian filter is space-invariant and separable: g(x) g(y) ∑ w ( s − s ) . I ( s) s 0 O( s ) = s * = ∑ w (s − s ) 0 s 0 s y x w s ( s ) = gauss( x ) . gauss( y ) Kernel center  Bilateral filter is NOT separable: I(s0)  Space-variant kernel due to local I(s) intensity dependency  However, even a highly non-linear filter like median filter is approximately separable (Narendra – PAMI 1981)
  • 10. Separable bilateral filtering result  Separable bilateral filtering is a good approximation of full kernel filtering: noisy Erika σnoise = 10 bilateral filtered in x-dimension followed by y-dimension filtering
  • 11. Separable bilateral filtering result  Separable bilateral filtering is a good approximation of full kernel filtering: noisy Erika σnoise = 10 bilateral filtered in x-dimension followed by y-dimension filtering Image size Brute-force Durrand 2002 Separable Aniso. diffusion 256x256 4.46 0.37 0.21 2.76 512x512 17.88 1.59 0.89 26.02 61x61x61 5.29 4.20 0.45 5.25 256x256x212 56 min Out of memory 50.3 Out of memory
  • 12. How separable bilateral filtering works?  Pros: extremely fast (fixed spatial weight + LUT for tonal weight)  Cons: effective filtering kernel is a slightly distorted
  • 13. Performance of separable bilateral filtering  Very good approximation of  Almost linear execution time the full-kernel filter per pixel
  • 14. MPEG-1 Foreman with bilateral preprocessing 16 without preprocessing with 3x3x3 full-kernel with 9x9x5 separable 15 RMSE 14 13 0 400 800 1200 bit-rate (K bits/s)  Better RMSE is achieved with separable bilateral filtering compared to full-kernel bilateral filtering with the same computation requirement
  • 15. Less artifact with Bilateral preprocessing
  • 16. Conclusions  Separable implementation of bilateral filtering:  Very good approximation of the original filter  Much faster than the original or other approximations  Applications in video preprocessing:  Improved quality of compressed video  Reduced processing time → real-time possibility
  • 17. Literature  C. Tomasi and R. Manduchi, Bilateral fitering for gray and color images, Proc. of ICCV, USA, 1998, 839-846.  T.Q. Pham and L.J. van Vliet, Separable bilateral filtering for fast video preprocessing, Proc. of ICME’05.  F. Durrand and J. Dorsey, Fast bilateral filtering for the display of high dynamic range images, Proc. of SIGGRAPH’02, 2002, 844-847.  P. Perona and J. Malik, Scale-space filtering and edge detection using anisotropic diffusion, PAMI, vol. 12, no. 7, 1990, 629-639.  R. v.d. Boomgaard and J. v.d. Weijer. On the equivalence of local- mode finding, robust estimation and mean-shift analysis as used in early vision tasks. In Proc. of ICPR, pages 927-930, 2002.