BM3D滤波算法+BM3D_image，可实现

BM3D去噪算法+图像实例，matlab源码，亲测可实现

function [PSNR y_est] = BM3D（y z sigma profile print_to_screen）

%  代码的模板基于http://www.cs.tut.fi/~foi/GCF-BM3D此网址提供的BM3D.m，
%  此文件将原文件中bm3d_thr和bm3d_wiener两个函数用matlab语言进行了实现，并附在最后，
%  由于原文件有些参数并未体现在论文思想中，这一情况同样体现在原文件的bm3d_thr和bm3d_wiener两个函数上，
%  又或者是我理解的偏差，导致编写的bm3d_thr和bm3d_wiener两个函数有部分参数并未用到，不过不影响函数使用，
%  另外，在调试过程中也为了提高PSNR，修改了部分参数，
%  此文件使用时，将预处理的图像与此文件放在同一个文件夹下，修改image_name（100行左右）使其与图像文件名一致，
%  可通过修改约位于128行的sigma修改方差
%  另外，也许是由于算法问题，此文件处理图像较慢，256x256需要2到3分钟，512x512需要8到15分钟，还请注意
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  BM3D is algorithm for attenuation of additive white Gaussian noise from 
%  grayscale images. This algorithm reproduces the results from the article:
%
%  [1] K. Dabov A. Foi V. Katkovnik and K. Egiazarian “Image denoising 
%   by sparse 3D transform-domain collaborative filtering“ 
%   IEEE Transactions on Image Processing December 2006 in review
%   preprint at http://www.cs.tut.fi/~foi/GCF-BM3D.
%
%  FUNCTION INTERFACE:
%
%  [PSNR y_est] = BM3D（y z sigma profile print_to_screen）
%
%  ! The function can work without any of the input arguments 
%   in which case the internal default ones are used !
% 
%  BASIC USAGE EXAMPLES:
%
%     Case 1） Using the default parameters （i.e. image name sigma etc.）
% 
%      [PSNR y_est] = BM3D;
% 
%     Case 2） Using an external noisy image:
%
%      % Read a grayscale image and scale its intensities in range [01]
%      y = im2double（imread（‘Cameraman256.png‘））; 
%      % Generate the same seed used in the experimental results of [1]
%      randn（‘seed‘ 0）;
%      % Standard deviation of the noise --- corresponding to intensity 
%      %  range [0255] despite that the input was scaled in [01]
%      sigma = 25;
%      % Add the AWGN with zero mean and standard deviation ‘sigma‘
%      z = y + （sigma/255）*randn（size（y））;
%      % Denoise ‘z‘. The denoised image is ‘y_est‘ and ‘NA = 1‘ because 
%      %  the true image was not provided
%      [NA y_est] = BM3D（1 z sigma）; 
%      % Compute the putput PSNR
%      PSNR = 10*log10（1/mean（（y（:）-y_est（:））.^2））
%      % show the noisy image ‘z‘ and the denoised ‘y_est‘
%      figure; imshow（z）;   
%      figure; imshow（y_est）;
% 
%     Case 3） If the original image y is provided as the first input 
%      argument then some additional information is printed （PSNRs 
%      figures etc.）. That is “[NA y_est] = BM3D（1 z sigma）;“ in the
%      above code should be replaced with:
% 
%      [PSNR y_est] = BM3D（y z sigma）;
% 
% 
%  INPUT ARGUMENTS （OPTIONAL）:
%
%     1） y （matrix M x N）: Noise-free image （needed for computing PSNR）
%                           replace with the scalar 1 if not available.
%     2） z （matrix M x N）: Noisy image （intensities in range [01] or [0255]）
%     3） sigma （double）  : Std. dev. of the noise （corresponding to intensities
%                          in range [0255] even if the range of z is [01]）
%     4） profile （char）  : ‘np‘ --> 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-04-08 20:21  BM3D\
     文件        6148  2019-04-08 20:21  BM3D\.DS_Store
     目录           0  2019-04-08 20:22  __MACOSX\
     目录           0  2019-04-08 20:22  __MACOSX\BM3D\
     文件         120  2019-04-08 20:21  __MACOSX\BM3D\._.DS_Store
     文件       32050  2019-04-03 20:02  BM3D\BM3D.m
     文件         232  2019-04-03 20:02  __MACOSX\BM3D\._BM3D.m
     目录           0  2019-04-03 20:00  BM3D\BM3D_images\
     文件      177762  2005-11-14 08:14  BM3D\BM3D_images\boat.png
     目录           0  2019-04-08 20:22  __MACOSX\BM3D\BM3D_images\
     文件         212  2005-11-14 08:14  __MACOSX\BM3D\BM3D_images\._boat.png
     文件       33216  2006-12-30 12:59  BM3D\BM3D_images\montage.png
     文件         212  2006-12-30 12:59  __MACOSX\BM3D\BM3D_images\._montage.png
     文件      429132  2005-12-08 11:21  BM3D\BM3D_images\image_F16_512rgb.png
     文件         212  2005-12-08 11:21  __MACOSX\BM3D\BM3D_images\._image_F16_512rgb.png
     文件      476195  2005-10-29 12:18  BM3D\BM3D_images\image_Lena512rgb.png
     文件         212  2005-10-29 12:18  __MACOSX\BM3D\BM3D_images\._image_Lena512rgb.png
     文件      185727  2005-11-14 08:14  BM3D\BM3D_images\barbara.png
     文件         212  2005-11-14 08:14  __MACOSX\BM3D\BM3D_images\._barbara.png
     文件      110937  2005-10-29 12:19  BM3D\BM3D_images\image_House256rgb.png
     文件         212  2005-10-29 12:19  __MACOSX\BM3D\BM3D_images\._image_House256rgb.png
     文件       37559  2006-04-17 11:20  BM3D\BM3D_images\Cameraman256.png
     文件         212  2006-04-17 11:20  __MACOSX\BM3D\BM3D_images\._Cameraman256.png
     文件       34985  2005-11-14 08:14  BM3D\BM3D_images\house.png
     文件         212  2005-11-14 08:14  __MACOSX\BM3D\BM3D_images\._house.png
     文件      153785  2006-04-17 11:20  BM3D\BM3D_images\Lena512.png
     文件         212  2006-04-17 11:20  __MACOSX\BM3D\BM3D_images\._Lena512.png
     文件      184561  2006-06-06 20:48  BM3D\BM3D_images\fingerprint.png
     文件         212  2006-06-06 20:48  __MACOSX\BM3D\BM3D_images\._fingerprint.png
     文件      209817  2006-12-30 12:53  BM3D\BM3D_images\man.png
     文件         212  2006-12-30 12:53  __MACOSX\BM3D\BM3D_images\._man.png
............此处省略13个文件信息