
Abstract |
Software |
Results |
People |
Related work |
Publications |
We propose a novel image denoising strategy based on an enhanced sparse representation in transform-domain. The enhancement of the sparsity is achieved by grouping similar 2D image fragments (e.g. blocks) into 3D data arrays which we call "groups".
Collaborative filtering is a special procedure developed to deal with these 3D groups. We realize it using the three successive steps: 3D transformation of 3D group, shrinkage of transform spectrum, and inverse 3D transformation. The result is a 3D estimate that consists of the jointly filtered grouped image blocks. By attenuating the noise, the collaborative filtering reveals even the finest details shared by grouped blocks and at the same time it preserves the essential unique features of each individual block. The filtered blocks are then returned to their original positions.
Because these blocks are overlapping, for each pixel we obtain many different estimates which need to be combined. Aggregation is a particular averaging procedure which is exploited to take advantage of this redundancy.
A significant improvement is obtained by a specially developed collaborative Wiener filtering.
We develop algorithms based on this novel denoising strategy. The experimental results presented here demonstrate that the developed methods achieve state-of-the-art denoising performance in terms of both peak signal-to-noise ratio and subjective visual quality.
before you proceed with downloading any of the files.
Download zipped BM3D MATLAB software
|
|
Download zipped test images |
|
Download zipped test sequences |
| σ¹/PSNR | House 256x256 |
Peppers 512x512 |
Lena 512x512 |
Baboon 512x512 |
F16 512x512 |
Kodak image 1 768x512 |
Kodak image 2 768x512 |
Kodak image 3 768x512 |
| 5 / 34.16 | 38.97 | 36.82 | 37.82 | 35.25 | 39.68 | 39.37 | 40.34 | 42.48 |
| 10 / 28.13 | 36.23 | 33.78 | 35.22 | 30.64 | 36.69 | 34.71 | 36.57 | 39.05 |
| 15 / 24.61 | 34.85 | 32.60 | 33.94 | 28.39 | 35.00 | 32.13 | 34.59 | 37.05 |
| 20 / 22.11 | 33.84 | 31.83 | 33.02 | 26.97 | 33.77 | 30.40 | 33.33 | 35.64 |
| 25 / 20.18 | 33.03 | 31.20 | 32.27 | 25.96 | 32.78 | 29.13 | 32.44 | 34.54 |
| 30 / 18.59 | 32.33 | 30.62 | 31.59 | 25.15 | 31.94 | 28.14 | 31.72 | 33.59 |
| 35 / 17.25 | 31.58 | 30.00 | 30.91 | 24.46 | 31.13 | 27.31 | 31.07 | 32.63 |
| 50 / 14.16 | 30.22 | 28.68 | 29.72 | 23.14 | 29.41 | 25.65 | 29.77 | 31.15 |
| 75 / 10.63 | 28.33 | 27.12 | 28.19 | 21.71 | 27.60 | 24.14 | 28.33 | 29.52 |
| σ¹/PSNR | Salesman 288x352@50 |
Tennis 240x352@150 |
Flower garden 240x352@150 |
Miss America 288x360@150 |
Costguard 144x176@300 |
Foreman 288x352@300 |
Bus 288x352@150 |
Bicycle 576x720@30 |
| 5 / 34.16 | 40.44 | 38.47 | 36.49 | 41.58 | 38.25 | 39.77 | 37.55 | 40.89 |
| 10 / 28.13 | 37.21 | 34.68 | 32.11 | 39.61 | 34.78 | 36.46 | 33.32 | 37.62 |
| 15 / 24.61 | 35.44 | 32.63 | 29.81 | 38.64 | 33.00 | 34.64 | 31.05 | 35.67 |
| 20 / 22.11 | 34.04 | 31.20 | 28.24 | 38.85 | 31.71 | 33.30 | 29.57 | 34.18 |
| 25 / 20.18 | 32.79 | 30.11 | 27.00 | 37.10 | 30.62 | 32.19 | 28.48 | 32.90 |
| 30 / 18.59 | 31.68 | 29.22 | 25.89 | 36.41 | 29.68 | 31.27 | 27.59 | 31.77 |
| 35 / 17.25 | 30.72 | 28.56 | 25.16 | 36.87 | 28.92 | 30.56 | 26.91 | 30.85 |
¹ The noisy images/videos were created by adding zero-mean white Gaussian noise with the following MATLAB commands:
randn('seed', 0);
noisy = original + sigma*randn(size(original));
| σ | Salesman 288x352@50 |
Tennis 240x352@150 |
Flower garden 240x352@150 |
|||
| Noisy | Denoised | Noisy | Denoised | Noisy | Denoised | |
| 15 | 24.71 | 35.44 | 24.78 | 32.61 | 24.88 | 29.78 |
| 20 | 22.32 | 34.07 | 22.33 | 31.18 | 22.46 | 28.21 |
Aram Danielyan
Kostadin Dabov
Alessandro Foi
Vladimir Katkovnik
Karen Egiazarian
Compressed Sensing Image Reconstruction, Image Upsampling, and Image/Video Super-Resolution via Recursive Spatially Adaptive Filtering
A preliminary version of BM3D using exclusively the DFT
Shape-Adaptive Transforms Filtering (Pointwise SA-DCT algorithm).
Katkovnik, V., A. Foi, K. Egiazarian, and J. Astola, “From local kernel to nonlocal multiple-model image denoising”, Int. J. Computer Vision, vol. 86, no. 1, pp. 1-32, January 2010.
doi:10.1007/s11263-009-0272-7
K. Dabov, A. Foi, V. Katkovnik, and K. Egiazarian, “BM3D Image Denoising with Shape-Adaptive Principal Component Analysis”, Proc. Workshop on Signal Processing with Adaptive Sparse Structured Representations (SPARS'09), Saint-Malo, France, April 2009.
Poster