Wavelet

Development Platform:

Matlab

den1_9.m：Code Content

							function X = den1_9(x, wname, n, thr, thr1)
% 基于综合阈值的图像去噪方法
% 将thr1和thr之间的线性处理方式改成指数的
% x' = a^(x-b) - 1
% a = (λ + 1)^(1 / (λ-λ1))
% b = λ1
% thr1 = 0.6 * thr;
a = (thr + 1) ^ (1/(thr-thr1));
b = thr1;
[C, S] = wavedec2(x, n, wname);                     % 对图像进行小波分解
dcoef = C( prod(S(1, :)) + 1 : end);                % 高频部分系数
ind = find( abs(dcoef) < thr1) + prod(S(1, :));     % 小于thr1的系数
C(ind) = 0;                                         %   直接置零
ind = find( abs(dcoef) >= thr1 & abs(dcoef) < thr )...
     + prod(S(1, :));                               % 大于thr1小于thr的系数
C(ind) = sign(C(ind)) .* ...
    (a .^( abs(C(ind)) - b) - 1);                   % 进行缩减
% 重构图像
X = waverec2(C, S, wname);

						