bm3dsharp.m

BM3D去噪算法的实现和相关文档

function [y_hat] = BM3DSHARP(z, sigma, alpha_sharp, profile, print_to_screen)
%
%  Joint sharpening and denoising with BM3D. This is implementation of the 
%  BM3D-SH3D sharpening method that is developed in:
%
%  [1] K. Dabov, A. Foi, V. Katkovnik, and K. Egiazarian, "Joint image
%   sharpening and denoising by 3D transform-domain collaborative filtering," 
%   Proc. 2007 Int. TICSP Workshop Spectral Meth. Multirate Signal Process.,
%   SMMSP 2007, Moscow, Russia, September 2007.
%
%  FUNCTION INTERFACE:
%
%  [ysharp] = BM3DSHARP(z, sigma, alpha_sharp, profile, print_to_screen)
%
%  The function can work without any of the input arguments, hence they are
%  optional!
%
%  INPUTS (OPTIONAL):
%
%        1) z (matrix, size MxN)       : Input image (noisy and with poor contrast)
%        2) sigma (double)             : Noise (IF ANY noise) standard deviation (signal assumed
%                                          in the range [0, 255])
%        3) alpha_sharp (double)       : Sharpening parameter (default: 1.5):
%                                         (1,inf) -> sharpen
%                                          1      -> no sharpening
%                                         (0,1)   -> de-sharpen
%        4) profile (char vector)      : 'lc' --> fast 
%                                        'np' --> normal (default)
%        5) print_to_screen (boolean)  : 0 --> do not print output
%                                          information (and do not plot figures)
%                                        1 --> print figures (default)
%
%   OUTPUTS:
%        1) ysharp (matrix, size MxN)  : Sharpened image (in the range [0,1])
%
%  BASIC USAGE EXAMPLES:     
%     
%     sigma = 10;
%     z = im2double(imread('cameraman.tif'));
%     z = z + (sigma/255)*randn(size(z));
%     alpha_sharp = 1.3;
%     [ysharp] = BM3DSHARP(z, sigma, alpha_sharp);
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Copyright